Nonlinear response of vibrational excitons: simulating the two-dimensional infrared spectrum of liquid water

Theoretical Two Dimensional Infrared Spectroscopy of Aqueous Solutions of tert-Butyl Alcohol: Variation of the Dynamics of Spectral Diffusion along the Percolation Transition

Binary mixtures of water and tert-butyl alcohol (TBA) are known to exhibit the so-called percolation transition where small clusters of TBA molecules span into large aggregates beyond a threshold concentration of the alcohol. In the present study, we have investigated the linear and two-dimensional infrared spectral features of aqueous solutions of TBA for varying concentration of the alcohol along the percolation transition. The percolation transition is characterized through calculations of intermolecular radial distribution functions and average size of the largest cluster of TBA molecules. It is found that, with variation of alcohol concentration, the radial distribution functions of the central carbon atoms of TBA molecules show a nonmonotonic change in the height of the first peak and also the size of the largest cluster of TBA molecules show a jump in the increase of its size for TBA mole fraction between 0.04 and 0.06 corresponding to a transition from smaller clusters to larger spanning aggregates. However, it is found that the linear infrared spectrum of water does not exhibit any noticeable changes on variation of TBA concentration along the percolation transition. Subsequently, two-dimensional infrared (2DIR) spectra and vibrational frequency time correlation function of water are calculated for all the TBA-water solutions considered in this study. The spectral diffusion of water calculated from 2DIR is found to slow down with increase of the TBA concentration. The time scales of spectral diffusion of water, as characterized by the relaxation of frequency time correlation function, 2DIR metric of central line slope, and also the hydrogen bond time correlation functions, are found to exhibit a noticeable jump along the percolation transition. The hydrophilic group of TBA is found to retard the water dynamics more effectively than the hydrophobic groups. Also, the jump in the dynamical slowdown along the percolation transition is found to be more significant for water molecules at the hydrophilic sites.

Simulation of two-dimensional infrared Raman spectroscopy with application to proteins

Two-dimensional infrared Raman spectroscopy is a powerful technique for studying the structure and interaction in molecular and biological systems. Here, we present a new implementation of the simulation of the two-dimensional infrared Raman signals. The implementation builds on the numerical integration of the Schrödinger equation approach. It combines the prediction of dynamics from molecular dynamics with a map-based approach for obtaining Hamiltonian trajectories and response function calculations. The new implementation is tested on the amide-I region for two proteins, where one is dominated by α-helices and the other by β-sheets. We find that the predicted spectra agree well with experimental observations. We further find that the two-dimensional infrared Raman spectra at least of the studied proteins are much less sensitive to the laser polarization used compared to conventional two-dimensional infrared experiments. The present implementation and findings pave the way for future applications for the interpretation of two-dimensional infrared Raman spectra.

Computational spectroscopy of complex systems

Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.

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.

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 Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

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.

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.

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.

Vibrational dynamics of confined supercooled water

The quest for a possible liquid-liquid coexistence line in supercooled water below its homogeneous nucleation temperature is faced by confining water within a porous silica substrate (MCM-41). This system is investigated by synchrotron radiation infrared spectroscopy, exploring both the intramolecular and the intermolecular vibrational dynamics, in the temperature range from ambient down to ∼120 K, along several isobaric paths between 0.7 kbar and 3.0 kbar. Upon lowering the temperature, the OH-stretching band shows that the intramolecular vibrational dynamics continuously evolves from predominantly liquidlike to predominantly icelike. An abrupt change in the line shape of the intermolecular vibrational band between 220 K and 240 K, depending on the pressure, is the signature of nucleation of ice within the MCM-41 pores. These findings do not support the presence of two liquid phases and provide evidence for the coexistence of liquid water and ice in water confined in MCM-41.

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.

Ballistic and diffusive vibrational energy transport in molecules

Energy transport in molecules is essential for many areas of science and technology. Strong covalent bonds of a molecular backbone can facilitate the involvement of the molecule's high-frequency modes in energy transport, which, under certain conditions, makes the transport fast and efficient. We discuss such conditions and describe various transport regimes in molecules, including ballistic, diffusive, directed diffusion, and intermediate regime cases, in light of recently developed experimental and theoretical approaches.

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.

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.

Water at surfaces with tunable surface chemistries

