Quasi-elastic neutron scattering study of dimethyl-sulfoxide–water mixtures: Probing molecular mobility in a nonideal solution

Confinement Effects on Glass-Forming Aqueous Dimethyl Sulfoxide Solutions

Combining broadband dielectric spectroscopy and nuclear magnetic resonance studies, we analyze the reorientation dynamics and the translational diffusion associated with the glassy slowdown of the eutectic aqueous dimethyl sulfoxide solution in nano-sized confinements, explicitly, in silica pores with different diameters and in ficoll and lysozyme matrices at different concentrations. We observe that both rotational and diffusive dynamics are slower and more heterogeneous in the confinements than in the bulk but the degree of these effects depends on the properties of the confinement and differs for the components of the solution. For the hard and the soft matrices, the slowdown and the heterogeneity become more prominent when the size of the confinement is reduced. In addition, the dynamics are more retarded for dimethyl sulfoxide than for water, implying specific guest-host interactions. Moreover, we find that the temperature dependence of the reorientation dynamics and of the translational diffusion differs in severe confinements, indicating a breakdown of the Stokes–Einstein–Debye relation. It is discussed to what extent these confinement effects can be rationalized in the framework of core-shell models, which assume bulk-like and slowed-down motions in central and interfacial confinement regions, respectively.

On the crystal structures and phase transitions of hydrates in the binary dimethyl sulfoxide–water system

Neutron powder diffraction data have been collected from a series of flash-frozen aqueous solutions of dimethyl sulfoxide (DMSO) with concentrations between 25 and 66.7 mol% DMSO. These reveal the existence of three stoichiometric hydrates, which crystallize on warming between 175 and 195 K. DMSO trihydrate crystallizes in the monoclinic space group P21/c, with unit-cell parameters at 195 K of a = 10.26619 (3), b = 7.01113 (2), c = 10.06897 (3) Å, β = 101.5030 (2)° and V = 710.183 (3) Å3 (Z = 4). Two of the symmetry-inequivalent water molecules form a sheet of tiled four- and eight-sided rings; the DMSO molecules are sandwiched between these sheets and linked along the b axis by the third water molecule to generate water–DMSO–water tapes. Two different polymorphs of DMSO dihydrate have been identified. The α phase is monoclinic (space group P21/c), with unit-cell parameters at 175 K of a = 6.30304 (4), b = 9.05700 (5), c = 11.22013 (7) Å, β = 105.9691 (4)° and V = 615.802 (4) Å3 (Z = 4). Its structure contains water–DMSO–water chains, but these are polymerized in such a manner as to form sheets of reniform eight-sided rings, with the methyl groups extending on either side of the sheet. On warming above 198 K, α-DMSO·2H2O undergoes a solid-state transformation to a mixture of DMSO·3H2O + anhydrous DMSO, and there is then a stable eutectic between these two phases at ∼203 K. The β-phase of DMSO dihydrate has been observed in a rapidly frozen eutectic melt and in very DMSO-rich mixtures. It is observed to be unstable with respect to the α-phase; above ∼180 K, β-DMSO·2H2O converts irreversibly to α-DMSO·2H2O. At 175 K, the lattice parameters of β-DMSO·2H2O are a = 6.17448 (10), b = 11.61635 (16), c = 8.66530 (12) Å, β = 101.663 (1)° and V = 608.684 (10) Å3 (Z = 4), hence this polymorph is just 1.16% denser than the α-phase under identical conditions. Like the other two hydrates, the space group appears likely, on the basis of systematic absences, to be P21/c, but the structure has not yet been determined. Our results reconcile 60 years of contradictory interpretations of the phase relations in the binary DMSO–water system, particularly between mole fractions of 0.25–0.50, and confirm empirical and theoretical studies of the liquid structure around the eutectic composition (33.33 mol% DMSO).

Picosecond self-diffusion in ethanol-water mixtures

