Perspective: Accurate ro-vibrational calculations on small molecules

Methane dimer rovibrational states and Raman transition moments

Benchmark-quality rovibrational data are reported for the methane dimer from variational nuclear motion computations using an ab initio intermolecular potential energy surface reported by [M. P. Metz et al., Phys. Chem. Chem. Phys., 2019, 21, 13504-13525]. A simple polarizability model is used to compute Raman transition moments that may be relevant for future direct observation of the intermolecular dynamics. Non-negligible ΔK ≠ 0 transition moments arise in this symmetric top system due to strong rovibrational couplings.

From the Automated Calculation of Potential Energy Surfaces to Accurate Infrared Spectra

Advances in the development of quantum chemical methods and progress in multicore architectures in computer science made the simulation of infrared spectra of isolated molecules competitive with respect to established experimental methods. Although it is mainly the multidimensional potential energy surface that controls the accuracy of these calculations, the subsequent vibrational structure calculations need to be carefully converged in order to yield accurate results. As both aspects need to be considered in a balanced way, we focus on approaches for molecules of up to 12-15 atoms with respect to both parts, which have been automated to some extent so that they can be employed in routine applications. Alternatives to machine learning will be discussed, which appear to be attractive, as long as local regions of the potential energy surface are sufficient. The automatization of these methods is still in its infancy, and the generalization to molecules with large amplitude motions or molecular clusters is far from trivial, but many systems relevant for astrophysical studies are already in reach.

Towards the generation of potential energy surfaces of weakly bound medium-sized molecular systems: the case of benzonitrile-He complex

Currently, the explicitly correlated coupled cluster method is used routinely to generate the multi-dimensional potential energy surfaces (mD-PESs) of van der Waals complexes of small molecular systems relevant for atmospheric, astrophysical and industrial applications. Although very accurate, this method is computationally prohibitive for medium and large molecules containing clusters. For instance, the recent detections of complex organic molecules (COMs) in the interstellar medium, such as benzonitrile, revealed the need to establish an accurate enough electronic structure approach to map the mD-PESs of these species interacting with the surrounding gases. As a benchmark, we have treated the case of the polar molecule benzonitrile interacting with helium, where we use post-Hartree-Fock and symmetry-adapted perturbation theory (SAPT) techniques. Accordingly, we show that MP2 and distinguishable-cluster approximation (DCSD) cannot be used for this purpose, whereas accurate enough PESs may be obtained using the corresponding explicitly correlated versions (MP2-F12 or DCSD-F12) with a reduction in computational costs. Alternatively, computations revealed that SAPT(DFT) is as performant as CCSD(T)-F12/aug-cc-pVTZ, making it the method of choice for mapping the mD-PESs of COMs containing clusters. Therefore, we have used this approach to generate the 3D-PES of the benzonitrile-He complex along the intermonomer Jacobi coordinates. As an application, we have incorporated the analytic form of this PES into quantum dynamical computations to determine the cross sections of the rotational (de-)excitation of benzonitrile colliding with helium at a collision energy of 10 cm-1.

Semiclassical and VSCF/VCI Calculations of the Vibrational Energies of trans- and gauche-Ethanol Using a CCSD(T) Potential Energy Surface

A recent full-dimensional Δ-Machine learning potential energy surface (PES) for ethanol is employed in semiclassical and vibrational self-consistent field (VSCF) and virtual-state configuration interaction (VCI) calculations, using MULTIMODE, to determine the anharmonic vibrational frequencies of vibration for both the trans and gauche conformers of ethanol. Both semiclassical and VSCF/VCI energies agree well with the experimental data. We find significant mixing between the VSCF basis states due to Fermi resonances between bending and stretching modes. The same effects are also accurately described by the full-dimensional semiclassical calculations. These are the first high-level anharmonic calculations using a PES, in particular a "gold-standard" CCSD(T) one.

Energy Landscape and Structural and Spectroscopic Characterization of Diazirine and Its Cyclic Isomers

