Infrared spectroscopy of small size‐selected water clusters

Spectroscopic and Theoretical Identifications of Two Structural Motifs of (H2O)10 Cluster

Precise characterization of archetypal systems of aqueous hydrogen-bonding networks is essential for developing accurate potential functions and universal models of water. The structures of water clusters (H2O)n (n = 2-9) have been verified recently through size-specific infrared spectroscopy with a vacuum ultraviolet free electron laser (VUV-FEL) and quantum chemical studies. For (H2O)10, the pentagonal prism and butterfly motifs were proposed to be important building blocks and were observed in previous experiments. Here we report the size-specific infrared spectra of (H2O)10 via a joint experimental and theoretical study. Well-resolved spectra provide a unique signature for the coexistence of pentagonal prism and butterfly motifs. These (H2O)10 motifs develop from the dominant structures of (H2O)n (n = 8, 9) clusters. This work provides an intriguing prelude to the diverse structure of liquid water and opens avenues for size-dependent measurement of larger systems to understand the stepwise formation mechanism of hydrogen-bonding networks.

Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods

Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological, and materials systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.

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.

Plausible photomolecular effect leading to water evaporation exceeding the thermal limit

We report in this work several unexpected experimental observations on evaporation from hydrogels under visible light illumination. 1) Partially wetted hydrogels become absorbing in the visible spectral range, where the absorption by both the water and the hydrogel materials is negligible. 2) Illumination of hydrogel under solar or visible-spectrum light-emitting diode leads to evaporation rates exceeding the thermal evaporation limit, even in hydrogels without additional absorbers. 3) The evaporation rates are wavelength dependent, peaking at 520 nm. 4) Temperature of the vapor phase becomes cooler under light illumination and shows a flat region due to breaking-up of the clusters that saturates air. And 5) vapor phase transmission spectra under light show new features and peak shifts. We interpret these observations by introducing the hypothesis that photons in the visible spectrum can cleave water clusters off surfaces due to large electrical field gradients and quadrupole force on molecular clusters. We call the light-induced evaporation process the photomolecular effect. The photomolecular evaporation might be happening widely in nature, potentially impacting climate and plants' growth, and can be exploited for clean water and energy technologies.

Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region

We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.

Water binding and hygroscopicity in π-conjugated polyelectrolytes

The presence of water strongly influences structure, dynamics and properties of ion-containing soft matter. Yet, the hydration of such matter is not well understood. Here, we show through a large study of monovalent π-conjugated polyelectrolytes that their reversible hydration, up to several water molecules per ion pair, occurs chiefly at the interface between the ion clusters and the hydrophobic matrix without disrupting ion packing. This establishes the appropriate model to be surface hydration, not the often-assumed internal hydration of the ion clusters. Through detailed analysis of desorption energies and O-H vibrational frequencies, together with OPLS4 and DFT calculations, we have elucidated key binding motifs of the sorbed water. Type-I water, which desorbs below 50 °C, corresponds to hydrogen-bonded water clusters constituting secondary hydration. Type-II water, which typically desorbs over 50-150 °C, corresponds to water bound to the anion under the influence of a proximal cation, or to a cation‒anion pair, at the cluster surface. This constitutes primary hydration. Type-III water, which irreversibly desorbs beyond 150 °C, corresponds to water kinetically trapped between ions. Its amount varies strongly with processing and heat treatment. As a consequence, hygroscopicity-which is the water sorption capacity per ion pair-depends not only on the ions, but also their cluster morphology.

Microhydration of the adamantane cation: intracluster proton transfer to solvent in [Ad(H2O)n=1-5]+ for n ≥ 3

