Path integral Liouville dynamics for thermal equilibrium systems

Area of Science:

• Quantum dynamics
• Computational chemistry
• Statistical mechanics

Background:

• Simulating condensed phase molecular systems requires accurate methods to capture quantum effects.
• Existing methods may struggle with strongly anharmonic potentials or nonlinear operators.
• The Wigner phase space offers a framework for developing new quantum dynamics methods.

Purpose of the Study:

• Introduce and validate a novel imaginary time path integral method: Path Integral Liouville Dynamics (PILD).
• Assess PILD's ability to conserve quantum canonical distributions.
• Evaluate PILD's capacity to recover exact thermal correlation functions.

Main Methods:

• Derivation of PILD from equilibrium Liouville dynamics in the Wigner phase space.
• Numerical testing of PILD using a simple Langevin thermostat.
• Application to two strongly anharmonic model problems.

Main Results:

• PILD successfully conserves the quantum canonical distribution.
• PILD accurately recovers exact thermal correlation functions for nonlinear operators.
• The method demonstrated effectiveness for anharmonic systems.

Conclusions:

• Path Integral Liouville Dynamics (PILD) is a promising new approach for condensed phase molecular simulations.
• PILD offers a general method applicable to various potential energy surfaces.
• The method shows potential for accurate simulations in classical, high-temperature, and harmonic limits.