A novel method for analyzing energy relaxation in condensed phases using nonequilibrium molecular dynamics simulations: application to the energy relaxation of intermolecular motions in liquid water

We present a novel method to investigate energy relaxation processes in condensed phases using nonequilibrium molecular dynamics simulations. This method can reveal details of the time evolution of energy relaxation like two-color third-order IR spectroscopy. Nonetheless, the computational cost of this method is significantly lower than that of third-order response functions. We apply this method to the energy relaxation of intermolecular motions in liquid water. We show that the intermolecular energy relaxation in water is characterized by four energy transfer processes. The structural changes of the liquid associated with the energy relaxation are also analyzed by the nonequilibrium molecular dynamics technique.     

Influence of hydration on protein dynamics: combining dielectric and neutron scattering spectroscopy data

Combining dielectric spectroscopy and neutron scattering data for hydrated lysozyme powders, we were able to identify several relaxation processes and follow protein dynamics at different hydration levels over a broad frequency and temperature range. We ascribe the main dielectric process to protein's structural relaxation coupled to hydration water and the slowest dielectric process to a larger scale protein's motions. Both relaxations exhibit a smooth, slightly super-Arrhenius temperature dependence between 300 and 180 K. The temperature dependence of the slowest process follows the main dielectric relaxation, emphasizing that the same friction mechanism might control both processes. No signs of a proposed sharp fragile-to-strong crossover at T approximately 220 K are observed in temperature dependences of these processes. Both processes show strong dependence on hydration: the main dielectric process slows down by six orders with a decrease in hydration from h approximately 0.37 (grams of water per grams of protein) to h approximately 0.05. The slowest process shows even stronger dependence on hydration. The third (fastest) dielectric relaxation process has been detected only in samples with high hydration ( h approximately 0.3 and higher). We ascribe it to a secondary relaxation of hydration water. The mechanism of the protein dynamic transition and a general picture of the protein dynamics are discussed.     

Ultrafast intermolecular dynamics of liquid water: a theoretical study on two-dimensional infrared spectroscopy

Physical and chemical properties of liquid water are dominated by hydrogen bond structure and dynamics. Recent studies on nonlinear vibrational spectroscopy of intramolecular motion provided new insight into ultrafast hydrogen bond dynamics. However, our understanding of intermolecular dynamics of water is still limited. We theoretically investigated the intermolecular dynamics of liquid water in terms of two-dimensional infrared (2D IR) spectroscopy. The 2D IR spectrum of intermolecular frequency region (<1000 cm(-1)) is calculated by using the equilibrium and nonequilibrium hybrid molecular dynamics method. We find the ultrafast loss of the correlation of the libration motion with the time scale of approximately 110 fs. It is also found that the energy relaxation from the libration motion to the low frequency motion takes place with the time scale of about 180 fs. We analyze the effect of the hindered translation motion on these ultrafast dynamics. It is shown that both the frequency modulation of libration motion and the energy relaxation from the libration to the low frequency motion significantly slow down in the absence of the hindered translation motion. The present result reveals that the anharmonic coupling between the hindered translation and libration motions is essential for the ultrafast relaxation dynamics in liquid water.     

Relation between solvent and protein dynamics as studied by dielectric spectroscopy

We present results obtained by dielectric spectroscopy in wide frequency (10(-2)-10(9) Hz) and temperature ranges on human hemoglobin in the three different solvents water, glycerol, and methanol, at a solvent level of 0.8 g of solvent/g of protein. In this broad frequency region, there are motions on several time-scales in the measured temperature range (110-370 K for water, 170-410 K for glycerol, and 110-310 K for methanol). For all samples, the dielectric data shows at least four relaxation processes, with frequency dependences that are well described by the Havriliak-Negami or Cole-Cole functions. The fastest and most pronounced process in the dielectric spectra of hemoglobin in glycerol and methanol solutions is similar to the alpha-relaxation of the corresponding bulk solvent (but shifted to slower dynamics due to surface interactions). For water solutions, however, this process corresponds to earlier results obtained for water confined in various systems and it is most likely due to a local beta-relaxation. The slowing down of the glycerol and methanol relaxations and the good agreement with earlier results on confined water show that this process is affected by the interaction with the protein surface. The second fastest process is attributed to motions of polar side groups on the protein, with a possible contribution from tightly bound solvent molecules. This process is shifted to slower dynamics with increasing solvent viscosity, and it shows a crossover in its temperature dependence from Arrhenius behavior at low temperatures to non-Arrhenius behavior at higher temperatures where there seems to be an onset of cooperativity effects. The origins of the two slowest relaxation processes (visible at high temperatures and low frequencies), which show saddlelike temperature dependences for the solvents water and methanol, are most likely due to motions of the polypeptide backbone and an even more global motion in the protein molecule.     

