Ultrafast Vibrational Dynamics of Hydrogen Bonds in the Condensed Phase

Observing vibronic coupling in a strongly hydrogen bonded system with coherent multidimensional vibrational-electronic spectroscopy

The coupled structural and electronic parameters of intramolecular hydrogen bonding play an important role in ultrafast chemical reactions, such as proton transfer processes. We perform one- and two-dimensional vibrational-electronic (1D and 2D VE) spectroscopy experiments to understand the couplings between vibrational and electronic coordinates in 10-Hydroxybenzo[h]quinoline, an ultrafast proton transfer system. The experiments reveal that the OH stretch (νOH) is strongly coupled to the electronic excitation, and Fourier analysis of the 1D data shows coherent oscillations from the low frequency backbone vibrational modes coupled to the νOH mode, resulting in an electronically detected vibronic signal. In-plane low-frequency vibrations at 242 and 386 cm-1 change the hydrogen bond distance and modulate the observed electronic signal in the polarization-selective 1D VE experiment through orientation-dependent coupling with the νOH mode. Resolution of the excitation frequency axis with 2D VE experiments reveals that excitation frequency, detection frequency, and experimental delay affect the frequency and strength of the vibronic transitions observed. Our results demonstrate evidence of direct coupling of the high frequency νOH mode with the S1 ← S0 electronic transition in 10-Hydroxybenzo[h]quinoline (HBQ), and orientation-dependent couplings of the low-frequency 242 and 386 cm-1 modes to the νOH mode and the electronic transition. This demonstration of multidimensional VE spectroscopy on HBQ reveals the potential of using 1D and 2D VE spectroscopy to develop a quantitative understanding of the role of vibronic coupling in hydrogen bonding and ultrafast proton transfer for complex systems.

Hydration structure and dynamics of phosphoric acid and its anions-Ultrafast 2D-IR spectroscopy and ab initio molecular dynamics simulations

The hydration shells of phosphate ions and phosphate groups of nucleotides and phospholipid membranes display markedly different structures and hydrogen-bond strengths. Understanding phosphate hydration requires insight into the spatial arrangements of water molecules around phosphates and in thermally activated structure fluctuations on ultrafast time scales. Femtosecond two-dimensional infrared spectroscopy of phosphate vibrations, particularly asymmetric stretching vibrations between 1000 and 1200 cm-1, and ab initio molecular dynamics (AIMD) simulations are combined to map and characterize dynamic local hydration structures and phosphate-water interactions. Phosphoric acid H3PO4 and its anions H2PO4-, HPO42-, and PO43- are studied in aqueous environments of different pH value. The hydration shells of phosphates providing OH donor groups in hydrogen bonds with the first water layer undergo ultrafast structural fluctuations, which induce a pronounced spectral diffusion of vibrational excitations on a sub-300 fs time scale. With a decreasing number of phosphate OH groups, the hydration shell becomes more ordered and rigid. The 2D-IR line shapes observed with hydrated PO43- ions display a pronounced inhomogeneous broadening, reflecting a distribution of hydration geometries without fast equilibration. The AIMD simulations allow for an in-depth characterization of the hydration geometries with different numbers of water molecules in the first hydration layer and different correlation functions of the fluctuating electric field that the water environment exerts on the vibrational phosphate oscillators.

Molecular dynamics simulation study of water structure and dynamics on the gold electrode surface with adsorbed 4-mercaptobenzonitrile

Understanding water dynamics at charged interfaces is of great importance in various fields, such as catalysis, biomedical processes, and solar cell materials. In this study, we implemented molecular dynamics simulations of a system of pure water interfaced with Au electrodes, on one side of which 4-mercaptobenzonitrile (4-MBN) molecules are adsorbed. We calculated time correlation functions of various dynamic quantities, such as the hydrogen bond status of the N atom of the adsorbed 4-MBN molecules, the rotational motion of the water OH bond, hydrogen bonds between 4-MBN and water, and hydrogen bonds between water molecules in the interface region. Using the Luzar-Chandler model, we analyzed the hydrogen bond dynamics between a 4-MBN and a water molecule. The dynamic quantities we calculated can be divided into two categories: those related to the collective behavior of interfacial water molecules and the H-bond interaction between a water molecule and the CN group of 4-MBN. We found that these two categories of dynamic quantities exhibit opposite trends in response to applied potentials on the Au electrode. We anticipate that the present work will help improve our understanding of the interfacial dynamics of water in various electrolyte systems.

Unified picture of vibrational relaxation of OH stretch at the air/water interface

