Isotope effects in water as investigated by neutron diffraction and path integral molecular dynamics

Pressure dependent structure of amorphous magnesium aluminosilicates: The effect of replacing magnesia by alumina at the enstatite composition

The effect of replacing magnesia by alumina on the pressure-dependent structure of amorphous enstatite was investigated by applying in situ high-pressure neutron diffraction with magnesium isotope substitution to glassy (MgO)0.375(Al2O3)0.125(SiO2)0.5. The replacement leads to a factor of 2.4 increase in the rate-of-change of the Mg-O coordination number with pressure, which increases from 4.76(4) at ambient pressure to 6.51(4) at 8.2 GPa, and accompanies a larger probability of magnesium finding bridging oxygen atoms as nearest-neighbors. The Al-O coordination number increases from 4.17(7) to 5.24(8) over the same pressure interval at a rate that increases when the pressure is above ∼3.5 GPa. On recovering the glass to ambient conditions, the Mg-O and Al-O coordination numbers reduce to 5.32(4) and 4.42(6), respectively. The Al-O value is in accordance with the results from solid-state 27Al nuclear magnetic resonance spectroscopy, which show the presence of six-coordinated aluminum species that are absent in the uncompressed material. These findings explain the appearance of distinct pressure-dependent structural transformation regimes in the preparation of permanently densified magnesium aluminosilicate glasses. They also indicate an anomalous minimum in the pressure dependence of the bulk modulus with an onset that suggests a pressure-dependent threshold for transitioning between scratch-resistant and crack-resistant material properties.

Structure of diopside, enstatite, and magnesium aluminosilicate glasses: A joint approach using neutron and x-ray diffraction and solid-state NMR

Neutron diffraction with magnesium isotope substitution, high energy x-ray diffraction, and 29Si, 27Al, and 25Mg solid-state nuclear magnetic resonance (NMR) spectroscopy were used to measure the structure of glassy diopside (CaMgSi2O6), enstatite (MgSiO3), and four (MgO) x(Al2O3) y(SiO2)1−x−y glasses, with x = 0.375 or 0.25 along the 50 mol. % silica tie-line (1 − x − y = 0.5) or with x = 0.3 or 0.2 along the 60 mol. % silica tie-line (1 − x − y = 0.6). The bound coherent neutron scattering length of the isotope 25Mg was remeasured, and the value of 3.720(12) fm was obtained from a Rietveld refinement of the powder diffraction patterns measured for crystalline 25MgO. The diffraction results for the glasses show a broad asymmetric distribution of Mg–O nearest-neighbors with a coordination number of 4.40(4) and 4.46(4) for the diopside and enstatite glasses, respectively. As magnesia is replaced by alumina along a tie-line with 50 or 60 mol. % silica, the Mg–O coordination number increases with the weighted bond distance as less Mg2+ ions adopt a network-modifying role and more of these ions adopt a predominantly charge-compensating role. 25Mg magic angle spinning (MAS) NMR results could not resolve the different coordination environments of Mg2+ under the employed field strength (14.1 T) and spinning rate (20 kHz). The results emphasize the power of neutron diffraction with isotope substitution to provide unambiguous site-specific information on the coordination environment of magnesium in disordered materials.

Long-Range Structures of Amorphous Solid Water

High-energy X-ray diffraction (XRD) and Fourier transform infrared spectroscopy (FTIR) of amorphous solid water (ASW) were studied during vapor deposition and the heating process. From the diffraction patterns, the oxygen–oxygen pair distribution functions (PDFs) were calculated up to the eighth coordination shell and an r = 23 Å. The PDF of ASW obtained both during vapor deposition at 80 K as well as the subsequent heating are consistent with that of low-density amorphous ice. The formation and temperature-induced collapse of micropores were observed in the XRD data and in the FTIR measurements, more specifically, in the OH stretch and the dangling mode. Above 140 K, ASW crystallizes into a stacking disordered ice, I_sd. It is observed that the fourth, fifth, and sixth peaks in the PDF, corresponding to structural arrangements between 8 and 12 Å, are the most sensitive to the onset of crystallization.

Structure and dynamics of aqueous NaCl solutions at high temperatures and pressures

