Pronounced non-Condon effects in the ultrafast infrared spectroscopy of water

Dynamics and Vibrational Spectroscopy of Alcohols in Ionic Liquids: Methanol and Ethanol

The structure, dynamics, and vibrational spectroscopy of dilute HOD, methanol, and ethanol in the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [emim][NTf₂], ionic liquid (IL) are investigated with molecular dynamics (MD) simulations. The structure of the ILs around the solutes is qualitatively similar, where the OD bond of the deuterated alcohols donates an interaction to an [NTf₂] anion and the [emim] cations interact with the oxygen atom of the OD group. The slowest time scale for the reorientational dynamics of the OD bond varied considerably for HOD, methanol, and ethanol (27, 71, and 87 ps, respectively). In contrast, the slowest time scales for spectral diffusion of the OD vibrational frequency were 11 ps for each of the three solutes, which indicates that the dynamics of the IL is relatively unchanged by the presence of the alcohols at dilute concentration. The theoretical results for the reorientational and spectral diffusion dynamics compare favorably with prior two-dimensional infrared (2D IR) spectroscopic measurements.

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.

Ab initio Modeling of the Vibrational Sum-Frequency Generation Spectrum of Interfacial Water

Understanding the structural and dynamical features of interfacial water is of greatest interest in physics, chemistry, biology, and materials science. Vibrational sum-frequency generation (SFG) spectroscopy, which is sensitive to the molecular orientation and dynamics on the surfaces or at the interfaces, allows one to study a wide variety of interfacial systems. The structural and dynamical features of interfacial water at the air/water interface have been extensively investigated by SFG spectroscopy. However, the interpretations of the spectroscopic features have been under intense debate. Here, we report a simulated SFG spectrum of the air/water interface based on ab initio molecular dynamics simulations, which covers the OH stretching, bending, and libration modes of interfacial water. Quantitative agreement between our present simulations and the most recent experimental studies ensures that ab initio simulations predict unbiased structural features and electrical properties of interfacial systems. By utilizing the kinetic energy spectral density (KESD) analysis to decompose the simulated spectra, the spectroscopic features can then be assigned to specific hydrogen-bonding configurations of interfacial water molecules.

Temperature dependence of the ultrafast vibrational echo spectroscopy of OD modes in liquid water from first principles simulations

The temperature dependence of the vibrational spectral diffusion of OD modes in liquid water is investigated through calculations of vibrational echo spectral observables from first principles molecular dynamics.

On the Hydrogen Bond Strength and Vibrational Spectroscopy of Liquid Water

In the present work, we introduce two new metrics i.e. hydrogen-bond strength and charge-transfer between the donor/acceptor water molecules as a measure of hydrogen-bond rearrangement dynamics. Further, we also employ a simple model based on energy flux through the donor-acceptor water pairs to quantify the extent of the local hydrogen-bond network reorganization. Most importantly, we report a linear relationship between the OH stretch frequency and the charge and energy transfer through donor-acceptor water pairs. We demonstrate that the vibrational frequency fluctuations, which are used to determine third-order non-linear spectroscopic observables like the short-time slope of three pulse photon echo, can be used as an analog of the fluctuations in the hydrogen-bond strength and charge-transfer. The timescales obtained from our hydrogen-bond strength correlation and charge-transfer correlation decay are in excellent agreement with the computed frequency-time correlation function, as well as with recent vibrational echo experiments.

Modeling Carbon Dioxide Vibrational Frequencies in Ionic Liquids: III. Dynamics and Spectroscopy