The elucidation of the energy dissipation process is crucial for understanding various phenomena occurring in nature. Yet, the vibrational relaxation and its timescale at the water interface, where the hydrogen-bonding network is truncated, are not well understood and are still under debate. In the present study, we focus on the OH stretch of interfacial water at the air/water interface and investigate its vibrational relaxation by femtosecond time-resolved, heterodyne-detected vibrational sum-frequency generation (TR-HD-VSFG) spectroscopy. The temporal change of the vibrationally excited hydrogen-bonded (HB) OH stretch band (ν=1→2 transition) is measured, enabling us to determine reliable vibrational relaxation (T1) time. The T1 times obtained with direct excitations of HB OH stretch are 0.2-0.4 ps, which are similar to the T1 time in bulk water and do not noticeably change with the excitation frequency. It suggests that vibrational relaxation of the interfacial HB OH proceeds predominantly with the intramolecular relaxation mechanism as in the case of bulk water. The delayed rise and following decay of the excited-state HB OH band are observed with excitation of free OH stretch, indicating conversion from excited free OH to excited HB OH (~0.9 ps) followed by relaxation to low-frequency vibrations (~0.3 ps). This study provides a complete set of the T1 time of the interfacial OH stretch and presents a unified picture of its vibrational relaxation at the air/water interface.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

The Role of H-Bonds in the Excited-State Properties of Multichromophoric Systems: Static and Dynamic Aspects

Given their importance, hydrogen bonds (H-bonds) have been the subject of intense investigation since their discovery. Indeed, H-bonds play a fundamental role in determining the structure, the electronic properties, and the dynamics of complex systems, including biologically relevant materials such as DNA and proteins. While H-bonds have been largely investigated for systems in their electronic ground state, fewer studies have focused on how the presence of H-bonds could affect the static and dynamic properties of electronic excited states. This review presents an overview of the more relevant progress in studying the role of H-bond interactions in modulating excited-state features in multichromophoric biomimetic complex systems. The most promising spectroscopic techniques that can be used for investigating the H-bond effects in excited states and for characterizing the ultrafast processes associated with their dynamics are briefly summarized. Then, experimental insights into the modulation of the electronic properties resulting from the presence of H-bond interactions are provided, and the role of the H-bond in tuning the excited-state dynamics and the related photophysical processes is discussed.

Direct detection of photo-induced reactions by IR: from Brook rearrangement to photo-catalysis

In situ IR detection of photoreactions induced by the light of LEDs at appropriate wavelengths provides a simple, cost-effective, and versatile method to get insight into mechanistic details. In particular, conversions of functional groups can be selectively followed. Overlapping UV-Vis bands or fluorescence from the reactants and products and the incident light do not obstruct IR detection. Compared with in situ photo-NMR, our setup does not require tedious sample preparation (optical fibers) and offers a selective detection of reactions, even at positions where ¹H-NMR lines overlap or ¹H resonances are not clear-cut. We illustrate the applicability of our setup following the photo-Brook rearrangement of (adamant-1-yl-carbonyl)-tris(trimethylsilyl)silane, address photo-induced α-bond cleavage (1-hydroxycyclohexyl phenyl ketone), study photoreduction using tris(bipyridine)ruthenium(II), investigate photo-oxygenation of double bonds with molecular oxygen and the fluorescent 2,4,6-triphenylpyrylium photocatalyst, and address photo-polymerization. With the LED/FT-IR combination, reactions can be qualitatively followed in fluid solution, (highly) viscous environments, and in the solid state. Viscosity changes during the reaction (e.g., during a polymerization) do not obstruct the method.

Water model for hydrophobic cavities: structure and energy from quantum-chemical calculations

This ab initio study aims to design a series of large water clusters having a hollow clathrate-like cage able to host hydrophobic solutes of various sizes. Starting from the (H₂O)_n (n = 18, 20, 24 and 28) hollow cages, water layers have been added in a stepwise manner in order to model the configuration of water molecules beyond the primary shell. The large (H₂O)₁₀₀, (H₂O)₁₂₀ and (H₂O)₁₄₀ clusters complete the hydrogen bonding network of the cage with optimal and regular tiling of the do-, tetra-decahedron and hexa-decahedron, respectively. This study is corroborated by an investigation of dense water clusters up to the (H₂O)₁₂₃ one, being highly consistent with experimental data on ice concerning the electronic and zero-point energies for aggregate formation at 0 K and enthalpy and entropy at 273 K. The cavity creation profoundly alters the orientation of water molecules compared with those found in dense clusters. Nevertheless, such a large reorganization is necessary to maximize the water-water attraction by making it similar to the one found in dense clusters. The cage formation is an endothermic process; however, the computed values are large compared with previous reports for hydrocarbon aqueous solutions. Larger clusters are required for a more fruitful comparison.

Vibrational relaxation by methylated xanthines in solution: Insights from 2D IR spectroscopy and calculations

