The effect of the counterion on water mobility in reverse micelles studied by molecular dynamics simulations

In this study, mobility and structure of water molecules in Aerosol OT (bis(2-ethylhexyl) sulfosuccinate, AOT) reverse micelles with water content w0 = 5 and Na+, K+, Cs+ counterions have been explored with molecular dynamics (MD) simulations. Using the Faeder/Ladanyi model (J. Phys. Chem. B, 2000, 104, 1033) of the reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function, FS(Q,t), for water hydrogen atoms that could be measured in a quasielastic neutron scattering experiment. Separate intermediate scattering functions FRS(Q,t) and FCMS(Q,t) were determined for rotational and translational motion. We find that the decay of FCMS(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior arises from decreased water mobility for molecules close to the interface and from confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay, which is consistent with relatively rapid restricted rotation and slower rotational relaxation over the full angular range. Rotational relaxation is anisotropic, with the O-H bond short-time rotational mobility considerably higher than that of the molecular dipole. This behavior is related to the decreased density of water-water hydrogen bonds in the vicinity of the interface compared to core or bulk water. We find that the interfacial mobility of water molecules is quite different for the three counterion types, but that the core mobility exhibits weak counterion dependence. Differences in interfacial mobility are strongly correlated with structural features, especially ion-water coordination, and the extent of disruption by the counterions of the water hydrogen bond network.     

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

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

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.     

An investigation of water dynamics in binary mixtures of water and dimethyl sulfoxide

The motion of water molecules in mixtures of water and d6-dimethyl sulfoxide (DMSO) has been explored through molecular dynamics (MD) simulations using the SPC/E water model (J. Chem. Phys. 1987, 91, 6269) and the P2 DMSO model (J. Chem. Phys. 1993, 98, 8160). We evaluate the self-intermediate scattering functions, FS(Q,t), which are related by a Fourier transform to the incoherent structure factors, S(Q,omega), measured in quasielastic neutron scattering (QNS) experiments. We compare our results to recent QNS experiments on these mixtures reported by Bordallo et al. (J. Chem. Phys. 2004, 121, 12457). In addition to comparing the MD data to the experimental signals, which correspond to a convolution of S(Q,omega) with a resolution function, we examine the rotational and translational components of FS(Q,t) and investigate to what extent simulation results for the single-molecule dynamics follow the dynamical models that are used in the analysis of the experimental data. We find that the agreement between the experimental signal and the MD data is quite good and that the portion of FS(Q,t) due to translational dynamics is well represented by the jump-diffusion model. The model parameters and their composition dependence are in reasonable agreement with experiment, exhibiting similar trends in water mobility with composition. Specifically, we find that water motion is less hindered in water-rich and water-poor mixtures than it is near equimolar composition. We find that the extent of coupling between rotational and translational motion contributing to FS(Q,t) increases as the equimolar composition of the mixture is approached. Thus, the decoupling approximation, which is used to extract information on rotational relaxation from QNS spectra at higher momentum transfer (Q) values, becomes less accurate than that in water-rich or DMSO-rich mixtures. We also find that rotational relaxation deviates quite strongly from the isotropic rotational diffusion model. We explore this issue further by investigating the behavior of orientational time correlations for different unit vectors and corresponding to Legendre polynomials of orders 1-4. We find that the rotational time correlations of water molecules behave in a way that is more consistent with the extended jump rotation model recently proposed by Laage and Hynes (Science 2006, 311, 832).     

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).     

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.     

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. 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.     

Quasielastic small-angle neutron scattering from heavy water solutions of cyclodextrins

We present a model for quasielastic neutron scattering (QENS) by an aqueous solution of compact and inflexible molecules. This model accounts for time-dependent spatial pair correlations between the atoms of the same as well as of distinct molecules and includes all coherent and incoherent neutron scattering contributions. The extension of the static theory of the excluded volume effect [A. K. Soper, J. Phys.: Condens. Matter 9, 2399 (1997)] to the time-dependent (dynamic) case allows us to obtain simplified model expressions for QENS spectra in the low Q region in the uniform fluid approximation. The resulting expressions describe the quasielastic small-angle neutron scattering (QESANS) spectra of D(2)O solutions of native and methylated cyclodextrins well, yielding in particular translational and rotational diffusion coefficients of these compounds in aqueous solution. Finally, we discuss the full potential of the QESANS analysis (that is, beyond the uniform fluid approximation), in particular, the information on solute-solvent interactions (e.g., hydration shell properties) that such an analysis can provide, in principle.     

