Is water one liquid or two?

The idea that water is a mixture of two distinct states is analyzed in some detail. It is shown that the known compressibility of water is in fact sufficiently small that for a volume of water of size 1 nm3, the density fluctuations are of order 4% of the average density. This is much smaller than the ≈25% density fluctuations that would be required for significant regions of high and low density water to occur on this volume scale. It is also pointed out that the density fluctuations in water are, if anything, smaller than those that occur in other common liquids which do not have the anomalous properties of water. It is shown that if the distribution of density fluctuations is unimodal, the system is in the one-phase region, and if bimodal, it is in the two-phase region. None of the liquid or amorphous phases of water explored in this work give any sign of being in the two-phase region. Existing neutron and X-ray scattering data on water in the amorphous phases, and in the stable liquid phases as a function pressure and temperature, are subject to a new set of empirical potential structure refinement simulations. These simulations are interrogated for their configurational entropy, using a spherical harmonic reconstruction of the full orientational pair correlation function. It is shown that the excess pair entropy derived from this function, plus the known perfect gas contributions, give a reasonable account of the total entropy of water, within the likely errors. This estimated entropy follows the expected declining trend with decreasing temperature. Evidence that higher density water will have higher entropy than lower density water emerges, in accordance with what is expected from the negative thermal expansion coefficient of water at low temperatures. However, this entropy increase is not large and goes through a maximum before declining at yet higher densities and pressures, in a manner reminiscent of what has been previously observed in the diffusion coefficient as a function of pressure. There is no evidence that ambient water can be regarded as patches of high density, high entropy and low density, low entropy liquid, as some have claimed, since high density water has a similar entropy to low density water. There is some evidence that the distinction between these two states will become more pronounced as the temperature is lowered. Extensive discussion of the use of order parameters to describe water structure is given, and it is pointed out that these indices generally cannot be used to infer two-state behavior.

Ice crystallization observed in highly supercooled confined water.

We investigate the state of water confined in the cylindrical pores of MCM-41 type mesoporous silica, with pore diameters of 2.8 nm and 4.5 nm, over the temperature range 160-290 K by combining small angle neutron scattering and wide angle diffraction. This allows us to observe simultaneously the intermolecular correlations in the local water structure (which shows up in a main water peak around Q = 1.7 Å-1), the two-dimensional hexagonal arrangement of water cylinders in the silica matrix (which gives rise to a pronounced Bragg peak around Q = 0.2 Å-1), and the so-called Porod scattering at smaller Q, which arises from larger scale interfacial scattering within the material. In the literature, the temperature evolution of the intensity of the Bragg peak has been interpreted as the signature of a density minimum in confined water at approximately 210 K. Here we show that, under the conditions of our experiment, a fraction of freezable water coexists with a layer of non-freezable water within the pore volume. The overall temperature dependence of our data in the different Q regions, as well as the comparison of the data for the two pore sizes, leads us to conclude that the observed variation in the intensity of the Bragg diffraction peak is actually caused by a liquid to ice transition in the freezable fraction of confined water.

Coarse-grained empirical potential structure refinement: Application to a reverse aqueous micelle.

Conventional atomistic computer simulations, involving perhaps up to 106atoms, can achieve length-scales on the order of a few 10s of nm. Yet many heterogeneous systems, such as colloids, nano-structured materials, or biological systems, can involve correlations over distances up 100s of nm, perhaps even 1µm in some instances. For such systems it is necessary to invoke coarse-graining, where single atoms are replaced by agglomerations of atoms, usually represented as spheres, in order for the simulation to be performed within a practical computer memory and time scale. Small angle scattering and reflectivity measurements, both X-ray and neutron, are routinely used to investigate structure in these systems, and traditionally the data have been interpreted in terms of discrete objects, such as spheres, sheets, and cylinders, and combinations thereof. Here we combine the coarse-grained computer simulation approach with neutron small angle scattering to refine the structure of a heterogeneous system, in the present case a reverse aqueous micelle of sodium-dioctyl sulfosuccinate (AOT) and iso-octane. The method closely follows empirical potential structure refinement and involves deriving an empirical interaction potential from the scattering data. As in traditional coarse-grained methods, individual atoms are replaced by spherical density profiles, which, unlike real atoms, can inter-penetrate to a significant extent. The method works over an arbitrary range of length-scales, but is limited to around 2 orders of magnitude in distance above a specified dimension. The smallest value for this dimension is of order 1nm, but the largest dimension is arbitrary. This article is part of a Special Issue entitled "Recent Advances in Bionanomaterials" Guest Editor: Dr. Marie-Louise Saboungi and Dr. Samuel D. Bader.