Two-dimensional infrared (2D IR) spectroscopy, infrared pump-infrared probe spectroscopy, and density functional theory calculations were used to study vibrational relaxation by ring and carbonyl stretching modes in a series of methylated xanthine derivatives in acetonitrile and deuterium oxide (heavy water). Isotropic signals from the excited symmetric and asymmetric carbonyl stretch modes decay biexponentially in both solvents. Coherent energy transfer between the symmetric and asymmetric carbonyl stretching modes gives rise to a quantum beat in the time-dependent anisotropy signals. The damping time of the coherent oscillation agrees with the fast decay component of the carbonyl bleach recovery signals, indicating that this time constant reflects intramolecular vibrational redistribution (IVR) to other solute modes. Despite their similar frequencies, the excited ring modes decay monoexponentially with a time constant that matches the slow decay component of the carbonyl modes. The slow decay times, which are faster in heavy water than in acetonitrile, approximately match the ones observed in previous UV pump-IR probe measurements on the same compounds. The slow component is assigned to intermolecular energy transfer to solvent bath modes from low-frequency solute modes, which are populated by IVR and are anharmonically coupled to the carbonyl and ring stretch modes. 2D IR measurements indicate that the carbonyl stretching modes are weakly coupled to the delocalized ring modes, resulting in slow exchange that cannot explain the common solvent-dependence. IVR is suggested to occur at different rates for the carbonyl vs ring modes due to differences in mode-specific couplings and not to differences in the density of accessible states.

Low-temperature vibronic spectroscopy of condensed chromophore exhibiting inhomogeneous distribution of vibrational frequencies in a mixed quantum-classical environment

This work has been motivated by the recent paper by the author [M. Toutounji, Phys. Chem. Chem. Phys., 2021, 23, 21981] whereby a mixed quantum-classical Liouville equation was used to probe the spectroscopy and dynamics of a spin-boson system. A mixed quantum-classical Liouville equation treats the system of interest quantum mechanically, the bath classically, and the coupling term mixed quantum-classical mechanically. This paper offers a two-fold advantage: correcting the treatment of the electronic transition decay (width in frequency domain) and assessing the local heterogeneous vibrational structure. The homogeneous linear absorption spectrum of a chromophore embedded in a mixed quantum-classical environment at low temperature is composed of a sharp peak called a zero-phonon line (ZPL) and a broad phonon sideband (PSB), whereby the ZPL and the PSB are assimilated by a Lorentzian function and Voigt profiles, respectively. The PSB, in this case, is characterized by a local heterogeneous structure due to a dispersive medium of vibrations, modeled by vibrational Gaussian distributions to represent the arising inhomogeneous broadening and Lorentzians to model the homogeneous vibrations. This description seems to model proteins and amorphous solids exhibiting a local heterogeneous structure as both electronic and vibrational inhomogeneous broadening seems to be large in these media. This work provides a derivation of linear absorption lineshape and vibronic transition dipole moment time correlation functions, both of which account for pure electronic dephasing (ZPL width) the Voigt profile description of the phonon profiles (PSB) in dispersive media.

Van der Waals interactions regulating the hydration of 2-methacryloyloxyethyl phosphorylcholine, the constructing monomer of biocompatible polymers

Van der Waals (VDW) interactions provide fantastic properties for biological systems that function at room temperature. The VDW interaction, which primarily contributes to weak hydrogen bonding, is expected to play a key role in regulating hydrophobic hydration to express the biologically inert biocompatible function of polymerized MPCs (2-methacryloyloxyethyl phosphorylcholine). This report explores at the molecular level the biologically inert function of polymerized MPCs through an array of vibrational spectroscopic and computational characterization of MPC monomers, as temperature-dependent change of intramolecular weak hydrogen bonding. Synchrotron Fourier transform infrared microspectroscopy and terahertz time-domain spectroscopy were used to investigate temperature-dependent spectral changes in the low frequency vibrations of the MPC over the temperature range from cryogenic to room temperature, and the results were analysed by highly reliable well-established density functional theory (DFT) calculations. Complicated spectral features in the low frequency energy region and the uncertain conformations of the MPC in the amorphous powder state are clearly resolved under a polarizable continuum model and dispersion correction to pure DFT calculations.

Charge Gradients around Dendritic Voids Cause Nanoscale Inhomogeneities in Liquid Water

Water is the matrix of life and serves as a solvent for numerous physical and chemical processes. The origins of the nature of inhomogeneities that exist in liquid water and the time scales over which they occur remains an open question. Here, we report femtosecond elastic second harmonic scattering (fs-ESHS) of liquid water in comparison to an isotropic liquid (CCl₄) and show that water is indeed a nonuniform liquid. The coherent fs-ESHS intensity was interpreted, using molecular dynamics simulations, as arising from charge density fluctuations with enhanced nanoscale polarizabilities around transient voids having an average lifetime of 300 fs. Although voids were also present in CCl₄, they were not characterized by hydrogen bond defects and did not show strong polarizability fluctuations, leading to fs-ESHS of an isotropic liquid. The voids increased in number at higher temperatures above room temperature, in agreement with the fs-ESHS results.

From Networked to Isolated: Observing Water Hydrogen Bonds in Concentrated Electrolytes with Two-Dimensional Infrared Spectroscopy

