Local hydrogen bonding dynamics and collective reorganization in water: ultrafast infrared spectroscopy of HOD/D(2)O

We present an investigation into hydrogen bonding dynamics and kinetics in water using femtosecond infrared spectroscopy of the OH stretching vibration of HOD in D(2)O. Infrared vibrational echo peak shift and polarization-selective pump-probe experiments were performed with mid-IR pulses short enough to capture all relevant dynamical processes. The experiments are self-consistently analyzed with a nonlinear response function expressed in terms of three dynamical parameters for the OH stretching vibration: the frequency correlation function, the lifetime, and the second Legendre polynomial dipole reorientation correlation function. It also accounts for vibrational-relaxation-induced excitation of intermolecular motion that appears as heating. The long time, picosecond behavior is consistent with previous work, but new dynamics are revealed on the sub-200 fs time scale. The frequency correlation function is characterized by a 50 fs decay and 180 fs beat associated with underdamped intermolecular vibrations of hydrogen bonding partners prior to 1.4 ps exponential relaxation. The reorientational correlation function observes a 50 fs librational decay prior to 3 ps diffusive reorientation. Both of these correlation functions compare favorably with the predictions from classical molecular dynamics simulations. The time-dependent behavior can be separated into short and long time scales by the 340 fs correlation time for OH frequency shifts. The fast time scales arise from dynamics that are mainly local: fluctuations in hydrogen bond distances and angles within relatively fixed intermolecular configurations. On time scales longer than the correlation time, dephasing and reorientations reflect collective reorganization of the liquid structure. Since the OH transition frequency and dipole are only weakly sensitive to these collective coordinates, this is a kinetic regime which gives an effective rate for exchange of intermolecular structures.     

Water-methanol mixtures: topology of hydrogen bonded network

Molecular dynamics simulation has been performed to study the structure of water-methanol mixtures. Besides the evaluation of partial radial distribution functions describing the hydrogen-bonded structure of the mixtures with different composition, the statistical analysis of configurations was introduced resulting in a new insight in the clustering properties and topology of hydrogen-bonded network. The results have shown that mixtures of methanol and water exhibit extended structures in solution. At low methanol concentration the water molecules form a percolated network, the methanol molecules are incorporated as monomers or short chains and together form a percolated system. In methanol-rich mixtures short water chains and longer methanol chains build up the hydrogen-bonded clusters in the system. On the basis of the statistical analysis of configurations obtained from molecular dynamics simulation it has been found that more methanol molecules are connected to non-cyclic entities, while more water molecules form rings that might have been predicted on the basis of the stoichiometry of the mixtures. This finding can be explained by the presence of microscopic configurational inhomogeneity in water-methanol mixtures.     

Molecular dynamics study of triosephosphate isomerase from Trypanosoma cruzi in water/decane mixtures

A comprehensive study of the triosephosphate isomerase from the parasite Trypanosoma cruzi (TcTIM) in water, in decane, and in three water/decane mixtures was performed using molecular dynamics (MD) simulations in a time scale of 40 ns. The structure and dynamics of the enzyme, as well as the solvent molecules' distribution and mobility, were analyzed in detail. In the presence of decane, the amplitudes of the most important internal motions of the enzyme backbone were observed to depend on the solvent concentration: the higher the water concentration, the greater the amplitudes. Contrary to this trend, the amplitudes of the TcTIM motions in pure water were similar to those of the simulation with the lowest water concentration. The enzyme was observed to be almost motionless in pure decane due to a sharp increase of the number of intramolecular hydrogen bonds. This caused a contraction of the enzyme structure accompanied by a loss of secondary structure and of a decrease of the hydrophilic solvent accessible surface. A similar behavior, although to a lesser extent, was observed in the simulation at the lowest water concentration. Our results suggest that the presence of decane molecules located at specific sites of the enzyme might accelerate its internal movements, although a minimum number of water molecules is needed for the protein to keep its structure and dynamics. Altogether, this work provides new insight into protein and water behavior in organic solvents as well as into the dynamics of TcTIM itself.     

Nonideality in diffusion of ionic and neutral solutes and hydrogen bond dynamics in dimethyl sulfoxide-chloroform mixtures of varying composition

