An official website of the United States government
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
Surface hopping, electron translation factors, electron rotation factors, momentum conservation, and size consistency
For a system without spin–orbit coupling, the (i) nuclear plus electronic linear momentum and (ii) nuclear plus orbital electronic angular momentum are good quantum numbers. Thus, when a molecular system undergoes a nonadiabatic transition, there should be no change in the total linear or angular momentum. Now, the standard surface hopping algorithm ignores the electronic momentum and indirectly equates the momentum of the nuclear degrees of freedom to the total momentum. However, even with this simplification, the algorithm still does not conserve either the nuclear linear or the nuclear angular momenta. Here, we show that one way to address these failures is to dress the derivative couplings (i.e., the hopping directions) in two ways: (i) we disallow changes in the nuclear linear momentum by working in a translating basis (which is well known and leads to electron translation factors) and (ii) we disallow changes in the nuclear angular momentum by working in a basis that rotates around the center of mass [which is not well-known and leads to a novel, rotationally removable component of the derivative coupling that we will call electron rotation factors below, cf. Eq. (96)]. The present findings should be helpful in the short term as far as interpreting surface hopping calculations for singlet systems (without spin) and then developing the new surface hopping algorithm in the long term for systems where one cannot ignore the electronic orbital and/or spin angular momentum.
more »
« less
- Award ID(s):
- 2102402
- PAR ID:
- 10519765
- Publisher / Repository:
- J Chem Pays
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 159
- Issue:
- 11
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Within the context of fewest-switch surface hopping (FSSH) dynamics, one often wishes to remove the angular component of the derivative coupling between states J and K. In a previous set of papers, Shu et al. [J. Phys. Chem. Lett. 11, 1135–1140 (2020)] posited one approach for such a removal based on direct projection, while we isolated a second approach by constructing and differentiating a rotationally invariant basis. Unfortunately, neither approach was able to demonstrate a one-electron operatorÔ whose matrix element JÔK was the angular component of the derivative coupling. Here, we show that a one-electron operator can, in fact, be constructed efficiently in a semi-local fashion. The present results yield physical insight into designing new surface hopping algorithms and are of immediate use for FSSH calculations.more » « less
-
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
-
We demonstrate that working with a correct phase-space electronic Hamiltonian captures electronic inertial effects. In particular, we show that phase space surface hopping dynamics do not suffer (at least to very high order) from non-physical non-adiabatic transitions between electronic eigenstates during the course of pure nuclear translational and rotational motion. This work opens up many new avenues for quantitatively investigating complex phenomena, including angular momentum transfer between chiral phonons and electrons as well as chiral-induced spin selectivity effects.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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0160965
