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.     

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

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.     

Reorientational dynamics of the dodecahydro-closo-dodecaborate anion in Cs2B12H12

Rapid reorientational motions of the B(12)H(12)(2-) icosahedral anion, a key intermediate in borohydride dehydrogenation, are revealed by quasielastic neutron scattering (QENS) measurements of Cs(2)B(12)H(12) between 430 and 530 K. At 430 K, over the range of momentum transfers collected, the elastic incoherent structure factor (EISF) is consistent with a model for reorientational jumps about a single molecular axis. At temperatures of 480 K and higher, however, the EISF suggests the emergence of multiaxis reorientation by dynamically similar, independent jumps about two axes, on average, preserving crystallographic order. Alternatively, if one assumes that the anions are undergoing temperature-dependent rotational trapping, then the EISF is also consistent with a jump model involving a temperature-dependent mobile fraction of anions statistically tumbling between discrete crystallographic sites. Although neutron vibrational spectra demonstrate that the anion torsional modes soften dramatically with increasing temperature, the QENS-derived activation energy of 333 meV for reorientation clearly shows that the anions are not undergoing isotropic rotational diffusion.     

Thermodynamic, structural, and dynamic properties of supercooled water confined in mesoporous MCM-41 studied with calorimetric, neutron diffraction, and neutron spin echo measurements

Thermodynamic, structural, and dynamic properties of heavy water (D(2)O) confined in mesoporous silica glass MCM-41 C10, C12, and C14 were investigated by differential scanning calorimetry, neutron diffraction, and neutron spin echo (NSE) measurements, respectively. The DSC data showed that no crystallization of D(2)O confined in C10 occurs in a temperature range between 298 and 180 K, and that crystalline ice is formed at 204 and 221 K for C12 and C14, respectively. For C10, the neutron radial distribution functions of confined D(2)O suggested a structural change in the supercooled state between 223 and 173 K. For C10 sample, it has been found that the tetrahedral-like water structure is partially enhanced in the central part of pores at 173 K. For all the samples, the intermediate scattering functions from the NSE measurements are fitted by the Kohlrausch-Williams-Watts stretched exponential function which implies that confined supercooled D(2)O exhibits a wide distribution of relaxation times. For C10, C12, and C14 samples, between 298 and 240 K, the relaxation times of supercooled D(2)O follow remarkably well the Vogel-Fulcher-Tamman equation; for C10 sample, below 240 K, the relaxation times of nonfreezing D(2)O show an Arrhenius type behavior. From the present experimental results on calorimetric, structural, and dynamic properties, it has been concluded that supercooled D(2)O confined in MCM-41 C10 experiences a transition from high-density to low-density hydrogen-bonded structure at around 229 K.     

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.     

Structural and dynamical analysis of the hydration of the Alzheimer's beta-amyloid peptide

An analysis of the water molecules in the first solvation shell obtained from the molecular dynamics simulation of the amyloid beta(10-35)NH2 peptide and the amyloid beta(10-35)NH2E22Q "Dutch" mutant peptide is presented. The structure, energetics, and dynamics of water in the hydration shell have been investigated using a variety of measures, including the hydrogen bond network, the water residence times for all the peptide residues, the diffusion constant, experimentally determined HN amide proton exchange, and the transition probabilities for water to move from one residue to another or into the bulk. The results of the study indicate that: (1) the water molecules at the peptide-solvent interface are organized in an ordered structure similar for the two peptide systems but different from that of the bulk, (2) the peptide structure inhibits diffusion perpendicular to the peptide surface by a factor of 3 to 5 relative to diffusion parallel to the peptide surface, which is comparable to diffusion of bulk water, (3) water in the first solvation shell shows dynamical relaxation on fast (1-2 ps) and slow (10-40 ps) time scales, (4) a novel solvent relaxation master equation is shown to capture the details of the fast relaxation of water in the peptide's first solvation shell, (5) the interaction between the peptide and the solvent is stronger in the wild type than in the E22Q mutant peptide, in agreement with earlier results obtained from computer simulations [Massi, F.; Straub, J. E. Biophys J 2001, 81, 697] correlated with the observed enhanced activity of the E22Q mutant peptide.     

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.     

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.     

The dynamics of glass-forming polymers in the microscopic-mesoscopic time scale. A quasielastic neutron scattering phenomenological approach