Identifying new nitrogenated hydrocarbon molecules in the interstellar medium (ISM) is challenging because of the lack of comprehensive spectroscopic data from experiments. In this computational work, we focus on investigating the structures, relative energies, spectroscopic constants, and energy landscape of the cyclic isomers of diazirine (c-CH2N2) using ab initio quantum chemical methods. Density functional theory (DFT) methods and coupled cluster theory with singles and doubles including perturbative triples [CCSD(T)] and CCSD(T) with the explicitly correlated F12b correction [CCSD(T)-F12b] were employed for this purpose along with large correlation consistent cc-pVTZ, cc-pVQZ, and cc-pV5Z basis sets. Harmonic vibrational frequencies, infrared vibrational intensities, rotational constants, and dipole moments are reported. Anharmonic vibrational fundamentals along with centrifugal distortion constants, and vibration-rotation interaction constants are also reported for all the cyclic isomers. The energies computed with the CCSD(T) and CCSD(T)-F12b methods were extrapolated to the one-particle complete basis set (CBS) limit following a three-point formula. At the CCSD(T)-F12b/CBS level of theory, the 3,3H-diazirine (c-CH2N2) is the lowest energy cyclic isomer followed by 1,3H-diazirine, (E)-1,2H-diazirine, and (Z)-1,2H-diazirine, which are 20.1, 47.8, and 51.3 kcal mol-1 above the 3,3H-diazirine, respectively. Accurate structures and spectroscopic constants that are reported here could be useful for future identification of these cyclic nitrogenated organic molecules in the interstellar medium or circumstellar disks.

Calculation of absolute Raman scattering cross-sections using vibrational self-consistent field/vibrational configuration interaction wave functions

In the present study, the differential scattering cross-sections, depolarization ratios and Raman shifts of small molecular systems are obtained from configuration iteration wave functions of vibrational self-consistent field (VSCF) states. The transition polarizabilities were modeled using the Placzek approximation, neglecting those contributions not arising from the electric dipole mechanism. This theoretical approach is considered a good approximation for samples that absorb in the UV range if the excitation radiation falls in the visible region, as is the case of the molecules selected for the present study, namely: water, methane, and acetylene. Potential energy and electronic polarizability surfaces are calculated by the CCSD(T) and CC3 methods with aug-cc-p(C)V(T,Q,5)Z basis sets. The vibrational Hamiltonian includes the vibrational angular momentum contribution of the Watson kinetic energy operator. As expected, due to the variational nature of the VSCF and vibrational configuration interaction (VCI) methods, the Raman transition wavenumbers are substantially improved over the harmonic predictions. Surprisingly, the scattering cross-sections obtained using the harmonic approximation or the VSCF method better agrees with the experimental values than those cross-sections predicted using VCI wave functions. The more significant deviations of the VCI results from the experimental reference may be related to the significant uncertainties of the measured cross-sections. Still, it may also indicate that the VCI Raman transition moments may require a more accurate description of the electronic polarizability surface. Finally, the depolarization ratios calculated for H₂ O and C₂ D₂ using harmonic and VCI wave functions have similar accuracy, whereas, for C₂ H₂ and C₂ HD, the VCI results are more accurate.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

General Perturb-Then-Diagonalize Model for the Vibrational Frequencies and Intensities of Molecules Belonging to Abelian and Non-Abelian Symmetry Groups

In this paper, we show that the standard second-order vibrational perturbation theory (VPT2) for Abelian groups can be used also for non-Abelian groups without employing specific equations for two- or threefold degenerate vibrations but rather handling in the proper way all the degeneracy issues and deriving the peculiar spectroscopic signatures of non-Abelian groups (e.g., l -doubling) by a posteriori transformations of the eigenfunctions. Comparison with the results of previous conventional implementations shows a perfect agreement for the vibrational energies of linear and symmetric tops, thus paving the route to the transparent extension of the equations already available for asymmetric tops to the energies of spherical tops and the infrared and Raman intensities of molecules belonging to non-Abelian symmetry groups. The whole procedure has been implemented in our general engine for vibro-rotational computations beyond the rigid rotor/harmonic oscillator model and has been validated on a number of test cases.