Aqueous interfaces are ubiquitous in natural environments, spanning atmospheric, geological, oceanographic, and biological systems, as well as in technical applications, such as fuel cells and membrane filtration. Where liquid water terminates at a surface, an interfacial region is formed, which exhibits distinct properties from the bulk aqueous phase. The unique properties of water are governed by the hydrogen-bonded network. The chemical and physical properties of the surface dictate the boundary conditions of the bulk hydrogen-bonded network and thus the interfacial properties of the water and any molecules in that region. Understanding the properties of interfacial water requires systematically characterizing the structure and dynamics of interfacial water as a function of the surface chemistry. In this review, we focus on the use of experimental surface-specific spectroscopic methods to understand the properties of interfacial water as a function of surface chemistry. Investigations of the air-water interface, as well as efforts in tuning the properties of the air-water interface by adding solutes or surfactants, are briefly discussed. Buried aqueous interfaces can be accessed with careful selection of spectroscopic technique and sample configuration, further expanding the range of chemical environments that can be probed, including solid inorganic materials, polymers, and water immiscible liquids. Solid substrates can be finely tuned by functionalization with self-assembled monolayers, polymers, or biomolecules. These variables provide a platform for systematically tuning the chemical nature of the interface and examining the resulting water structure. Finally, time-resolved methods to probe the dynamics of interfacial water are briefly summarized before discussing the current status and future directions in studying the structure and dynamics of interfacial water.

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.

Coupling between intra- and intermolecular motions in liquid water revealed by two-dimensional terahertz-infrared-visible spectroscopy

The interaction between intramolecular and intermolecular degrees of freedom in liquid water underlies fundamental chemical and physical phenomena such as energy dissipation and proton transfer. Yet, it has been challenging to elucidate the coupling between these different types of modes. Here, we report on the direct observation and quantification of the coupling between intermolecular and intramolecular coordinates using two-dimensional, ultra-broadband, terahertz-infrared-visible (2D TIRV) spectroscopy and molecular dynamics calculations. Our study reveals strong coupling of the O-H stretch vibration, independent of the degree of delocalization of this high-frequency mode, to low-frequency intermolecular motions over a wide frequency range from 50 to 250 cm-1, corresponding to both the intermolecular hydrogen bond bending (≈ 60 cm-1) and stretching (≈ 180 cm-1) modes. Our results provide mechanistic insights into the coupling of the O-H stretch vibration to collective, delocalized intermolecular modes.

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.

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.

Two-dimensional infrared spectroscopy of neat ice Ih

The assignment of the distinct peaks observed in the OH stretch lineshape of ice Ih is controversial. Recent two-dimensional infrared spectroscopic measurements provided new data. The spectra are, however, challenging to interpret and here we provide simulations that help overcome experimental issues as thermal signals and finite pulse duration. We find good agreement with experiment and the difference between H2O and D2O ices is well accounted for. The overall dynamics is demonstrated to be faster than observed for the corresponding liquid water. We find that excitonic cross peaks exist between the dominant exciton peaks. This leads us to conclude that the cross peaks arise due to the formation of delocalized exciton states, which have essentially no directional correlation between their transition dipoles as opposed to what is commonly seen, for example, in isolated water, where the transition dipoles of the eigenstates are perpendicular to each other.

Ultrafast Vibrational Echo Spectroscopy of Liquid Water from First-Principles Simulations

Vibrational echo spectroscopy has become a powerful technique to study vibrational spectral diffusion in water and aqueous solutions. The dynamics of vibrational spectral diffusion is intimately related to the hydrogen bond fluctuations in liquid water and other hydrogen bonded liquids. Earlier theoretical calculations of vibrational echo spectroscopy of aqueous systems were based on classical molecular dynamics simulations involving empirical force fields of water. In the current work, we have employed the method of ab initio molecular dynamics simulation to calculate the spectral observables of vibrational echo and two-dimensional infrared (2D-IR) spectroscopy of liquid water at room temperature under Condon and cumulant approximations. The time scales extracted from the temporal decay of the frequency-time correlation function (FTCF), short-time slope of three pulse photon echo (SP3E), dynamic line width (DLW), and the slope of nodal line of 2D-IR spectra are found to be in reasonably close agreement with each other which reinforces the assertion that signatures of FTCF can be captured using three pulse photon echo and 2D-IR spectroscopy.

Collective vibrations of water-solvated hydroxide ions investigated with broadband 2DIR spectroscopy

The infrared spectra of aqueous solutions of NaOH and other strong bases exhibit a broad continuum absorption for frequencies between 800 and 3500 cm(-1), which is attributed to the strong interactions of the OH(-) ion with its solvating water molecules. To provide molecular insight into the origin of the broad continuum absorption feature, we have performed ultrafast transient absorption and 2DIR experiments on aqueous NaOH by exciting the O-H stretch vibrations and probing the response from 1350 to 3800 cm(-1) using a newly developed sub-70 fs broadband mid-infrared source. These experiments, in conjunction with harmonic vibrational analysis of OH(-)(H2O)n (n = 17) clusters, reveal that O-H stretch vibrations of aqueous hydroxides arise from coupled vibrations of multiple water molecules solvating the ion. We classify the vibrations of the hydroxide complex by symmetry defined by the relative phase of vibrations of the O-H bonds hydrogen bonded to the ion. Although broad and overlapping spectral features are observed for 3- and 4-coordinate ion complexes, we find a resolvable splitting between asymmetric and symmetric stretch vibrations, and assign the 2850 cm(-1) peak infrared spectra of aqueous hydroxides to asymmetric stretch vibrations.

Local orientational order in liquids revealed by resonant vibrational energy transfer

