Spectra of water dimer from a new ab initio potential with flexible monomers

Infrared spectra of partially deuterated water dimers in the fundamental O-D stretch region

Spectra of mixed H/D water dimers are studied in the mid-infrared region of the O-D stretch fundamental (2630-2800 cm-1) using a pulsed supersonic slit jet and a tunable optical parametric oscillator infrared source. Over 30 bands, belonging to nine of the ten possible isotopologues (only H2O-HOD is missed), are observed and analyzed. The implications for excited state tunneling splittings, lifetime effects, and vibrational shifts are discussed. These are the first significant new experimental results on (gas phase) mixed water dimers in over 25 years, and they are valuable for testing water dimer theoretical calculations, a field which continues to be of lively current interest.

Tunneling splittings in the vibrationally excited states of water trimer

Tunneling splitting (TS) patterns in vibrationally excited states of the water trimer are calculated, taking into account six tunneling pathways that describe the flips of free OH bonds and five bifurcation mechanisms that break and reform hydrogen bonds in the trimer ring. A tunneling matrix (TM) model is used to derive the energy shifts due to tunneling in terms of the six distinct TM elements in symbolic form. TM elements are calculated using the recently-developed modified WKB (Wentzel-Kramers-Brillouin) method in full dimensionality. Convergence was achieved for the lowest six excited vibrational modes. Bifurcation widths of the pseudorotational quartets are shown to be of comparable size to the ground-state widths, obtained using instanton theory, or increased for some particular modes of vibration. The largest increase is obtained for the excited out-of-phase flip of two adjacent water monomers with free OH bonds pointing in opposite directions relative to the ring plane. Bifurcation widths in (D2O)3 are found to be two orders of magnitude smaller than in (H2O)3. Geometrical arguments were used to explain the order of states in some TS multiplets in vibrationally excited water trimers.

Spectra of the D2O dimer in the O-D fundamental stretch region: The acceptor symmetric stretch fundamental and new combination bands

The O-D stretch fundamental region of the deuterated water dimer, (D2O)2, is further studied using a pulsed supersonic slit jet and a tunable optical parametric oscillator infrared source. The previously unobserved acceptor symmetric O-D stretch fundamental vibration is detected, with Ka = 0 ← 0 and 1 ← 0 sub-bands at about 2669 and 2674 cm-1, respectively. The analysis indicates that the various water dimer tunneling splittings generally decrease in the excited vibrational state, similar to the three other previously observed O-D stretch fundamentals. Two new (D2O)2 combination bands are observed, giving information on intermolecular vibrations in the excited O-D stretch states. The likely vibrational assignments for these and a previously observed combination band are discussed.

Microhydration of Tertiary Amines: Robust Resonances in Red-Shifted Water

Tertiary amines are strong hydrogen bond acceptors. When a water molecule donates one of the OH groups, its in-phase stretching vibration wavenumber is decreased to such an extent that it comes close to the water bending overtone. This gives rise to anharmonic phenomena such as classical Fermi resonance, resonance with multiple-quantum dark states, or combination transitions with low-frequency intermolecular modes. These effects, which contribute to the spectral breadth of room-temperature hydrogen-bonded amine complexes, are disentangled by Fourier transform infrared spectroscopy in pulsed supersonic slit jet expansions. Monohydrates of the amines quinuclidine, N-methylpyrrolidine, N-methylpiperidine, and dimethylcyclohexylamine exhibit systematic mode coupling signatures. These suggest relatively fast energy flow out of the excited OH stretching fundamental into intra- and intermolecular degrees of freedom of the hydrogen-bonded water molecule. Trimeric complexes are spectroscopically separated from the amine monohydrates.

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.

Computing excited OH stretch states of water dimer in 12D using contracted intermolecular and intramolecular basis functions