MULTIMODE Calculations of Vibrational Spectroscopy and 1d Interconformer Tunneling Dynamics in Glycine Using a Full-Dimensional Potential Energy Surface

A full-dimensional, permutationally invariant polynomial potential energy surface for glycine recently reported (R. Conte et al., J. Chem. Phys. 2020, 153, 244301) is used with the code MULTIMODE to determine the IR absorption spectra for Conformers I and II using a new separable dipole moment function. The calculated spectra agree well with the experimental ones. The full-dimensional nature of the potential allows us also to examine dynamical results, such as tunneling rates. Remarkably, using a one-dimensional path based on the potential energy surface to estimate the tunneling rate from Conformer VI to Conformer I, good agreement is found with the recent experimental measurement. Finally a brief comparison of our potential energy surface with a recently reported sGDML one is made.

On the synergy of matrix-isolation infrared spectroscopy and vibrational configuration interaction computations

The key feature of matrix-isolation infrared (MI-IR) spectroscopy is the isolation of single guest molecules in a host system at cryogenic conditions. The matrix mostly hinders rotation of the guest molecule, providing access to pure vibrational features. Vibrational self-consistent field (VSCF) and configuration interaction computations (VCI) on ab initio multimode potential energy surfaces (PES) give rise to anharmonic vibrational spectra. In a single-sourced combination of these experimental and computational approaches, we have established an iterative spectroscopic characterization procedure. The present article reviews the scope of this procedure by highlighting the strengths and limitations based on the examples of water, carbon dioxide, methane, methanol, and fluoroethane. An assessment of setups for the construction of the multimode PES on the example of methanol demonstrates that CCSD(T)-F12 level of theory is preferable to compute (a) accurate vibrational frequencies and (b) equilibrium or vibrationally averaged structural parameters. Our procedure has allowed us to uniquely assign unknown or disputed bands and enabled us to clarify problematic spectral regions that are crowded with combination bands and overtones. Besides spectroscopic assignment, the excellent agreement between theory and experiment paves the way to tackle questions of rather fundamental nature as to whether or not matrix effects are systematic, and it shows the limits of conventional notations used by spectroscopists.

Decomposing anharmonicity and mode-coupling from matrix effects in the IR spectra of matrix-isolated carbon dioxide and methane

A combined experimental and computational approach revealed similarities and differences in the vibrational signature of matrix-isolated carbon dioxide and methane.

General mathematical formulation of scattering processes in atom-diatomic collisions in the RmatReact methodology

Accurately modelling cold and ultracold reactive collisions occurring over deep potential wells, such as [Formula: see text], requires the development of new theoretical and computational methodologies. One potentially useful framework is the R-matrix method adopted widely for electron-molecule collisions which has more recently been applied to non-reactive heavy-particle collisions such as Ar-Ar. The existing treatment of non-reactive elastic and inelastic scattering needs to be substantially extended to enable modelling of reactive collisions: this is the subject of this paper. Herein, we develop the general mathematical formulation for non-reactive elastic and inelastic scattering, photoassociation, photodissociation, charge exchange and reactive scattering using the R-matrix method. Of particular note is that the inner region, of central importance to calculable R-matrix methodologies, must be finite in all scattering coordinates rather than a single scattering coordinate as for non-reactive scattering. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H₃⁺, H₅⁺ and beyond'.

The smallest proton-bound dimer H5+: theoretical progress

The protonated hydrogen dimer, H5+, is the smallest system including proton transfer, and has been of long-standing interest since its first laboratory observation in 1962. H5+ and its isotopologues are the intermediate complexes in deuterium fractionation reactions, and are of central importance in molecular astrophysics. The recently recorded infrared spectra of both H5+ and D5+ reveal a rich vibrational dynamics of the cations, which presents a challenge for standard theoretical approaches. Although H5+ is a four-electron ion, which makes highly accurate electronic structure calculations tractable, the construction of ab initio-based potential energy and dipole moment surfaces has proved a hard task. In the same vein, the difficulties in treating the nuclear motion could also become cumbersome due to their high dimensionality, floppiness and/or symmetry. These systems are prototypical examples for studying large-amplitude motions, as they are highly delocalized, interconverting between equivalent minima through internal rotation and proton transfer motions requiring state-of-the-art treatments. Recent advances in the computational vibrational spectroscopy of the H5+ cation and its isotopologues are reported from full quantum spectral simulations, providing important information in a rigorous manner, and open perspectives for further future investigations. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

