The water dimer I: Experimental characterization

A theoretical study of the photochemistry of 1,3-cyclopentadiene and its cyano derivatives bound to a water dimer: Assessing reactivity of ionized clusters and possible photoproducts

We theoretically investigate the photoionization scenarios of molecular complexes involving cyclopentadiene and cyanocyclopentadiene bound to water dimers. Using electronic structure calculations within density-functional theory (DFT) and time dependent DFT (TD-DFT), we explore the potential photochemical pathways following ionization, and determine the charge transfer excitations related to the possible subsequent reactions. Our findings suggest that the investigated photochemical pathways of the hydrated complexes take place in two well-defined ultraviolet regions: (i) 8.2-9.5 eV for the cyclic compounds and (ii) 11.2-11.4 eV for the bound water dimer. We quantify how H-bonding effects can influence the photoionization channels. Before forming possible photoproducts, we also examine the regiospecificity of OH addition to 1,3-cyclopentadiene and its cyano derivatives We analyze our results in light of photoionization studies of jet-cooled molecular complexes and possible implications in astrochemical environments.

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.

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.

Effect of an Electric Field on the Structure and Stability of Atmospheric Clusters

We study the influence of an applied electric field on the structure and stability of some common bimolecular clusters that are found in the atmosphere. These clusters play an important role in new particle formation (NPF). For low values of the electric field (i.e., |E| ≤ 0.01 V Å-1), we demonstrate that the field response of the clusters can be predicted from simply calculating the dipole moment of the cluster and the dipole moments of the constituent molecules and that the influence on the association energy of the cluster is minimal (i.e., <0.5 kcal mol-1). For higher field strengths |E| > 0.2 V Å-1, there can be more dramatic effects on both structure and energetics, as the induced dipole, charge transfer, and geometric distortion play a larger role. Although such large fields are not very relevant in the atmosphere, they do exist in some situations of experimental interest, such as near interfaces and in intense laser fields.

The stability of oxygen-centered radicals and its response to hydrogen bonding interactions

The stability of various alkoxy/aryloxy/peroxy radicals, as well as TEMPO and triplet dioxygen (3 O2 ) has been explored at a variety of theoretical levels. Good correlations between RSEtheor and RSEexp are found for hybrid DFT methods, for compound schemes such as G3B3-D3, and also for DLPNO-CCSD(T) calculations. The effects of hydrogen bonding interactions on the stability of oxygen-centered radicals have been probed by addition of a single solvating water molecule. While this water molecule always acts as a H-bond donor to the oxygen-centered radical itself, it can act as a H-bond donor or acceptor to the respective closed-shell parent.

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.

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.

Ortho-para interconversion of nuclear states of H2O through replica transition state: prospect of quantum entanglement at homodromic Bjerrum defect site

CONTEXT

From a nuclear spin prospective, water exists as para and ortho nuclear spin isomers (isotopomers). Spin interconversions in isolated molecules of water are forbidden, but many recent reports have shown them to happen in bulk, through dynamic proton exchanges happening between interconnected networks of a large array of water molecules. In this contribution, a possible explanation for an unexpected slow or delayed interconversion of ortho-para water in ice observed in an earlier reported experiment is provided. Using the results of quantum mechanical investigations, we have discussed the roles played by Bjerrum defects in the dynamic proton exchanges and ortho-para spin state interconversions. We guess that at the sites of the Bjerrum defects, there are possibilities of quantum entanglements of states, through pairwise interactions. Based on the perfectly correlated exchange happening via a replica transition state, we speculate that it can have significant influences on ortho-para interconversions of water. We also conjecture that the overall ortho-para interconversion is not a continuous process, rather can be imagined to be happening serendipitously, but within the boundary of the rules of quantum mechanics.

METHODS

All computations were performed with Gaussian 09 program. B3LYP/6-31++G(d,p) methodology was used to compute all the stationary points. Further energy corrections were computed using CCSD(T)/aug-cc-pVTZ methodology. Intrinsic reaction coordinate (IRC) path computations were carried out for the transition states.

Analysis of the Electron Density of a Water Molecule Encapsulated by Two Cholic Acid Residues