Radical cations of diamondoids are important intermediates in their functionalization reactions in polar solvents. To explore the role of the solvent at the molecular level, we characterize herein microhydrated radical cation clusters of the parent molecule of the diamondoid family, adamantane (C₁₀H₁₆, Ad), by infrared photodissociation (IRPD) spectroscopy of mass-selected [Ad(H₂O)_n=1-5]⁺ clusters. IRPD spectra of the cation ground electronic state recorded in the CH/OH stretch and fingerprint ranges reveal the first steps of this fundamental H-substitution reaction at the molecular level. Analysis of size-dependent frequency shifts with dispersion-corrected density functional theory calculations (B3LYP-D3/cc-pVTZ) provides detailed information about the acidity of the proton of Ad⁺ as a function of the degree of hydration, the structure of the hydration shell, and the strengths of the CH⋯O and OH⋯O hydrogen bonds (H-bonds) of the hydration network. For n = 1, H₂O strongly activates the acidic C-H bond of Ad⁺ by acting as a proton acceptor in a strong CH⋯O ionic H-bond with cation-dipole configuration. For n = 2, the proton is almost equally shared between the adamantyl radical (C₁₀H₁₅, Ady) and the (H₂O)₂ dimer in a strong C⋯H⋯O ionic H-bond. For n ≥ 3, the proton is completely transferred to the H-bonded hydration network. The threshold for this size-dependent intracluster proton transfer to solvent is consistent with the proton affinities of Ady and (H₂O)_n and confirmed by collision-induced dissociation experiments. Comparison with other related microhydrated cations reveals that the acidity of the CH proton of Ad⁺ is in the range of strongly acidic phenol⁺ but lower than for cationic linear alkanes such as pentane⁺. Significantly, the presented IRPD spectra of microhydrated Ad⁺ provide the first spectroscopic molecular-level insight of the chemical reactivity and reaction mechanism of the important class of transient diamondoid radical cations in aqueous solution.

Microhydration of the Pyrrole Cation (Py+) Revealed by IR Spectroscopy: Ionization-Induced Rearrangement of the Hydrogen-Bonded Network of Py+(H2O)2

Microhydration of heterocyclic aromatic molecules can be an appropriate fundamental model to shed light on intermolecular interactions and functions of macromolecules and biomolecules. We characterize herein the microhydration process of the pyrrole cation (Py⁺) by infrared photodissociation (IRPD) spectroscopy and dispersion-corrected density functional theory calculations (B3LYP-D3/aug-cc-pVTZ). Analysis of IRPD spectra of mass-selected Py⁺(H₂O)₂ and its cold Ar-tagged cluster in the NH and OH stretch range combined with geometric parameters of intermolecular structures, binding energies, and natural atomic charge distribution provides a clear picture of the growth of the hydration shell and cooperativity effects. Py⁺(H₂O)₂ is formed by stepwise hydration of the acidic NH group of Py⁺ by a hydrogen-bonded (H₂O)₂ chain with NH···OH···OH configuration. In this linear H-bonded hydration chain, strong cooperativity, mainly arising from the positive charge, strengthens both the NH···O and OH···O H-bonds with respect to those of Py⁺H₂O and (H₂O)₂, respectively. The linear chain structure of the Py⁺(H₂O)₂ cation is discussed in terms of the ionization-induced rearrangement of the hydration shell of the neutral Py(H₂O)₂ global minimum characterized by the so-called "σ-π bridge structure" featuring a cyclic NH···OH···OH···π H-bonded network. Emission of the π electron from Py by ionization generates a repulsive interaction between the positive π site of Py⁺ and the π-bonded OH hydrogen of (H₂O)₂, thereby breaking this OH···π hydrogen bond and driving the hydration structure toward the linear chain motif of the global minimum on the cation potential.

Nature and Strength of Sulfur-Centered Hydrogen Bond in Methanethiol Aqueous Solutions

Methanethiol (M) and water (W) clusters like dimers (M₁W₁, M₂, and W₂), trimers (M₁W₂, M₂W₁, M₃, and W₃), and tetramers (M₁W₃, M₂W₂, M₃W₁, M₄, and W₄) were studied to assess the strength of sulfur-centered hydrogen bonding using different levels of theories, viz, HF, MP2, MP3, MP4, B3LYP, B3LYP-D3, CCSD, CCSD(T)-F12, and CCSD(T) along with aug-cc-pVNZ (where N = D, T, and Q) basis sets. Interaction energies were found to be in the range of -3.3 to -5.3 kcal/mol for the dimers, -8.0 to -16.7 kcal/mol for the trimers, and -13.5 to -29.5 kcal/mol for the tetramers at the B3LYP-D3/CBS limit level of theory. Normal modes of vibrations computed at the B3LYP/cc-pVDZ level of theory were seen to be in good agreement with the experimental values. Local energy decomposition calculations using the DLPNO-CCSD(T) level of theory indicated the domination of electrostatic interactions' contribution to the interaction energy in all cluster systems. Furthermore, atoms in molecules and natural bond orbital calculations both carried out at the B3LYP-D3/aug-cc-pVQZ level of theory aided in visualizing the hydrogen bonds besides proving a rationale for the strength of the hydrogen bonds and thereby the stability of these cluster systems.

