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.

Mapping the structural trends in zinc aluminosilicate glasses.

The structure of zinc aluminosilicate glasses with the composition (ZnO)x(Al2O3)y(SiO2)1-x-y, where 0 ≤ x < 1, 0 ≤ y < 1, and x + y < 1, was investigated over a wide composition range by combining neutron and high-energy x-ray diffraction with 27Al magic angle spinning nuclear magnetic resonance spectroscopy. The results were interpreted using an analytical model for the composition-dependent structure in which the zinc ions do not act as network formers. Four-coordinated aluminum atoms were found to be in the majority for all the investigated glasses, with five-coordinated aluminum atoms as the main minority species. Mean Al-O bond distances of 1.764(5) and 1.855(5) Å were obtained for the four- and five-coordinated aluminum atoms, respectively. The coordination environment of zinc was not observed to be invariant. Instead, it is dependent on whether zinc plays a predominantly network-modifying or charge-compensating role and, therefore, varies systematically with the glass composition. The Zn-O coordination number and bond distance were found to be 4.36(9) and 2.00(1) Å, respectively, for the network-modifying role vs 5.96(10) and 2.08(1) Å, respectively, for the charge-compensating role. The more open coordination environment of the charge-compensator is related to an enhanced probability of zinc finding bridging oxygen atoms as nearest-neighbors, reflecting a change in the connectivity of the glass network comprising four-coordinated silicon and aluminum atoms as the alumina content is increased.

Structure of amorphous materials in the NASICON system Na1+xTi2SixP3-xO12.

The structure of glasses in the sodium (Na) super-ionic conductor (NASICON) system Na1+xTi2SixP3-xO12withx= 0.8 andx= 1.0 was explored by combining neutron and high-energy x-ray diffraction with29Si,31P and23Na solid-state nuclear magnetic resonance (NMR) spectroscopy. The29Si magic angle spinning (MAS) NMR spectra reveal that the silica component remains fully polymerized in the form of Si4units, i.e. the silicon atoms are bound to four bridging oxygen atoms. The31P{23Na} rotational echo adiabatic passage double resonance (REAPDOR) NMR data suggest that the31P MAS NMR line shape originates from four-coordinated Pnunits, wheren= 1, 2 or 3 is the number of bridging oxygen atoms per phosphorus atom. These sites differ in their31P-23Na dipolar coupling strengths. The results support an intermediate range order scenario of a phosphosilicate mixed network-former glass in which the phosphate groups selectively attract the Na+modifier ions. Titanium takes a sub-octahedral coordination environment with a mean Ti-O coordination number of 5.17(4) forx= 0.8 and 4.86(4) forx= 1.0. A mismatch between the P-O and Si-O bond lengths of 8% is likely to inhibit the incorporation of silicon into the phosphorus sites of the NASICON crystal structure.

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.

Structure of molten NaCl and the decay of the pair-correlations.

The structure of molten NaCl is investigated by combining neutron and x-ray diffraction with molecular dynamics simulations that employed interaction potentials with either rigid or polarizable ions. Special attention is paid to the asymptotic decay of the pair-correlation functions, which is related to the small-k behavior of the partial structure factors, where k denotes the magnitude of the scattering vector. The rigid-ion approach gives access to an effective restricted primitive model in which the anion and cation have equal but opposite charges and are otherwise identical. For this model, the decay of the pair-correlation functions is in qualitative agreement with simple theory. The polarizable ion approach gives a good account of the diffraction results and yields thermodynamic parameters (density, isothermal compressibility, Debye screening length, and heat capacity) in accord with experiment. The longest decay length for the partial pair-distribution functions is a factor of ≃2.5 times greater than the nearest-neighbor distance. The results are commensurate with the decay lengths found for the effective restricted primitive model, which are much shorter than those found in experiments on concentrated electrolytes or ionic liquids using surface force apparatus.

Structural model for amorphous aluminosilicates.

