Gaussian resolutions for equilibrium density matrices

Abstract

A Gaussian resolution method for the computation of equilibrium density matrices rho(T) for a general multidimensional quantum problem is presented. The variational principle applied to the ``imaginary time'' Schroedinger equation provides the equations of motion for Gaussians in a resolution of rho(T) described by their width matrix, center and scale factor, all treated as dynamical variables. The method is computationally very inexpensive, has favorable scaling with the system size and is surprisingly accurate in a wide temperature range, even for cases involving quantum tunneling. Incorporation of symmetry constraints, such as reflection or particle statistics, is also discussed.

4 pages

References in corpus (1)

• Gaussian Wavepacket Dynamics: semiquantal and semiclassical phase space formalism

Cited by in corpus (7)

• Thermodynamics and equilibrium structure of Ne_38 cluster: Quantum Mechanics versus Classical
• Path integral approach to the Wigner representation of canonical density operators for discrete systems coupled to harmonic baths
• Imaginary time Gaussian dynamics of the Ar_3 cluster
• A fast Variational Gaussian Wave-packet method: Size-induced structural transitions in large neon clusters
• Simulation of open quantum systems
• First-order corrections to semiclassical Gaussian partition functions for clusters of atoms
• Structural information about the ${\text{Ar}}_6$ cluster with the frozen Gaussian imaginary time propagator