On the Debye-Waller factor of hexagonal ice: a computer simulation study

We investigate by molecular dynamics (MD) simulations the temperature dependence of the Debye-Waller (DW) factor of hexagonal ice with 25 different proton-disordered configurations. Each initial configuration is composed of 288 water molecules with no net dipole moment. The intermolecular interaction of water is described by TIP4P potential. Each production run of the simulation is 15 ns or longer. We observe a change in slope of the DW factor around 200 K, which cannot be explained within the framework of either classical or quantum harmonic approximation. Configurations generated by MD simulations are subjected to the steepest descent energy minimization. Analysis of the local energy minimum structures reveals that water molecules above 200 K jump to other lattice sites via some local energy minimum structures which contain some water molecules sitting on the locations other than the lattice sites. As time evolves, these defect molecules move back and forth to the lattice sites yielding defect-free structures. Those motions are responsible for the unusual increase in the DW factor at high temperatures. In making a transition from an energy-minimum structure to another one, a small number of water molecules are involved in a highly cooperative fashion. The larger DW factor at higher temperature arises from jump-like motions of water molecules among these locally stable configurations which may or may not be a family of the proton-disordered ice forms satisfying the "ice rule".     

Crystallization, melting, and structure of water nanoparticles at atmospherically relevant temperatures

Water nanoparticles play an important role in atmospheric processes, yet their equilibrium and nonequilibrium liquid-ice phase transitions and the structures they form on freezing are not yet fully elucidated. Here we use molecular dynamics simulations with the mW water model to investigate the nonequilibrium freezing and equilibrium melting of water nanoparticles with radii R between 1 and 4.7 nm and the structure of the ice formed by crystallization at temperatures between 150 and 200 K. The ice crystallized in the particles is a hybrid form of ice I with stacked layers of the cubic and hexagonal ice polymorphs in a ratio approximately 2:1. The ratio of cubic ice to hexagonal ice is insensitive to the radius of the water particle and is comparable to that found in simulations of bulk water around the same temperature. Heating frozen particles that contain multiple crystallites leads to Ostwald ripening and annealing of the ice structures, accompanied by an increase in the amount of ice at the expense of the liquid water, before the particles finally melt from the hybrid ice I to liquid, without a transition to hexagonal ice. The melting temperatures T(m) of the nanoparticles are not affected by the ratio of cubic to hexagonal layers in the crystal. T(m) of the ice particles decreases from 255 to 170 K with the particle size and is well described by the Gibbs-Thomson equation, T(m)(R) = T(m)(bulk) - K(GT)/(R - d), with constant K(GT) = 82 ± 5 K·nm and a premelted liquid of width d = 0.26 ± 0.05 nm, about one monolayer. The freezing temperatures also decrease with the particles' radii. These results are important for understanding the composition, freezing, and melting properties of ice and liquid water particles under atmospheric conditions.     

Molecular dynamics with quantum transitions study of the vibrational relaxation of the HOD bend fundamental in liquid D2O

The molecular dynamics with quantum transitions method is used to study the vibrational relaxation of the HOD bend fundamental in liquid D(2)O. All of the vibrational bending degrees of freedom of the HOD and D(2)O molecules are described by quantum mechanics, while the remaining translational and rotational degrees of freedom are described classically. The effect of the coupling between the rotational and vibrational degrees of freedom of the deuterated water molecules is analyzed. A kinetic mechanism based on three steps is proposed in order to interpret the dynamics of the system. It is shown that intermolecular vibrational energy transfer plays an important role in the relaxation process and also that the transfer of energy into the rotational degrees of freedom is favored over the transfer of energy into the translational motions. The thermalization of the system after the relaxation is reached in a shorter time scale than that of the recovery of the hydrogen bond network. The relaxation and equilibration times obtained compare well with experimental and previous theoretical results.     

High frequency dynamics and structural relaxation process in liquid ammonia

