Vibrational energy relaxation of a diatomic molecule in a room-temperature ionic liquid

Vibrational energy relaxation (VER) dynamics of a diatomic solute in ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate (EMI(+)PF(6) (-)) are studied via equilibrium and nonequilibrium molecular dynamics simulations. The time scale for VER is found to decrease markedly with the increasing solute dipole moment, consonant with many previous studies in polar solvents. A detailed analysis of nonequilibrium results shows that for a dipolar solute, dissipation of an excess solute vibrational energy occurs almost exclusively via the Lennard-Jones interactions between the solute and solvent, while an oscillatory energy exchange between the two is mainly controlled by their electrostatic interactions. Regardless of the anharmonicity of the solute vibrational potential, VER becomes accelerated as the initial vibrational energy increases. This is attributed primarily to the enhancement in variations of the solvent force on the solute bond, induced by large-amplitude solute vibrations. One interesting finding is that if a time variable scaled with the initial excitation energy is employed, dissipation dynamics of the excess vibrational energy of the dipolar solute tend to show a universal behavior irrespective of its initial vibrational state. Comparison with water and acetonitrile shows that overall characteristics of VER in EMI(+)PF(6) (-) are similar to those in acetonitrile, while relaxation in water is much faster than the two. It is also found that the Landau-Teller theory predictions for VER time scale obtained via equilibrium simulations of the solvent force autocorrelation function are in reasonable agreement with the nonequilibrium results.     

Non-monotonic, distance-dependent relaxation of water in reverse micelles: propagation of surface induced frustration along hydrogen bond networks

Layer-wise, distance-dependent orientational relaxation of water confined in reverse micelles (RM) is studied using theoretical and computational tools. We use both a newly constructed "spins on a ring" (SOR) Ising-type model (with Shore-Zwanzig rotational dynamics) and atomistic simulations with explicit water. Our study explores the effect of reverse micelle size and role of intermolecular correlations, compromised by the presence of a highly polar surface, on the distance (from the interface) dependence of water relaxation. The "spins on a ring" model can capture some aspects of distance dependence of relaxation, such as acceleration of orientational relaxation at intermediate layers. In atomistic simulations, layer-wise decomposition of hydrogen bond formation pattern clearly reveals that hydrogen bond arrangement of water at a certain distance away from the surface can remain frustrated due to the interaction with the polar surface head groups. This layer-wise analysis also reveals the presence of a non-monotonic slow relaxation component which can be attributed to this frustration effect and which is accentuated in small to intermediate size RMs. For large size RMs, the long time component decreases monotonically from the interface to the interior of the RMs with slowest relaxation observed at the interface.     

Theoretical studies of solute vibrational energy relaxation at liquid interfaces

Recent advances in the theoretical understanding of solute vibrational energy relaxation at liquid interfaces and surfaces are described. Non-equilibrium molecular dynamics simulations of the relaxation of an initially excited solute molecule are combined with equilibrium force autocorrelation calculations to gain insight into the factors that influence the vibrational relaxation rate. Diatomic and triatomic nonpolar, polar, and ionic solute molecules adsorbed at the liquid/vapor interface of several liquids as well as at the water/CCl(4) liquid/liquid interface are considered. In general, the vibrational relaxation rate is significantly slower (a factor of 3 to 4) at the liquid/vapor and liquid/liquid interface than in the bulk due to the reduced density, which gives rise to a reduced contribution of the repulsive solvent-solute forces on the vibrational mode. The surface effects on the ionic solutes are much smaller (50% or less slower relaxation relative to the bulk). This is due to the fact that ionic solutes at the interface are able to keep part of their solvation shell to a degree that depends on their size. Thus, a significant portion of the repulsive forces is maintained. A high degree of correlation is found between the peak height of the solvent-solute radial distribution function and the vibrational relaxation rate. The relaxation rate at the liquid/liquid interface strongly depends on the location of the solute across the interface and correlates with the change in the density and polarity profile of the interface.     

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.     

Structure, dynamics, and the free energy of solute adsorption at liquid-vapor interfaces of simple dipolar systems: molecular dynamics results for pure and mixed Stockmayer fluids

