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.     

Diffusive solvent dynamics in a polymer gel electrolyte studied by quasielastic neutron scattering

A quasielastic neutron scattering study has been performed on a polymer gel electrolyte consisting of lithium perchlorate dissolved in ethylene carbonate/propylene carbonate and stabilized with poly(methyl methacrylate). The dynamics of the solvent, which is crucial for the ion conduction in this system, was probed using the hydrogen/deuterium contrast variation method with nondeuterated solvent and a deuterated polymer matrix. Two relaxation processes of the solvent were studied in the 10-400 microeV range at different temperatures. From analysis of the momentum transfer dependence of the processes we conclude that the faster process ( approximately 100 microeV) is related to rotational diffusion of the solvent and the slower process ( approximately 10 microeV) to translational diffusion of the solvent. The translational diffusion is found to be similar to the diffusion in the corresponding liquid electrolyte at short distances, but geometrically constrained by the polymer matrix at distances beyond approximately 5 A. The study indicates that the hindered diffusion of the solvent on a length scale of the polymer network interchain distance ( approximately 5-20 A) is sufficient to explain the reduced macroscopic diffusivity and ion conductivity of the gel electrolyte compared to the liquid electrolyte.     

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.     

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.     

Self-diffusion in molecular liquids: medium-chain n-alkanes and coenzyme Q10 studied by quasielastic neutron scattering

A systematic time-of-flight quasielastic neutron scattering (TOF-QENS) study on diffusion of n-alkanes in a melt is presented for the first time. As another example of a medium-chain molecule, coenzyme Q(10) is investigated in the same way. The data were evaluated both in the frequency and in the time domain. TOF-QENS data can be satisfactorily described by different models, and it turned out that the determined diffusion coefficients are largely independent of the applied model. The derived diffusion coefficients are compared with values measured by pulsed-field gradient nuclear magnetic resonance (PFG-NMR). With increasing chain length, an increasing difference between the TOF-QENS diffusion coefficient and the PFG-NMR diffusion coefficient is observed. This discrepancy in the diffusion coefficients is most likely due to a change of the diffusion mechanism on a nanometer length scale for molecules of medium-chain length.     

Coherent dynamics of meta-toluidine investigated by quasielastic neutron scattering

The coherent dynamics of a typical fragile glass former, meta-toluidine, was investigated at the molecular level using quasielastic neutron scattering, with time-of-flight and neutron spin echo spectrometers. It is well known that the static structure factor of meta-toluidine shows a prepeak originating from clustering of the molecules through hydrogen bonding between the amine groups. The dynamics of meta-toluidine was measured for several values of the wavevector transfer Q, which is equivalent to an inverse length scale, in a range encompassing the prepeak and the structure factor peak. Data were collected in the temperature range corresponding to the liquid and supercooled states, down to the glass transition. At least two dynamical processes were identified. This paper focuses on the slowest relaxation process in the system, the α-relaxation, which was found to scale with the macroscopic shear viscosity at all the investigated Q values. No evidence of "de Gennes" narrowing associated with the prepeak was observed, in contrast with what happens at the Q value corresponding to the interparticle distance. Moreover, using partially deuterated samples, the dynamics of the clusters was found to be correlated to the single-particle dynamics of the meta-toluidine molecules.     

In-plane pyridinium cation reorientation in bis-thiourea chloride, bromide and iodide: quasielastic neutron scattering combined with molecular dynamics simulations

We have studied the dynamics of bis-thiourea pyridinium chloride and bromide by means of quasielastic neutron scattering (QENS). The QENS data allow describing the geometry of the in-plane motion of the pyridinium cation and reveal that it is similar to the motion previously observed in bis-thiourea pyridinium iodide. Molecular dynamics (MD) simulations have been performed to investigate the cation dynamics on the high temperature phase of the full series of compounds: bis-thiourea pyridinium chloride, bromide and iodide. Three different models of intermolecular potential have been tested and the agreement between the simulated and experimental elastic incoherent structure factors (EISFs) is used to select the more realistic one. The detailed analysis of the MD results indicates that Coulombic interactions together with the formation of hydrogen bonds between the pyridinium cation and the host sublattice influence strongly the geometry of the in-plane cation reorientation.     

Dynamics of surface water in ZrO2 studied by quasielastic neutron scattering

A quasielastic neutron scattering experiment has revealed the dynamics of surface water in a high surface area zirconium oxide in the temperature range of 300-360 K. The characteristic times of the rotational (picoseconds) and translational (tens of picoseconds) components of diffusion motion are well separated. The rotational correlation time shows an Arrhenius-type behavior with an activation energy of 4.48 kJ/mol, which is lower compared to bulk water. The rotational diffusion at room temperature is slower by about a factor of 2 compared to bulk water, whereas the translational diffusion slows down by a factor of 40. In contrast to bulk water, the translational correlation time exhibits an Arrhenius-type temperature dependence with an activation energy of 11.38 kJ/mol. Comparison of different models for jump diffusion processes suggests that water molecules perform two-dimensional jumps at a well-defined, almost temperature-independent distance of 4.21-4.32 A. Such a large jump distance indicates a low molecular density of the layer of diffusing molecules. We argue that undissociated water molecules on an average form two hydrations layers on top of the surface layer of hydroxyl groups, and all the layers have similar molecular density. Quasielastic neutron scattering experiment assesses the dynamics of the outermost hydration layer, whereas slower motion of the water molecules in the inner hydration layer contributes to the elastic signal.     

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

