Derivation of rotation-vibration Hamiltonians that satisfy the Casimir condition

Umbrella inversion of ammonia redux

Umbrella inversion of ammonia is a prototypical example of large-amplitude vibrational motion, described with a symmetric double-well potential. The transition state of the latter corresponds to a planar (D_3h) molecular geometry, whereas the two equilibrium configurations are equivalent (C_3v) pyramidal structures, with the nitrogen atom being either 'above' or 'below' the plane of the hydrogen atoms. As commonly understood, inversion motion of ammonia corresponds to the coherent, anharmonic, vibrational motion of the molecule, which shuttles back and forth between the two potential wells; that is, oscillation of the nitrogen atom from one side of the H₃ plane to the other, via coherent tunneling. However, this intuitively appealing view of umbrella inversion results from a reduced description of the dynamics, which includes only the inversion vibrational coordinate and fully neglects all the other molecular degrees of freedom. As such, this textbook picture of inversion motion ignores the fact that the two equilibrium structures of ammonia are superimposable, and can only be distinguished by labelling the identical hydrogen nuclei. A correct description of umbrella inversion, which incorporates nuclear permutations, requires the inclusion of other molecular modes. Indeed, it is well known that the quantum symmetrization postulate engenders entanglement between ammonia's nuclear-spin, inversion, and rotation. Using the explicit expressions of the corresponding zeroth-order eigenstates, we clearly show that the inversion density of any multilevel wavepacket of ammonia, including the case of perfectly aligned molecules, is symmetrically distributed between the two potential wells, at all times. This follows from a rigorous demonstration based on the evaluation of the expectation values of the inversion coordinate or equivalent projection operators. However, provided that these wavepackets involve inversion-rotation levels with opposite parity, the inversion density may exhibit dynamical spatial localization. In the latter case, the space-fixed inversion density or, equivalently, the expectation values of the projections of the inversion coordinate on the space-fixed axes, may oscillate between opposite directions in the space-fixed frame. Nevertheless, in all cases, localization of ammonia in a single potential well is impossible, even partially or transiently. This is equivalent to saying that the nitrogen atom has the same probability (one-half) to be on either side of the H₃ plane, for any wavepacket of the molecule and at all times-a conclusion which is in perfect accord with the principle of the indistinguishability of identical particles (nuclei).

On neglecting Coriolis and related couplings in first-principles rovibrational spectroscopy: Considerations of symmetry, accuracy, and simplicity. II. Case studies for H2O isotopologues, H3+, O3, and NH3

For centuries, it has been known that vibrational and rotational degrees of freedom are in general not separable. Nevertheless, surprisingly little is known about the best strategies for approximately separating these degrees of freedom in practice-even in the case of semirigid molecules, where the separation is most meaningful. There is also some confusion in the literature about the proper way to quantify the magnitude of the Coriolis (i.e., rotation-vibration) coupling in rovibrational Hamiltonians or its effect on the rovibrational eigenenergies. In this study, a vibrational-coordinate-independent metric is proposed to quantify the magnitude of the Coriolis contribution to the rovibrational Hamiltonian. The impact of Coriolis coupling on the rovibrational eigenenergies is computed numerically exactly, using both full and various truncated Hamiltonians. The role played by the choice of the vibrational coordinate system-and especially by the choice of "embedding" or body-fixed frame-is examined extensively, both numerically and analytically. This investigation targets several molecular prototypes, all of which serve as important benchmarks for the high-resolution spectroscopic community. Most of these are triatomic molecules, including water (H₂¹⁶O), its deuterated isotopologues (D₂¹⁶O and HD¹⁶O), H₃⁺, and ozone (¹⁶O₃), but the tetratomic ammonia molecule (¹⁴NH₃) is also investigated. These studies provide important insight into the nature of Coriolis coupling under various circumstances. The findings of this study also have significant practical ramifications, vis-à-vis the use of simplifying numerical approximation techniques in nuclear-motion computations.

On neglecting Coriolis and related couplings in first-principles rovibrational spectroscopy: considerations of symmetry, accuracy, and simplicity

The rotation-vibration (Coriolis) coupling contribution to variationally computed rovibrational energy levels is investigated, employing triatomic AB[Formula: see text] molecules as models. In particular, calculations are performed for H[Formula: see text][Formula: see text]O, across a range of vibrational and rotational excitations, both with and without the Coriolis contribution. A variety of different embedding choices are considered, together with a hierarchy of increasingly severe approximations culminating in a generalized version of the so-called "centrifugal sudden" method. Several surprising and remarkable conclusions are found, including that the Eckart embedding is not the best embedding choice.

Eckart-Sayvetz conditions revisited

It is shown that vibrational displacements satisfying the Eckart-Sayvetz conditions can be constructed by projection of unconstrained displacements. This result has a number of interesting direct and indirect ramifications: (i) The normal coordinates corresponding to an electronic state or an isotopologue of a molecule are transformed to those of another state or isotopologue by a linear and, in general, non-orthogonal transformation. (ii) Novel interpretation of axis switching. (iii) One may enhance the separation of rotational-large-amplitude internal motions and the vibrational motions beyond that offered by the standard use of the Eckart-Sayvetz conditions. (iv) The rotational-vibrational Hamiltonian given in terms of curvilinear internal coordinates may be derived with elementary mathematical tools while taking into account the Eckart conditions with or without enhancement.

Numerically constructed internal-coordinate Hamiltonian with Eckart embedding and its application for the inversion tunneling of ammonia

It is shown that the use of an Eckart-frame embedding with a kinetic energy operator expressed in curvilinear internal coordinates becomes feasible and straightforward to implement for arbitrary molecular compositions and internal coordinates if the operator is defined numerically over a (discrete variable representation) grid. The algorithm proposed utilizes the transformation method of Dymarsky and Kudin to maintain the rotational Eckart condition. In order to demonstrate the applicability and flexibility of our approach the non-rigid ammonia molecule is considered and the corresponding rotational-vibrational energy levels and wave functions are computed using kinetic energy operators with three different embeddings. Two of them fulfill the Eckart conditions corresponding to a trigonal pyramidal (C3v) and a trigonal planar (D3h) reference structure and the third one is a non-Eckart frame. The computed energy levels are, of course, identical, and the structure of the three different wave-function representations are analyzed in terms of the rigid rotor functions for a symmetric top. The possible advantages of one frame representation over another are discussed concerning the interpretation of the rovibrational states in terms of the traditional rigid rotor labels.

Combined perturbative-variational investigation of the vibrations of CHBr3 and CDBr3

A full dimensional vibrational treatment of CHBr(3) and CDBr(3) using Van Vleck perturbation theory followed by a variational calculation is presented. The calculation of a force field, and its adjustment for better match with experiment, is discussed. The computed eigenstates and spectral features are compared to experiment. Changes in intensities of the nu(1) and 2nu(4) bands upon simple alterations of the dipole moment expansion are described.

Body frames in the separation of collective angles in quantum N-body problems

The application of the concept of body-fixed reference frames, proposed by C. Eckart [Phys. Rev. 47, 552 (1935)], to the problem of the separation of three collective angles in quantum N-body problems is analyzed based on the technique recently developed by Meremianin and Briggs [Phys. Rep. 384, 121 (2003)]. Special attention is paid to the body frame defined by the "second Eckart condition" which minimizes vibro-rotational couplings near the equilibrium position. The important case of the Eckart frame for three-body systems is considered in detail. The connection of the basis vectors of the Eckart frame with Jacobi vectors is derived. All results of this work are valid for an arbitrary choice of internal (body-frame) coordinates.