Anharmonic exciton dynamics and energy dissipation in liquid water from two-dimensional infrared spectroscopy

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.

Generation and Implementation of Continuum Infrared Pulses for Broadband Detection in 2D IR Spectroscopy

In recent years, new methods of generating continuum mid-infrared pulses through filamentation in gases have been developed for ultrafast time-resolved infrared vibrational spectroscopy. The generated infrared pulses can have thousands of wavenumbers of bandwidth, spanning the entire mid-IR region while retaining pulse length below 100 fs. This technology has had a significant impact on problems involving ultrafast structural dynamics in congested spectra with broad features, such as those found in aqueous solutions and molecules with strong intermolecular interactions. This study describes the recent advances in generating and characterizing these pulses and the practical aspects of implementing these sources for broadband detection in transient absorption and 2D IR spectroscopy.

Systematic analysis of electronic barrier heights and widths for concerted proton transfer in cyclic hydrogen bonded clusters: (HF)n, (HCl)n and (H2O)n where n = 3, 4, 5

The MP2 and CCSD(T) methods are paired with correlation consistent basis sets as large as aug-cc-pVQZ to optimize the structures of the cyclic minima for (HF)n, (HCl)n and (H2O)n where n = 3-5, as well as the corresponding transition states (TSs) for concerted proton transfer (CPT). MP2 and CCSD(T) harmonic vibrational frequencies confirm the nature of each minimum and TS. Both conventional and explicitly correlated CCSD(T) computations are employed to assess the electronic dissociation energies and barrier heights for CPT near the complete basis (CBS) limit for all 9 clusters. Results for (HF)n are consistent with prior studies identifying Cnh and Dnh point group symmetry for the minima and TSs, respectively. Our computations also confirm that CPT proceeds through Cs TS structures for the C1 minima of (H2O)3 and (H2O)5, whereas the process goes through a TS with D2d symmetry for the S4 global minimum of (H2O)4. This work corroborates earlier findings that the minima for (HCl)3, (HCl)4 and (HCl)5 have C3h, S4 and C1 point group symmetry, respectively, and that the Cnh structures are not minima for n = 4 and 5. Moreover, our computations show the TSs for CPT in (HCl)3, (HCl)4 and (HCl)5 have D3h, D2d, and C2 point group symmetry, respectively. At the CCSD(T) CBS limit, (HF)4 and (HF)5 have the smallest electronic barrier heights for CPT (≈15 kcal mol-1 for both), followed by the HF trimer (≈21 kcal mol-1). The barriers are appreciably higher for the other clusters (around 27 kcal mol-1 for (H2O)4 and (HCl)3; roughly 30 kcal mol-1 for (H2O)3, (H2O)5 and (HCl)4; up to 38 kcal mol-1 for (HCl)5). At the CBS limit, MP2 significantly underestimates the CCSD(T) barrier heights (e.g., by ca. 2, 4 and 7 kcal mol-1 for the pentamers of HF, H2O and HCl, respectively), whereas CCSD overestimates these barriers by roughly the same magnitude. Scaling the barrier heights and dissociation energies by the number of fragments in the cluster reveals strong linear relationships between the two quantities and with the magnitudes of the imaginary vibrational frequency for the TSs.

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.

Investigating Intramolecular H Atom Transfer Dynamics in β-Diketones with Ultrafast Infrared Spectroscopies and Theoretical Modeling

The vibrational signatures and ultrafast dynamics of the intramolecular H-bond in a series of β-diketones are investigated with 2D IR spectroscopy and computational modeling. The chosen β-diketones exhibit a range of H atom donor-acceptor distances and asymmetry along the H atom transfer coordinate that tunes the intramolecular H-bond strength. The species with the strongest H-bonds are calculated to have very soft H atom potentials, resulting in highly red-shifted OH stretch fundamental frequencies and dislocation of the H atom upon vibrational excitation. These soft potentials lead to significant coupling to the other normal mode coordinates and give rise to the very broad vibrational signatures observed experimentally. The 2D IR spectra in both the OH and OD stretch regions of the light and deuterated isotopologues reveal broadened and long-lived ground-state bleach signatures of the vibrationally hot molecules. Polarization-sensitive transient absorption measurements in the OH and OD stretch regions reveal notable isotopic differences in orientational dynamics. Orientational relaxation was measured to occur on ∼600 fs and ∼2 ps time scales for the light and deuterated isotopologues, respectively. The orientational dynamics are interpreted in terms of activated H/D atom transfer events driven by collective intramolecular structural rearrangements.