Due to the ubiquity and importance of water, water dimer has been intensively studied. Computing the (ro-)vibrational spectrum of water dimer is challenging. The potential has eight wells separated by low barriers, which makes harmonic approximations of limited utility. A variational approach is imperative, but difficult because there are 12 coupled vibrational coordinates. In this paper, we use a product contracted basis whose functions are products of intramolecular and intermolecular functions computed using an iterative eigensolver. An intermediate matrix F facilitates calculating matrix elements. Using F, it is possible to do calculations on a general potential without storing the potential on the full quadrature grid. We find that surprisingly many intermolecular functions are required. This is due to the importance of coupling between inter- and intra-molecular coordinates. The full G₁₆ symmetry of water dimer is exploited. We calculate, for the first time, monomer excited stretch states and compare P(1) transition frequencies with their experimental counterparts. We also compare with experimental vibrational shifts and tunneling splittings. Surprisingly, we find that the largest tunneling splitting, which does not involve the interchange of the two monomers, is smaller in the asymmetric stretch excited state than in the ground state. Differences between levels we compute and those obtained with a [6+6]D adiabatic approximation [Leforestier et al. J. Chem. Phys. 137 014305 (2012)] are ∼0.6 cm^-1 for states without monomer excitation, ∼4 cm^-1 for monomer excited bend states, and as large as ∼10 cm^-1 for monomer excited stretch states.

Δ-Machine Learned Potential Energy Surfaces and Force Fields

There has been great progress in developing machine-learned potential energy surfaces (PESs) for molecules and clusters with more than 10 atoms. Unfortunately, this number of atoms generally limits the level of electronic structure theory to less than the "gold standard" CCSD(T) level. Indeed, for the well-known MD17 dataset for molecules with 9-20 atoms, all of the energies and forces were obtained with DFT calculations (PBE). This Perspective is focused on a Δ-machine learning method that we recently proposed and applied to bring DFT-based PESs to close to CCSD(T) accuracy. This is demonstrated for hydronium, N-methylacetamide, acetyl acetone, and ethanol. For 15-atom tropolone, it appears that special approaches (e.g., molecular tailoring, local CCSD(T)) are needed to obtain the CCSD(T) energies. A new aspect of this approach is the extension of Δ-machine learning to force fields. The approach is based on many-body corrections to polarizable force field potentials. This is examined in detail using the TTM2.1 water potential. The corrections make use of our recent CCSD(T) datasets for 2-b, 3-b, and 4-b interactions for water. These datasets were used to develop a new fully ab initio potential for water, termed q-AQUA.

Vibrational Tunneling Spectra of Molecules with Asymmetric Wells: A Combined Vibrational Configuration Interaction and Instanton Approach

A combined approach that uses the vibrational configuration interaction (VCI) and semiclassical instanton theory was developed to study vibrational tunneling spectra of molecules with multiple wells in full dimensionality. The method can be applied to calculate low-lying vibrational states in the systems with an arbitrary number of minima, which are not necessarily equal in energy or shape. It was tested on a two-dimensional double-well model system and on malonaldehyde, and the calculations reproduced the exact quantum mechanical (QM) results with high accuracy. The method was subsequently applied to calculate the vibrational spectrum of the asymmetrically deuterated malonaldehyde with nondegenerate vibrational frequencies in the two wells. The spectrum is obtained at a cost of single-well VCI calculations used to calculate the local energies. The interactions between states of different wells are computed semiclassically using the instanton theory at a comparatively negligible computational cost. The method is particularly suited to systems in which the wells are separated by large potential barriers and tunneling splittings are small, for example, in some water clusters, when the exact QM methods come at a prohibitive computational cost.

Reduced-dimensional vibrational models of the water dimer

A model based on the finite-basis representation of a vibrational Hamiltonian expressed in internal coordinates is developed. The model relies on a many-mode, low-order expansion of both the kinetic energy operator and the potential energy surface (PES). Polyad truncations and energy ceilings are used to control the size of the vibrational basis to facilitate accurate computations of the OH stretch and HOH bend intramolecular transitions of the water dimer (H₂ ¹⁶O)₂. Advantages and potential pitfalls of the applied approximations are highlighted. The importance of choices related to the treatment of the kinetic energy operator in reduced-dimensional calculations and the accuracy of different water dimer PESs are discussed. A range of different reduced-dimensional computations are performed to investigate the wavenumber shifts in the intramolecular transitions caused by the coupling between the intra- and intermolecular modes. With the use of symmetry, full 12-dimensional vibrational energy levels of the water dimer are calculated, predicting accurately the experimentally observed intramolecular fundamentals. It is found that one can also predict accurate intramolecular transition wavenumbers for the water dimer by combining a set of computationally inexpensive reduced-dimensional calculations, thereby guiding future effective-Hamiltonian treatments.