Cholic acid is a trihydroxy bile acid with a nice peculiarity: the average distance between the oxygen atoms (O7 and O12) of the hydroxy groups located at C7 and C12 carbon atoms is 4.5 Å, a value which perfectly matches with the O/O tetrahedral edge distance in Ih ice. In the solid phase, they are involved in the formation of hydrogen bonds with other cholic acid units and solvents. This fact was satisfactorily used for designing a cholic dimer which encapsulates one single water molecule between two cholic residues, its oxygen atom (Ow) being exactly located at the centroid of a distorted tetrahedron formed by the four steroid hydroxy groups. The water molecule participates in four hydrogen bonds, with the water simultaneously being an acceptor from the 2 O12 (hydrogen lengths are 2.177 Å and 2.114 Å) and a donor towards the 2 O7 (hydrogen bond lengths are 1.866 Å and 1.920 Å). These facts suggest that this system can be a nice model for the theoretical study of the formation of ice-like structures. These are frequently proposed to describe the water structure found in a plethora of systems (water interfaces, metal complexes, solubilized hydrophobic species, proteins, and confined carbon nanotubes). The above tetrahedral structure is proposed as a reference model for those systems, and the results obtained from the application of the atoms in molecules theory are presented here. Furthermore, the structure of the whole system allows a division into two interesting subsystems in which water is the acceptor of one hydrogen bond and the donor of another. The analysis of the calculated electron density is performed through its gradient vector and the Laplacian. The calculation of the complexation energy used correction of the basis set superposition error (BSSE) with the counterpoise method. As expected, four critical points located in the H…O bond paths were identified. All calculated parameters obey the proposed criteria for hydrogen bonds. The total energy for the interaction in the tetrahedral structure is 54.29 kJ/mol, while the summation obtained of the two independent subsystems and the one between the alkyl rings without water is only 2.5 kJ/mol higher. This concordance, together with the calculated values for the electron density, the Laplacian of the electron density, and the lengths of the oxygen atom and the hydrogen atom (involved in the formation of each hydrogen bond) to the hydrogen bond critical point, suggests that each pair of hydrogen bonds can be considered independent of each other.

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.

Formation of heterogeneous clusters in superfluid helium nanodroplets: phthalocyanine and water

Clusters consisting of a single phthalocyanine molecule and a single water molecule are investigated by means of electronic spectroscopy in superfluid helium droplets. A recent spectroscopic study of those clusters [J. Fischer, F. Schlaghaufer, E.-M. Lottner, A. Slenczka, L. Christiansen, H. Stapelfeldt, M. Karra, B. Friedrich, T. Mullan, M. Schütz and D. Usvyat, J. Phys. Chem. A, 2019, 123, 10057-10064] which all exhibit a water induced electronic shift to the red is now complemented by the corresponding clusters exhibiting a water induced shift to the blue. These clusters will be analyzed by means of fluorescence excitation spectra, dispersed emission spectra, and additional experimental observations made feasible by helium droplets as cryogenic reactor. Together with the blue shifted clusters a total number of at least 6 isomeric variants could be identified in helium droplets. This contrasts to a number of only three isomeric variants obtained from quantum chemical calculations [J. Fischer, F. Schlaghaufer, E.-M. Lottner, A. Slenczka, L. Christiansen, H. Stapelfeldt, M. Karra, B. Friedrich, T. Mullan, M. Schütz and D. Usvyat, J. Phys. Chem. A, 2019, 123, 10057-10064] disregarding the helium environment and to a single isomer identified in a molecular beam experiment [J. Menapace and E. Bernstein, J. Chem. Phys., 1987, 87, 6877-6889]. The discrepancy in the number of isomers provides evidence of a profound involvement of helium in clustering. Moreover, the discrepancies between the gas phase experiment and quantum chemical calculations similarly reveal the influence of the dynamics of cluster formation on the population of global and local minima that are accessible as isomeric variants.

Water: new aspect of hydrogen bonding in the solid state

