Atomistic Simulation of Adiabatic Reactive Processes Based on Multi-State Potential Energy Surfaces

The adiabatic reactive molecular dynamics (ARMD) method provides a framework to study chemical reactions using molecular dynamics simulations with minimal computational overhead. Here, ARMD is generalized to an arbitrary reactive process between two states in which reactants and products can be treated by an atomistic force field. The implementation is described, and the method is applied to two systems: the kinetics of NO rebinding to myoglobin (Mb) as a validation system and the conformational transition in neuroglobin (Ngb) which explores the full functionality of ARMD. For MbNO, the nonexponential kinetics observed both in experiment and earlier ARMD studies is reproduced. Furthermore, the sensitivity of the results with respect to the asymptotic separation between the two potential energy surfaces (NO bound and unbound) is studied.