Segregated water observed in a putative fish embryo cryopreservative.

Development of new cryopreservation strategies has major potential in medicine and agriculture and is critical to the conservation of endangered species that currently cannot be preserved. A critical property of any potential cryopreservative solution is its ability to prevent cell-damaging ice formation during cooling and subsequent heating. This study focuses on the freezing behaviour of promising model cryoprotective solutions. We perform neutron scattering analysis, combined with computer modelling, of the water structure after quench cooling these solutions. It is found that water in this solution forms nano-clusters encapsulated by the surrounding matrix of cryoprotectant solute molecules. We posit that these small volumes inhibit ice formation, because water does not have space for the structural relaxation required to crystallize on the timescale of the cooling process.

Low-Density Water Structure Observed in a Nanosegregated Cryoprotectant Solution at Low Temperatures from 285 to 238 K.

The structure of liquid water is defined by its molecular association through hydrogen bonding. Two different structures have been proposed for liquid water at low temperatures: low-density liquid (LDL) and high-density liquid (HDL) water. Here, we demonstrate a platform that can be exploited to experimentally probe the structure of liquid water in equilibrium at temperatures down to 238 K. We make use of a cryoprotectant molecule, glycerol, that, when mixed with water, lowers the freezing temperature of the solution nonmonotonically with glycerol concentration. We use a combination of neutron diffraction measurements and computational modeling to examine the structure of water in glycerol-water liquid mixtures at low temperatures from 285 to 238 K. We confirm that the mixtures are nanosegregated into regions of glycerol-rich and water-rich clusters. We examine the water structure and reveal that, at the temperatures studied here, water forms a low-density water structure that is more tetrahedral than the structure at room temperature. We postulate that nanosegregation allows water to form a low-density structure that is protected by an extensive and encapsulating glycerol interface.

Polar stacking of molecules in liquid chloroform.

Using neutron diffraction and the isotopic substitution technique we have investigated the local structure of liquid chloroform. A strong tendency for polar stacking of molecules with collinear alignment of dipole moments is found. We speculate that these polar stacks contribute to the performance of chloroform as a solvent.

Neutron diffraction study of aqueous Laponite suspensions at the NIMROD diffractometer.

The process of dynamical arrest, leading to formation of different arrested states such as glasses and gels, along with the closely related process of aging, is central for both basic research and technology. Here we report on a study of the time-dependent structural evolution of two aqueous Laponite clay suspensions at different weight concentrations. Neutron diffraction experiments have been performed with the near and intermediate range order diffractometer (NIMROD) that allows studies of the structure of liquids and disordered materials over a continuous length scale ranging from 1 to 300 Å, i.e., from the atomistic to the mesoscopic scales. NIMROD is presently a unique diffractometer, bridging the length scales traditionally investigated by small angle neutron scattering or small angle x-ray scattering with that accessible by traditional diffractometers for liquids. Interestingly, we have unveiled a signature of aging of both suspensions in the length scale region of NIMROD. This phenomenon, ascribed to sporadic contacts between Laponite platelets at long times, has been observed with the sample arrested as gel or as repulsive glass. Moreover, water molecules within the layers closest to Laponite platelets surface show orientational and translational order, which maps into the crystalline structure of Laponite.

Empirical potential structure refinement of semi-crystalline polymer systems: polytetrafluoroethylene and polychlorotrifluoroethylene.

Empirical potential structure refinement (EPSR) simulations are performed on total neutron scattering data from powder samples of polytetrafluoroethylene (PTFE) and polychlorotrifluoroethylene (PCTFE), both at 300 K. Starting from single strands of polymer consisting of between 30 and 60 monomers of tetrafluoroethylene and chlorotrifluoroethylene in each case, hexagonal simulation cells are constructed consisting of an array 25 (5×5) such strands placed on a hexagonal lattice. Allowed simulation moves are polymer translation moves along all three Cartesian axes, whole polymer rotations about the polymer axis, and individual atom moves within each polymer. For PTFE a number of Bragg peaks are visible in the scattering data and these are found to be consistent with a lattice spacing a(=b) = 5.69(1) Å with a dihedral angle along the (helical) chain of 166° which gives a repeat distance along the chain (c-axis) of ~19.6 Å. The positions of the Bragg peaks are well reproduced by this model, although there is a mismatch in the amplitudes of some of the higher order reflections between simulation and data. For PCTFE there is only one visible Bragg peak (100) which is well reproduced by a hexagonal lattice of atactic parallel polymers with a spacing of a(=b) = 6.37(1) Å. In this case the absence of distinct reflections along the polymer c-axis makes characterization of the internal dihedral angle difficult, but a model with a dihedral angle of 166° was less successful at fitting the diffuse scattering than a model where this angle was set to 180°, giving a nearly straight trans (zig-zag) structure. For PCTFE little change in structure could be discerned when the material was heated to 550 K, apart from a slight increase in lattice spacing. In both cases there is substantial diffuse scattering between the Bragg peaks, and this is correctly replicated by the EPSR simulations.

