Sub-Doppler Infrared Spectra and Torsion-Rotation Energy Manifold of Methanol in the CH-Stretch Fundamental Region

Exact tunneling splittings from symmetrized path integrals

We develop a new simulation technique based on path-integral molecular dynamics for calculating ground-state tunneling splitting patterns from ratios of symmetrized partition functions. In particular, molecular systems are rigorously projected onto their J = 0 rotational state by an "Eckart spring" that connects two adjacent beads in a ring polymer. Using this procedure, the tunneling splitting can be obtained from thermodynamic integration at just one (sufficiently low) temperature. Converged results are formally identical to the values that would have been obtained by solving the full rovibrational Schrödinger equation on a given Born-Oppenheimer potential energy surface. The new approach is showcased with simulations of hydronium and methanol, which are in good agreement with wavefunction-based calculations and experimental measurements. The method will be of particular use for the study of low-barrier methyl rotations and other floppy modes, where instanton theory is not valid.

On the separability of large-amplitude motions in anharmonic frequency calculations

Nuclear vibrational theories based upon the Watson Hamiltonian are ubiquitous in quantum chemistry, but are generally unable to model systems in which the wavefunction can delocalise over multiple energy minima, i.e. molecules that have low-energy torsion and inversion barriers. In a 2019 Chemical Reviews article, Puzzarini et al. note that a common workaround is to simply decouple these problematic modes from all other vibrations in the system during anharmonic frequency calculations. They also point out that this approximation can be "ill-suited", but do not quantify the errors introduced. In this work, we present the first systematic investigation into how separating out or constraining torsion and inversion vibrations within potential energy surface (PES) expansions affects the accuracy of computed fundamental wavenumbers for the remaining vibrational modes, using a test set of 19 tetratomic molecules for which high quality analytic potential energy surfaces and fully-coupled anharmonic reference fundamental frequencies are available. We find that the most effective and efficient strategy is to remove the mode in question from the PES expansion entirely. This introduces errors of up to +10 cm-1 in stretching fundamentals that would otherwise couple to the dropped mode, and ±5 cm-1 in all other fundamentals. These errors are approximately commensurate with, but not necessarily additional to, errors due to the choice of electronic structure model used in constructing spectroscopically accurate PES.

Diffusion Monte Carlo Calculations of Zero-Point Energies of Methanol and Deuterated Methanol

Diffusion Monte Carlo (DMC) simulations have been used to obtain quantum zero-point energies of methanol and all its isotopologs and isotopomers, using a new, accurate semi-global potential energy surface. This potential energy surface is a precise, permutationally invariant fit to 6676 ab initio energies, obtained at the CCSD(T)-F12b/aug-cc-pVDZ level of theory. Quantum zero-point energies of deuterated methanol isotopomers are very close to each other and so a simple statistical argument can be used to estimate the populations of each isotopomer at very low-temperatures. The DMC simulations also indicate that there is virtually zero probability for H/D exchange in the zero-point state. © 2019 Wiley Periodicals, Inc.

An extended E⊗e Jahn-Teller Hamiltonian for large-amplitude motion: Application to vibrational conical intersections in CH3SH and CH3OH

An extended E⊗e Jahn-Teller Hamiltonian is presented for the case where the (slow) nuclear motion extends far from the symmetry point and may be described approximately as motion on a sphere. Rather than the traditional power series expansion in the displacement from the C_3v symmetry point, an expansion in the spherical harmonics is employed. Application is made to the vibrational Jahn-Teller effect in CH₃XH, with X = S, O, where the equilibrium CXH angles are 83° and 72°, respectively. In addition to the symmetry-required conical intersection (CI) at the C_3v symmetry point, ab initio calculations reveal sets of six symmetry-allowed vibrational CIs in each molecule. The CIs for each molecule are arranged differently in the large-amplitude space, and that difference is reflected in the infrared spectra. The CIs in CH₃SH are found in both eclipsed and staggered geometries, whereas those for CH₃OH are found only in the eclipsed geometry near the torsional saddle point. This difference between the two molecules is reflected in the respective high-resolution spectra in the CH stretch fundamental region.