Discretized hierarchical equations of motion in mixed Liouville-Wigner space for two-dimensional vibrational spectroscopies of liquid water

A model of a bulk water system describing the vibrational motion of intramolecular and intermolecular modes is constructed, enabling analysis of its linear and nonlinear vibrational spectra as well as the energy transfer processes between the vibrational modes. The model is described as a system of four interacting anharmonic oscillators nonlinearly coupled to their respective heat baths. To perform a rigorous numerical investigation of the non-Markovian and nonperturbative quantum dissipative dynamics of the model, we derive discretized hierarchical equations of motion in mixed Liouville-Wigner space, with Lagrange-Hermite mesh discretization being employed in the Liouville space of the intramolecular modes and Lagrange-Hermite mesh discretization and Hermite discretization in the Wigner space of the intermolecular modes. One-dimensional infrared and Raman spectra and two-dimensional terahertz-infrared-visible and infrared-infrared-Raman spectra are computed as demonstrations of the quantum dissipative description provided by our model.

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.

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.

Exchange-Mediated Transport in Battery Electrolytes: Ultrafast or Ultraslow?

Understanding the mechanisms of charge transport in batteries is important for the rational design of new electrolyte formulations. Persistent questions about ion transport mechanisms in battery electrolytes are often framed in terms of vehicular diffusion by persistent ion-solvent complexes versus structural diffusion through the breaking and reformation of ion-solvent contacts, i.e., solvent exchange events. Ultrafast two-dimensional (2D) IR spectroscopy can probe exchange processes directly via the evolution of the cross-peaks on picosecond time scales. However, vibrational energy transfer in the absence of solvent exchange gives rise to the same spectral signatures, hiding the desired processes. We employ 2D IR on solvent resonances of a mixture of acetonitrile isotopologues to differentiate chemical exchange and energy-transfer dynamics in a comprehensive series of Li⁺, Mg²⁺, Zn²⁺, Ca²⁺, and Ba²⁺ bis(trifluoromethylsulfonyl)imide electrolytes from the dilute to the superconcentrated regime. No exchange phenomena occur within at least 100 ps, regardless of the ion identity, salt concentration, and presence of water. All of the observed spectral dynamics originate from the intermolecular energy transfer. These results place the lower experimental boundary on the ion-solvent residence times to several hundred picoseconds, much slower than previously suggested. With the help of MD simulations and conductivity measurements on the Li⁺ and Zn²⁺ systems, we discuss these results as a continuum of vehicular and structural modalities that vary with concentration and emphasize the importance of collective electrolyte motions to ion transport. These results hold broadly applicable to many battery-relevant ions and solvents.

Ultrafast Structure and Vibrational Dynamics of a Cyano-Containing Non-Fullerene Acceptor for Organic Solar Cells Revealed by Two-Dimensional Infrared Spectroscopy