An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2-25

We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2-11, 16-17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D_e's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2-25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83-119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D_e's) throughout the entire cluster range (n = 2-25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ∼4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6-3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.

Microhydrated clusters of a pharmaceutical drug: infrared spectra and structures of amantadineH+(H2O)n

Solvation of pharmaceutical drugs has an important effect on their structure and function. Analysis of infrared photodissociation spectra of amantadineH⁺(H₂O)_n=1-4 clusters in the sensitive OH, NH, and CH stretch range by quantum chemical calculations (B3LYP-D3/cc-pVTZ) provides a first impression of the interaction of this pharmaceutically active cation with water at the molecular level. The size-dependent frequency shifts reveal detailed information about the acidity of the protons of the NH₃⁺ group of N-protonated amantadineH⁺ (AmaH⁺) and the strength of the NH⋯O and OH⋯O hydrogen bonds (H-bonds) of the hydration network. The preferred cluster growth begins with sequential hydration of the NH₃⁺ group by NH⋯O ionic H-bonds (n = 1-3), followed by the extension of the solvent network through OH⋯O H-bonds. However, smaller populations of cluster isomers with an H-bonded solvent network and free N-H bonds are already observed for n ≥ 2, indicating the subtle competition between noncooperative ion hydration and cooperative H-bonding. Interestingly, cyclic water ring structures are identified for n ≥ 3, each with two NH⋯O and two OH⋯O H-bonds. Despite the increasing destabilization of the N-H proton donor bonds upon gradual hydration, no proton transfer to the (H₂O)_n solvent cluster is observed up to n = 4. In addition to ammonium cluster ions, a small population of microhydrated iminium isomers is also detected, which is substantially lower for the hydrophilic H₂O than for the hydrophobic Ar environment.

Spectroscopic snapshot for neutral water nonamer (H2O)9: Adding a H2O onto a hydrogen bond-unbroken edge of (H2O)8

Structural characterization of neutral water clusters is crucial to understanding the structures and properties of water, but it has been proven to be a challenging experimental target due to the difficulty in size selection. Here, we report the size-specific infrared spectra of confinement-free neutral water nonamer (H₂O)₉ based on threshold photoionization, using a tunable vacuum ultraviolet free-electron laser. Distinct OH stretch vibrational fundamentals in the 3200-3350 cm^-1 region are observed, providing unique spectral signatures for the formation of an unprecedented (H₂O)₉ structure evolved by adding a ninth water molecule onto a hydrogen bond-unbroken edge of the (H₂O)₈ octamer with D_2d symmetry. This nonamer structure coexists with the five previously identified structures that can be viewed as derived by inserting a ninth water molecule into a hydrogen bond-broken edge of the D_2d/S₄ octamer. These findings provide key microscopic information for systematic understanding of the formation and growth mechanism of dynamical hydrogen-bonding networks that are responsible for the structure and properties of condensed-phase water.

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.

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.

Raman study of water deposited in solid argon matrix

New Raman data are presented concerning H₂O and D₂O water aggregation in argon matrix having the ratio of number of argon atoms to water molecules close to 40:1. Experiments were conducted at temperatures from 8 K to 34 K allowing observation of OH and OD stretching vibrations of water monomers, dimers, trimers and higher multimers, as well as broad bands corresponding to solid amorphous water. Molecular dynamics simulations were performed for thirteen or sometimes fourteen water molecules dispersed among 500 argon atoms. Resulting final configurations included dimers, trimers, tetramers and pentamers, all in open chain configurations which upon optimization resulted in mostly cyclic conformations. Observed OH stretching vibrations were assigned by comparing calculated normal modes in harmonic approximation at the B3LYP/aug-cc-pVDZ and PBEPBE1/aug-cc-pVDZ level of theory with our data and previously observed bands from infrared matrix isolation studies and Raman jet cooled experiments. Raman bands assigned to water multimers in argon matrix are shifted 20 to 25 cm^-1 towards lower wavenumbers with respect to the positions of OH stretching vibrations of almost free water clusters.

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.

Bonding analysis of water clusters using quasi-atomic orbitals

The quasi-atomic orbital (QUAO) bonding analysis introduced by Ruedenberg and co-workers is used to develop an understanding of the hydrogen bonds in small water clusters, from the dimer through the hexamer (bag, boat, book, cyclic, prism and cage conformers). Using kinetic bond orders as a metric, it is demonstrated that as the number of waters in simple cyclic clusters increases, the hydrogen bonds strengthen, from the dimer through the cyclic hexamer. However, for the more complex hexamer isomers, the strength of the hydrogen bonds varies, depending on whether the cluster contains double acceptors and/or double donors. The QUAO analysis also reveals the three-center bonding nature of hydrogen bonds in water clusters.

In Situ Characterization of Dehydration during Ion Transport in Polymeric Nanochannels

The transport of hydrated ions across nanochannels is central to biological systems and membrane-based applications, yet little is known about their hydrated structure during transport due to the absence of in situ characterization techniques. Herein, we report experimentally resolved ion dehydration during transmembrane transport using modified in situ liquid ToF-SIMS in combination with MD simulations for a mechanistic reasoning. Notably, complete dehydration was not necessary for transport to occur across membranes with sub-nanometer pores. Partial shedding of water molecules from ion solvation shells, observed as a decrease in the average hydration number, allowed the alkali-metal ions studied here (lithium, sodium, and potassium) to permeate membranes with pores smaller than their solvated size. We find that ions generally cannot hold more than two water molecules during this sterically limited transport. In nanopores larger than the size of the solvation shell, we show that ionic mobility governs the ion hydration number distribution. Viscous effects, such as interactions with carboxyl groups inside the membrane, preferentially hinder the transport of the mono- and dihydrates. Our novel technique for studying ion solvation in situ represents a significant technological leap for the nanofluidics field and may enable important advances in ion separation, biosensing, and battery applications.

Ordered Water Layer on the Macroscopically Hydrophobic Fluorinated Polymer Surface and Its Ultrafast Vibrational Dynamics

Hydrophobic-like water monolayers have been predicted at the metal and some polar surfaces by theoretical simulations. However, direct experimental evidence for the presence of this water layer at surfaces, particularly at biomolecule and polymer surfaces, is yet to be validated at room temperature. Here we observe experimentally that an ordered molecular water layer is present at the hydrophobic fluorinated polymer such as polytetrafluoroethylene (PTFE) surface by using sum frequency generation vibrational spectroscopy. The macroscopic hydrophobicity of PTFE surface is actually hydrophilic at the molecular level. The macroscopically hydrophobic character of PTFE is indeed resulting from the hydrophobicity of the ordered two-dimension (2D) water layer, in which cyclic water tetramer structure is found. The water layer at humidity of ≤40% has a vibrational relaxation time of 550 ± 60 fs. The vibrational relaxation time in the frequency range of 3200-3400 cm^-1 shows remarkable difference from the interfacial water at the air/H₂O interface and the lipid/H₂O interface. No discernible frequency dependence of the vibrational relaxation time is observed, indicating the homogeneous dynamics of OH groups in the water layer. These insights into the water layer at the macroscopically hydrophobic surface may contribute to a better understanding of the hydrophobic interaction and interfacial water dynamics.

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.

Microhydration of substituted diamondoid radical cations of biological relevance: infrared spectra of amantadine+-(H2O)n = 1-3 clusters

Hydration of biomolecules and pharmaceutical compounds has a strong impact on their structure, reactivity, and function. Herein, we explore the microhydration structure around the radical cation of the widespread pharmaceutical drug amantadine (C16H15NH2, Ama) by infrared photodissociation (IRPD) spectroscopy of mass-selected Ama+Wn = 1-3 clusters (W = H2O) recorded in the NH, CH, and OH stretch range of the cation ground electronic state. Analysis of the size-dependent frequency shifts by dispersion-corrected density functional theory calculations (B3LYP-D3/cc-pVTZ) provides detailed information about the acidity of the protons of the NH2 group of Ama+ and the structure and strength of the NHO and OHO hydrogen bonds (H-bonds) of the hydration network. The preferred sequential cluster growth begins with hydration of the two acidic NH protons of the NH2 group (n = 1-2) and continues with an extension of the H-bonded hydration network by forming an OHO H-bond of the third W to one ligand in the first hydration subshell (n = 3), like in the W2 dimer. For n = 2, a minor population corresponds to Ama+W2 structures with a W2 unit attached to Ama+via a NHW2 H-bond. Although the N-H proton donor bonds are progressively destabilized by gradual microhydration, no proton transfer to the Wn solvent cluster is observed in the investigated size range (n ≤ 3). Besides the microhydration structure, we also obtain a first impression of the structure and IR spectrum of bare Ama+, as well as the effects of both ionization and hydration on the structure of the adamantyl cage. Comparison of Ama+ with aliphatic and aromatic primary amine radical cations reveals differences in the acidity of the NH2 group and the resulting interaction with W caused by substitution of the cycloalkyl cage.

Water accelerated self-healing of hydrophobic copolymers

Previous studies have shown that copolymer compositions can significantly impact self-healing properties. This was accomplished by enhancement of van der Waals (vdW) forces which facilitate self-healing in relatively narrow copolymer compositional range. In this work we report the acceleration of self-healing in alternating/random hydrophobic acrylic-based copolymers in the presence of confined water molecules. Under these conditions competing vdW interactions do not allow H₂O-diester H-bonding, thus forcing nBA side groups to adapt L-shape conformations, generating stronger dipole-dipole interactions resulting in shorter inter-chain distances compared to 'key-and-lock' associations without water. The perturbation of vdW forces upon mechanical damage in the presence of controllable amount of confined water is energetically unfavorable leading the enhancement of self-healing efficiency of hydrophobic copolymers by a factor of three. The concept may be applicable to other self-healing mechanisms involving reversible covalent bonding, supramolecular chemistry, or polymers with phase-separated morphologies.

Infrared spectroscopic study of hydrogen bonding topologies in the smallest ice cube

The water octamer with its cubic structure consisting of six four-membered rings presents an excellent cluster system for unraveling the cooperative interactions driven by subtle changes in the hydrogen-bonding topology. Despite prediction of many distinct structures, it has not been possible to extract the structural information encoded in their vibrational spectra because this requires size-selectivity of the neutral clusters with sufficient resolution to identify the contributions of the different isomeric forms. Here we report the size-specific infrared spectra of the isolated cold, neutral water octamer using a scheme based on threshold photoionization using a tunable vacuum ultraviolet free electron laser. A plethora of sharp vibrational bands features are observed. Theoretical analysis of these patterns reveals the coexistence of five cubic isomers, including two with chirality. The relative energies of these structures are found to reflect topology-dependent, delocalized multi-center hydrogen-bonding interactions. These results demonstrate that even with a common structural motif, the degree of cooperativity among the hydrogen-bonding network creates a hierarchy of distinct species. The implications of these results on possible metastable forms of ice are speculated.

Spectroscopic studies of water clusters provide insight into the hydrogen bond structure of water and ice. The authors measure infrared spectra of neutral water octamers using a threshold photoionization technique based on a tunable vacuum-UV free electron laser, identifying two cubic isomers in addition to those previously observed.

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.

Microhydration of protonated biomolecular building blocks: protonated pyrimidine

Protonation and hydration of biomolecules govern their structure, conformation, and function. Herein, we explore the microhydration structure in mass-selected protonated pyrimidine-water clusters (H⁺Pym-W_n, n = 1-4) by a combination of infrared photodissociation spectroscopy (IRPD) between 2450 and 3900 cm^-1 and density functional theory (DFT) calculations at the dispersion-corrected B3LYP-D3/aug-cc-pVTZ level. We further present the IR spectrum of H⁺Pym-N₂ to evaluate the effect of solvent polarity on the intrinsic molecular parameters of H⁺Pym. Our combined spectroscopic and computational approach unequivocally shows that protonation of Pym occurs at one of the two equivalent basic ring N atoms and that the ligands in H⁺Pym-L (L = N₂ or W) preferentially form linear H-bonds to the resulting acidic NH group. Successive addition of water ligands results in the formation of a H-bonded solvent network which increasingly weakens the NH group. Despite substantial activation of the N-H bond upon microhydration, no intracluster proton transfer occurs up to n = 4 because of the balance of relative proton affinities of Pym and W_n and the involved solvation energies. Comparison to neutral Pym-W_n clusters reveals the drastic effects of protonation on microhydration with respect to both structure and interaction strength.

Infrared spectroscopy of neutral water clusters at finite temperature: Evidence for a noncyclic pentamer

Infrared spectroscopic study of neutral water clusters is crucial to understanding of the hydrogen-bonding networks in liquid water and ice. Here we report infrared spectra of size-selected neutral water clusters, (H₂O) _n (n = 3-6), in the OH stretching vibration region, based on threshold photoionization using a tunable vacuum ultraviolet free-electron laser. Distinct OH stretch vibrational fundamentals observed in the 3,500-3,600-cm^-1 region of (H₂O)₅ provide unique spectral signatures for the formation of a noncyclic pentamer, which coexists with the global-minimum cyclic structure previously identified in the gas phase. The main features of infrared spectra of the pentamer and hexamer, (H₂O) _n (n = 5 and 6), span the entire OH stretching band of liquid water, suggesting that they start to exhibit the richness and diversity of hydrogen-bonding networks in bulk water.

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.

Facile assembly of 2D α-zirconium phosphate supported silver nanoparticles: superior and recyclable catalysis

A novel, efficient and durable two-dimensional ZrP@PDA/Ag nanocatalyst for the reduction of 4-nitrophenol.

Experimental Investigation of Hydrophobically Modified α-ZrP Nanosheets for Enhancing Oil Recovery in Low-Permeability Sandstone Cores

Highly crystalline α-zirconium phosphate (α-ZrP) nanoparticles were synthesized and exfoliated into nanosheets, and then the hydrophilic nanosheets were modified into hydrophobic nanosheets with octadecyltrichlorosilane (OTS). Scanning electron microscopy, Fourier transform infrared spectroscopy, X-ray diffraction, and thermogravimetric analysis were applied to confirm the morphology and chemical structure of the nanosheets. Contact angle measurement was conducted to explore the wettability alteration of the hydrophobically modified α-ZrP nanosheets, and the result showed that the wettability of the core was changed into hydrophobicity. When ZrP-OTS nanosheets were injected during water-flooding, a Pickering emulsion will be formed. The droplet diameters and viscosities of the Pickering emulsions were measured. The hydrophobically modified α-ZrP nanosheets were applied in low-permeability sandstone cores and various concentrations were tested. The injectivity of the hydrophobically modified α-ZrP nanosheets was also studied and the result indicated that the nanosheets exhibit good injectivity. The mechanisms for enhancing oil recovery by utilizing hydrophobic α-ZrP nanosheets were analyzed: forming Pickering emulsions and increasing the viscosity of the displacing phase. Forming emulsions and increasing the viscosity of the flooding phase can enhance the microdisplacement efficiency, while good injectivity can also enhance the macrodisplacement efficiency. The result indicated the possibility of using hydrophobically modified α-ZrP nanosheets for enhancing oil recovery in a low-permeability reservoir.

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.

Drastic Change in Electronic Transition upon Hydrogen Bond Network Switching in 3-Aminopyridine-(H2O) n Clusters

The hydration structures of 3-aminopyridine (3AP)-(H₂O) _n ( n = 2-4) in supersonic jets have been investigated by measuring the electronic and vibrational spectra with the aid of quantum chemical calculations. The S₁-S₀ electronic transition of 3AP-(H₂O)₂ is observed at a slightly red-shifted position from 3AP-(H₂O)₁, while further hydration induces drastic red shifts and complicated vibrational structures. We assign the cluster structures of 3AP-(H₂O)₂ as a cyclic structure composed of the homodromic hydrogen bond (H-bond) chain connecting the pyridyl CH bond acting as the proton donor toward a pyridyl nitrogen acceptor. For 3AP-(H₂O) _n ( n = 3, 4), on the other hand, the initial donor site in the H-bond network changes from a pyridyl CH group to an amino group. The observed red shift derived from H-bond network switching can be reproduced very well with the S₁-S₀ origin band estimation obtained by applying geometry optimization and subsequent harmonic vibrational analysis of (TD-)DFT calculations to each electronic state of the isomer structure. It is suggested that the drastic red shift of the electronic transition upon H-bond network switching is due to a much more stabilized "quinoid-like" structure in the ππ* state by the H-bond formation of an amino group.

Relaxation of the H2O Overtone Bending Vibration in the Water Dimer···Hydroxyl Radical Complex

The relaxation mechanism of the overtone bending vibration in the collision of the water dimer with the vibrationally excited hydroxyl radical is studied by use of trajectory procedures. The transfer of the OH(v = 1) energy to the dimer stretches is followed by a near-resonant first overtone transition to the donor monomer. Nearly a quarter of the trajectories undergo a complex-mode collision, forming the (H₂O)₂···OH complex bound by a hydrogen bond with the lifetime ranging from a subpicosecond scale to >100 ps. The overtone vibration relaxes to the ground state, transferring approximately half of its energy to the dimer hydrogen-bonding (H₂O···H₂O) and the remaining half to the complex hydrogen-bonding (H₂O)₂···OH, via near-resonant pathways, each consisting of a series of intermolecular low-frequency vibrations.

Energy Transfer to the Hydrogen Bond in the (H2O)2 + H2O Collision

Trajectory procedures are used to study the collision between the vibrationally excited H₂O and the ground-state (H₂O)₂ with particular reference to energy transfer to the hydrogen bond through the inter- and intramolecular pathways. In nearly 98% of the trajectories, energy transfer processes occur on a subpicosecond scale (≤0.7 ps). The H₂O transfers approximately three-quarters of its excitation energy to the OH stretches of the dimer. The first step of the intramolecular pathway in the dimer involves a near-resonant first overtone transition from the OH stretch to the bending mode. The energy transfer probability in the presence of the 1:2 resonance is 0.61 at 300 K. The bending mode then redistributes its energy to low-frequency intermolecular vibrations in a series of small excitation steps, with the pathway which results in the hydrogen-bonding modes gaining most of the available energy. The hydrogen bonding in ∼50% of the trajectories ruptures on vibrational excitation, leaving one quantum in the bend of the monomer fragment. In a small fraction of trajectories, the duration of collision is longer than 1 ps, during which the dimer and H₂O form a short-lived complex through a secondary hydrogen bond, which undergoes large amplitude oscillations.

Microhydration of PAH+ cations: evolution of hydration network in naphthalene+-(H2O) n clusters (n ≤ 5)

The interaction of polycyclic aromatic hydrocarbon molecules with water (H2O = W) is of fundamental importance in chemistry and biology. Herein, size-selected microhydrated naphthalene cation nanoclusters, Np+-W n (n ≤ 5), are characterized by infrared photodissociation (IRPD) spectroscopy in the C-H and O-H stretch range to follow the stepwise evolution of the hydration network around this prototypical PAH+ cation. The IRPD spectra are highly sensitive to the hydration structure and are analyzed by dispersion-corrected density functional theory calculations (B3LYP-D3/aug-cc-pVTZ) to determine the predominant structural isomers. For n = 1, W forms a bifurcated CH···O ionic hydrogen bond (H-bond) to two acidic CH protons of the bicyclic ring. For n ≥ 2, the formation of H-bonded solvent networks dominates over interior ion solvation, because of strong cooperativity in the former case. For n ≥ 3, cyclic W n solvent structures are attached to the CH protons of Np+. However, while for n = 3 the W3 ring binds in the CH···O plane to Np+, for n ≥ 4 the cyclic W n clusters are additionally stabilized by stacking interactions, leading to sandwich-type configurations. No intracluster proton transfer from Np+ to the W n solvent is observed in the studied size range (n ≤ 5), because of the high proton affinity of the naphthyl radical compared to W n . This is different from microhydrated benzene+ clusters, (Bz-W n )+, for which proton transfer is energetically favorable for n ≥ 4 due to the much lower proton affinity of the phenyl radical. Hence, because of the presence of polycyclic rings, the interaction of PAH+ cations with W is qualitatively different from that of monocyclic Bz+ with respect to interaction strength, structure of the hydration shell, and chemical reactivity. These differences are rationalized and quantified by quantum chemical analysis using the natural bond orbital (NBO) and noncovalent interaction (NCI) approaches.