In recent years, interest in carbon capture and sequestration has led to numerous investigations of the ability of ionic liquids to act as recyclable CO₂-sorbent materials. Herein, we investigate the structure and dynamics of a model physisorbing ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate ([C₄C₁Im][PF₆]), from the perspective of CO₂ using two-dimensional (2D) IR spectroscopy and molecular dynamics simulations. A direct comparison of experimentally measured and calculated 2D IR line shapes confirms the validity of the simulations and spectroscopic calculations. Taken together, the simulations and experiments reveal new insights into the interactions of a CO₂ solute with the surrounding ionic liquid and how these interactions manifest in the 2D IR spectra. In particular, higher CO₂ asymmetric stretch vibrational frequencies are associated with softer, less populated solvent cages and lower frequencies are associated with tighter, more highly populated solvent cages. The CO₂ interacts most strongly with the anions, and these interactions persist for more than 1 ns. The second strongest interactions are with the imidazolium cation ring that last 100 ps, and the weakest interactions are with the cation butyl tail that persist for 10 ps. The principal contributors to spectral diffusion of the CO₂ asymmetric stretch vibrational frequency due to the dynamical evolution of the solvent are through Lennard-Jones interactions at short times and electrostatics at long times.

Hydrogen bonding and vibrational energy relaxation of interfacial water: A full DFT molecular dynamics simulation

The air-water interface has been a subject of extensive theoretical and experimental studies due to its ubiquity in nature and its importance as a model system for aqueous hydrophobic interfaces. We report on the structure and vibrational energy transfer dynamics of this interfacial water system studied with equilibrium and non-equilibrium molecular dynamics simulations employing a density functional theory -based description of the system and the kinetic energy spectral density analysis. The interfacial water molecules are found to make fewer and weaker hydrogen (H)-bonds on average compared to those in the bulk. We also find that (i) the H-bonded OH groups conjugate to the free OH exhibit rather low vibrational frequencies (3000-3500 cm^-1); (ii) the presence of a significant fraction (>10%) of free and randomly oriented water molecules at the interface ("labile water"), neither of whose OH groups are strong H-bond donors; (iii) the inertial rotation of free OH groups, especially from the labile water, contribute to the population decay of excited free OH groups with comparable rate and magnitude as intramolecular energy transfer between the OH groups. These results suggest that the labile water, which might not be easily detectable by the conventional vibrational sum frequency generation method, plays an important role in the surface water dynamics.

Spectral line shapes in linear absorption and two-dimensional spectroscopy with skewed frequency distributions

The effect of Gaussian dynamics on the line shapes in linear absorption and two-dimensional correlation spectroscopy is well understood as the second-order cumulant expansion provides exact spectra. Gaussian solvent dynamics can be well analyzed using slope line analysis of two-dimensional correlation spectra as a function of the waiting time between pump and probe fields. Non-Gaussian effects are not as well understood, even though these effects are common in nature. The interpretation of the spectra, thus far, relies on complex case to case analysis. We investigate spectra resulting from two physical mechanisms for non-Gaussian dynamics, one relying on the anharmonicity of the bath and the other on non-linear couplings between bath coordinates. These results are compared with outcomes from a simpler log-normal dynamics model. We find that the skewed spectral line shapes in all cases can be analyzed in terms of the log-normal model, with a minimal number of free parameters. The effect of log-normal dynamics on the spectral line shapes is analyzed in terms of frequency correlation functions, maxline slope analysis, and anti-diagonal linewidths. A triangular line shape is a telltale signature of the skewness induced by log-normal dynamics. We find that maxline slope analysis, as for Gaussian dynamics, is a good measure of the solvent dynamics for log-normal dynamics.

Nuclear quantum effects on the vibrational dynamics of liquid water

Based on quantum-mechanical path-integral molecular dynamics simulations, the impact of nuclear quantum effects on the vibrational and hydrogen bond dynamics in liquid water is investigated. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time-series analysis, while the time scales of the vibrational spectral diffusion are determined from frequency-time correlation functions, joint probability distributions, and the slope of three-pulse photon echo. We find that the inclusion of nuclear quantum effects leads not only to a redshift of the vibrational frequency distribution by around 130 cm^-1 but also to an acceleration of the vibrational dynamics by as much as 30%. In addition, quantum fluctuations also entail a significantly faster decay of correlation in the initial diffusive regime, which is in agreement with recent vibrational echo experiments.

