Ultrafast Energy Equilibration in Hydrogen-Bonded Liquids

Interplay between anion-receptor and anion-solvent interactions in halide receptor complexes characterized with ultrafast infrared spectroscopies

The competition between host-guest binding and solvent interactions is a crucial factor in determining the binding affinities and selectivity of molecular receptor species. The interplay between these competing interactions, however, have been difficult to disentangle. In particular, the development of molecular-level descriptions of solute-solvent interactions remains a grand experimental challenge. Herein, we investigate the prototypical halide receptor meso-octamethylcalix[4]pyrrole (OMCP) complexed with either chloride or bromide anions in both dichloromethane (DCM) and chloroform (trichloromethane, TCM) solvent using ultrafast infrared transient absorption and 2D IR spectroscopies. OMCP·Br- complexes in both solvents display slower vibrational relaxation dynamics of the OMCP pyrrole NH stretches, consistent with weaker H-bonding interactions with OMCP compared to chloride and less efficient intermolecular relaxation to the solvent. Further, OMCP·Br- complexes show nearly static spectral diffusion dynamics compared to OMCP·Cl-, indicating larger structural fluctuations occur within chloride complexes. Importantly, distinct differences in the vibrational spectra and dynamics are observed between DCM and TCM solutions. The data are consistent with stronger and more perturbative solvent effects in TCM compared to DCM, despite DCM's larger dielectric constant and smaller reported OMCP·X- binding affinities. These differences are attributed to the presence of weak H-bond interactions between halides and TCM, in addition to competing interactions from the bulky tetrabutylammonium countercation. The data provide important experimental benchmarks for quantifying the role of solvent and countercation interactions in anion host-guest complexes.

The Role of Sum-Frequency Generation Spectroscopy in Understanding On-Surface Reactions and Dynamics in Atmospheric Model-Systems

Surfaces, both water/air and solid/water, play an important role in mediating a multitude of processes central to atmospheric chemistry, particularly in the aerosol phase. However, the study of both static and dynamic properties of surfaces is highly challenging from an experimental standpoint, leading to a lack of molecular level information about the processes that take place at these systems and how they differ from bulk. One of the few techniques that has been able to capture ultrafast surface phenomena is time-resolved sum-frequency generation (SFG) spectroscopy. Since it is both surface-specific and chemically sensitive, the extension of this spectroscopic technique to the time domain makes it possible to study dynamic processes on the femtosecond time scale. In this Perspective, we will explore recent advances made in the field both in terms of studying energy dissipation as well as chemical reactions and the role the surface geometry plays in these processes.

Deciphering Structural and Dynamical Properties of Hydrated Cobalt Porphyrins via Ab Initio Quantum Mechanical Charge Field Molecular Dynamics Simulation

The present study successfully implemented the ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism for the investigation of structural and dynamical properties of hydrated cobalt-porphyrin complexes. Considering the significance of cobalt ions in biological systems (for instance, vitamin B12), which reportedly incorporate cobalt ions in a d6, low spin, +3 state chelated in the corrin ring, an analog of porphyrin, the current study is focused on cobalt in the oxidation states +2 and +3 bound to the parent porphyrin lead structures embedded in an aqueous solution. These cobalt-porphyrin complexes were investigated in terms of their structural and dynamical properties at the quantum chemical level. The structural attributes of these hydrated complexes revealed the contrasting features of the water binding to these solutes, including a detailed evaluation of the associated dynamics. The study also yielded notable findings in regard to the respective electronic configurations vs coordination, which suggested that Co(II)-POR possesses a 5-fold square pyramidal coordination geometry in an aqueous solution containing the metal ion coordinating to four nitrogen atoms of the porphyrin ring and one axial water as the fifth ligand. On the other hand, high-spin Co(III)-POR was hypothesized to be more stable due to the smaller size-to-charge ratio of the cobalt ion, but the high-spin complex demonstrated unstable structural and dynamical behavior. However, the corresponding properties of the hydrated Co(III)_LS-POR revealed a stable structure in an aqueous solution, thus suggesting the Co(III) ion to be in a low-spin state when bound to the porphyrin ring. Moreover, the structural and dynamical data were augmented by computing the free energy of water binding to the cobalt ions and the solvent-accessible surface area, which provide further information on thermochemical properties of the metal-water interaction and the hydrogen bonding potential of the porphyrin ring in these hydrated systems.

