Structural dynamics of supercooled water from quasielastic neutron scattering and molecular simulations

One of the outstanding challenges presented by liquid water is to understand how molecules can move on a picosecond time scale despite being incorporated in a three-dimensional network of relatively strong H-bonds. This challenge is exacerbated in the supercooled state, where the dramatic slowing down of structural dynamics is reminiscent of the, equally poorly understood, generic behavior of liquids near the glass transition temperature. By probing single-molecule dynamics on a wide range of time and length scales, quasielastic neutron scattering (QENS) can potentially reveal the mechanistic details of water's structural dynamics, but because of interpretational ambiguities this potential has not been fully realized. To resolve these issues, we present here an extensive set of high-quality QENS data from water in the range 253-293 K and a corresponding set of molecular dynamics (MD) simulations to facilitate and validate the interpretation. Using a model-free approach, we analyze the QENS data in terms of two motional components. Based on the dynamical clustering observed in MD trajectories, we identify these components with two distinct types of structural dynamics: picosecond local (L) structural fluctuations within dynamical basins and slower interbasin jumps (J). The Q-dependence of the dominant QENS component, associated with J dynamics, can be quantitatively rationalized with a continuous-time random walk (CTRW) model with an apparent jump length that depends on low-order moments of the jump length and waiting time distributions. Using a simple coarse-graining algorithm to quantitatively identify dynamical basins, we map the newtonian MD trajectory on a CTRW trajectory, from which the jump length and waiting time distributions are computed. The jump length distribution is gaussian and the rms jump length increases from 1.5 to 1.9 A? as the temperature increases from 253 to 293 K. The rms basin radius increases from 0.71 to 0.75 A? over the same range. The waiting time distribution is exponential at all investigated temperatures, ruling out significant dynamical heterogeneity. However, a simulation at 238 K reveals a small but significant dynamical heterogeneity. The macroscopic diffusion coefficient deduced from the QENS data agrees quantitatively with NMR and tracer results. We compare our QENS analysis with existing approaches, arguing that the apparent dynamical heterogeneity implied by stretched exponential fitting functions results from the failure to distinguish intrabasin (L) from interbasin (J) structural dynamics. We propose that the apparent dynamical singularity at ～220 K corresponds to freezing out of J dynamics, while the calorimetric glass transition corresponds to freezing out of L dynamics.     

Dynamics of water sorbed in reverse osmosis polyamide membrane

Dynamical motion of water sorbed in reverse osmosis polyamide membrane (ROPM) material is reported as studied by quasielastic neutron scattering (QENS) technique. The ROPM studied here has pore size of 4.4 Å as determined by positron annihilation lifetime spectroscopy. Analysis of the QENS data showed that diffusion behavior of the water within the membrane is describable by random jump diffusion model. A much longer residence time is found as compared to bulk water. Positive shift of the freezing point as observed in differential scanning calorimetry indicates presence of strong attractive interaction corroborating the slower diffusivity as observed in QENS.     

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.     

Dynamics of a triphenylene discotic molecule,HAT6, in the columnar and isotropic liquid phases

Discotic molecules have planar, disklike polyaromatic cores that can self-assemble into "molecular wires". Highly anisotropic charge transfer along the wires arises when there is sufficient intermolecular overlap of the pi-orbitals of the molecular cores. Discotic materials can be applied in molecular electronics, field-effect transistors, and-recently with record quantum efficiencies-photovoltaics (Schmidt-Mende, L.; Fechtenk?tter, A.; Müllen, K.; Moons, E.; Frien, R. H.; MacKenzie, J. D. Science 2001, 293, 1119). A combination of quasielastic neutron scattering (QENS) measurements with molecular dynamics simulations on the discotic molecule hexakis(n-hexyloxy)triphenylene (HAT6) shows that the dynamics of the cores and tails of discotic molecules are strongly correlated. Core and tail dynamics are not separated, the system being characterized by overall in-plane motion, on a time scale of 0.2 ps, and softer out-of-plane motions at 7 ps. Because charge transfer between the molecules is on similar time scales, these motions are relevant for the conducting properties of the materials. Both types of motion are dominated by van der Waals interactions. Small-amplitude in-plane motions in which the disks move over each other are almost entirely determined by tail/tail interactions, these also playing an important role in the out-of-plane motion. The QENS measurements reveal that these motions are little changed by passing from the columnar phase to the isotropic liquid phase, just above the clearing temperature. The model of four HAT6 molecules in a column reproduces the measured QENS spectrum of the liquid phase, suggesting that correlations persist within the liquid phase over about this number of disks.     

Equation of state of water under negative pressure