Molecular insight into the hydrogen bonding and micro-segregation of a cryoprotectant molecule.

Glycerol-water liquid mixtures are intriguing hydrogen-bonded systems and essential in many fields of chemistry, ranging from basic molecular research to widespread use in industrial and biomedical applications as cryoprotective solutions. Despite much research on these mixtures, the details of their microscopic structure are still not understood. One common notion is that glycerol acts to diminish the hydrogen bonding ability of water, a recurring hypothesis that remains untested by direct experimental approaches. The present work characterizes the structure of glycerol-water mixtures, across the concentration range, using a combination of neutron diffraction experiments and computational modeling. Contrary to previous expectations, we show that the hydrogen bonding ability of water is not diminished in the presence of glycerol. We show that glycerol-water hydrogen bonds effectively take the place of water-water hydrogen bonds, allowing water to maintain its full hydrogen bonding capacity regardless of the quantity of glycerol in the environment. We provide a quantitative measurement of all hydrogen bonding in the system and reveal a concentration range where a microsegregated, bipercolating liquid mixture exists in coexistence with a considerable interface region. This work highlights the role of hydrogen bonding connectivity rather than water structuring/destructuring effects in these important cryoprotective systems.

The hydrogen-bonding ability of the amino acid glutamine revealed by neutron diffraction experiments.

Hydrogen bonding between glutamine residues has been identified as playing an important role in the intermolecular association and aggregation of proteins. To establish the molecular mechanisms of glutamine interactions, neutron diffraction coupled with hydrogen/deuterium isotopic substitution in combination with computational modeling has been used to investigate the structure and hydration of glutamine in aqueous solution. The final structures obtained are consistent with the experimental data and provide insight into the hydrogen-bonding ability of glutamine. We find that the backbone of glutamine is able to coordinate more water molecules than the side chain, suggesting that charged groups on the glutamine molecule are more successful in attracting water than the dipole in the side chain. In both the backbone and the side chain, we find that the carbonyl groups interact more readily with water molecules than the amine groups. We find that glutamine-glutamine interactions are present, despite their low concentration in this dilute solution. This is evidenced through the occurrence of dimers of glutamine molecules in the solution, demonstrating the effective propensity of this molecule to associate through backbone-backbone, backbone-side chain, and side chain-side chain hydrogen bond interactions. The formation of dimers of glutamine molecules in such a dilute solution (30 mg/mL glutamine) may have implications in the aggregation of glutamine-rich proteins in neurological diseases where aggregation is prevalent.

Preference for isolated water molecules in a concentrated glycerol-water mixture.

Neutron diffraction coupled with hydrogen/deuterium isotopic substitution has been used to investigate the structure of a concentrated glycerol water (4:1 mole fraction) solution. The neutron diffraction data were used to constrain a three-dimensional computational model that is experimentally relevant using the empirical potential structure refinement technique. From interrogation of this model, we find that glycerol-glycerol hydrogen bonding is largely unperturbed by the presence of water in the solution. We find that glycerol-water hydrogen bonding is prevalent, suggesting that water molecules effectively take the place of glycerol molecules in this concentrated solution. In contrast, we find that water-water hydrogen bonding is significantly perturbed. While the first coordination shell of water in the concentrated solution remains similar to that of pure water, water-water hydrogen bonding is greatly diminished beyond the first neighbor distance. Interestingly, the majority of water molecules exist as single monomers in the concentrated glycerol solution. The preference of isolated water molecules results in a solution that is well mixed with optimal glycerol-water hydrogen bonding. These results highlight the importance of preferential hydrogen bonding in aqueous solutions and suggest a mechanism for cryoprotection by which glycerol effectively hydrogen bonds with water, resulting in a disrupted hydrogen-bonded water network.

The structure of glycerol in the liquid state: a neutron diffraction study.

Neutron diffraction coupled with hydrogen/deuterium isotopic substitution has been used to investigate the structure of the pure cryoprotectant glycerol in the liquid state at 298 K and 1 atm. The neutron diffraction data were used to constrain a 3 dimensional computational model that is experimentally relevant using the empirical potential structure refinement (EPSR) technique. These simulations lead to a model structure of the glycerol molecule that is consistent with the experimental data. Interestingly, from interrogation of this structure, it is found that the number of hydrogen bonds per molecule is larger than had previously been suggested. Furthermore, converse to previous work, no evidence for intra-molecular hydrogen bonds is found. These results highlight the importance and relevance of using experimental data to inform computational modelling of even simple liquid systems.