All water-water contacts in the crystal structures from the Cambridge Structural Database with d OO ≤ 4.0 Å have been found. These contacts were analysed on the basis of their geometries and interaction energies from CCSD(T)/CBS calculations. The results show 6729 attractive water-water contacts, of which 4717 are classical hydrogen bonds (d OH ≤ 3.0 Å and α ≥ 120°) with most being stronger than -3.3 kcal mol-1. Beyond the region of these hydrogen bonds, there is a large number of attractive interactions (2062). The majority are antiparallel dipolar interactions, where the O-H bonds of two water molecules lying in parallel planes are oriented antiparallel to each other. Developing geometric criteria for these antiparallel dipoles (β1, β2 ≥ 160°, 80 ≤ α ≤ 140° and T HOHO > 40°) yielded 1282 attractive contacts. The interaction energies of these antiparallel oriented water molecules are up to -4.7 kcal mol-1, while most of the contacts have interaction energies in the range -0.9 to -2.1 kcal mol-1. This study suggests that the geometric criteria for defining attractive water-water interactions should be broader than the classical hydrogen-bonding criteria, a change that may reveal undiscovered and unappreciated interactions controlling molecular structure and chemistry.

Assessing the Experimental Hydrogen Bonding Energy of the Cyclic Water Dimer Transition State with a Cyclooctatetraene-Based Molecular Balance

We have conducted an experimental and computational study of cyclooctatetraene-1,4/1,6-dimethanol (1,4 and 1,6) as a molecular balance with the goal in mind to determine the otherwise inaccessible hydrogen bonding energy (HBE) of the cyclic water dimer, which constitutes a transition state. The 1,4/1,6 folding equilibrium is governed by an intramolecular hydrogen bond in the folded 1,6-isomer, in which the OH groups adopt a cyclic planar geometry, akin to the structure of the cyclic water dimer transition state. We characterized hydrogen bonding in 1,6 and reference complexes utilizing SAPT2 + (3)δMP2/aug-cc-pVTZ and selected quantum theory of atoms in molecule descriptors at M06-2XD3(0)/ma-def2-TZVPP. Additionally, we computed HBEs at the DLPNO-CCSD(T)/aug-cc-pVQZ level of theory. We find that hydrogen bonding in 1,6 is very similar to the interaction in the C_i symmetric cyclic water dimer TS, both in magnitude and character. We experimentally determined the Gibbs free energy of the folding process (ΔG_eq) in a variety of organic solvents via nuclear magnetic resonance spectroscopy measurements at room temperature. By combining experimentally obtained ΔG_eq values with corrections derived from accurate computational methods, we provide estimates for the HBE of cyclic water dimers and the cyclic water dimer TS, as the most stable cyclic water dimer.

Equilibrium and Dynamical Characteristics of Hydrogen Bond Bifurcations in Water-Water and Water-Ammonia Dimers: A Path Integral Molecular Dynamics Study

We present path integral molecular dynamics results that describe the effects of nuclear quantum fluctuations on equilibrium and dynamical characteristics pertaining to bifurcation pathways in hydrogen bonded dimers combining water and ammonia, at cryogenic temperatures of the order of 20 K. Along these isomerizations, the hydrogen atoms in the molecules acting as hydrogen-bond donors interchange their original dangling/connective characters. Our results reveal that the resulting quantum transition paths comprise three stages: the initial and final ones involve overall rotations during which the two protons retain their classical-like characteristics. Effects from quantum fluctuation are clearly manifested in the changes operated at the intermediate passages over transition states, as the spatial extents of the protons stretch over typical lengths comparable to the distances between connective and dangling basins of attractions. Consequently, the classical over-the-hill path is replaced by a tunneling controlled mechanism which, within the path integral perspective, can be cast in terms of concerted inter-basin migrations of polymer beads from dangling-to-connective and from connective-to-dangling, at practically no energy costs. We also estimated the characteristic timescales describing such interconversions within the approximate ring polymer rate theory. Effects derived from full and partial deuteration are also discussed.

A Hessian-Free Method to Prevent Zero-Point Energy Leakage in Classical Trajectories