We report on the simultaneous measurements of the speed of sound and the density in liquid water under negative pressure. Application of a focused acoustic wave to the bulk liquid is able to generate negative pressures before nucleation of the vapor phase occurs. A method for time-resolved Brillouin scattering is developed to measure the speed of sound during the passage of a 1 MHz ultrasonic wave. This is coupled with a fiber optic probe hydrophone which allows the determination of the density. Together, these methods give an ambient temperature equation of state of metastable liquid water down to the acoustic cavitation threshold. Empirical equations of state of water are based on experimental data at positive pressure; the validity of their extrapolation to negative pressures had been tested only indirectly or with very weakly metastable liquid. We provide thermodynamic data that prove the fidelity of recent equations of state down to -26 MPa. However, this raises questions regarding the nature of the cavitation threshold observed in acoustic experiments, which is far less negative than expected.     

Interlayer water molecules in vanadium pentoxide hydrate. 8. Dynamic properties by quasi-elastic neutron scattering

The dynamics of water molecules in the layered vanadium pentoxide hydrate, V(2)O(5).nH(2)O, were studied by quasi-elastic neutron scattering (QENS) measurements. Heterogeneity of the dynamic properties was confirmed by alpha-relaxation model analysis. Translational diffusion of monolayer and double-layer water molecules is by site-to-site diffusion and is reduced relative to that of bulk water. Water molecules lose their mobility markedly and solidify with decreasing temperature. However, mobile water remains at 253 K. Rotational diffusion coefficients are unaffected by confinement and are very similar to the bulk values determined at temperatures in the range 253-298 K. The dynamic speed characterized by QENS is much faster than that expected from the data determined by deuterium NMR (DNMR) measurements at low temperatures.     

Influence of isotherm inflection on the loading dependence of the diffusivities of n-hexane and n-heptane in MFI zeolite. Quasi-elastic neutron scattering experiments supplemented by molecular simulations

Quasi-Elastic Neutron Scattering (QENS) experiments were carried out to determine (a) Fick diffusivity, D (b) self-diffusivity, Dself, and (c) 1/Gamma, the inverse of the thermodynamic correction factor, for n-hexane (nC6) and n-heptane (nC7) in MFI zeolite (all silica silicalite-1) at 300 K for a variety of loadings. These experimental results are compared with configurational-bias Monte Carlo (CBMC) and molecular dynamics (MD) simulations of, respectively, the adsorption isotherms and diffusivities. For n-hexane, the CBMC simulated isotherm shows a slight inflection at a loading=4 molecules per unit cell; this inflection manifests, also, in the loading dependence of 1/Gamma, obtained from QENS. The trend in the loading dependence of the Fick D and Dself of nC6 obtained from QENS matches the MD simulation results. For nC7 the CBMC simulated isotherm shows a strong inflection at a loading=4 molecules per unit cell. At this loading=4, 1/Gamma tends to zero and there is a very good match between QENS and molecular simulations for the loading dependence of 1/Gamma. Both MD simulations and QENS data on the Fick diffusivity shows a sharp maximum at a loading in the region of=4. For both nC6 and nC7 the simulated values of diffusivity are about an order of magnitude higher than those determined from QENS.     

Dynamics of water in a molecular sieve by quasielastic neutron scattering

We have investigated the dynamics of water confined in a molecular sieve, with a cylindrical pore diameter of 10 A, by means of quasielastic neutron scattering (QENS). Both the incoherent and coherent intermediate scattering functions I(Q,t) were determined by time-of-flight QENS and the neutron spin-echo technique, respectively. The results show that I(Q,t) is considerably more stretched in time with a slightly larger average relaxation time in the case of coherent scattering. From the Q dependence of I(Q,t) it is clear that the observed dynamics is almost of an ordinary translational nature. A comparison with previous dielectric measurements suggests a possible merging of the alpha and beta relaxations of the confined water at T=185 K, although the alpha relaxation cannot be directly observed at lower temperatures due to the severe confinement. The present results are discussed in relation to previous results for water confined in a Na-vermiculite clay, where the average relaxation time from spin-echo measurements was found to be slower than in the present system (particularly at low temperatures).     

Water motion in reverse micelles studied by quasielastic neutron scattering and molecular dynamics simulations