Simple Quantum Dynamics with Thermalization

In this paper, we introduce two simple quantum dynamics methods. One is based on the popular surface-hopping method, and the other is based on rescaling of the propagation on the bath ground-state potential surface. The first method is special, as it avoids specific feedback from the simulated quantum system to the bath and can be applied for precalculated classical trajectories. It is based on the equipartition theorem to determine if hops between different potential energy surfaces are allowed. By comparing with the formally exact Hierarchical Equations Of Motion approach for four model systems we find that the method generally approximates the quantum dynamics toward thermal equilibrium very well. The second method is based on rescaling of the nonadiabatic coupling and also neglect the effect of the state of the quantum system on the bath. By the nature of the approximations, they cannot reproduce the effect of bath relaxation following excitation. However, the methods are both computationally more tractable than the conventional fewest switches surface hopping, and we foresee that the methods will be powerful for simulations of quantum dynamics in systems with complex bath dynamics, where the system-bath coupling is not too strong compared to the thermal energy.

Decomposition of the Experimental Raman and Infrared Spectra of Acidic Water into Proton, Special Pair, and Counterion Contributions

Textbooks describe excess protons in liquid water as hydronium (H₃O⁺) ions, although their true structure remains lively debated. To address this question, we have combined Raman and infrared (IR) multivariate curve resolution spectroscopy with ab initio molecular dynamics and anharmonic vibrational spectroscopic calculations. Our results are used to resolve, for the first time, the vibrational spectra of hydrated protons and counterions and reveal that there is little ion-pairing below 2 M. Moreover, we find that isolated excess protons are strongly IR active and nearly Raman inactive (with vibrational frequencies of ∼1500 ± 500 cm^-1), while flanking water OH vibrations are both IR and Raman active (with higher frequencies of ∼2500 ± 500 cm^-1). The emerging picture is consistent with Georg Zundel's seminal work, as well as recent ultrafast dynamics studies, leading to the conclusion that protons in liquid water are primarily hydrated by two flanking water molecules, with a broad range of proton hydrogen bond lengths and asymmetries.

Electrostatic frequency maps for amide-I mode of β-peptide: Comparison of molecular mechanics force field and DFT calculations

The spectroscopy of amide-I vibrations has been widely utilized for the understanding of dynamical structure of polypeptides. For the modeling of amide-I spectra, two frequency maps were built for β-peptide analogue (N-ethylpropionamide, NEPA) in a number of solvents within different schemes (molecular mechanics force field based, GM map; DFT calculation based, GD map), respectively. The electrostatic potentials on the amide unit that originated from solvents and peptide backbone were correlated to the amide-I frequency shift from gas phase to solution phase during map parameterization. GM map is easier to construct with negligible computational cost since the frequency calculations for the samples are purely based on force field, while GD map utilizes sophisticated DFT calculations on the representative solute-solvent clusters and brings insight into the electronic structures of solvated NEPA and its chemical environments. The results show that the maps' predicted amide-I frequencies present solvation environmental sensitivities and exhibit their specific characters with respect to the map protocols, and the obtained vibrational parameters are in satisfactory agreement with experimental amide-I spectra of NEPA in solution phase. Although different theoretical schemes based maps have their advantages and disadvantages, the present maps show their potentials in interpreting the amide-I spectra for β-peptides, respectively.

A "Universal" Spectroscopic Map for the OH Stretching Mode in Alcohols

Empirical maps are presented for the OH stretching vibrations in neat alcohols in which the relevant spectroscopic quantities are expressed in terms of the electric field exerted on the hydrogen atom by the surrounding liquid. It is found, by examination of the four lowest linear alcohols, methanol, ethanol, n-propanol, and n-butanol, that a single map can be used for alcohols with different alkyl groups. This "universal" map is in very good agreement with maps optimized for the individual alcohols but differs from those previously developed for water. This suggests that one map can be used for all alcohols, perhaps even those not examined in the present study. The universal map gives IR lineshapes in good agreement with measured spectra for isotopically dilute methanol and ethanol, while the two-dimensional IR photon echo spectra give results that differ from experiments. The role of non-Condon effects, reorientation dynamics, hydrogen bonding, and spectral diffusion is discussed.