We have performed molecular dynamics simulation studies of the structural, thermodynamic, and dynamical properties of liquid-vapor interfaces of pure and binary Stockmayer fluids of different polarity. The density profiles, the width of the liquid-vapor interface, and the orientational structure of the interfaces are calculated to characterize the structural aspects of the interfaces. Among the thermodynamic properties, we have computed the surface tension and also the free energy of transfer of a charged solute across the liquid-vapor interface for both pure and mixed fluids. Among the dynamical properties of the interfaces, we have calculated the time dependence of the velocity and angular velocity autocorrelation functions, continuous and intermittent survival probabilities, mean square displacements, diffusion coefficients, and also the dipole correlation functions and orientational relaxation times of interfacial solvent molecules. It is found that the width of the interfaces decreases with increase of concentration of the more polar component. The dipole vectors of the interfacial molecules tend to align parallel to the surfaces and this alignment is enhanced with increasing dipole moment of the fluid molecules. Also, the surface tension shows an increasing trend with increase of dipole moment of the molecules. The dynamical properties of the interfaces are found to be different from those of the corresponding bulk liquid phases. In general, the molecules at the interfaces are found to rotate and translate in the parallel direction at a somewhat faster rate than the bulk molecules. Also, on increase of concentration of the more polar component, the diffusion and orientational relaxation of interfacial molecules are found to show a weaker slowing down than those of the bulk molecules, which can be attributed to the preferential presence of the more polar component in the bulk liquid regions. The temporal behavior of the interfacial survival probabilities reveals a decrease of the survival times with increasing polarity, which can be attributed to a corresponding decrease in the interfacial thickness. Results are presented for both continuous and intermittent survival times and the origins of their differences are discussed. The free energy calculations reveal no minimum at the interfaces for adsorption of a charged solute, which shows that the ions would prefer to stay in the interior of the liquid phases, rather than at interfaces, for these model dipolar systems.     

Manifestation of nonequilibrium initial conditions in molecular rotation: the generalized J-diffusion model

In order to adequately describe molecular rotation far from equilibrium, we have generalized the J-diffusion model by allowing the rotational relaxation rate to be angular momentum dependent. The calculated nonequilibrium rotational correlation functions (CFs) are shown to decay much slower than their equilibrium counterparts, and orientational CFs of hot molecules exhibit coherent behavior, which persists for several rotational periods. As distinct from the results of standard theories, rotational and orientational CFs are found to dependent strongly on the nonequilibrium preparation of the molecular ensemble. We predict the Arrhenius energy dependence of rotational relaxation times and violation of the Hubbard relations for orientational relaxation times. The standard and generalized J-diffusion models are shown to be almost indistinguishable under equilibrium conditions. Far from equilibrium, their predictions may differ dramatically.     

A first principles simulation study of fluctuations of hydrogen bonds and vibrational frequencies of water at liquid-vapor interface

We present a theoretical study of the structure and dynamics of water-vapor interface by means of ab initio molecular dynamics simulations. The inhomogeneous density, hydrogen bond and orientational profiles, voids and vibrational frequency distributions are investigated. We have also studied various dynamical properties of the interface such as diffusion, orientational relaxation, hydrogen bond dynamics and vibrational frequency fluctuations. The diffusion and orientational relaxation of water molecules are found to be faster at the interface which can be correlated with the voids present in the system. The hydrogen bond dynamics, however, is found to be slightly slower at the interface than that in bulk water. The correlations of hydrogen bond relaxation with the dynamics of vibrational frequency fluctuations are also discussed.     

Orientational motions of vibrational chromophores in molecules at the air/water interface with time-resolved sum frequency generation

The first time-resolved experiments in which interfacial molecules are pumped to excited electronic states and probed by vibrational sum frequency generation (SFG) are reported. This method was used to measure the out-of-plane rotation dynamics, i.e. time dependent changes in the polar angle, of a vibrational chromophore of an interfacial molecule. The chromophore is the carbonyl group, the rotation observed is that of the -C=O bond axis, with respect to the interfacial normal, and the interfacial molecule is coumarin 314 (C314) at the air/water interface. The orientational relaxation time was found to be 220+/-20 ps, which is much faster than the orientational relaxation time of the permanent dipole moment axis of C314 at the same interface, as obtained from pump-second harmonic probe experiments. Possible effects on the rotation of the -C=O bond axis due to the carbonyl group hydrogen bonding with interfacial water are discussed. From the measured equilibrium orientation of the permanent dipole moment axis and the carbonyl axis, and knowledge of their relative orientation in the molecule, the absolute orientation of C314 at the air/water interface is obtained.     