The problem associated with the zero-point energy (ZPE) leak in classical trajectory calculations is well known. Since ZPE is a manifestation of the quantum uncertainty principle, there are no restrictions on energy during the classical propagation of nuclei. This phenomenon can lead to unphysical results, such as forming products without the ZPE in the internal vibrational degrees of freedom (DOFs). The ZPE leakage also permits reactions below the quantum threshold for the reaction. We have developed a new Hessian-free method, inspired by the Lowe-Andersen thermostat model, to prevent energy dipping below a threshold in the local-pair (LP) vibrational DOFs. The idea is to pump the leaked energy to the corresponding local vibrational mode taken from the other vibrational DOFs. We have applied the new correction protocol on the ab-initio ground-state molecular dynamics simulation of the water dimer (H₂O)₂, which dissociates due to unphysical ZPE spilling from high-frequency OH modes. The LP-ZPE method has been able to prevent the ZPE spilling of the OH stretching modes by pumping back the leaked energy into the corresponding modes, while this energy is taken from the other modes of the dimer itself, keeping the system as a microcanonical ensemble.

Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses results obtained with the real-time equation-of-motion CC cumulant approach (RT-EOM-CC), and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the CCSD(T) formulation as a difference of energies for N and N-1 electron systems. All methods show good agreement with each other. They also agree well with experiment, with errors usually below 0.1 eV for the ionization potentials.We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

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.

Hydrogen Bonding in the Dimer and Monohydrate of 2-Adamantanol: A Test Case for Dispersion-Corrected Density Functional Methods

Weakly-bound intermolecular clusters constitute reductionist physical models for non-covalent interactions. Here we report the observation of the monomer, the dimer and the monohydrate of 2-adamantanol, a secondary alcohol with a bulky ten-carbon aliphatic skeleton. The molecular species were generated in a supersonic jet expansion and characterized using broadband chirped-pulse microwave spectroscopy in the 2–8 GHz frequency region. Two different gauche-gauche O-H···O hydrogen-bonded isomers were observed for the dimer of 2-adamantanol, while a single isomer was observed for the monomer and the monohydrate. The experimental rotational parameters were compared with molecular orbital calculations using density functional theory (B3LYP-D3(BJ), B2PLYP-D3(BJ), CAM-B3LYP-D3(BJ), ωB97XD), additionally providing energetic and electron density characterization. The shallow potential energy surface makes the dimer an interesting case study to benchmark dispersion-corrected computational methods and conformational search procedures.

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.

Confined Water Vapor in ZIF-8 Nanopores

Metal-organic frameworks (MOFs) possess an ordered and size-controllable porous structure, making them an interesting heterogeneous confining environment for water. Herein, molecular dynamics simulations are employed to investigate the structure of confined water vapor in zeolitic imidazolate framework-8 (ZIF-8) nanopores. Water dimers, which are rarely observed in liquid or water vapor, can form in ZIF-8 at room temperature. The six-ring-member gate is the main location of a water dimer in ZIF-8. The terminal methyl and CH groups of the imidazole linker interact with the water dimer by relatively weak hydrogen bonding. The above-presented findings provide a foundation for the elucidation of water confined in ZIF-8 and demonstrate the potential of obtaining low-order clusters of water by MOFs.

Nonlinear rotational spectroscopy reveals many-body interactions in water molecules

Because of their central importance in chemistry and biology, water molecules have been the subject of decades of intense spectroscopic investigations. Rotational spectroscopy of water vapor has yielded detailed information about the structure and dynamics of isolated water molecules, as well as water dimers and clusters. Nonlinear rotational spectroscopy in the terahertz regime has been developed recently to investigate the rotational dynamics of linear and symmetric-top molecules whose rotational energy levels are regularly spaced. However, it has not been applied to water or other lower-symmetry molecules with irregularly spaced levels. We report the use of recently developed two-dimensional (2D) terahertz rotational spectroscopy to observe high-order rotational coherences and correlations between rotational transitions that were previously unobservable. The results include two-quantum (2Q) peaks at frequencies that are shifted slightly from the sums of distinct rotational transitions on two different molecules. These results directly reveal the presence of previously unseen metastable water complexes with lifetimes of 100 ps or longer. Several such peaks observed at distinct 2Q frequencies indicate that the complexes have multiple preferred bimolecular geometries. Our results demonstrate the sensitivity of rotational correlations measured in 2D terahertz spectroscopy to molecular interactions and complexation in the gas phase.

Superphane: a new lantern-like receptor for encapsulation of a water dimer

A new superphane, featuring an aesthetically pleasing structure, was successfully obtained via one-pot synthesis of a hexakis-amine and m-phthalaldehyde in a [2+6] manner. It proved capable of entrapping a water dimer within its cavity as inferred from the mass spectroscopy, crystallographical analysis, NMR spectroscopy, and theoretical calculations.