Incoherent quasielastic neutron scattering study of molecular dynamics of 4-n-octyl-4'-cyanobiphenyl

We report incoherent quasielastic neutron scattering experiments on the thermotropic liquid crystal 4-n-octyl-4'-cyanobiphenyl. The combination of time-of-flight and backscattering data allows analysis of the intermediate scattering function over about three decades of relaxation times. Translational diffusion and uniaxial molecular rotations are clearly identified as the major relaxation processes in, respectively, the nanosecond and picosecond time scales.     

Correlated atomic motions in liquid deuterium fluoride studied by coherent quasielastic neutron scattering

The collective dynamics of liquid deuterium fluoride are studied by means of high-resolution quasielastic and inelastic neutron scattering over a range of four decades in energy transfer. The spectra show a low-energy coherent quasielastic component which arises from correlated stochastic motions as well as a broad inelastic feature originating from overdamped density oscillations. While these results are at variance with previous works which report on the presence of propagating collective modes, they are fully consistent with neutron diffraction, nuclear magnetic resonance, and infrared/Raman experiments on this prototypical hydrogen-bonded fluid.     

Hydration of tricalcium and dicalcium silicate mixtures studied using quasielastic neutron scattering

Quasielastic neutron scattering was used to study the hydration reaction of tricalcium and dicalcium silicate mixtures by following the fixation of hydrogen into the reaction products, and by applying hydration models to the data. The reaction kinetics were well-described by an Avrami-derived model for the nucleation and growth regime during early hydration times and a diffusion-limited model for later periods. This study showed that the hydration reaction is not a simple linear combination of the reactions for the individual components. Compressive strength tests correlated with the neutron scattering data, suggesting that the details of the interaction affect the microstructure and therefore the strength of the product. Results suggest that favorable reaction mechanics provide optimal strength when an 80-95% tricalcium silicate and 20-5% dicalcium silicate mixture is used.     

Patchy worm-like micelles: solution structure studied by small-angle neutron scattering

Triblock terpolymers exhibit a rich self-organization behavior including the formation of fascinating cylindrical core-shell structures with a phase separated corona. After crystallization-induced self-assembly of polystyrene-block-polyethylene-block-poly(methyl methacrylate) triblock terpolymers (abbreviated as SEMs = Styrene-Ethylene-Methacrylates) from solution, worm-like core-shell micelles with a patchy corona of polystyrene and poly(methyl methacrylate) were observed by transmission electron microscopy. However, the solution structure is still a matter of debate. Here, we present a method to distinguish in situ between a Janus-type (two faced) and a patchy (multiple compartments) configuration of the corona. To discriminate between both models the scattering intensity must be determined mainly by one corona compartment. Contrast variation in small-angle neutron scattering enables us to focus on one compartment of the worm-like micelles. The results validate the existence of the patchy structure also in solution.     

Deswelling of microfibril bundles in drying wood studied by small-angle neutron scattering and molecular dynamics

Cellulose - Structural changes of cellulose microfibrils and microfibril bundles in unmodified spruce cell wall due to drying in air were investigated using time-resolved small-angle neutron...     

Solvent and lipid dynamics of hydrated lipid bilayers by incoherent quasielastic neutron scattering

The microscopic dynamics of the planar, multilamellar lipid bilayer system 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) has been investigated using quasielastic neutron scattering. The DMPC was hydrated to a level corresponding to approximately nine water molecules per lipid molecule. Selective deuteration has been used to separately extract the dynamics of the water, the acyl chains, and the polar head groups from the strong incoherent scattering of the remaining hydrogen atoms. Furthermore, the motions parallel and perpendicular to the bilayers were probed by using two different sample orientations relative to the incident neutron beam. For both sample orientations, the results showed an onset of water motions at 260 K on the experimental time scale of about 100 ps. From lack of wave-vector dependence of the onset temperature for water motions, it is evident that the observed water dynamics is of mainly rotational character at such low temperatures. At 290 K, i.e., slightly below the gel-to-liquid transition around 295 K, the nature of the water dynamics had changed to a more translational character, well described by a jump-diffusion model. On the limited experimental time and length (about 10 A) scales, this jump-diffusion process was isotropic, despite the very anisotropic system. The acyl chains exhibited a weak onset of anharmonic motions already at 120 K, probably due to conformational changes (trans-gauche and/or syn-anti) in the plane of the lipid bilayers. Other anharmonic motions were not observed on the experimental time scale until temperature had been reached above the gel-to-liquid transition around 295 K, where the acyl chains start to show more substantial motions.     

