Probing intermolecular couplings in liquid water with two-dimensional infrared photon echo spectroscopy

Simulating the isotropic Raman spectra of O-H stretching mode in liquid H2O based on a machine learning potential: the influence of vibrational couplings

In this study, we trained a deep potential (DP) for H2O, an accurate machine learning (ML) potential. We performed molecular dynamics (MD) simulations of liquid water using the DP model (or DeePMD simulations). Our results showed that the DP model exhibits DFT-level accuracy, and the DeePMD simulation is a promising approach for modeling the structural properties of liquid water. Based on the DeePMD simulation trajectories, we calculated the isotropic Raman spectra of the O-H stretching mode using the surface-specific velocity-velocity correlation function (ssVVCF), showing that the DeePMD/ssVVCF approach can correctly capture the bimodal characteristics of the experimental Raman spectra, with one peak located near 3400 cm-1 and the other near 3250 cm-1. The success of the DeePMD/ssVVCF approach should be credited to (1) the DFT-level accuracy of the DP model for H2O, (2) the ssVVCF formulation considering the coupling between vibrational modes, and (3) non-Condon effects. Furthermore, the DeePMD simulations revealed that the anharmonic interactions between the coupled water molecules in the first and second hydration shells should play an essential role in the strong mixing of the H-O-H bending mode and the O-H stretching mode, leading to the delocalization of the O-H stretching band. In particular, increasing the strength of hydrogen bonds would enhance the bend-stretch coupling, leading to the red-shifting of the O-H vibrational spectra and the increase in the intensity of the shoulder around 3250 cm-1.

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.

Role of Intermolecular Charge Fluxes in the Hydrogen-Bond-Induced Frequency Shifts of the OH Stretching Mode of Water

The relation between the vibrational properties and the electrostatic situations of the vibrating functional group is useful to predict vibrational spectroscopic features based on, for example, classical molecular dynamics of liquids or biomolecular systems, but to pursue its generality or the extent of applicability, it is required to understand the mechanisms giving rise to it. Here such an analysis is carried out for the OH stretching mode of water. By examining the correlations among various (structural, vibrational, and electrostatic) properties and by analyzing the spatial characteristics of the behavior of electrons occurring upon the vibration, it is shown that the dependence of the vibrational frequency and the dipole derivative of the OH stretching mode on the electric field is not of purely electrostatic origin, and the delocalized electronic motions occurring with this mode, called intermolecular charge fluxes, related to both the dipole first and second derivatives play important roles.

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.

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.

Vibration Spectral Dynamics of Weakly Coordinating Water Molecules near an Anion: FPMD Simulations of an Aqueous Solution of Tetrafluoroborate

The extent to which the ions affect the nearby water molecules will decide the structure-making or breaking nature of those ions in aqueous solutions. The effects of a weakly coordinating anion on the structure, dynamics, and vibrational properties of water molecules are not so significant as compared to an anion capable of making strong ion-water hydrogen bonds. The present work deals with the first-principles molecular dynamics study of an aqueous solution of such a weakly coordinating anion, tetrafluoroborate (BF₄^-), using dispersion-corrected DFT-based first-principles molecular dynamics (FPMD) simulations. Various structural, dynamical, and spectral properties, such as radial distribution functions (RDFs), rotational dynamics, vibrational density of states (VDOS), hydrogen bond as well as dangling OH autocorrelation functions, and residence dynamics, were calculated to investigate the effects of the anion on nearby water molecules. The process of spectral diffusion was assessed through a time series wavelet transformation of trajectories obtained from FPMD simulations. The first ion-water solvation shell extends up to 5.5 Å, containing around 20 water molecules. The lifetime of the ion-water hydrogen bond is found to be 1.19 ps, whereas the water-water hydrogen bond lifetime is found to be 1.13 ps. Inside the solvation shell, the persistence time of dangling OH chromophores and the average frequency of OH modes inside the solvation shell are found to be more compared to bulk. Three time scales are found for solvation shell OH modes from the frequency-frequency correlation function. A very short time scale is found for the intact ion-water interaction; the short time scale is for the ion-water hydrogen bond, and the long time scale is for escape dynamics of water molecules from the ion solvation shell. From the mean squared displacement, it is found that solvation water molecules diffuse slower than the bulk. However, solvation shell water molecules show faster relaxation from the analysis of rotational anisotropy. Within the longer time scale of spectral diffusion, this process (which is related to various dynamics of the molecules) is not yet complete, as compared to fast anisotropic decay. This fact is similar to the experimental finding of spectral diffusion and anisotropy time scales in the aqueous solution of borohydride anion. The calculated results are also compared with available experimental data wherever possible.

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.

