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.

Transition State Theory within Constrained Nuclear-Electronic Orbital Framework: Accurate Hydrogen Atom Transfer Reaction Rates with Direct Incorporation of Quantum Nuclear Delocalization Effects

Constrained Nuclear-Electronic Orbital Transition State Theory Using Energy Surfaces with Nuclear Quantum Effects

Assessment of Electron-Proton Correlation Functionals for Vibrational Spectra of Shared-Proton Systems by Constrained Nuclear-Electronic Orbital Density Functional Theory

Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory for Multiple Proton Transfer Processes

Nonadiabatic Hydrogen Tunneling Dynamics for Multiple Proton Transfer Processes with Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory

Constrained Nuclear-Electronic Orbital QM/MM Approach for Simulating Complex Systems with Quantum Nuclear Delocalization Effects Incorporated

Nuclear-Electronic Orbital Time-Dependent Configuration Interaction Method

Quantum Suppression of Intramolecular Deuterium Kinetic Isotope Effects in a Pericyclic Hydrogen Transfer Reaction.

Computational Prediction for Singlet- and Triplet-Transition Energies of Charge-Transfer Compounds.

Elucidating the NuclearQuantum Dynamics of Intramolecular Double Hydrogen Transfer in Porphycene

Computation of Rate Coefficients in Solutions Based on Transition State Theory Combined with a Heuristically Corrected Polarizable Continuum Model: Intermolecular Diels-Alder Reactions as Case Studies

Nuclear Quantum Effects in Proton or Hydrogen Transfer

Heavy Atom Tunneling in Organic Reactions at Coupled Cluster Potential Accuracy with a Parallel Implementation of Anharmonic Constant Calculations and Semiclassical Transition State Theory

Accurate transition state generation with an object-aware equivariant elementary reaction diffusion model

Reaction mechanism and isotope effects derived from centroid transition state theory in intramolecular proton transfer reactions

Quantum Monte Carlo for Noncovalent Interactions: An Efficient Protocol Attaining Benchmark Accuracy

Nuclear-Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Direct Dynamics with Nuclear-Electronic Orbital Density Functional Theory

Correlation in Quantum Chemical Calculation and Its Effect on the Uncertainty of Theoretically Predicted Rate Coefficients and Branching Ratios

Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited State Geometry Optimizations and Adiabatic Excitation Energies

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix