Quasiclassical Correlation Functions from the Wigner Density Using the Stability Matrix.

Area of Science:

• Quantum mechanics
• Statistical mechanics
• Computational chemistry

Background:

• Classical trajectory calculations often require sampling from quantized phase space distributions.
• The Weyl-Wigner transform is commonly used but computationally challenging to construct and sample.
• Accurate initial conditions are crucial for calculating time correlation functions.

Purpose of the Study:

• To develop a more efficient method for sampling quantized phase space distributions.
• To simplify the calculation of time correlation functions in classical trajectory simulations.
• To reduce the computational cost associated with quantum corrections in classical dynamics.

Main Methods:

• Transferring operator dependence from the phase space distribution to the dynamics.
• Augmenting classical equations of motion with stability matrix differential equations.
• Employing a local harmonic approximation for dynamical derivatives.

Main Results:

• Enables sampling from the simpler Wigner phase space density instead of the Weyl-Wigner transform.
• Significantly reduces computational cost for nonlinear operator correlation functions.
• Local harmonic approximation quantitatively captures quasiclassical results for dissipative systems.

Conclusions:

• The proposed method offers a computationally efficient alternative for incorporating zero-point energy in classical trajectory calculations.
• The local harmonic approximation shows promise, particularly for systems in dissipative environments.
• This approach facilitates more accessible and accurate quantum-classical simulations.