Ab initio intermolecular potential energy surfaces of the water-rare gas atom complexes

An intramolecular vibrationally excited intermolecular potential energy surface and predicted 2OH overtone spectroscopy of H2O-Kr

A new five-dimensional potential energy surface (PES) for H2O-Kr which explicitly includes the intramolecular 2OH overtone state of the H2O monomer is presented. The intermolecular potential energies were evaluated using explicitly correlated coupled cluster theory [CCSD(T)-F12] with a large basis set. Four vibrationally averaged analytical intermolecular PESs for H2O-Kr with H2O molecules in its |00+〉, |02+〉, |02-〉, and |11+〉 states are obtained by fitting to the multi-dimensional Morse/Long-Range potential function form. Each vibrationally averaged PES fitted to 578 points has root-mean-square (RMS) deviations smaller than 0.14 cm-1 and requires only 58 parameters. The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to calculate the rovibrational energy levels for |00+〉, |02+〉, |02-〉, and |11+〉 states of the H2O-Kr complexes. The calculated |02-〉Πf/e(101) ← |00+〉Σe(000) and |02+〉Πf/e(110) ← |00+〉Σe(101) infrared transitions are in excellent agreement with the experimental values with RMS discrepancies being only 0.007 and 0.016 cm-1, respectively. These analytical PESs can be used to provide reliable theoretical guidance for future infrared overtone spectroscopy of H2O-Kr.

Computational investigations of stable multiple-cage-occupancy He clathrate-like hydrostructures

One of the several possibilities offered by the interesting clathrate hydrates is the opportunity to encapsulate several atoms or molecules, in such a way that more efficient storage materials could be explored or new molecules that otherwise do not exist could be created. These types of applications are receiving growing attention from technologists and chemists, given the future positive implications that they entail. In this context, we investigated the multiple cage occupancy of helium clathrate hydrates, to establish stable novel hydrate structures or ones similar to those predicted previously by experimental and theoretical studies. To this purpose, we analyzed the feasibility of including an increased number of He atoms inside the small (D) and large (H) cages of the sII structure through first-principles properly assessed density functional approaches. On the one hand, we have computed energetic and structural properties, in which we examined the guest-host and guest-guest interactions in both individual and two-adjacent clathrate-like sII cages by means of binding and evaporation energies. On the other hand, we have carried out a thermodynamical analysis on the stability of such He-containing hydrostructures in terms of changes in enthalpy, ΔH, Gibbs free energy, ΔG, and entropy, ΔS, during their formation process at various temperature and pressure values. In this way, we have been able to make a comparison with experiments, reaffirming the ability of computational DFT approaches to describe such weak guest-host interactions. In principle, the most stable structure involves the encapsulation of one and four He atoms inside the D and H sII cages, respectively; however, more He atoms could be entrapped under lower temperature and/or higher pressure thermodynamic conditions. We foresee such accurate computational quantum chemistry approaches contributing to the current emerging machine-learning model development.

Vibrationally excited intermolecular potential energy surfaces and the predicted near infrared overtone (vOH = 2 ← 0) spectra of a H2O-Ne complex

The ab initio intra- and inter-molecular potential energy surfaces (PESs) for the H₂O-Ne system that explicitly incorporate the intramolecular overtone state (v_OH = 2) of H₂O are presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster theory [CCSD(T)-F12] level with an augmented correlation-consistent triple zeta basis set and an additional bond function. The vibrationally averaged three-dimensional intermolecular potentials for |00⁺〉, |02⁺〉, |02^-〉 and |11⁺〉 are obtained analytically by fitting to the multi-dimensional Morse/Long-range potential function form. These fits to 46 980 points have a root-mean-square (RMS) discrepancy of 0.12 cm^-1 for interaction energies less than 1000.0 cm^-1. With the vibrationally averaged PESs for the H₂O-Ne, we employed the combined radial discrete variable representation/angular finite basis representation method and Lanczos algorithm to calculate rovibrational energy levels (J = 0-10, n_s ≤ 2). The predicted infrared transitions and intensities of the para- and ortho-H₂O-Ne complex are in good agreement with the available experimental data for |02^-〉 ← |00⁺〉, |02⁺〉 ← |00⁺〉 transitions. In particular, the RMS discrepancy for |02^-〉∑^e(0₀₀,0) ← |00⁺〉∑^e(1₀₁,0), including P and R branch patterns, is only 0.045 cm^-1, which is comparable with the experimental values. These results will provide reliable theoretical guidance for the future infrared overtone spectroscopy of clusters.

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.

