Intermolecular potential and second virial coefficient of the water–helium complex

A close coupling study of the bending relaxation of H2O by collision with He

We present a close coupling study of the bending relaxation of H₂O by collision with He, taking explicitly into account the bending-rotation coupling within the rigid-bender close-coupling method. A 4D potential energy surface is developed based on a large grid of ab initio points calculated at the coupled-cluster single double triple level of theory. The bound states energies of the He-H₂O complex are computed and found to be in excellent agreement with previous theoretical calculations. The dynamics results also compare very well with the rigid-rotor results available in the Basecol database and with experimental data for both rotational transitions and bending relaxation. The bending-rotation coupling is also demonstrated to be very efficient in increasing bending relaxation when the rotational excitation of H₂O increases.

Homogeneous water nucleation: Experimental study on pressure and carrier gas effects

Homogeneous nucleation of water is investigated in argon and in nitrogen at about 240 K and 0.1 MPa, 1 MPa, and 2 MPa by means of a pulse expansion wave tube. The surface tension reduction at high pressure qualitatively explains the observed enhancement of the nucleation rate of water in argon as well as in nitrogen. The differences in nucleation rates for the two mixtures at high pressure are consistent with the differences in adsorption behavior of the different carrier gas molecules. At low pressure, there is not enough carrier gas available to ensure the growing clusters are adequately thermalized by collisions with carrier gas molecules so that the nucleation rate is lower than under isothermal conditions. This reduction depends on the carrier gas, pressure, and temperature. A qualitative agreement between experiments and theory is found for argon and nitrogen as carrier gases. As expected, the reduction in the nucleation rates is more pronounced at higher temperatures. For helium as the carrier gas, non-isothermal effects appear to be substantially stronger than predicted by theory. The critical cluster sizes are determined experimentally and theoretically according to the Gibbs-Thomson equation, showing a reasonable agreement as documented in the literature. Finally, we propose an empirical correction of the classical nucleation theory for the nucleation rate calculation. The empirical expression is in agreement with the experimental data for the analyzed mixtures (water-helium, water-argon, and water-nitrogen) and thermodynamic conditions (0.06 MPa-2 MPa and 220 K-260 K).

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.

Zero- and high-pressure mechanisms in the complex forming reactions of OH with methanol and formaldehyde at low temperatures

A recent Ring Polymer Molecular Dynamics study of the reactions of OH with methanol and formaldehyde, at zero pressure and below 100 K, has shown the formation of long lived complexes, with long lifetimes, longer than 100 ns for the lower temperatures studied, 20-100 K (del Mazo-Sevillano et al., 2019). These long lifetimes support the existence of multi collision events with the He buffer-gas atoms under experimental conditions, as suggested by several transition state theory studies of these reactions. In this work we study these secondary collisions, as a dynamical approach to study pressure effects on these reactions. For this purpose, the potential energy surfaces of He with H2CO, OH, H2O and HCO are calculated at highly accurate ab initio level. The stability of some of the complexes is studied using Path Integral Molecular dynamics techniques, determining that OH-H2CO complexes can be formed up to 100 K or higher temperatures, while the weaker He-H2CO complexes dissociate at approximately 50 K. The predicted IR intensity spectra shows new features which could help the identification of the OH-H2CO complex. Finally, the He-H2CO + OH and OH-H2CO + He collisions are studied using quassi-classical trajectories, finding that the cross section to produce HCO + H2O products increases with decreasing collision energy, and that it is ten times higher in the He-H2CO + OH case.

Collision energy dependence of state-to-state differential cross sections for rotationally inelastic scattering of H2O by He

