An official website of the United States government
We present a phase-space electronic Hamiltonian ĤPS (parameterized by both nuclear position X and momentum P) that boosts each electron into the moving frame of the nuclei that are closest in real space. The final form for the phase space Hamiltonian does not assume the existence of an atomic orbital basis, and relative to standard Born–Oppenheimer theory, the newly proposed one-electron operators can be expressed directly as functions of electronic and nuclear positions and momentum. We show that (i) quantum–classical dynamics along such a Hamiltonian maintains momentum conservation and that (ii) diagonalizing such a Hamiltonian can recover the electronic momentum and electronic current density reasonably well. In conjunction with other reports in the literature that such a phase-space approach can also recover vibrational circular dichroism spectra, we submit that the present phase-space approach offers a testable and powerful approach to post-BO electronic structure theory. Moreover, the approach is inexpensive and can be immediately applied to simulations of chiral induced spin selectivity experiments (where the transfer of angular momentum between nuclei and electrons is considered critical).
more »
« less
Practical phase-space electronic Hamiltonians for ab initio dynamics
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 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, Ĥ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
- PAR ID:
- 10520189
- Publisher / Repository:
- jcp
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 160
- Issue:
- 12
- ISSN:
- 0021-9606
- 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 block-diagonalize 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 large-amplitude motions and is not dependent on small-amplitude 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 statesWe 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.more » « less
-
We present an ab initio method for computing vibro-polariton and phonon-polariton spectra of molecules and solids coupled to the photon modes of optical cavities. We demonstrate that if interactions of cavity photon modes with both nuclear and electronic degrees of freedom are treated on the level of the cavity Born–Oppenheimer approximation, spectra can be expressed in terms of the matter response to electric fields and nuclear displacements, which are readily available in standard density functional perturbation theory implementations. In this framework, results over a range of cavity parameters can be obtained without the need for additional electronic structure calculations, enabling efficient calculations on a wide range of parameters. Furthermore, this approach enables results to be more readily interpreted in terms of the more familiar cavity-independent molecular electric field response properties, such as polarizability and Born effective charges, which enter into the vibro-polariton calculation. Using corresponding electric field response properties of bulk insulating systems, we are also able to obtain the Γ point phonon-polariton spectra of two dimensional (2D) insulators. Results for a selection of cavity-coupled molecular and 2D crystal systems are presented to demonstrate the method.more » « less
-
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 phase-space 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 chiral-induced spin selectivity effect.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0192084