Molecular dynamics simulations of charged and neutral solutes in dimethyl sulfoxide (DMSO)-chloroform mixtures reveal pronounced nonideality in the solute diffusion with changes of composition of the mixtures. The diffusion coefficient of the anionic solute first decreases, passes through a minimum at DMSO mole fraction of about 0.50, and then increases to reach its value for pure DMSO. The diffusion coefficients of the cationic and neutral solutes are found to decrease with increase in DMSO content of the solvent mixture. The extent of nonideality in the diffusion and orientational relaxation of solvent molecules is found to be somewhat stronger than that in diffusion of the anionic solute in these mixtures. We have also calculated the relaxation of hydrogen bonds formed between DMSO and chloroform molecules. The lifetimes of DMSO-chloroform hydrogen bonds are found to increase monotonically with increase in DMSO concentration. The average number of hydrogen bonds and their average energies are also computed. It is found that an increase in DMSO concentration causes a decrease in the number of DMSO-chloroform hydrogen bonds per DMSO or chloroform molecules but increases the strength of these hydrogen bonds.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

Vibrational spectral diffusion and hydrogen bond dynamics in heavy water from first principles

We present a first-principles theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in heavy water without using any empirical model potentials. The calculations are based on ab initio molecular dynamics simulations for trajectory generation and a time series analysis using the wavelet method for frequency calculations. It is found that, in deuterated water, although a one-to-one relation does not exist between the instantaneous frequency of an OD bond and the distance of its associated hydrogen bond, such a relation does hold on average. The dynamics of spectral diffusion is investigated by means of frequency-time correlation and spectral hole dynamics calculations. Both of these functions are found to have a short-time decay with a time scale of approximately 100 fs corresponding to dynamics of intact hydrogen bonds and a slower long-time decay with a time constant of approximately 2 ps corresponding to lifetimes of hydrogen bonds. The connection of the slower time scale to the dynamics of local structural relaxation is also discussed. The dynamics of hydrogen bond making is shown to have a rather fast time scale of approximately 100 fs; hence, it can also contribute to the short-time dynamics of spectral diffusion. A damped oscillation is also found at around 150-200 fs, which is shown to have come from underdamped intermolecular vibrations of a hydrogen-bonded water pair. Such assignments are confirmed by independent calculations of power spectra of intermolecular motion and hydrogen bond kinetics using the population correlation function formalism. The details of the time constants of frequency correlations and spectral shifts are found to depend on the frequencies of chosen OD bonds and are analyzed in terms of the dynamics of hydrogen bonds of varying strengths and also of free non-hydrogen-bonded OD groups.     

Long-range influence of carbohydrates on the solvation dynamics of water--answers from terahertz absorption measurements and molecular modeling simulations

We present new terahertz (THz) spectroscopic measurements of solvated sugars and compare the effect of two disaccharides (trehalose and lactose) and one monosaccharide (glucose) with respect to the solute-induced changes in the sub-picosecond network dynamics of the hydration water. We found that the solute affects the fast collective network motions of the solvent, even beyond the first solvation layer. For all three carbohydrates, we find an increase of 2-4% in the THz absorption coefficient of the hydration water in comparison to bulk water. Concentration-dependent changes in the THz absorption between 2.1 and 2.8 THz of the solute-water mixture were measured with a precision better than 1% and were used to deduce a dynamical hydration shell, which extends from the surface up to 5.7 +/- 0.4 and 6.5 +/- 0.9 A for the disaccharides lactose and trehalose, respectively, and 3.7 +/- 0.9 A for the glucose. This exceeds the values for the static hydration shell as determined, for example, by scattering, where the long-range structure was found to be not significantly affected by the solute beyond the first hydration shell. When comparing all three carbohydrates, we found that the solute-induced change in the THz absorption depends on the product of molar concentration of the solute and the number of hydrogen bonds between the carbohydrate and water molecules. We can conclude that the long-range influence on the sub-picosecond collective water network motions of the hydration water is directly correlated with the average number of hydrogen bonds between the molecule and adjacent water molecules for carbohydrates. This implies that monosaccharides have a smaller influence on the surrounding water molecules than disaccharides. This could explain the bioprotection mechanism of sugar-water mixtures, which has been found to be more effective for disaccharides than for monosaccharides.     