We report the self-diffusion in ethanol-water mixtures as a function of the water-ethanol ratio measured at different temperatures using quasi-elastic neutron spectroscopy (QENS). For our protiated samples, QENS is mainly sensitive to the dominant ensemble-averaged incoherent scattering from the hydrogen atoms of the liquid mixtures. The energy range and resolution render our experiment sensitive to the picosecond time scale and nanometer length scale. These observation scales complement different scales accessible by nuclear magnetic resonance techniques. Subsequent to testing different models, we find that a simple jump-diffusion model averaging over both types of molecules, water and ethanol, best fits our data.

Molecular dynamics study of unfolding of lysozyme in water and its mixtures with dimethyl sulfoxide

All-atom explicit solvent molecular dynamics was used to study the process of unfolding of hen egg white lysozyme in water and mixtures of water with dimethyl sulfoxide at different compositions. We have determined the kinetic parameters of unfolding at a constant temperature 450K. For each run, the time of disruption of the tertiary structure of lysozyme t_u was defined as the moment when a certain structural criterion computed from the trajectory reaches its critical value. A good agreement is observed between the results obtained using several different criteria. The secondary structure according to DSSP calculations is found to be partially unfolded to the moment of disruption of tertiary structure, but some of its elements keep for a long time after that. The values of t_u averaged over ten 30ns-long trajectories for each solvent composition are shown to decrease very rapidly with addition of dimethyl sulfoxide, and rather small amounts of dimethyl sulfoxide are found to change the pathway of unfolding. In pure water, despite the loss of tertiary contacts and disruption of secondary structure elements, the protein preserves its compact globular state at least over 130ns of simulation, while even at 5mol percents of dimethyl sulfoxide it loses its compactness within 30ns. The proposed methodology is a generally applicable tool to quantify the rate of protein unfolding in simulation studies.

Local lateral environment of the molecules at the surface of DMSO-water mixtures

Molecular dynamics simulations of the liquid-vapour interface of dimethyl sulphoxide (DMSO)-water mixtures of 11 different compositions, including two neat systems are performed on the canonical (N, V, T) ensemble at 298 K. The molecules constituting the surface layer of these systems are selected by means of the identification of the truly interfacial molecules (ITIM) method, and their local lateral environment at the liquid surface is investigated by performing Voronoi analysis. The obtained results reveal that both molecules prefer to be in a mixed local environment, consisting of both kinds of molecules, at the liquid surface, and this preference is even stronger here than in the bulk liquid phase. Neat-like patches, in which a molecule is surrounded by like neighbours, are not found. However, vacancies that are surrounded solely by water molecules are observed at the liquid surface. Our results show that strongly hydrogen bonded DMSO·H2O complexes, known to exist in the bulk phase of these mixtures, are absent from the liquid surface.

The local environment of the molecules in water–DMSO mixtures, as seen from computer simulations and Voronoi polyhedra analysis

The local structure of DMSO–water mixtures is studied by computer simulation and Voronoi analysis.

Dielectric constant of polyhydric alcohol-DMSO mixture solution at the microwave frequency

Dielectric spectrum of polyhydric alcohol (1,2,3)-DMSO (dimethyl sulfoxide) mixtures, at full concentration, have been determined by the dielectric relaxation spectroscopy (DRS) method at frequencies from 20 MHz to 20 GHz at room temperature. The mixture behavior is described according to four Davidson-Cole terms whose evolution with composition is analyzed. The binding energies and hydrogen bond (HB) numbers between solute-solute and solvent-solute pairs are obtained from the permittivity using the Luzar model. The average HBs number associated with DMSO decreases with increasing mole fraction of DMSO. The binding energy of solute-solute (E11) and solvent-solute (E12) interaction decreases with the increased numbers of carbon atoms in the alcohols.

A model for predicting the permeation of dimethyl sulfoxide into articular cartilage, and its application to the liquidus-tracking method