Differences in the Vibrational Dynamics of H(2)O and D(2)O: Observation of Symmetric and Antisymmetric Stretching Vibrations in Heavy Water

Water's ability to donate and accept hydrogen bonds leads to unique and complex collective dynamical phenomena associated with its hydrogen-bond network. It is appreciated that the vibrations governing liquid water's molecular dynamics are delocalized, with nuclear motion evolving coherently over the span of several molecules. Using two-dimensional infrared spectroscopy, we have found that the nuclear motions of heavy water, D2O, are qualitatively different than those of H2O. The nonlinear spectrum of liquid D2O reveals distinct O-D stretching resonances, in contrast to H2O. Furthermore, our data indicates that condensed-phase O-D vibrations have a different character than those in the gas phase, which we understand in terms of weakly delocalized symmetric and antisymmetric stretching vibrations. This difference in molecular dynamics reflects the shift in the balance between intra- and intermolecular couplings upon deuteration, an effect which can be understood in terms of the anharmonicity of the nuclear potential energy surface.

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

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

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.

Manipulating Excited-State Dynamics of Individual Light-Harvesting Chromophores through Restricted Motions in a Hydrated Nanoscale Protein Cavity

Manipulating the photophysical properties of light-absorbing units is a crucial element in the design of biomimetic light-harvesting systems. Using a highly tunable synthetic platform combined with transient absorption and time-resolved fluorescence measurements and molecular dynamics simulations, we interrogate isolated chromophores covalently linked to different positions in the interior of the hydrated nanoscale cavity of a supramolecular protein assembly. We find that, following photoexcitation, the time scales over which these chromophores are solvated, undergo conformational rearrangements, and return to the ground state are highly sensitive to their position within this cavity and are significantly slower than in a bulk aqueous solution. Molecular dynamics simulations reveal the hindered translations and rotations of water molecules within the protein cavity with spatial specificity. The results presented herein show that fully hydrated nanoscale protein cavities are a promising way to mimic the tight protein pockets found in natural light-harvesting complexes. We also show that the interplay between protein, solvent, and chromophores can be used to substantially tune the relaxation processes within artificial light-harvesting assemblies in order to significantly improve the yield of interchromophore energy transfer and extend the range of excitation transport. Our observations have implications for other important, similarly sized bioinspired materials, such as nanoreactors and biocompatible targeted delivery agents.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. III. Water clustering and vibrational energy transfer

Water clustering and connectivity around lipid bilayers strongly influences the properties of membranes and is important for functions such as proton and ion transport. Vibrational anisotropic pump-probe spectroscopy is a powerful tool for understanding such clustering, as the measured anisotropy depends upon the time-scale and degree of intra- and intermolecular vibrational energy transfer. In this article, we use molecular dynamics simulations and theoretical vibrational spectroscopy to help interpret recent experimental measurements of the anisotropy of water in lipid multi-bilayers as a function of both lipid hydration level and isotopic substitution. Our calculations are in satisfactory agreement with the experiments of Piatkowski, Heij, and Bakker, and from our simulations we can directly probe water clustering and connectivity. We find that at low hydration levels, many water molecules are in fact isolated, although up to 70% of hydration water forms small water clusters or chains. At intermediate hydration levels, water forms a wide range of cluster sizes, while at higher hydration levels, the majority of water molecules are part of a large, percolating water cluster. Therefore, the size, number, and nature of water clusters are strongly dependent on lipid hydration level, and the measured anisotropy reflects this through its dependence on intermolecular energy transfer.

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.

Fluctuations and Relaxation Dynamics of Liquid Water Revealed by Linear and Nonlinear Spectroscopy

Many efforts have been devoted to elucidating the intra- and intermolecular dynamics of liquid water because of their important roles in many fields of science and engineering. Nonlinear spectroscopy is a powerful tool to investigate the dynamics. Because nonlinear response functions are described by more than one time variable, it is possible to analyze static and dynamic mode couplings. Here we review the intra- and intermolecular dynamics of liquid water revealed by recent linear and nonlinear spectroscopic experiments and computer simulations. In particular, we discuss the population relaxation, anisotropy decay, and spectral diffusion of the intra- and intermolecular motions of water and their temperature dependence, which play important roles in ultrafast dynamics and relaxations in water.

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.