Single particle and pair dynamics in water-formic acid mixtures containing ionic and neutral solutes: nonideality in dynamical properties

A series of molecular dynamics simulations of water-formic acid mixtures containing either an ionic solute or a neutral hydrophobic solute has been performed to study the extent of nonideality in the dynamics of these solutes for varying composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to show nonideal behavior with variation of composition, and similar nonideality is also observed for the diffusion and orientational relaxation of solvent molecules in these mixtures. The diffusion coefficient of a neutral hydrophobic solute, however, decreases monotonically with increase in water concentration. We have also investigated some of the pair dynamical properties such as water-water and water-formic acid hydrogen bond relaxation and residence dynamics of water molecules in water and formic acid hydration shells. The lifetimes of water-water hydrogen bonds are found to be longer than those between formic acid carbonyl oxygen-water hydrogen bonds, whereas the lifetimes of formic acid hydroxyl hydrogen-water hydrogen bonds are longer than those of water-water hydrogen bonds. In general, the hydrogen bond lifetimes for both water-water and water-formic acid hydrogen bonds are found to decrease with increase in water concentration. Residence times of water molecules also show the same trend with increase in formic acid concentration. Interestingly, these pair dynamical properties show a monotonic dependence on composition without any maximum or minimum and behave almost ideally with respect to changes in the composition of the mixtures. The present calculations are performed with fixed-charge nonpolarizable models of the solvent and solute molecules without taking into account many-body polarization effects in an explicit manner.     

Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations

A molecular dynamics simulation study is presented for the relaxation of the polarizability anisotropy in liquid mixtures of formamide and water, using a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors of the molecules, and the dipole-quadrupole polarizability of water species. The long time diffusive decay of the collective polarizability anisotropy correlations exhibits a substantial slowing down as the formamide mole fraction increases in the mixture. The diffusive times for the polarizability relaxation obtained from the authors' simulations are in good agreement with optical Kerr effect experimental data, and they are found to correlate nearly linearly with the estimated mean lifetimes of the hydrogen bonds within the mixture, suggesting that the relaxation of the hydrogen bond network is responsible to some extent for the collective relaxation of the polarizability anisotropy of the mixture. The short time behavior of the polarizability anisotropy relaxation was investigated by computing the nuclear response function, R(t), which is very rapidly dominated by the formamide contribution as it is added to water, due to the much larger polarizability anisotropy of formamide molecules compared to that of water. Several contributions to the Raman spectrum were also analyzed as a function of composition, and the dynamical origin of the different bands was determined.     

Molecular mechanisms of cryoprotection in aqueous proline: light scattering and molecular dynamics simulations

Free proline amino acid is a natural cryoprotectant expressed by numerous organisms under low-temperature stress. Previous reports have suggested that complex assemblies underlie its functional properties. We investigate here aqueous proline solutions as a function of temperature using combinations of Raman spectroscopy, Rayleigh-Brillouin light scattering, and molecular dynamics simulations with the view to revealing the molecular origins of the mixtures' functionality as a cryoprotectant. The evolution of the Brillouin frequency shifts and line widths with temperature shows that, above a critical proline concentration, the water-like dynamics is suppressed and viscoelastic behavior emerges: Here, the Landau-Placzek ratio also shows a temperature-independent maximum arising from concentration fluctuations. Molecular dynamics simulations reveal that the water-water correlations in the mixtures depend much more weakly on temperature than does bulk water. By contrast, the water OH Raman bands exhibit strong red-shifts on cooling similar to those seen in ices; however, no evidence of ice lattice phonons is observed in the low-frequency spectrum. We attribute this primarily to enhanced proline-water hydrogen bonding. In general, the picture that emerges is that aqueous proline is a heterogeneous mixture on molecular length scales (characterized by significant concentration fluctuations rather than well-defined aggregates). Simulations reveal that proline also appears to suppress the normal dependence of water structure on temperature and preserves the ambient-temperature correlations even in very cold solutions. The water structure in cold proline solutions therefore appears to be similar to that at a higher effective temperature. This, coupled with the emergence of glassy dynamics offers a molecular explanation for the functional properties of proline as a cryoprotectant without the need to invoke previously proposed complex aggregates.     

Microscale liquid dynamics and the effect on macroscale propagation in pillar arrays