The inelastic scattering of H₂O by He as a function of collision energy in the range 381 cm^-1 to 763 cm^-1 at an energy interval of approximately 100 cm^-1 has been investigated in a crossed beam experiment using velocity map imaging. Change in collision energy was achieved by varying the collision angle between the H₂O and He beam. We measured the state-to-state differential cross section (DCS) of scattered H₂O products for the final rotational states J_KaKc = 1₁₀, 1₁₁, 2₂₁ and 4₁₄. Rotational excitation of H₂O is probed by (2 + 1) resonance enhanced multiphoton ionization (REMPI) spectroscopy. DCS measurements over a wide range of collision energies allowed us to probe the H₂O-He potential energy surface (PES) with greater detail than in previous work. We found that a classical approximation of rotational rainbows can predict the collision energy dependence of the DCS. Close-coupling quantum mechanical calculations were used to produce DCS and partial cross sections. The forward-backward ratio (FBR), is introduced here to compare the experimental and theoretical DCS. Both theory and experiments suggest that an increase in the collision energy is accompanied with more forward scattering.

Helium Diffraction as a Probe of Structure and Proton Order on Model Ice Surfaces

Helium diffraction has the potential to reveal the degree of proton order at an ice surface, and has been used in the past to benchmark theoretical work. We demonstrate that previous calculations do not represent the diffraction experiment to a sufficient degree of accuracy. By combining a realistic helium-water potential with quantum calculations using exact close-coupling methods we demonstrate that the scattering is strongly energy dependent. Proton order may be inferred best from selective adsorption resonances of the helium atom, which involve multiple scattering. We use the results to discuss the validity of the latest assumptions for the ice Ih surface with respect to proton ordering.

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.

Dipole oscillator strength distributions, properties, and dispersion energies for ethylene, propene, and 1-butene

A recommended isotropic dipole oscillator strength distribution (DOSD) has been constructed for the ethylene molecule through the use of quantum mechanical constraint techniques and experimental dipole oscillator strength (DOS) data; the DOS data employed are recent experimental results not available at the time of the original constrained DOSD analysis of this molecule. The constraints are furnished by molar refractivity data and the Thomas–Reiche–Kuhn sum rule. The DOSD is used to evaluate a variety of isotropic dipole oscillator strength sums, logarithmic dipole oscillator strength sums, and mean excitation energies for ethylene. Pseudo-DOSDs for this molecule, and for propene and 1–butene, which are based on an earlier constrained DOSD analysis for these molecules, are developed. They are used to obtain reliable results for the isotropic dipole–dipole dispersion-energy coefficients C6, for the interactions of the alkenes with each other and with 47 other species, and the triple-dipole dispersion-energy coefficients C9 for interactions involving any triple of molecules taken from ethylene, propene, and 1–butene.Key words: alkenes, dipole properties, pseudo-states, dipole–dipole and triple-dipole dispersion energies, long-range additive, non-additive interaction energies.

Close-coupling study of rotational energy transfer and differential scattering in H2O collisions with He atoms

Quantum close-coupling scattering calculations of rotational energy transfer (RET) of rotationally excited H(2)O due to collisions with He are presented for collision energies between 10(-6) and 1000 cm(-1) with para-H(2)O initially in levels 1(1,1), 2(0,2), 2(1,1), and 2(2,0) and ortho-H(2)O in levels 1(1,0), 2(1,2), and 2(2,1). Quenching cross sections and rate coefficients for state-to-state RET were computed. Both elastic and inelastic differential cross sections are also calculated and compared with relative experimental results giving generally good agreement in all cases, but less so for inelastic results. Significant differences in the computed collisional parameters, obtained on three different potential energy surfaces (PESs), were found particularly in the ultracold regime. In the thermal regime, the rate coefficients calculated on each of the surfaces are generally in better agreement and comparable, but typically larger, than those obtained in a previous calculation. Unfortunately, a lack of absolute differential or integral inelastic experimental data prevents firm determination of a preferred PES.

Intermolecular potential and second virial coefficient of the water-nitrogen complex

The authors construct a rigid-body (five-dimensional) potential energy surface for the water-nitrogen complex using the systematic intermolecular potential extrapolation routine. The intermolecular potential is then extrapolated to the limit of a complete basis set. An analytic fit of this surface is obtained, and, using this, the global minimum energy is found. The minimum is located in an arrangement in which N2 is near the H atom of H2O, almost collinear with the OH bond. The best estimate of the binding energy is 441 cm-1 (1 cm-1 approximately 1.986 43x10(-23) J). The extrapolated potential is then used to calculate the second cross virial coefficient over a wide temperature range (100-3000 K). These calculated second virial coefficients are generally consistent with experimental data, but for the most part the former have smaller uncertainties.