Ultrafast vibrational dynamics of the free OD at the air/water interface: Negligible isotopic dilution effect but large isotope substitution effect

Vibrational relaxation dynamics of the OH stretch of water at the air/water interface has been a subject of intensive research, facilitated by recent developments in ultrafast interface-selective nonlinear spectroscopy. However, a reliable determination of the vibrational relaxation dynamics in the OD stretch region at the air/D₂O interface has not been yet achieved. Here, we report a study of the vibrational relaxation of the free OD carried out by time-resolved heterodyne-detected vibrational sum frequency generation spectroscopy. The results obtained with the aid of singular value decomposition analysis indicate that the vibrational relaxation (T₁) time of the free OD at the air/D₂O interface and air/isotopically diluted water (HOD-H₂O) interfaces show no detectable isotopic dilution effect within the experimental error, as in the case of the free OH in the OH stretch region. Thus, it is concluded that the relaxation of the excited free OH/OD predominantly proceeds with their reorientation, negating a major contribution of the intramolecular energy transfer. It is also shown that the T₁ time of the free OD is substantially longer than that of the free OH, further supporting the reorientation relaxation mechanism. The large difference in the T₁ time between the free OD and the free OH (factor of ∼2) may indicate the nuclear quantum effect on the diffusive reorientation of the free OD/OH because this difference is significantly larger than the value expected for a classical rotational motion.

Vibrational Dynamics of the Intramolecular H-Bond in Acetylacetone Investigated with Transient and 2D IR Spectroscopy

Acetylacetone (AcAc) has proven to be a fruitful but highly challenging model system for the experimental and computational interrogation of strong intramolecular hydrogen bonds. Key questions remain, however, regarding the identity of the minimum-energy structure of AcAc and the dynamics of intramolecular proton transfer. Here, we investigate the OH/OD stretch and bend regions of the enol tautomer of AcAc and its deuterated isotopologue with transient absorption and 2D IR spectroscopy. The OH bend region reveals a single dominant diagonal transition near 1625 cm^-1 with intense cross peaks to lower-frequency modes, demonstrating highly mixed fingerprint transitions that contain OH bend character. The anharmonic coupling of the OH bend results in a highly elongated OH bend excited-state absorption transition that indicates a large manifold of OH bend overtone/combination bands in the OH stretch region that leads to strong bend-stretch Fermi resonance interactions. The OH and OD stretch regions consist of broad ground-state bleach signals, but there is no clear evidence of ω₂₁ excited-state absorptions due to rapid population relaxation arising from strong intramolecular coupling to bending, fingerprint, and low-frequency H-bond modes. Orientational relaxation dynamics persist for timescales longer than the vibrational lifetimes, with polarization anisotropy components decaying within approximately 2 and 10 periods of the O-O oscillation for the OH and OD stretch, respectively. The significant isotopic dependence of the orientational dynamics is discussed in the context of intramolecular mode coupling, diffusional processes, and contributions from proton/deuteron transfer dynamics.

Least-Squares Fitting of Multidimensional Spectra to Kubo Line-Shape Models

We report a comprehensive study of the efficacy of least-squares fitting of multidimensional spectra to generalized Kubo line-shape models and introduce a novel least-squares fitting metric, termed the scale invariant gradient norm (SIGN), that enables a highly reliable and versatile algorithm. The precision of dephasing parameters is between 8× and 50× better for nonlinear model fitting compared to that for the centerline-slope (CLS) method, which effectively increases data acquisition efficiency by 1–2 orders of magnitude. Whereas the CLS method requires sequential fitting of both the nonlinear and linear spectra, our model fitting algorithm only requires nonlinear spectra but accurately predicts the linear spectrum. We show an experimental example in which the CLS time constants differ by 60% for independent measurements of the same system, while the Kubo time constants differ by only 10% for model fitting. This suggests that model fitting is a far more robust method of measuring spectral diffusion than the CLS method, which is more susceptible to structured residual signals that are not removable by pure solvent subtraction. Statistical analysis of the CLS method reveals a fundamental oversight in accounting for the propagation of uncertainty by Kubo time constants in the process of fitting to the linear absorption spectrum. A standalone desktop app and source code for the least-squares fitting algorithm are freely available, with example line-shape models and data. We have written the MATLAB source code in a generic framework where users may supply custom line-shape models. Using this application, a standard desktop fits a 12-parameter generalized Kubo model to a 10⁶ data-point spectrum in a few minutes.