The structure of a concentrated solution of NaCl in D₂O was investigated by in situ high-pressure neutron diffraction with chlorine isotope substitution to give site-specific information on the coordination environment of the chloride ion. A broad range of densities was explored by first increasing the temperature from 323 to 423 K at 0.1 kbar and then increasing the pressure from 0.1 to 33.8 kbar at 423 K, thus mapping a cyclic variation in the static dielectric constant of the pure solvent. The experimental work was complemented by molecular dynamics simulations using the TIP4P/2005 model for water, which were validated against the measured equation of state and diffraction results. Pressure-induced anion ordering is observed, which is accompanied by a dramatic increase in the Cl-O and O-O coordination numbers. With the aid of bond-distance resolved bond-angle maps, it is found that the increased coordination numbers do not originate from a sizable alteration to the number of either Cl⋯D-O or O⋯D-O hydrogen bonds but from the appearance of non-hydrogen-bonded configurations. Increased pressure leads to a marked decrease in the self-diffusion coefficients but has only a moderate effect on the ion-water residence times. Contact ion pairs are observed under all conditions, mostly in the form of charge-neutral NaCl⁰ units, and coexist with solvent-separated Na⁺-Na⁺ and Cl^--Cl^- ion pairs. The exchange of water molecules with Na⁺ adopts a concerted mechanism under ambient conditions but becomes non-concerted as the state conditions are changed. Our findings are important for understanding the role of extreme conditions in geochemical processes.

Modeling Liquid Water by Climbing up Jacob's Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn-Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob's ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric-isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

Experimental Determination of Relationship between Intramolecular O-D Bond Length and Its Stretching Vibrational Frequency of D2O Molecule in the Liquid State

Experimental evidence has been obtained for the structure-spectra relationship of hydrogen bonds in aqueous solutions. Intramolecular O-D distance, r_OD, has been determined by the least-squares fitting analysis of the neutron interference term in the high-Q region observed for pure D₂O and concentrated aqueous solutions. The average O-D stretching frequency, ν_OD, has been obtained from the position of the center of gravity of the observed ATR-IR O-D stretching band. The linear relationship between r_OD and ν_OD has been confirmed in the liquid state. The slope of dν_OD/dr_OD is evaluated to be -21 000 ± 1000 cm^-1 Å^-1.

Structure of crystalline and amorphous materials in the NASICON system Na1+xAlxGe2-x(PO4)3

The structure of crystalline and amorphous materials in the sodium (Na) super-ionic conductor system Na_1+xAl_xGe_2-x(PO₄)₃ with x = 0, 0.4, and 0.8 was investigated by combining (i) neutron and x-ray powder diffraction and pair-distribution function analysis with (ii) ²⁷Al and ³¹P magic angle spinning (MAS) and ³¹P/²³Na double-resonance nuclear magnetic resonance (NMR) spectroscopy. A Rietveld analysis of the powder diffraction patterns shows that the x = 0 and x = 0.4 compositions crystallize into space group-type R3̄, whereas the x = 0.8 composition crystallizes into space group-type R3̄c. For the as-prepared glass, the pair-distribution functions and ²⁷Al MAS NMR spectra show the formation of sub-octahedral Ge and Al centered units, which leads to the creation of non-bridging oxygen (NBO) atoms. The influence of these atoms on the ion mobility is discussed. When the as-prepared glass is relaxed by thermal annealing, there is an increase in the Ge and Al coordination numbers that leads to a decrease in the fraction of NBO atoms. A model is proposed for the x = 0 glass in which super-structural units containing octahedral Ge⁽⁶⁾ and tetrahedral P⁽³⁾ motifs are embedded in a matrix of tetrahedral Ge⁽⁴⁾ units, where superscripts denote the number of bridging oxygen atoms. The super-structural units can grow in size by a reaction in which NBO atoms on the P⁽³⁾ motifs are used to convert Ge⁽⁴⁾ to Ge⁽⁶⁾ units. The resultant P⁽⁴⁾ motifs thereby provide the nucleation sites for crystal growth via a homogeneous nucleation mechanism.

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF-Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born-Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. Here, using tetrahydrofuran (THF)-water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D₈THF to examine the molecular origins of the isotope-dependent loss of the THF-water miscibility gap.

Opposing Electronic and Nuclear Quantum Effects on Hydrogen Bonds in H2 O and D2 O

