Development of semiclassical molecular dynamics simulation method.

Area of Science:

• Quantum mechanics
• Chemical dynamics
• Computational chemistry

Background:

• Quantum mechanical effects like tunneling are vital in chemical and biological processes.
• Molecular dynamics (MD) methods traditionally rely on classical mechanics, neglecting these quantum effects.

Purpose of the Study:

• To develop a method for incorporating quantum mechanical tunneling into classical molecular dynamics (MD).
• To enhance the accuracy of MD simulations for chemical reactions by including quantum effects.

Main Methods:

• Detecting caustics along classical trajectories to identify classically forbidden regions.
• Determining the optimal tunneling path with minimum action starting from caustics.
• Estimating the real phase associated with quantum tunneling.

Main Results:

• A numerical demonstration using the collinear reaction O + HCl → OH + Cl successfully illustrates the proposed method.
• The method provides a way to incorporate tunneling effects into classical MD simulations.
• The approach is generalizable to on-the-fly ab initio calculations.

Conclusions:

• The developed method offers a significant advancement in simulating chemical reactions with quantum tunneling.
• Future work can involve developing new semiclassical MD methods by incorporating nonadiabatic transitions using Zhu-Nakamura theory.