Full-dimensional, ab initio potential energy surface for CH3OH → CH3+OH

Variational Vibrational States of Methanol (12D)

Full-dimensional (12D) vibrational states of the methanol molecule (CH3OH) have been computed using the GENIUSH-Smolyak approach and the potential energy surface from Qu and Bowman (2013). All vibrational energies are converged better than 0.5 cm-1 with respect to the basis and grid size up to the first overtone of the CO stretch, ca. 2000 cm-1 beyond the zero-point vibrational energy. About 70 torsion-vibration states are reported and assigned. The computed vibrational energies agree with the available experimental data within less than a few cm-1 in most cases, which confirms the good accuracy of the potential energy surface. The computations are carried out using curvilinear normal coordinates with the option of path-following coefficients, which minimize the coupling of the small- and large-amplitude motions. It is important to ensure tight numerical fulfillment of the C3v(M) molecular symmetry for every geometry and coefficient set used to define the curvilinear normal coordinates along the torsional coordinate to obtain a faithful description of degeneracy in this floppy system. The reported values may provide a computational reference for fundamental spectroscopy, astrochemistry, and for the search of the proton-to-electron mass ratio variation using the methanol molecule.

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.

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.

Torsional Potentials of Glyoxal, Oxalyl Halides, and Their Thiocarbonyl Derivatives: Challenges for Popular Density Functional Approximations

The reliability of popular density functionals was studied for the description of torsional profiles of 36 molecules: glyoxal, oxalyl halides, and their thiocarbonyl derivatives. HF and 18 functionals of varying complexity, from local density to range-separated hybrid approximations and double-hybrid, have been considered and benchmarked against CCSD(T)-level rotational profiles. For molecules containing heavy halogens, most functionals fail to reproduce barrier heights accurately and a number of functionals introduce spurious minima. Dispersion corrections show no improvement. Calibrated torsion-corrected atom-centered potentials rectify the shortcomings of PBE and also improve on σ-hole based intermolecular binding in dimers and crystals.

Stereoisomers of hydroxymethanes: Probing structural and spectroscopic features upon substitution

Ab initio studies on CH_x(OH)_4-x (x = 0-3) polyols are carried out to derive their structural and spectroscopic features. Several stereoisomers (both equilibrium structures and transition states) are found. Some are predicted here for the first time. We determined hence their geometrical parameters, vibrational frequencies, electronic excitation energies for the singlet manifold, and IR spectra. While the IR spectra for all polyols present similar shapes, their UV spectra exhibit however distinct band origin that are specific to each polyol and more interestingly to each diasteroisomer. Stereoelectronic effects are also noticed and discussed. It is suggested that UV spectroscopy is an efficient probe to experimentally identify polyols in mixtures involving polyols.