Search

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.

The classical equations of motion of a charged particle in a spherically symmetric distribution of magnetic monopoles can be transformed into a system of linear equations, thereby providing a type of integrability. In the case of a single monopole, the solution was given long ago by Poincar\'e. In the case of a uniform distribution of monopoles the solution can be expressed in terms of parabolic cylinder functions (essentially the eigenfunctions of an inverted harmonic oscillator). This solution is relevant to recent studies of nonassociative star products, symplectic lifts of twisted Poisson structures and fluids and plasmas of electric and magnetic charges., Comment: 20 pages 1 figure

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.

The Wigner $3j$ symbols of the quantum angular momentum theory are related to the vector coupling or Clebsch-Gordan coefficients and to the Hahn and dual Hahn polynomials of the discrete orthogonal hyperspherical family, of use in discretization approximations. We point out the important role of the Regge symmetries for defining the screen where images of the coefficients are projected, and for discussing their asymptotic properties and semiclassical behavior. Recursion relationships are formulated as eigenvalue equations, and exploited both for computational purposes and for physical interpretations., Comment: 14 pages, 8 figures, presented at ICCSA 2014, 14th International Conference on Computational Science and Applications

The Schl\"afli identity, which is important in Regge calculus and loop quantum gravity, is examined from a symplectic and semiclassical standpoint in the special case of flat, 3-dimensional space. In this case a proof is given, based on symplectic geometry. A series of symplectic and Lagrangian manifolds related to the Schl\"afli identity, including several versions of a Lagrangian manifold of tetrahedra, are discussed. Semiclassical interpretations of the various steps are provided. Possible generalizations to 3-dimensional spaces of constant (nonzero) curvature, involving Poisson-Lie groups and q-deformed spin networks, are discussed., Comment: 40 pages, 8 figures

This article presents and discusses in detail the results of extensive exact calculations of the most basic ingredients of spin networks, the Racah coefficients (or Wigner 6j symbols), exhibiting their salient features when considered as a function of two variables - a natural choice due to their origin as elements of a square orthogonal matrix - and illustrated by use of a projection on a square "screen" introduced recently. On these screens, shown are images which provide a systematic classification of features previously introduced to represent the caustic and ridge curves (which delimit the boundaries between oscillatory and evanescent behaviour according to the asymptotic analysis of semiclassical approaches). Particular relevance is given to the surprising role of the intriguing symmetries discovered long ago by Regge and recently revisited; from their use, together with other newly discovered properties and in conjunction with the traditional combinatorial ones, a picture emerges of the amplitudes and phases of these discrete wavefunctions, of interest in wide areas as building blocks of basic and applied quantum mechanics., Comment: 16 pages, 13 figures, presented at ICCSA 2013 13th International Conference on Computational Science and Application

Maslov indices are integers that appear in semiclassical wave functions and quantization conditions. They are often notoriously difficult to compute. We present methods of computing the Maslov index that rely only on typically elementary Poisson brackets and simple linear algebra. We also present a singular differential form, whose integral along a curve gives the Maslov index of that curve. The form is closed but not exact, and transforms by an exact differential under canonical transformations. We illustrate the method with the $6j$-symbol, which is important in angular momentum theory and in quantum gravity., Comment: 5 pages

We derive a new asymptotic formula for the Wigner $9j$-symbol, in the limit of one small and eight large angular momenta, using a novel gauge-invariant factorization for the asymptotic solution of a set of coupled wave equations. Our factorization eliminates the geometric phases completely, using gauge-invariant non-canonical coordinates, parallel transports of spinors, and quantum rotation matrices. Our derivation generalizes to higher $3nj$-symbols. We display without proof some new asymptotic formulas for the $12j$-symbol and the $15j$-symbol in the appendices. This work contributes a new asymptotic formula of the Wigner $9j$-symbol to the quantum theory of angular momentum, and serves as an example of a new general method for deriving asymptotic formulas for $3nj$-symbols., Comment: 18 pages, 16 figures. To appear in Phys. Rev. A