Superconcentrated electrolytes have emerged as a promising class of materials for energy storage devices, with evidence that high voltage performance is possible even with water as the solvent. Here, we study the changes in the water hydrogen bonding network induced by the dissolution of lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) in concentrations ranging from the dilute to the superconcentrated regimes. Using time-resolved two-dimensional infrared spectroscopy, we observe the progressive disruption of the water-water hydrogen bond network and the appearance of isolated water molecules interacting only with ions, which can be identified and spectroscopically isolated through the intermolecular cross-peaks between the water and the TFSI^- ions. Analyzing the vibrational relaxation of excitations of the H₂O stretching mode, we observe a transition in the dominant relaxation path as the bulk-like water vanishes and is replaced by ion-solvation water with the rapid single-step relaxation of delocalized stretching vibrations into the low frequency modes being replaced by multistep relaxation through the intramolecular H₂O bend and into the TFSI^- high frequency modes prior to relaxing to the low frequency structural degrees of freedom. These results definitively demonstrate the absence of vibrationally bulk-like water in the presence of high concentrations of LiTFSI and especially in the superconcentrated regime, while additionally revealing aspects of the water hydrogen bond network that have been difficult to discern from the vibrational spectroscopy of the neat liquid.

Field-Induced Electron Generation in Water: Solvation Dynamics and Many-Body Interactions

The solvated electron represents an elementary quantum system in a liquid environment. Electrons solvated in water have raised strong interest because of their prototypical properties, their role in radiation chemistry, and their relevance for charge separation and transport. Nonequilibrium dynamics of photogenerated electrons in water occur on ultrafast time scales and include charge transfer, localization, and energy dissipation processes. We present new insight into the role of fluctuating electric fields of the liquid for generating electrons in the presence of an external terahertz field and address polaronic many-body properties of solvated electrons. This Perspective combines a review of recent results from experiment and theory with a discussion of basic electric interactions of electrons in water.

Characterization of Acetonitrile Isotopologues as Vibrational Probes of Electrolytes

Acetonitrile has emerged as a solvent candidate for novel electrolyte formulations in metal-ion batteries and supercapacitors. It features a bright local C≡N stretch vibrational mode whose infrared (IR) signature is sensitive to battery-relevant cations (Li⁺, Mg²⁺, Zn²⁺, Ca²⁺) both in pure form and in the presence of water admixture across a full possible range of concentrations from the dilute to the superconcentrated regime. Stationary and time-resolved IR spectroscopy thus emerges as a natural tool to study site-specific intermolecular interactions from the solvent perspective without introducing an extrinsic probe that perturbs solution morphology and may not represent the intrinsic dynamics in these electrolytes. The metal-coordinated acetonitrile, water-separated metal-acetonitrile pair, and free solvent each have a distinct vibrational signature that allows their unambiguous differentiation. The IR band frequency of the metal-coordinated acetonitrile depends on the ion charge density. To study the ion transport dynamics, it is necessary to differentiate energy-transfer processes from structural interconversions in these electrolytes. Isotope labeling the solvent is a necessary prerequisite to separate these processes. We discuss the design principles and choice of the CD₃¹³CN label and characterize its vibrational spectroscopy in these electrolytes. The Fermi resonance between ¹³C≡N and C-D stretches complicates the spectral response but does not prevent its effective utilization. Time-resolved two-dimensional (2D) IR spectroscopy can be performed on a mixture of acetonitrile isotopologues and much can be learned about the structural dynamics of various species in these formulations.

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.

Two-dimensional Infrared Spectroscopy Reveals Better Insights of Structure and Dynamics of Protein

Proteins play an important role in biological and biochemical processes taking place in the living system. To uncover these fundamental processes of the living system, it is an absolutely necessary task to understand the structure and dynamics of the protein. Vibrational spectroscopy is an established tool to explore protein structure and dynamics. In particular, two-dimensional infrared (2DIR) spectroscopy has already proven its versatility to explore the protein structure and its ultrafast dynamics, and it has essentially unprecedented time resolutions to observe the vibrational dynamics of the protein. Providing several examples from our theoretical and experimental efforts, it is established here that two-dimensional vibrational spectroscopy provides exceptionally more information than one-dimensional vibrational spectroscopy. The structural information of the protein is encoded in the position, shape, and strength of the peak in 2DIR spectra. The time evolution of the 2DIR spectra allows for the visualisation of molecular motions.

Energy relaxation dynamics of hydrogen-bonded OH vibration conjugated with free OH bond at an air/water interface