Solvent dynamics following charge transfer at the liquid-liquid interface

The dynamics of solvent reorganization following charge transfer at the interface between two immiscible liquids, one polar and the other non-polar are investigated by molecular dynamics. Both charge separation and charge recombination processes in two different orientations at the surface are considered and are compared with the same processes in the bulk polar solvent. The relaxation at the surface is significantly slower than in the bulk, in contrast to predictions of continuum models. The linear response approximation works reasonably well in the bulk but fails at the interface. The Langevin model is in poor agreement with the exact equilibrium correlations, especially at short time, but can qualitatively account for the overall rate.     

Molecular dynamics study of the liquid-vapor interface of acetonitrile: equilibrium and dynamical properties

The equilibrium and dynamical properties of the liquid-vapor interface of pure acetonitrile are studied by means of molecular dynamics simulations. Both nonpolarizable and polarizable models are employed in the present study. For the nonpolarizable model, the simulations are carried out for two different system sizes and at two different temperatures whereas the simulation with the polarizable model is done for a single system. The inhomogeneous density, anisotropic orientational profile, the width of the interface, and also the surface tension are calculated at room temperature and also at a lower temperature of 273 K. The dynamical aspects of the interface are investigated in terms of the single-particle dynamical properties such as the relaxation of velocity autocorrelation and the translational diffusion coefficients along the perpendicular and parallel directions and the dipole orientational relaxation of the interfacial acetonitrile molecules. The results of the interfacial dynamics are compared with those of the corresponding bulk phases at both temperatures. The convergence of the calculated results with respect to the length of simulation runs and the system size are also discussed.     

Structure and dynamics of water confined in silica nanopores

We report the results of molecular simulation of water in silica nanopores at full hydration and room temperature. The model systems are approximately cylindrical pores in amorphous silica, with diameters ranging from 20 to 40 ?. The filled pores are prepared using grand canonical Monte Carlo simulation and molecular dynamics simulation is used to calculate the water structure and dynamics. We found that water forms two distinct molecular layers at the interface and exhibits uniform, but somewhat lower than bulk liquid, density in the core region. The hydrogen bond density profile follows similar trends, with lower than bulk density in the core and enhancements at the interface, due to hydrogen bonds between water and surface non-bridging oxygens and OH groups. Our studies of water dynamics included translational mean squared displacements, orientational time correlations, survival probabilities in interfacial shells, and hydrogen bond population relaxation. We found that the radial-axial anisotropy in translational motion largely follows the predictions of a model of free diffusion in a cylinder. However, both translational and rotational water mobilities are strongly dependent on the proximity to the interface, with pronounced slowdown in layers near the interface. Within these layers, the effects of interface curvature are relatively modest, with only a small increase in mobility in going from the 20 to 40 ? diameter pore. Hydrogen bond population relaxation is nearly bulk-like in the core, but considerably slower in the interfacial region.     

Orientational dynamics of water trapped between two nanoscopic hydrophobic solutes: A molecular dynamics simulation study

We investigate thoroughly the effect of confinement and solute topology on the orientational dynamics of water molecule in the interplate region between two nanoscopic hydrophobic paraffinlike plates. Results are obtained from molecular dynamics simulations of aqueous solutions of paraffinlike plates in the isothermal-isobaric ensemble. An analysis of survival time auto correlation function shows that the residence time of the water molecule in the confined region between two model nanoscopic hydrophobic plates depends on solute surface topology (intermolecular distance within the paraffinlike plate). As expected, the extent of confinement also changes the residence time of water molecules considerably. Orientational dynamics was analyzed along three different directions, viz., dipole moment, HH, and perpendicular to molecular plane vectors. It has been demonstrated that the rotational dynamics of the confined water does not follow the Debye rotational diffusion model, and surface topology of the solute plate and the extent of confinement have considerable effect on the rotational dynamics of the confined water molecules.     

