Bridging Theory with Experiment: A Benchmark Study of Thermally Averaged Structural and Effective Spectroscopic Parameters of the Water Molecule

Intermolecular Bending States and Tunneling Splittings of Water Trimer from Rigorous 9D Quantum Calculations: I. Methodology, Energy Levels, and Low-Frequency Spectrum

We present the computational methodology that enables the first rigorous nine-dimensional (9D) quantum calculations of the intermolecular bending states of the water trimer, as well as its low-frequency spectrum for direct comparison with experiment. The water monomers, treated as rigid, have their centers of mass (cm's) at the corners of an equilateral triangle, and the intermonomer cm-to-cm distance is set to a value slightly larger than that in the equilibrium geometry of the trimer. The remaining nine strongly coupled large-amplitude bending (angular) degrees of freedom (DOFs) enter the 9D bend Hamiltonian of the three coupled 3D rigid-water hindered rotors. Its 9D eigenstates encompass excited librational vibrations of the trimer, as well as their torsional and bifurcation tunneling splittings, which have been the subject of much interest. The calculations of these eigenstates are extremely demanding, and a sophisticated computational scheme is developed that exploits the molecular symmetry group of the water trimer, G48, in order to make them feasible in a reasonable amount of time. The spectrum of the low-frequency vibrations of the water trimer simulated using the eigenstates of the 9D bend Hamiltonian agrees remarkably well with the experimentally observed far-infrared (FIR) spectrum of the trimer in helium nanodroplets over the entire frequency range of the measurements from 70 to 620 cm-1. This shows that most peaks in the experimental FIR spectrum are associated with the intermolecular bending vibrations of the trimer. Moreover, the ground-state torsional tunneling splittings from the present 9D calculations are in excellent agreement with the spectroscopic data. These results demonstrate the high quality of the ab initio 2 + 3-body PES employed for the DOFs included in the bound-state calculations.

Solving the Puzzle of the Carbonic Acid Vibrational Spectrum - an Anharmonic Story

Against the general belief that carbonic acid is too unstable for synthesis, it was possible to synthesize the solid[1,2] as well as gas-phase carbonic acid.[3] It was suggested that solid carbonic acid might exist in Earth's upper troposphere and in the harsh environments of other solar bodies,[4] where it undergoes a cycle of synthesis, decomposition, and dimerization.[5] To provide spectroscopic data for probing the existence of extraterrestrial carbonic acid,[2,6] matrix-isolation infrared (MI-IR) spectroscopy has shown to be essential.[3,4,6-8] However, early assignments within the harmonic approximation using scaling factors impeded a full interpretation of the rather complex MI-IR spectrum of H2CO3. Recently, carbonic acid was detected in the Galactic center molecular cloud,[9] triggering new interest in the anharmonic spectrum.[10] In this regard, we substantially reassign our argon MI-IR spectra based on accurate anharmonic calculations. We calculate a four-mode potential energy surface (PES) at the explicitly correlated coupled-cluster theory using up to triple-zeta basis sets, i. e., CCSD(T)-F12/cc-pVTZ-F12. On this PES, we perform vibrational self-consistent field and configuration interaction (VSCF/VCI) calculations to obtain accurate vibrational transition frequencies and resonance analysis of the fundamentals, first overtones, and combination bands. In total, 12 new bands can be assigned, extending the spectral data for carbonic acid and thus simplifying detection in more complex environments. Furthermore, we clarify disputed assignments between the cc- and ct-conformer.

Quantum Definition of Molecular Structure

Molecular structure, a key concept of chemistry, has remained elusive from the perspective of all-particle quantum mechanics, despite many efforts. Viewing molecular structure as a manifestation of strong statistical correlation between nuclear positions, we propose a practical method based on Markov chain Monte Carlo sampling and unsupervised machine learning. Application to the D3+ molecule unambiguously shows that it possesses an equilateral triangular structure. These results provide a major step forward in our understanding of the molecular structure from fundamental quantum principles.

Comprehensive quantum calculation of the first dielectric virial coefficient of water

We present a complete calculation, fully accounting for quantum effects and for molecular flexibility, of the first dielectric virial coefficient of water and its isotopologues. The contribution of the electronic polarizability is computed from a state-of-the-art intramolecular potential and polarizability surface from the literature, and its small temperature dependence is quantified. The dipolar polarizability is calculated in a similar manner with an accurate literature dipole-moment surface; it differs from the classical result both due to the different molecular geometries sampled at different temperatures and due to the quantization of rotation. We calculate the dipolar contribution independently from spectroscopic information in the HITRAN2020 database and find that the two methods yield consistent results. The resulting first dielectric virial coefficient provides a complete description of the dielectric constant at low density that can be used in humidity metrology and as a boundary condition for new formulations for the static dielectric constant of water and heavy water.

Comprehensive quantum chemical analysis of the (ro)vibrational spectrum of thiirane and its deuterated isotopologue

The (ro)vibrational spectra of thiirane, c-C2H4S, and its fully deuterated isotopologue, c-C2D4S, have been studied by means of vibrational configuration interaction theory, VCI, its incremental variant, iVCI, and subsequent variational rovibrational calculations, RVCI, which rely on multidimensional potential energy surfaces of coupled-cluster quality including up to four-mode coupling terms. Accurate geometrical parameters, fundamental vibrational transitions and first overtones, rovibrational spectra and rotational spectroscopic constants have been determined from these calculations and were compared with experimental results whenever available. A number of tentative misassignments in the vibrational spectra could be resolved and most results for the deuterated thiirane are high-level predictions, which may guide experiments to come. Besides this, a new implementation of infrared intensities within the iVCI framework has been tested for the transitions of the title compounds and are compared with results obtained from standard VCI calculations.

Quantum Chemical Investigation of the Cold Water Dimer Spectrum in the First OH-Stretching Overtone Region Provides a New Interpretation

Intramolecular vibrational transition wavenumbers and intensities were calculated in the fundamental HOH-bending, fundamental OH-stretching, first OH-stretching-HOH-bending combination, and first OH-stretching overtone (ΔvOH = 2) regions of the water dimer's spectrum. Furthermore, the rotational-vibrational spectrum was calculated in the ΔvOH = 2 region at 10 K, corresponding to the temperature of the existing jet-expansion experiments. The calculated spectrum was obtained by combining results from a full-dimensional (12D) vibrational and a reduced-dimensional vibrational-rotational-tunneling model. The ΔvOH = 2 spectral region is rich in features due to contributions from multiple vibrational-rotational-tunneling sub-bands. Origins of the experimental vibrational bands depend on the assignment of the observed sub-bands. Based on our calculations, we assign the observed sub-bands, and our reassignment leads to new values for the vibrational band origins of the free donor and antisymmetric acceptor OH-stretching first overtones of ∼7227 and ∼7238 cm-1, respectively. The observed bands with origins at 7192.34 and ∼7366 cm-1 are assigned to the symmetric acceptor OH-stretching first overtone and the OH-stretching combination of the donor, respectively.

H2O-HCN complex: A new potential energy surface and intermolecular rovibrational states from rigorous quantum calculations

In this work the H2O-HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O-HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O-HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O-HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former.

Water as a Lévy Rotor

A probability density function describing the angular evolution of a fixed-length atom-atom vector as a Lévy rotor is derived containing just two dynamical parameters: the Lévy parameter α and a rotational time constant τ. A Lévy parameter α<2 signals anomalous (non-Brownian) motion. Molecular dynamics simulation of water at 298 K validates the probability density function for the intramolecular ^{1}H─^{1}H dynamics. The rotational dynamics of water is found to be approximately Brownian at subpicosecond time intervals, becomes increasingly anomalous at longer time intervals due to hydrogen-bond breaking and reforming, before becoming indistinguishable from Brownian dynamics beyond about 25 ps. The Lévy rotor model is used to estimate the intramolecular contribution to the longitudinal nuclear-magnetic-resonance (NMR) relaxation rate R_{1,intra}. It is found that R_{1,intra} contributes 65%±7% to the overall relaxation rate of water at room temperature.

Three- and four-site models for heavy water: SPC/E-HW, TIP3P-HW, and TIP4P/2005-HW

Heavy water or deuterium oxide, D₂O, is used as a solvent in various biophysical and chemical experiments. To model such experiments with molecular dynamics simulations, effective pair potentials for heavy water are required, which reproduce the well-known physicochemical differences relative to light water. We present three effective pair potentials for heavy water, denoted SPC/E-HW, TIP3P-HW, and TIP4P/2005-HW. The models were parameterized by modifying the widely used three- and four-site models for light water, with the aim of maintaining the specific characteristics of the light water models. At room temperature, SPC/E-HW and TIP3P-HW capture the modulations relative to light water of the mass and electron densities, heat of vaporization, diffusion coefficient, and water structure. TIP4P/2005-HW captures, in addition, the density of heavy water over a wide temperature range.

High-Level Rovibrational Calculations on Ketenimine

From an astrochemical point of view ketenimine (CH₂CNH) is a complex organic molecule (COM) and therefore likely to be a building block for biologically relevant molecules. Since it has been detected in the star-forming region Sagittarius B2(N), it is of high relevance in this field. Although experimental data are available for certain bands, for some energy ranges such as above 1200 cm^-1 reliable data virtually do not exist. In addition, high-level ab initio calculations are neither reported for ketenimine nor for one of its deuterated isotopologues. In this paper, we provide for the first time data from accurate quantum chemical calculations and a thorough analysis of the full rovibrational spectrum. Based on high-level potential energy surfaces obtained from explicitly correlated coupled-cluster calculations including up to 4-mode coupling terms, the (ro)vibrational spectrum of ketenimine has been studied in detail by variational calculations relying on rovibrational configuration interaction (RVCI) theory. Strong Fermi resonances were found for all isotopologues. Rovibrational infrared intensities have been obtained from dipole moment surfaces determined from the distinguishable cluster approximation. A comparison of the spectra of the CH₂CNH molecule with experimental data validates our results, but also reveals new insight about the system, which shows very strong Coriolis coupling effects.

Potential energy surface and bound states of the H2O-HF complex

We present the first global five-dimensional potential energy surface for the H₂O-HF dimer, a prototypical hydrogen bonded complex. Large scale ab initio calculations were carried out using the explicitly correlated coupled cluster approach with single- and double-excitations together with non-iterative perturbative treatment of triple excitations with the augmented correlation-consistent triple zeta basis sets, in which the water and hydrogen fluoride monomers were frozen at their vibrationally averaged geometries. The ab initio data points were fitted to obtain a global potential energy surface for the complex. The equilibrium geometry of the complex corresponds to the formation of a hydrogen bond with water acting as a proton acceptor and a binding energy of D_e = 3059 cm^-1 (8.75 kcal/mol). The energies and wavefunctions of the lowest bound states of the complex were computed using a variational approach, and the dissociation energies of both ortho-H₂O-HF (D₀ = 2089.4 cm^-1 or 5.97 kcal/mol) and para-H₂O-HF (D₀ = 2079.6 cm^-1 or 5.95 kcal/mol) were obtained. The rotational constant of the complex was found to be in good agreement with the available experimental data.

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.

Toward Elimination of Discrepancies between Theory and Experiment: Anharmonic Rotational-Vibrational Spectrum of Water in Solid Noble Gas Matrices

Rotational–vibrational spectroscopy of water in solid noble gas matrices has been studied for many decades. Despite that, discrepancies persist in the literature about the assignment of specific bands. We tackle the involved rotational–vibrational spectrum of the water isotopologues H₂¹⁶O, HD¹⁶O, and D₂¹⁶O with an unprecedented combination of experimental high-resolution matrix isolation infrared (MI-IR) spectroscopy and computational anharmonic vibrational spectroscopy by vibrational configuration interaction (VCI) on high-level ab initio potential energy surfaces. With VCI, the average deviation to gas-phase experiments is reduced from >100 to ≈1 cm^–1 when compared to harmonic vibrational spectra. Discrepancies between MI-IR and VCI spectra are identified as matrix effects rather than missing anharmonicity in the theoretical approach. Matrix effects are small in Ne (≈1.5 cm^–1) and a bit larger in Ar (≈10 cm^–1). Controversial assignments in Ne MI-IR spectra are resolved, for example, concerning the ν₃ triad in HDO. We identify new transitions, for example, the ν₂ 1₀₁ ← 1₁₀ transition in D₂O and H₂O or the ν₃ 0₀₀ ← 1₀₁ transition in D₂O, and reassign bands, for example, the band at 3718.9 cm^–1 that is newly assigned as the 1₁₀ ← 1₁₁ transition. The identification and solution of discrepancies for a well-studied benchmark system such as water prove the importance of an iterative and one-hand combination of theory and experiment in the field of high-resolution infrared spectroscopy of single molecules. As the computational costs involved in the VCI approach are reasonably low, such combined experimental/theoretical studies can be extended to molecules larger than triatomics.

Soft experimental constraints for soft interactions: a spectroscopic benchmark data set for weak and strong hydrogen bonds

An experimental benchmark data base on rotational constants, vibrational properties and energy differences for weakly and more strongly hydrogen-bonded complexes and their constituents from the spectroscopic literature is assembled. It is characterized in detail and finally contracted to a more compact, discriminatory set (ENCH-51, for Experimental Non-Covalent Harmonic with 51 entries). The meeting points between theory and experiment consist of equilibrium rotational constants and harmonic frequencies and energies, which are back-corrected from experimental observables and are very easily accessible by quantum chemical calculations. The relative performance of B3LYP-D3, PBE0-D3 and M06-2X density functional theory predictions with a quadruple-zeta basis set is used to illustrate systematic errors, error compensation and selective performance for structural, vibrational and energetical observables. The current focus is on perspectives and different benchmarking methodologies, rather than on a specific theoretical method or a specific class of compounds. Extension of the data base in chemical, observable and quantum chemical method space is encouraged.

Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states

Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H₂O)₂, CH₄·H₂O, and (CH₄)₂.