Vibrational energy relaxation dynamics of the excited hydrogen-bonded (H-bonded) OH conjugated with free OH (OD) at an air/water (for both pure water and isotopically diluted water) interface are elucidated via non-equilibrium ab initio molecular dynamics (NE-AIMD) simulations. The calculated results are compared with those of the excited H-bonded OH in bulk liquid water reported previously. In the case of pure water, the relaxation timescale (vibrational lifetime) of the excited H-bonded OH at the interface is T₁ = 0.13 ps, which is slightly larger than that in the bulk (T₁ = 0.11 ps). Conversely, in the case of isotopically diluted water, the relaxation timescale of T₁ = 0.74 ps in the bulk decreases to T₁ = 0.26 ps at the interface, suggesting that the relaxation dynamics of the H-bonded OH are strongly dependent on the surrounding H-bond environments particularly for the isotopically diluted conditions. The relaxation paths and their rates are estimated by introducing certain constraints on the vibrational modes except for the target path in the NE-AIMD simulation to decompose the total energy relaxation rate into contributions to possible relaxation pathways. It is found that the main relaxation pathway in the case of pure water is due to intermolecular OH⋯OH vibrational coupling, which is similar to the relaxation in the bulk. In the case of isotopically diluted water, the main pathway is due to intramolecular stretch and bend couplings, which show more efficient relaxation than in the bulk because of strong H-bonding interactions specific to the air/water interface.

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.

Infrared spectra of hydrogen bond network in lamellar perfluorocarboxylic acid monohydrates

The infrared spectra of the long-chain perfluorocarboxylic acid monohydrates differ markedly from those of the anhydrous dimers. Consequently, the structure of the solid perfluorocarboxylic acid monohydrates must differ from any known dimer-containing carboxylic acid crystals. Consideration of the significant features of the infrared spectra of the long-chain perfluorocarboxylic acid monohydrates, supplemented by their Raman spectra, and comparison with the spectra of auxiliary substances have led us to conclude that the rather strong neutral carboxyl-hydroxyl to water bonding can best explain the observations. The infrared spectra indicate the presence of fairly short hydrogen bonds connecting the water molecules to the carbonyl groups. In the construction of the hydrogen bonding pattern of the perfluorocarboxylic acid monohydrates, the oxalic acid dihydrate plays the key role. The striking similarity between the infrared spectra of the oxalic acid dihydrates and the perfluorocarboxylic acid monohydrates in the regions characteristic of water and OH⋯O vibration suggests that the structure of the hydrated carboxyl groups is the same in both crystals. These regions are characterized by the sharp doublet at 3539 cm^-1 and 3464 cm^-1, which is due to the H₂O ν₁ and ν₃ stretching vibrations, respectively, and the broad absorption between 3000 cm^-1 and 1500 cm^-1 with the intense band at 1970 cm^-1, both associated with the vibration of the OH⋯O group. The later peak consists of two band components at near 1980 cm^-1 and 2020 cm^-1. These band components show different behaviour when the temperature, polarization or deuteration is changed. In general, the infrared spectra of long-chain perfluorocarboxylic acids represent the system with very short hydrogen bonds connecting the water molecules to the carboxylates. This hydrogen bond pattern should be very similar to that found in the crystals of α-oxalic acid dihydrate.

Direct observation of coherent femtosecond solvent reorganization coupled to intramolecular electron transfer

It is well known that the solvent plays a critical role in ultrafast electron-transfer reactions. However, solvent reorganization occurs on multiple length scales, and selectively measuring short-range solute-solvent interactions at the atomic level with femtosecond time resolution remains a challenge. Here we report femtosecond X-ray scattering and emission measurements following photoinduced charge-transfer excitation in a mixed-valence bimetallic (FeⁱⁱRuⁱⁱⁱ) complex in water, and their interpretation using non-equilibrium molecular dynamics simulations. Combined experimental and computational analysis reveals that the charge-transfer excited state has a lifetime of 62 fs and that coherent translational motions of the first solvation shell are coupled to the back electron transfer. Our molecular dynamics simulations identify that the observed coherent translational motions arise from hydrogen bonding changes between the solute and nearby water molecules upon photoexcitation, and have an amplitude of tenths of ångströms, 120-200 cm^-1 frequency and ~100 fs relaxation time. This study provides an atomistic view of coherent solvent reorganization mediating ultrafast intramolecular electron transfer.

Temperature-Dependent Low-Frequency Vibrations of Thiamine Crystal Containing Hydrated Ions

Low-frequency vibrations of crystalline molecules are very sensitive to the local environment in which the molecules, for example, hydrated ions captured in crystals, find themselves. We present low-temperature X-ray crystallographic measurements on the harvested thiamine crystal containing hydrated ions and its temperature-dependent terahertz spectra and synchrotron infrared microspectra. It is found from the X-ray structure that the hydrated ions and hydration water are in a similar environment to liquid, although those are captured in crystals. The vibrationally resolved THz spectra of two states in the present organic crystals containing hydrated ions are well explained by the difference in the hydrogen-bonded pattern. Peak assignments were performed based on highly accurate first-principles calculations incorporating relativistic effects and dispersion corrections. The temperature dependences are observed for the vibrations around the chloride ions and hydration water due to the loose binding of chloride ions, the bond elongation with increasing temperature, and the cleavage of weak hydrogen bonds.

