Search

According to a result of Wigner and von Neumann [1], real symmetric matrices with a doubly degenerate lowest eigenvalue form a submanifold of codimension 2 within the space of all real symmetric matrices. This mathematical result has important consequences for chemistry. First, it implies that degeneracies do not occur within generic one-parameter families of real symmetric matrices - this is the famous non-crossing rule, and is responsible for the phenomenon of avoided crossings in the energy levels of diatomic molecules. Second, it implies that energy levels are expected to cross in polyatomic molecules, with crossings taking place on a submanifold of nuclear configuration space which is codimension 2 - this submanifold is the famous conical intersection seam, of central importance in nonadiabatic chemistry. In this paper we extend the analysis of Wigner and von Neumann to include symmetry. We introduce a symmetry group, and consider parametrised families of matrices which respect an action of that symmetry group on both the parameter space and on the space of matrices. A concrete application is to molecules, for which the relevant symmetry group is generated by permutations of atomic nuclei combined with spatial reflections and rotations. In the presence of this extra symmetry, we find that energy level crossings do not typically occur on codimension 2 submanifolds, and connect our findings with the discovery of confluences of conical intersection seams in the chemical literature. We give a classification of confluences for triatomic molecules and planar molecules, unifying the previous literature on this topic, and predict several new types of confluence.

We show that following the standard mantra of quantum chemistry and diagonalizing the Born-Oppenheimer (BO) Hamiltonian $\hat H_{\rm BO}(\bm R)$ is not the optimal means to construct potential energy surfaces. A better approach is to diagonalize a phase-space electronic Hamiltonian, $\hat H_{\rm PS}(\bm R,\bm P)$, which is parameterized by both nuclear position $\bm R$ and nuclear momentum $\bm P$. The foundation of such a non-perturbative phase-space electronic Hamiltonian can be made rigorous using a partial Wigner transform and the method has exactly the same cost as BO for a semiclassical calculation (and only a slight increase in cost for a quantum nuclear calculation). For a three-particle system, with two heavy particles and one light particle, numerical results show that a phase space electronic Hamiltonian produces not only meaningful electronic momenta (which are completely ignored by BO theory) but also far better vibrational energies. As such, for high level results and/or systems with degeneracies and spin degrees of freedom, we anticipate that future electronic structure and quantum chemistry packages will need to take as input not just the positions of the nuclei but also their momenta.

We show that a phase space electronic Hamiltonian $\hat{H}_{PS}(\mathbf{X},\mathbf{P})$, parameterized by both nuclear position $\mathbf{X}$ and momentum $\mathbf{P}$, can recover not just experimental vibrational circular dichroism (VCD) signals, but also a meaningful electronic current density that explains the features of the VCD rotatory 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 we have isolated a meaningful alternative approach to electronic structure theory, one that entirely avoids Born-Oppenheimer theory and frozen nuclei. While the dynamical implications of such a phase space electronic Hamiltonian are not yet known, we hypothesize that, by offering classical trajectories the conserve 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.