Photoelectron Spectroscopy of the Water Dimer Reveals Unpredicted Vibrational Structure

Hydrogen bonds and proton transfer reactions can be considered as being at the very heart of aqueous chemistry and of utmost importance for many processes of biological relevance. Nevertheless, these processes are not yet well understood, even in seemingly simple model systems like small water clusters. We present a study of the photoelectron spectrum of the water dimer, revealing previously unresolved vibrational structure with 10-30 meV (80-242 cm^-1) typical splitting, in disagreement with a previous theoretical photoionization study predicting an apparent main vibrational progression with an ∼130 meV spacing [Kamarchik et al.; J. Chem. Phys. 2010, 132, 194311]. The observed vibrational structure and its deviation from the theoretical prediction is discussed in terms of known difficulties with calculations of strongly coupled anharmonic systems involving large amplitude motions. Potential contributions of the nonzero vibrational energy of the neutral water dimer at a finite experimental internal temperature are addressed. The internal temperature is estimated from the breakdown diagram associated with the dissociative ionization of the water dimer to be around to 130 K. This analysis also provides two additional, independently measured values for the 0 K appearance energy of the hydronium ion (H₃O⁺) from dissociative ionization of the water dimer.

New insights in chemical reactivity from quantum chemical topology

Based on the quantum chemical topology of the modified electron localization function ELF_x , an efficient and robust mechanistic methodology designed to identify the favorable reaction pathway between two reactants is proposed. We first recall and reshape how the supermolecular interaction energy can be evaluated from only three distinct terms, namely the intermolecular coulomb energy, the intermolecular exchange-correlation energy and the intramolecular energies of reactants. Thereafter, we show that the reactivity between the reactants is driven by the first-order variation in the coulomb intermolecular energy defined in terms of the response to changes in the number of electrons. Illustrative examples with the formation of the dative bond B-N involved in the BH₃ NH₃ molecule and the typical formation of the hydrogen bond in the canonical water dimer are presented. For these selected systems, our approach unveils a noticeable mimicking of E_dual onto the DFT intermolecular interaction energy surface calculated between the both reactants. An automated reaction-path algorithm aimed to determine the most favorable relative orientations when the two molecules approach each other is also outlined.

Pauli Repulsion Enhances Mobility of Ultraconfined Water