The effect of hydrogen bonds on the ultrafast relaxation dynamics of a BODIPY dimer

The influence of hydrogen bonds (H-bonds) in the structure, dynamics, and functionality of biological and artificial complex systems is the subject of intense investigation. In this broad context, particular attention has recently been focused on the ultrafast H-bond dependent dynamical properties in the electronic excited state because of their potentially dramatic consequences on the mechanism, dynamics, and efficiency of photochemical reactions and photophysical processes of crucial importance for life and technology. Excited-state H-bond dynamics generally occur on ultrafast time scales of hundreds of femtoseconds or less, making the characterization of associated mechanisms particularly challenging with conventional time-resolved techniques. Here, 2D electronic spectroscopy is exploited to shed light on this still largely unexplored dynamic mechanism. An H-bonded molecular dimer prepared by self-assembly of two boron-dipyrromethene dyes has been specifically designed and synthesized for this aim. The obtained results confirm that upon formation of H-bonds and the dimer, a new ultrafast relaxation channel is activated in the ultrafast dynamics, mediated by the vibrational motions of the hydrogen donor and acceptor groups. This relaxation channel also involves, beyond intra-molecular relaxations, an inter-molecular transfer process. This is particularly significant considering the long distance between the centers of mass of the two molecules. These findings suggest that the design of H-bonded structures is a particularly powerful tool to drive the ultrafast dynamics in complex materials.

Molecular Mechanism of Acceleration and Retardation of Collective Orientation Relaxation of Water Molecules in Aqueous Solutions

The collective orientation relaxation (COR) of water molecules in aqueous solutions is faster or slower with an increase in the concentration of the solutions than that in pure water; for example, acceleration (deceleration) of the COR is observed in a solution of sodium chloride (tetramethylammonium chloride) with increasing concentration. However, the molecular mechanism of the solution and concentration dependence of the relaxation time of the COR has not yet been clarified. We theoretically investigate the concentration dependence of the COR of water molecules in solutions of tetramethylammonium chloride (TMACl), guanidinium chloride (GdmCl), and sodium chloride (NaCl). Based on the Mori-Zwanzig equation, we identify two opposing factors that determine the COR of water molecules in any aqueous solution: the correlation of dipole moments and the single-molecule orientation relaxation. We reveal the molecular mechanism of the concentration dependence of the relaxation time of the COR in the TMACl, GdmCl, and NaCl solutions in terms of these two factors.

A unified quantum model susceptible to elucidate the dissimilarity of IR spectral density of dicarboxylic acid crystals: Phthalic and terephthalic acid crystals cases

Over the last decades, several approaches have been developed for elucidating the infrared spectral density of dicarboxylic acid crystals, which has been served as prototype for determining hydrogen bonds dynamics. These approaches differ in how accurately the simulated spectra can superimpose the experimental ones. In this study, we present a superdimer quantum approach susceptible to elucidate the infrared spectral properties of some particular dicarboxylic acid crystals using a newly proposed algorithm, which favors the rule of Davydov coupling in the generation of the spectra. The approach, which is herein effectively applied to terephthalic and phthalic acid dimer crystals, ascribes the non-conventional IR spectral properties of these particular acid crystals to the existence of superdimer structure in their lattices. In this superdimer structure, a strong vibronic coupling mechanism, namely Davydov coupling, takes place between the proton stretching vibrations in the (COOH)₂ cycles. This strong coupling exciton, generated by the resonance arising in the two coupled (COOH)₂ cycles of the aromatic rings of the superdimer, in conjunction with the strong anharmonic coupling between the fast and slow modes of each hydrogen bonds provide a strong support basis for a common explanation of the physical properties of these two different crystalline systems. The numerical simulations, involving the implications of the superdimer model, are systematically correlated with the experimental spectra. A decent agreement between the evaluated spectra and the experimental bandshapes of terephthalic and phthalic dicarboxylic acid crystals was obtained using a set of physically sound parameters as inputs in the theoretical formulation. The superdimer quantum approach thereby underscore the potential of the dynamical cooperative interactions between "Davydov coupling" and "strong anharmonic coupling" mechanisms in the generation of the spectral features of terephthalic and phthalic dicarboxylic acid crystals, suggesting that the congregated effects of these two mechanisms can be considered as the most reliable source of the non-conventional IR spectral properties observed. It is therefore expected that this novel algorithm reduces the discrepancies between the simulated spectra compared to the experimental one and simplify the computation of spectra in more complex hydrogen bonded systems.

Spatially dependent H-bond dynamics at interfaces of water/biomimetic self-assembled lattice materials

