Molecular Dynamics Simulations on Relaxed Reduced-Dimensional Potential Energy Surfaces

Area of Science:

• Computational Chemistry
• Chemical Dynamics
• Theoretical Chemistry

Background:

• High-level ab initio calculations enable accurate potential energy surface (PES) construction for molecular dynamics (MD) simulations.
• Full-dimensional PESs are computationally expensive for larger molecules, limiting reaction dynamics studies.
• Reduced-dimensional PESs offer a computationally feasible alternative by focusing on essential reactive coordinates.

Purpose of the Study:

• To develop and present a new molecular dynamics (MD) implementation for utilizing relaxed reduced-dimensional potential energy surfaces (PESs).
• To enable standard microcanonical (NVE) and canonical (NVT) MD simulations on these reduced PESs.
• To validate the accuracy of the new MD method for reaction dynamics.

Main Methods:

• Developed a novel MD algorithm to generate trajectories on relaxed reduced-dimensional PESs.
• Applied the method to study the pyramidal inversion of ammonia (NH3).
• Validated results against ab initio MD simulations, full-dimensional PES MD simulations, and experimental data.

Main Results:

• The new MD implementation successfully utilizes relaxed reduced-dimensional PESs for NVE and NVT simulations.
• Simulations of ammonia's pyramidal inversion on a 3D PES showed good agreement with ab initio MD, full-dimensional PES MD, and experimental results.
• The method demonstrates the feasibility of studying reaction dynamics on reduced-dimensional PESs.

Conclusions:

• The presented MD implementation provides a computationally efficient and accurate approach for reaction dynamics studies of larger molecular systems.
• Reduced-dimensional PESs, when used with appropriate MD algorithms, can effectively model complex chemical processes.
• This work bridges the gap between theoretical PES construction and practical MD simulations for systems where full-dimensional approaches are intractable.