Rovibrational quantum dynamics of the vinyl radical and its deuterated isotopologues

Rotational–vibrational states up to 3200 cm⁻¹, beyond the highest-lying stretching fundamental, are computed variationally for the vinyl radical (VR), H₂C_βC_αH, and the following deuterated isotopologues of VR: CH₂CD, CHDCH, and CD₂CD.

Fully quantum calculation of the second and third virial coefficients of water and its isotopologues from ab initio potentials

Path-Integral Monte Carlo methods were applied to calculate the second, B(T), and the third, C(T), virial coefficients for water. A fully quantum approach and state-of-the-art flexible-monomer pair and three-body potentials were used. Flexible-monomer potentials allow calculations for any isotopologue; we performed calculations for both H2O and D2O. For B(T) of H2O, the quantum effect contributes 25% of the value at 300 K and is not entirely negligible even at 1000 K, in accordance with recent literature findings. The effect of monomer flexibility, while not as large as some claims in the literature, is significant compared to the experimental uncertainty. It is of opposite sign to the quantum effect, smaller in magnitude than the latter below 500 K, and varies from 2% at 300 K to 10% at 700 K. When monomer flexibility is accounted for, results from the CCpol-8sf pair potential are in excellent agreement with the available experimental data and provide reliable B(T) values at temperatures outside the range of experimental data. The flexible-monomer MB-pol pair potential yields B(T) values that are slightly too high compared to experiment. For C(T), our calculations confirm earlier findings that the use of three-body potential is necessary for meaningful predictions. However, due to various uncertainties of the potentials used, especially the three-body ones, we were not able to establish benchmark values of C(T), although our results are in qualitative agreement with available experimental data. The quantum effect, never before included for water, reduces the magnitude of the classical value for H2O by a factor of 2.5 at 300 K and is not entirely negligible even at 1000 K.

