Search

Modern electronic structure theory is built around the Born–Oppenheimer approximation and the construction of an electronic Hamiltonian H ̂ e l (X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the well-known theory of electron translation (Γ′) and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phase-space electronic Hamiltonian that depends on both nuclear position and momentum, H ̂ P S (X , P). While classical Born–Oppenheimer dynamics that run along the eigensurfaces of the operator H ̂ e l (X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of H ̂ P S (X , P) can better capture both nuclear and electronic properties (including the elusive electronic momentum studied by Nafie). Moreover, only the latter (as opposed to the former) conserves the total linear and angular momentum in general. [ABSTRACT FROM AUTHOR]

We show that standard Ehrenfest dynamics does not conserve linear and angular momentum when using a basis of truncated adiabatic states. However, we also show that previously proposed effective Ehrenfest equations of motion [M. Amano and K. Takatsuka, "Quantum fluctuation of electronic wave-packet dynamics coupled with classical nuclear motions," J. Chem. Phys. 122, 084113 (2005) and V. Krishna, "Path integral formulation for quantum nonadiabatic dynamics and the mixed quantum classical limit," J. Chem. Phys. 126, 134107 (2007)] involving the non-Abelian Berry force do maintain momentum conservation. As a numerical example, we investigate the Kramers doublet of the methoxy radical using generalized Hartree–Fock with spin–orbit coupling and confirm that angular momentum is conserved with the proper equations of motion. Our work makes clear some of the limitations of the Born–Oppenheimer approximation when using ab initio electronic structure theory to treat systems with unpaired electronic spin degrees of freedom, and we demonstrate that Ehrenfest dynamics can offer much improved, qualitatively correct results. [ABSTRACT FROM AUTHOR]

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems. [ABSTRACT FROM AUTHOR]

The nuclear–electronic orbital (NEO) method is a multicomponent quantum chemistry theory that describes electronic and nuclear quantum effects simultaneously while avoiding the Born–Oppenheimer approximation for certain nuclei. Typically specified hydrogen nuclei are treated quantum mechanically at the same level as the electrons, and the NEO potential energy surface depends on the classical nuclear coordinates. This approach includes nuclear quantum effects such as zero-point energy and nuclear delocalization directly into the potential energy surface. An extended NEO potential energy surface depending on the expectation values of the quantum nuclei incorporates coupling between the quantum and classical nuclei. Herein, theoretical methodology is developed to optimize and characterize stationary points on the standard or extended NEO potential energy surface, to generate the NEO minimum energy path from a transition state down to the corresponding reactant and product, and to compute thermochemical properties. For this purpose, the analytic coordinate Hessian is developed and implemented at the NEO Hartree–Fock level of theory. These NEO Hessians are used to study the SN2 reaction of ClCH3Cl− and the hydride transfer of C4H9+. For each system, analysis of the single imaginary mode at the transition state and the intrinsic reaction coordinate along the minimum energy path identifies the dominant nuclear motions driving the chemical reaction. Visualization of the electronic and protonic orbitals along the minimum energy path illustrates the coupled electronic and protonic motions beyond the Born–Oppenheimer approximation. This work provides the foundation for applying the NEO approach at various correlated levels of theory to a wide range of chemical reactions. [ABSTRACT FROM AUTHOR]