Langevin integration for isothermal-isobaric condition with a large time step

Area of Science:

• Computational chemistry
• Molecular dynamics simulations
• Statistical mechanics

Background:

• Langevin integration is crucial for molecular dynamics (MD) simulations.
• Existing methods face challenges in accurate temperature and pressure evaluation, especially with large time steps.
• Inaccuracies in momentum space affect pressure calculations in conventional Langevin dynamics.

Purpose of the Study:

• To develop an accurate method for evaluating temperature and pressure in Langevin integration.
• To address inaccuracies in pressure evaluation inherent in current Langevin dynamics approaches.
• To improve the stability and reliability of isothermal-isobaric MD simulations.

Main Methods:

• Building upon the Leimkuhler and Matthews approach for Langevin integration.
• Modifying the calculation of kinetic energy for pressure evaluation using half-time step momentum.
• Performing molecular dynamics simulations with the refined pressure evaluation method.

Main Results:

• The proposed method significantly improves the quality of configuration space compared to other Langevin dynamics methods.
• Calculating kinetic energy with half-time step momentum reduces numerical errors in pressure evaluation.
• Achieved improved accuracy and stability in isothermal-isobaric MD simulations with a time step of 5 fs.

Conclusions:

• The refined pressure evaluation method enhances the accuracy and stability of isothermal-isobaric molecular dynamics simulations.
• Utilizing half-time step momentum effectively mitigates errors associated with large time steps.
• This advancement allows for more reliable simulations using longer time steps, improving computational efficiency.