Probing the distribution of water molecules hydrating lipid membranes with ultrafast Förster vibrational energy transfer

We determine the relative positioning of water molecules in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) membranes by measuring the rate of vibrational resonant (Förster) energy transfer between the water hydroxyl stretch vibrations. The rate of Förster energy transfer is strongly distance dependent and thus gives detailed information on the relative positioning of the water molecules. We determine the rate of intermolecular Förster energy by measuring the anisotropy dynamics of excited O-D stretch vibrations of HDO and D(2)O molecules with polarization-resolved femtosecond mid-infrared spectroscopy. We study the dynamics for deuterium fractions between 0.1 and 1 and for hydration levels between 2 and 12 water molecules per DOPC lipid molecule. We find that most of the water molecules hydrating the membrane are contained in nanoclusters and have an average intermolecular distance of 3.4 Å. The density of these nanoclusters increases with increasing hydration level of the DOPC membranes.

Frequency distribution of the amide-I vibration sorted by residues in amyloid fibrils revealed by 2D-IR measurements and simulations

The infrared optical response of amyloid fibrils Aβ(1-40) is investigated. Simulations of two models corresponding to different protonation states are compared with experiment. The simulations reveal that vibrational frequency distributions inside the fibrils are dominated by side chain fluctuations. We further confirm earlier suggestions based on 2D-IR measurements that water molecules can be trapped inside the fibrils.

Decelerated water dynamics and vibrational couplings of hydrated DNA mapped by two-dimensional infrared spectroscopy

Double-stranded DNA oligomers containing 23 alternating adenine-thymine base pairs are studied at different hydration levels by femtosecond two-dimensional (2D) infrared spectrosopy. Coupled NH stretching modes of the A-T pairs and OH stretching excitations of the water shell are discerned in the 2D spectra. Limited changes of NH stretching frequencies and line shapes with increasing hydration suggest spectral dynamics governed by DNA rather than water fluctuations. In contrast, OH stretching excitations of the water shell around fully hydrated DNA undergo spectral diffusion on a ~500 fs time scale. The center line slopes of the 2D spectra of hydrated DNA demonstrate a slower decay of the frequency-time correlation function (TCF) than that in neat water, as is evident from a comparison with 2D spectra of neat H(2)O and theoretical TCFs. We attribute this behavior to reduced structural fluctuations of the water shell and a reduced rate of resonant OH stretching energy transfer.

Reduced coupling of water molecules near the surface of reverse micelles

We report on vibrational dynamics of water near the surface of AOT reverse micelles studied by narrow-band excitation, mid-IR pump-probe spectroscopy. Evidence of OH-stretch frequency splitting into the symmetric and asymmetric modes is clearly observed for the interfacial H(2)O molecules. The polarization memory of interfacial waters is preserved over an exceptionally extended >10 ps timescale which is a factor of 100 longer than in bulk water. These observations point towards negligibly small intermolecular vibrational coupling between the water molecules as well as strongly reduced water rotational mobility within the interfacial water layer.

Strong intermolecular coupling between the HF stretching and H2O bending vibrations in HF:H2O binary amorphous solids: breakdown of the electrostatic description of the hydrogen bond

The coupling mechanism between the HF stretching and H(2)O bending vibrations observed in the infrared spectra of HF:H(2)O binary amorphous solids is analyzed using a simple cluster model. The intermolecular vibrational coupling derived from electrostatic potentials is one order of magnitude smaller, and of the opposite sign, than that obtained from electronic structure-based potentials. This highlights the distinctively covalent character of strong H-bonds and unveils fundamental weaknesses of electrostatic descriptions of vibrational energy transfer in liquid water and aqueous solutions.

Intra- and intermolecular charge fluxes induced by the OH stretching mode of water and their effects on the infrared intensities and intermolecular vibrational coupling

Intra- and intermolecular charge fluxes induced by the OH stretching mode of water and their effects on the infrared (IR) intensities and intermolecular vibrational coupling are analyzed theoretically. Density functional theoretical (DFT) calculations are carried out for the clusters consisting of 28 and 30 water molecules, and the change in the electron density induced by the OH stretching mode, i.e., the electron density derivative ∂ρ((el))(r)/∂Q(OH), is calculated and its effect is analyzed. It is shown that the electron density derivative extends to the spatial region of the water molecule that accepts the hydrogen bond donated by the vibrating OH bond, indicating that intermolecular charge flux is induced. Both the amount of vibration-induced change in the electric charge (the charge derivative) and the distance between charge centers are large for the intermolecular charge flux, giving rise to a significant enhancement of the dipole derivative, and hence the IR intensity. After showing quantitatively that it is possible to understand the intermolecular vibrational coupling between the OH stretching oscillators with the picture of the "electric field modulation by one vibration" and the "response by the other vibration to the electric field modulation", the effect of the intermolecular charge flux on the intermolecular vibrational coupling is analyzed. It is shown that the intermolecular charge flux plays an important role also in this relation, but the two factors determining the dipole derivative (the charge derivative and the distance between charge centers) should be carefully taken into account.

