Intermolecular Vibrations of the Water Trimer, a Matrix Isolation Study

Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory

This work systematically examines the interactions between a single argon atom and the edges and faces of cyclic H2O clusters containing three-five water molecules (Ar(H2O)n=3-5). Full geometry optimizations and subsequent harmonic vibrational frequency computations were performed using MP2 with a triple-ζ correlation consistent basis set augmented with diffuse functions on the heavy atoms (cc-pVTZ for H and aug-cc-pVTZ for O and Ar; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the two-body-many-body (2b:Mb) and three-body-many-body (3b:Mb) techniques; here, high-level CCSD(T) computations capture up through the two-body or three-body contributions from the many-body expansion, respectively, while less demanding MP2 computations recover all higher-order contributions. Five unique stationary points have been identified in which Ar binds to the cyclic water trimer, along with four for (H2O)4 and three for (H2O)5. To the best of our knowledge, eleven of these twelve structures have been characterized here for the first time. Ar consistently binds more strongly to the faces than the edges of the cyclic (H2O)n clusters, by as much as a factor of two. The 3b:Mb electronic energies computed with the haTZ basis set indicate that Ar binds to the faces of the water clusters by at least 3 kJ mol-1 and by nearly 6 kJ mol-1 for one Ar(H2O)5 complex. An analysis of the interaction energies for the different binding motifs based on symmetry-adapted perturbation theory (SAPT) indicates that dispersion interactions are primarily responsible for the observed trends. The binding of a single Ar atom to a face of these cyclic water clusters can induce perturbations to the harmonic vibrational frequencies on the order of 5 cm-1 for some hydrogen-bonded OH stretching frequencies.

Wetting of a Hydrophobic Surface: Far-IR Action Spectroscopy and Dynamics of Microhydrated Naphthalene

The interaction of water and polycyclic aromatic hydrocarbons is of fundamental importance in areas as diverse as materials science and atmospheric and interstellar chemistry. The interplay between hydrogen bonding and dipole-π interactions results in subtle dynamics that are challenging to describe from first principles. Here, we employ far-IR action vibrational spectroscopy with the infrared free-electron laser FELIX to investigate naphthalene with one to three water molecules. We observe diffuse bands associated with intermolecular vibrational modes that serve as direct probes of the loose binding of water to the naphthalene surface. These signatures are poorly reproduced by static DFT or Møller-Plesset computations. Instead, a rationalization is achieved through Born-Oppenheimer Molecular Dynamics simulations, revealing the active mobility of water over the surface, even at low temperatures. Therefore, our work provides direct insights into the wetting interactions associated with shallow potential energy surfaces while simultaneously demonstrating a solid experimental-computational framework for their investigation.

The effect of alkylation on the micro-solvation of ethers revealed by highly localized water librational motion

The specific far-infrared spectral signatures associated with highly localized large-amplitude out-of-plane librational motion of water molecules have recently been demonstrated to provide sensitive spectroscopic probes for the micro-solvation of organic molecules [Mihrin et al., Phys. Chem. Chem. Phys. 21(4), 1717 (2019)]. The present work employs this direct far-infrared spectroscopic approach to investigate the non-covalent intermolecular forces involved in the micro-solvation of a selection of seven ether molecules with systematically varied alkyl substituents: dimethyl ether, diethyl ether, diisopropyl ether, ethyl methyl ether, t-butyl methyl ether, and t-butyl ethyl ether. The ranking of the observed out-of-plane water librational band signatures for this selected series of ether–water complexes embedded in inert neon matrices at 4 K reveals information about the interplay of directional intermolecular hydrogen bond motifs and non-directional and long-range dispersion interactions for the micro-solvated structures. These far-infrared observables differentiate minor subtle effects introduced by specific alkyl substituents and serve as rigorous experimental benchmarks for modern quantum chemical methodologies of various levels of scalability, which often fail to accurately predict the structural variations and corresponding vibrational signatures of the closely related systems. The accurate interaction energies of the series of ether–water complexes have been predicted by the domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations, followed by a local energy decomposition analysis of the energy components. In some cases, the secondary dispersion forces are in direct competition with the primary intermolecular hydrogen bonds as witnessed by the specific out-of-plane librational signatures.

A Rather Universal Vibrational Resonance in 1:1 Hydrates of Carbonyl Compounds

When the lower frequency OH stretching fundamental of a water molecule is shifted to the 3500 cm^-1 spectral range by the solvation of a carbonyl compound, in this case a ketone, its infrared intensity is shared with a dark state. It is shown by chemical and isotope substitution for more than a dozen systems that the location of this resonance is remarkably substitution-independent. Harmonic and anharmonic model calculations support its assignment to a combination of the water bending overtone and in-plane water libration. This previously unrecognized intramolecular-intermolecular coupling in single solvent water has a strength of 7-10 cm^-1. It may have been sporadically observed before in a few other carbonyl compounds such as amides, without any previous exploration of its potential universality. The resulting generic picosecond energy redistribution channel for aqueous solutions may represent a slow counterpart and doorway model of what happens on a subpicosecond time scale when the hydrogen bonds become stronger, such as in carboxylic acid dimers or protonated water clusters.

Lactic Acid Spectroscopy: Intra- and Intermolecular Interactions

Lactic acid, a relevant molecule in biology and the environment, is an α-hydroxy acid with a high propensity to form hydrogen bonds, both internally and to other hydrogen-bond-accepting molecules. This work includes the novel recording of infrared spectra of gas-phase lactic acid using Fourier transform infrared spectroscopy, and the vibrational absorption features of lactic acid are assigned with the aid of computationally simulated vibrational spectra with anharmonic corrections. Theoretical chemistry methods are used to relate intramolecular hydrogen-bond strengths to the relative stability of lactic acid conformers. The formation of hydrogen-bonded lactic acid dimers and 1:1 water complexes is investigated by simulated vibrational spectra and calculated thermodynamic parameters for the lactic acid monomer and dimer and its water complex in the gas phase. The results of this study are discussed in the context of environmental chemistry with an emphasis on indoor environments.

Terahertz and mid-infrared spectroscopy of matrix-isolated clusters and matrix-sublimation ice of D2O

We have established an apparatus for terahertz and mid-infrared spectroscopy in an ultrahigh vacuum and have measured absorption spectra of D₂O clusters trapped in solid Ar. To assign terahertz absorption peaks due to the D₂O dimer, trimer, and tetramer, the dependence of the spectrum on the annealing temperature and D₂O dilution was analyzed. The assignment was also examined by ab initio calculations with the use of the "our own N-layered integrated molecular orbital and molecular mechanics" method, where the flexibility of surrounding Ar atoms was systematically incorporated. We identified all the intermolecular fundamentals of the dimer and those with significant intensities of the trimer and tetramer, whose structural symmetries were revealed to be broken down. After isolating the D₂O clusters in solid Ar, we sublimated only Ar atoms to leave behind matrix-sublimation ice, which was found to be amorphous- or crystal-like depending on the formation conditions: the dilution and sublimation temperature. The crystallinity of matrix-sublimation ice was determined by decomposing its terahertz spectrum into the spectra of amorphous and crystalline ices. Since the crystallinity got higher by raising the dilution and sublimation temperature, the diffusion of the D₂O monomer on the surface of sublimating solid Ar was found to be crucial to the crystallization of the sublimation ice.

The water trimer reaction OH + (H2O)3→ (H2O)2OH + H2O

All important stationary points on the potential energy surface (PES) for the reaction OH + (H₂O)₃→ (H₂O)₂OH + H₂O have been fully optimized using the "gold standard" CCSD(T) method with the large Dunning correlation-consistent cc-pVQZ basis sets. Three types of pathways were found. For the pathway without hydrogen abstraction, the barrier height of the transition state (TS1) is predicted to lie 5.9 kcal mol^-1 below the reactants. The two major complexes (H₂O)₃OH (CP1 and CP2a) are found to lie 6.3 and 11.0 kcal mol^-1, respectively, below the reactants [OH + (H₂O)₃]. For one of the H-abstraction pathways the lowest classical barrier height is predicted to be much higher, 6.1 kcal mol^-1 (TS2a) above the reactants. For the other H-abstraction pathway the barrier height is even higher, 15.0 (TS3) kcal mol^-1. Vibrational frequencies and the zero-point vibrational energies connected to the PES are also reported. The energy barriers for the H-abstraction pathways are compared with those for the OH + (H₂O)₂ and OH + H₂O reactions, and the effects of the third water on the energetics are usually minor (0.2 kcal mol^-1).

Highly localized H2O librational motion as a far-infrared spectroscopic probe for microsolvation of organic molecules

The most prominent spectroscopic observable for the hydrogen bonding between individual molecules in liquid water is the broad absorption band detected in the spectral region between 300 and 900 cm-1. The present work demonstrates how the associated large-amplitude out-of-plane OH librational motion of H2O molecules also directly reflects the microsolvation of organic compounds. This highly localized OH librational motion of the first solvating H2O molecule causes a significant change of dipole moment and gives rise to a strong characteristic band in the far-infrared spectral region, which is correlated quantitatively with the complexation energy. The out-of-plane OH librational band origins ranging from 324.5 to 658.9 cm-1 have been assigned experimentally for a series of four binary hydrogen-bonded H2O complexes embedded in solid neon involving S-, O- and N-containing compounds with increasing hydrogen bond acceptor capability. The hydrogen bond energies for altogether eight binary H2O complexes relative to the experimental value of 13.2 ± 0.12 kJ mol-1 for the prototypical (H2O)2 system [Rocher-Casterline et al., J. Chem. Phys., 2011, 134, 211101] are revealed directly by these far-infrared spectroscopic observables. The far-infrared spectral signatures are able to capture even minor differences in the hydrogen bond acceptor capability of O atoms with slightly different alkyl substituents in the order H-O-C(CH3)3 > CH3-O-CH3 > H-O-CH(CH3)2 > H-O-CH2CH3.

Conformational Landscape of the 1/1 Diacetyl/Water Complex Investigated by Infrared Spectroscopy and ab Initio Calculations

The complexes of diacetyl with water have been studied experimentally by Fourier transform infrared (FTIR) spectroscopy coupled to solid neon matrix and supersonic jet, and anharmonic ab initio calculations. The vibrational analysis of neon matrix spectra over the 100-7500 cm^-1 infrared range confirms the existence of two nearly isoenergetic one-to-one (1/1) diacetyl-water S₁ and S₂ isomers already evidenced in a previous argon matrix study. A third form (S₃) predicted slightly less stable ( J. Mol. Mod. 2015 , 21 , 214 ) is not observed. The correct agreement obtained between neon matrix and anharmonic calculated vibrational frequencies is exploited in several cases to derive band assignments for the vibrational modes of a specific isomer. Thereafter, theoretical x_ij anharmonic coupling constants are used for the attribution of combination bands and overtones relative to the 1/1 dimer. Finally, the most stable isomer of the one-to-two (1/2) diacetyl-water complex is identified in the OH stretching region of water on the grounds of comparison of experimental and calculated vibrational shifts between water dimer and the three most stable 1/2 isomers.

Fast crystalline ice formation at extremely low temperature through water/neon matrix sublimation

Crystalline ice formation requires water molecules to be sufficiently mobile to find and settle on the thermodynamically most stable site.

Spectroscopic identification of ethanol-water conformers by large-amplitude hydrogen bond librational modes

The far-infrared absorption spectra have been recorded for hydrogen-bonded complexes of water with ethanol embedded in cryogenic neon matrices at 2.8 K. The partial isotopic H/D-substitution of the ethanol subunit enabled by a dual inlet deposition procedure enables the observation and unambiguous assignment of the intermolecular high-frequency out-of-plane and the low-frequency in-plane donor OH librational modes for two different conformations of the mixed binary ethanol/water complex. The resolved donor OH librational bands confirm directly previous experimental evidence that ethanol acts as the O⋯HO hydrogen bond acceptor in the two most stable conformations. In the most stable conformation, the water subunit forces the ethanol molecule into its less stable gauche configuration upon dimerization owing to a cooperative secondary weak O⋯HC hydrogen bond interaction evidenced by a significantly blue-shift of the low-frequency in-plane donor OH librational band origin. The strong correlation between the low-frequency in-plane donor OH librational motion and the secondary intermolecular O⋯HC hydrogen bond is demonstrated by electronic structure calculations. The experimental findings are further supported by CCSD(T)-F12/aug-cc-pVQZ calculations of the conformational energy differences together with second-order vibrational perturbation theory calculations of the large-amplitude donor OH librational band origins.

Calculations of the IR spectra of bend fundamentals of (H2O)n=3,4,5 using the WHBB_2 potential and dipole moment surfaces

Stimulated by new experiments from the Havenith group, we report IR spectra of the bend fundamentals of (H₂O)_n=3,4,5, using anharmonic, coupled-mode VSCF/VCI calculations, done in a subspace of modes consisting of all the monomer bends plus the hydrogen-bonded OH stretches.

Infrared spectroscopy of the ν2 band of the water monomer and small water clusters (H2O)n=2,3,4 in helium droplets

Ro-vibrational transitions of water and water clusters (H₂O)_n=1,2,3,4 for the ν₂ bending vibration were observed and assigned to distinct structures.

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.

The influence of large-amplitude librational motion on the hydrogen bond energy for alcohol–water complexes

The intermolecular large-amplitude OH librational modes for mixed hydrogen-bonded complexes of water with methanol and t-butanol are unambiguously assigned for the first time.

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.

Determination of Critical Indices by “Slow” Spectroscopy: NMR Shifts by Statistical Thermodynamics and Density Functional Theory Calculations

The temperature dependencies of NMR shifts in the critical region of two coexisting phases have been simulated using statistical thermodynamics and graph-theory consideration of equilibrium processes of molecular association. Microparameters of magnetic screening of various water and water/pyridine structures used in the statistical averaging have been evaluated by density functional theory calculations (PBE1PBE and B3PW91 functionals in the 6-311++G** basis set). The gauge-including atomic orbital (GIAO) approach has been applied to ensure gauge invariance of the results. Solvent effects were taken into account by a polarized continuum model (PCM). NMR shifts "order parameters" (Deltadelta = |delta+ - delta-|) and "diameters" (phidelta = |(delta+ + delta-)/2 - deltaC|, where delta+, delta-, and deltaC are the chemical shifts of coexisting phases and at the critical point respectively) have been calculated in each case close to the lower critical solution point (TL) and processed using linear regression analysis of Deltadelta approximately |T - TL| and phidelta approximately |T - TL| in the log-log plot. It has been shown that the critical index beta can be evaluated with high precision from the slope of Deltadelta = f(T - TL) at any realistic set of model input parameters. The slope of diameter has been found to depend on both input beta and alpha values. The obtained phidelta slopes (0.58-0.63) are very close to 2beta values. The results are discussed within the concept of complete scaling. Results of simulation are compared and supported by experimental NMR data for water/2,6-lutidine, acetic anhydride/n-heptane, and acetic anhydride/cyclohexane systems.

Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.