Vibrational Spectroscopy of the Water Dimer at Jet-Cooled and Atmospheric Temperatures

The vibrational spectroscopy of the water dimer provides an understanding of basic hydrogen bonding in water clusters, and with about one water dimer for every 1,000 water molecules, it plays a critical role in atmospheric science. Here, we review how the experimental and theoretical progress of the past decades has improved our understanding of water dimer vibrational spectroscopy under both cold and warm conditions. We focus on the intramolecular OH-stretching transitions of the donor unit, because these are the ones mostly affected by dimer formation and because their assignment has proven a challenge. We review cold experimental results from early matrix isolation to recent mass-selected jet expansion techniques and, in parallel, the improvements in the theoretical anharmonic models. We discuss and illustrate changes in the vibrational spectra of complexes upon increasing temperature, and the difficulties in recording and calculating these spectra. In the atmosphere, water dimer spectra at ambient temperature are crucial. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

GPU-Accelerated Neural Network Potential Energy Surfaces for Diffusion Monte Carlo

Diffusion Monte Carlo (DMC) provides a powerful method for understanding the vibrational landscape of molecules that are not well-described by conventional methods. The most computationally demanding step of these calculations is the evaluation of the potential energy. In this work, a general approach is developed in which a neural network potential energy surface is trained by using data generated from a small-scale DMC calculation. Once trained, the neural network can be evaluated by using highly parallelizable calls to a graphics processing unit (GPU). The power of this approach is demonstrated for DMC simulations on H₂O, CH₅⁺, and (H₂O)₂. The need to include permutation symmetry in the neural network potentials is explored and incorporated into the molecular descriptors of CH₅⁺ and (H₂O)₂. It is shown that the zero-point energies and wave functions obtained by using the neural network potentials are nearly identical to the results obtained when using the potential energy surfaces that were used to train the neural networks at a substantial savings in the computational requirements of the simulations.

Are "Bright-State" Models Appropriate for Analyzing Fermi-Coupled Bands in Molecular Vibrational Spectra?

Bright-state models are often applied to "deperturb" Fermi-coupled bands in molecular vibrational spectra, in cases where a harmonically forbidden transition "borrows" intensity from an energetically nearby allowed transition. However, forbidden transitions can also acquire intensity through anharmonic couplings on the potential energy surface ("mechanical anharmonicity") or dipole moment surface ("electrical anharmonicity") that are not accounted for within the bright-state model. In this work, we compare deperturbation shifts obtained by analysis of experimental data with those predicted using the bright-state model, for a series of discrete encapsulated chloride hydrate isotopomers. Predicted band center shifts and Fermi coupling matrix elements obtained using the bright-state model are larger than those estimated from experimental data.

Tunneling splittings of vibrationally excited states using general instanton paths

A multidimensional semiclassical method for calculating tunneling splittings in vibrationally excited states of molecules using Cartesian coordinates is developed. It is an extension of the theory by Mil'nikov and Nakamura [J. Chem. Phys. 122, 124311 (2005)] to asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. Additionally, new terms are introduced in the description of the semiclassical wavefunction that drastically improves the splitting estimates for certain systems. The method is based on the instanton theory and builds the semiclassical wavefunction of the vibrationally excited states from the ground-state instanton wavefunction along the minimum action path and its harmonic neighborhood. The splittings of excited states are thus obtained at a negligible added numerical effort. The cost is concentrated, as for the ground-state splittings, in the instanton path optimization and the hessian evaluation along the path. The method can thus be applied without modification to many mid-sized molecules in full dimensionality and in combination with on-the-fly evaluation of electronic potentials. The tests were performed on several model potentials and on the water dimer.