We show empirically that a phase-space non-Born-Oppenheimer electronic Hamiltonian approach to quantum chemistry (where the electronic Hamiltonian is parameterized by both nuclear position and momentum, (H(R,P)) is both a practical and accurate means to recover vibrational circular dichroism spectra. We further hypothesize that such a phase space approach may lead to very new dynamical physics beyond spectroscopy circular dichroism, with potential implications for understanding chiral induced spin selectivity (CISS), noting that classical phase space approaches conserve the total nuclear plus electronic momentum, whereas classical Born-Oppenheimer approaches do not (they conserve only the nuclear momentum), Comment: 46 pages, 5 figures

Modern electronic structure theory is built around the Born-Oppenheimer approximation and the construction of an electronic Hamiltonian H_{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 (Gamma') and rotational (Gamma'') 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_{PS}(X,P). While classical Born-Oppenheimer dynamics that run along the eigensurfaces of the operator H_{el}(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of H_{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.

Within the context of FSSH dynamics, one often wishes to remove the angular component of the derivative coupling between states $\left|J\right>$ and $\left|K\right>$. In a set of previous papers, Truhlar {\em et al.} posited one approach for such a removal based on direct projection, while we isolated a second approach by constructing and differentiating rotationally invariant basis. Unfortunately, neither approach was able to demonstrate a {\em one-electron operator} $\hat{O}$ whose matrix element $\left$ 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 be of immediate use for FSSH calculations.

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. The result is a power series in $\kappa^2$, where $\kappa =(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 method was applied to Born-Oppenheimer theory by Littlejohn and Weigert (1993), in a treatment that notably produced the correction $K_{22}$ to the Born-Oppenheimer Hamiltonian (see {\em infra}). Recently Matyus and Teufel (2019) have applied an improved and more elegant version of the method to Born-Oppenheimer theory, and have calculated the Born-Oppenheimer Hamiltonian for multiple potential energy surfaces to order $\kappa^6$. One of the shortcomings of earlier methods is that 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 need to be a part of the discussion. This article presents an independent derivation of the Moyal expansion in molecular Born-Oppenheimer theory, with special attention to the translational degrees of freedom. 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., Comment: 88 pages; 1 figure

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[Amano2005,Krishna2007] 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 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.

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 [ETFs]); 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 [ERFs] 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 new surface hopping algorithm in the long term for systems where one cannot ignore the electronic orbital and/or spin angular momentum.

We prove both analytically and numerically that the total angular momentum of a molecular system undergoing adiabatic Born-Oppenheimer dynamics is conserved only when pseudo-magnetic Berry forces are taken into account. This finding sheds light on the nature of Berry forces for molecular systems with spin-orbit coupling and highlights how ab initio Born-Oppenheimer molecular dynamics simulations can successfully capture the entanglement of spin and nuclear degrees of freedom as modulated by electronic interactions.

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 $\kappa=(m/M)^{1/4}$, and what the correlation is between angular momentum quantum numbers in the various representations. These questions are addressed both in problems involving a single potential energy surface, and those with multiple, strongly coupled surfaces; and both in 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..., Comment: 75 pages, 4 figures

This article concerns the use of parallel transport to create a diabatic basis. The advantages of the parallel-transported basis include the facility with which Taylor series expansions can be carried out in the neighborhood of a point or a manifold such as a seam (the locus of degeneracies of the electronic Hamiltonian), and the close relationship between the derivative couplings and the curvature in this basis. These are important for analytic treatments of the nuclear Schr\"odinger equation in a neighborhood of degeneracies. The parallel-transported basis bears a close relationship to the singular-value basis; in this article both are expanded in power series about a reference point and they are shown to agree through second order but not beyond. Taylor series expansions are effected through the projection operator, whose expansion does not involve energy denominators or any type of singularity, and in terms of which both the singular-value basis and the parallel-transported basis can be expressed. The parallel-transported basis is a version of Poincar\'e gauge, well known in electromagnetism, which provides a relationship between the derivative couplings and the curvature and which, along with a formula due to Mead, affords an efficient method for calculating Taylor series of the basis states and the derivative couplings. The case in which fine structure effects are included in the electronic Hamiltonian is covered., Comment: 44 pages, 2 figures

Nuclear Berry curvature effects emerge from electronic spin degeneracy and canlead to non-trivial spin-dependent (nonadiabatic) nuclear dynamics. However, such effects are completely neglected in all current mixed quantum-classical methods such as fewest switches surface-hopping. In this work, we present a phase-space surface-hopping (PSSH) approach to simulate singlet-triplet intersystem crossing dynamics. We show that with a simple pseudo-diabatic ansatz, a PSSH algorithm can capture the relevant Berry curvature effects and make predictions in agreement with exact quantum dynamics for a simple singlet-triplet model Hamiltonian. Thus, this approach represents an important step towards simulating photochemical and spin processes concomitantly, as relevant to intersystem crossing and spin-lattice relaxation dynamics.

Chemical relaxation phenomena, including photochemistry and electron transfer processes, form a vigorous area of research in which nonadiabatic dynamics plays a fundamental role. Here, we show that for nonadiabatic dynamics with two electronic states and a complex-valued Hamiltonian that does not obey time-reversal symmetry, the optimal semiclassical approach is to run surface hopping dynamics on a set of phase-space adiabatic surfaces. In order to generate such phase-adiabats, one must isolate a proper set of diabats and apply a phase gauge transformation, before eventually diagonalizing the total Hamiltonian (which is now parameterized by both R and P). The resulting algorithm is valid in both the adiabatic and nonadiabatic limits, incorporates all Berry curvature effects, and allows for the study of semiclassical nonadiabatic dynamics in the presence of spin-orbit coupling and/or external magnetic fields.

Another manifestation of the Aharonov-Bohm effect is introduced to chemistry, in fact to nuclear dynamics and high-resolution molecular spectroscopy. As demonstrated, the overall rotation of a symmetric-top molecule influences the dynamics of an internal vibrational motion in a way that is analogous to the presence of a solenoid carrying magnetic flux. To a good approximation, the low-energy rovibrational energy-level structure of the quasistructural molecular ion H5+ can be understood entirely in terms of this effect.

Three categories of algebras with morphisms generalising the usual set of algebra homomorphisms are described. The Sweedler product provides a hom-tensor equivalence relating these three categories, and a tool enabling the universal measuring coalgebra to be calculated in small cases. A parallel theory for modules is presented.

A new one-dimensional model is proposed for the low-energy vibrational quantum dynamics of CH5+ based on the motion of an effective particle confined to a 60-vertex graph ${\Gamma}_{60}$ with a single edge length parameter. Within this model, the quantum states of CH5+ are obtained in analytic form and are related to combinatorial properties of ${\Gamma}_{60}$. The bipartite structure of ${\Gamma}_{60}$ gives a simple explanation for curious symmetries observed in numerically exact variational calculations on CH5+.

Starting with the exact factorization of the molecular wavefunction, this paper presents the results from the numerical implementation in nonadiabatic molecular dynamics of the recently proposed bohmion method. Within the context of quantum hydrodynamics, we introduce a regularized nuclear Bohm potential admitting solutions comprising a train of $\delta$-functions which provide a finite-dimensional sampling of the hydrodynamic flow paths. The bohmion method inherits all the basic conservation laws from its underlying variational structure and captures electronic decoherence. After reviewing the general theory, the method is applied to the well-known Tully models, which are used here as benchmark problems. In the present case of study, we show that the new method accurately reproduces both electronic decoherence and nuclear population dynamics.

While the treatment of conical intersections in molecular dynamics generally requires nonadiabatic approaches, the Born-Oppenheimer adiabatic approximation is still adopted as a valid alternative in certain circumstances. In the context of Mead-Truhlar minimal coupling, this paper presents a new closure of the nuclear Born-Oppenheimer equation, thereby leading to a molecular dynamics scheme capturing geometric phase effects. Specifically, a semiclassical closure of the nuclear Ehrenfest dynamics is obtained through a convenient prescription for the nuclear Bohmian trajectories. The conical intersections are suitably regularized in the resulting nuclear particle motion and the associated Lorentz force involves a smoothened Berry curvature identifying a loop-dependent geometric phase. In turn, this geometric phase rapidly reaches the usual topological index as the loop expands away from the original singularity. This feature reproduces the phenomenology appearing in recent exact nonadiabatic studies, as shown explicitly in the Jahn-Teller problem for linear vibronic coupling. Likewise, a newly proposed regularization of the diagonal correction term is also shown to reproduce quite faithfully the energy surface presented in recent nonadiabatic studies., Comment: Third version with minor changes. To appear in Phys. Rev. A