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.     

Diffusion of water in clays on the microscopic scale: modeling and experiment

Diffusion of water in montmorillonite clays at low hydration has been studied on the microscopic scale by two quasi-elastic neutron scattering techniques, neutron spin-echo (NSE) and time-of-flight (TOF), and by classical microscopic simulation. Experiment and simulation are compared both directly on the level of intermediate scattering functions, I(Q, t), and indirectly on the level of relaxation times after a model of atomic motion is applied. Regarding the dynamics of water in Na- and Cs-monohydrated montmorillonite samples, the simulation and NSE results show a very good agreement, both indicating diffusion coefficients of the order of (1-3) x 10(-10) m(2) s(-1). The TOF technique significantly underestimates water relaxation times (therefore overestimates water dynamics), by a factor of up to 3 and 7 in the two systems, respectively, primarily due to insufficiently long correlation times being probed. In the case of the Na-bihydrated system, the TOF results are in closer agreement with the other two techniques (the techniques differ by a factor of 2-3 at most), giving diffusion coefficients of (5-10) x 10(-10) m(2) s(-1). Attention has been also paid to the elastic incoherent structure factor, EISF(Q). Simulation has played a key role in understanding the various contributions to EISF(Q) in clay systems and in clearly distinguishing the signatures of "apparent" and true confinement. Indirectly, simulation highlights the difficulty in interpreting the EISF(Q) signal from powder clay samples used in experiments.     

Reduced mobility of di-propylene glycol methylether in its aqueous mixtures by quasielastic neutron scattering

The hydrogen (H-) bonding interplay between water and other organic molecules is important both in nature and in a wide range of technological applications. Structural relaxation and, thus, diffusion in aqueous mixtures are generally dependent on both the strength and the structure of the H-bonds. To investigate diffusion in H-bonding mixtures, we present a quasielastic neutron scattering study of di-propylene glycol methylether (2PGME) mixed with H(2)O (or D(2)O) over the concentration range 0-90 wt.% water. We observe a nonmonotonic behavior of the dynamics with a maximum in average relaxation time for the mixture with 30 wt.% water, which is more than a factor 2 larger compared to that of either of the pure constituents. This is a result in qualitative agreement with previous calorimetric studies and the behavior of aqueous mixtures of simple mono-alcohols. More surprisingly, we notice that the dynamics of the 2PGME molecules in the mixture is slowed down by more than a factor 3 at 30 wt.% water but that the water dynamics indicates an almost monotonous behavior. Furthermore, in the low momentum transfer (Q) range of the 2PGME, where the intermediate scattering function I(Q,t) is considerably stretched in time (i.e., the stretching parameter β ? 1), it is evident for the 2PGME-D(2)O samples that the Q-dependence of the inverse average relaxation time, <τ>(-1), is greater than 2. This implies that the relaxation dynamics is partly homogenously stretched, i.e., the relaxation of each relaxing unit is somewhat intrinsically stretched in time.     

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.     

Dynamics of propylene glycol and its oligomers confined to a single molecular layer

The dynamics of propylene glycol (PG) and its oligomers 7-PG and poly-propylene glycol (PPG), with M(w) = 4000 (approximately 70 monomers), confined in a Na-vermiculite clay have been investigated by quasielastic neutron scattering. The liquids are confined to single molecular layers between clay platelets, giving a true two-dimensional liquid. Data from three different spectrometers of different resolutions were Fourier transformed to S(Q,t) and combined to give an extended dynamical time range of 0.3-2000 ps. An attempt was made to distinguish the diffusive motion from the methyl group rotation and a fast local motion of hydrogen in the polymer backbone. The results show that the average relaxation time tau(d) of this diffusive process is, as expected, larger than the relaxation time tau averaged over all dynamical processes observed in the experimental time window. More interesting, it is evident that the severe confinement has a relatively small effect on tau(d) at T = 300 K, this holds particularly for the longest oligomer, PPG. The most significant difference is that the chain-length dependence of tau(d) is weaker for the confined liquids, although the slowing down in bulk PG due to the formation of a three-dimensional network of OH-bonded end groups reduces this difference. The estimated average relaxation time tau at Q = 0.92 Angstroms(-1) for all the observed processes is in excellent agreement with the previously reported dielectric alpha relaxation time in the studied temperature range of 260-380 K. The average relaxation time tau (as well as the dielectric alpha relaxation time) is also almost unaffected by the confinement to a single molecular layer, suggesting that the interaction with the clay surfaces is weak and that the reduced dimensionality has only a weak influence on the time scale of all the dynamical processes observed in this study.     

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

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

Plasticizer effect on the dynamics of polyvinylchloride studied by dielectric spectroscopy and quasielastic neutron scattering

We have studied the influence of plasticization on the microscopic dynamics of a glass-forming polymer. For this purpose we studied polyvinylchloride (PVC) with and without the commercially used plasticizer dioctylphthalate (DOP). We used dielectric spectroscopy and inelastic neutron scattering employing the neutron spin echo (NSE) technique. For both kinds of spectra the alpha relaxation could be consistently described by a model involving a distribution of individual relaxations of the Kohlrausch type. In contrast to earlier studies it turned out that an asymmetric distribution is necessary to fit the data at the lower temperatures investigated here. The shape parameters of the distribution (width, skewness) for PVC and PVC/DOP turned out to coincide when the characteristic relaxation times were the same. This means that the plasticizer only induces a remapping of the temperature dependence of the alpha relaxation. Comparison of NSE spectra S(Q,t)S(Q) at different scattering vectors Q gave the result that the slowing down at the structure factor peak Q(max) is surprisingly small for PVC while it is in the normal range for PVC/DOP.     

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.     