Ab initio molecular dynamics study on energy relaxation path of hydrogen-bonded OH vibration in bulk water

The vibrational energy relaxation paths of hydrogen-bonded (H-bonded) OH excited in pure water and in isotopically diluted (deuterated) water are elucidated via non-equilibrium ab initio molecular dynamics (NE-AIMD) simulations. The present study extends the previous NE-AIMD simulation for the energy relaxation of an excited free OH vibration at an air/water interface [T. Ishiyama, J. Chem. Phys. 154, 104708 (2021)] to the energy relaxation of an excited H-bonded OH vibration in bulk water. The present simulation shows that the excited OH vibration in pure water dissipates its energy on a timescale of 0.1 ps, whereas that in deuterated water relaxes on a timescale of 0.7 ps, consistent with the experimental observations. To decompose these relaxation energies into the components due to intramolecular and intermolecular couplings, constraints are introduced on the vibrational modes except for the target path in the NE-AIMD simulation. In the case of pure water, 80% of the total relaxation is attributed to the pathway due to the resonant intermolecular OH⋯OH stretch coupling, and the remaining 17% and 3% are attributed to intramolecular couplings with the bend overtone and with the conjugate OH stretch, respectively. This result strongly supports a significant role for the Förster transfer mechanism of pure water due to the intermolecular dipole-dipole interactions. In the case of deuterated water, on the other hand, 36% of the total relaxation is due to the intermolecular stretch coupling, and all the remaining 64% arises from coupling with the intramolecular bend overtone.

Slow Proton Transfer in Nanoconfined Water

The transport of protons in nanoconfined environments, such as in nanochannels of biological or artificial proton conductive membranes, is essential to chemistry, biology, and nanotechnology. In water, proton diffusion occurs by hopping of protons between water molecules. This process involves the rearrangement of many hydrogen bonds and as such can be strongly affected by nanoconfinement. We study the vibrational and structural dynamics of hydrated protons in water nanodroplets stabilized by a cationic surfactant using polarization-resolved femtosecond infrared transient absorption spectroscopy. We determine the time scale of proton hopping in the center of the water nanodroplets from the dynamics of the anisotropy of the transient absorption signals. We find that in small nanodroplets with a diameter <4 nm, proton hopping is more than 10 times slower than in bulk water. Even in relatively large nanodroplets with a diameter of ∼7 nm, we find that the rate of proton hopping is slowed by ∼4 times compared with bulk water.

Interfacial Vibrational Dynamics of Ice Ih and Liquid Water

Insights into energy flow dynamics at ice surfaces are essential for understanding chemical dynamics relevant to atmospheric and geographical sciences. Here, employing ultrafast surface-specific spectroscopy, we report the interfacial vibrational dynamics of ice I_h. A comparison to liquid water surfaces reveals accelerated vibrational energy relaxation and dissipation at the ice surface for hydrogen-bonded OH groups. In contrast, free-OH groups sticking into the vapor phase exhibit substantially slower vibrational dynamics on ice. The acceleration and deceleration of vibrational dynamics of these different OH groups at the ice surface are attributed to enhanced intermolecular coupling and reduced rotational mobility, respectively. Our results highlight the unique properties of free-OH groups on ice, putatively linked to the high catalytic activities of ice surfaces.

Hidden Isolated OH at the Charged Hydrophobic Interface Revealed by Two-Dimensional Heterodyne-Detected VSFG Spectroscopy