The dynamic structure factor S(Q,omega) of liquid ammonia has been measured by inelastic x-ray scattering in the terahertz frequency region as a function of the temperature in the range of 220-298 K at a pressure P=85 bars. The data have been analyzed using the generalized hydrodynamic formalism with a three term memory function to take into account the thermal, the structural, (alpha) and the microscopic (mu) relaxation processes affecting the dynamics of the liquid. This allows to extract the temperature dependence of the structural relaxation time (tau(alpha)) and strength (Delta(alpha)). The former quantity follows an Arrhenius behavior with an activation energy E(a)=2.6+/-0.2 kcal/mol, while the latter is temperature independent suggesting that there are no changes in the interparticle potential and arrangement with T. The obtained results, compared with those already existing in liquid water and liquid hydrogen fluoride, suggest the strong influence of the connectivity of the molecular network on the structural relaxation.     

Water dynamics in n-propylene glycol aqueous solutions

The relaxation dynamics of dipropylene glycol and tripropylene glycol (nPG-n=2,3) water solutions on the nPG-rich side has been studied by broadband dielectric spectroscopy and differential scanning calorimetry in the temperature range of 130-280 K. Two relaxation processes are observed for all the hydration levels; the slower process (I) is related to the alpha relaxation of the solution whereas the faster one (II) is associated with the reorientation of water molecules in the mixture. Dielectric data for process (II) at temperatures between 150 and 200 K indicate the existence of a critical water concentration (x(c)) below which water mobility is highly restricted. Below x(c), nPG-water domains drive the dielectric signal whereas above x(c), water-water domains dominate the dielectric response at low temperatures. The results also show that process (II) at low temperatures is due to local motions of water molecules in the glassy frozen matrix. Additionally, we will show that the glass transition temperatures (T(g)) for aqueous PG, 2PG, and 3PG solutions do not extrapolate to approximately 136 K, regardless of the extrapolation method. Instead, we find that the extrapolated T(g) value for water from these solutions lies in the neighborhood of 165 K.     

No fragile-to-strong crossover in LiCl-H2O solution

Dynamics of water, especially in the temperature range of the "no man's land", remain a mystery. We present detailed study of dynamics in aqueous LiCl solution that is often considered as a model for bulk water. We employ broadband dielectric and light scattering spectroscopy in a broad frequency and temperature range. Our analysis reveals no sign of the fragile-to-strong crossover (FSC) neither in structural relaxation nor in translational motions. Our experimental results combined with a large selection of literature data lead to the clear conclusion-there is no FSC in dynamics of aqueous solutions at T ～ 200-230 K. Instead, our analysis reveals appearance of the so-called excess wing at the high frequency tail of the structural relaxation peak. We discuss the localized nature of the relaxation process that contributes to the excess wing.     

Further evidence that interfacial water is the main "driving force" of protein dynamics: a neutron scattering study on perdeuterated C-phycocyanin

The fundamental role of hydration water (also called interfacial water) is widely recognized in protein flexibility, especially in the existence of the so-called protein "dynamical transition" at around 220 K. In the present study, we take advantage of perdeuterated C-phycocyanin (CPC) and elastic incoherent neutron scattering (EINS) to distinguish between protein dynamics and interfacial water dynamics. Powders of hydrogenated (hCPC) and perdeuterated (dCPC) CPC protein have been hydrated, respectively, with D(2)O or H(2)O and measured by EINS to separately probe protein dynamics (hCPC/D(2)O) and water dynamics (dCPC/H(2)O) at different time- and length-scales. We find that "fast" (<20 ps) local mean-square displacements (MSD) of both protein and interfacial water coincide all along the temperature range, with the same dynamical transition temperature at ~220 K. On higher resolution (<400 ps), two different types of motions can be separated: (i) localized motions with the same amplitude for CPC and hydration water and two transitions at ~170 and ~240 K for both; (ii) large scale fluctuations exhibiting for both water molecules and CPC protein a single transition at ~240 K, with a significantly higher amplitude for the interfacial water than for CPC. Moreover, by comparing these motions with bulk water MSD measured under the same conditions, we show no coupling between bulk water dynamics and protein dynamics all along the temperature range. These results show that interfacial water is the main "driving force" governing both local and large scale motions in proteins.     

Temperature dependence of solvation dynamics of probe molecules in methanol-doped ice and in liquid ethanol