Solvent-dependent spectral diffusion in a hydrogen bonded "vibrational aggregate"

Two-dimensional infrared spectroscopy (2DIR) is used to measure the viscosity-dependent spectral diffusion of a model vibrational probe, Mn(2)(CO)(10) (dimanganese decacarbonyl, DMDC), in a series of alcohols with time scales ranging from 2.67 ps in methanol to 5.33 ps in 1-hexanol. Alcohol-alkane solvent mixtures were found to produce indistinguishable linear IR spectra, while still demonstrating viscosity-dependent spectral diffusion. Using a vibrational exciton model to characterize the inhomogeneous energy landscape, several analogies emerge with multichromophoric electronic systems, such as J-aggregates and light-harvesting protein complexes. An excitonic, local vibrational mode Hamiltonian parametrized to reproduce the vibrational structure of DMDC serves as a starting point from which site energies (i.e., local carbonyl frequencies) are given Gaussian distributed disorder. The model gives excellent agreement with both the linear IR spectrum and the inhomogeneous widths extracted from 2DIR, indicating the system can be considered to be a "vibrational aggregate." This model naturally leads to exchange narrowing due to disorder-induced exciton localization, producing line widths consistent with our 1D and 2D measurements. Further, the diagonal disorder alone effectively reduces the molecular symmetry, leading to the appearance of Raman bands in the IR spectrum in accord with the measurements. Here, we show that the static inhomogeneity of the excitonic model with disorder successfully captures the essential details of the 1D spectrum while predicting the degree of IR activity of forbidden modes as well as the inhomogeneous widths and relative magnitudes of the transition moments.

Two-dimensional terahertz correlation spectra of electronic excitations in semiconductor quantum wells

We discuss a novel approach for nonlinear two-dimensional (2D) spectroscopy in the terahertz (THz) frequency range which is based on a collinear interaction geometry of a sequence of THz pulses with the sample. The nonlinear polarization is determined by a phase-resolved measurement of the electric field transmitted through the sample as a function of the delay τ between two phase-locked pulses and the "real" time t. The information provided by a single 2D scan along the τ and t axes is equivalent to that from a noncollinear photon-echo setup equipped with four local oscillators, each interacting with a different diffracted order. We address basic concepts of collinear 2D THz spectroscopy, in particular data analysis and phasing issues. Different rephasing and nonrephasing contributions to the third-order response are separated and 2D correlation spectra derived. As a prototype application, 2D correlation spectra of intersubband excitations of electrons in semiconductor quantum wells are presented.

Effects of nonadditive interactions on ion solvation at the water/vapor interface: a molecular dynamics study

The solvation of halide ions at the water/vapor interface is investigated by using molecular dynamics simulations with nonpolarizable molecular mechanical (MM), polarizable MM, and quantum mechanical (QM)/MM methods. The free energy profile of the ion solvation is decomposed into the energy and the entropic contributions along the ion displacement from inside to the surface of water. It is found that the surface affinity of the ion, relative to the bulk value, is determined by a subtle balance between the energetic destabilization and the entropic stabilization with the ion displacement. The amount of energetic destabilization is found to be reduced when nonadditive interactions are included, as in the polarizable MM and QM/MM models. The structure of water around the ion at the interface is also largely modified when the higher order effects are considered. For example, the induced dipole effect enhances the solvation structure around the ion at the interface significantly and thus reduces the amount of entropic stabilization at the interface, relative to in the bulk. It is found that this induced dipole effect causes the slowing in the ion-water hydrogen bond dynamics at the interface. On the other hand, the higher order induced multipole effects in the QM/MM method suppress both the excessive enhancement of the solvation structure and the slowing of the ion-water hydrogen bond dynamics at the interface. The present study demonstrates that not only the induced dipole moment but also the higher order induced multipole moments, which are neglected in standard empirical models, are essential for the correct description of the ion solvation at the water/vapor interface.