Separable cooperative and localized translational motions of water confined by a chemically heterogeneous environment

We report quasi-elastic neutron scattering experiments at two resolutions that probe timescales of picoseconds to nanoseconds for the hydration dynamics of water, confined in a concentrated solution of N-acetyl-leucine-methylamide (NALMA) peptides in water over a temperature range of 248 K to 288 K. The two QENS resolutions used allow for a clean separation of two observable translational components, and ultimately two very different relaxation processes, that become evident when analyzed under a combination of the jump diffusion model and the relaxation cage model. The first translational motion is a localized beta-relaxation process of the bound surface water, and exhibits an Arrhenius temperature dependence and a large activation energy of approximately 8 kcal mol(-1). The second non-Arrhenius translational component is a dynamical signature of the alpha-relaxation of more fluid water, exhibiting a glass transition temperature of approximately 116 K when fit to the Volger Fulcher Tamman functional form. These peptide solutions provide a novel experimental system for examining confinement in order to understand the dynamical transition in bulk supercooled water by removing the unwanted interface of the confining material on water dynamics.     

High-resolution neutron-scattering study of slow dynamics of surface water molecules in zirconium oxide

We have performed a quasielastic neutron-scattering experiment on backscattering spectrometer with sub-mueV resolution to investigate the slow dynamics of surface water in zirconium oxide using the sample studied previously with a time-of-flight neutron spectrometer [E. Mamontov, J. Chem. Phys. 121, 9087 (2004)]. The backscattering measurements in the temperature range of 240-300 K have revealed a translational dynamics slower by another order of magnitude compared to the translational dynamics of the outer hydration layer observed in the time-of-flight experiment. The relaxation function of this slow motion is described by a stretched exponential with the stretch factors between 0.8 and 0.9, indicating a distribution of the relaxation times. The temperature dependence of the average residence time is non-Arrhenius, suggesting that the translational motion studied in this work is more complex than surface jump diffusion previously observed for the molecules of the outer hydration layer. The observed slow dynamics is ascribed to the molecules of the inner hydration layer that form more hydrogen bonds compared to the molecules of the outer hydration layer. Despite being slower by two orders of magnitude, the translational motion of the molecules of the inner hydration layer may have more in common with bulk water compared to the outer hydration layer, the dynamics of which is slower than that of bulk water by just one order of magnitude.     

Entropy and dynamics of water in hydration layers of a bilayer

