Abstract:Direct dynamics simulations of chemical reactions typically require the selection of a method for generating the potential energy surfaces and a method for the dynamical propagation of the nuclei on these surfaces. The nuclear-electronic orbital (NEO) framework avoids this Born-Oppenheimer separation by treating specified nuclei on the same level as the electrons with wave function methods or density functional theory (DFT). The NEO approach is particularly applicable to proton, hydride, and proton-coupled electron transfer reactions, where the transferring proton(s) and all electrons are treated quantum mechanically. In this manner, the zero-point energy, density delocalization, and anharmonicity of the transferring protons are inherently and efficiently included in the energies, optimized geometries, and dynamics.This Account describes how various NEO methods can be used for direct dynamics simulations on electron-proton vibronic surfaces. The strengths and limitations of these approaches are discussed, and illustrative examples are presented. The NEO-DFT method can be used to simulate chemical reactions on the ground state vibronic surface, as illustrated by the application to hydride transfer in C4H9+. The NEO multistate DFT (NEO-MSDFT) method is useful for simulating ground state reactions in which the proton density becomes bilobal during the dynamics, a characteristic of hydrogen tunneling, as illustrated by proton transfer in malonaldehyde. The NEO time-dependent DFT (NEO-TDDFT) method produces excited electronic, vibrational, and vibronic surfaces. The application of linear-response NEO-TDDFT to H2 and H3+, as well as the partially and fully deuterated counterparts, shows that this approach produces accurate fundamental vibrational excitation energies when all nuclei and all electrons are treated quantum mechanically. Moreover, when only specified nuclei are treated quantum mechanically, this approach can be used to optimize geometries on excited state vibronic surfaces, as illustrated by photoinduced single and double proton transfer systems, and to conduct adiabatic dynamics on these surfaces. The real-time NEO-TDDFT method provides an alternative approach for simulating nonequilibrium nuclear-electronic dynamics of such systems. These various NEO methods can be combined with nonadiabatic dynamics methods such as Ehrenfest and surface hopping dynamics to include the nonadiabatic effects between the quantum and classical subsystems. The real-time NEO-TDDFT Ehrenfest dynamics simulation of excited state intramolecular proton transfer in o-hydroxybenzaldehyde illustrates the power of this type of combined approach. The field of multicomponent quantum chemistry is in the early stages, and the methods discussed herein provide the foundation for a wide range of promising future directions to be explored. An appealing future direction is the expansion of the real-time NEO-TDDFT method to describe the dynamics of all nuclei and electrons on the same level. Direct dynamics simulations using NEO wave function methods such as equation-of-motion coupled cluster or multiconfigurational approaches are also attractive but computationally expensive options. The further development of NEO direct dynamics methods will enable the simulation of the nuclear-electronic dynamics for a vast array of chemical and biological processes that extend beyond the Born-Oppenheimer approximation.

Time-Resolved Chemical Bonding Structure Evolution by Direct-Dynamics Chemical Simulations

Direct Chemical Dynamics Simulations of H3+ + CO Bimolecular Reaction

Towards simultaneous imaging of ultrafast nuclear and electronic dynamics in small molecules

Molecular dynamics simulation of a chemical reaction in solution

Tabletop imaging of structural evolutions in chemical reactions

Modeling chemical reactions in classical molecular dynamics simulations

From Ab Initio Potential Energy Surfaces to State-Resolved Reactivities: X + H2O ↔ HX + OH [X = F, Cl, and O(3P)] Reactions

Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms using Molecular Configuration Transformer

Halogen Bond Structure and Dynamics from Molecular Simulations

Chemical bond reorganization in intramolecular proton transfer revealed by ultrafast X-ray photoelectron spectroscopy

Spectroscopic and Structural Probing of Excited-State Molecular Dynamics with Time-Resolved Photoelectron Spectroscopy and Ultrafast Electron Diffraction

Bond dissociation dynamics of single molecules on Ag(111)

Direct Dynamics with Nuclear-Electronic Orbital Density Functional Theory

Ultrafast chemical kinetics: Elementary chemical act

Above-Threshold Ionization and Laser-Induced Electron Diffraction in Diatomic Molecules

Dynamics and Kinetics Study of "In-Water" Chemical Reactions by Enhanced Sampling of Reactive Trajectories.

Electronic Force Density Fields: Insights into Partial Bonds, Transition States, and Chemical Structure Evolution

Femtosecond intramolecular rearrangement of the CH3NCS radical cation

Dynamic hetero-metallic bondings visualized by sequential atom imaging

Real-time Observation of Ultrafast Molecular Rotation in Weakly Bound Dimers

Dynamic crystal structure of a molecular framework