Electronic currents and Born-Oppenheimer molecular dynamics

A basis-free phase space electronic Hamiltonian that recovers beyond Born-Oppenheimer electronic momentum and current density

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).

An electronic phase-space Hamiltonian approach for electronic current density and vibrational circular dichroism

The Born-Oppenheimer framework stipulates that chemistry and physics occur on potential energy surfaces VBO(X) parameterized by a nuclear coordinate X, which are built by diagonalizing a BO Hamiltonian ĤBO(X). However, such a framework cannot recover many measurable chemical and physical features, including vibrational circular dichroism spectra. In this article, we show that a phase-space electronic Hamiltonian ĤPS(X,P), parameterized by both nuclear position X and momentum P, with a similar computational cost as solving ĤBO(X), can recover not just experimental vibrational circular dichroism signals but also a meaningful electronic current density that explains the features of the vibrational circular dichroism rotational strengths. Combined with earlier demonstrations that such Hamiltonians can also recover qualitatively correct electronic momenta with electronic densities that approximately satisfy a continuity equation, the data would suggest that, if one looks closely enough, chemistry in fact occurs on potential energy surfaces parameterized by both X and P, EPS(X, P). While the dynamical implications of such a phase-space electronic Hamiltonian are not yet known, we hypothesize that, by offering classical trajectories that explicitly offer nonzero electronic momentum while also conserving the total angular momentum (unlike Born-Oppenheimer theory), this new phase-space electronic structure Hamiltonian may well explain some fraction of the chiral-induced spin selectivity effect.

A simple one-electron expression for electron rotational factors

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.

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.

Diagonalizing the Born-Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom

This article describes a method for calculating higher order or nonadiabatic corrections in Born-Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner-Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ2, where κ = (m/M)1/4 is the usual Born-Oppenheimer parameter. The lowest order term is the usual Born-Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born-Oppenheimer theory. We show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born-Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.

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 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.

Ultrafast Molecular Frame Quantum Tomography

We develop and experimentally demonstrate a methodology for a full molecular frame quantum tomography (MFQT) of dynamical polyatomic systems. We exemplify this approach through the complete characterization of an electronically nonadiabatic wave packet in ammonia (NH_{3}). The method exploits both energy and time-domain spectroscopic data, and yields the lab frame density matrix (LFDM) for the system, the elements of which are populations and coherences. The LFDM fully characterizes electronic and nuclear dynamics in the molecular frame, yielding the time- and orientation-angle dependent expectation values of any relevant operator. For example, the time-dependent molecular frame electronic probability density may be constructed, yielding information on electronic dynamics in the molecular frame. In NH_{3}, we observe that electronic coherences are induced by nuclear dynamics which nonadiabatically drive electronic motions (charge migration) in the molecular frame. Here, the nuclear dynamics are rotational and it is nonadiabatic Coriolis coupling which drives the coherences. Interestingly, the nuclear-driven electronic coherence is preserved over longer timescales. In general, MFQT can help quantify entanglement between electronic and nuclear degrees of freedom, and provide new routes to the study of ultrafast molecular dynamics, charge migration, quantum information processing, and optimal control schemes.

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.

Adiabatic electronic flux in molecules and in condensed matter

The theory of adiabatic electron transport in a correlated condensed-matter system is rooted in a seminal paper by Niu and Thouless [J. Phys A {\bf 17}, 2453 (1984)]; I adopt here an analogous logic in order to { retrieve the known expression for the adiabatic electronic flux in a molecular system [L. A. Nafie, J. Chem. Phys. {\bf 79}, 4950 (1983)]. Its derivation here is considerably simpler than those available in the current quantum-chemistry literature; it also explicitly identifies the adiabaticity parameter, in terms of which the adiabatic flux and the electron density are both exact to first order. It is shown that continuity equation is conserved to the same order. For the sake of completeness,} I also briefly outline the relevance of the macroscopic electronic flux to the physics of solids and liquids.

On the molecular electronic flux: Role of nonadiabaticity and violation of conservation

