The complex of N2with H2O, D2O, and HDO: A combinedabinitioand diffusion Monte Carlo study

Spectroscopic study of the tunneling dynamics in N2-water observed in the O-D stretch region

The O-D stretch rovibrational spectra of N₂-D₂O and N₂-DOH were measured and analyzed. A combination band involving the in-plane N₂ bending vibration was also observed. These bands were recorded using a pulsed-slit supersonic jet expansion and a mid-infrared tunable optical parametric oscillator. The spectra were analyzed by considering the feasible tunneling motions, and transitions were fitted to independent asymmetric rotors for each tunneling state. The rotational constants of the four tunneling components of N₂-D₂O were retrieved for the excited vibrational states. A two order of magnitude increase in the tunneling splittings is observed for the asymmetric O-D stretch (ν₃ in D₂O) excitation compared to the symmetric stretch (ν₁ in D₂O) and to the ground vibrational state. This last finding indicates that the ν₃ vibrational state is likely perturbed by a combination state that includes ν₁. Finally, the observation of a local perturbation in the ν₃ vibrational band, affecting the positions of few rovibrational levels, provides an experimental lower limit of the dissociation energy of the complex, D₀ > 120 cm^-1.

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.

Structure and water attachment rates of ice in the atmosphere: role of nitrogen

In this work we perform computer simulations of the ice surface in order to elucidate the role of nitrogen in the crystal growth rates and crystal habits of snow in the atmosphere. In pure water vapor at temperatures typical of ice crystal formation in cirrus clouds, we find that basal and primary prismatic facets exhibit a layer of premelted ice, with thickness in the subnanometer range. For partial pressures of 1 bar, well above the expected values in the troposphere, we find that only small amounts of nitrogen are adsorbed. The adsorption takes place onto the premelted surface, and hardly any nitrogen dissolves within the premelting film. The premelting film thickness does not change either. We quantify the resulting change of the ice/vapor surface tension to be in the hundredth of mN m^-1 and find that the structure of the pristine ice surface is not changed in a significant manner. We perform a trajectory analysis of colliding water molecules, and find that the attachment rates from direct ballistic collision are very close to unity irrespective of the nitrogen pressure. Nitrogen is however at sufficient density to deflect a fraction of trajectories with smaller distance than the mean free path. Our results show explicitly that the reported differences in growth rates measured in pure water vapor and a controlled nitrogen atmosphere are not related to a significant disruption of the ice surface due to nitrogen adsorption. On the contrary, we show clearly from our trajectory analysis that nitrogen slows down the crystal growth rates due to collisions between water molecules with bulk nitrogen gas. This clarifies the long standing controversy of the role of inert gases on crystal growth rates and demonstrates their influence is solely related to the diffusion limited flow of water vapor across the gas phase.

The vibration-rotation-tunneling levels of N2-H2O and N2-D2O

In this paper, we report vibration-rotation-tunneling levels of the van der Waals clusters N2-H2O and N2-D2O computed from an ab initio potential energy surface. The only dynamical approximation is that the monomers are rigid. We use a symmetry adapted Lanczos algorithm and an uncoupled product basis set. The pattern of the cluster's levels is complicated by splittings caused by H-H exchange tunneling (larger splitting) and N-N exchange tunneling (smaller splitting). An interesting result that emerges from our calculation is that whereas in N2-H2O, the symmetric H-H tunnelling state is below the anti-symmetric H-H tunnelling state for both K = 0 and K = 1, the order is reversed in N2-D2O for K = 1. The only experimental splitting measurements are the D-D exchange tunneling splittings reported by Zhu et al. [J. Chem. Phys. 139, 214309 (2013)] for N2-D2O in the v2 = 1 region of D2O. Due to the inverted order of the split levels, they measure the sum of the K = 0 and K = 1 tunneling splittings, which is in excellent agreement with our calculated result. Other splittings we predict, in particular those of N2-H2O, may guide future experiments.

H(2)O--N(2) collision-induced absorption band intensity in the region of the N(2) fundamental: ab initio investigation of its temperature dependence and comparison with laboratory data

