Gaussian Variational Schemes on Bounded and Unbounded Domains