Water around hydrophobic groups mediates hydrophobic interactions that play key roles in many chemical and biological processes. Thus, the molecular-level elucidation of the properties of water in the vicinity of hydrophobic groups is important. We report on the structure and dynamics of water at two oppositely charged hydrophobic ion/water interfaces, that is, the tetraphenylborate-ion (TPB^- )/water and tetraphenylarsonium-ion (TPA⁺ )/water interfaces, which are clarified by two-dimensional heterodyne-detected vibrational sum-frequency generation (2D HD-VSFG) spectroscopy. The obtained 2D HD-VSFG spectra of the anionic TPB^- interface reveal the existence of distinct π-hydrogen bonded OH groups in addition to the usual hydrogen-bonded OH groups, which are hidden in the steady-state spectrum. In contrast, 2D HD-VSFG spectra of the cationic TPA⁺ interface only show the presence of usual hydrogen-bonded OH groups. The present study demonstrates that the sign of the interfacial charge governs the structure and dynamics of water molecules that face the hydrophobic region.

Hydration of Closely Related Manganese and Magnesium Porphyrins in Aqueous Solutions: Ab Initio Quantum Mechanical Charge Field Molecular Dynamics Simulation Study

To the best of our knowledge, the current study based on ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) is the first to explore the difference in the hydration behavior between Mn(II)- and Mg(II)-associated porphyrins (Mn(II)-POR and Mg(II)-POR) in aqueous solution. The simulation study highlights similar and dissimilar characteristics of the structural, dynamical, and thermodynamical properties of these closely related metals bound to porphyrins in aqueous solution. The structural analysis is based on radial and angular distribution functions, coordination number distributions, and angular-radial distributions. Both hydrated systems demonstrate similar pentacoordinated structures formed via the axial coordination of one water molecule to the metal ion in addition to the four nitrogen atoms of the porphyrin ring. However, in the case of Mn(II)-POR, the formation of a distorted square pyramidal geometry was observed. It was envisaged as a weak coordination of the water molecule to the Mn(II) atom and thus higher atomic fluctuation for all atoms in contrast to that for the hydrated Mg(II)-POR. The dynamical data in terms of the mean residence times, velocity autocorrelation function, free energy, and other parameters revealed the difference in the metal binding effect because the Mn(II) atom was observed to inhibit H-bond formation more than the presence of Mg(II) atoms in the core of the porphyrin. The current study thus highlights the significant differences in the structural and dynamical properties of Mn(II)- and Mg(II)-associated porphyrin systems.

Energy Relaxation and Thermal Diffusion in Infrared Pump-Probe Spectroscopy of Hydrogen-Bonded Liquids

Infrared pump-probe spectroscopy provides detailed information about the dynamics of hydrogen-bonded liquids. Due to dissipation of the absorbed pump pulse energy, thermal equilibration dynamics also contributes to the observed signal. Disentangling this contribution from the molecular response remains a challenge. By performing non-equilibrium molecular dynamics simulations of liquid-deuterated methanol, we show that faster molecular vibrational relaxation and slower heat diffusion are decoupled and occur on different length scales. Transient structures of the hydrogen bonding network influence thermal relaxation by affecting thermal diffusivity over a length scale of several nanometers.

Ultrafast energy relaxation dynamics of amide I vibrations coupled with protein-bound water molecules

The influence of hydration water on the vibrational energy relaxation in a protein holds the key to understand ultrafast protein dynamics, but its detection is a major challenge. Here, we report measurements on the ultrafast vibrational dynamics of amide I vibrations of proteins at the lipid membrane/H₂O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy. We find that the relaxation time of the amide I mode shows a very strong dependence on the H₂O exposure, but not on the D₂O exposure. This observation indicates that the exposure of amide I bond to H₂O opens up a resonant relaxation channel and facilitates direct resonant vibrational energy transfer from the amide I mode to the H₂O bending mode. The protein backbone motions can thus be energetically coupled with protein-bound water molecules. Our findings highlight the influence of H₂O on the ultrafast structure dynamics of proteins.

Energy Relaxation and Structural Dynamics of Protons in Water/DMSO Mixtures

We investigate the structure and dynamics of proton solvation structures in mixed water/dimethyl sulfoxide (DMSO) solvents using two-color mid-infrared femtosecond pump–probe spectroscopy. At a water fraction below 20%, protons are mainly solvated as (DMSO-H)⁺ and (DMSO-H)⁺-H₂O structures. We find that excitation of the OH-stretch vibration of the proton in (DMSO-H)⁺-H₂O structures leads to an ultrafast contraction of the hydrogen bond between (DMSO-H)⁺ and H₂O. This excited state relaxes rapidly with T₁ = 95 ± 10 fs and leads in part to a strong local heating effect and in part to predissociation of the protonated cluster into (DMSO-H)⁺ and water monomers.

