skip to main content


Title: Representation and conservation of angular momentum in the Born–Oppenheimer theory of polyatomic molecules
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
Award ID(s):
2102402
NSF-PAR ID:
10418876
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
158
Issue:
10
ISSN:
0021-9606
Page Range / eLocation ID:
104302
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate the Polariton induced conical intersection (PICI) created from coupling a diatomic molecule with the quantized photon mode inside an optical cavity, and the corresponding Berry Phase effects. We use the rigorous Pauli–Fierz Hamiltonian to describe the quantum light-matter interactions between a LiF molecule and the cavity, and use the exact quantum propagation to investigate the polariton quantum dynamics. The molecular rotations relative to the cavity polarization direction play a role as the tuning mode of the PICI, resulting in an effective CI even within a diatomic molecule. To clearly demonstrate the dynamical effects of the Berry phase, we construct two additional models that have the same Born–Oppenheimer surface, but the effects of the geometric phase are removed. We find that when the initial wavefunction is placed in the lower polaritonic surface, the Berry phase causes a π phase-shift in the wavefunction after the encirclement around the CI, indicated from the nuclear probability distribution. On the other hand, when the initial wavefunction is placed in the upper polaritonic surface, the geometric phase significantly influences the couplings between polaritonic states and therefore, the population dynamics between them. These BP effects are further demonstrated through the photo-fragment angular distribution. PICI created from the quantized radiation field has the promise to open up new possibilities to modulate photochemical reactivities. 
    more » « less
  2. null (Ed.)
    The exact energy and angular momentum conservation laws are derived by the Noether method for the Hamiltonian and symplectic representations of the gauge-free electromagnetic gyrokinetic Vlasov–Maxwell equations. These gyrokinetic equations, which are solely expressed in terms of electromagnetic fields, describe the low-frequency turbulent fluctuations that perturb a time-independent toroidally-axisymmetric magnetized plasma. The explicit proofs presented here provide a complete picture of the transfer of energy and angular momentum between the gyrocentres and the perturbed electromagnetic fields, in which the crucial roles played by gyrocentre polarization and magnetization effects are highlighted. In addition to yielding an exact angular momentum conservation law, the gyrokinetic Noether equation yields an exact momentum transport equation, which might be useful in more general equilibrium magnetic geometries. 
    more » « less
  3. The effects of hyperfine structure on ultracold molecular collisions in external fields are largely unexplored due to major computational challenges associated with rapidly proliferating hyperfine and rotational channels coupled by highly anisotropic intermolecular interactions. We explore a new basis set for incorporating the effects of hyperfine structure and external magnetic fields in quantum scattering calculations on ultracold molecular collisions. The basis is composed of direct products of the eigenfunctions of the total {\it rotational} angular momentum (TRAM) of the collision complex Jr and the electron/nuclear spin basis functions of the collision partners. The separation of the rotational and spin degrees of freedom ensures rigorous conservation of Jr even in the presence of external magnetic fields and isotropic hyperfine interactions. The resulting block-diagonal structure of the scattering Hamiltonian enables coupled-channel calculations on highly anisotropic atom-molecule and molecule-molecule collisions to be performed independently for each value of Jr, with an added advantage of eliminating the unphysical states present in the total angular momentum representation. We illustrate the efficiency of the TRAM basis by calculating state-to-state cross sections for ultracold He + YbF collisions in a magnetic field. The size of the TRAM basis required to reach numerical convergence is 8 times smaller than that of the uncoupled basis used previously, providing a computational gain of three orders of magnitude. The TRAM basis is therefore well suited for rigorous quantum scattering calculations on ultracold molecular collisions in the presence of hyperfine interactions and external magnetic fields. 
    more » « less
  4. null (Ed.)
    Abstract The convergence characteristics of three geometrically accurate spatial finite elements (FEs) are examined in this study using an eigenvalue analysis. The spatial beam, plate, and solid elements considered in this investigation are suited for both structural and multibody system (MBS) applications. These spatial elements are based on geometry derived from the kinematic description of the absolute nodal coordinate formulation (ANCF). In order to allow for an accurate reference-configuration geometry description, the element shape functions are formulated using constant geometry coefficients defined using the position-vector gradients in the reference configuration. The change in the position-vector gradients is used to define a velocity transformation matrix that leads to constant element inertia and stiffness matrices in the case of infinitesimal rotations. In contrast to conventional structural finite elements, the elements considered in this study can be used to describe the initial geometry with the same degree of accuracy as B-spline and nonuniform rational B-spline (NURBS) representations, widely used in the computer-aided design (CAD). An eigenvalue analysis is performed to evaluate the element convergence characteristics in the case of different geometries, including straight, tapered, and curved configurations. The frequencies obtained are compared with those obtained using a commercial FE software and analytical solutions. The stiffness matrix is obtained using both the general continuum mechanics (GCM) approach and the newly proposed strain split method (SSM) in order to investigate its effectiveness as a locking alleviation technique. 
    more » « less
  5. The translation equivariance of convolutional layers enables convolutional neural networks to generalize well on image problems. While translation equivariance provides a powerful inductive bias for images, we often additionally desire equivariance to other transformations, such as rotations, especially for non-image data. We propose a general method to construct a convolutional layer that is equivariant to transformations from any specified Lie group with a surjective exponential map. Incorporating equivariance to a new group requires implementing only the group exponential and logarithm maps, enabling rapid prototyping. Showcasing the simplicity and generality of our method, we apply the same model architecture to images, ball-and-stick molecular data, and Hamiltonian dynamical systems. For Hamiltonian systems, the equivariance of our models is especially impactful, leading to exact conservation of linear and angular momentum. 
    more » « less