He Inclusion in Ice-like and Clathrate-like Frameworks: A Benchmark Quantum Chemistry Study of Guest-Host Interactions

Energetics and structural properties of selected type and size He@hydrate frameworks, e.g., from regular structured ice channels to clathrate-like cages, are presented from first-principles quantum chemistry methods. The scarcity of information on He@hydrates makes such complexes challenging targets, while their computational study entails an interesting and arduous task. Some of them have been synthesized in the laboratory, which motivates further investigations on their stability. Hence, the main focus is to examine the performance and accuracy of different wave function-based electronic structure methods, such as MP2, CCSD(T), their explicitly correlated (F12) and domain-based local pair-natural orbital (DLPNO) analogs, as well as modern and conventional density functional theory (DFT) approaches, and analytical model potentials available. Different structures are considered, starting from the "simplest system" formed by a noble gas atom (such as He) and one water molecule, followed by the study of the "fundamental units" present in all ice-like and clathrate-like frameworks (such as pentamers and hexamers) and finally the description of interactions in the "building blocks" of three-dimensional (3D) ice channels (e.g., horizontal and perpendicular ice II and I_h) and clathrate-like cages, such as the 5¹² present in the most common sI, sII, and sH clathrate-hydrate structures. The idea is to provide well-converged DLPNO-CCSD(T) and DFMP2/CBS reference datasets that in turn are used to validate how DFT functionals (in total, 29 approaches from generalized-gradient approximation (GGA), meta-GGA, to hybrid and range-separated functionals, including dispersion correction treatments, were checked) and analytical semiempirical/ab initio-based potentials perform compared with high-level alternatives. Within all tested approaches, those best-performing were identified and classified. Most of the DFT/DFT-D functionals, as well as available analytical pairwise model potentials, face difficulties in describing both hydrogen-bonded water frameworks and dispersion bound He-water interactions. Including dispersion corrections yields an overall well-balanced performance for LCωPBE-D3BJ and PBE0-D4 functionals. Such benchmark datasets can benefit research into the development of new cheminformatics models, as can serve to guide and cross-check methodologies, lending increased predicted power to future molecular simulations for investigating the role of structures and phase transitions from nanoscale clusters to macroscopic crystalline structures.

Nuclear spin/parity dependent spectroscopy and predissociation dynamics in vOH = 2 ← 0 overtone excited Ne-H2O clusters: Theory and experiment

Vibrationally state selective overtone spectroscopy and state- and nuclear spin-dependent predissociation dynamics of weakly bound ortho- and para-Ne-H₂O complexes (D_0(ortho) = 34.66 cm^-1 and D_0(para) = 31.67 cm^-1) are reported, based on near-infrared excitation of van der Waals cluster bands correlating with v_OH = 2 ← 0 overtone transitions (|02^-〉 and |02⁺〉) out of the ortho (1₀₁) and para (0₀₀) internal rotor states of the H₂O moiety. Quantum theoretical calculations for nuclear motion on a high level potential energy surface [CCSD(T)/VnZf12 (n = 3, 4)], corrected for basis set superposition error and extrapolated to the complete basis set (CBS) limit, are employed to successfully predict and assign Π-Σ, Σ-Σ, and Σ-Π infrared bands in the spectra, where Σ or Π represent approximate projections of the body-fixed H₂O angular momentum along the Ne-H₂O internuclear axis. IR-UV pump-probe experimental capabilities permit real-time measurements of the vibrational predissociation dynamics, which indicate facile intramolecular vibrational energy transfer from the H₂O v_OH = 2 overtone vibrations into the VdWs (van der Waals) dissociation coordinate on the τ_prediss = 15-25 ns time scale. Whereas all predicted strong transitions in the ortho-Ne-H₂O complexes are readily detected and assigned, vibrationally mediated photolysis spectra for the corresponding para-Ne-H₂O bands are surprisingly absent despite ab initio predictions of Q-branch intensities with S/N > 20-40. Such behavior signals the presence of highly selective nuclear spin ortho-para predissociation dynamics in the upper state, for which we offer a simple mechanism based on Ne-atom mediated intramolecular vibrational relaxation in the H₂O subunit (i.e., |02^±〉 → {|01^±〉; v₂ = 2}), which is confirmed by the ab initio energy level predictions and the nascent OH rotational (N), spin orbit (Π_1/2,_3/2), and lambda doublet product distributions.

Experimental and theoretical investigations of H2O-Ar

