Quantum Thermal Bath for Path Integral Molecular Dynamics Simulation.

Area of Science:

• Computational Chemistry and Physics
• Quantum Mechanics
• Materials Science

Background:

• Standard molecular dynamics (MD) simulations often neglect nuclear quantum effects.
• Quantum thermal bath (QTB) methods approximate nuclear quantum effects within MD.
• Path integral molecular dynamics (PIMD) provides exact quantum results but is computationally expensive.

Purpose of the Study:

• To combine QTB and PIMD methods for improved accuracy and efficiency in simulating quantum systems.
• To develop a more accurate QTB-PIMD approach for strongly anharmonic systems.
• To validate the new method on benchmark systems and a relevant material.

Main Methods:

• Development of a novel power spectral density for the random force in the QTB method.
• Introduction of a modified centroid-virial estimator for kinetic energy, tailored for QTB-PIMD.
• Application of the combined QTB-PIMD method to a 1D double-well system, a ferroelectric transition, and hydrogen atom distribution in fuel cell materials.

Main Results:

• The developed QTB-PIMD method accurately captures nuclear quantum effects.
• The approach demonstrates significant improvements for strongly anharmonic systems compared to standalone QTB-MD.
• The method achieves exact results with reduced computational cost compared to standard PIMD.

Conclusions:

• The combined QTB-PIMD method offers a computationally efficient and accurate way to include nuclear quantum effects.
• This hybrid approach overcomes limitations of both QTB-MD and PIMD for anharmonic systems.
• QTB-PIMD is a promising technique for accurate simulations in condensed matter physics and materials science.