The effect of extending the O-H bond length(s) in water on the hydrogen-bonding strength has been investigated using static ab initio molecular orbital calculations. The "polar flattening" effect that causes a slight σ-hole to form on hydrogen atoms is strengthened when the bond is stretched, so that the σ-hole becomes more positive and hydrogen bonding stronger. In opposition to this electronic effect, path-integral ab initio molecular-dynamics simulations show that the nuclear quantum effect weakens the hydrogen bond in the water dimer. Thus, static electronic effects strengthen the hydrogen bond in H2 O relative to D2 O, whereas nuclear quantum effects weaken it. These quantum fluctuations are stronger for the water dimer than in bulk water.

X-ray Scattering and O-O Pair-Distribution Functions of Amorphous Ices

The structure factor and oxygen–oxygen pair-distribution functions of amorphous ices at liquid nitrogen temperature (T = 77 K) have been derived from wide-angle X-ray scattering (WAXS) up to interatomic distances of r = 23 Å, where local structure differences between the amorphous ices can be seen for the entire range. The distances to the first coordination shell for low-, high-, and very-high-density amorphous ice (LDA, HDA, VHDA) were determined to be 2.75, 2.78, and 2.80 Å, respectively, with high accuracy due to measurements up to a large momentum transfer of 23 Å^–1. Similarities in pair-distribution functions between LDA and supercooled water at 254.1 K, HDA and liquid water at 365.9 K, and VHDA and high-pressure liquid water were found up to around 8 Å, but beyond that at longer distances, the similarities were lost. In addition, the structure of the high-density amorphous ices was compared to high-pressure crystalline ices IV, IX , and XII, and conclusions were drawn about the local ordering.

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.

Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs

Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.

Family of Oxygen-Oxygen Radial Distribution Functions for Water

In a typical X-ray diffraction experiment, the elastically scattered intensity, I(Q), is the experimental observable. I(Q) contains contributions from both intramolecular as well as intermolecular correlations embodied in the scattering factors, HOO(Q) and HOH(Q), with negligible contributions from HHH(Q). Thus, to accurately define the oxygen-oxygen radial distribution function, gOO(r), a model of the electron density is required to accurately weigh the HOO(Q) component relative to the intramolecular and oxygen-hydrogen correlations from the total intensity observable. In this work, we carefully define the electron density model and its underlying assumptions and more explicitly utilize two restraints on the allowable gOO(r) functions, which must conform to both very low experimental errors at high Q and the need to satisfy the isothermal compressibility at low Q. Although highly restrained by these conditions, the underdetermined nature of the problem is such that we present a family of gOO(r) values that provide equally good agreement with the high-Q intensity and compressibility restraints and with physically correct behavior at small r.

The structure of water around the compressibility minimum

Here we present diffraction data that yield the oxygen-oxygen pair distribution function, g(OO)(r) over the range 254.2-365.9 K. The running O-O coordination number, which represents the integral of the pair distribution function as a function of radial distance, is found to exhibit an isosbestic point at 3.30(5) Å. The probability of finding an oxygen atom surrounding another oxygen at this distance is therefore shown to be independent of temperature and corresponds to an O-O coordination number of 4.3(2). Moreover, the experimental data also show a continuous transition associated with the second peak position in g(OO)(r) concomitant with the compressibility minimum at 319 K.

Strong isotope effects on melting dynamics and ice crystallisation processes in cryo vitrification solutions

The nucleation and growth of crystalline ice during cooling, and further crystallization processes during re-warming are considered to be key processes determining the success of low temperature storage of biological objects, as used in medical, agricultural and nature conservation applications. To avoid these problems a method, termed vitrification, is being developed to inhibit ice formation by use of high concentration of cryoprotectants and ultra-rapid cooling, but this is only successful across a limited number of biological objects and in small volume applications. This study explores physical processes of ice crystal formation in a model cryoprotective solution used previously in trials on vitrification of complex biological systems, to improve our understanding of the process and identify limiting biophysical factors. Here we present results of neutron scattering experiments which show that even if ice crystal formation has been suppressed during quench cooling, the water molecules, mobilised during warming, can crystallise as detectable ice. The crystallisation happens right after melting of the glass phase formed during quench cooling, whilst the sample is still transiting deep cryogenic temperatures. We also observe strong water isotope effects on ice crystallisation processes in the cryoprotectant mixture. In the neutron scattering experiment with a fully protiated water component, we observe ready crystallisation occurring just after the glass melting transition. On the contrary with a fully deuteriated water component, the process of crystallisation is either completely or substantially supressed. This behaviour might be explained by nuclear quantum effects in water. The strong isotope effect, observed here, may play an important role in development of new cryopreservation strategies.

