This content will become publicly available on March 28, 2025
Modern electronic structure theory is built around the Born–Oppenheimer approximation and the construction of an electronic Hamiltonian Ĥel(X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the wellknown theory of electron translation (Γ′) and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phasespace electronic Hamiltonian that depends on both nuclear position and momentum, ĤPS(X,P). While classical Born–Oppenheimer dynamics that run along the eigensurfaces of the operator Ĥel(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of ĤPS(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.
more » « less Award ID(s):
 2102402
 NSFPAR ID:
 10520189
 Publisher / Repository:
 jcp
 Date Published:
 Journal Name:
 The Journal of Chemical Physics
 Volume:
 160
 Issue:
 12
 ISSN:
 00219606
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

This paper concerns the representation of angular momentum operators in the Born–Oppenheimer theory of polyatomic molecules and the various forms of the associated conservation laws. Topics addressed include the question of whether these conservation laws are exactly equivalent or only to some order of the Born–Oppenheimer parameter κ = ( m/ M) 1/4 and what the correlation is between angular momentum quantum numbers in the various representations. These questions are addressed in both problems involving a single potential energy surface and those with multiple, strongly coupled surfaces and in both the electrostatic model and those for which fine structure and electron spin are important. The analysis leads to an examination of the transformation laws under rotations of the electronic Hamiltonian; of the basis states, both adiabatic and diabatic, along with their phase conventions; of the potential energy matrix; and of the derivative couplings. These transformation laws are placed in the geometrical context of the structures in the nuclear configuration space that are induced by rotations, which include the rotational orbits or fibers, the surfaces upon which the orientation of the molecule changes but not its shape, and the section, an initial value surface that cuts transversally through the fibers. Finally, it is suggested that the usual Born–Oppenheimer approximation can be replaced by a dressing transformation, that is, a sequence of unitary transformations that blockdiagonalize the Hamiltonian. When the dressing transformation is carried out, we find that the angular momentum operator does not change. This is a part of a system of exact equivalences among various representations of angular momentum operators in Born–Oppenheimer theory. Our analysis accommodates largeamplitude motions and is not dependent on smallamplitude expansions about an equilibrium position. Our analysis applies to noncollinear configurations of a polyatomic molecule; this covers all but a subset of measure zero (the collinear configurations) in the nuclear configuration space.more » « less

Total angular momentum conservation in Ehrenfest dynamics with a truncated basis of adiabatic states
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 wavepacket 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 nonAbelian 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.

We demonstrate that, for systems with spin–orbit coupling and an odd number of electrons, the standard fewest switches surface hopping algorithm does not conserve the total linear or angular momentum. This lack of conservation arises not so much from the hopping direction (which is easily adjusted) but more generally from propagating adiabatic dynamics along surfaces that are not time reversible. We show that one solution to this problem is to run along eigenvalues of phasespace electronic Hamiltonians H(R, P) (i.e., electronic Hamiltonians that depend on both nuclear position and momentum) with an electronic–nuclear coupling Γ · P [see Eq. (25)], and we delineate the conditions that must be satisfied by the operator Γ. The present results should be extremely useful as far as developing new semiclassical approaches that can treat systems where the nuclear, electronic orbital, and electronic spin degrees of freedom altogether are all coupled together, hopefully including systems displaying the chiralinduced spin selectivity effect.

Within a tightbinding approximation, we numerically determine the time evolution of graphene electronic states in the presence of classically vibrating nuclei. There is no reliance on the Born–Oppenheimer approximation within the porbital tightbinding basis, although our approximation is “atomically adiabatic”: the basis porbitals are taken to follow nuclear positions. Our calculations show that the strict adiabatic Born–Oppenheimer approximation fails badly. We find that a diabatic (lazy electrons responding weakly to nuclear distortions) Born–Oppenheimer model provides a much more accurate picture and suggests a generalized manybody Bloch orbitalnuclear basis set for describing electron–phonon interactions in graphene.more » « less

Within the nuclearelectronic orbital (NEO) framework, the realtime NEO timedependent density functional theory (RTNEOTDDFT) approach enables the simulation of coupled electronicnuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of longtime nuclear quantum dynamics. Herein, the electronic Born–Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the realtime nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an orderofmagnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RTNEOTDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RTNEOEhrenfest dynamics and its BO counterpart can describe proton delocalization during the realtime nuclear quantum dynamics. Thus, the BO RTNEO approach provides the foundation for a wide range of chemical and biological applications.