Computational analysis of vibrational frequencies and rovibrational spectroscopic constants of hydrogen sulfide dimer using MP2 and CCSD(T)

Previous studies have shown that the weakly bonded H₂S dimer demands high level quantum chemical calculations to reproduce experimental values. We investigated the hydrogen bonding of H₂S dimer using MP2 and CCSD(T) levels of theory in combination with aug-cc-pV(D,T,Q)Z basis sets. More precisely, the binding energies, potential energy curves, rovibrational spectroscopic constants, decomposition lifetime, and normal vibrational frequencies were calculated. In addition, we introduced the local mode analysis of Konkoli-Cremer to quantify the hydrogen bonding in the H₂S dimer as well as providing for the first time the comprehensive decomposition of normal vibrational modes into local modes contributions, and a decomposition lifetime based on rate constant. The local mode force constant of the H₂S dimer hydrogen bond is smaller than that of the water dimer, in accordance with the weaker hydrogen bonding in the H₂S dimer.

Explicitly correlated ab initio potential energy surface and predicted rovibrational spectra for H2O-N2 and D2O-N2 complexes

An ab initio intermolecular potential energy surface (PES) for the van der Waals complex of H₂O-N₂ that explicitly incorporates the intramolecular Q₂ bending normal mode of the H₂O monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster theory [CCSD(T)-F12] with an augmented correlation-consistent triple zeta basis set and an additional bond function. Analytic five-dimensional intermolecular PESs for ν₂(H₂O) = 0 and 1 are obtained by fitting to the multi-dimensional Morse/long-range potential function form. These fits to 40 890 points have the root-mean-square (rms) discrepancy of 0.88 cm^-1 for interaction energies less than 2000.0 cm^-1. The resulting vibrationally averaged PESs provide good representations of the experimental microwave and infrared data: for microwave transitions of H₂O-N₂, the rms discrepancy is only 0.0003 cm^-1, and for infrared transitions of the A₁ symmetry of the H₂O(ν₂ = 1 ← 0)-N₂, the rms discrepancy is 0.001 cm^-1. The calculated infrared band origin shifts associated with the ν₂ bending vibration of water are 2.210 cm^-1 and 1.323 cm^-1 for H₂O-N₂ and D₂O-N₂, respectively, in good agreement with the experimental values of 2.254 cm^-1 and 1.266 cm^-1. The benchmark tests and comparisons of the predicted spectral properties are carried out between CCSD(T)-F12a and CCSD(T)-F12b approaches.

Observation of the Low-Frequency Spectrum of the Water Trimer as a Sensitive Test of the Water-Trimer Potential and the Dipole-Moment Surface

Intermolecular interactions in bulk water are dominated by pairwise and non‐pairwise cooperative interactions. While accurate descriptions of the pairwise interactions are available and can be tested by precise low‐frequency spectra of the water dimer up to 550 cm⁻¹, the same does not hold for the three‐body interactions. Here, we report the first comprehensive spectrum of the water trimer in the frequency region from 70 to 620 cm⁻¹ using helium‐nanodroplet isolation and free‐electron lasers. By comparison to accompanying high‐level quantum calculations, the experimentally observed intermolecular bands can be assigned. The transition frequencies of the degenerate translation, the degenerate in‐plane and the non‐degenerate out‐of‐plane libration, as well as additional bands of the out‐of‐plane librational mode are reported for the first time. These provide a benchmark for state‐of‐the‐art water potentials and dipole‐moment surfaces, especially with respect to three‐body interactions.

Automatic Generation of Flexible-Monomer Intermolecular Potential Energy Surfaces