We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme [two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (TΔS) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules.     

Dynamics of Water in NaxCoO2.yH2O

Incoherent inelastic neutron scattering experiments were performed on Na0.7CoO2 and Na0.28CoO2.1.3H2O in order to understand how the dynamics of the hydrogen-bond network of water is modified in the triangular crystalline lattice NaxCoO2.yH2O. Using quasi-elastic neutron scattering (QENS), we were able to differentiate between two types of proton dynamics: a fast process (due to water strongly bound into the sodium cobalt oxyhydrate structure during the hydration process) and a slow process (likely attributable to a collective motion). High-resolution QENS experiments, carried out on Na0.28CoO2.1.3H2O, show that, at temperatures above 310 K, the water dynamics can be well-described by a random jump diffusion model characterized by a diffusion constant equal to 0.9 x 10(-9)m2/s, which is significantly lower than the rate of diffusion for bulk water. Furthermore, our results indicate that, at room temperature, the sodium ions have no influence on the rotational dynamics of the "fast" water molecules.     

Monolayer of aerosol-OT surfactants adsorbed at the air/water interface: an atomistic computer simulation study

An atomistic molecular dynamics (MD) simulation has been carried out to investigate the structural and dynamical properties of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature and at a surface coverage corresponding to that at its critical micelle concentration (78 A(2)/molecule). The estimated thickness of the adsorbed layer is in good agreement with neutron reflection data. The study shows that the surfactants exhibit diffusive motion in the plane of the interface. It is observed that the surfactant monolayer has a strong influence in restricting both the translational and reorientational motions of the water molecules close to the interface. A drastic difference in the dipolar reorientational motion of water molecules in the aqueous layer is observed with a small variation of the distance from the surfactant headgroups. It has been observed that the water molecules in the first hydration layer (region 1) form strong hydrogen bonds with surfactant headgoups. This results in the slower structural relaxation of water-water hydrogen bonds in the first hydration layer compared to that in the pure bulk water. Most interestingly, we notice that the water molecules present in the layer immediately after the first hydration layer form weaker hydrogen bonds and thus relax faster than even pure bulk water.     

Hydrogen bond breaking mechanism and water reorientational dynamics in the hydration layer of lysozyme

The mechanism and the rate of hydrogen bond-breaking in the hydration layer surrounding an aqueous protein are important ingredients required to understand the various aspects of protein dynamics, its function, and stability. Here, we use computer simulation and a time correlation function technique to understand these aspects in the hydration layer of lysozyme. Water molecules in the layer are found to exhibit three distinct bond-breaking mechanisms. A large angle orientational jump of the donor water molecule is common among all of them. In the most common ( approximately 80%) bond-breaking event in the layer, the new acceptor water molecule comes from the first coordination shell (initially within 3.5 A of the donor), and the old acceptor water molecule remains within the first coordination shell, even after the bond-breaking. This is in contrast to that in bulk water, in which both of the acceptor molecules involve the second coordination shell. Additionally, the motion of the incoming and the outgoing acceptor molecules involved is not diffusive in the hydration layer, in contrast to their observed diffusive motion in the bulk. The difference in rotational dynamics between the bulk and the hydration layer water molecules is clearly manifested in the calculated time-dependent angular van Hove self-correlation function ( G(theta, t)) which has a pronounced two-peak structure in the layer, and this can be traced to the constrained translational motion in the layer. The longevity of the surrounding hydrogen bond network is found to be significantly enhanced near a hydrophilic residue.     

Influence of pyrogenic silica hydration on water diffusion in surface layers

The surface hydration of pyrogenic silica (aerosil) has been studied by slow neutron scattering. It was shown that the rate of diffusion changes with the thickness of the layer monotonically. The mean square displacement of water molecules from the equilibrium position in aerosil hydration shell is smaller than in bulk water, but twice larger than in ice.     

Water dynamics in graphite oxide investigated with neutron scattering

Graphite oxide is an inorganic multilayer system that preserves the layered structure of graphite but not the conjugated bond structure. In the past few years, detailed studies of the static structure of graphite oxide were carried out. This was mainly done by NMR investigations and led to a new structural model of graphite oxide. The layer distance of graphite oxide increases with increasing humidity level, giving rise to different spacings of the carbon layers in the range from 6 to 12 A. As a consequence, different types of motions of water and functional groups appear. Information about the mobility of the water molecules is not yet complete but is crucial for the understanding of the structure of the carbon layers as well as the intercalation process. In this paper, the hydration- and temperature-dependent dynamic behavior of graphite oxide will be investigated by quasielastic neutron scattering using the time-of-flight spectrometer NEAT at the Hahn-Meitner-Institut Berlin. The character of the embedded water does not change over a wide range of hydration levels. Especially the interlayer water remains tightly bound and does not show any translational motion. In samples with excess water, however, the water is also distributed in noninterlayer voids, leading to the observation of additional motions of bulklike or confined water. The dynamic behavior of hydrated graphite oxide can be described by a consistent model that combines two two-site jump motions for the motions of the water molecules and the motions of OH groups.     

Dynamics of trehalose molecules in confined solutions

The dynamics of trehalose molecules in aqueous solutions confined in silica gel have been studied by quasielastic neutron scattering (QENS). Small-angle neutron scattering measurements confirmed the absence of both sugar clustering and matrix deformation of the gels, indicating that the results obtained are representative of homogeneous trehalose solutions confined in a uniform matrix. The pore size in the gel is estimated to be 18 nm, comparable to the distances in cell membranes. For the QENS measurements, the gel was prepared from D2O in order to accentuate the scattering from the trehalose. Values for the translational diffusion constant and effective jump distance were derived from model fits to the scattering function. Comparison with QENS and NMR results in the literature for bulk trehalose shows that confinement on a length scale of 18 nm has no significant effect on the translational diffusion of trehalose molecules.     

Low temperature properties of acetonitrile confined in MCM-41

The effect of confinement on the phase changes and dynamics of acetonitrile in mesoporous MCM-41 was studied by use of adsorption, FT-IR, DSC, and quasi-elastic neutron scattering (QENS) measurements. Acetonitrile molecules in a monolayer interact strongly with surface hydroxyls to be registered and perturb the triple bond in the C[triple bond]N group. Adsorbed molecules above the monolayer through to the central part of the cylindrical pores are capillary condensed molecules (cc-acetonitrile), but they do not show the hysteresis loop in adsorption-desorption isotherms, i.e., second order capillary condensation. FT-IR measurements indicated that the condensed phase is very similar to the bulk liquid. The cc-acetonitrile freezes at temperatures that depend on the pore size of the MCM-41 down to 29.1 A (C14), below which it is not frozen. In addition, phase changes between alpha-type and beta-type acetonitriles were observed below the melting points. Application of the Gibbs-Thomson equation, assuming the unfrozen layer thickness to be 0.7 nm, gave the interface free energy differences between the interfaces, i.e., Deltagamma(l/alpha) = 22.4 mJ m(-2) for the liquid/pore surface (ps) and alpha-type/ps, and Deltagamma(alpha/beta) = 3.17 mJ m(-2) for alpha-type/ps and beta-type/ps, respectively. QENS experiments substantiate the differing behaviors of monolayer acetonitrile and cc-acetonitrile. The monolayer acetonitrile molecules are anchored so as not to translate. The two Lorentzian analysis of QENS spectra for cc-acetonitriles showed translational motion but markedly slowed. However, the activation energy for cc-acetonitrile in MCM-41 (C18) is 7.0 kJ mol(-1) compared to the bulk value of 12.7 kJ mol(-1). The relaxation times for tumbling rotational diffusion of cc-acetonitrile are similar to bulk values.     

Solid-state ligand dynamics in interpenetrating Mn[N(CN)(2)](2)(pyrazine): a neutron spectroscopy study

We have used quasielastic neutron scattering to probe the solid-state ligand dynamics in the coordination polymer Mn[N(CN)(2)](2)(pyz) [pyz = pyrazine] which has double-interpenetrating 3D lattices. A reversible structural phase transition occurs at 410 K as shown by neutron spectroscopy and differential scanning calorimetry. The origin of this transition is linked to rotational dynamics associated with the bridging pyz ligands. At 425 K, the pyrazine ring motion can be solely regarded as a 180 degrees reorientational jump about the axis defined by the Mn-N coordinative bonds, occurring with a correlation time of approximately 70 ps. This model can be extended to the 200-410 K temperature region using high-resolution backscattering spectroscopy to measure an identical motion on the time scale of nanoseconds with an activation energy of 24 +/- 2 kJ mol(-1). In contrast, no quasielastic scattering is seen for the 2D layered variant beta-Cu[N(CN)(2)](2)(pyz), owing to its more compact layer packing motif. Importantly, this work represents the very first study of solid-state rotational dynamics in an interpenetrating lattice structure.     

Water dynamics in silica nanopores: the self-intermediate scattering functions

The dynamics of water molecules confined in approximately cylindrical silica nanopores is investigated using molecular simulation. The model systems are pores of diameter varying between 20 and 40 ? containing water at room temperature and at full hydration, prepared using grand canonical Monte Carlo simulation. Water dynamics in these systems is studied via molecular dynamics simulation. The results of the basic characterization of these systems have been reported in A. A. Milischuk and B. M. Ladanyi [J. Chem. Phys. 135, 174709 (2011)]. The main focus of the present study is the self-intermediate scattering function (ISF), F(S)(Q, t), of water hydrogens, the observable in quasi-elastic neutron scattering experiments. We investigate how F(S)(Q, t) depends on the pore diameter, the direction and magnitude of the momentum transfer Q, and the proximity of water molecules to the silica surface. We also study the contributions to F(S)(Q, t) from rotational and translational motions of water molecules and the extent of rotation-translation coupling present in F(S)(Q, t). We find that F(S)(Q, t) depends strongly on the pore diameter and that this dependence is due mainly to the contributions to the ISF from water translational motion and can be attributed to the decreased mobility of water molecules near the silica surface. The relaxation rate depends on the direction of Q and is faster for Q in the axial than in the radial direction. As the magnitude of Q increases, this difference diminishes but does not disappear. We find that its source is mainly the anisotropy in translational diffusion at low Q and in molecular reorientation at higher Q values.     

Quasi-elastic neutron scattering study of dimethyl-sulfoxide-water mixtures: probing molecular mobility in a nonideal solution

The translational and rotational motions of water and dimethyl sulfoxide, [DMSO, (CH(3))(2)SO] have been investigated using quasi-elastic neutron scattering. Water-DMSO mixtures at five DMSO mole fractions, chi(DMSO), ranging from 0 to 0.75, were measured. Hydrogen-deuterium substitution was used to extract independently the water proton dynamics (d-DMSO-H(2)O), the DMSO methyl proton dynamics (h-DMSO-D(2)O) and to obtain background corrections (d-DMSO-D(2)O). The translational diffusion of water slows down significantly compared to bulk water at all chi(DMSO)>0. The rotational time constant for water exhibits a maximum at chi(DMSO)=0.33 that corresponds to the observed maximum of the viscosity of the mixture. Data for DMSO can be analyzed in terms of a relatively slow tumbling of the molecule about its center-of-mass in conjunction with random translational diffusion. The rotational time constant for this motion exhibits some dependence on chi(DMSO), while the translational diffusion constant shows no clear variation for chi(DMSO)>0. The results presented reinforce the idea that due to the stronger associative nature of DMSO, DMSO-water aggregates are formed over the whole composition range, disturbing the tetrahedral natural arrangement of the water molecules. As a consequence adding DMSO to water causes a drastic slowing down of the dynamics of the water molecule, and vice versa.     

Water hydrogen bond analysis on hydrophilic and hydrophobic biomolecule sites

Elastic and quasielastic neutron scattering experiments have been used to investigate the hydrogen bonding network dynamics of hydration water on hydrophilic and hydrophobic sites. To this end the evolution of hydration water dynamics of a prototypical hydrophobic amino acid with polar backbone, N-acetyl-leucine-methylamide (NALMA), and hydrophilic amino acid, N-acetyl-glycine-methylamide (NAGMA), has been investigated as a function of the molecular ratio water : peptide. The results suggest that the dynamical contribution of the intrinsic and low hydration molecules of water is characteristic of pure librational/rotational movement. The water molecule remains attached to the hydrophilic site with only the possibility of hindered rotations that eventually break the bond with the peptide and reform it immediately after. A gradual evolution from librational motions to hindered rotations is observed as a function of temperature. When the hydration increases, we observe (together with the hindered rotations of hydrogen bonds) a slow diffusion of water molecules on the surface of the peptides.     

High-density hydration layer of lysozymes: molecular dynamics decomposition of solution scattering data

A characterization of the physical properties of protein hydration water is critical for understanding protein structure and function. Recent small-angle X-ray and neutron scattering data indicate that the density of water on the surface of lysozyme is significantly higher than in bulk water. Here, we provide an interpretation of the scattering results using a molecular dynamics simulation, which allows us to make quantitative predictions about density variations in the first hydration shell. The perturbation relative to bulk water involves statistically significant changes in the average water structure in the first hydration layer. The water density in the first hydration shell is increased by 5% with respect to the bulk. In regions of higher water density, the water dipoles align more parallel to each other and the number of hydrogen bonds per water molecule is higher. Increased water density is found for water molecules interacting with hydrogen and carbon atoms in the backbone or with nonpolar or negatively charged side-chain groups.     

Structure and dynamics of a fluid phase bilayer on a solid support as observed by a molecular dynamics computer simulation

Simulations of a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine lipid bilayer interacting with a solid surface of hydroxylated nanoporous amorphous silica have been carried out over a range of lipid-solid substrate distances. The porous solid surface allowed the water layer to dynamically adjust its thickness, maintaining equal pressures above and below the membrane bilayer. Qualitative estimates of the force between the surfaces leads to an estimated lipid-silicon distance in very good agreement with the results of neutron scattering experiments. Detailed analysis of the simulation at the separation suggested by experiment shows that for this type of solid support the water layer between surfaces is very narrow, consisting only of bound waters hydrating the lipid head groups and hydrophilic silica surface. The reduced hydration, however, has only minor effects on the head group hydration, the orientation of water molecules at the interface, and the membrane dipole potential. Whereas these structural properties were not sensitive to the presence of the solid substrate, the calculated diffusion coefficient for translation of the lipid molecules was altered significantly by the silica surface.     

Experimental evidence of fragile-to-strong dynamic crossover in DNA hydration water

We used high-resolution quasielastic neutron scattering spectroscopy to study the single-particle dynamics of water molecules on the surface of hydrated DNA samples. Both H(2)O and D(2)O hydrated samples were measured. The contribution of scattering from DNA is subtracted out by taking the difference of the signals between the two samples. The measurement was made at a series of temperatures from 270 down to 185 K. The relaxing-cage model was used to analyze the quasielastic spectra. This allowed us to extract a Q-independent average translational relaxation time of water molecules as a function of temperature. We observe clear evidence of a fragile-to-strong dynamic crossover (FSC) at T(L)=222+/-2 K by plotting log versus T. The coincidence of the dynamic transition temperature T(c) of DNA, signaling the onset of anharmonic molecular motion, and the FSC temperature T(L) of the hydration water suggests that the change of mobility of the hydration water molecules across T(L) drives the dynamic transition in DNA.     

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.     

Adsorption dynamics of water on Pt{110}-(1 x 2)

The dynamics of H(2)O adsorption on Pt{110}-(1 x 2) is studied using supersonic molecular beam and temperature programed desorption techniques. The sticking probabilities are measured using the King and Wells method at a surface temperature of 165 K. The absolute initial sticking probability s(0) of H(2)O is 0.54+/-0.03 for an incident kinetic energy of 27 kJmol. However, an unusual molecular beam flux dependence on s(0) is also found. At low water coverage (theta<1), the sticking probability is independent of coverage due either to diffusion in an extrinsic precursor state formed above bilayer islands or to incorporation into the islands. We define theta=1 as the water coverage when the dissociative sticking probability of D(2) on a surface predosed with water has dropped to zero. The slow falling H(2)O sticking probability at theta>1 results from compression of the bilayer and the formation of multilayers. Temperature programed desorption of water shows fractional order kinetics consistent with hydrogen-bonded islands on the surface. A remarkable dependence of the initial sticking probability on the translational (1-27 kJ/mol) and internal energies of water is observed: s(0) is found to be essentially a step function of translational energy, increasing fivefold at a threshold energy of 5 kJ/mol. The threshold migrates to higher energies with increasing nozzle temperature (300-700 K). We conclude that both rotational state and rotational alignment of the water molecules in the seeded supersonic expansion are implicated in dictating the adsorption process.     

Probing surface properties and glass-liquid transition of amorphous solid water: temperature-programmed TOF-SIMS and TPD studies of adsorption/desorption of hexane

The interaction of hexane with amorphous solid water has been investigated in terms of the surface diffusion, hydrogen bond imperfections, hydrophobic hydration, crystallization, and glass-liquid transition. The hexane exhibits two main peaks in temperature-programmed desorption: one is ascribed to a complex formed at the surface or subsurface sites (135 K) and the other is caused by a bulk complex (165 K). The latter is associated with the presence of frozen-in imperfections in hydrogen bonds and formed provided that the annealing temperature of the film is below 130 K, whereas the former is created even when the film is annealed up to 150 K. Thus, the hexane-water interaction is hardly characterized by simple physisorption. The hexane is incorporated in the bulk during reorganization of hydrogen bonds due to rotational and translational diffusions of water molecules above 120-140 K, whereas the surface complex is formed even below 120 K due to the surface diffusion of molecules. The film undergoes abrupt dewetting at 165 K as a consequence of the glass-liquid transition. The slow evolution of the fluidity in the supercooled liquid phase may be responsible for the delay of the structural relaxation (165 K) relative to the onset of the translational molecular diffusion (135-140 K).