We have studied the solvation statics and dynamics of coumarin 343 and a strong photoacid (pK* approximately 0.7) 2-naphthol-6, 8-disulfonate (2N68DS) in methanol-doped ice (1% molar concentration of methanol) and in cold liquid ethanol in the temperature range of 160-270 K. Both probe molecules show a relatively fast solvation dynamics in ice, ranging from a few tens of picoseconds at about 240 K to nanoseconds at about 160 K. At about 160 K in doped ice, we observe a sharp decrease of the dynamic Stokes shift of both coumarin 343 and 2N68DS. Its value is approximately only 200 cm-1 at approximately 160 K compared to about 1100 cm-1 at T >/= 200 K (at times longer than t > 10 ps). We find a good correlation between the inefficient and slow excited-state proton-transfer rate at low-temperature ice, T < 180 K, and the dramatic decrease of the solvation energy, as measured by the dynamic band shift, at these low temperatures. We find that the average solvation rate in ice is similar to its value in liquid ethanol at all given temperatures in the range of 200-250 K. The surprisingly fast solvation rate in ice is explained by the relatively large freedom of the water hydrogen rotation in ice Ih.     

NMR investigation of the porosity of hypercrosslinked polystyrene and the properties of water confined in its nanopores

The melting of water frozen preliminarily at 180 K in a free internal volume of water-swollen hypercrosslinked polystyrene networks with degrees of crosslinking ranging from 43 to 500% is studied by NMR. It is found that ice melts within a narrow range of low temperatures, 195?C225 K, demonstrating that the pores in the networks are small and uniform in size. It is, however, impossible to calculate the pore size via the Gibbs-Thomson equation, since the structure of water (and hence its properties) depend on a sample??s rate of freezing. We conclude that the activation of the orientational motions of water molecules starts to manifest itself already at 200 K; upon reaching 220?C230 K, the total mobility of water molecules is due largely to translational motions with an activation energy of 27?C28 kJ/mol.     

Observation of fragile-to-strong liquid transition in surface water in CeO2

A quasielastic neutron-scattering experiment carried out on a backscattering spectrometer with sub-micro eV resolution in the temperature range of 200-250 K has revealed the dynamics of surface water in cerium oxide on the time scale of hundreds of picoseconds. This slow dynamics is attributed to the translational mobility of the water molecules in contact with the surface hydroxyl groups. The relaxation function of this slow motion can be described by a slightly stretched exponential with the stretch factor exceeding 0.9, which indicates almost a Debye-type dynamics. Down to about 220 K, the temperature dependence of the residence time for water molecules follows a Vogel-Fulcher-Tamman law with the glass transition temperature of 181 K. At lower temperatures, the residence time behavior abruptly changes, indicating a fragile-to-strong liquid transition in surface water at about 215 K.     

Vibrational dynamics of ice in reverse micelles

The ultrafast vibrational dynamics of HDO:D(2)O ice at 180 K in anionic reverse micelles is studied by midinfrared femtosecond pump-probe spectroscopy. Solutions containing reverse micelles are cooled to low temperatures by a fast-freezing procedure. The heating dynamics of the micellar solutions is studied to characterize the micellar structure. Small reverse micelles with a water content up to approximately 150 water molecules contain an amorphous form of ice that shows remarkably different vibrational dynamics compared to bulk hexagonal ice. The micellar amorphous ice has a much longer vibrational lifetime than bulk hexagonal ice and micellar liquid water. The vibrational lifetime is observed to increase linearly from 0.7 to 4 ps with the resonance frequency ranging from 3100 to 3500 cm(-1). From the pump dependence of the vibrational relaxation the homogeneous linewidth of the amorphous ice is determined (55+/-5 cm(-1)).     

Multiple relaxation processes versus the fragile-to-strong transition in confined water

Broadband dielectric spectroscopy data on water confined in three different environments, namely at the surface of a globular protein or inside the small pores of two silica substrates, in the temperature range 140 K ≤ T ≤ 300 K, are presented and discussed in comparison with previous results from different techniques. It is found that all samples show a fast relaxation process, independently of the hydration level and confinement size. This relaxation is well known in the literature and its cross-over from Arrhenius to non-Arrhenius temperature behavior is the object of vivid debate, given its claimed relation to the existence of a second critical point of water. We find such a cross-over at a temperature of ~180 K, and assign the relaxation process to the layer of molecules adjacent and strongly interacting with the substrate surface. This is the water layer known to have the highest density and slowest translational dynamics compared to the average: its apparent cross-over may be due to the freezing of some degree of freedom and survival of very localized motions alone, to the onset of finite size effects, or to the presence of a calorimetric glass transition of the hydration shell at ~170 K. Another relaxation process is visible in water confined in the silica matrices: this is slower than the previous one and has distinct temperature behaviors, depending on the size of the confining volume and consequent ice nucleation.     

Ultrafast vibrational dynamics of interfacial water