We have used continuous-wave cavity ring-down spectroscopy to record the spectrum of H₂O-Ar in the 2OH excitation range of H₂O. 24 sub-bands have been observed. Their rotational structure (T_rot = 12 K) is analyzed and the lines are fitted separately for ortho and para species together with microwave and far infrared data from the literature, with a unitless standard deviation σ=0.98 and 1.31, respectively. Their vibrational analysis is supported by a theoretical input based on an intramolecular potential energy surface obtained through ab initio calculations and computation of the rotational energy of sub-states of the complex with the water monomer in excited vibrational states up to the first hexad. For the ground and (010) vibrational states, the theoretical results agree well with experimental energies and rotational constants in the literature. For the excited vibrational states of the first hexad, they guided the assignment of the observed sub-bands. The upper state vibrational predissociation lifetime is estimated to be 3 ns from observed spectral linewidths.

Confinement Effects on the Nuclear Spin Isomer Conversion of H2O

The mechanism for interconversion between the nuclear spin isomers (NSI) of H₂O remains shrouded in uncertainties. The temperature dependence displayed by NSI interconversion rates for H₂O isolated in an argon matrix provides evidence that confinement effects are responsible for the dramatic increase in their kinetics with respect to the gas phase, providing new pathways for o-H₂O↔p-H₂O conversion in endohedral compounds. This reveals intramolecular aspects of the interconversion mechanism which may improve methodologies for the separation and storage of NSI en route to applications ranging from magnetic resonance spectroscopy and imaging to interpretations of spin temperatures in the interstellar medium.

Impact of metal ion's charge on the interatomic Coulombic decay widths in microsolvated clusters

Interatomic Coulombic decay (ICD) is an efficient electronic decay mechanism of electronically excited atoms and molecules embedded in an environment. For the series of isoelectronic Na(+), Mg(2+), and Al(3+) ions in aqueous solution, ultrashort ICD lifetimes of 3.1 fs, 1.5 fs, and 0.9 fs, respectively, were observed experimentally. The magnitude of the ICD lifetimes and their variation within the series were qualitatively explained by shortening metal-oxygen equilibrium distances and the increasing polarization of the water molecules as the metal's charge grows. We carried out an extensive ab initio investigation of the variation of the ICD widths with the metal-oxygen distances and the number of water neighbors in Na(+)-(H2O)m (m = 1-4) and Mg(2+)-(H2O)n (n = 1-6) clusters including and excluding polarization effects in the decaying state. We demonstrated that the effect of the induced polarization of the water ligand and the equilibrium cation-oxygen distance are equally important in determining the ordering and ratios of the ICD lifetimes in the series. Moreover, we showed that the induced polarization of the water molecules leads to a slower than linear growth of ICD width with the number of equivalent water neighbors; the non-linearity is stronger for Mg(2+). The ab initio ICD widths in microsolvated Na(+)-(H2O)4 and Mg(2+)-(H2O)6 clusters are found to be in good agreement with the experimental values.

Computational investigations of the thermodynamic properties of size-selected water and Ar-water clusters: high-pressure transitions

Classical parallel-tempering Monte Carlo simulations in the isothermal-isobaric ensemble were carried out for the (H2O)20 and Ar(H2O)20 clusters, over a wide range of temperatures (30-1000 K) and pressures (3 kPa-10 GPa) in order to study their thermodynamic properties and structural changes. The TIP4P/ice water model is employed for the water-water interactions, while both semiempirical and ab initio-based potentials are used to model the interaction between the rare-gas atoms and the water molecules. Temperature-pressure phase diagrams for these cluster systems were constructed by employing a two-dimensional multiple-histogram method. Structural changes were detected by analyzing the heat capacity landscape and the Pearson correlation coefficient profile for the interaction energy and volume. Those at high pressure correspond to solid-to-solid transitions and are found to be related to clathrate-like cages around the Ar atom. It is also shown that the formation and thermodynamic stability of such structures are determined by the intermolecular interaction between the rare-gas atoms and the host water molecules.

Water clusters in an argon matrix: infrared spectra from molecular dynamics simulations with a self-consistent charge density functional-based tight binding/force-field potential