Understanding hydrogen-bond interactions in self-assembled lattice materials is crucial for preparing such materials, but the role of hydrogen bonds (H bonds) remains unclear. To gain insight into H-bond interactions at the materials' intrinsic spatial scale, we investigated ultrafast H-bond dynamics between water and biomimetic self-assembled lattice materials (composed of sodium dodecyl sulfate and β-cyclodextrin) in a spatially resolved manner. To accomplish this, we developed an infrared pump, vibrational sum-frequency generation (VSFG) probe hyperspectral microscope. With this hyperspectral imaging method, we were able to observe that the primary and secondary OH groups of β-cyclodextrin exhibit markedly different dynamics, suggesting distinct H-bond environments, despite being separated by only a few angstroms. We also observed another ultrafast dynamic reflecting a weakening and restoring of H bonds between bound water and the secondary OH of β-cyclodextrin, which exhibited spatial uniformity within self-assembled domains, but heterogeneity between domains. The restoration dynamics further suggest heterogeneous hydration among the self-assembly domains. The ultrafast nature and meso- and microscopic ordering of H-bond dynamics could contribute to the flexibility and crystallinity of the material--two critically important factors for crystalline lattice self-assemblies--shedding light on engineering intermolecular interactions for self-assembled lattice materials.

Icelike Vibrational Properties of Strong Hydrogen Bonds in Hydrated Lithium Nitrate

The hydrogen bond network accounts for many of the extraordinary physical properties of liquid water and ice. Its vibrational dynamics are quite complex in their entirety but can be accessed in detail by investigating small groups of only a few water molecules. Here, aqueous salt hydrates turned out to be an exceptional model system for water molecules arranged in well-defined geometrical structures that can be accessed by means of femtosecond spectroscopy of the OH stretching vibration. In this study, we find striking resemblance between the vibrational properties of three water molecules connected via strong hydrogen bonds in the trihydrate of LiNO₃ and those of ordinary ice I_h. As in ice, the vibrations of the hydrate water molecules show ultrafast excited state dynamics that are strongly accelerated when proceeding from deuterated to neat H₂O samples. The latter is analyzed by means of an additional relaxation channel that is due to Fermi resonance between the OH stretching vibration and the bend overtone accompanied by delocalization of the vibration over neighboring water molecules in the H₂O species. Moreover, in the hydrate and ice samples severe spectral broadening is examined when comparing fundamental and excited state absorption bands. Here, proton delocalization along the strong hydrogen bonds is given as a possible underlying mechanism.

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.

The nano-structural inhomogeneity of dynamic hydrogen bond network of TIP4P/2005 water

A method for studying the time dependence of the short-range molecular order of water has been proposed. In the present study, water is considered as a dynamic network between molecules at distances not exceeding 3.2 Å. The instantaneous configurations obtained with the molecular dynamics method have been sequentially analyzed. The mutual orientation of each molecule with its neighboring molecules has been studied and the interaction energy of each pair of neighbor molecules has been calculated. The majority of mutual orientation angles between molecules lie in the interval [0°; 20°]. More than 85% of the molecular pairs in each instantaneous configuration form H-bonds and the H-bond network includes all water molecules in the temperature range 233-293 K. The number of H-bonds fluctuates near the mean value and increases with decreasing temperature, and the energy of the vast majority of such bonds is much higher than the thermal energy. The interaction energy of 80% of the H-bonding molecular pairs lies in the interval [-7; -4] kcal/mol. The interaction energy of pairs that do not satisfy the H-bond angle criterion lies in the interval [-5; 4] kcal/mol; the number of such bonds does not exceed 15% and decreases with decreasing temperature. For the first time it has been found that in each instantaneous configuration the H-bond network contains built-in nanometric structural heterogeneities formed by shorter H-bonds. The fraction of molecules involved in the structural heterogeneities increases from 40% to 60% with a temperature decrease from 293 K to 233 K. Each heterogeneity has a finite lifetime and changeable structure, but they are constantly present during the entire simulation time.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study

We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.

Intramolecular hydrogen bonding protects the hydroxyl group from attack by fluctuating solvent forces

Ultrafast spectroscopy of molecular systems involving hydrogen- (H-) bonding has been at the forefront of fundamental chemical and physical research for several decades. Among the spectroscopic observables of the ultrafast dynamics is the pure dephasing of vibrationally excited molecules. Using third-order nonlinear vibrational spectroscopy, including polarization-selective transient grating measurements of vibrational lifetime and orientational diffusion as well as two-dimensional infrared spectroscopy, we determined different individual line shape components of hydroxyl stretching (νOH) excitations in a homologous series of chlorophenols and obtained the corresponding pure dephasing rates. The pure dephasing rates are correlated with vibrational anharmonicity of the νOH mode, which is tuned remotely from the hydroxyl site by changing the position of the chlorine substituents on the phenol ring. We found that in molecules where the hydroxyl group is in its free form, the pure dephasing rates scale linearly with the mode's anharmonicity such that assuming it is dominated by the third-order diagonal term, the ultrafast dynamics follow the prediction of the Kubo-Oxtoby theory. However, in the intramolecularly H-bonded ortho-chlorophenols, this trend is reversed, and the pure dephasing slows down by ∼50% for an increase in anharmonicity of only a few wavenumbers. Because the νOH mode's anharmonicity is known to reflect the H-bonding strength, our results suggest that intramolecular H-bonding can serve as a mechanism of protection from fluctuating forces exerted by the solvent. Such an effect can be relevant for ultrafast dynamics in biomolecules, where H-bonding plays a central role.