Entropy of water in the hydration layer of major and minor grooves of DNA

Transport properties (translational and rotational) of water in the two grooves of the B-DNA duplex are known to be different from those in the bulk. Here, we use a recently developed theoretical scheme to compute the entropies of water molecules in both of the grooves of DNA and compare them with that in the bulk. The scheme requires as input both translational and rotational velocity autocorrelation function (C(V)(t) and C(omega)(t), respectively) data. These velocity autocorrelation functions were computed from an atomistic MD simulation of a B-DNA duplex (36 base pairs long) in explicit water (TIP3P). The average values of the entropy of water at 300 K in both of the grooves of DNA (the TS value in the major groove is 6.71 kcal/mol and that in the minor groove is 6.41 kcal/mol) are found to be significantly lower than that in bulk water (the TS value is 7.27 kcal/mol). Thus, the entropic contribution to the free energy change (TDeltaS) of transferring a minor groove water molecule to the bulk is 0.86 kcal/mol and of transferring a major groove water to the bulk is 0.56 kcal/mol at 300 K, which is to be compared with 1.44 kcal/mol for melting of ice at 273 K. We also calculate the energy of interaction of each water molecule with the rest of the atoms in the system and hence calculate the chemical potential (Helmholtz free energy per water molecule, A = E - TS) in the different domains. The identical free energy value of water molecules in the different domains proves the robustness of the scheme. We propose that the configurational entropy of water in the grooves can be used as a measure of the mobility (or microviscosity) of water molecules in a given domain.     

Ion pairs of crystal violet in sodium bis(2-ethylhexyl)sulfosuccinate reverse micelles

The interfacial localization and the ion pair formation of the positively charged dye crystal violet (CV) in sodium bis(2-ethylhexyl)sulfosuccinate reverse micelles (AOT RMs) were studied by several structural and spectroscopic techniques and by quantum chemical calculations. The size and shape of the AOT RMs in the presence of CV were investigated by small-angle X-ray scattering, showing that CV does not significantly change the RM structure. CV localization as a function of the water to surfactant molar ratio (w(0)) was characterized by H(1) and (13)C NMR, indicating the close proximity of CV to the sulfosuccinate group of AOT at small and large w(0) values. These results were confirmed by calculation of magnetic shielding constants using the gauge-independent atomic orbital method with the HF/6-31G(d) basis set. Two different types of ion pairs between AOT and CV, i.e., contact ion pair (CIPs) and solvent-separated ion pair (SSIPs), were characterized by UV-vis spectroscopy and quantum chemical calculations using the semiempirical ZINDO-CI method. In nonpolar isotropic solvents CIPs are formed with an association constant (K(ASSOC)) of 2 x 10(4) mol(-1) L in isooctane and 750 mol(-1) L in chloroform. In AOT RMs at low w(0), CV-AOT CIPs are also formed. By increasing w(0), there is a sharp decrease in the CIP association free energy, and SSIPs are formed. (CV(+))(H(2)O)(AOT(-)) SSIPs are stable in the AOT RM up to the largest w(0) tested (w(0) = 33).     

Molecular simulation study of water mobility in aerosol-OT reverse micelles