A method is developed for automatic generation of nonreactive intermolecular two-body potential energy surfaces (PESs) including intramonomer degrees of freedom. This method, called flex-autoPES, is an extension of the autoPES method developed earlier, which assumes rigid monomers. In both cases, the whole PES development proceeds without any human intervention. The functional form used is a sum of products of site-site functions (both atomic and off-atomic sites can be used). The leading terms with sites involving different monomers are of physically motivated form. The long-range part of a PES is computed from monomer properties without using any dimer information. The close-range part is fitted to dimer interaction energies computed using electronic structure methods. Virtually any method can be used in such calculations, but the use of symmetry-adapted perturbation theory provides a seamless connection to the long-range part of the PES. The performance of the flex-autoPES code was tested by developing a full-dimensional PES for the water dimer and PESs including only some soft intramonomer degrees of freedom for the ethylene glycol dimer and for the ethylene glycol-water dimer. In the case of the water dimer, the root-mean-square error (RMSE) of the PES from the data points with negative total energies is 0.03 kcal/mol, and we expect this PES to be more accurate than any previously published PES of this type. For the ethylene glycol dimer and the ethylene glycol-water dimers, the analogous RMSEs are 0.25 and 0.1 kcal/mol, respectively.

Instanton theory of ground-state tunneling splittings with general paths

We derive a multidimensional instanton theory for calculating ground-state tunneling splittings in Cartesian coordinates for general paths. It is an extension of the method by Mil'nikov and Nakamura [J. Chem. Phys. 115, 6881 (2001)] to include asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. The approach avoids multiple expensive matrix diagonalizations to converge the fluctuation prefactor in the ring-polymer instanton (RPI) method, and instead replaces them by an integration of a Riccati differential equation. When combined with the string method for locating instantons, we avoid the need to converge the calculation with respect to the imaginary time period of the semiclassical orbit, thereby reducing the number of convergence parameters of the optimized object to just one: the number of equally spaced system replicas used to represent the instanton path. The entirety of the numerical effort is thus concentrated in optimizing the shape of the path and evaluating hessians along the path, which is a dramatic improvement over RPI. In addition to the standard instanton approximations, we neglect the coupling of vibrational modes to external rotations. The method is tested on the model potential of malonaldehyde and on the water dimer and trimer, giving close agreement with RPI at a much-reduced cost.

Infrared Spectroscopy of Neutral Water Dimer Based on a Tunable Vacuum Ultraviolet Free Electron Laser

Infrared (IR) spectroscopy provides detailed structural and dynamical information on clusters at the fingerprint level. Herein, we demonstrate the capability of a tunable vacuum ultraviolet free electron laser (VUV-FEL) for selective detection of a wide variety of neutral water clusters and for recording the size-dependent IR spectra. The present technique does not require the presence of an ultraviolet chromophore or a dipole moment and is generally applicable for IR spectroscopy of neutral clusters free from confinement. To show the features of our technique, we report here the IR spectra of neutral water dimer in the OH stretch region, providing benchmarks for theoretical study of the accurate description of hydrogen bonding structures involved in liquid water and ice. Quantum mechanical calculations on a 12-dimensional ab initio potential energy surface are utilized to simulate the anharmonic vibrational spectra of water dimer. These results help to resolve the controversy of the exact vibrational assignment of each band feature of the water dimer.

Solvent- and concentration-induced self-assembly of an amphiphilic perylene dye

PDA-CA self-assembly behavior depends on temperature, concentration, and solvent.

Observation of the Low-Frequency Spectrum of the Water Dimer as a Sensitive Test of the Water Dimer Potential and Dipole Moment Surfaces

Using the helium nanodroplet isolation setup at the ultrabright free-electron laser source FELIX in Nijmegen (BoHeNDI@FELIX), the intermolecular modes of water dimer in the frequency region from 70 to 550 cm^-1 were recorded. Observed bands were assigned to donor torsion, acceptor wag, acceptor twist, intermolecular stretch, donor torsion overtone, and in-plane and out-of-plane librational modes. This experimental data set provides a sensitive test for state-of-the-art water potentials and dipole moment surfaces. Theoretical calculations of the IR spectrum are presented using high-level quantum and approximate quasiclassical molecular dynamics approaches. These calculations use the full-dimensional ab initio WHHB potential and dipole moment surfaces. Based on the experimental data, a considerable increase of the acceptor switch and a bifurcation tunneling splitting in the librational mode is deduced, which is a consequence of the effective decrease in the tunneling barrier.

Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements

The rovibrational energy levels and intensities of the CH₃F-H₂ dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H₂O-H₂ dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H₂ are more strongly bound than dimers containing para-H₂. Energies and wavefunctions of the discrete rovibrational levels of CH₃F-paraH₂ complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH₃F-paraH₂ complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH₃F-orthoH₂ dimer has been predicted for the first time.

Toward breaking the curse of dimensionality in (ro)vibrational computations of molecular systems with multiple large-amplitude motions

Methodological progress is reported in the challenging direction of a black-box-type variational solution of the (ro)vibrational Schrödinger equation applicable to floppy, polyatomic systems with multiple large-amplitude motions. This progress is achieved through the combination of (i) the numerical kinetic-energy operator (KEO) approach of Mátyus et al. [J. Chem. Phys. 130, 134112 (2009)] and (ii) the Smolyak nonproduct grid method of Avila and Carrington, Jr. [J. Chem. Phys. 131, 174103 (2009)]. The numerical representation of the KEO makes it possible to choose internal coordinates and a body-fixed frame best suited for the molecular system. The Smolyak scheme reduces the size of the direct-product grid representation by orders of magnitude, while retaining some of the useful features of it. As a result, multidimensional (ro)vibrational states are computed with system-adapted coordinates, a compact basis- and grid-representation, and an iterative eigensolver. Details of the methodological developments and the first numerical applications are presented for the CH₄·Ar complex treated in full (12D) vibrational dimensionality.

Spectra of the D2O dimer in the O-D fundamental stretch region: Vibrational dependence of tunneling splittings and lifetimes

The fundamental O-D stretch region (2600-2800 cm^-1) of the fully deuterated water dimer (D₂O)₂ is studied using a pulsed supersonic slit jet source and a tunable optical parametric oscillator source. Relatively high spectral resolution (0.002 cm^-1) enables all six dimer tunneling components to be observed, in most cases, for the acceptor asymmetric O-D stretch, the donor free O-D stretch, and the donor bound O-D stretch vibrations. The dominant acceptor switching tunneling splittings are observed to decrease moderately in the excited O-D stretch states, to roughly 75% of their ground state values, whereas the smaller donor-acceptor interchange splittings show more dramatic and irregular decreases. Excited state predissociation lifetimes, as determined from the observed line broadening, show large variations (0.2 ≤ τ ≤ 5 ns) depending on the vibrational state, K-value, and tunneling symmetry. Another very weak band is tentatively assigned to a combination mode involving an intramolecular O-D stretch plus an intermolecular twist overtone. Asymmetric O-D stretch bands of the mixed isotopologue dimers D₂O-DOH and D₂O-HOD are also observed and analyzed.

Computing vibration-rotation-tunnelling levels of HOD dimer

Using an accurate 6D water dimer potential energy surface, we compute vibration-rotation-tunnelling levels of HOD dimer, by assuming that the two monomers are rigid. HOD dimer has two isomers, a D-bonded isomer and an H-bonded isomer, and the wavefunctions of both isomers have amplitude in four wells. HOD dimer is important because, unlike the case of H2O dimer or D2O dimer, it is possible to measure the largest tunnelling splitting. Results for HOD dimer, therefore facilitate the testing of H2O dimer potentials. In J. Chem. Phys., 1995, 102, 1114, experimental results were interpreted in terms of 1D models. Experimental splittings of both isomers, obtained by fitting an energy level equation to spectra, are in good agreement with those we compute.

The water–carbon monoxide dimer: new infrared spectra, ab initio rovibrational energy level calculations, and an interesting in-termolecular mode

Bound state rovibrational energy level calculations using a high-level intermolecular potential surface are reported for H₂O–CO and D₂O–CO.

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₄)₂.

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.

Designer Metallic Acceptor-Containing Halogen Bonds: General Strategies