An analytical model is developed for the composition-dependent structure of the amorphous aluminosilicate materials (M2O)x(Al2O3)y(SiO2)1-x-y and (MO)x(Al2O3)y(SiO2)1-x-y, where 0 ≤ x ≤ 1 and 0 ≤ y ≤ 1. The model is based on a simple set of reactions and contains a single adjustable parameter p (0 ≤ p ≤ 1). The latter is found from 27Al solid-state nuclear magnetic resonance (NMR) experiments in the regime where R = x/y ≥ 1, aided by new experiments on the magnesium and zinc aluminosilicate systems. The parameter p decreases linearly as the cation field strength of M+ or M2+ increases, as per the observation previously made for the degree of aluminum avoidance [Lee et al., J. Phys. Chem. C 120, 737 (2016)]. The results indicate that as the cation field strength increases, there are less fourfold coordinated aluminum atoms to contribute toward the glass network, and Al-O-Al bonds become more prevalent in a progressive breakdown of Loewenstein's aluminum avoidance rule. The model gives a good account of the composition-dependent fraction of non-bridging oxygen (NBO) atoms for R ≥ 1, as assessed from the results obtained from solid-state NMR experiments. An extension of the model to (M2O3)x(Al2O3)y(SiO2)1-x-y glasses leads, however, to an excess of NBO atoms, the proportion of which can be reduced by invoking network-forming fivefold coordinated Al atoms and/or oxygen triclusters. The model provides a benchmark for predicting the structure-related properties of aluminosilicate materials and a starting point for predicting the evolution in the structure of these materials under the extreme conditions encountered in the Earth's interior or in processes such as sharp-contact loading.

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

The structure of a concentrated solution of NaCl in D2O 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 NaCl0 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.

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 Na1+xAlxGe2-x(PO4)3 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) 27Al and 31P magic angle spinning (MAS) and 31P/23Na 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 27Al 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(6) and tetrahedral P(3) motifs are embedded in a matrix of tetrahedral Ge(4) 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(3) motifs are used to convert Ge(4) to Ge(6) units. The resultant P(4) motifs thereby provide the nucleation sites for crystal growth via a homogeneous nucleation mechanism.

Persistent homology in two-dimensional atomic networks.

The topology of two-dimensional network materials is investigated by persistent homology analysis. The constraint of two dimensions allows for a direct comparison of key persistent homology metrics (persistence diagrams, cycles, and Betti numbers) with more traditional metrics such as the ring-size distributions. Two different types of networks are employed in which the topology is manipulated systematically. In the first, comparatively rigid networks are generated for a triangle-raft model, which are representative of materials such as silica bilayers. In the second, more flexible networks are generated using a bond-switching algorithm, which are representative of materials such as graphene. Bands are identified in the persistence diagrams by reference to the length scales associated with distorted polygons. The triangle-raft models with the largest ordering allow specific bands Bn (n = 1, 2, 3, …) to be allocated to configurations of atoms separated by n bonds. The persistence diagrams for the more disordered network models also display bands albeit less pronounced. The persistent homology method thereby provides information on n-body correlations that is not accessible from structure factors or radial distribution functions. An analysis of the persistent cycles gives the primitive ring statistics, provided the level of disorder is not too large. The method also gives information on the regularity of rings that is unavailable from a ring-statistics analysis. The utility of the persistent homology method is demonstrated by its application to experimentally-obtained configurations of silica bilayers and graphene.

Structure of As-Se glasses by neutron diffraction with isotope substitution.

The method of neutron diffraction with selenium isotope substitution is used to measure the structure of glassy As0.30Se0.70, As0.35Se0.65, and As0.40Se0.60. The method delivers three difference functions for each sample in which the As-As, As-Se, or Se-Se correlations are eliminated. The measured coordination numbers are consistent with the "8-N" rule and show that the As0.30Se0.70 network is chemically ordered, a composition near to which there is a minimum in the fragility index and a boundary to the intermediate phase. Chemical ordering in glassy As0.35Se0.65 and As0.40Se0.60 is, however, broken by the appearance of As-As bonds, the fraction of which increases with the arsenic content of the glass. For the As0.40Se0.60 material, a substantial fraction of As-As and Se-Se defect pairs (∼11%) is frozen into the network structure on glass formation.

Many-body effects at the origin of structural transitions in B2O3.

The structural properties of glassy diboron trioxide, g-B2O3, are investigated from ambient to high pressure conditions using two types of atomic force-field models that account for many-body effects. These models are parameterized by a dipole- and force-fitting procedure of reference datasets created via first-principles calculations on a series of configurations. The predictions of the models are tested against experimental data, where particular attention is paid to the structural transitions in g-B2O3 that involve changes to both the short- and medium-range order. The models outperform those previously devised, where improvement originates from the incorporation of two key physical ingredients, namely, (i) the polarizability of the oxide ion and (ii) the ability of an oxide ion to change both size and shape in response to its coordination environment. The results highlight the importance of many-body effects for accurately modeling this challenging system.

Structure of semiconducting versus fast-ion conducting glasses in the Ag-Ge-Se system.