In this work, we present results from molecular dynamics simulations on the single-molecule relaxation of water within reverse micelles (RMs) of different sizes formed by the surfactant aerosol-OT (AOT, sodium bis(2-ethylhexyl)sulfosuccinate) in isooctane. Results are presented for RM water content w(0) = [H(2)O]/[AOT] in the range from 2.0 to 7.5. We show that translational diffusion of water within the RM can, to a good approximation, be decoupled from the translation of the RM through the isooctane solvent. Water translational mobility within the RM is restricted by the water pool dimensions, and thus, the water mean-squared displacements (MSDs) level off in time. Comparison with models of diffusion in confined geometries shows that a version of the Gaussian confinement model with a biexponential decay of correlations provides a good fit to the MSDs, while a model of free diffusion within a sphere agrees less well with simulation results. We find that the local diffusivity is considerably reduced in the interfacial region, especially as w(0) decreases. Molecular orientational relaxation is monitored by examining the behavior of OH and dipole vectors. For both vectors, orientational relaxation slows down close to the interface and as w(0) decreases. For the OH vector, reorientation is strongly affected by the presence of charged species at the RM interface and these effects are especially pronounced for water molecules hydrogen-bonded to surfactant sites that serve as hydrogen-bond acceptors. For the dipole vector, orientational relaxation near the interface slows down more than that for the OH vector due mainly to the influence of ion-dipole interactions with the sodium counterions. We investigate water OH and dipole reorientation mechanisms by studying the w(0) and interfacial shell dependence of orientational time correlations for different Legendre polynomial orders.     

Experimental evidences for molecular origin of low-Q peak in neutron/x-ray scattering of 1-alkyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide ionic liquids

Short- and long-range liquid structures of [C(n)mIm(+)][TFSA(-)] with n = 2, 4, 6, 8, 10, and 12 have been studied by high-energy x-ray diffraction (HEXRD) and small-angle neutron scattering (SANS) experiments with the aid of MD simulations. Observed x-ray structure factor, S(Q), for the ionic liquids with the alkyl-chain length n > 6 exhibited a characteristic peak in the low-Q range of 0.2-0.4 A?(-1), indicating the heterogeneity of their ionic liquids. SANS profiles I(H)(Q) and I(D)(Q) for the normal and the alkyl group deuterated ionic liquids, respectively, showed significant peaks for n = 10 and 12 without no form factor component for large spherical or spheroidal aggregates like micelles in solution. The peaks for n = 10 and 12 evidently disappeared in the difference SANS profiles ΔI(Q) [=I(D)(Q) - I(H)(Q)], although that for n = 12 slightly remained. This suggests that the long-range correlations originated from the alkyl groups hardly contribute to the low-Q peak intensity in SANS. To reveal molecular origin of the low-Q peak, we introduce here a new function; x-ray structure factor intensity at a given Q as a function of r, S(Q) (peak)(r). The S(Q) (peak)(r) function suggests that the observed low-Q peak intensity depending on n is originated from liquid structures at two r-region of 5-8 and 8-15 A? for all ionic liquids examined except for n = 12. Atomistic MD simulations are consistent with the HEXRD and SANS experiments, and then we discussed the relationship between both variations of low-Q peak and real-space structure with lengthening the alkyl group of the C(n)mIm.     

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 solvent and rotational relaxation of glycerol in the nanocavity of reverse micelles

The dynamics of solvent and rotational relaxation of Coumarin 480 and Coumarin 490 in glycerol containing bis-2-ethyl hexyl sulfosuccinate sodium salt (AOT) reverse micelles have been investigated with steady-state and time-resolved fluorescence spectroscopy. We observed slower solvent relaxation of glycerol confined in the nanocavity of AOT reverse micelles compared to that in pure glycerol. However, the slowing down in the solvation time on going from neat glycerol to glycerol confined reverse micelles is not comparable to that on going from pure water or acetonitrile to water or acetonitrile confined AOT reverse micellar aggregates. While solvent relaxation times were found to decrease with increasing glycerol content in the reverse micellar pool, rotational relaxation times were found to increase with increase in glycerol content.     

Reduced coupling of water molecules near the surface of reverse micelles

We report on vibrational dynamics of water near the surface of AOT reverse micelles studied by narrow-band excitation, mid-IR pump-probe spectroscopy. Evidence of OH-stretch frequency splitting into the symmetric and asymmetric modes is clearly observed for the interfacial H(2)O molecules. The polarization memory of interfacial waters is preserved over an exceptionally extended >10 ps timescale which is a factor of 100 longer than in bulk water. These observations point towards negligibly small intermolecular vibrational coupling between the water molecules as well as strongly reduced water rotational mobility within the interfacial water layer.     

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.     

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 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.