Long-term storage of articular cartilage (AC) has excited great interest due to the practical surgical significance of this tissue. The liquidus-tracking (LT) method developed by Pegg et al. (2006) [29] for vitreous preservation of AC achieved reasonable survival of post-warming chondrocytes in situ, but the design of the entire procedure was more dependent on trial and error. Mathematical modeling would help to better understand the LT process, and thereby make possible improvements to attain higher cell survival. Mass transfer plays a dominant role in the LT process. In the present study, a diffusion model based on the free-volume theory and the Flory-Huggins thermodynamics theory was developed to predict the permeation of dimethyl sulfoxide (Me2SO) into AC. A comparison between the predicted mean concentration of Me2SO in the AC disc and the experimental data over wide temperature and concentration ranges [-30 to 37 °C, 10 to 64.5% (w/w)] shows that the developed model can accurately describe the permeation of Me2SO into AC [coefficient of determination (R(2)): 0.951-1.000, mean relative error (MRE): 0.8-12.8%]. With this model, the spatial and temporal distribution of Me2SO in the AC disc during a loading/unloading process can be obtained. Application of the model to Pegg et al.'s LT procedure revealed that the liquidus line is virtually not followed for the center part of the AC disc. The presently developed model will be a useful tool in the analysis and design of the LT method for vitreous preservation of AC.

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.

Concrete and Cement Paste Studied by Quasi-Elastic Neutron Scattering

In a world where the effects of climate change on weather patterns is accepted as real and serious, the problem of decreasing the production of carbon dioxide is perceived as increasingly important. The cement industry produces 5–7% of the world’s carbon dioxide emission and its survival will depend on improvements in the production of concrete which will be both more durable and require less carbon dioxide per unit of manufacture than the currently produced concrete. The durability of concrete is related to its ability to limit fluid transmission and knowledge of how to reduce the rate at which water will be transmitted through cement paste is critical to improving durability. However, because of the complex chemical and physical nature of cement pastes, understanding water mobility is a great challenge. Many techniques are not applicable simply because they are not sensitive to the range of size from angstroms to microns and the extent of water interaction with the cement where water can either be chemically bound at hydroxyls or physically free in large pores. In this review paper, we present the most up to date results on the physical chemistry of the water/ cement paste interactions studied by quasi-elastic neutron scattering. These results bring new insight to the mobility of water in the gel pores, the small pores (radius less than 50Å) that control the rate of water transmission in the cement pastes from which high quality concrete will be made.

Hindered water motions in hardened cement pastes investigated over broad time and length scales

We investigated the dynamics of confined water in different hydrated cement pastes with minimized contributions of capillary water. It was found that the water motions are extremely reduced compared to those of bulk water. The onset of water mobility, which was modified by the local environment, was investigated with elastic temperature scans using the high-resolution neutron backscattering instrument SPHERES. Using a Cauchy-Lorenz distribution, the quasi-elastic signal observed in the spectra obtained by the backscattering spectrometer was analyzed, leading to the identification of rotational motions with relaxation times of 0.3 ns. Additionally, neutron spin echo (NSE) spectroscopy was used to measure the water diffusion over the local network of pores. The motions observed in the NSE time scale were characterized by diffusion constants ranging from 0.6 to 1.1 x 10(-9) m(2) s(-1) most likely related to water molecules removed from the interface. In summary, our results indicate that the local diffusion observed in the gel pores of hardened cement pastes is on the order of that found in deeply supercooled water. Finally, the importance of the magnetic properties of cement pastes were discussed in relation to the observation of a quasi-elastic signal on the dried sample spectra measured using the time-of-flight spectrometer.

Noncooperative dimethyl sulfoxide-induced dissection of insulin fibrils: toward soluble building blocks of amyloid

