skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Incorporating reduced axis system embedding into ab initio tunnelling-rotation Hamiltonians with curvilinear vibrational Møller–Plesset perturbation theory
Large-amplitude vibrational motion influences the rovibrational structure of molecules that tunnel between multiple wells. Reaction path (RP) Hamiltonians, and curvilinear coordinates more gen- erally, are useful for modelling pure vibrational motion in these systems and provide a practical framework for calculating accurate ab initio anharmonic vibrational energies and tunnelling split- tings with perturbation theory. These computational tools also offer the means to address rotation- vibration coupling associated with large-amplitude motion in rotating molecules. In this paper, we incorporate the reduced axis system (RAS) frameembeddingwithRPHamiltoniansandsecond-order vibrational Møller-Plesset perturbation theory (VMP2). Because the RP-RAS Hamiltonian eliminates rotation-vibration momentumcoupling everywhere along a one-dimensional reaction path, it is well suited for rovibrational VMP2 methods, the convergence of which relies critically on approximate vibration-vibration and vibration-rotation separability. The accuracy of this combined RP-RAS-VMP2 scheme is demonstrated by comparisons with numerically exact variational calculations and VMP2 parameters based on traditional Eckart embeddings for reduced-dimension models of torsional tunnelling in hydrogenperoxideandinversion tunnelling in cyclopropyl radical. Thefavourablecom- putational scaling ofVMP2makes it a promising strategy for calculating accurate tunnelling-rotation parameters for medium-sized and larger molecules in full dimensionality.  more » « less
Award ID(s):
2110489
PAR ID:
10528750
Author(s) / Creator(s):
Publisher / Repository:
Taylor & Francis
Date Published:
Journal Name:
Molecular Physics
Volume:
2024
ISSN:
0026-8976
Page Range / eLocation ID:
e2304098
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Quantum state-resolved spectroscopy was recently achieved for C60 molecules when cooled by buffer gas collisions and probed with a midinfrared frequency comb. This rovibrational quantum state resolution for the largest molecule on record is facilitated by the remarkable symmetry and rigidity of C60, which also present new opportunities and challenges to explore energy transfer between quantum states in this many-atom system. Here we combine state-specific optical pumping, buffer gas collisions, and ultrasensitive intracavity nonlinear spectroscopy to initiate and probe the rotation-vibration energy transfer and relaxation. This approach provides the first detailed characterization of C60 collisional energy transfer for a variety of collision partners, and determines the rotational and vibrational inelastic collision cross sections. These results compare well with our theoretical modeling of the collisions, and establish a route towards quantum state control of a new class of unprecedentedly large molecules. 
    more » « less
  2. Context.Accurate42Ti(p,γ)43V reaction rates are crucial for understanding the nucleosynthesis path of the rapid capture process (rpprocess) that occurs in X-ray bursts. Aims.We aim to improve the thermonuclear rates of42Ti(p,γ)43V based on more complete resonance information and a more accurate direct component, together with the recently released nuclear masses data. We also explore the impact of the newly obtained rates on therpprocess. Methods.We reevaluated the reaction rate of42Ti(p,γ)43V by the sum of the isolated resonance contribution instead of the Hauser-Feshbach statistical model. We used a Monte Carlo method to derive the associated uncertainties of new rates. The nucleosynthesis simulations were performed via the NuGrid post-processing code ppn. Results.The new rates differ from previous estimations due to the use of a series of updated resonance parameters and a direct S factor. Compared with the previous results from the Hauser-Feshbach statistical model, which assumes compound nucleus43V with a sufficiently high-level density in the energy region of astrophysical interest, large differences exist over the entire temperature region ofrp-process interest, up to two orders of magnitude. We consistently calculated the photodisintegration rate using our new nuclear masses via the detailed balance principle, and found the discrepancies among the different reverse rates are much larger than those for the forward rate, up to ten orders of magnitude at the temperature of 108K. Using a trajectory with a peak temperature of 1.95×109K, we performed therp-process nucleosynthesis simulations to investigate the impact of the new rates. Our calculations show that the adoption of the new forward and reverse rates result in abundance variations for Sc and Ca of 128% and 49%, respectively, compared to the variations for the statistical model rates. On the other hand, the overall abundance pattern is not significantly affected. The results of using new rates also confirm that therp-process path does not bypass the isotope43V. Conclusions.Our study found that the Hauser-Feshbach statistical model is inappropriate to the reaction rate evaluation for42Ti(p,γ)43V. The adoption of the new rates confirms that the reaction path of42Ti(p,γ)43V(p,γ)44Cr(β+)44V is a key branch of therpprocess in X-ray bursts. 
    more » « less
  3. In this work, we present the spectral analysis of 1 H- and 2 H-1,2,3-triazole vibrationally excited states alongside provisional and practical computational predictions of the excited-state quartic centrifugal distortion constants. The low-energy fundamental vibrational states of 1 H-1,2,3-triazole and five of its deuteriated isotopologues ([1- 2 H]-, [4- 2 H]-, [5- 2 H]-, [4,5- 2 H]-, and [1,4,5- 2 H]-1 H-1,2,3-triazole), as well as those of 2 H-1,2,3-triazole and five of its deuteriated isotopologues ([2- 2 H]-, [4- 2 H]-, [2,4- 2 H]-, [4,5- 2 H]-, and [2,4,5- 2 H]-2 H-1,2,3-triazole), are studied using millimeter-wave spectroscopy in the 130–375 GHz frequency region. The normal and [2- 2 H]-isotopologues of 2 H-1,2,3-triazole are also analyzed using high-resolution infrared spectroscopy, determining the precise energies of three of their low-energy fundamental states. The resulting spectroscopic constants for each of the vibrationally excited states are reported for the first time. Coupled-cluster vibration–rotation interaction constants are compared with each of their experimentally determined values, often showing agreement within 500 kHz. Newly available coupled-cluster predictions of the excited-state quartic centrifugal distortion constants based on fourth-order vibrational perturbation theory are benchmarked using a large number of the 1,2,3-triazole tautomer isotopologues and vibrationally excited states studied. 
    more » « less
  4. The vibrational spectra of condensed and gas-phase systems are influenced by thequantum-mechanical behavior of light nuclei. Full-dimensional simulations of approximate quantum dynamics are possible thanks to the imaginary time path-integral (PI) formulation of quantum statistical mechanics, albeit at a high computational cost which increases sharply with decreasing temperature. By leveraging advances in machine-learned coarse-graining, we develop a PI method with the reduced computational cost of a classical simulation. We also propose a simple temperature elevation scheme to significantly attenuate the artifacts of standard PI approaches as well as eliminate the unfavorable temperature scaling of the computational cost. We illustrate the approach, by calculating vibrational spectra using standard models of water molecules and bulk water, demonstrating significant computational savings and dramatically improved accuracy compared to more expensive reference approaches. Our simple, efficient, and accurate method has prospects for routine calculations of vibrational spectra for a wide range of molecular systems - with an explicit treatment of the quantum nature of nuclei. 
    more » « less
  5. null (Ed.)
    Scattering resonances above dissociation threshold are computed for four isotopically substituted ozone species: 16 O 18 O 16 O, 16 O 16 O 18 O, 18 O 16 O 18 O and 16 O 18 O 18 O, using a variational method with accurate treatment of the rotation–vibration coupling terms (Coriolis effect) for all values of the total angular momentum J from 0 to 4. To make these calculations numerically affordable, a new approach was developed which employs one vibrational basis set optimized for a typical rotational excitation ( J , Λ ), to run coupled rotation–vibration calculations at several desired values of J . In order to quantify the effect of Coriolis coupling, new data are contrasted with those computed using the symmetric-top rotor approximation, where the rotation–vibration coupling terms are neglected. It is found that, overall, the major properties of scattering resonances (such as their lifetimes, the number of these states, and their cumulative partition function Q ) are all influenced by the Coriolis effect and this influence grows as the angular momentum J is raised. However, it is found that the four isotopically substituted ozone molecules are affected roughly equally by the Coriolis coupling. When the ratio η of partition functions for asymmetric over symmetric ozone molecules is computed, the Coriolis effect largely cancels, and this cancelation seems to occur for all values of J . Therefore, it does not seem grounded to attribute any appreciable mass-independent symmetry-driven isotopic fractionation to the Coriolis coupling effect. 
    more » « less