The semiclassical mechanics of the Wigner 6j-symbol is examined from the standpoint of WKB theory for multidimensional, integrable systems, to explore the geometrical issues surrounding the Ponzano-Regge formula. The relations among the methods of Roberts and others for deriving the Ponzano-Regge formula are discussed, and a new approach, based on the recoupling of four angular momenta, is presented. A generalization of the Yutsis-type of spin network is developed for this purpose. Special attention is devoted to symplectic reduction, the reduced phase space of the 6j-symbol (the 2-sphere of Kapovich and Millson), and the reduction of Poisson bracket expressions for semiclassical amplitudes. General principles for the semiclassical study of arbitrary spin networks are laid down; some of these were used in our recent derivation of the asymptotic formula for the Wigner 9j-symbol., Comment: 64 pages, 50 figures

We present the asymptotic formula for the Wigner 9j-symbol, valid when all quantum numbers are large, in the classically allowed region. As in the Ponzano-Regge formula for the 6j-symbol, the action is expressed in terms of lengths of edges and dihedral angles of a geometrical figure, but the angles require care in definition. Rules are presented for converting spin networks into the associated geometrical figures. The amplitude is expressed as the determinant of a 2x2 matrix of Poisson brackets. The 9j-symbol possesses caustics associated with the fold and elliptic and hyperbolic umbilic catastrophes. The asymptotic formula obeys the exact symmetries of the 9j-symbol., Comment: 17 pages, 7 figures

A new uniform asymptotic approximation for the Wigner $6j$ symbol is given in terms of Wigner rotation matrices ($d$-matrices). The approximation is uniform in the sense that it applies for all values of the quantum numbers, even those near caustics. The derivation of the new approximation is not given, but the geometrical ideas supporting it are discussed and numerical tests are presented, including comparisons with the exact $6j$-symbol and with the Ponzano-Regge approximation., Comment: 44 pages plus 20 figures

We analyze the asymptotics of the Wigner $3j$-symbol as a matrix element connecting eigenfunctions of a pair of integrable systems, obtained by lifting the problem of the addition of angular momenta into the space of Schwinger's oscillators. A novel element is the appearance of compact Lagrangian manifolds that are not tori, due to the fact that the observables defining the quantum states are noncommuting. These manifolds can be quantized by generalized Bohr-Sommerfeld rules and yield all the correct quantum numbers. The geometry of the classical angular momentum vectors emerges in a clear manner. Efficient methods for computing amplitude determinants in terms of Poisson brackets are developed and illustrated., Comment: 7 figure files

The internal space for a molecule, atom, or other n-body system can be conveniently parameterised by 3n-9 kinematic angles and three kinematic invariants. For a fixed set of kinematic invariants, the kinematic angles parameterise a subspace, called a kinematic orbit, of the n-body internal space. Building on an earlier analysis of the three- and four-body problems, we derive the form of these kinematic orbits (that is, their topology) for the general n-body problem. The case n=5 is studied in detail, along with the previously studied cases n=3,4., Comment: 38 pages, submitted to J. Phys. A

For a three-body system, a quantum wave function $\Psi^\ell_m$ with definite $\ell$ and $m$ quantum numbers may be expressed in terms of an internal wave function $\chi^\ell_k$ which is a function of three internal coordinates. This article provides necessary and sufficient constraints on $\chi^\ell_k$ to ensure that the external wave function $\Psi^\ell_m$ is analytic. These constraints effectively amount to boundary conditions on $\chi^\ell_k$ and its derivatives at the boundary of the internal space. Such conditions find similarities in the (planar) two-body problem where the wave function (to lowest order) has the form $r^{|m|}$ at the origin. We expect the boundary conditions to prove useful for constructing singularity free three-body basis sets for the case of nonvanishing angular momentum., Comment: 41 pages, submitted to Phys. Rev. A

The rovibrational kinetic energy for an arbitrary number of rigid molecules is computed. The result has the same general form as the kinetic energy in the molecular rovibrational Hamiltonian, although certain quantities are augmented to account for the rotational energy of the monomers. No specific choices of internal coordinates or body frame are made in order to accommodate the large variety of such conventions. However, special attention is paid to how key quantities transform when these conventions are changed. An example system is explicitly analysed as an illustration of the formalism., Comment: 20 pages, submitted to Molecular Physics