Non-fullerene molecules, such as ITIC (3,9-bis(2-methylene-(3-(1,1-dicyanomethylene) indanone)-5,5,11,11-tetrakis(4-hexylphenyl)dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']-dithiophene), are among the most promising non-fullerene acceptors for organic solar cells (OSCs). In this work, using the cyano stretching mode as a vibrational marker, the structural and vibrational energy dynamics of ITIC were examined on an ultrafast time scale with two-dimensional infrared spectroscopy. Two IR-active modes studied here mainly correspond to two anti-symmetric combinations of symmetric and asymmetric stretching vibrations of two C≡N modes originating from two -C(CN)₂ chromophores that are located across the ITIC system, which were found to have significantly larger off-diagonal anharmonicity than their corresponding diagonal anharmonicities. This indicates strong anharmonic vibrational coupling between the two modes, which is supported by ab initio anharmonic frequency computations. Transient IR results indicate picosecond intramolecular vibrational energy transfer between the two C≡N modes upon excitation. The structural basis for these vibrational and energetic features is the conjugating molecular frame that is composed of a network of single/double bonds connecting the two -C(CN)₂ chromophores and may enable efficient vibration delocalization, in addition to its well-known electron delocalization capability. The importance of the results for the OSC applications is discussed.

Machine Learning Approach for Describing Water OH Stretch Vibrations

A machine learning approach employing neural networks is developed to calculate the vibrational frequency shifts and transition dipole moments of the symmetric and antisymmetric OH stretch vibrations of a water molecule surrounded by water molecules. We employed the atom-centered symmetry functions (ACSFs), polynomial functions, and Gaussian-type orbital-based density vectors as descriptor functions and compared their performances in predicting vibrational frequency shifts using the trained neural networks. The ACSFs perform best in modeling the frequency shifts of the OH stretch vibration of water among the types of descriptor functions considered in this paper. However, the differences in performance among these three descriptors are not significant. We also tried a feature selection method called CUR matrix decomposition to assess the importance and leverage of the individual functions in the set of selected descriptor functions. We found that a significant number of those functions included in the set of descriptor functions give redundant information in describing the configuration of the water system. We here show that the predicted vibrational frequency shifts by trained neural networks successfully describe the solvent-solute interaction-induced fluctuations of OH stretch frequencies.

Structural Characterization of Protonated Water Clusters Confined in HZSM-5 Zeolites

A molecular description of the structure and behavior of water confined in aluminosilicate zeolite pores is a crucial component for understanding zeolite acid chemistry under hydrous conditions. In this study, we use a combination of ultrafast two-dimensional infrared (2D IR) spectroscopy and ab initio molecular dynamics (AIMD) to study H₂O confined in the pores of highly hydrated zeolite HZSM-5 (∼13 and ∼6 equivalents of H₂O per Al atom). The 2D IR spectrum reveals correlations between the vibrations of both terminal and H-bonded O-H groups and the continuum absorption of the excess proton. These data are used to characterize the hydrogen-bonding network within the cluster by quantifying single-, double-, and non-hydrogen-bond donor water molecules. These results are found to be in good agreement with the statistics calculated from an AIMD simulation of an H⁺(H₂O)₈ cluster in HZSM-5. Furthermore, IR spectral assignments to local O-H environments are validated with DFT calculations on clusters drawn from AIMD simulations. The simulations reveal that the excess charge is detached from the zeolite and resides near the more highly coordinated water molecules in the cluster. When they are taken together, these results unambiguously assign the complex IR spectrum of highly hydrated HZSM-5, providing quantitative information on the molecular environments and hydrogen-bonding topology of protonated water clusters under extreme confinement.

Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport

The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and provide an outlook for prioritizing strategic research directions.

Computing the frequency fluctuation dynamics of highly coupled vibrational transitions using neural networks

The description of frequency fluctuations for highly coupled vibrational transitions has been a challenging problem in physical chemistry. In particular, the complexity of their vibrational Hamiltonian does not allow us to directly derive the time evolution of vibrational frequencies for these systems. In this paper, we present a new approach to this problem by exploiting the artificial neural network to describe the vibrational frequencies without relying on the deconstruction of the vibrational Hamiltonian. To this end, we first explored the use of the methodology to predict the frequency fluctuations of the amide I mode of N-methylacetamide in water. The results show good performance compared with the previous experimental and theoretical results. In the second part, the neural network approach is used to investigate the frequency fluctuations of the highly coupled carbonyl stretch modes for the organic carbonates in the solvation shell of the lithium ion. In this case, the frequency fluctuation predicted by the neural networks shows a good agreement with the experimental results, which suggests that this model can be used to describe the dynamics of the frequency in highly coupled transitions.

Vibrational couplings and energy transfer pathways of water's bending mode

Coupling between vibrational modes is essential for energy transfer and dissipation in condensed matter. For water, different O-H stretch modes are known to be very strongly coupled both within and between water molecules, leading to ultrafast dissipation and delocalization of vibrational energy. In contrast, the information on the vibrational coupling of the H-O-H bending mode of water is lacking, even though the bending mode is an essential intermediate for the energy relaxation pathway from the stretch mode to the heat bath. By combining static and femtosecond infrared, Raman, and hyper-Raman spectroscopies for isotopically diluted water with ab initio molecular dynamics simulations, we find the vibrational coupling of the bending mode differs significantly from the stretch mode: the intramode intermolecular coupling of the bending mode is very weak, in stark contrast to the stretch mode. Our results elucidate the vibrational energy transfer pathways of water. Specifically, the librational motion is essential for the vibrational energy relaxation and orientational dynamics of H-O-H bending mode.

Vibrational energy transfer in water involves intermolecular coupling of O-H stretching modes, but much less is known about the role of the bending modes. Here the authors, combining static and femtosecond infrared, Raman, and hyper-Raman spectroscopy and ab initio molecular dynamics simulations, provide insight into the energy dynamics of the bend vibrations.

The Bending Mode of Water: A Powerful Probe for Hydrogen Bond Structure of Aqueous Systems

Insights into the microscopic structure and dynamics of the water's hydrogen-bonded network are crucial to understand the role of water in biology, atmospheric and geochemical processes, and chemical reactions in aqueous systems. Vibrational spectroscopy of water has provided many such insights, in particular using the O-H stretch mode. In this Perspective, we summarize our recent studies that have revealed that the H-O-H bending mode can be an equally powerful reporter for the microscopic structure of water and provides more direct access to the hydrogen-bonded network than the conventionally studied O-H stretch mode. We discuss the fundamental vibrational properties of the water bending mode, such as the intermolecular vibrational coupling, and its effects on the spectral lineshapes and vibrational dynamics. Several examples of static and ultrafast bending mode spectroscopy illustrate how the water bending mode provides an excellent window on the microscopic structure of both bulk and interfacial water.

Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations

The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H₅O₂ ⁺ and Eigen H₉O₄ ⁺ ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure-spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm^-1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. Finally, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.

Vibrational coupling to hydration shell - Mechanism to performance enhancement of qualitative analysis in NIR spectroscopy of carbohydrates in aqueous environment

Vibrational coupling between carbohydrates and the hydration shell is unveiled as the underlying mechanism that improves wavenumber-selectively the carbohydrate discrimination performance by near-infrared (NIR) spectroscopy. The investigation is based on measurement of six carbohydrates (fructose, glucose, mannose, ribose, xylose and sorbitol) in aqueous solution in different concentration levels (5 mg/L, ~0.03 mmol/dm^-3 and 20 mg/L, ~0.1 mmol/dm^-3). The results of multivariate classification are interpreted by quantum mechanical NIR spectra simulations. The simulation unveils that the phenomenon is vibration-selective and thus wavenumber-selective, and leads to an enhancement of the qualitative information contained in the specific spectral regions. The location of these regions and the related performance correspond fully to the appearance and magnitude of the unveiled cooperative vibration effect.

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.

Kinetic Isotope Effects and Hydrogen Tunnelling in PCET Oxidations of Ascorbate: New Insights into Aqueous Chemistry?

Recent experimental studies of kinetic isotope effects (KIE-s) and hydrogen tunnelling comprising three proton-coupled electron transfer (PCET) oxidations of ascorbate monoanion, (a) in aqueous reaction solutions, (b) in the mixed water-organic cosolvent systems, (c) in aqueous solutions of various salts and (d) in fairly diluted aqueous solutions of the various partial hydrophobes are reviewed. A number of new insights into the wealth of the kinetic isotope phenomena in the PCET reactions have been obtained. The modulation of KIE-s and hydrogen tunnelling observed when partially hydrophobic solutes are added into water reaction solution, in the case of fairly diluted solutions is revealed as the strong linear correlation of the isotopic ratios of the Arrhenius prefactors Ah/Ad and the isotopic differences in activation energies ΔE_a (D,H). The observation has been proposed to be a signature of the involvement of the collective intermolecular excitonic vibrational dynamics of water in activation processes and aqueous chemistry.

Hydrogen Tunnelling as a Probe of the Involvement of Water Vibrational Dynamics in Aqueous Chemistry?

Our study of tunnelling in proton-coupled electron transfer (PCET) oxidation of ascorbate with hexacyanoferrate(III) follows the insights obtained from ultrafast 2D IR spectroscopy and theoretical studies of the vibrational water dynamics that led to the proposal of the involvement of collective intermolecular excitonic vibrational water dynamics in aqueous chemistry. To test the proposal, the hydrogen tunnelling modulation observed in the PCET reaction studied in the presence of low concentrations of various partial hydrophobic solutes in the water reaction system has been analyzed in terms of the proposed involvement of the collective intermolecular vibrational water dynamics in activation process in the case. The strongly linear correlation between common tunnelling signatures, isotopic values of Arrhenius prefactor ratios ln AH/AD and isotopic differences in activation enthalpies ΔΔH‡ (H,D) observed in the process in fairly diluted water solutions containing various partial hydrophobic solutes (such as dioxane, acetonitrile, ethanol, and quaternary ammonium ions) points to the common physical origin of the phenomenon in all the cases. It is suggested that the phenomenon can be rooted in an interplay of delocalized collective intermolecular vibrational dynamics of water correlated with vibrations of the coupled transition configuration, where the donor-acceptor oscillations, the motions being to some degree along the reaction coordinate, lead to modulation of hydrogen tunnelling in the reaction.

Capturing intrinsic site-dependent spectral signatures and lifetimes of isolated OH oscillators in extended water networks

The extremely broad infrared spectrum of water in the OH stretching region is a manifestation of how profoundly a water molecule is distorted when embedded in its extended hydrogen-bonding network. Many effects contribute to this breadth in solution at room temperature, which raises the question as to what the spectrum of a single OH oscillator would be in the absence of thermal fluctuations and coupling to nearby OH groups. We report the intrinsic spectral responses of isolated OH oscillators embedded in two cold (~20 K), hydrogen-bonded water cages adopted by the Cs⁺·(HDO)(D₂O)₁₉ and D₃O⁺·(HDO)(D₂O)₁₉ clusters. Most OH oscillators yield single, isolated features that occur with linewidths that increase approximately linearly with their redshifts. Oscillators near 3,400 cm^-1, however, occur with a second feature, which indicates that OH stretch excitation of these molecules drives low-frequency, phonon-type motions of the cage. The excited state lifetimes inferred from the broadening are considered in the context of fluctuations in the local electric fields that are available even at low temperature.

Entropic barriers in the kinetics of aqueous proton transfer

Aqueous proton transport is uniquely rapid among aqueous processes, mediated by fluctuating hydrogen bond reorganization in liquid water. In a process known as Grotthuss diffusion, the excess charge diffuses primarily by sequential proton transfers between water molecules rather than standard Brownian motion, which explains the anomalously high electrical conductivity of acidic solutions. Employing ultrafast IR spectroscopy, we use the orientational anisotropy decay of the bending vibrations of the hydrated proton complex to study the picosecond aqueous proton transfer kinetics as a function of temperature, concentration, and counterion. We find that the orientational anisotropy decay exhibits Arrhenius behavior, with an apparent activation energy of 2.4 kcal/mol in 1M and 2M HCl. Interestingly, acidic solutions at high concentration with longer proton transfer time scales display corresponding decreases in activation energy. We interpret this counterintuitive trend by considering the entropic and enthalpic contributions to the activation free energy for proton transfer. Halide counteranions at high concentrations impose entropic barriers to proton transfer in the form of constraints on the solution's collective H-bond fluctuations and obstruction of potential proton transfer pathways. The corresponding proton transfer barrier decreases due to weaker water-halide H-bonds in close proximity to the excess proton, but the entropic effects dominate and result in a net reduction in the proton transfer rate. We estimate the activation free energy for proton transfer as ∼1.0 kcal/mol at 280 K.

Deconstructing water's diffuse OH stretching vibrational spectrum with cold clusters

The diffuse vibrational envelope displayed by water precludes direct observation of how different hydrogen-bond topologies dictate the spectral response of individual hydroxy group (OH) oscillators. Using cold, isotopically labeled cluster ions, we report the spectral signatures of a single, intact water (H₂O) molecule embedded at various sites in the clathrate-like cage structure adopted by the Cs⁺·(D₂O)₂₀ ion. These patterns reveal the site-dependent correlation between the frequencies of the two OH groups on the same water molecule and establish that the bound OH companion of the free OH group exclusively accounts for bands in the lower-energy region of the spectrum. The observed multiplet structures reveal the homogeneous linewidths of the fundamentals and quantify the anharmonic contributions arising from coupling to both the intramolecular bending and intermolecular soft modes.

Intensified C≡C Stretching Vibrator and Its Potential Role in Monitoring Ultrafast Energy Transfer in 2D Carbon Material by Nonlinear Vibrational Spectroscopy

In this work, an intensity-enhanced C≡C stretching infrared (IR) absorption is observed in hexakis[(trimethylsilyl)ethynyl]benzene (HTEB), whose IR transition dipole magnitude becomes comparable to that of a typical C═O stretch, and the enhancement is believed to be due to a joint effect of π-π conjugation and hyperconjugation associated with a terminal trimethylsilyl group. Using dynamical time-dependent two-dimensional infrared (2D IR) spectroscopy, a picosecond intramolecular energy redistribution process is observed between two nondegenerate C≡C stretching modes, whose symmetry breaking is attributed to a noncovalent halogen-bonding interaction between HTEB and solvent CH₂Cl₂. The rigid structure of HTEB and limited structural dynamics are also inferred from the insignificant initial spectral diffusion value extracted from the 2D IR spectra. This work provides the first nonlinear infrared investigation of the conventionally weak C≡C stretch. The methods outlined are particularly important for detailed understanding of the structure-related processes such as vibrational energy transfer in novel C≡C species containing materials such as graphdiyne.

Ion Hydration and Ion Pairing as Probed by THz Spectroscopy

Ion hydration is of pivotal importance for many fundamental processes. Various spectroscopic methods are used to study the retardation of the hydration bond dynamics in the vicinity of anions and cations. Here we introduce THz-FTIR spectroscopy as a powerful method to answer the open questions. We show through dissection of THz spectra that we can pinpoint characteristic absorption features that can be attributed to the rattling modes of strongly hydrating ions within their hydration cages as well as vibrationally induced charge fluctuations in the case of weakly hydrating ions. Further analysis yields information on anion-cation cooperativity, the size of the dynamic hydration shell, as well as the lifetimes of these collective ion-hydration water modes and their connecting thermal bath states. Our study provides evidence for a non-additive behavior, thus questioning the simplified Hofmeister model. THz spectroscopy enables ion pairing to be observed and quantified at a high salt concentration.

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.

IR spectral assignments for the hydrated excess proton in liquid water

The local environmental sensitivity of infrared (IR) spectroscopy to a hydrogen-bonding structure makes it a powerful tool for investigating the structure and dynamics of excess protons in water. Although of significant interest, the line broadening that results from the ultrafast evolution of different solvated proton-water structures makes the assignment of liquid-phase IR spectra a challenging task. In this work, we apply a normal mode analysis using density functional theory of thousands of proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multistate empirical valence bond proton model (MS-EVB 3.2). These calculations are used to obtain a vibrational density of states and IR spectral density, which are decomposed on the basis of solvated proton structure and the frequency dependent mode character. Decompositions are presented on the basis of the proton sharing parameter δ, often used to distinguish Eigen and Zundel species, the stretch and bend character of the modes, the mode delocalization, and the vibrational mode symmetry. We find there is a wide distribution of vibrational frequencies spanning 1200-3000 cm^-1 for every local proton configuration, with the region 2000-2600 cm^-1 being mostly governed by the distorted Eigen-like configuration. We find a continuous red shift of the special-pair O⋯H⁺⋯O stretching frequency, and an increase in the flanking water bending intensity with decreasing δ. Also, we find that the flanking water stretch mode of the Zundel-like species is strongly mixed with the flanking water bend, and the special pair proton oscillation band is strongly coupled with the bend modes of the central H₅O2+moiety.

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.

Interpreting Quasi-Thermal Effects in Ultrafast Spectroscopy of Hydrogen-Bonded Systems

Vibrational excitation of molecules in the condensed phase relaxes through vibrational modes of decreasing energy to ultimately generate an equilibrium state in which the energy is distributed among low-frequency modes. In ultrafast vibrational spectroscopy, changes in the vibrational features of hydrogen-bonded NH and OH stretch modes are typically observed to persist long after these high-frequency vibrations have relaxed. Due to the resemblance to the spectral changes caused by heating the sample, these features are typically described as arising from a hot ground state. However, these spectral features appear on ultrafast time scales that are much too fast to result from a true thermal state, and significant differences between the thermal difference spectrum and the induced quasi-thermal changes in ultrafast spectroscopy are often observed. Here, we examine and directly compare the thermal and quasi-thermal responses of the hydrogen-bonded homodimer of 7-azaindole with temperature-dependent FTIR spectroscopy and ultrafast mid-IR continuum spectroscopy. We find that the thermal difference spectra contain contributions from both dissociation of the hydrogen bonds and from frequency shifts due to changes in the thermal population of low-frequency modes. The transient spectra in ultrafast vibrational spectroscopy are also found to contain two contributions: initial frequency shifts over 2.3 ± 0.11 ps associated with equilibration of the initial excitation, and frequency shifts associated with the excitation of several fingerprint modes, which decay over 21.8 ± 0.11 ps, giving rise to a quasi-thermal response caused by a distribution of fingerprint modes being excited within the sample ensemble. This resembles the thermal frequency shifts due to population changes of low-frequency modes, but not the overall thermal spectrum, which is dominated by features caused by dimer dissociation. These findings provide insight into the changes in the vibrational spectrum from different origins and are important for assigning, analyzing, and comparing features in thermal and ultrafast vibrational spectroscopy of hydrogen-bonded complexes.

Picosecond Proton Transfer Kinetics in Water Revealed with Ultrafast IR Spectroscopy

Aqueous proton transport involves the ultrafast interconversion of hydrated proton species that are closely linked to the hydrogen bond dynamics of water, which has been a long-standing challenge to experiments. In this study, we use ultrafast IR spectroscopy to investigate the distinct vibrational transition centered at 1750 cm^-1 in strong acid solutions, which arises from bending vibrations of the hydrated proton complex. Broadband ultrafast two-dimensional IR spectroscopy and transient absorption are used to measure vibrational relaxation, spectral diffusion, and orientational relaxation dynamics. The hydrated proton bend displays fast vibrational relaxation and spectral diffusion timescales of 200-300 fs; however, the transient absorption anisotropy decays on a remarkably long 2.5 ps timescale, which matches the timescale for hydrogen bond reorganization in liquid water. These observations are indications that the bending vibration of the aqueous proton complex is relatively localized, with an orientation that is insensitive to fast hydrogen bonding fluctuations and dependent on collective structural relaxation of the liquid to reorient. We conclude that the orientational relaxation is a result of proton transfer between configurations that are well described by a Zundel-like proton shared between two flanking water molecules.

Controlled modification of biomolecules by ultrashort laser pulses in polar liquids

Targeted chemical modification of peptides and proteins by laser pulses in a biologically relevant environment, i.e. aqueous solvent at room temperature, allows for accurate control of biological processes. However, the traditional laser methods of control of chemical reactions are applicable only to a small class of photosensitive biomolecules because of strong and ultrafast perturbations from biomolecule-solvent interactions. Here, we report excitation of harmonics of vibration modes of solvent molecules by femtosecond laser pulses to produce controlled chemical modifications of non-photosensitive peptides and proteins in polar liquids under room conditions. The principal modifications included lysine formylation and methionine sulfoxidation both of which occur with nearly 100% yield under atmospheric conditions. That modification occurred only if the laser irradiance exceeded certain threshold level. The threshold, type, and extent of the modifications were completely controlled by solvent composition, laser wavelength, and peak irradiance of ultrashort laser pulses. This approach is expected to assist in establishing rigorous control over a broad class of biological processes in cells and tissues at the molecular level.

Direct local solvent probing by transient infrared spectroscopy reveals the mechanism of hydrogen-bond induced nonradiative deactivation

The hydrogen-bond induced quenching of an excited chromophore is visualised by probing O–H vibrations of the interacting solvent molecules.

The fluorescence quenching of organic dyes via H-bonding interactions is a well-known phenomenon. However, the mechanism of this Hydrogen-Bond Induced Nonradiative Deactivation (HBIND) is not understood. Insight into this process is obtained by probing in the infrared the O–H stretching vibration of the solvent after electronic excitation of a dye with H-bond accepting cyano groups. The fluorescence lifetime of this dye was previously found to decrease from 1.5 ns to 110 ps when going from an aprotic solvent to the strongly protic hexafluoroisopropanol (HFP). Prompt strengthening of the H-bond with the dye was identified by the presence of a broad positive O–H band of HFP, located at lower frequency than the O–H band of the pure solvent. Further strengthening occurs within a few picoseconds before the excited H-bonded complex decays to the ground state in 110 ps. The latter process is accompanied by the dissipation of energy from the dye to the solvent and the rise of a characteristic hot solvent band in the transient spectrum. Polarization-resolved measurements evidence a collinear alignment of the nitrile and hydroxyl groups in the H-bonded complex, which persists during the whole excited-state lifetime. Measurements in other fluorinated alcohols and in chloroform/HFP mixtures reveal that the HBIND efficiency depends not only on the strength of the H-bond interactions between the dye and the solvent but also on the ability of the solvent to form an extended H-bond network. The HBIND process can be viewed as an enhanced internal conversion of an excited complex consisting of the dye molecule connected to a large H-bond network.

Dimension of discrete variable representation for mixed quantum/classical computation of three lowest vibrational states of OH stretching in liquid water

Three low-lying vibrational states of molecular systems are responsible for the signals of linear and third-order nonlinear vibrational spectroscopies. Theoretical studies based on mixed quantum/classical calculations provide a powerful way to analyze those experiments. A statistically meaningful result can be obtained from the calculations by solving the vibrational Schrödinger equation over many numbers of molecular configurations. The discrete variable representation (DVR) method is a useful technique to calculate vibrational eigenstates subject to an arbitrary anharmonic potential surface. Considering the large number of molecular configurations over which the DVR calculations are repeated, the calculations are desired to be optimized in balance between the cost and accuracy. We determine a dimension of the DVR method which appears to be optimum for the calculations of the three states of molecular vibrations with anharmonic strengths often found in realistic molecular systems. We apply the numerical technique to calculate the local OH stretching frequencies of liquid water, which are well known to be widely distributed due to the inhomogeneity in molecular configuration, and found that the frequencies of the 0-1 and 1-2 transitions are highly correlated. An empirical relation between the two frequencies is suggested and compared with the experimental data of nonlinear IR spectroscopies.