Molecular simulation study of water mobility in aerosol-OT reverse micelles

In this work, we present results from molecular dynamics simulations on the single-molecule relaxation of water within reverse micelles (RMs) of different sizes formed by the surfactant aerosol-OT (AOT, sodium bis(2-ethylhexyl)sulfosuccinate) in isooctane. Results are presented for RM water content w(0) = [H(2)O]/[AOT] in the range from 2.0 to 7.5. We show that translational diffusion of water within the RM can, to a good approximation, be decoupled from the translation of the RM through the isooctane solvent. Water translational mobility within the RM is restricted by the water pool dimensions, and thus, the water mean-squared displacements (MSDs) level off in time. Comparison with models of diffusion in confined geometries shows that a version of the Gaussian confinement model with a biexponential decay of correlations provides a good fit to the MSDs, while a model of free diffusion within a sphere agrees less well with simulation results. We find that the local diffusivity is considerably reduced in the interfacial region, especially as w(0) decreases. Molecular orientational relaxation is monitored by examining the behavior of OH and dipole vectors. For both vectors, orientational relaxation slows down close to the interface and as w(0) decreases. For the OH vector, reorientation is strongly affected by the presence of charged species at the RM interface and these effects are especially pronounced for water molecules hydrogen-bonded to surfactant sites that serve as hydrogen-bond acceptors. For the dipole vector, orientational relaxation near the interface slows down more than that for the OH vector due mainly to the influence of ion-dipole interactions with the sodium counterions. We investigate water OH and dipole reorientation mechanisms by studying the w(0) and interfacial shell dependence of orientational time correlations for different Legendre polynomial orders.     

Molecular rotation at negatively charged surfactant/aqueous interfaces

The effect of charge on the rotational dynamics of the molecular probe coumarin 314 (C314) at air/water interfaces covered with the negatively charged surfactant sodium dodecyl sulfate (SDS) was investigated using femtosecond time-resolved second harmonic spectroscopy. The out-of-plane orientational time constant at the highest SDS surface coverage of 100 A2 per molecule is 383 +/- 9 ps. The rotational dynamics is slower than at the air/water interface where the out-of-plane reorientational time constant is 336 +/- 6 ps. At the air/water interface the rotational dynamics is over three times slower than the bulk orientational diffusion time of 100 ps. The relatively small effect of the surfactant charge density on the C314 rotation time constant is surprising, considering the marked dependence of the C314 orientation, spectra, and surfactant phase diagram on surfactant density.     

Distribution, structure, and dynamics of cesium and iodide ions at the H2O-CCl4 and H2O-vapor interfaces

Molecular dynamics simulations utilizing many-body potentials of H2O-CCl4 and H2O-vapor interfaces were carried out at different cesium and iodide ion concentrations to compare ion distribution, interfacial orientational and structural properties, and dynamics. It was found that cesium was repelled by both interfaces, and iodide was active at both interfaces, but to a much greater degree at the H2O-vapor interface. At the interface, the iodide dipole was strongly induced, orienting perpendicular to the interface for both systems, leading to stronger hydrogen bonds with water. For the H2O-CCl4 interface, though, there was a compensation between these strong hydrogen bonds and short to moderate ranged repulsion between iodide and CCl4. Hydrogen bond distance and angular distributions showed weaker water-water hydrogen bonds at both interfaces, but generally stronger water-iodide hydrogen bonds. Both translational and rotational dynamics of water were faster at the interface, while for CCl4, its translational dynamics was slower, but rotational dynamics faster at the interface. For many of the studied systems and species, translational diffusion was found to be anisotropic at both interfacial and bulk regions.     

Orientational dynamics for an amphiphilic-solvent solution