Liquid dynamics in micropillar arrays have received significant fundamental interest and have offered opportunities for the development of advanced microfluidic, thermal management, and energy-harvesting devices. However, a comprehensive understanding of complex liquid behavior and the effect on macroscopic propagation rates in micropillar arrays is needed. In this work, we investigated the microscopic sweeping behavior of the liquid front along the spreading direction in micropillar arrays where the sweeping distance scales with the one-fifth power of time. We explain the scaling with a simplified model that captures the capillary pressure gradient at the liquid front. Furthermore, we show that such microscopic dynamics is the mechanism that decreases the macroscopic propagation rate. This effect is a result of the reduction in the interfacial energy difference used to generate the capillary pressure, which is explained with an energy-based model and corroborated with experiments. The results indicate the importance of accounting for the microscopic dynamics of the liquid on microstructured surfaces, particularly in sparse geometries.     

Molecular dynamics simulations of aqueous solutions of ethanolamines

We report on molecular dynamics simulations performed at constant temperature and pressure to study ethanolamines as pure components and in aqueous solutions. A new geometric integration algorithm that preserves the correct phase space volume is employed to study molecules having up to three ethanol chains. The most stable geometry, rotational barriers, and atomic charges were obtained by ab initio calculations in the gas phase. The calculated dipole moments agree well with available experimental data. The most stable conformation, due to intramolecular hydrogen bonding interactions, has a ringlike structure in one of the ethanol chains, leading to high molecular stability. All molecular dynamics simulations were performed in the liquid phase. The interaction parameters are the same for the atoms in the ethanol chains, reducing the number of variables in the potential model. Intermolecular hydrogen bonding is also analyzed, and it is shown that water associates at low water mole fractions. The force field reproduced (within 1%) the experimental liquid densities at different temperatures of pure components and aqueous solutions at 313 K. The excess and partial molar volumes are analyzed as a function of ethanolamine concentration.     

Experimental and kinetic study of shock initiated ignition in homogeneous methane–hydrogen–air mixtures at engine‐relevant conditions

The ignition delay time of two stoichiometric methane/hydrogen/air mixtures has been measured in a shock tube facility at pressures from 16 to 40 atm and temperatures from 1000 to 1300 K. Overall, the observed reduction in ignition delay with some methane replaced by hydrogen is relatively small given the large concentration of hydrogen involved in the current study. With a high hydrogen mole fraction (35% of the total fuel), a reduction of the ignition‐promoting effect was observed with reduced temperature. A detailed chemical kinetic mechanism was used to simulate ignitions of test mixtures behind reflected shocks. An analysis of the mechanism indicates that at higher temperatures, the rapid decomposition of hydrogen molecules leads to a quick formation of H radical pools, which promote the chain branching through H + O₂ ? O + OH. At lower temperatures, the branching efficiency of hydrogen is low; a weak effect of hydrogen on methane ignition could be result from the reaction between H₂ and methylperoxy CH₃O₂, which contributes extra H radicals to the reaction system. The effects of hydrogen also decrease with increasing pressure; this is related to the negative pressure dependence of hydrogen at the second ignition limit. © 2006 Wiley Periodicals, Inc. Int J Chem Kinet 38: 221–233, 2006     

Molecular dynamics investigation with the time resolved optical Kerr effect on the CS2-C6H6 mixtures

An investigation of the molecular dynamics in pure liquids and in mixtures through the technique of time resolved optical Kerr effect is performed. The samples studied were the mixtures of carbon disulfide (CS(2)) with benzene (C(6)H(6)). The molecular dynamics of the pure liquids is briefly discussed while the main results are obtained for the mixtures. A slow dynamics component is observed for the optical heterodyne detected optical Kerr effect transient decaying exponentially with time constants on picosecond time scale. The fast subpicosend time relaxations are analyzed in terms of the nondiffusive component of the spectral response that is associated with the molecular dynamics. The modifications of the spectrum are quantified, and the explanation of the observed changes is given in terms of the structural interaction configurations that produced changes in the intermolecular potential within which the molecules execute librational motions.     

Rotational dynamics of UVITEX-OB in alkanes, alcohols and binary mixtures