Aqueous solvation from the water perspective

The response of water re-solvating a charge-transfer dye (deprotonated Coumarin 343) after photoexcitation has been measured by means of transient THz spectroscopy. Two steps of increasing THz absorption are observed, a first ∼10 ps step on the time scale of Debye relaxation of bulk water and a much slower step on a 3.9 ns time scale, the latter of which reflecting heating of the bulk solution upon electronic relaxation of the dye molecules from the S₁ back into the S₀ state. As an additional reference experiment, the hydroxyl vibration of water has been excited directly by a short IR pulse, establishing that the THz signal measures an elevated temperature within ∼1 ps. This result shows that the first step upon dye excitation (10 ps) is not limited by the response time of the THz signal; it rather reflects the reorientation of water molecules in the solvation layer. The apparent discrepancy between the relatively slow reorientation time and the general notion that water is among the fastest solvents with a solvation time in the sub-picosecond regime is discussed. Furthermore, non-equilibrium molecular dynamics simulations have been performed, revealing a close-to-quantitative agreement with experiment, which allows one to disentangle the contribution of heating to the overall THz response from that of water orientation.

Persistent optical hole-burning spectroscopy of nano-confined dye molecules in liquid at room temperature: Spectral narrowing due to a glassy state and extraordinary relaxation in a nano-cage

Persistent optical hole-burning spectroscopy has been conducted for a dye molecule within a very small (∼1 nm) reverse micelle at room temperature. The spectra show a spectral narrowing due to site-selective excitation. This definitely demonstrates that the surroundings of the dye molecule are in a glassy state regardless of a solution at room temperature. On the other hand, the hole-burning spectra exhibit large shifts from excitation frequencies, and their positions are almost independent of excitation frequencies. The hole-burning spectra have been theoretically calculated by taking account of a vibronic absorption band of the dye molecule under the assumption that the surroundings of the dye molecule are in a glassy state. The calculated results agree with the experimental ones that were obtained for the dye molecule in a polymer glass for comparison, where it has been found that the ratio of hole-burning efficiencies of vibronic- to electronic-band excitations is quite high. On the other hand, the theoretical results do not explain the large spectral shift from the excitation frequency and small spectral narrowing observed in the hole-burning spectra measured for the dye-containing reverse micelle. It is thought that the spectral shift and broadening occur within the measurement time owing to the relaxation process of the surroundings that are hot with the thermal energy deposited by the dye molecule optically excited. Furthermore, the relaxation should be temporary because the cooling of the inside of the reverse micelle takes place with the dissipation of the excess thermal energy to the outer oil solvent, and so the surroundings of the dye molecule return to the glassy state and do not attain the thermal equilibrium. These results suggest that a very small reverse micelle provides a unique reaction field in which the diffusional motion can be controlled by light in a glassy state.

Delocalization and stretch-bend mixing of the HOH bend in liquid water

Liquid water's rich sub-picosecond vibrational dynamics arise from the interplay of different high- and low-frequency modes evolving in a strong yet fluctuating hydrogen bond network. Recent studies of the OH stretching excitations of H₂O indicate that they are delocalized over several molecules, raising questions about whether the bending vibrations are similarly delocalized. In this paper, we take advantage of an improved 50 fs time-resolution and broadband infrared (IR) spectroscopy to interrogate the 2D IR lineshape and spectral dynamics of the HOH bending vibration of liquid H₂O. Indications of strong bend-stretch coupling are observed in early time 2D IR spectra through a broad excited state absorption that extends from 1500 cm^-1 to beyond 1900 cm^-1, which corresponds to transitions from the bend to the bend overtone and OH stretching band between 3150 and 3550 cm^-1. Pump-probe measurements reveal a fast 180 fs vibrational relaxation time, which results in a hot-ground state spectrum that is the same as observed for water IR excitation at any other frequency. The fastest dynamical time scale is 80 fs for the polarization anisotropy decay, providing evidence for the delocalized or excitonic character of the bend. Normal mode analysis conducted on water clusters extracted from molecular dynamics simulations corroborate significant stretch-bend mixing and indicate delocalization of δ_HOH on 2-7 water molecules.