Anomalously large isotope effect in the glass transition of water

We present the discovery of an unusually large isotope effect in the structural relaxation and the glass transition temperature Tg of water. Dielectric relaxation spectroscopy of low-density as well as of vapor-deposited amorphous water reveal Tg differences of 10 ± 2 K between H2O and D2O, sharply contrasting with other hydrogen-bonded liquids for which H/D exchange increases Tg by typically less than 1 K. We show that the large isotope effect and the unusual variation of relaxation times in water at low temperatures can be explained in terms of quantum effects. Thus, our findings shed new light on water's peculiar low-temperature dynamics and the possible role of quantum effects in its structural relaxation, and possibly in dynamics of other low-molecular-weight liquids.

Mode-distribution analysis of quasielastic neutron scattering and application to liquid water

A quasielastic neutron scattering (QENS) experiment is a particular technique that endeavors to define a relationship between time and space for the diffusion dynamics of atoms and molecules. However, in most cases, analyses of QENS data are model dependent, which may distort attempts to elucidate the actual diffusion dynamics. We have developed a method for processing QENS data without a specific model, wherein all modes can be described as combinations of the relaxations based on the exponential law. By this method, we can obtain a distribution function B(Q,Γ), which we call the mode-distribution function (MDF), to represent the number of relaxation modes and distributions of the relaxation times in the modes. The deduction of MDF is based on the maximum entropy method and is very versatile in QENS data analysis. To verify this method, reproducibility was checked against several analytical models, such as that with a mode of distributed relaxation time, that with two modes closely located, and that represented by the Kohlrausch-Williams-Watts function. We report the first application to experimental data of liquid water. In addition to the two known modes, the existence of a relaxation mode of water molecules with an intermediate time scale has been discovered. We propose that the fast mode might be assigned to an intermolecular motion and the intermediate motion might be assigned to a rotational motion of the water molecules instead of to the fast mode.

Hyperfine coupling of the hydrogen atom in high temperature water

The hyperfine coupling constant of the hydrogen atom has been measured in pressurized liquid water up to 300 °C. The reduced constant A(water)∕A(vacuum) is 0.9939 at room temperature, and decreases to a minimum of 0.9918 at 240 °C. The reduced constant then increases at higher temperature. The g-factor is 2.002244(10) at room temperature and decreases to 2.00221(1) at 240 °C. The change in g-factor is proportional to the change in hyperfine coupling. The behavior below 110 °C is in excellent agreement with a previously proposed model in which the H atom is confined to a harmonic solvent cage, and vibrations within the cage mix "p-type" character into the wavefunction, resulting inA(water)∕A(vacuum) < 1. The harmonic model breaks down above 130 °C. We demonstrate that a classical binary collision model using approximate partial molar volume information can recover the observed minima near 240 °C.

The bound coherent neutron scattering lengths of the oxygen isotopes

The technique of neutron interferometry was used to measure the bound coherent neutron scattering length b(coh) of the oxygen isotopes (17)O and (18)O. From the measured difference in optical path between two water samples, either H(2)(17)O or H(2)(18)O versus H(2)(nat)O, where nat denotes the natural isotopic composition, we obtain b(coh,(17)O) = 5.867(4) fm and b(coh,(18)O) = 6.009(5) fm, based on the accurately known value of b(coh,(nat)O) = 5.805(4) fm which is equal to b(coh,(16)O) within the experimental uncertainty. Our results for b(coh,(17)O) and b(coh,(18)O) differ appreciably from the standard tabulated values of 5.6(5) fm and 5.84(7) fm, respectively. In particular, our measured scattering-length contrast of 0.204(3) fm between (18)O and (nat)O is nearly a factor of 6 greater than the tabulated value, which renders feasible neutron diffraction experiments using (18)O isotope substitution and thereby offers new possibilities for measuring the partial structure factors of oxygen-containing compounds, such as water.