The present paper aims at ab initio and laboratory evaluation of the N(2) collision-induced absorption band intensity arising from interactions between N(2) and H(2)O molecules at wavelengths of around 4 μm. Quantum chemical calculations were performed in the space of five intermolecular coordinates and varying N--N bond length using Møller-Plesset perturbation and CCSD(T) methods with extrapolation of the electronic energy to the complete basis set. This made it possible to construct the intermolecular potential energy surface and to define the surface of the N--N dipole derivative with respect to internal coordinate. The intensity of the nitrogen fundamental was then calculated as a function of temperature using classical integration. Experimental spectra were recorded with a BOMEM DA3-002 FTIR spectrometer and 2 m base-length multipass White cell. Measurements were conducted at temperatures of 326, 339, 352 and 363 K. The retrieved water-nitrogen continuum significantly deviates from the MT_CKD model because the relatively strong nitrogen absorption induced by H(2)O was not included in this model. Substantial uncertainties in the measurements of the H(2)O-N(2) continuum meant that quantification of any temperature dependence was not possible. The comparison of the integrated N(2) fundamental band intensity with our theoretical estimates shows reasonably good agreement. Theory indicates that the intensity as a function of temperature has a minimum at approximately 500 K.

Infrared spectra and intensities of the H2O and N2 complexes in the range of the nu1- and nu3-bands of water

The IR spectra of complexes of water with nitrogen molecules in the range of the symmetric (nu(1)) and antisymmetric (nu(3)) bands of H(2)O have been studied in helium droplets. The infrared intensities of the nu(3) and nu(1) modes of H(2)O were found to be larger by factors of 1.3 and 2, respectively, in the N(2)-H(2)O complexes. These factors are smaller than those obtained in recent theoretical calculations. The conformation of the N(2)-H(2)O complex was estimated. Spectra and IR intensities of the (N(2))(2)-H(2)O and N(2)-(H(2)O)(2) complexes were also obtained and their structures are discussed.

Vibrational Spectra and Structure of CH3Cl:(H2O)2and CH3Cl:(D2O)2Complexes. IR Matrix Isolation and ab Initio Calculations

The infrared spectra of CH3Cl + H2O isolated in solid neon at low temperature have been investigated. High concentration studies of water (0.01%-4%) and subsequent annealing lead to the formation of the ternary CH3Cl:(H2O)2 complex. Detailed vibrational assignments were made on the observed spectra of water and deuterated water engaged in the complex. In parallel, structural, energetic, and vibrational properties of the complex have been studied at the second-order Møller-Plesset perturbation theory using several basis sets. Anaharmonic correction to the vibrational frequencies has been done with the standard second-order perturbation approach. It was shown that the ground state of the complex has a cyclic form for which the nonadditive three-body contribution was found to be around 10% of the interaction energy.

Stereographic projection diffusion monte carlo (SPDMC) algorithms for molecular condensed matter

We develop and test three algorithms for diffusion Monte Carlo simulations in non-Euclidean manifolds. The methods are based on the construction of the "velocity" distribution by rejection techniques and are capable of functioning in a broad class of non-Euclidean spaces generated by holonomic constraints. The formulation of the propagator for non-Euclidean manifolds avoids the use of Lagrange multipliers; it is derived instead from the Feynman quantization in manifolds proposed by DeWitt. The manifolds are mapped onto Rd by using stereographic projection coordinates. Numerical tests are conducted for the particle in a ring of unit radius subjected to a sinusoidal potential, for the electron in the field of an infinitely massive proton, and for a water molecule modeled as an asymmetric top subjected to an external field.

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.

Vibrational Spectra and Structures of H2O−NO, HDO−NO, and D2O−NO Complexes. An IR Matrix Isolation and DFT Study

The IR spectra of H2O+NO, HDO+NO, and D2O+NO, isolated in solid neon at low temperature have been investigated. Concentration effects and detailed vibrational analysis of deuterated and partially deuterated species allowed identification of three 1:1 HDO-NO species, two 1:1 D2O-NO species, and only one 1:1 H2O-NO complex. From comparison between the experimental spectra and the results of DFT calculations, it appeared that two different types of weakly bound complexes between water and nitric oxide can be formed in a neon matrix. The first species is a 1:1 complex where bonding occurs between water hydrogen and nitric oxide nitrogen, in which OH-N and OD-N intermolecular bonds are engaged. For this complex only DOD-NO, HOD-NO, and DOH-NO isotopic species have been experimentally detected and no IR bands of HOH-NO were observed. This result could be explained by the fact that the dissociation energy of HOH-NO is lower than those of DOD-NO, HOD-NO and DOH-NO. For the second detected 1:1 H2O-NO complex and its isotopic variants, the H2O-NO potential surface was explored systematically at the B3LYP level, but no stable species corresponding to the complex could be calculated. The structure of the second observed 1:1 H2O-NO complex results from columbic attractions between water and nitric oxide and could be stabilized only in matrix, probably by interaction between NO, water and (Ne)n.