Structure from Dynamics: Vibrational Dynamics of Interfacial Water as a Probe of Aqueous Heterogeneity

The structural heterogeneity of water at various interfaces can be revealed by time-resolved sum-frequency generation spectroscopy. The vibrational dynamics of the O–H stretch vibration of interfacial water can reflect structural variations. Specifically, the vibrational lifetime is typically found to increase with increasing frequency of the O–H stretch vibration, which can report on the hydrogen-bonding heterogeneity of water. We compare and contrast vibrational dynamics of water in contact with various surfaces, including vapor, biomolecules, and solid interfaces. The results reveal that variations in the vibrational lifetime with vibrational frequency are very typical, and can frequently be accounted for by the bulk-like heterogeneous response of interfacial water. Specific interfaces exist, however, for which the behavior is less straightforward. These insights into the heterogeneity of interfacial water thus obtained contribute to a better understanding of complex phenomena taking place at aqueous interfaces, such as photocatalytic reactions and protein folding.

Hydration of iron-porphyrins: ab initio quantum mechanical charge field molecular dynamics simulation study

The ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) simulation approach was successfully applied to Fe²⁺-P and Fe³⁺-P in water to evaluate their structural, dynamical and energetic properties. Based on the structural data, it was found that Fe²⁺-P accommodates one water molecule in the first coordination sphere of the Fe²⁺ ion including the four nitrogen atoms of the porphyrin system coordinating with central metal species. On the other hand, two water molecules were coordinated to Fe³⁺-P, thus forming a hexa-coordinated species. Comparison of dynamical properties such as the vibrational power spectrum and ligand mean residence times to other metal-free porphyrin systems demonstrate the ions' influence on the hydration structure, enabling a characterisation of the strong interaction of the ions which greatly reduces the hydrogen bonding potential of the complex. The association of water molecules with the metal ions in both solutes was quantified by computing the free energy of binding obtained via the potential of mean force. This further confirmed the strong association of water to the metal ions which was conversely weak as inferred from the energetic data for the Fe²⁺-P system.

Zinc- and copper-porphyrins in aqueous solution - two similar complexes with strongly contrasting hydration

We present detailed analysis of the hydration behavior of zinc and copper bound porphyrins treated via ab initio quantum mechanical charge field molecular dynamics which agrees well with available experimental data. The computed metal-water coordination in the case of zinc bound porphyrin demonstrates a strong association of water with zinc compared to the copper-water interaction which correlates well with the calculated free energy of binding. The H-bond dynamics in these hydrated systems yield weaker H-bond interactions as compared to that observed in the case of metal-free porphyrin; nevertheless, the effect of metal association with porphyrin resulted in shifts in the vibrational frequencies. These characteristic data suggest a contrasting behavior between these metalloporphyrins in solution which could also serve to correlate with the properties of biological systems.

Quantum Interference in the Vibrational Relaxation of the O-H Stretch Overtone of Liquid H2O

Using femtosecond two-color infrared pump-probe spectroscopy, we study the vibrational relaxation of the O-H stretch vibrations of liquid H2O after excitation of the overtone transition. The overtone transition has its maximum at 6900 cm(-1) (1.45 μm), which is a relatively high frequency in view of the central frequency of 3400 cm(-1) of the fundamental transition. The excitation of the overtone leads to a transient induced absorption of two-exciton states of the O-H stretch vibrations. When the overtone excitation frequency is tuned from 6600 to 7200 cm(-1), the vibrational relaxation time constant of the two-exciton states increases from 400 ± 50 fs to 540 ± 40 fs. These values define a limited range of relatively long relaxation time constants compared to the range of relaxation time constants of 250-550 fs that we recently observed for the one-exciton O-H stretch vibrational state of liquid H2O ( S. T. van der Post et al., Nature Comm. 2015 , 6 , 8384 ). We explain the high central frequency and the limited range of relatively long relaxation time constants of the overtone transition from the destructive quantum interference of the mechanical and electrical anharmonic contributions to the overtone transition probability. As a result of this destructive interference, the overtone transition of liquid H2O is dominated by molecules of which the O-H groups donate relatively weak hydrogen bonds to other H2O molecules.