Water is ubiquitous on Earth and dominates chemical and biological processes in daily life. However, how water behaves under some critical conditions is not fully understood. In this paper, we employed quantum first-principles calculations and dynamics simulations to reveal the unexpectedly high mobility of water molecules in ultraconfined spaces. The water molecules rotated more freely in the (4, 4) carbon nanotube than in the (5, 5) carbon nanotube, which is induced by the Pauli repulsion from the wall of the narrower channel when reducing the size of the channel from general confinement to ultraconfinement. Moreover, this quantum effect facilitates the transport of water molecules into the space within their van der Waals diameter easily, which is in contrast to the general understanding. Thus, the conventional concept that the tighter the confined space, the more difficult the motion of the confined object is not always correct. This quantum-induced enhancement of water mobility by Pauli repulsion calls us to pay more attention to the existence and the function of water in neglected ultraconfined spaces (e.g., cells and the Earth's crust) in the future.

Switching Hydrogen Bonding to π-Stacking: The Thiophenol Dimer and Trimer

We used jet-cooled broadband rotational spectroscopy to explore the balance between π-stacking and hydrogen-bonding interactions in the self-aggregation of thiophenol. Two different isomers were detected for the thiophenol dimer, revealing dispersion-controlled π-stacked structures anchored by a long S-H···S sulfur hydrogen bond. The weak intermolecular forces allow for noticeable internal dynamics in the dimers, as tunneling splittings are observed for the global minimum. The large-amplitude motion is ascribed to a concerted inversion motion between the two rings, exchanging the roles of the proton donor and acceptor in the thiol groups. The determined torsional barrier of B₂ = 250.3 cm^-1 is consistent with theoretical predictions (290-502 cm^-1) and the monomer barrier of 277.1(3) cm^-1. For the thiophenol trimer, a symmetric top structure was assigned in the spectrum. The results highlight the relevance of substituent effects to modulate π-stacking geometries and the role of the sulfur-centered hydrogen bonds.

Revisiting the Spectroscopy of Water Dimer in Jets

Laser spectroscopy in jets is one of the main sources of structural data from molecular aggregates. Consequently, numerous and sophisticated experimental systems have been developed to extract precise information, which is usually interpreted in the light of quantum mechanical calculations. However, even with the most sophisticated experiments, it is sometimes difficult to interpret the experimental results. We present here the example of water dimer and how after almost 70 years, the assignment of its mass-resolved IR spectrum still generates controversy that extends toward the mechanism of ionization of water aggregates.

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.

A Step-by-Step Process-Induced Unidirectional Oriented Water Wire in the Nanotube

The orientation of water molecules is a key requirement for the fast transport of water in nanotubes. It has been accepted that the flip of the water chain follows a concerted mechanism, which has led to the view that bidirectional water flux in nanotubes can be transformed into unidirectional transport when the orientation of water molecules is maintained in long nanotubes under the external field. In this Letter, on the basis of molecular dynamics simulations and first-principles calculations, we confirmed that the flip of the water chain is a step-by-step process, which is different from the perceived concerted mechanism. Further analysis indicated that without an external field, it needed more than 20 water molecules to maintain the unidirectional single-file water flow in a carbon nanotube at a duration time of seconds. Considering that the thickness of the cell membrane (normally 5-10 nm) is larger than the length threshold of the unidirectional water wire, this study suggested that it may not require the external field to maintain the unidirectional flow in the water channel at the macroscopic time scale.

Evidence of a new phase in gypsum-anhydrite transformations under microwave heating by in situ dielectric analysis and Raman spectroscopy

Mineral transformations of the gypsum-anhydrite system under microwave heating have been studied using in situ dielectric thermal analysis (MW-DETA) and Raman spectroscopy simultaneously. The dielectric properties of samples that were measured under microwave heating provided thorough information about the dynamics of the gypsum-anhydrite system transformations and its significance from the mineralogical point of view. In particular, the MW-DETA technique revealed a new intermediate phase with a γ-anhydrite structure. This phase corresponds to the soluble stage of γ-anhydrite, and it is characterized by a high ionic charge inside the crystal channels. The complete sequence is gypsum → 0.625-subhydrate → bassanite → hydro γ-anhydrite → anhydrous γ-anhydrite → β-anhydrite. The transformations were also assessed using DSC, TG, DTA and dielectric measurements at room temperature, as well as other techniques including X-ray powder diffraction (XRPD) and high-temperature XRD (HT-XRD). Correlations between the dielectric properties with temperature and the rest of the techniques elucidated the heating mechanisms of this material under microwave energy during the different stages. The in situ combination of the MW-DETA and the Raman analysis appears to be a powerful technique, providing new insights about the mechanisms which govern the volumetric heating of this and other materials.

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.

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.

Exploring optimization strategies for improving explicit water models: Rigid n-point model and polarizable model based on Drude oscillator

Rigid n-point water models are widely used in atomistic simulations, but have known accuracy drawbacks. Increasing the number of point charges, as well as adding electronic polarizability, are two common strategies for accuracy improvements. Both strategies come at considerable computational cost, which weighs heavily against modest possible accuracy improvements in practical simulations. In an effort to provide guidance for model development, here we have explored the limiting accuracy of "electrostatically globally optimal" n-point water models in terms of their ability to reproduce properties of water dimer-a mimic of the condensed state of water. For a given n, each model is built upon a set of reference multipole moments (e.g. ab initio) and then optimized to reproduce water dimer total dipole moment. The models are then evaluated with respect to the accuracy of reproducing the geometry of the water dimer. We find that global optimization of the charge distribution alone can deliver high accuracy of the water model: for n = 4 or n = 5, the geometry of the resulting water dimer can be almost within 50 of the ab initio reference, which is half that of the experimental error margin. Thus, global optimization of the charge distribution of classical n-point water models can lead to high accuracy models. We also find that while the accuracy improvement in going from n = 3 to n = 4 is substantial, the additional accuracy increase in going from n = 4 to n = 5 is marginal. Next, we have explored accuracy limitations of the standard practice of adding electronic polarizability (via a Drude particle) to a "rigid base"-pre-optimization rigid n-point water model. The resulting model (n = 3) shows a relatively small improvement in accuracy, suggesting that the strategy of merely adding the polarizability to an inferior accuracy water model used as the base cannot fix the defects of the latter. An alternative strategy in which the parameters of the rigid base model are globally optimized along with the polarizability parameter is much more promising: the resulting 3-point polarizable model out-performs even the 5-point optimal rigid model by a large margin. We suggest that future development efforts consider 3- and 4-point polarizable models where global optimization of the "rigid base" is coupled to optimization of the polarizability to deliver globally optimal solutions.

Temperature Dependence of Intramolecular Vibrational Bands in Small Water Clusters

Cyclic water clusters are pivotal for understanding atmospheric reactions as well as liquid water, yet the temperature (T) dependence of their dynamics and spectroscopy is poorly studied. The development of highly accurate water potentials, such as MB-pol, partly rectifies this. It remains to account for the quantum nuclear effects (NQE), because quantum nuclear dynamics become increasingly inaccurate at low temperatures. From a practical point of view, we find that NQE can be accounted for simply by subtracting a constant from the frequencies obtained from the velocity autocorrelation functions (VACF) of classical nuclear dynamics, resulting in unprecedented agreement with experiment, mostly within 5 cm^-1. We have performed classical simulations of (H₂O)_n clusters (n = 2-5) from 20 K and up to their melting temperature, calculating both all-atom and partial VACF, thus generating the temperature dependence of the vibrational frequencies (IR and Raman bands). Focusing on the hydrogen-bonded (HBed) OH stretch and HOH bend, we find opposing T dependencies. The HBed OH modes blue shift linearly with T, attributed to ring expansion rather than any specific conformational change. The lowest-frequency Raman concerted mode is predicted to show the largest such shift. In contrast, the HOH bend undergoes a red-shift, with the highest frequency concerted band undergoing the largest red-shift. These results can be explained by a coupled-oscillator model for n hydrogen atoms on a ring, constrained to move either tangentially (stretch) or perpendicularly (bend) to the ring. With increasing temperature and weakening of HBs, the intrinsic force constant increases (stretch) or remains constant (bend), while the nearest-neighbor coupling constant decreases, and this results in the interesting behavior revealed herein. T-dependent Raman studies are required for testing some of these predictions.

Toward Elimination of Discrepancies between Theory and Experiment: Anharmonic Rotational-Vibrational Spectrum of Water in Solid Noble Gas Matrices

Rotational–vibrational spectroscopy of water in solid noble gas matrices has been studied for many decades. Despite that, discrepancies persist in the literature about the assignment of specific bands. We tackle the involved rotational–vibrational spectrum of the water isotopologues H₂¹⁶O, HD¹⁶O, and D₂¹⁶O with an unprecedented combination of experimental high-resolution matrix isolation infrared (MI-IR) spectroscopy and computational anharmonic vibrational spectroscopy by vibrational configuration interaction (VCI) on high-level ab initio potential energy surfaces. With VCI, the average deviation to gas-phase experiments is reduced from >100 to ≈1 cm^–1 when compared to harmonic vibrational spectra. Discrepancies between MI-IR and VCI spectra are identified as matrix effects rather than missing anharmonicity in the theoretical approach. Matrix effects are small in Ne (≈1.5 cm^–1) and a bit larger in Ar (≈10 cm^–1). Controversial assignments in Ne MI-IR spectra are resolved, for example, concerning the ν₃ triad in HDO. We identify new transitions, for example, the ν₂ 1₀₁ ← 1₁₀ transition in D₂O and H₂O or the ν₃ 0₀₀ ← 1₀₁ transition in D₂O, and reassign bands, for example, the band at 3718.9 cm^–1 that is newly assigned as the 1₁₀ ← 1₁₁ transition. The identification and solution of discrepancies for a well-studied benchmark system such as water prove the importance of an iterative and one-hand combination of theory and experiment in the field of high-resolution infrared spectroscopy of single molecules. As the computational costs involved in the VCI approach are reasonably low, such combined experimental/theoretical studies can be extended to molecules larger than triatomics.

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.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

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.