Infrared water vapor continuum absorption at atmospheric temperatures

We have used a continuous-wave carbon dioxide laser in a single-mode realization of cavity ring-down spectroscopy to measure absorption coefficients of water vapor at 944 cm(-1) for several temperatures in the range 270-315 K. The conventional description of water vapor infrared absorption is applied, in which the absorption is modeled in two parts consisting of local line absorption and the remaining residual absorption, which has become known as the water vapor continuum. This water vapor continuum consists of distinct water-water, water-nitrogen, and water-oxygen continua. The water-water continuum absorption coefficient is found to have a magnitude of C(s)(296 K) = (1.82+/-0.02) x 10(-22) cm(2) molecule(-1) atm(-1), and the water-nitrogen coefficient has a magnitude of C(n)(296 K) = (7.3 +/- 0.4) x 10(-25) cm(2) molecule(-1) atm(-1). The temperature dependences of both the water-water and the water-nitrogen continua are shown to be well represented by a model describing the expected behavior of weakly bound binary complexes. Using this model, our data yield dissociation energies of D(e) = (-15.9 +/- 0.3) kJ/mole for the water dimer and D(e) = (-3.2 +/- 1.7) kJ/mole for the water-nitrogen complex. These values are in excellent agreement with recent theoretical predictions of D(e) = -15.7 kJ/mole (water dimer) and D(e) = -2.9 kJ/mole (water-nitrogen complex), as well as the experimentally determined value of D(e) = (-15.3 +/- 2.1) kJ/mole for the water dimer obtained by investigators employing a thermal conductivity technique. Although there is reasonably good agreement with the magnitude of the continuum absorption coefficients, the agreement on temperature dependence is less satisfactory. While our results are suggestive of the role played by water dimers and water complexes in producing the infrared continuum, the uncertain spectroscopy of the water dimer in this spectral region prevents us from making a firm conclusion. In the meantime, empirical models of water vapor continuum absorption, essential for atmospheric radiative transfer calculations, should be refined to give better agreement with our low-uncertainty continuum absorption data.

A 'first principles' potential energy surface for liquid water from VRT spectroscopy of water clusters

We present results of gas phase cluster and liquid water simulations from the recently determined VRT(ASP-W)III water dimer potential energy surface (the third fitting of the Anisotropic Site Potential with Woermer dispersion to vibration-rotation-tunnelling data). VRT(ASP-W)III is shown to not only be a model of high 'spectroscopic' accuracy for the water dimer, but also makes accurate predictions of vibrational ground-state properties for clusters up through the hexamer. Results of ambient liquid water simulations from VRT(ASP-W)III are compared with those from ab initio molecular dynamics, other potentials of 'spectroscopic' accuracy and with experiment. The results herein represent the first time to the authors' knowledge that a 'spectroscopic' potential surface is able to correctly model condensed phase properties of water.

Elucidating the role of many-body forces in liquid water. I. Simulations of water clusters on the VRT(ASP-W) potential surfaces

We test two new potentials for water, fit to vibration-rotation tunneling (VRT) data by employing diffusion quantum Monte Carlo simulations to calculate the vibrational ground-state properties of water clusters. These potentials, VRT(ASP-W)II and VRT(ASP-W)III, are fits of the highly detailed ASP-W (anisotropic site potential with Woermer dispersion) ab initio potential to (D(2)O)(2) microwave and far-infrared data, and along with the SAPT5s (five-site symmetry adapted perturbation theory) potentials, are the most accurate water dimer potential surfaces in the literature. The results from VRT(ASP-W)II and III are compared to those from the original ASP-W potential, the SAPT5s family of potentials, and several bulk water potentials. Only VRT(ASP-W)III and the spectroscopically "tuned" SAPT5st (with N-body induction included) accurately reproduce the vibrational ground-state structures of water clusters up to the hexamer. Finally, the importance of many-body induction and three-body dispersion are examined, and it is shown that the latter can have significant effects on water cluster properties despite its small magnitude.