Methane-water cross second virial coefficient with quantum corrections from an ab initio potential

We present our calculations of the cross second virial coefficient (B12) and of a related quantity, phi 12 = B12-TdB12/dT, for the methane-water system in the temperature range T = 200-1000 K. These calculations were performed using one of the ab initio potentials developed in previous work [Akin-Ojo and Szalewicz, J. Chem. Phys. 123, 134311 (2005)]. Quantum corrections of order variant Planck's over 2pi(2) were added to the computed classical values. We have estimated the uncertainties in our computed B12 and phi 12(T). This allowed evaluation of the quality of the experimental data to which we compare our results. We also provide an analytical expression for B12(T) as a function of the temperature T obtained by fitting the calculated values. This formula also predicts values of phi12(T) consistent with the directly calculated values.

Dipole Oscillator Strength Distributions and Properties for Methanol, Ethanol and Propan-1-ol and Related Dispersion Energies

Recommended isotropic dipole oscillator strength distributions (DOSDs) have been constructed for the methanol and ethanol molecules through the use of quantum mechanical constraint techniques and experimental dipole oscillator strength (DOS) data; the DOS data employed are recent experimental results not available at the time of the original constrained DOSD analysis of these molecules. The constraints are furnished by molar refractivity data and the Thomas-Reiche-Kuhn sum rule. The DOSDs are used to evaluate a variety of isotropic dipole oscillator strength sums, logarithmic dipole oscillator strength sums, and mean excitation energies for the molecules. Pseudo-DOSDs for these molecules, and for propan-1-ol based on an earlier constrained DOSD analysis for this molecule, are also presented. They are used to obtain reliable results for the isotropic dipole-dipole dispersion energy coefficients C6, for the interactions of the alcohols with each other and with 36 other species, and the triple-dipole dispersion energy coefficients C9for interactions involving any triple of molecules involving methanol, ethanol and propan-1-ol.

Intermolecular potential and second virial coefficient of the water–hydrogen complex

We construct a rigid-body (five-dimensional) potential-energy surface for the water-hydrogen complex using scaled perturbation theory (SPT). An analytic fit of this surface is obtained, and, using this, two minima are found. The global minimum has C2v symmetry, with the hydrogen molecule acting as a proton donor to the oxygen atom on water. A local minimum with Cs symmetry has the hydrogen molecule acting as a proton acceptor to one of the hydrogen atoms on water, where the OH bond and H2 are in a T-shaped configuration. The SPT global minimum is bound by 1097 microEh (Eh approximately 4.359744 x 10(-18) J). Our best estimate of the binding energy, from a complete basis set extrapolation of coupled-cluster calculations, is 1076.1 microEh. The fitted surface is used to calculate the second cross virial coefficient over a wide temperature range (100-3000 K). Three complementary methods are used to quantify quantum statistical mechanical effects that become significant at low temperatures. We compare our results with experimental data, which are available over a smaller temperature range (230-700 K). Generally good agreement is found, but the experimental data are subject to larger uncertainties.

An environmental pseudopotential approach to molecular interactions: Implementation in MOLPRO

We present the implementation into the MOLPRO package of a model for the interaction of a central system with its surrounding environment. The properties of a target system enclosed by a noncovalently bound environment or solvent are modeled as those of a system embedded into the effective pseudopotential arising from the exact electrostatic Coulomb potential and the approximated exchange-repulsion potential. For the latter we use the charge-density overlap model, which relates the exchange-repulsion interaction energy between two species with the overlap of their ground-state electron charge densities. The solutions of the modified Hartree-Fock equations for the target system are obtained self-consistently. This way the exchange-induction effects arising from the converged electron-charge density of the embedded system are implicitly included. Inclusion of the correlation effects is provided by the use of post-Hartree-Fock and density-functional techniques available in the MOLPRO package. The computational and conceptual advantages provided by this approach are shown in the calculation of the dipole polarizabilities of halide and chalcogenide anions in different environments.