The Opposite Effect of Metal Ions on Short-/Long-Range Water Structure: A Multiple Characterization Study

Inorganic electrolyte solutions are very important in our society as they dominate many biochemical and geochemical processes. Herein, an in-depth study was performed to illustrate the ion-induced effect on water structure by coupling NMR, viscometer, Raman and Molecular Dynamic (MD) simulations. The NMR coefficient (B_NMR) and diffusion coefficient (D) from NMR, and viscosity coefficient (B_vis) from a viscometer all proved that dissolved metal ions are capable of enhancing the association degree of adjacent water molecules, and the impact on water structure decreased in the order of Cr³⁺ > Fe³⁺ > Cu²⁺ > Zn²⁺. This regularity was further evidenced by Raman analysis; however, the deconvoluted Raman spectrum indicated the decrease in high association water with salt concentration and the increase in low association water before 200 mmol·L⁻¹. By virtue of MD simulations, the opposite changing manner proved to be the result of the opposite effect on short-/long-range water structure induced by metal ions. Our results may help to explain specific protein denaturation induced by metal ions.

Vibrational dynamics of aqueous hydroxide solutions probed using broadband 2DIR spectroscopy

We employed ultrafast transient absorption and broadband 2DIR spectroscopy to study the vibrational dynamics of aqueous hydroxide solutions by exciting the O-H stretch vibrations of the strongly hydrogen-bonded hydroxide solvation shell water and probing the continuum absorption of the solvated ion between 1500 and 3800 cm(-1). We observe rapid vibrational relaxation processes on 150-250 fs time scales across the entire probed spectral region as well as slower vibrational dynamics on 1-2 ps time scales. Furthermore, the O-H stretch excitation loses its frequency memory in 180 fs, and vibrational energy exchange between bulk-like water vibrations and hydroxide-associated water vibrations occurs in ∼200 fs. The fast dynamics in this system originate in strong nonlinear coupling between intra- and intermolecular vibrations and are explained in terms of non-adiabatic vibrational relaxation. These measurements indicate that the vibrational dynamics of the aqueous hydroxide complex are faster than the time scales reported for long-range transport of protons in aqueous hydroxide solutions.

Ultrafast 2D IR spectroscopy of the excess proton in liquid water

Despite decades of study, the structures adopted to accommodate an excess proton in water and the mechanism by which they interconvert remain elusive. We used ultrafast two-dimensional infrared (2D IR) spectroscopy to investigate protons in aqueous hydrochloric acid solutions. By exciting O–H stretching vibrations and detecting the spectral response throughout the mid-IR region, we observed the interaction between the stretching and bending vibrations characteristic of the flanking waters of the Zundel complex, [H(H2O)2]+, at 3200 and 1760 cm−1, respectively. From time-dependent shifts of the stretch-bend cross peak, we determined a lower limit on the lifetime of this complex of 480 femtoseconds. These results suggest a key role for the Zundel complex in aqueous proton transfer.

Strong frequency dependence of vibrational relaxation in bulk and surface water reveals sub-picosecond structural heterogeneity

Because of strong hydrogen bonding in liquid water, intermolecular interactions between water molecules are highly delocalized. Previous two-dimensional infrared spectroscopy experiments have indicated that this delocalization smears out the structural heterogeneity of neat H2O. Here we report on a systematic investigation of the ultrafast vibrational relaxation of bulk and interfacial water using time-resolved infrared and sum-frequency generation spectroscopies. These experiments reveal a remarkably strong dependence of the vibrational relaxation time on the frequency of the OH stretching vibration of liquid water in the bulk and at the air/water interface. For bulk water, the vibrational relaxation time increases continuously from 250 to 550 fs when the frequency is increased from 3,100 to 3,700 cm(-1). For hydrogen-bonded water at the air/water interface, the frequency dependence is even stronger. These results directly demonstrate that liquid water possesses substantial structural heterogeneity, both in the bulk and at the surface.

Polarization-modulation setup for ultrafast infrared anisotropy experiments to study liquid dynamics