Being electrostatic interactions in nature, hydrogen bonding (HB) and halogen bonding (XB) are considered to be two parallel worlds. In principle, all the applications that HB have could also be applied to XB. However, there has been no report on an anionic, metallic XB acceptor, but metal anions have been observed to be good HB acceptors. This missing mosaic piece of XB theory is because common metal anions are reactive for XB donors. In view of this, two strategies are proposed for designing metallic acceptor-containing XB using ab initio calculations. The first one is to utilize a metal cluster anion with a high electron detachment energy, such as the superatom, Al₁₃^- as the XB acceptor. The second strategy is to design a ligand-passivated/protected metal core that can maintain the negative charge; several exotic clusters, such as PtH₅^- , PtZnH₅^- , and PtMgH₅^- , are used as examples. Based on these two strategies, it is anticipated that more metallic acceptor-containing XBs will be discovered.

Single Mutations Reshape the Structural Correlation Network of the DMXAA-Human STING Complex

Subtle changes in protein sequences are able to alter ligand-protein interactions. Unraveling the mechanism of such phenomena is important for understanding ligand-protein interactions, including the DMXAA-STING interaction. DMXAA specifically binds to mouse STING instead of human STING. However, the S162A mutation and a newly discovered E260I mutation endow human STING^AQ with DMXAA sensitivity. Through molecular dynamics simulations, we revealed how these single mutations alter the DMXAA-STING interaction. Compared to mutated systems, structural correlations in the interaction of STING^AQ with DMXAA are stronger, and the correlations are cross-protomers in the dimeric protein. Analyses on correlation coefficients lead to the identification of two key interactions that mediate the strong cross-protomer correlation in the DMXAA-STING^AQ interaction network: DMXAA-267T-162S* and 238R-260E*. These two interactions are partially and totally interrupted by the S162A and E260I mutations, respectively. Moreover, a smaller number of water molecules are displaced upon DMXAA binding to STING^AQ than that on binding to its mutants, leading to a larger entropic penalty for the former. Considering the sensitivity of STING^AQ and two of its mutants to DMXAA, a strong structural correlation appears to discourage DMXAA-STING binding. Such an observation suggests that DMXAA derivatives, which are deprived of hydrogen-bond interaction with both 162S* and 267T, are potential agonists of human STING.

All-dimensional H2-CO potential: Validation with fully quantum second virial coefficients

We use a new high-accuracy all-dimensional potential to compute the cross second virial coefficient B₁₂(T) between molecular hydrogen and carbon monoxide. The path-integral method is used to fully account for quantum effects. Values are calculated from 10 K to 2000 K and the uncertainty of the potential is propagated into uncertainties of B₁₂. Our calculated B₁₂(T) are in excellent agreement with most of the limited experimental data available, but cover a much wider range of temperatures and have lower uncertainties. Similar to recently reported findings from scattering calculations, we find that the reduced-dimensionality potential obtained by averaging over the rovibrational motion of the monomers gives results that are a good approximation to those obtained when flexibility is fully taken into account. Also, the four-dimensional approximation with monomers taken at their vibrationally averaged bond lengths works well. This finding is important, since full-dimensional potentials are difficult to develop even for triatomic monomers and are not currently possible to obtain for larger molecules. Likewise, most types of accurate quantum mechanical calculations, e.g., spectral or scattering, are severely limited in the number of dimensions that can be handled.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

On the importance of full-dimensionality in low-energy molecular scattering calculations

Scattering of H₂ on CO is of great importance in astrophysics and also is a benchmark system for comparing theory to experiment. We present here a new 6-dimensional potential energy surface for the ground electronic state of H₂-CO with an estimated uncertainty of about 0.6 cm⁻¹ in the global minimum region, several times smaller than achieved earlier. This potential has been used in nearly exact 6-dimensional quantum scattering calculations to compute state-to-state cross-sections measured in low-energy crossed-beam experiments. Excellent agreement between theory and experiment has been achieved in all cases. We also show that the fully 6-dimensional approach is not needed with the current accuracy of experimental data since an equally good agreement with experiment was obtained using only a 4-dimensional treatment, which validates the rigid-rotor approach widely used in scattering calculations. This finding, which disagrees with some literature statements, is important since for larger systems full-dimensional scattering calculations are currently not possible.