Role of Solvent Water in the Temperature-Induced Self-Assembly of a Triblock Copolymer

Water-soluble triblock copolymers have received much attention in industrial applications and scientific fields. We here show that femtosecond mid-IR pump-probe spectroscopy is useful to study the role of water in the temperature-induced self-assembly of triblock copolymers. Our experimental results suggest two distinct subpopulations of water molecules: those that interact with other water molecules and those involved in the hydration of a triblock copolymer surface. We find that the vibrational dynamics of bulk-like water is not affected by either micellation or gelation of triblock copolymers. The increased population of water interacting with ether oxygen atoms of the copolymer during the unimer to micelle phase transition is important evidence for the entropic role of water in temperature-induced micelle formation at a low copolymer concentration. In contrast, at the critical gelation temperature and beyond, the population of surface-associated water molecules interacting with ether oxygen atoms decreases, which indicates important enthalpic control by water. The present study on the roles of water in the two different phase transitions of triblock copolymers sheds new light on the underlying mechanisms of temperature-induced self-aggregation behaviors of amphiphiles that are ubiquitous in nature.

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.

Vibrational Probes: From Small Molecule Solvatochromism Theory and Experiments to Applications in Complex Systems

The vibrational frequency of a chosen normal mode is one of the most accurately measurable spectroscopic properties of molecules in condensed phases. Accordingly, infrared absorption and Raman scattering spectroscopy have provided valuable information on both distributions and ensemble-average values of molecular vibrational frequencies, and these frequencies are now routinely used to investigate structure, conformation, and even absolute configuration of chemical and biological molecules of interest. Recent advancements in coherent time-domain nonlinear vibrational spectroscopy have allowed the study of heterogeneous distributions of local structures and thermally driven ultrafast fluctuations of vibrational frequencies. To fully utilize IR probe functional groups for quantitative bioassays, a variety of biological and chemical techniques have been developed to site-specifically introduce vibrational probe groups into proteins and nucleic acids. These IR-probe-labeled biomolecules and chemically reactive systems are subject to linear and nonlinear vibrational spectroscopic investigations and provide information on the local electric field, conformational changes, site-site protein contacts, and/or function-defining features of biomolecules. A rapidly expanding library of data from such experiments requires an interpretive method with atom-level chemical accuracy. However, despite prolonged efforts to develop an all-encompassing theory for describing vibrational solvatochromism and electrochromism as well as dynamic fluctuations of instantaneous vibrational frequencies, purely empirical and highly approximate theoretical models have often been used to interpret experimental results. They are, in many cases, based on the simple assumption that the vibrational frequency of an IR reporter is solely dictated by electric potential or field distribution around the vibrational chromophore. Such simplified description of vibrational solvatochromism generally referred to as vibrational Stark effect theory has been considered to be quite appealing and, even in some cases, e.g., carbonyl stretch modes in amide, ester, ketone, and carbonate compounds or proteins, it works quantitatively well, which makes it highly useful in determining the strength of local electric field around the IR chromophore. However, noting that the vibrational frequency shift results from changes of solute-solvent intermolecular interaction potential along its normal coordinate, Pauli exclusion repulsion, polarization, charge transfer, and dispersion interactions, in addition to the electrostatic interaction between distributed charges of both vibrational chromophore and solvent molecules, are to be properly included in the theoretical description of vibrational solvatochromism. Since the electrostatic and nonelectrostatic intermolecular interaction components have distinctively different distance and orientation dependences, they affect the solvatochromic vibrational properties in a completely different manner. Over the past few years, we have developed a systematic approach to simulating vibrational solvatochromic data based on the effective fragment potential approach, one of the most accurate and rigorous theories on intermolecular interactions. We have further elucidated the interplay of local electric field with the general vibrational solvatochromism of small IR probes in either solvents or complicated biological systems, with emphasis on contributions from non-Coulombic intermolecular interactions to vibrational frequency shifts and fluctuations. With its rigorous foundation and close relation to quantitative interpretation of experimental data, this and related theoretical approaches and experiments will be of use in studying and quantifying the structure and dynamics of biomolecules with unprecedented time and spatial resolution when combined with time-resolved vibrational spectroscopy and chemically sensitive vibrational imaging techniques.