The transition from a semiconductor to a fast-ion conductor with increasing silver content along the Ag x (Ge0.25Se0.75)(100-x) tie line (0≤x≤25) was investigated on multiple length scales by employing a combination of electric force microscopy, X-ray diffraction, and neutron diffraction. The microscopy results show separation into silver-rich and silver-poor phases, where the Ag-rich phase percolates at the onset of fast-ion conductivity. The method of neutron diffraction with Ag isotope substitution was applied to the x=5 and x=25 compositions, and the results indicate an evolution in structure of the Ag-rich phase with change of composition. The Ag-Se nearest-neighbours are distributed about a distance of 2.64(1) Å, and the Ag-Se coordination number increases from 2.6(3) at x=5 to 3.3(2) at x=25. For x=25, the measured Ag-Ag partial pair-distribution function gives 1.9(2) Ag-Ag nearest-neighbours at a distance of 3.02(2) Å. The results show breakage of Se-Se homopolar bonds as silver is added to the Ge0.25Se0.75 base glass, and the limit of glass-formation at x≃28 coincides with an elimination of these bonds. A model is proposed for tracking the breakage of Se-Se homopolar bonds as silver is added to the base glass.

Structure of rare-earth chalcogenide glasses by neutron and x-ray diffraction.

The method of neutron diffraction with isomorphic substitution was used to measure the structure of the rare-earth chalcogenide glasses [Formula: see text](Ga2 X 3)0.33(GeX 2)0.60 with [Formula: see text] or Ce and [Formula: see text] or Se. X-ray diffraction was also used to measure the structure of the sulphide glass. The results are consistent with networks that are built from GeX 4 and GaX 4 tetrahedra, and give R-S and R-Se coordination numbers of 8.0(2) and 8.5(4), respectively. The minimum nearest-neighbour R-R distance associated with rare-earth clustering is discussed.

Structure of amorphous GeSe9 by neutron diffraction and first-principles molecular dynamics: Impact of trajectory sampling and size effects.

The structure of glassy GeSe9 was investigated by combining neutron diffraction with density-functional-theory-based first-principles molecular dynamics. In the simulations, three different models of N = 260 atoms were prepared by sampling three independent temporal trajectories, and the glass structures were found to be substantially different from those obtained for models in which smaller numbers of atoms or more rapid quench rates were employed. In particular, the overall network structure is based on Sen chains that are cross-linked by Ge(Se4)1/2 tetrahedra, where the latter are predominantly corner as opposed to edge sharing. The occurrence of a substantial proportion of Ge-Se-Se connections does not support a model in which the material is phase separated into Se-rich and GeSe2-rich domains. The appearance of a first-sharp diffraction peak in the Bhatia-Thornton concentration-concentration partial structure factor does, however, indicate a non-uniform distribution of the Ge-centered structural motifs on an intermediate length scale.

Networks under pressure: the development of in situ high-pressure neutron diffraction for glassy and liquid materials.

The pressure-driven collapse in the structure of network-forming materials will be considered in the gigapascal (GPa) regime, where the development of in situ high-pressure neutron diffraction has enabled this technique to obtain new structural information. The improvements to the neutron diffraction methodology are discussed, and the complementary nature of the results is illustrated by considering the pressure-driven structural transformations for several key network-forming materials that have also been investigated by using other experimental techniques such as x-ray diffraction, inelastic x-ray scattering, x-ray absorption spectroscopy and Raman spectroscopy. A starting point is provided by the pressure-driven network collapse of the prototypical network-forming oxide glasses B2O3, SiO2 and GeO2. Here, the combined results help to show that the coordination number of network-forming structural motifs in a wide range of glassy and liquid oxide materials can be rationalised in terms of the oxygen-packing fraction over an extensive pressure and temperature range. The pressure-driven network collapse of the prototypical chalcogenide glass GeSe2 is also considered where, as for the case of glassy GeO2, site-specific structural information is now available from the method of in situ high-pressure neutron diffraction with isotope substitution. The application of in situ high-pressure neutron diffraction to other structurally disordered network-forming materials is also summarised. In all of this work a key theme concerns the rich diversity in the mechanisms of network collapse, which drive the changes in physico-chemical properties of these materials. A more complete picture of the mechanisms is provided by molecular dynamics simulations using theoretical schemes that give a good account of the experimental results.

High-pressure transformation of SiO2 glass from a tetrahedral to an octahedral network: a joint approach using neutron diffraction and molecular dynamics.

A combination of in situ high-pressure neutron diffraction at pressures up to 17.5(5) GPa and molecular dynamics simulations employing a many-body interatomic potential model is used to investigate the structure of cold-compressed silica glass. The simulations give a good account of the neutron diffraction results and of existing x-ray diffraction results at pressures up to ~60 GPa. On the basis of the molecular dynamics results, an atomistic model for densification is proposed in which rings are "zipped" by a pairing of five- and/or sixfold coordinated Si sites. The model gives an accurate description for the dependence of the mean primitive ring size ⟨n⟩ on the mean Si-O coordination number, thereby linking a parameter that is sensitive to ordering on multiple length scales to a readily measurable parameter that describes the local coordination environment.