Interaction of H2O with CO: potential energy surface, bound states and scattering calculations

We present the first scattering calculations for the H₂O–CO system based on a high accuracy potential energy surface.

Rovibrational quantum dynamical computations for deuterated isotopologues of the methane–water dimer

Rovibrational states of methane–water isotopologues are computed in a variational procedure and the wave functions are analyzed in terms of the rigid-rotor and coupled-rotors models.

Rovibrational transitions of the methane–water dimer from intermolecular quantum dynamical computations

Quantum dynamical computations of the methane–water dimer explain the far-infrared spectrum of this important prototype of a water–hydrocarbon interaction.

A systematic development of a polarizable potential of water

Based on extensive studies of existing potentials we propose a new molecular model for water. The new model is rigid and contains three Gaussian charges. Contrary to other models, all charges take part in the polarization of the molecule. They are connected by harmonic springs to their gas-phase positions: the negative charge to a prescribed point on the main axis of the molecule; the positive charges to the hydrogens. The mechanical equilibrium between the electrostatic forces and the spring forces determines the polarization of the molecule which is established by iteration at every timestep. The model gives excellent estimates for ambient liquid properties and reasonably good results from high-pressure solids to gas-phase clusters. We present a detailed description of the development of this model and a large number of calculated properties compared to the estimates of the nonpolarizable TIP4P∕2005 [J. L. F. Abascal and C. Vega, J. Chem. Phys. 123, 234505 (2005)], the polarizable GCPM [P. Paricaud, M. Predota, A. A. Chialvo, and P. T. Cummings, J. Chem. Phys. 122, 244511 (2005)], and our earlier BKd3 model [P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)]. The best overall performance is shown by the new model.