The dynamics of three glass-forming polymers, PVC, PB and PI, has been investigated by time of flight (TOF) neutron scattering in a time scale from 10⁻¹³s to 10¹¹s looking for the crossover from microscopic dynamics to segmental dynamics (α relaxation). A new analysis procedure has been applied to TOF data in order to separate harmonic vibrational and relaxational contributions. Due to the involved assumptions, this procedure can be considered only as a first approximation adequate for the case of “fragile” systems (in the Angell's meaning) like the polymers here investigated. The behaviour obtained was the same for the three polymers studied. The intermediate scattering function corresponding to the relaxational dynamics, I_o(Q,t), shows two different dynamical regimes separated by a crossover time t_c (≈ 2 ps), which hardly depends on Q (momentum transfer) and temperature. At t < t<_c, I_o(Q,t) displays a Debye-like behaviour (exponential decay). The activation energy found for the relaxation time corresponding to this regime was in the range of 2–5 Kcal/mol, i.e., in the range of the activation energy for local conformational transitions in isolated macromolecular chains. At t > t_c and, at least at high temperature, I_o(Q,t) shows a Kohlrausch-Williams-Watts (KWW) behaviour similar to the obtained one by means of backscattering neutron techniques in the mesoscopic time scale (10⁻¹¹s to 10⁻⁷s) and dielectric measurements in the macroscopic time scale (10⁻⁷s to 10°s). This KWW regime can be associated to the segmental dynamics involved in the α relaxation. A phenomenological interpretation is outlined. In this framework, the Debye-like regime is interpreted to be the segmental dynamics free from intermolecular hindrances. Therefore, t_c should be the time at which intermolecular interactions start to play a significant role concerning to the segmental dynamics. This interpretation recalls some of the basic ideas of the so called “Coupling-Model” proposed a long time ago by Ngai.     

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.     

Temperature-dependent dynamics of water confined in nafion membranes

We performed a neutron scattering study to investigate the dynamical behavior of water absorbed in Nafion at low hydration level as a function of temperature in the range 200-300 K. To single out the spectral contribution of the confined water, the measurements were done on samples hydrated with both H(2)O and D(2)O. Due to the strong incoherent scattering cross section of hydrogen atoms with respect to deuterium, in the difference spectra, the contribution from the Nafion membrane is subtracted out and the signal originates essentially from protons in the liquid phase. The main quantities we extracted as a function of the momentum transfer are the elastic incoherent structure factor (EISF) and the line width of the quasielastic component. Their trend suggests that the motion of hydrogen atoms can be schematized as a random jumping inside a confining region, which can be related to the boundaries of the space where water molecules move in the cluster they form around the sulfonic acid site. Through the calculated EISF, we obtained information on the size of such a region, which increases up to 260 K and then attains a constant value. Above this temperature, the number of water protons that are dynamically activated in the accessible time window increases with a faster rate. The jump diffusion dynamics is characterized by a typical jumping time which is stable at 5.3 ps up to approximately 260 K and then gradually decreases. The ensemble of the findings indicates that, within the limits of the energy resolution of the present experiment, water absorbed in the Nafion membrane undergoes a dynamical transition at around 260 K. We discuss the possible relationship of this dynamical onset with the behavior of the electrical conductivity of the membrane as a function of the temperature.     

Investigation of isothermal water infiltration into fired clay brick by scattered neutrons

A method based on neutron scattering was proposed to investigate isothermal water infiltration in porous media. Two different kinds of fired clay bricks were investigated. While the sample absorb water, scattered neutrons from the different wetted regions, along the flow direction were continuously recorded. The results were discussed in terms of the theory of water infiltration in unsaturated porous media as well as by an anomalous diffusion approach. It was shown that the infiltration process in the Canadian clay brick (CCB) is Fickian and the water diffusivity was analytically determined, while it is non-Fickian in the Egyptian clay brick (ECB). The infiltration process in ECB can be modeled as a two stage Fickian process, at the low and high absorption times. The anomalous diffusion approach failed to describe the diffusion process in the ECB at all water contents.     

Dynamic properties of microemulsions modified with homopolymers and diblock copolymers: the determination of bending moduli and renormalization effects

The properties of bicontinuous microemulsions, consisting of water, oil, and a surfactant, can be modified by the addition of diblock copolymers (boosting effect) and homopolymers (inverse boosting effect) or a combination of both. Here, the influence of the addition of homopolymers (PEP(X) and PEO(X), X=5k or 10k molecular weight) on the dynamics of the surfactant layer is studied with neutron spin echo spectroscopy (NSE). Combining the results with the previous findings for diblock copolymers allows for a better separation of viscosity and bending modulus effects. With the addition of homopolymers, a significant increase of the relaxation rate compared to the pure microemulsion has been observed. The influence on the bending rigidity kappa is measured with NSE experiments. Homopolymer addition reduces kappa by up to Deltakappa approximately -0.5k(B)T, whereas the diblock copolymer yields an increase of kappa by approximately 0.3k(B)T. Comparison of the bending moduli that are obtained by analysis of the dynamics to those obtained from small angle neutron scattering (SANS) sheds light on the different renormalization length scales for NSE and SANS. Variation of the surfactant concentration at otherwise constant conditions of homopolymer or diblock-copolymer concentration shows that NSE results are leading to the pure bending rigidity, while the renormalized one is measured with SANS.     

Phenylene ring dynamics in bisphenol-A-polysulfone by neutron scattering

We have investigated the dynamics of phenylene rings in a glassy polysulfone (bisphenol-A-polysulfone) by means of quasielastic neutron scattering. Nowadays it is well known that these molecular motions are directly connected with the mechanical properties of engineering thermoplastics in general. The particular system investigated by us has the advantage that by selective deuteration of the methyl groups, the neutron scattering measured is dominated by the incoherent contribution from the protons in the phenylene rings. In this way, the dynamics of such molecular groups can be experimentally isolated. Two different types of neutron spectrometers: time of flight and backscattering, were used in order to cover a wide dynamic range, which extends from microscopic (10(-13) s) to mesoscopic (10(-9) s) times. Moreover, neutron diffraction experiments with polarization analysis were also carried out in order to characterize the structural features of the sample investigated. Fast oscillations of increasing amplitude with temperature and pi-flips are identified for phenylene rings motions. Due to the structural disorder characteristic of the amorphous state, both molecular motions display a broad distribution of relaxation times, which spreads over several orders of magnitude. Based on the results obtained, we propose a model for phenylene rings dynamics, which combines the two kinds of molecular motions identified. This model nicely describes the neutron scattering results in the whole dynamic range investigated.     

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