Search

The ab initio molecular dynamics (AIMD) method provides a computational route for the real-time simulation of reactive chemistry. An often-overlooked capability of this approach is the opportunity to examine the electronic evolution of a chemical system. For AIMD trajectories based on Hartree-Fock or density functional theory methods, the real-time evolution of single-particle molecular orbitals (MOs) can provide detailed insights into the time-dependent electronic structure of molecules. The evolving electronic Hamiltonians at each MD step pose problems for tracking and visualizing a given MO's character, ordering, and associated phase throughout an MD trajectory, however. This report presents and assesses a simple algorithm for correcting these deficiencies by exploiting similarity projections of the electronic structure between neighboring MD steps. Two aspects bring this analysis beyond a simple step-to-step projection scheme. First, the challenging case of coincidental orbital degeneracies is resolved via a quadrupole-field perturbation that nonetheless rigorously preserves energy conservation. Second, the resulting orbitals are shown to evolve adiabatically, in spite of the "preservation of character" concept that undergirds a projection of neighboring steps' MOs. The method is tested on water clusters, which exhibit considerable dynamic degeneracies, as well as a classic organic nucleophilic substitution reaction, in which the adiabatic evolution of the bonding orbitals clarifies textbook interpretations of the electronic structure during this reactive collision.

The ab initio molecular dynamics (AIMD) method provides a computational route for the real-time simulation of reactive chemistry. An often-overlooked capability of this approach is the opportunity to examine the electronic evolution of a chemical system. For AIMD trajectories based on Hartree-Fock or density functional theory methods, the real-time evolution of single-particle molecular orbitals (MOs) can provide detailed insights into the time-dependent electronic structure of molecules. The evolving electronic Hamiltonians at each MD step pose problems for tracking and visualizing a given MO's character, ordering, and associated phase throughout an MD trajectory, however. This report presents and assesses a simple algorithm for correcting these deficiencies by exploiting similarity projections of the electronic structure between neighboring MD steps. Two aspects bring this analysis beyond a simple step-to-step projection scheme. First, the challenging case of coincidental orbital degeneracies is resolved via a quadrupole-field perturbation that nonetheless rigorously preserves energy conservation. Second, the resulting orbitals are shown to evolve adiabatically, in spite of the "preservation of character" concept that undergirds a projection of neighboring steps' MOs. The method is tested on water clusters, which exhibit considerable dynamic degeneracies, as well as a classic organic nucleophilic substitution reaction, in which the adiabatic evolution of the bonding orbitals clarifies textbook interpretations of the electronic structure during this reactive collision.