Quasielastic and inelastic neutron scattering on hydrated calcium silicate pastes

Using the inverse geometry spectrometer QENS at the Intense Pulsed Neutron Source of the Argonne National Laboratory, we collected quasielastic and inelastic neutron scattering spectra of hydrated tricalcium and dicalcium silicate, the main components of ordinary Portland cement. Data were obtained at different curing time, from a few hours to several months. Both the quasielastic and inelastic spectra have been analyzed at the same time according to the relaxing cage model, which is a model developed to describe the dynamics of water at supercooled temperatures. Short-time and long-time dynamics of hydration water in hydrated cement pastes as a function of the curing time have been simultaneously obtained. The results confirm the findings reported in previous experiments showing that it is possible to fit consistently the quasielastic and inelastic spectra giving insights on the effect of the curing time on the short-time vibrational dynamics of hydration water.     

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.     

Neutron scattering investigation of a diluted blend of poly(ethylene oxide) in polyethersulfone

By using quasielastic neutron scattering (QENS) with isotopic labeling we have investigated the component dynamics in a miscible blend of polyethersulfone (PES) and poly(ethylene oxide) (PEO) with 75% content in weight of PES. Due to the large difference in the glass-transition temperatures, T(g)'s, of the two polymers (T(g) (PEO) approximately equal to 220 K, T(g) (PES) approximately equal to 382 K) the dynamic asymmetry in the system dramatically increases when approaching the average T(g) of the blend, . For the fast (PEO) component, this leads to a behavior which hints a crossover from typical glass-forming liquidlike dynamics at high temperatures to confined dynamics close to induced by the freezing of the segmental motions of the slow PES. The features of the confined PEO motion observed by QENS are similar to those of the secondary gamma-relaxation detected for pure (semicrystalline) PEO. A neutron diffraction study of the short-range order of the homopolymers and the blend suggests that this coincidence could be due to similarities in the intermolecular packing of PEO and PES polymers.     

Dynamics of small-molecule glass formers confined in nanopores

We report a comparative neutron scattering study of the molecular mobility and nonexponential relaxation of three structurally similar glass-forming liquids, isopropanol, propylene glycol, and glycerol, both in bulk and confined in porous Vycor glass. Confinement reduces molecular mobility in all three liquids, and suppresses crystallization in isopropanol. High-resolution quasielastic neutron scattering spectra were fit to Fourier transformed Kohlrausch functions exp[-(t∕τ)(β)], describing the α-relaxation processes in these liquids. The stretching parameter β is roughly constant with wavevector Q and over the temperature range explored in bulk glycerol and propylene glycol, but varies both with Q and temperature in confinement. Average relaxation times <τ(Q)> are longer at lower temperatures and in confinement. They obey a power law <τ(Q)> ∝ Q(-γ), where the exponent γ is modified by confinement. Comparison of the bulk and confined liquids lends support to the idea that structural and∕or dynamical heterogeneity underlies the nonexponential relaxation of glass formers, as widely hypothesized in the literature.     

Confined polymer dynamics on clay platelets

The structure and dynamics of poly(ethylene oxide) adsorbed on dispersed clay platelets are investigated by small-angle neutron scattering and neutron spin-echo spectroscopy. The intermediate scattering function has a mobile contribution described by the Zimm theory and an immobile contribution that is constant within the time window. The immobile fraction as a function of the scattering vector Q is described by a Lorentz function, from which a localization length is determined. The relaxation rates grow with polymer concentration in agreement with dielectric measurements but contrary to pure polymer gels.     

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.     

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.     

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.     

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.     

Dynamic behavior of water within a polymer electrolyte fuel cell membrane at low hydration levels

Protonic conduction across the membrane of a polymer electrolyte fuel cell is intimately related to the dynamic behavior of water present within the membrane. To further the understanding of water dynamics in these materials, quasielastic neutron scattering (QENS) has been used to investigate the picosecond dynamic behavior of water within a perfluorosulfonated ionomer (PFSI) membrane under increasing hydration levels from dry to saturation. Evaluation of the elastic incoherent structure factor (EISF) reveals an increase in the characteristic length-scale of confinement as the number of water molecules in the membrane increases, tending to an asymptotic value at saturation. The fraction of elastic incoherent scattering observed at high Q over all hydration levels is well fit by a simple model that assumes a single, nondiffusing hydronium ion per membrane sulfonic acid site. The quasielastic component of the fitted data indicates confined dynamic behavior for scattering vectors less than 0.7 A(-1). As such, the dynamic behavior was interpreted using continuous diffusion confined within a sphere at Q < 0.7 A(-1) and random unconstrained jump diffusion at Q > 0.7 A(-1). As the number of water molecules in the membrane increases, the characteristic residence times obtained from both models is reduced. The increased dynamical frequency is further reflected in the diffusion coefficients predicted by both models. Between low hydration (2 H2O/SO3H) and saturation (16 H2O/SO3H), the continuous spherical diffusion coefficient changes from 0.46 +/- 0.12 to 1.04 +/- 0.12 (10(-5) cm2/s) and jump diffusion indicates an increase from 1.21 +/- 0.03 to 2.14 +/- 0.08 (10(-5) cm2/s). Overall, the dynamic behavior of water has been quantified over different length scale regimes, the results of which may be rationalized on the basis of the formation of water clusters in the hydrophilic domain that expand toward an asymptotic upper limit with increased hydration.     

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.