Water dynamics in hardened ordinary Portland cement paste or concrete: from quasielastic neutron scattering

Portland cement reacts with water to form an amorphous paste through a chemical reaction called hydration. In concrete the formation of pastes causes the mix to harden and gain strength to form a rock-like mass. Within this process lies the key to a remarkable peculiarity of concrete: it is plastic and soft when newly mixed, strong and durable when hardened. These qualities explain why one material, concrete, can build skyscrapers, bridges, sidewalks and superhighways, houses, and dams. The character of the concrete is determined by the quality of the paste. Creep and shrinkage of concrete specimens occur during the loss and gain of water from cement paste. To better understand the role of water in mature concrete, a series of quasielastic neutron scattering (QENS) experiments were carried out on cement pastes with water/cement ratio varying between 0.32 and 0.6. The samples were cured for about 28 days in sealed containers so that the initial water content would not change. These experiments were carried out with an actual sample of Portland cement rather than with the components of cement studied by other workers. The QENS spectra differentiated between three different water interactions: water that was chemically bound into the cement paste, the physically bound or "glassy water" that interacted with the surface of the gel pores in the paste, and unbound water molecules that are confined within the larger capillary pores of cement paste. The dynamics of the "glassy" and "unboud" water in an extended time scale, from a hundred picoseconds to a few nanoseconds, could be clearly differentiated from the data. While the observed motions on the picosecond time scale are mainly stochastic reorientations of the water molecules, the dynamics observed on the nanosecond range can be attributed to long-range diffusion. Diffusive motion was characterized by diffusion constants in the range of (0.6-2) 10(-9) m(2)/s, with significant reduction compared to the rate of diffusion for bulk water. This reduction of the water diffusion is discussed in terms of the interaction of the water with the calcium silicate gel and the ions present in the pore water.     

Investigation of chiral recognition by molecular micelles with molecular dynamics simulations

Molecular dynamics simulations were used to characterize the binding of the chiral drugs chlorthalidone and lorazepam to the molecular micelle poly-(sodium undecyl-(L)-leucine-valine). The project’s goal was to characterize the nature of chiral recognition in capillary electrophoresis separations that use molecular micelles as the chiral selector. The shapes and charge distributions of the chiral molecules investigated, their orientations within the molecular micelle chiral binding pockets, and the formation of stereoselective intermolecular hydrogen bonds with the molecular micelle were all found to play key roles in determining where and how lorazepam and chlorthalidone enantiomers interacted with the molecular micelle.     

Molecular dynamics of glucose in solution: a quasielastic neutron scattering study

The molecular dynamics of glucose dissolved in heavy water have been investigated at 280 K by the technique of quasielastic neutron scattering. The scattering was described by a dynamic structure factor that accounts for decoupled diffusive jumps and free rotational motions of the glucose molecules. With increasing glucose concentration, the diffusion constant decreases by a factor five and the time between jumps increases considerably. Our observations validate theoretical predictions concerning the impact of concentration on the environment of a glucose molecule and the formation of cages made by neighboring glucose molecules at higher concentrations.     

Core-shell structure of degradable, thermosensitive polymeric micelles studied by small-angle neutron scattering

The structure of assemblies of block copolymers composed of thermosensitive, biodegradable poly(N-(2-hydroxypropyl) methacrylamide-dilactate) and poly(ethylene glycol) (pHPMAmDL-b-PEG) has been studied by small-angle neutron scattering (SANS). Three amphiphilic copolymers with a fixed PEG of 5 kDa and a partially deuterated pHPMAmDL(d) block of 6700, 10400, or 21200 Da were used to form micelles in aqueous media by heating the polymeric solution from below to above the cloud point temperature (around 10 degrees C) of the thermosensitive block. Simultaneous and quantitative analysis of the scattering cross sections obtained at three different solvent contrasts is expedited using core-shell model, which assumed a homogeneous core of uniform scattering length density. The mean core radius increased from 13 to 18.5 nm with the molecular weight of the pHPMAmDL(d) block, while the thickness of the stabilizing PEG layer was around 8 nm for the three investigated assemblies. In addition, the volume fraction values of the stabilizing PEG chains in the shell are low and decreased from 31% to 14% with increasing the size of pHPMAmDL(d) block which shows that the shell layer of the assemblies is highly hydrated. The corresponding PEG chain grafting densities decreased from 0.22 to 0.11 nm-2 and the distance between PEG chains on the nanoparticles surface increased from 2.4 to 3.4 nm. The pHPMAmDL-b-PEG micelles showed a controlled instability due to hydrolysis of the lactic acid side groups in the thermosensitive block; that is, an increase of the degradation time leads to an increase of the size of the core which becomes less hydrophobic and consequently more hydrated. Neutron experiments supplied accurate information on how the size of the core and the micelle's aggregation number changed with the incubation time. This feature and the initially small size and dense structure in aqueous solution make the polymeric micelles suitable as carriers for hydrophobic drugs.