Steady-state and time resolved fluorescence anisotropy measurements were carried out to study the rotational diffusion dynamics of UVITEX-OB (U-OB) in series of alcohols, alkanes and binary mixtures of toluene and butanol at room temperature. The experimentally measured rotational reorientation times were compared with those estimated by the hydrodynamic and molecular models developed for microscopic friction. The experimental results are in good agreement with theoretical slip hydrodynamics and a deviation towards subslip behavior is noted. Also a faster rotation of the probe in binary mixture of toluene and butanol is noted as compared to that in alcohols and alkanes.     

Coupling between lysozyme and glycerol dynamics: microscopic insights from molecular-dynamics simulations

We explore possible molecular mechanisms behind the coupling of protein and solvent dynamics using atomistic molecular-dynamics simulations. For this purpose, we analyze the model protein lysozyme in glycerol, a well-known protein-preserving agent. We find that the dynamics of the hydrogen bond network between the solvent molecules in the first shell and the surface residues of the protein controls the structural relaxation (dynamics) of the whole protein. Specifically, we find a power-law relationship between the relaxation time of the aforementioned hydrogen bond network and the structural relaxation time of the protein obtained from the incoherent intermediate scattering function. We demonstrate that the relationship between the dynamics of the hydrogen bonds and the dynamics of the protein appears also in the dynamic transition temperature of the protein. A study of the dynamics of glycerol as a function of the distance from the surface of the protein indicates that the viscosity seen by the protein is not the one of the bulk solvent. The presence of the protein suppresses the dynamics of the surrounding solvent. This implies that the protein sees an effective viscosity higher than the one of the bulk solvent. We also found significant differences in the dynamics of surface and core residues of the protein. The former is found to follow the dynamics of the solvent more closely than the latter. These results allowed us to propose a molecular mechanism for the coupling of the solvent-protein dynamics.     

Communication: Unusual dynamics of hybrid nanoparticles and their binary mixtures

We present the results on the evolution of microscopic dynamics of hybrid nanoparticles and their binary mixtures as a function of temperature and wave vector. We find unexpectedly a nonmonotonic dependence of the structural relaxation time of the nanoparticles as a function of the morphology. In binary mixtures of two of the largest nanoparticles studied, we observe re-entrant vitrification as a function of the volume fraction of the smaller nanoparticle, which is unusual for such high diameter ratio. Possible explanation for the observed behavior is provided.     

Quasi-Hamiltonian equations of motion for internal coordinate molecular dynamics of polymers

Conventional molecular dynamics simulations of macromolecules require long computational times because the most interesting motions are very slow compared to the fast oscillations of bond lengths and bond angles that limit the integration time step. Simulation of dynamics in the space of internal coordinates, that is, with bond lengths, bond angles, and torsions as independent variables, gives a theoretical possibility of eliminating all uninteresting fast degrees of freedom from the system. This article presents a new method for internal coordinate molecular dynamics simulations of macromolecules. Equations of motion are derived that are applicable to branched chain molecules with any number of internal degrees of freedom. Equations use the canonical variables and they are much simpler than existing analogs. In the numerical tests the internal coordinate dynamics are compared with the traditional Cartesian coordinate molecular dynamics in simulations of a 56 residue globular protein. For the first time it was possible to compare the two alternative methods on identical molecular models in conventional quality tests. It is shown that the traditional and internal coordinate dynamics require the same time step size for the same accuracy and that in the standard geometry approximation of amino acids, that is, with fixed bond lengths, bond angles, and rigid aromatic groups, the characteristic step size is 4 fs, which is 2 times higher than with fixed bond lengths only. The step size can be increased up to 11 fs when rotation of hydrogen atoms is suppressed. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18 : 1354–1364, 1997     

Relaxation dynamics in glycerol-water mixtures. 2. Mesoscopic feature in water rich mixtures

The relaxation dynamics of water-rich glycerol-water mixtures is studied by broadband dielectric spectroscopy (BDS) at 173-323 K and differential scanning calorimetry (DSC) at 138-313 K. These data indicate the existence of the critical concentration of 40 mol % glycerol. In the studied temperature range for water-rich glycerol mixtures, the two states of water (ice and interfacial water) are observed in addition to water in the mesoscopic 40 mol % glycerol-water domains. The possible kinetics of water exchange between different water states is discussed in order to explain the mechanism of the broad melting behavior observed by DSC.