NIMROD: The Near and InterMediate Range Order Diffractometer of the ISIS second target station.

NIMROD is the Near and InterMediate Range Order Diffractometer of the ISIS second target station. Its design is optimized for structural studies of disordered materials and liquids on a continuous length scale that extends from the atomic, upward of 30 nm, while maintaining subatomic distance resolution. This capability is achieved by matching a low and wider angle array of high efficiency neutron scintillation detectors to the broad band-pass radiation delivered by a hybrid liquid water and liquid hydrogen neutron moderator assembly. The capabilities of the instrument bridge the gap between conventional small angle neutron scattering and wide angle diffraction through the use of a common calibration procedure for the entire length scale. This allows the instrument to obtain information on nanoscale systems and processes that are quantitatively linked to the local atomic and molecular order of the materials under investigation.

Water and trehalose: how much do they interact with each other?

The observation made by early naturalists that some organisms could tolerate extreme environmental condisions and "enjoy the advantage of real resurrection after death" [ Spallanzani , M. Opuscules de Physique Animale et Vegetale 1776 (translated from Italian by Senebier , J. Opuscules de Physique Animale et Vegetale 1787 , 2 , 203 - 285 )] stimulated research that still continues to this day. Cryptobiosis, the ability of an organism to tolerate adverse environments, such as dehydration and low temperatures, still represents an unsolved and fascinating problem. It has been shown that many sugars play an important role as bioprotectant agents, and among the best performers is the disaccharide trehalose. The current hypothesis links the efficiency of its protective role to strong modifications of the tetrahedral arrangement of water molecules in the sugar hydration shell, with trehalose forming many hydrogen bonds with the solvent. Here, we show, by means of state-of-the-art neutron diffraction experiments combined with EPSR simulations, that trehalose solvation induces very minor modifications of the water structure. Moreover, the number of water molecules hydrogen-bonded to the sugar is surprisingly small.

Multiscale approach to the structural study of water confined in MCM41.

We present a protocol for simultaneous structural characterization of a confined fluid and the confining substrate, along with the extraction of site-site pair correlation functions of the liquid of interest. This is based on neutron diffraction experiments, exploiting where feasible the isotopic substitution technique, analyzed through numerical coarse graining calculations and atomistic simulations. All of the subtleties of the experimental procedure, the needed ancillary measurements, and the recipe for tailoring the numerical codes to the real experiment and sample are described in the case of water confined in MCM41-S-15. In particular the excluded volume effects and the relevance of liquid-substrate cross-correlation terms in the neutron cross section are quantitatively discussed. The results obtained for the microscopic structure of water evidence a non-homogeneous distribution of molecules within the pore, with the presence of water-substrate hydrogen bonds, and a strong distortion of the water-water radial distribution functions with respect to those of bulk water extending at least up to three hydration layers.

Influence of concentration and anion size on hydration of H+ ions and water structure.

Neutron diffraction experiments with hydrogen isotope substitution on aqueous solutions of HCl and HBr have been performed at concentrations ranging from 1:17 to 1:83 solute per water molecules, at ambient conditions. Data are analyzed using the empirical potential structure refinement technique in order to extract information on both the ion hydration shells and the microscopic structure of the solvent. It is found that the influence of these solutes on the water structure is less concentration dependent than that of salts or hydroxides. Moreover protons readily form a strong H-bond with a water molecule upon solvation, at all proportions. The majority of them is also bonded via a longer bond to another water molecule, giving a prepeak in the g(OwOw). At high solute concentration, the second water molecule may be substituted by the counterion. In particular at solute concentrations of the order of 1:17 or higher, all protons have an anion within a distance of 4.5 A.

Percolation and three-dimensional structure of supercritical water.

It is well established that at ambient and supercooled conditions water can be described as a percolating network of H bonds. This work is aimed at identifying, by neutron diffraction experiments combined with computer simulations, a percolation line in supercritical water, where the extension of the H-bond network is in question. It is found that in real supercritical water liquidlike states are observed at or above the percolation threshold, while below this threshold gaslike water forms small, sheetlike configurations. Inspection of the three-dimensional arrangement of water molecules suggests that crossing of this percolation line is accompanied by a change of symmetry in the first neighboring shell of molecules from trigonal below the line to tetrahedral above.