Distinctive behavior and two-dimensional vibrational dynamics of water molecules inside glycine solvation shell

We present a first principles molecular dynamics study of a deuterated aqueous solution of a single glycine moiety to explore the structure, dynamics, and two-dimensional infrared spectra of water molecules found in the solvation shell of glycine.

Vibrational echo spectroscopy of aqueous sodium bromide solutions from first principles simulations

A theoretical study of the time-dependent vibrational echo spectroscopy of sodium bromide solutions in deuterated water at two different concentrations of 0.5 and 5.0 M and at temperatures of 300 and 350 K is presented using the method of ab initio molecular dynamics simulations. The instantaneous fluctuations in frequencies of local OD stretch modes are calculated using time-series analysis of the simulated trajectories. The third-order polarization and intensities of three pulse photon-echo are calculated from ab initio simulations. The timescales of vibrational spectral diffusion are determined from the frequency time correlation functions (FTCF) and short-time slope of three pulse photon echo (S3PE) calculated within the second-order cumulant and Condon approximations. It is found that under ambient conditions, the rate of vibrational spectral diffusion becomes slower with increase in ionic concentration. Decay of S3PE calculated for different systems give timescales, which are in close agreement with those of FTCF and also with the results of experimental time-dependent vibrational spectroscopic experiments. © 2019 Wiley Periodicals, Inc.

Superheating and Homogeneous Melting Dynamics of Bulk Ice

Homogeneous melting of crystals is a complex multistep process involving the formation of transient states at temperatures considerably higher than the melting point. The nature and persistence of these metastable structures are intimately connected to the melting process, and a precise definition of the temporal boundaries of these phenomena is not yet available. We set up a specifically designed experiment to probe by transient infrared absorption spectroscopy the entire dynamics, ranging from tens of picoseconds to microseconds, of superheating and melting of an ice crystal. In spite of a large excess of energy provided, only about 30% of the micrometric crystal liquefies in the first 20-25 ns because of the long persistence of the superheated metastable phase that extends for more than 100 ns. This behavior is ascribed to the population of low-energy states that trap a large amount of energy, favoring the formation of a metastable, likely plastic, ice phase.

Effect of hydrogen bond on solvation dynamics of coumarin153 in cyclohexane-phenol solvent mixtures by time-resolved optical Kerr fluorescence

Time-resolved optical Kerr fluorescence system was used to investigate time-resolved red-shift of coumarin 153 in different solvent mixtures. The mixtures included four mole fractions of phenol-cyclohexane solvents (0, 0.013, 0.08, and 0.3), and anisole-cyclohexane solvents with the mole fraction of 0.3. The measured time-resolved fluorescence showed that, in the solvent mixtures containing phenol, the time-dependent frequency shift accelerated with the increase in the mole fraction of phenol-cyclohexane mixtures. However, the time-dependent red-shift in the fluorescence was not observed in the anisole-cyclohexane mixture, the solvent polarity could not influence the spectral Stokes shift compared with phenol. The results indicated that coumarin 153 formed an excited hydrogen bond with phenol, and the excited hydrogen bond was strengthened with an increase in the mole fraction of phenol. And, these processes also suggested that the increase in the phenol ratio improves a large number of hydrogen bond formed between phenol and carbonyl group of coumarin 153, the charge distribution will be faster towards lower the free energy of the system due to the stronger dipole moment. Therefore, the corresponding solvation response in phenol-cyclohexane mixtures with higher mole fractions decays very rapidly.

Vibrational Relaxation Dynamics of the Core and Outer Part of Proton-Hydration Clusters

We study the ultrafast relaxation dynamics of hydrated proton clusters in acetonitrile using femtosecond mid-infrared pump-probe spectroscopy. We observe a strong dependence of transient absorption dynamics on the frequency of excitation. When we excite the OH vibrations with frequencies ≤3100 cm^–1, we observe an ultrafast energy relaxation that leads to the heating of the local environment of the proton. This response is assigned to the OH vibrations of the water molecules in the core of the hydrated proton cluster. When we excite with frequencies ≥3200 cm^–1, we observe a relatively slow vibrational relaxation with a T₁ time constant ranging from 0.22 ± 0.04 ps at ν_ex = 3200 cm^–1 to 0.37 ± 0.02 ps at ν_ex = 3520 cm^–1. We assign this response to water molecules in the outer part of the hydrated proton cluster.