Identifying and characterising the different structural length scales in liquids and glasses: an experimental approach.

The structure of several network-forming liquids and glasses is considered, where a focus is placed on the detailed information that is made available by using the method of neutron diffraction with isotope substitution (NDIS). In the case of binary network glass-forming materials with the MX2 stoichiometry (e.g. GeO2, GeSe2, ZnCl2), two different length scales at distances greater than the nearest-neighbour distance manifest themselves by peaks in the measured diffraction patterns. The network properties are influenced by a competition between the ordering on these "intermediate" and "extended" length scales, which can be manipulated by changing the chemical identity of the atomic constituents or by varying state parameters such as the temperature and pressure. The extended-range ordering, which describes the decay of the pair-correlation functions at large-r, can be represented by making a pole analysis of the Ornstein-Zernike equations, an approach that can also be used to describe the large-r behaviour of the pair-correlation functions for liquid and amorphous metals where packing constraints are important. The first applications are then described of the NDIS method to measure the detailed structure of aerodynamically-levitated laser-heated droplets of "fragile" glass-forming liquid oxides (CaAl2O4 and CaSiO3) at high-temperatures (~2000 K) and the structure of a "strong" network-forming glass (GeO2) under pressures ranging from ambient to ~8 GPa. The high-temperature experiments show structural changes on multiple length scales when the oxides are vitrified. The high-pressure experiment offers insight into the density-driven mechanisms of network collapse in GeO2 glass, and parallels are drawn with the high-pressure behaviour of silica glass. Finally, the hydrogen-bonded network of water is considered, where the first application of the method of oxygen NDIS is used to measure the structures of light versus heavy water and a difference of approximately equal to 0.5% is found between the O-D and O-H intra-molecular bond lengths. The experimental data are best matched by using path integral molecular dynamics simulations with a flexible anharmonic water model, and the results support a competing quantum effects model for water in which its structural and dynamical properties are governed by an offset between intra-molecular and inter-molecular quantum contributions.

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.

Mechanisms of network collapse in GeO2 glass: high-pressure neutron diffraction with isotope substitution as arbitrator of competing models.

The structure of the network forming glass GeO(2) is investigated by making the first application of the method of in situ neutron diffraction with isotope substitution at pressures increasing from ambient to 8 GPa. Of the various models, the experimental results are in quantitative agreement only with molecular dynamics simulations made using interaction potentials that include dipole-polarization effects. When the reduced density ρ/ρ(0) > or approximately equal to 1.16, where ρ(0) is the value at ambient pressure, network collapse proceeds via an interplay between the predominance of distorted square pyramidal GeO(5) units versus octahedral GeO(6) units as they replace tetrahedral GeO(4) units. This replacement necessitates the formation of threefold coordinated oxygen atoms and leads to an increase with density in the number of small rings, where a preference is shown for sixfold rings when ρ/ρ(0) = 1 and fourfold rings when ρ/ρ(0) = 1.64.

Density-driven structural transformations in network forming glasses: a high-pressure neutron diffraction study of GeO2 glass up to 17.5 GPa.

The structure of GeO(2) glass was investigated at pressures up to 17.5(5) GPa using in situ time-of-flight neutron diffraction with a Paris-Edinburgh press employing sintered diamond anvils. A new methodology and data correction procedure were developed, enabling a reliable measurement of structure factors that are largely free from diamond Bragg peaks. Calibration curves, which are important for neutron diffraction work on disordered materials, were constructed for pressure as a function of applied load for both single and double toroid anvil geometries. The diffraction data are compared to new molecular-dynamics simulations made using transferrable interaction potentials that include dipole-polarization effects. The results, when taken together with those from other experimental methods, are consistent with four densification mechanisms. The first, at pressures up to approximately equal 5 GPa, is associated with a reorganization of GeO(4) units. The second, extending over the range from approximately equal 5 to 10 GPa, corresponds to a regime where GeO(4) units are replaced predominantly by GeO(5) units. In the third, as the pressure increases beyond ~10 GPa, appreciable concentrations of GeO(6) units begin to form and there is a decrease in the rate of change of the intermediate-range order as measured by the pressure dependence of the position of the first sharp diffraction peak. In the fourth, at about 30 GPa, the transformation to a predominantly octahedral glass is achieved and further densification proceeds via compression of the Ge-O bonds. The observed changes in the measured diffraction patterns for GeO(2) occur at similar dimensionless number densities to those found for SiO(2), indicating similar densification mechanisms for both glasses. This implies a regime from about 15 to 24 GPa where SiO(4) units are replaced predominantly by SiO(5) units, and a regime beyond ~24 GPa where appreciable concentrations of SiO(6) units begin to form.