Motion of water molecules in Aerosol OT [sodium bis(2-ethylhexyl) sulfosuccinate, AOT] reverse micelles with water content w(0) ranging from 1 to 5 has been explored both experimentally through quasielastic neutron scattering (QENS) and with molecular dynamics (MD) simulations. The experiments were performed at the energy resolution of 85 microeV over the momentum transfer (Q) range of 0.36-2.53 A(-1) on samples in which the nonpolar phase (isooctane) and the AOT alkyl chains were deuterated, thereby suppressing their contribution to the QENS signal. QENS results were analyzed via a jump-diffusion/isotropic rotation model, which fits the results reasonably well despite the fact that confinement effects are not explicitly taken into account. This analysis indicates that in reverse micelles with low-water content (w(0)=1 and 2.5) translational diffusion rate is too slow to be detected, while for w(0)=5 the diffusion coefficient is much smaller than for bulk water. Rotational diffusion coefficients obtained from this analysis increase with w(0) and are smaller than for bulk water, but rotational mobility is less drastically reduced than translational mobility. Using the Faeder/Ladanyi model [J. Phys. Chem. B 104, 1033 (2000)] of reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function F(S)(Q,t) for water hydrogens. Comparison of the time Fourier transform of this F(S)(Q,t) with the QENS dynamic structure factor S(Q,omega), shows good agreement between the model and experiment. Separate intermediate scattering functions F(S) (R)(Q,t) and F(S) (CM)(Q,t) were determined for rotational and translational motion. Consistent with the decoupling approximation used in the analysis of QENS data, the product of F(S) (R)(Q,t) and F(S) (CM)(Q,t) is a good approximation to the total F(S)(Q,t). We find that the decay of F(S) (CM)(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior is due to lower water mobility close to the interface and to confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay. However, rotational mobility of O-H bond vectors in the interfacial region remains fairly high due to the lower density of water-water hydrogen bonds in the vicinity of the interface.     

正丙醇水溶液准弹性中子散射光谱的分子动力学模拟

准弹性中子散射(QENS)光谱是获取溶液分子动力学性质的重要方法,但光谱解析模型的有效性和去耦合近似的合理性仍存在争议.本文利用分子动力学模拟方法获取纯水和正丙醇水溶液中羟基氢原子的自相关中间散射函数FS(Q,t)和去耦合近似函数FP(Q,t),以及相关性质来评价它们的合理性.结果表明,在低动量转移范围内平-转去耦合近似相对合理,水分子的平-转耦合贡献较小,混合溶液中水分子的平-转耦合项和转动项随动量转移Q值增大而增大,二者显现相互抵消趋势.对于混合溶液中的正丙醇羟基氢原子,由于FS(Q,t)和质心自相关中间散射函数FCM(Q,t)偏差较大,利用实验光谱直接拟合分子平动扩散系数是不合适的.三种平动模型获取的纯水和正丙醇水溶液分子平动扩散系数与实验结果一致,略高于Einstein均方位移方法所得结果.水分子在纯水和混合溶液中表现为跳跃转动,而不是连续转动.正丙醇分子存在转动各向异性,羟基氢原子沿羟基向量为跳跃转动,沿相对质心向量可近似为连续转动.模拟结果显示,高动量转移范围平-转耦合项贡献较大,直接拟合实验光谱获取分子转动扩散系数或弛豫时间是不合适的.鉴于低动量转移范围内转动和平转耦合贡献较小,以及二者的抵消作用,在此范围内获取水分子平动信息是现实可行的.     

Quasielastic neutron scattering study of hydrogen motions in an aqueous poly(vinyl methyl ether) solution

We present a quasielastic neutron scattering (QENS) investigation of the component dynamics in an aqueous Poly(vinyl methyl ether) (PVME) solution (30% water content in weight). In the glassy state, an important shift in the Boson peak of PVME is found upon hydration. At higher temperatures, the diffusive-like motions of the components take place with very different characteristic times, revealing a strong dynamic asymmetry that increases with decreasing T. For both components, we observe stretching of the scattering functions with respect to those in the bulk and non-Gaussian behavior in the whole momentum transfer range investigated. To explain these observations we invoke a distribution of mobilities for both components, probably originated from structural heterogeneities. The diffusive-like motion of PVME in solution takes place faster and apparently in a more continuous way than in bulk. We find that the T-dependence of the characteristic relaxation time of water changes at T ? 225 K, near the temperature where a crossover from a low temperature Arrhenius to a high temperature cooperative behavior has been observed by broadband dielectric spectroscopy (BDS) [S. Cerveny, J. Colmenero and A. Alegri?a, Macromolecules, 38, 7056 (2005)]. This observation might be a signature of the onset of confined dynamics of water due to the freezing of the PVME dynamics, that has been selectively followed by these QENS experiments. On the other hand, revisiting the BDS results on this system we could identify an additional "fast" process that can be attributed to water motions coupled with PVME local relaxations that could strongly affect the QENS results. Both kinds of interpretations, confinement effects due to the increasing dynamic asymmetry and influence of localized motions, could provide alternative scenarios to the invoked "strong-to-fragile" transition.     

Inverse freezing in α-cyclodextrin solutions probed by quasi elastic neutron scattering

Solutions of α-cyclodextrin, 4-methylpyridine and water undergo a reversible liquid–solid phase transition on heating (“inverse freezing”). In this paper quasi elastic neutron scattering (QENS) measurements are reported, from which the diffusive dynamics of different components in the liquid and solid phases is determined. Results imply that, in solution, α-cyclodextrin is contained in a solvation shell of 4-methylpyridine molecules, while in the crystal phase, the majority of the cell contents are immobilized on the nanosecond timescale. This information will be important in evaluating the entropy of the system in its different phases and in understanding the unusual phase transition.     

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.     

Single particle and collective hydration dynamics for hydrophobic and hydrophilic peptides

We have conducted extensive molecular dynamics simulations to study the single particle and collective dynamics of water in solutions of N-acetyl-glycine-methylamide, a model hydrophilic protein backbone, and N-acetyl-leucine-methylamide, a model (amphiphilic) hydrophobic peptide, as a function of peptide concentration. Various analytical models commonly used in the analysis of incoherent quasielastic neutron scattering (QENS), are tested against the translational and rotational intermediate scattering function, the mean square displacement of the water molecule center of mass, and fits to the second-order rotational correlation function of water evaluated directly from the simulation data. We find that while the agreement between the model-free analysis and analytical QENS models is quantitatively poor, the qualitative feature of dynamical heterogeneity due to caging is captured well by all approaches. The center of mass collective and single particle intermediate scattering functions of water calculated for these peptide solutions show that the crossover from collective to single particle-dominated motions occurs at a higher value of Q for high concentration solutions relative to low concentration because of the greater restriction in movement of water molecules due to confinement. Finally, we have shown that at the same level of confinement of the two peptides, the aqueous amphiphilic amino acid solution shows the strongest deviation between single particle and collective dynamics relative to the hydrophilic amino acid, indicating that chemical heterogeneity induces even greater spatial heterogeneity in the water dynamics.     

Interlayer water molecules in vanadium pentoxide hydrate. IX. Anisotropic translational diffusion leading to anisotropic ac conductivity

The anisotropy of the dynamic properties of interlayer water molecules along the a and b axes of vanadium pentoxide hydrate, orthorhombic V2O5.nH2O, was studied using quasielastic neutron scattering (QENS) in relation to the anisotropy of the ac conductivity. The QENS spectra were analyzed using a stretched exponential function and a Lorentzian function. Both methods showed that the double-layer water molecules along the b axis are more mobile than those along the a axis. The difference in mobility between the two axes is more pronounced using a Lorentzian function analysis. These facts suggest that the diffusion coefficient of water molecules along the b axis is larger than that along the a axis, which is closely related to the ac conductivity originating from proton hopping. The anisotropy of the dynamic motion of water molecules can be attributed to the shorter b-axis length (b=3.60 A), with respect to the longer and less regular repetition of the atomic arrangements along the a axis (42.34 A).     

Molecular view of water dynamics near model peptides

Incoherent quasi-elastic neutron scattering (QENS) has been used to measure the dynamics of water molecules in solutions of a model protein backbone, N-acetyl-glycine-methylamide (NAGMA), as a function of concentration, for comparison with results for water dynamics in aqueous solutions of the N-acetyl-leucine-methylamide (NALMA) hydrophobic peptide at comparable concentrations. From the analysis of the elastic incoherent structure factor, we find significant fractions of elastic intensity at high and low concentrations for both solutes, which corresponds to a greater population of protons with rotational time scales outside the experimental resolution (>13 ps). The higher-concentration solutions show a component of the elastic fraction that we propose is due to water motions that are strongly coupled to the solute motions, while for low-concentration solutions an additional component is activated due to dynamic coupling between inner and outer hydration layers. An important difference between the solute types at the highest concentration studied is found from stretched exponential fits to their experimental intermediate scattering functions, showing more pronounced anomalous diffusion signatures for NALMA, including a smaller stretched exponent beta and a longer structural relaxation time tau than those found for NAGMA. The more normal water diffusion exhibited near the hydrophilic NAGMA provides experimental support for an explanation of the origin of the anomalous diffusion behavior of NALMA as arising from frustrated interactions between water molecules when a chemical interface is formed upon addition of a hydrophobic side chain, inducing spatial heterogeneity in the hydration dynamics in the two types of regions of the NALMA peptide. We place our QENS measurements on model biological solutes in the context of other spectroscopic techniques and provide both confirming as well as complementary dynamic information that attempts to give a unifying molecular view of hydration dynamics signatures near peptides and proteins.     

Quasielastic neutron scattering in soft matter

The general trend in soft matter is to study systems of increasing complexity which are more technologically and biologically relevant. This is facilitated by the capability of quasielastic neutron scattering (QENS) to selectively probe spatially resolved dynamical modes at a molecular level. The large number of recent publications using QENS for investigating complex and multi-component soft matter systems, serves as recognition of the suitability of this technique by the scientific community. Exploiting its complementarity with molecular dynamics (MD) simulations and other experimental techniques is the basis of a successful methodology for this scientific challenge. We illustrate the potential of QENS with three kinds of soft materials whose structural units increase in size/complexity: lipids, polymers and biomolecules.