Molecular modeling and assignment of IR spectra of the hydrated excess proton in isotopically dilute water

Infrared (IR) spectroscopy of the water O-H stretch has been widely used to probe both the local hydrogen-bonding structure and dynamics of aqueous systems. Although of significant interest, the IR spectroscopy of excess protons in water remains difficult to assign as a result of extensive and strong intermolecular interactions in hydrated proton complexes. As an alternate approach, we develop a mixed quantum-classical model for the vibrational spectroscopy of the excess proton in isotopically dilute water that draws on frozen proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multi-state empirical valence bond proton model (MS-EVB 3.2). A semi-empirical single oscillator spectroscopic map for the instantaneous transition frequency and transition dipole moment is constructed using potential energy surfaces for the O-H stretch coordinate of the excess proton using electronic structure calculations. Calculated spectra are compared with experimental spectra of dilute H⁺ in D₂O obtained from double-difference FTIR to demonstrate the validity of the map. The model is also used to decompose IR spectra into contributions from different aqueous proton configurations. We find that the O-H transition frequency continuously decreases as the oxygen-oxygen length for a special pair proton decreases, shifting from Eigen- to Zundel-like configurations. The same shift is accompanied by a shift of the flanking water stretches of the Zundel complex to higher frequency than the hydronium O-H vibrations.

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.

Probing the dynamics of N-methylacetamide in methanol via ab initio molecular dynamics

Two-dimensional infrared (2D IR) spectroscopy of amide 1 vibrational bands provides a valuable probe of proteins as well as molecules such as N-methylacetamide (NMA), which present peptide-like H-bonding possibilities to a solvent.

Observation of isolated ionic liquid cations and water molecules in an inert solvent

1-Octyl-3-methyl imidazolium halides ([OMIM]I and [OMIM]Cl) were loaded on top of CCl₄, and an in situ inclusion process was monitored from the CCl₄ phase as time elapses by infrared absorption spectroscopy. Absorption from the bands corresponding to the C(2)-H and C(4,5)-H stretch modes in the imidazolium cation was reduced significantly compared to the bulk IL spectra. This indicates that (1) the [OMIM] cation exists in CCl₄ as a monomer, dissociated from the anion and other cations, and (2) hydrogen bonding between the anion and the cation increases the dipole strength of the CH moieties in the imidazolium ring. [OMIM]I was found to transfer into the CCl₄ phase much faster than [OMIM]Cl, and this instigated us to compare the transfer of aqueous solutions of ionic liquids, 1-butyl-3-methyl imidazolium halides ([BMIM]I and [BMIM]Cl) into the CCl₄ matrix. Not only [BMIM]I but also water molecules transferred faster compared to those in the [BMIM]Cl aqueous solution. Water molecules in CCl₄ were shown to form clusters in [BMIM]I; presumably, I^- anions work as nucleation centers of water clusters.

Vibrational Relaxation of the Aqueous Proton in Acetonitrile: Ultrafast Cluster Cooling and Vibrational Predissociation