An infrared pump-probe setup using rapid polarization modulation has been developed to perform time-resolved vibrational anisotropy measurements. A photo-elastic modulator is used as a rapidly switchable half-wave plate, enabling the measurement of transient absorptions for parallel and perpendicular polarizations of the pump and probe pulses on a shot-to-shot basis. In this way, infrared intensity fluctuations are nearly completely canceled, significantly enhancing the accuracy of the transient-anisotropy measurement. The method is tested on the OD-stretch vibration of HDO in H₂O, for which the signal-to-noise ratio is found to be 4 times better than with conventional methods.

Hydration properties determining the reactivity of nitrite in aqueous solution

The knowledge of the hydration properties of the nitrite ion is key to understanding its reaction mechanism controlled by solvent effects. Here, ab initio quantum mechanical charge field molecular dynamics was performed to obtain the structural and dynamical properties of the hydration shell in an aqueous solution of nitrite ions, elucidated by data analysis using a molecular approach and an extended quantitative analysis of all superimposed trajectories with three-dimensional alignment (density map). The pattern of the power spectra corresponded to the experimental data, indicating the suitability of the Hartree-Fock method coupled with double-ζ plus polarization and diffuse functional basis sets to study this system. The density maps revealed the structure of the hydration shell, that presented a higher density in the N-O bond direction than in the axis vertical to the molecular plane, whereas the atomic and molecular radial distribution functions provided vague information. The number of actual contacts indicated 4.6 water molecules interacting with a nitrite ion, and 1.5 extra water molecules located in the molecular hydration shell, forming a H-bonding network with the bulk water. The mean residence times for the water ligands designated the strength of the hydration spheres for the oxygen sites, whilst the results for the nitrogen sites over-estimated the number of water molecules from other sites and indicated a weak structure. These results show the influence of the water molecules surrounding the nitrite ion creating an anisotropic hydration shell, suggesting that the reactive sites are situated above and below the molecular plane with a lower water density.

Changes of water hydrogen bond network with different externalities

It is crucial to uncover the mystery of water cluster and structural motif to have an insight into the abundant anomalies bound to water. In this context, the analysis of influence factors is an alternative way to shed light on the nature of water clusters. Water structure has been tentatively explained within different frameworks of structural models. Based on comprehensive analysis and summary of the studies on the response of water to four externalities (i.e., temperature, pressure, solutes and external fields), the changing trends of water structure and a deduced intrinsic structural motif are put forward in this work. The variations in physicochemical and biological effects of water induced by each externality are also discussed to emphasize the role of water in our daily life. On this basis, the underlying problems that need to be further studied are formulated by pointing out the limitations attached to current study techniques and to outline prominent studies that have come up recently.

Hydrogen bond dynamics in primary alcohols: a femtosecond infrared study

Hydrogen-bonded liquids are excellent solvents, in part due to the highly dynamic character of the directional interaction associated with the hydrogen bond. Here we study the vibrational and reorientational dynamics of deuterated hydroxyl groups in various primary alcohols using polarization-resolved femtosecond infrared spectroscopy. We show that the relaxation of the OD stretch vibration is similar for ethanol and its higher homologues (∼0.9 ps), while it is appreciably faster for methanol (∼0.75 ps). The fast relaxation for methanol is attributed to strong coupling of the OD stretch vibration to the overtone of the CH3 rocking mode. Subsequent to excited state relaxation, the dissipation of the excess energy leads to structural relaxation of the alcohol liquid structure. We show that this relaxation of the H-bonded network depends on the alkyl chain length. We find that the anisotropy of the excitation decays by both thermal diffusion from excited OD groups to nonexcited molecules and reorientational motion. The reorientation is described well by a model employing two relaxation times that increase linearly with increasing alcohol size. The short reorientation time is assigned to the partial reorientation of molecules within the alcohol cluster, while the long reorientation times can be attributed to breaking and reforming of hydrogen bonds.

Probing ultrafast temperature changes of aqueous solutions with coherent terahertz pulses

We introduce an infrared pump-terahertz probe technique to measure the thermalization dynamics of aqueous solutions with a time resolution <200  fs. This technique makes use of the sensitivity of the terahertz absorption to the temperature of the hydrogen bond network. The thermalization dynamics of different aqueous solutions are measured and compared to the dynamics inferred from ultrafast infrared pump-infrared probe measurements on the intramolecular stretch vibration of water. This technique can shed new light on important aspects of energy transfer and heat dynamics and is applicable to a wide range of systems.