We demonstrate that local orientational ordering in a liquid can be observed in the decay of the vibrational anisotropy caused by resonant transfer of vibrational excitations between its constituent molecules. We show that the functional form of this decay is determined by the (distribution of) angles between the vibrating bonds of the molecules between which energy transfer occurs, and that the initial drop in the decay reflects the average angle between nearest neighbors. We use this effect to observe the difference in local orientational ordering in the two hydrogen-bonded liquids ethanol and N-methylacetamide.

OH-stretching in synthetic hydrogen-bonded chains

We study hydrogen bond dynamics in stereoselectively synthesized polyalcohols by combining linear and two-dimensional (2D) infrared spectroscopy experiments with simulations. We consider two variants of the polyalcohols: the all-syn and all-anti tetrol, which because of their different stereochemistry of the hydroxyl groups form a linear hydrogen-bonded chain that is stable for tens of picoseconds or a system where hydrogen bonds are formed and broken on a picosecond timescale, respectively. The differences in structure and hydrogen bond dynamics gives rise to significant differences in the linear spectra for the two compounds. Furthermore, we show that the stronger hydrogen bonding for the all-syn variant leads to faster fluctuations of the site frequencies than for the all-anti one, which is reflected in the higher degree of homogeneous broadening in the 2D spectra. Because of the different stereochemistry, the coupling in the all-syn molecule is stronger than for the all-anti one, which leads to a faster delocalization of a local excitation. This explains the previously observed pump-frequency independent vibrational lifetime for the all-syn variant, since the excitation loses the memory of the pump frequency before relaxation. For the all-anti form, the coupling is weak and the excitation remains in the initially excited state, maintaining the memory of the pump frequency.

Effects of excluded volume and correlated molecular orientations on Förster resonance energy transfer in liquid water

Förster theory for the survival probability of excited chromophores is generalized to include the effects of excluded volume and orientation correlation in the molecular distribution. An analytical expression for survival probability was derived and written in terms of a few simple elementary functions. Because of the excluded volume, the survival probability exhibits exponential decay at early times and stretched exponential decay at later times. Experimental schemes to determine the size of the molecular excluded volume are suggested. With the present generalization of theory, we analyzed vibrational resonance energy transfer kinetics in neat water. Excluded volume effects prove to be important and slow down the kinetics at early times. The majority of intermolecular resonance energy transfer was found to occur with exponential kinetics, as opposed to the stretched exponential behavior predicted by Förster theory. Quantum yields of intra-molecular vibrational relaxation, intra-, and intermolecular energy transfer were calculated to be 0.413, 0.167, and 0.420, respectively.

Ultrafast dynamics of liquid water: frequency fluctuations of the OH stretch and the HOH bend

Frequency fluctuations of the OH stretch and the HOH bend in liquid water are reported from the third-order response function evaluated using the TTM3-F potential for water. The simulated two-dimensional infrared spectra of the OH stretch are similar to previously reported theoretical results. The present study suggests that the frequency fluctuation of the HOH bend is faster than that of the OH stretch. The ultrafast loss of the frequency correlation of the HOH bend is due to the strong couplings with the OH stretch as well as the intermolecular hydrogen bond bend.

Slow hydrogen-bond switching dynamics at the water surface revealed by theoretical two-dimensional sum-frequency spectroscopy

Using our newly developed explicit three-body (E3B) water model, we simulate the surface of liquid water. We find that the timescale for hydrogen-bond switching dynamics at the surface is about three times slower than that in the bulk. In contrast, with this model rotational dynamics are slightly faster at the surface than in the bulk. We consider vibrational two-dimensional (2D) sum-frequency generation (2DSFG) spectroscopy as a technique for observing hydrogen-bond rearrangement dynamics at the water surface. We calculate the nonlinear susceptibility for this spectroscopy for two different polarization conditions, and in each case we see the appearance of cross-peaks on the timescale of a few picoseconds, signaling hydrogen-bond rearrangement on this timescale. We thus conclude that this 2D spectroscopy will be an excellent experimental technique for observing slow hydrogen-bond switching dynamics at the water surface.

An Efficient N(3)-Scaling Propagation Scheme for Simulating Two-Dimensional Infrared and Visible Spectra

In this paper, we develop and test a new approximate propagation scheme for calculating two-dimensional infrared and visible spectra. The new scheme scales one order more efficiently with the system size than the existing schemes. A Trotter type of approximation is used for the matrix exponent that describes the time evolution of the quantum system. This is needed for calculating the third-order response functions governing the two-dimensional infrared and visible spectra. The method is tested on a model dimer system, the amide I spectrum of the Gramicidin A antimicrobial peptide, the spectrum of the OH stretching vibration of bulk water, and a homogeneous J-aggregate. Errors due to the approximation are hardly observable in the calculated spectra. Scaling simulations with different system sizes are used to demonstrate that the new scheme is indeed scaling with the system size to the third power, one order faster than the existing methods.