Computing thermal Wigner densities with the phase integration method
- Département de Chimie, Ecole Normale Supérieure-PSL Research University, 24, rue Lhomond, 75005 Paris (France)
- Dipartimento di Fisica, Università di Roma “La Sapienza” Piazzale A. Moro 5, 00185 Rome (Italy)
We discuss how the Phase Integration Method (PIM), recently developed to compute symmetrized time correlation functions [M. Monteferrante, S. Bonella, and G. Ciccotti, Mol. Phys. 109, 3015 (2011)], can be adapted to sampling/generating the thermal Wigner density, a key ingredient, for example, in many approximate schemes for simulating quantum time dependent properties. PIM combines a path integral representation of the density with a cumulant expansion to represent the Wigner function in a form calculable via existing Monte Carlo algorithms for sampling noisy probability densities. The method is able to capture highly non-classical effects such as correlation among the momenta and coordinates parts of the density, or correlations among the momenta themselves. By using alternatives to cumulants, it can also indicate the presence of negative parts of the Wigner density. Both properties are demonstrated by comparing PIM results to those of reference quantum calculations on a set of model problems.
- OSTI ID:
- 22419817
- Journal Information:
- Journal of Chemical Physics, Vol. 141, Issue 8; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Path integral Liouville dynamics for thermal equilibrium systems
Disentangling covariant Wigner functions for chiral fermions