We report investigations of the vibrational dynamics of water molecules at the water–air and at the water–lipid interface. Following vibrational excitation with an intense femtosecond infrared pulse resonant with the O–H stretch vibration of water, we follow the subsequent relaxation processes using the surface-specific spectroscopic technique of sum frequency generation. This allows us to selectively follow the vibrational relaxation of the approximately one monolayer of water molecules at the interface. Although the surface vibrational spectra of water at the interface with air and lipids are very similar, we find dramatic variations in both the rates and mechanisms of vibrational relaxation. For water at the water–air interface, very rapid exchange of vibrational energy occurs with water molecules in the bulk, and this intermolecular energy transfer process dominates the response. For membrane-bound water at the lipid interface, intermolecular energy transfer is suppressed, and intramolecular relaxation dominates. The difference in relaxation mechanism can be understood from differences in the local environments experienced by the interfacial water molecules in the two different systems.     

Host-assisted intramolecular vibrational relaxation at low temperatures: OH in an argon cage

The vibrational relaxation of hydroxyl radicals in the A (2)Sigma(+) (v=1) state has been studied using the semiclassical perturbation treatment at cryogenic temperatures. The radical is considered to be trapped in a closest packed cage composed of the 12 nearest argon atoms and undergoes local translation and hindered rotation around the cage center. The primary relaxation pathway is towards local translation, followed by energy transfer to rotation through hindered-to-free rotational transitions. Free-to-free rotational transitions are found to be unimportant. All pathways are accompanied by the propagation of energy to argon phonon modes. The deexcitation probability of OH(v=1) is 1.3 x 10(-7) and the rate constant is 4.7 x 10(5) s(-1) between 4 and 10 K. The negligible temperature dependence is attributed to the presence of intermolecular attraction (>kT) in the guest-host encounter, which counteracts the T(2) dependence resulting from local translation. Calculated relaxation time scales are much shorter than those of homonuclear molecules, suggesting the importance of the hindered and free motions of OH and strong guest-host interactions.     

The low frequency dynamics of supercooled LiBr, 6H2O

We present results of a series of experiments performed on LiBr, 6H(2)0 from room temperature down to 172 K ≈ 1.2T(g). These ultrasound, Brillouin and depolarized light scattering, and transient grating experiments show that, above 215 K, this solution behaves like supercooled water: its zero frequency sound velocity C(0) continuously decreases with decreasing temperature, and the reorientational dynamics of the water molecules can be directly detected at some temperatures of this domain. Conversely, below 215 K, a new regime sets in, where the apparent C(0) is practically temperature independent and where a β, Arrenhius like, relaxation process coexists with the usual, Vogel-Fulcher like, α relaxation process of the supercooled liquid. These results are similar to those recently obtained in LiCl, 6H(2)O. The onset of the new regime is possibly due to an increase of the interaction of the water molecules with a neighboring Li(+) ion when lowering the temperature. We also compare our results with published dielectric data on water solutions of glass forming polyalcohols. Some of them present a low temperature splitting of their relaxation time similar to what is found in LiBr, 6H(2)O.     

Thermal signature of hydrophobic hydration dynamics

Hydrophobic hydration, the perturbation of the aqueous solvent near an apolar solute or interface, is a fundamental ingredient in many chemical and biological processes. Both bulk water and aqueous solutions of apolar solutes behave anomalously at low temperatures for reasons that are not fully understood. Here, we use (2)H NMR relaxation to characterize the rotational dynamics in hydrophobic hydration shells over a wide temperature range, extending down to 243 K. We examine four partly hydrophobic solutes: the peptides N-acetyl-glycine-N'-methylamide and N-acetyl-leucine-N'-methylamide, and the osmolytes trimethylamine N-oxide and tetramethylurea. For all four solutes, we find that water rotates with lower activation energy in the hydration shell than in bulk water below 255 +/- 2 K. At still lower temperatures, water rotation is predicted to be faster in the shell than in bulk. We rationalize this behavior in terms of the geometric constraints imposed by the solute. These findings reverse the classical "iceberg" view of hydrophobic hydration by indicating that hydrophobic hydration water is less ice-like than bulk water. Our results also challenge the "structural temperature" concept. The two investigated osmolytes have opposite effects on protein stability but have virtually the same effect on water dynamics, suggesting that they do not act indirectly via solvent perturbations. The NMR-derived picture of hydrophobic hydration dynamics differs substantially from views emerging from recent quasielastic neutron scattering and pump-probe infrared spectroscopy studies of the same solutes. We discuss the possible reasons for these discrepancies.