Analysis of electron flux within and in between molecules is crucial in the study of real-time dynamics of molecular electron wavepacket evolution such as those in attosecond laser chemistry and ultrafast chemical reaction dynamics. We here address two mutually correlated issues on the conservation law of molecular electronic flux, which serves as a key consistency condition for electron dynamics. The first one is about a close relation between "weak" nonadiabaticity and the electron dynamics in low-energy chemical reactions. We show that the electronic flux in adiabatic reactions can be consistently reproduced by taking account of nonadiabaticity. Such nonadiabaticity is usually weak in the sense that it does not have a major effect on nuclear dynamics, whereas it plays an important role in electronic dynamics. Our discussion is based on a nonadiabatic extension of the electronic wavefunction similar in idea to the complete adiabatic formalism developed by Nafie [J. Chem. Phys. 79, 4950 (1983)], which has also recently been reformulated by Patchkovskii [J. Chem. Phys. 137, 084109 (2012)]. We give straightforward proof of the theoretical assertion presented by Nafie using a time-dependent mixed quantum-classical framework and a standard perturbation expansion. Explicitly taking account of the flux conservation, we show that the nonadiabatically induced flux realizes the adiabatic time evolution of the electronic density. In other words, the divergence of the nonadiabatic flux equals the time derivative of the electronic density along an adiabatic time evolution of the target molecule. The second issue is about the accurate computationability of the flux. The calculation of flux needs an accurate representation of the (relative) quantum phase, in addition to the amplitude factor, of a total wavefunction and demands special attention for practical calculations. This paper is the first one to approach this issue directly and show how the difficulties arise explicitly. In doing so, we reveal that a number of widely accepted truncation techniques for static property calculations are potential sources of numerical flux non-conservation. We also theoretically propose alternative strategies to realize better flux conservation.

Probing Electronic Fluxes via Time-Resolved X-Ray Scattering

The current flux density is a vector field that can be used to describe theoretically how electrons flow in a system out of equilibrium. In this work, we unequivocally demonstrate that the signal obtained from time-resolved x-ray scattering does not only map the time evolution of the electronic charge distribution, but also encodes information about the associated electronic current flux density. We show how the electronic current flux density qualitatively maps the distribution of electronic momenta and reveals the underlying mechanism of ultrafast charge migration processes, while also providing quantitative information about the timescales of electronic coherences.

Concerted Electronic and Nuclear Fluxes During Coherent Tunnelling in Asymmetric Double-Well Potentials

The quantum theory of concerted electronic and nuclear fluxes (CENFs) during coherent periodic tunnelling from reactants (R) to products (P) and back to R in molecules with asymmetric double-well potentials is developed. The results are deduced from the solution of the time-dependent Schrödinger equation as a coherent superposition of two eigenstates; here, these are the two states of the lowest tunnelling doublet. This allows the periodic time evolutions of the resulting electronic and nuclear probability densities (EPDs and NPDs) as well as the CENFs to be expressed in terms of simple sinusodial functions. These analytical results reveal various phenomena during coherent tunnelling in asymmetric double-well potentials, e.g., all EPDs and NPDs as well as all CENFs are synchronous. Distortion of the symmetric reference to a system with an asymmetric double-well potential breaks the spatial symmetry of the EPDs and NPDs, but, surprisingly, the symmetry of the CENFs is conserved. Exemplary application to the Cope rearrangement of semibullvalene shows that tunnelling of the ideal symmetric system can be suppressed by asymmetries induced by rather small external electric fields. The amplitude for the half tunnelling, half nontunnelling border is as low as 0.218 × 10(-8) V/cm. At the same time, the delocalized eigenstates of the symmetric reference, which can be regarded as Schrödinger's cat-type states representing R and P with equal probabilities, get localized at one or the other minima of the asymmetric double-well potential, representing either R or P.

Staircase patterns of nuclear fluxes during coherent tunneling in excited doublets of symmetric double well potentials

Step-by-step flux for one-by-one transfers of the lobes of the density, from the reactant (left) to the product (right) in the excited tunneling doublet.

Quantum theory of concerted electronic and nuclear fluxes associated with adiabatic intramolecular processes

Example of concerted electronic (right) and nuclear (left) fluxes: isomerization of B₄.