In this work, we performed Monte Carlo simulations on a lattice model for spontaneous amphiphilic aggregation, in order to study the orientational and hydrogen-bonding dynamics of water on different regions inside the micellar solution. We employed an associating lattice gas model that mimics the aqueous solvent, which presents a rich phase diagram with first- and second-order transition lines. Even though this is a simplified model, it makes possible to investigate the orientational dynamics of water in an equilibrium solution of amphiphiles, as well as the influence of the different phases of the solvent in the interfacial and bulk water dynamics. By means of extensive simulations, we showed that, at high temperatures, the behavior of the orientational relaxation and hydrogen bonding of water molecules in the bulk, first, and second hydration shells are considerable different. We observe the appearance of a very slow component for water molecules in the first hydration shell of micelles when the system reaches a high-density phase, consistent with previous theoretical and experimental studies concerning biological water. Also, at high temperatures, we find that water molecules in the second hydration shell of micelles have an orientational decay similar to that of bulk water, but with a generally slower dynamics. Otherwise, at low temperatures, we have two components for the orientational relaxation of bulk water in the low density liquid phase, and only a single component in the high density liquid (HDL) phase, which reflect the symmetry properties of the different phases of the solvent model. In the very dense region of water molecules in the first hydration shell of micelles at low temperatures, we find two components for the orientational relaxation on both liquid phases, one of them much slower than that in the single component of bulk water in the HDL phase. This happens even though our model does not present any hindrance to the water rotational freedom caused by the presence of the amphiphiles.     

Hydrogen bonded structure and dynamics of liquid-vapor interface of water-ammonia mixture: an ab initio molecular dynamics study

We have carried out ab initio molecular dynamics simulations of a liquid-vapor interfacial system consisting of a mixture of water and ammonia molecules. We have made a detailed analysis of the structural and dynamical properties of the bulk and interfacial regions of the mixture. Among structural properties, we have looked at the inhomogeneous density profiles of water and ammonia molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions of bulk and interfacial molecules. It is found that the interfacial molecules show preference for specific orientations so as to form water-ammonia hydrogen bonds at the interface with ammonia as the acceptor. The structure of the system is also investigated in terms of inter-atomic voids present in the system. Among the dynamical properties, we have calculated the diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion in bulk and interfacial regions. It is found that the diffusion and orientation relaxation of the interfacial molecules are faster than those of the bulk. However, the hydrogen bond lifetimes are longer at the interface which can be correlated with the time scales found from the decay of frequency time correlations.     

Water dynamics--the effects of ions and nanoconfinement

Hydrogen bond dynamics of water in highly concentrated NaBr salt solutions and reverse micelles are studied using ultrafast 2D-IR vibrational echo spectroscopy and polarization-selective IR pump-probe experiments performed on the OD hydroxyl stretch of dilute HOD in H(2)O. The vibrational echo experiments measure spectral diffusion, and the pump-probe experiments measure orientational relaxation. Both experimental observables are directly related to the structural dynamics of water's hydrogen bond network. The measurements performed on NaBr solutions as a function of concentration show that the hydrogen bond dynamics slow as the NaBr concentration increases. The most pronounced change is in the longest time scale dynamics which are related to the global rearrangement of the hydrogen bond structure. Complete hydrogen bond network randomization slows by a factor of approximately 3 in approximately 6 M NaBr solution compared to that in bulk water. The hydrogen bond dynamics of water in nanoscopically confined environments are studied by encapsulating water molecules in ionic head group (AOT) and nonionic head group (Igepal CO 520) reverse micelles. Water dynamics in the nanopools of AOT reverse micelles are studied as a function of size by observing orientational relaxation. Orientational relaxation dynamics deviate significantly from bulk water when the size of the reverse micelles is smaller than several nm and become nonexponential and slower as the size of the reverse micelles decreases. In the smallest reverse micelles, orientational relaxation (hydrogen bond structural randomization) is almost 20 times slower than that in bulk water. To determine if the changes in dynamics from bulk water are caused by the influence of the ionic head groups of AOT or the nanoconfinement, the water dynamics in 4 nm nanopools in AOT reverse micelles (ionic) and Igepal reverse micelles (nonionic) are compared. It is found that the water orientational relaxation in the 4 nm diameter nanopools of the two types of reverse micelles is almost identical, which indicates that confinement by an interface to form a nanoscopic water pool is a primary factor governing the dynamics of nanoscopic water rather than the presence of charged groups at the interface.