We study the ultrafast O-H stretch vibrational relaxation dynamics of protonated water clusters embedded in a matrix of deuterated acetonitrile, using polarization-resolved mid-IR femtosecond spectroscopy. The clusters are produced by mixing triflic (trifluoromethanesulfonic) acid and H2O in molar ratios of 1:1, 1:2, and 1:3, thus varying the degree of hydration of the proton. At all hydration levels the excited O-H stretch vibration of the hydrated proton shows an ultrafast vibrational relaxation with a time constant T1 < 100 fs, leading to an ultrafast local heating of the protonated water cluster. This excess thermal energy, initially highly localized to the region of the excited proton, first re-distributes over the aqueous cluster and then dissipates into the surrounding acetonitrile matrix. For clusters with a triflic acid to H2O ratio of 1:3 these processes occur with time constants of 320 ± 20 fs and 1.4 ± 0.1 ps, respectively. The cooling of the clusters reveals a long-living, underlying transient absorption change with high anisotropy. We argue that this feature stems from the vibrational predissociation of a small fraction of the proton hydration structures, directly following the ultrafast infrared excitation.

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

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

Theoretical study of vibrational energy transfer of free OH groups at the water-air interface

Recent experimental studies have shown that the vibrational dynamics of free OH groups at the water-air interface is significantly different from that in bulk water. In this work, by performing molecular dynamics simulations and mixed quantum/classical calculations, we investigate different vibrational energy transfer pathways of free OH groups at the water-air interface. The calculated intramolecular vibrational energy transfer rate constant and the free OH bond reorientation time scale agree well with the experiment. It is also found that, due to the small intermolecular vibrational couplings, the intermolecular vibrational energy transfer pathway that is very important in bulk water plays a much less significant role in the vibrational energy relaxation of the free OH groups at the water-air interface.

Hydrogen bond dynamics in bulk alcohols

Hydrogen-bonded liquids play a significant role in numerous chemical and biological phenomena. In the past decade, impressive developments in multidimensional vibrational spectroscopy and combined molecular dynamics-quantum mechanical simulation have established many intriguing features of hydrogen bond dynamics in one of the fundamental solvents in nature, water. The next class of a hydrogen-bonded liquid--alcohols--has attracted much less attention. This is surprising given such important differences between water and alcohols as the imbalance between the number of hydrogen bonds, each molecule can accept (two) and donate (one) and the very presence of the hydrophobic group in alcohols. Here, we use polarization-resolved pump-probe and 2D infrared spectroscopy supported by extensive theoretical modeling to investigate hydrogen bond dynamics in methanol, ethanol, and isopropanol employing the OH stretching mode as a reporter. The sub-ps dynamics in alcohols are similar to those in water as they are determined by similar librational and hydrogen-bond stretch motions. However, lower density of hydrogen bond acceptors and donors in alcohols leads to the appearance of slow diffusion-controlled hydrogen bond exchange dynamics, which are essentially absent in water. We anticipate that the findings herein would have a potential impact on fundamental chemistry and biology as many processes in nature involve the interplay of hydrophobic and hydrophilic groups.

On the interplay of the potential energy and dipole moment surfaces in controlling the infrared activity of liquid water

Infrared vibrational spectroscopy is a valuable tool for probing molecular structure and dynamics. However, obtaining an unambiguous molecular-level interpretation of the spectral features is made difficult, in part, due to the complex interplay of the dipole moment with the underlying vibrational structure. Here, we disentangle the contributions of the potential energy surface (PES) and dipole moment surface (DMS) to the infrared spectrum of liquid water by examining three classes of models, ranging in complexity from simple point charge models to accurate representations of the many-body interactions. By decoupling the PES from the DMS in the calculation of the infrared spectra, we demonstrate that the PES, by directly modulating the vibrational structure, primarily controls the width and position of the spectroscopic features. Due to the dependence of the molecular dipole moment on the hydration environment, many-body electrostatic effects result in a ∼100 cm(-1) redshift in the peak of the OH stretch band. Interestingly, while an accurate description of many-body collective motion is required to generate the correct (vibrational) structure of the liquid, the infrared intensity in the OH stretching region appears to be a measure of the local structure due to the dominance of the one-body and short-ranged two-body contributions to the total dipole moment.