Highly accurate experimentally determined energy levels of H3. J Chem Phys 2019;150:214303. [PMID: 31176322 DOI: 10.1063/1.5099454]

A sub-Doppler rovibrational spectroscopic survey of H₃ ⁺ has been conducted which included 36 transitions in the ν₂ ← 0 fundamental band, 15 transitions in the 2ν₂ ²←ν₂ hot band, and 7 transitions in the 2ν₂ ²←0 overtone band, improving the uncertainties of most transitions by more than an order of magnitude to ∼4 MHz. Combination differences were used to determine relative energy levels and forbidden rotational transitions up to J = 6. A fit of the ground state to an effective Hamiltonian was used to connect ortho and para states, and to determine the absolute energy levels relative to the forbidden (0, 0) state. Ultimately, 62 rovibrational energy levels in the ground, ν₂, and 2ν₂ ² states were determined with ∼10 MHz uncertainty. Comparing the experimentally determined energy levels with ab initio calculations revealed an unexpected dependence of the residuals on the quantum number G.

Theoretical rotation-vibration spectroscopy of cis- and trans-diphosphene (P2H2) and the deuterated species P2HD

Growing astronomical interest in phosphorous (P) chemistry is stimulating the search for new interstellar P-bearing molecules, a task requiring detailed knowledge of the microwave and infrared molecular spectrum. In this work, we present comprehensive rotation-vibration line lists of the cis- and trans-isomers of diphosphene (P₂H₂). The line lists have been generated using robust, first-principles methodologies based on newly computed, high-level ab initio potential energy and dipole moment surfaces. Transitions are considered between states with energies up to 8000 cm^-1 and total angular momentum J ≤ 25. These are the first-ever line lists to be reported for P₂H₂, and they should significantly facilitate future spectroscopic characterization of this system. The deuterated species trans-P₂HD and the effect of its dynamic dipole moment on the rovibrational spectrum are also discussed.

The Hunt for Elusive Molecules: Insights from Joint Theoretical and Experimental Investigations

Rotational spectroscopy is an invaluable tool to unambiguously determine the molecular structure of a species, and sometimes even to establish its very existence. This article illustrates how experimental and theoretical state-of-the-art tools can be used in tandem to investigate the rotational structure of molecules, with particular emphasis on those that have long remained elusive. The examples of three emblematic species-gauche-butadiene, disilicon carbide, and germanium dicarbide-highlight the close, mutually beneficial interaction between high-level theoretical calculations and sensitive microwave measurements. Prospects to detect other elusive molecules of chemical and astronomical interest are discussed.

Quantum approaches to vibrational dynamics and spectroscopy: is ease of interpretation sacrificed as rigor increases?

The subject of this Perspective is quantum approaches, beyond the harmonic approximation, to vibrational dynamics and IR spectroscopy.

The rotation–vibration spectrum of methyl fluoride from first principles

State-of-the-art theoretical methods are used to investigate the rotation–vibration spectrum of methyl fluoride.

Non-adiabatic mass correction to the rovibrational states of molecules: Numerical application for the H 2 + molecular ion

General transformation expressions of the second-order non-adiabatic Hamiltonian of the atomic nuclei, including the kinetic-energy correction terms, are derived upon the change from laboratory-fixed Cartesian coordinates to general curvilinear coordinate systems commonly used in rovibrational computations. The kinetic-energy or so-called "mass-correction" tensor elements are computed with the stochastic variational method and floating explicitly correlated Gaussian functions for the H 2 + molecular ion in its ground electronic state. {Further numerical applications for the ⁴ He 2 + molecular ion are presented in the forthcoming paper, Paper II [E. Mátyus, J. Chem. Phys. 149, 194112 (2018)]}. The general, curvilinear non-adiabatic kinetic energy operator expressions are used in the examples, and non-adiabatic rovibrational energies and corrections are determined by solving the rovibrational Schrödinger equation including the diagonal Born-Oppenheimer as well as the mass-tensor corrections.