The present theoretical study aims at investigating the effects of an argon matrix on the structures, energetics, dynamics, and infrared (IR) spectra of small water clusters (H2O)n (n = 1-6). The potential energy surface is obtained from a hybrid self-consistent charge density functional-based tight binding/force-field approach (SCC-DFTB/FF) in which the water clusters are treated at the SCC-DFTB level and the matrix is modeled at the FF level by a cluster consisting of ∼340 Ar atoms with a face centered cubic (fcc) structure, namely (H2O)n/Ar. With respect to a pure FF scheme, this allows a quantum description of the molecular system embedded in the matrix, along with all-atom geometry optimization and molecular dynamics (MD) simulations of the (H2O)n/Ar system. Finite-temperature IR spectra are derived from the MD simulations. The SCC-DFTB/FF scheme is first benchmarked on (H2O)Arn clusters against correlated wave function results and DFT calculations performed in the present work, and against FF data available in the literature. Regarding (H2O)n/Ar systems, the geometries of the water clusters are found to adapt to the fcc environment, possibly leading to intermolecular distortion and matrix perturbation. Several energetical quantities are estimated to characterize the water clusters in the matrix. In the particular case of the water hexamer, substitution and insertion energies for the prism, bag, and cage are found to be lower than that for the 6-member ring isomer. Finite-temperature MD simulations show that the water monomer has a quasifree rotation motion at 13 K, in agreement with experimental data. In the case of the water dimer, the only large-amplitude motion is a distortion-rotation intermolecular motion, whereas only vibration motions around the nuclei equilibrium positions are observed for clusters with larger sizes. Regarding the IR spectra, we find that the matrix environment leads to redshifts of the stretching modes and almost no shift of the bending modes. This is in agreement with experimental data. Furthermore, in the case of the water monomer and dimer, the magnitudes of the computed shifts are in fair agreement with the experimental values. The complex case of the water hexamer, which presents several low-energy isomers, is discussed.

Rovibrational energies and spectroscopic constants for H2O-Ng complexes

In this work, rovibrational energies and spectroscopic constants for the water -Ng complexes (Ng = He, Ne, Ar, Kr and Xe) were calculated through two different approaches (by solving the Nuclear Schrödinger equation and by applying the Dunham's method) and using two different potential energy curves (PEC). These PEC were determined using potential parameters obtained through molecular beam scattering experiments and accurate theoretical calculation, respectively. It was found that the theoretical rovibrational energies are in a good agreement (only for the lowest numbers of vibrational states) with those obtained through experimental PEC. Another important conclusions was regarding the calculated first two rovibrational energies for the H 2 O-Ar system, that are in a good agreement with the experimental data.

A new ab initio intermolecular potential energy surface and predicted rotational spectra of the Ne-H2S complex

We report a new three-dimensional ab initio intermolecular potential energy surface for the Ne-H(2)S complex with H(2)S monomer fixed at its experimental average structure. Using the supermolecular approach, the intermolecular potential energies were evaluated at CCSD(T) (coupled cluster with single and double and perturbative triple excitations) level with large basis sets including bond functions. The full counterpoise procedure was employed to correct the basis set superposition error. The planar T-shaped global minimum is located at the intermolecular distance of 3.51 Å with a well depth of 71.57 cm(-1). An additional planar local minimum was found to be separated from the global minimum with an energy barrier of 23.11 cm(-1). In addition, two first-order and one second-order saddle points were also located. The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to evaluate the rovibrational energy levels for eight isotopic species of the Ne-H(2)S complexes. The rotational transition frequencies for the eight isotopomers were also determined for the ground and first vibrational excited states, which are all in very good agreement with the available experimental values.

Water as a solute: competitive shell formation in (He,Ne) mixed microdroplets

Quantum Monte Carlo (QMC) stochastic calculations are carried out for a series of mixed rare gas clusters containing He and Ne which further include one H(2)O molecule as a single dopant. The ab initio potentials employed in the calculations, and the structural details provided by the QMC results, clearly reveal that the differences in the interaction forces which exist between the two solvent adatoms and the molecular solute are causing strongly competing environments that generate preferential shell structuring when surrounding the water molecule. The different behaviour of the two solvents, and the energetics of mixing, are analyzed in detail for small aggregates and for larger mixtures, revealing structural effects which originate from the different networking between solvent adatoms.

Ab initio characterization of the Ca-HCl van der Waals complex

The equilibrium structure and three-dimensional potential energy surface of the Ca-HCl van der Waals complex in its ground electronic state have been determined from accurate ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with basis sets of quadruple- and quintuple-zeta quality. The core-electron correlation, high-order valence-electron correlation, and scalar relativistic effects were investigated. The Ca-HCl complex was confirmed to be linear at equilibrium, with the vibrationless dissociation energy (into Ca and HCl) D(e) of 287 cm(-1). The vibration-rotation energy levels of various Ca-HCl isotopomers were predicted using the variational method. The predicted spectroscopic constants can be useful in a further analysis of high-resolution vibration-rotation spectra of the Ca-HCl complex.