The enormous molecular weight complicates detailed structural studies of amyloid fibrils and obscures identification of biologically active forms of protein aggregates in amyloid-related diseases. Here we show that aqueous solutions of dimethyl sulfoxide (DMSO) solubilize insulin fibrils while maintaining their beta-pleated structure. This is accompanied by a marked decrease in the fluorescence of thioflavin T. According to atomic force microscopy images and dynamic light scattering measurements, the partial DMSO-induced dissection of insulin fibrils favors formation of smaller soluble oligomers, which retain a limited capacity to induce daughter generation of fibrils through seeding to the native insulin, as well as the ability to reassemble into fibrils upon removal of DMSO through dialysis against water. These findings suggest that the DMSO-induced ensembles of insulin molecules are closely related to elementary building blocks of amyloid fibrils.

Dynamical properties of confined water nanoclusters: Simulation study of hydrated zeolite NaA: structural and vibrational properties

Water nanoclusters confined to zeolitic cavities have been extensively investigated by various experimental techniques. We report a series of molecular dynamics simulations at different temperatures and for water nanoclusters of different sizes in order to attempt an atomistic interpretation of the properties of these systems. The cavities of zeolite NaA are spherical in shape and about 1 nm in diameter and can host nanoclusters of water containing nearly up to 24 water molecules. A modified interaction potential, yielding a better reproduction of experimental hydration energy and water diffusivity across a number of different zeolites, is proposed. Molecular dynamics simulations reproduce the known experimental structural features obtained by X-ray diffraction. Variations of simulated vibrational IR and IINS spectra with temperature and size of nanoclusters are in good agreement with experiment. The simulated water nanoclusters in zeolite NaA are found to be too small to crystallize and, at low temperature, behave as amorphous ice, in agreement with recent experimental results for similar water nanoclusters in reverse micelles.

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

Vibrational Spectroscopy and Dynamics of Azide Ion in Ionic Liquid and Dimethyl Sulfoxide Water Mixtures

Steady-state and time-resolved infrared spectroscopy of the azide (N(3)-) anion has been used to characterize aqueous mixtures both with the ionic liquid (IL) 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF(4)]) and with dimethyl sulfoxide (DMSO). In the DMSO-water mixtures, two anion vibrational bands are observed for low water mole fractions (0 > X(w) > 0.25), which indicates a heterogeneous ion solvation environment. The band at 2000 cm(-1) observed for neat DMSO does not shift but decreases in amplitude as the amount of water is increased. Another band appears at slightly higher frequency at low X(w) (=0.05). As the amount of water is increased, this band shifts to higher frequency and becomes stronger and is attributed to azide with an increasing degree of hydration. At intermediate and high X(w), a single band is observed that shifts almost linearly with water mole fraction toward the bulk water value. The heterogeneity is evident from the infrared pump-probe studies in which the decay times depend on probe frequency at low mole fraction. For the azide spectra in IL-water mixtures, a single azide band is observed for each mole fraction mixture. The azide band shifts almost linearly with mole fraction, indicating nearly ideal mixing behavior. As with the DMSO-water mixtures, the time-resolved IR decay times are probe-frequency-dependent at low mole fraction, again indicating heterogeneous solvation. In both the DMSO and IL mixtures with water, the relaxation times are slower than would be expected from ideal mixing, suggesting that vibrational relaxation of azide is more sensitive than its vibrational frequency to the solvent structure. The results are discussed in terms of preferential solvation and the degree to which the azide shift and vibrational relaxation depend on the degree of water association in the mixtures.

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.

Study of the dynamical properties of water in disaccharide solutions

This work presents quasi-elastic neutron scattering (QENS) and neutron spin echo (NSE) findings on homologous disaccharide (i.e. trehalose, maltose and sucrose)/water solutions as a function of temperature. The dynamical properties of these systems are investigated by QENS, which, on the picosecond scale, allows for the characterisation of the diffusion of both solutes and solvent. On the other hand, NSE investigates the dynamics on the nanosecond scale, allowing for the relaxation times of the disaccharide/water systems to be evaluated. The experimental data highlight a strong slowing down of water in the presence of disaccharides. The whole set of findings indicates, therefore, a noticeable disaccharide-water interaction, which is more intense in the case of trehalose. This feature can justify its higher bioprotective effectiveness.