Characterization of Azirine and Its Structural Isomers

The structures and spectroscopic properties of azirine (C₂H₃N), a nitrogen-containing three-membered cyclic molecule, and its isomers were studied with state-of-the-art ab initio quantum chemical methods. Azirine is isomeric with methyl cyanide (CH₃CN) and methyl isocyanide (CH₃NC)-both observed in the star-forming regions of Sgr B2. In this study, we characterize the stationary points on the potential energy surface, relative energies, dipole moments, rotational constants, and harmonic vibrational frequencies of the 2 H-azirine ( a), 1 H-2,2 H-azirine ( b, carbene isomer), and 1 H-azirine ( c) cyclic isomers. The CCSD(T) method and density functional theory (DFT), using the ωB97-X functional, along with Dunning's cc-pVXZ (X = T and Q) basis sets were used to optimize molecular geometries and calculate vibrational frequencies. The 2 H-azirine, an imine isomer ( a), was found to be the lowest in energy among the cyclic isomers, followed by the carbene isomer ( b), and last the 1 H-azirine, an enamine isomer ( c). All three cyclic isomers have a C _s symmetry equilibrium structure. Azirines, if identified (three linear C₂H₃N isomers are already identified in the same source toward the galactic center, Sgr B2), would be the first nitrogen-containing cyclic molecules identified in an astronomical observation.

Quantum and classical IR spectra of (HCOOH)2, (DCOOH)2 and (DCOOD)2 using ab initio potential energy and dipole moment surfaces

Full-dimensional (24 modes) quantum calculation of the IR spectrum of (DCOOD)₂, and comparison with classical MD one.

Optical Measurement of Radiocarbon below Unity Fraction Modern by Linear Absorption Spectroscopy

High-precision measurements of radiocarbon (¹⁴C) near or below a fraction modern ¹⁴C of 1 (F¹⁴C ≤ 1) are challenging and costly. An accurate, ultrasensitive linear absorption approach to detecting ¹⁴C would provide a simple and robust benchtop alternative to off-site accelerator mass spectrometry facilities. Here we report the quantitative measurement of ¹⁴C in gas-phase samples of CO₂ with F¹⁴C < 1 using cavity ring-down spectroscopy in the linear absorption regime. Repeated analysis of CO₂ derived from the combustion of either biogenic or petrogenic sources revealed a robust ability to differentiate samples with F¹⁴C < 1. With a combined uncertainty of ¹⁴C/¹²C = 130 fmol/mol (F¹⁴C = 0.11), initial performance of the calibration-free instrument is sufficient to investigate a variety of applications in radiocarbon measurement science including the study of biofuels and bioplastics, illicitly traded specimens, bomb dating, and atmospheric transport.

Investigation of Nonadiabatic Effects for the Vibrational Spectrum of a Triatomic Molecule: The Use of a Single Potential Energy Surface with Distance-Dependent Masses for H3. J Phys Chem A 2017;121:7016-7030. [PMID: 28820589 DOI: 10.1021/acs.jpca.7b04703]

On the basis of first-principles, the influence of nonadiabatic effects on the vibrational bound states of H₃⁺ has been investigated using distance-dependent reduced masses and only one single potential energy surface. For these new vibrational calculations, potentials based on explicitly correlated wave functions are used where, in addition, adiabatic corrections and relativistic contributions are taken into account. For the first time, several different fully distance-dependent reduced mass surfaces in three dimensions have been incorporated in the vibrational calculations.

Advances in spectroscopy and dynamics of small and medium sized molecules and clusters

Investigations of the spectroscopy and dynamics of small- and medium-sized molecules and clusters represent a hot topic in atmospheric chemistry, biology, physics, atto- and femto-chemistry and astrophysics.