A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat

Area of Science:

• Computational Chemistry and Physics
• Statistical Mechanics
• Molecular Dynamics Simulations

Background:

• Path Integral Molecular Dynamics (PIMD) is a powerful quantum simulation technique.
• Efficient thermostatting is crucial for accurate PIMD simulations.
• Existing Langevin equation thermostats in PIMD can limit simulation efficiency.

Purpose of the Study:

• To introduce a novel, simple algorithm for thermostatting PIMD using the Langevin equation.
• To improve the efficiency and accuracy of PIMD simulations.
• To compare the performance of the new algorithm against existing methods.

Main Methods:

• Development of a new thermostatting algorithm for PIMD based on the Langevin equation.
• Utilizing the staging transformation of path integral beads.
• Employing optimum friction coefficients derived from the free particle limit.
• Performing PIMD simulations on realistic systems: H2O molecule, liquid para-hydrogen, and liquid water.

Main Results:

• The new thermostatting algorithm significantly increases the accessible time interval in PIMD simulations by a factor of 4-6 or more for equivalent accuracy.
• The algorithm demonstrates improved efficiency compared to the standard path integral Langevin equation thermostat.
• Error analysis confirms the robustness of the new algorithm across different systems.

Conclusions:

• The novel simple algorithm offers a substantial improvement in the efficiency of path integral molecular dynamics simulations.
• This method provides a more accurate and computationally feasible approach for studying quantum systems.
• The findings are applicable to a wide range of molecular systems, including liquids and molecules.