Automation of an "Aculight" continuous-wave optical parametric oscillator

We report the automation of a continuous-wave, singly resonant, optical parametric oscillator (Lockheed-Martin Aculight ARGOS 2400-SF-15). This commercially available optical parametric oscillator (OPO) is capable of producing >1 W of continuously tunable idler output between 2.2 and 4.6 μm. An algorithm based on the feedback from a high accuracy wavemeter is implemented to synchronize three separate OPO tuning elements; the translation of a fan-out type periodically poled lithium niobate crystal, the rotation of an intracavity etalon, and the continuous tuning of the pump and idler wavelengths via piezoelectric strain of the tunable fiber pump laser. This allows for several hundred wavenumbers of efficient, automatic, continuous tuning of the idler wave. Continuous feedback from the wavemeter limits the absolute frequency accuracy to ±20 MHz. The broad, automatic tuning of the OPO is demonstrated via its implementation as a probe laser for the infrared action spectroscopy of methanol solvated in helium nanodroplets. LabVIEW virtual instruments for the automation of this OPO laser system are reported, along with detailed schematics of the associated hardware developed at the University of Georgia.

Vibrational Levels of Methanol Calculated by the Reaction Path Version of MULTIMODE, Using an ab initio, Full-Dimensional Potential

An accurate potential energy surface has been determined for methanol from ab initio potential data at the CCSD(T) level of theory with an aug-cc-pVTZ basis. The resulting potential function is valid over all twelve vibrational degrees of freedom for all near-equilibrium and torsional configurations. A torsional reaction path has been derived for this potential, from which the low-lying vibrational levels of methanol have been calculated by the reaction path version of MULTIMODE. Comparisons with experiment and other calculations are made.

Torsion-vibration coupling in methanol: The adiabatic approximation and intramolecular vibrational redistribution scaling

The four-dimensional model Hamiltonian of Wang and Perry [J. Chem. Phys. 109, 10795 (1998)] is used to compare the approximate adiabatic separation of the torsion and CH stretches in methanol to an exact solution of the same Hamiltonian. The adiabatic approximation accounts for the pattern of the energy levels in the lowest torsional states, including the inverted tunneling splittings, but does not account for the pattern of systematic two- and four-fold near degeneracies at high torsional excitation. In the adiabatic basis, the nonadiabatic couplings mix the torsional and vibrational degrees of freedom and hence are a source for intramolecular vibrational redistribution (IVR). These IVR matrix elements are found to decrease by only a factor of 2 or 3 with each higher coupling order, in agreement with the results of Pearman and Gruebele [Z. Phys. Chem. Munich 214, 1439 (2000)]. This gentle scaling behavior, which contrasts with a steeper falloff with coupling order in more rigid molecules, points to a more important role for direct high-order couplings in torsional molecules. In this model, the scaling behavior derives from a single coupling term that is low order in the torsional angular momentum in combination with one-dimensional torsional functions that include contributions from many torsional angular momenta.

Internal Rotation and Spin Conversion of CH 3 OH in Solid para-Hydrogen

The quantum solid para-hydrogen (p-H2) has recently proven useful in matrix isolation spectroscopy. Spectral lines of compounds embedded in this host are unusually narrow, and several species have been reported to rotate in p-H2. We found that a p-H2 matrix inhibits rotation of isolated methanol (CH3OH) but still allows internal rotation about the C-O bond, with splittings of the E/A torsional doublet in internal rotation-coupled vibrational modes that are qualitatively consistent with those for CH3OH in the gaseous phase. This simplified high-resolution spectrum further revealed the slow conversion of nuclear spin symmetry from species E to species A in the host matrix, offering potential insight into nuclear spin conversion in astrophysical sources.

Vibrational overtone spectroscopy of jet-cooled methanol from 5000 to 14 000 cm−1

Spectra of jet-cooled methanol in the overtone and combination region from 5000 to 14 000 cm(-1) have been obtained by means of infrared laser-assisted photofragment spectroscopy. Many of the observed features are assigned to combination bands of the type nnu(1)+nu(6), nnu(1)+nu(8), and nnu(1)+nu(6)+nu(8) (n=1,2,3), where nu(1) is the OH stretch, nu(6) is the OH bend, and nu(8) is the CO stretch. These bands show sharp torsion-rotation structure with features as narrow as 0.1 cm(-1). We also observe CH stretch overtones that are weaker than the OH containing combination bands and lack distinct torsion-rotation structure above v(CH)=2. The extent of observed structure on these bands allows us to place limits on the intramolecular vibrational energy redistribution decay rates in the upper vibrational states. We report a global fit of the observed band centers to a simple expression involving low-order anharmonicity constants.

Organic Composition of C/1999 S4 (LINEAR): A Comet Formed Near Jupiter?

In the current paradigm, Oort cloud comets formed in the giant planets' region of the solar nebula, where temperatures and other conditions varied greatly. The measured compositions of four such comets (Halley, Hyakutake, Hale-Bopp, and Lee) are consistent with formation from interstellar ices in the cold nebular region beyond Uranus. The composition of comet C/1999 S4 (LINEAR) differs greatly, which suggests that its ices condensed from processed nebular gas, probably in the Jupiter-Saturn region. Its unusual organic composition may require reevaluation of the prebiotic organic material delivered to the young Earth by comets.

Torsional Splittings in Small-Amplitude Vibrational Fundamental States of Methanol-Type Molecules

Group-theoretical methods are used to show that inverted torsional splittings in fundamental levels of small-amplitude vibrations of methanol-like molecules can be parameterized and understood in terms of the energy level patterns induced when a pair of high-barrier torsionally split components of given v(t) and (t)A+(t)E symmetry species in the molecular symmetry group G(6) is allowed to interact with small-amplitude vibrational modes of symmetry (v)E. Such doubly degenerate (v)E vibrational modes arise rather naturally in G(6) (isomorphic with the point-group C(3v)) for those methyl-group vibrations in point-group-C(s) asymmetric tops such as CH(3)-CHO that are analogs of the degenerate methyl-group stretch, bend, and rocking vibrations in point-group-C(3v) symmetric tops such as CH(3)-C identical withC-H. The present group-theoretical treatment is somewhat different from, but (as a comparison of model parameters shows) still fundamentally similar to, the recent local mode explanation of inverted torsional splittings in the C-H stretching fundamental region in methanol. Copyright 2001 Academic Press.

Sub-Doppler Infrared Spectra of the OH-Stretch Fundamental of (13)C-Methanol

High-resolution infrared spectra of the (13)C-methanol OH stretching band, nu(1), were obtained by slit-jet absorption spectroscopy at a rotational temperature of 12-15 K. Twenty-eight subbands were assigned within the range of upper state quantum numbers, K' = 0 to 2 and J' = 0 to 10. The upper state energy levels are heavily perturbed; about half of the assigned subbands were found to be split by perturbations with matrix elements in the range 1-3 cm(-1). The doubled lines were deperturbed and together with the "unperturbed" lines were fitted to a global torsion-rotation Hamiltonian with root-mean-square deviation of 0.41 cm(-1) to yield a torsional barrier height of approximately 405 cm(-1) in the vibrationally excited state. Copyright 2000 Academic Press.

Dark state illuminated: infrared spectrum and inverted torsional structure of the nu(11) out-of-plane CH3-rocking mode of methanol

The nu(11) out-of-plane CH3-rocking band of methanol, one of the last fundamentals remaining en route to full vibrational description of this prototype internal rotor, has been identified in the Fourier transform spectrum of CH3OH in the 1150 cm(-1) region. The nu(11) torsional energy pattern is found to be inverted, the first such discovery for a bending state and an important pointer to potential generality of this phenomenon for the whole class of threefold internal rotor molecules. Level-crossing resonances giving transfer channels for intramolecular vibrational energy redistribution (IVR) have been found, and new insights are reported for far-infrared laser emission involving nu(11) levels.