Plutonium oxide melt structure and covalency.

Advances in nuclear power reactors include the use of mixed oxide fuel, containing uranium and plutonium oxides. The high-temperature behaviour and structure of PuO2-x above 1,800 K remain largely unexplored, and these conditions must be considered for reactor design and planning for the mitigation of severe accidents. Here, we measure the atomic structure of PuO2-x through the melting transition up to 3,000 ± 50 K using X-ray scattering of aerodynamically levitated and laser-beam-heated samples, with O/Pu ranging from 1.57 to 1.76. Liquid structural models consistent with the X-ray data are developed using machine-learned interatomic potentials and density functional theory. Molten PuO1.76 contains some degree of covalent Pu-O bonding, signalled by the degeneracy of Pu 5f and O 2p orbitals. The liquid is isomorphous with molten CeO1.75, demonstrating the latter as a non-radioactive, non-toxic, structural surrogate when differences in the oxidation potentials of Pu and Ce are accounted for. These characterizations provide essential constraints for modelling pertinent to reactor safety design.

Modeling the Structure of Ketoprofen-Poly(vinylpyrrolidone) Amorphous Solid Dispersions with Empirical Potential Structure Refinements of X-ray Scattering Data.

The metastability of amorphous formulations poses barriers to their safe and widespread commercialization. The propensity of amorphous solid dispersions (ASDs) to crystallize is directly linked to their molecular structure. Amorphous structures are inherently complex and thus difficult to fully characterize by experiments, which makes structural simulations an attractive route for investigating which structural characteristics correlate with ASD stability. In this study, we use empirical potential structure refinement (EPSR) to create molecular models of ketoprofen-poly(vinylpyrrolidone) (KTP/PVP) ASDs with 0-75 wt % drug loading. The EPSR technique uses X-ray total scattering measurements as constraints, yielding models that are consistent with the X-ray data. We perform several simulations to assess the sensitivity of the EPSR approach to input parameters such as intramolecular bond rotations, PVP molecule length, and PVP tacticity. Even at low drug loading (25 wt %), ∼40% of KTP molecules participate in KTP-KTP hydrogen bonding. The extent of KTP-PVP hydrogen bonding does not decrease significantly at higher drug loadings. However, the models' relative uncertainties are too large to conclude whether ASDs' lower stabilities at high drug loadings are due to changes in drug-excipient hydrogen bonding or a decrease in steric hindrance of KTP molecules. This study illustrates how EPSR, combined with total scattering measurements, can be a powerful tool for investigating structural characteristics in amorphous formulations and developing ASDs with improved stability.

Probing Molecular Packing of Amorphous Pharmaceutical Solids Using X-ray Atomic Pair Distribution Function and Solid-State NMR.

The structural investigation of amorphous pharmaceuticals is of paramount importance in comprehending their physicochemical stability. However, it has remained a relatively underexplored realm primarily due to the limited availability of high-resolution analytical tools. In this study, we utilized the combined power of X-ray pair distribution functions (PDFs) and solid-state nuclear magnetic resonance (ssNMR) techniques to probe the molecular packing of amorphous posaconazole and its amorphous solid dispersion at the molecular level. Leveraging synchrotron X-ray PDF data and employing the empirical potential structure refinement (EPSR) methodology, we unraveled the existence of a rigid conformation and discerned short-range intermolecular C-F contacts within amorphous posaconazole. Encouragingly, our ssNMR 19F-13C distance measurements offered corroborative evidence supporting these findings. Furthermore, employing principal component analysis on the X-ray PDF and ssNMR data sets enabled us to gain invaluable insights into the chemical nature of the intermolecular interactions governing the drug-polymer interplay. These outcomes not only furnish crucial structural insights facilitating the comprehension of the underlying mechanisms governing the physicochemical stability but also underscore the efficacy of synergistically harnessing X-ray PDF and ssNMR techniques, complemented by robust modeling strategies, to achieve a high-resolution exploration of amorphous structures.

Molecular interactions in short-chain perfluoroalkyl carboxylic acids and aqueous solutions.

The presence of short-chain per- and polyfluoroalkyl substances in water poses a major health and environmental challenge. Here, we have performed high-energy small- and wide-angle X-ray scattering measurements on CF3[CF2]nCOOH (where n = 1, 2, 3 represents the chain length) and their aqueous solutions at 10% mole concentrations to characterize their molecular interactions at the atomic and nanometer length scales. The experimental wide-angle structure factors have been modelled using Empirical Potential Structural Refinement. The oxygen-oxygen partial X-ray pair distribution functions show that the coordination number between the hydroxyl oxygen on the acid and surrounding oxygen water molecules increases significantly with acid chain length, rising from 3.2 for n = 1 to 4.1 for n = 3. The small-angle scattering is dominated by a sharp, high-intensity peak at Q1 ∼ 0.2 Å-1 and a smaller peak at Q2 = 1.2 Å-1 for n = 3, both of which decrease with decreasing chain length. The Q2 peak is attributed to groups of adjacent non-bonded acid molecules, and Q1 has contributions from both correlations between acid molecules and water-water interactions. In all cases, the models show nanoscale aggregation occurs in the form of denser channels of winding hydrogen-bonded chains, approximately 20 water molecules in length, surrounding clusters of acid molecules. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 2)'.

The structure of molten calcium ferrite under various redox conditions.

Laser-heated melts based on the 43CaO-57Fe2O3-x eutectic, close to the calcium ferrite (CF) composition, were measured with high-energy X-ray diffraction using aerodynamic levitation over a range of redox states controlled by CO/CO2 gas atmospheres. The iron-oxygen coordination number was found to rise from 4.4 ± 0.3 at 15% Fe3+ to 5.3 ± 0.3 at 87% Fe3+. Empirical potential structure refinement modelling was used to obtain the ferric and ferrous partial pair distribution functions. It was found that the Fe2+ iron-oxygen coordination number is consistently approximately 10% higher in CF than in pure iron oxide, while Fe3+ is essentially identical in all but the most oxygen-rich environments (where it is higher in CF compared with FeOx). The model also shows calcium octahedra to be the dominant species across all redox environments, although the population of CaO7 increases with the availability of oxygen at the expense of CaO4 and CaO5. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

The local ordering of polar solvents around crystalline carbon nitride nanosheets in solution.

The crystalline graphitic carbon nitride, poly-triazine imide (PTI) is highly unusual among layered materials since it is spontaneously soluble in aprotic, polar solvents including dimethylformamide (DMF). The PTI material consists of layers of carbon nitride intercalated with LiBr. When dissolved, the resulting solutions consist of dissolved, luminescent single to multilayer nanosheets of around 60-125 nm in diameter and Li+ and Br- ions originating from the intercalating salt. To understand this unique solubility, the structure of these solutions has been investigated by high-energy X-ray and neutron diffraction. Although the diffraction patterns are dominated by inter-solvent correlations there are clear differences between the X-ray diffraction data of the PTI solution and the solvent in the 4-6 Å-1 range, with real space differences persisting to at least 10 Å. Structural modelling using both neutron and X-ray datasets as a constraint reveal the formation of distinct, dense solvation shells surrounding the nanoparticles with a layer of Br-close to the PTI-solvent interface. This solvent ordering provides a configuration that is energetically favourable underpinning thermodynamically driven PTI dissolution. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 2)'.

AL4GAP: Active learning workflow for generating DFT-SCAN accurate machine-learning potentials for combinatorial molten salt mixtures.

Machine learning interatomic potentials have emerged as a powerful tool for bypassing the spatiotemporal limitations of ab initio simulations, but major challenges remain in their efficient parameterization. We present AL4GAP, an ensemble active learning software workflow for generating multicomposition Gaussian approximation potentials (GAP) for arbitrary molten salt mixtures. The workflow capabilities include: (1) setting up user-defined combinatorial chemical spaces of charge neutral mixtures of arbitrary molten mixtures spanning 11 cations (Li, Na, K, Rb, Cs, Mg, Ca, Sr, Ba and two heavy species, Nd, and Th) and 4 anions (F, Cl, Br, and I), (2) configurational sampling using low-cost empirical parameterizations, (3) active learning for down-selecting configurational samples for single point density functional theory calculations at the level of Strongly Constrained and Appropriately Normed (SCAN) exchange-correlation functional, and (4) Bayesian optimization for hyperparameter tuning of two-body and many-body GAP models. We apply the AL4GAP workflow to showcase high throughput generation of five independent GAP models for multicomposition binary-mixture melts, each of increasing complexity with respect to charge valency and electronic structure, namely: LiCl-KCl, NaCl-CaCl2, KCl-NdCl3, CaCl2-NdCl3, and KCl-ThCl4. Our results indicate that GAP models can accurately predict structure for diverse molten salt mixture with density functional theory (DFT)-SCAN accuracy, capturing the intermediate range ordering characteristic of the multivalent cationic melts.

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.

Structures of glasses created by multiple kinetic arrests.

X-ray scattering has been used to characterize glassy itraconazole (ITZ) prepared by cooling at different rates. Faster cooling produces ITZ glasses with lower (or zero) smectic order with more sinusoidal density modulation, larger molecular spacing, and shorter lateral correlation between the rod-like molecules. We find that each glass is characterized by not one, but two fictive temperatures Tf (the temperature at which a chosen order parameter is frozen in the equilibrium liquid). The higher Tf is associated with the regularity of smectic layers and lateral packing, while the lower Tf with the molecular spacings between and within smectic layers. This indicates that different structural features are frozen on different timescales. The two timescales for ITZ correspond to its two relaxation modes observed by dielectric spectroscopy: the slower δ mode (end-over-end rotation) is associated with the freezing of the regularity of molecular packing and the faster α mode (rotation about the long axis) with the freezing of the spacing between molecules. Our finding suggests a way to selectively control the structural features of glasses.

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.

Experimentally Driven Automated Machine-Learned Interatomic Potential for a Refractory Oxide.

Understanding the structure and properties of refractory oxides is critical for high temperature applications. In this work, a combined experimental and simulation approach uses an automated closed loop via an active learner, which is initialized by x-ray and neutron diffraction measurements, and sequentially improves a machine-learning model until the experimentally predetermined phase space is covered. A multiphase potential is generated for a canonical example of the archetypal refractory oxide, HfO_{2}, by drawing a minimum number of training configurations from room temperature to the liquid state at ∼2900 °C. The method significantly reduces model development time and human effort.

Structure of ice confined in silica nanopores.

Observed anomalous thermodynamic properties of confined water such as deviations in the melting point and freezing point motivate the determination of the structure of confined water as a function of pore size and temperature. In this study, we investigate the dynamic evolution of the structure of confined ice in SBA-15 porous materials with pore diameters of 4 nm, 6 nm, and 8 nm at temperatures ranging from 183 K to 300 K using in operando Wide-Angle X-Ray Scattering (WAXS) measurements, X-Ray Partial Distribution Function (PDF) measurements, and classical Molecular Dynamics (MD) simulations. Formation of hexagonal ice structures is noted in all the three pore sizes. In silica nanopores with diameters of 4 nm, cubic ice formation is noted in addition to hexagonal ice. Longer lasting hydrogen bonds and longer residence times of the water molecules in the first coordination shell contribute to observed crystalline organization of ice in confinement. Self-diffusion coefficients of confined liquid water, predicted from classical MD simulations, are four orders of magnitude higher compared to ice formed in confinement. These experimental and simulation results provide comprehensive insights underlying the organization of confined water and ice in silica nanopores and the underlying physico-chemical interactions that contribute to the observed structures.

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol.

Synchrotron x-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol, and D-sorbitol) from above the glass transition temperatures Tg to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently, reflecting the different natures of bonding. Upon cooling from 400 K, the O⋯O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrests at Tg; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above Tg, glycerol loses the fewest HBs per OH for a given temperature rise scaled by Tg, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C⋯C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above Tg causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol being the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

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.

In Situ High-Temperature Synchrotron Diffraction Studies of (Fe,Cr,Al)3O4 Spinels.

The modeling of a loss-of-coolant-accident scenario involving nuclear fuels with FeCrAl cladding materials in consideration to replace a Zircaloy requires knowledge of the thermodynamics of oxidized structures. At temperatures higher than 1500 °C, oxidation of FeCrAl alloys forms (Fe,Cr,Al)3O4 spinels. In situ high-energy X-ray diffraction in a conical nozzle levitator installed at beamline 6-ID-D of the APS was used to study the structural evolution of the oxides as a function of the temperature. Single-phase (spinel) and multiphase (spinel-corundum-FeAlO3) regions are mapped as a function of the temperature for three different compositions of FeCrAl oxidation products. The thermal expansion coefficients and cation distribution in the spinel structure have been refined. The temperature at which complete melting of the fuel cladding is expected has been determined by the liquidus temperatures of the oxidized products to be between 1657 and 1834 °C in a 20% O2/Ar atmosphere using the cooling trace method. The liquidus temperature increases with increasing Al and Cr content in the spinel phase.

Diffusive dynamics during the high-to-low density transition in amorphous ice.

Water exists in high- and low-density amorphous ice forms (HDA and LDA), which could correspond to the glassy states of high- (HDL) and low-density liquid (LDL) in the metastable part of the phase diagram. However, the nature of both the glass transition and the high-to-low-density transition are debated and new experimental evidence is needed. Here we combine wide-angle X-ray scattering (WAXS) with X-ray photon-correlation spectroscopy (XPCS) in the small-angle X-ray scattering (SAXS) geometry to probe both the structural and dynamical properties during the high-to-low-density transition in amorphous ice at 1 bar. By analyzing the structure factor and the radial distribution function, the coexistence of two structurally distinct domains is observed at T = 125 K. XPCS probes the dynamics in momentum space, which in the SAXS geometry reflects structural relaxation on the nanometer length scale. The dynamics of HDA are characterized by a slow component with a large time constant, arising from viscoelastic relaxation and stress release from nanometer-sized heterogeneities. Above 110 K a faster, strongly temperature-dependent component appears, with momentum transfer dependence pointing toward nanoscale diffusion. This dynamical component slows down after transition into the low-density form at 130 K, but remains diffusive. The diffusive character of both the high- and low-density forms is discussed among different interpretations and the results are most consistent with the hypothesis of a liquid-liquid transition in the ultraviscous regime.

X-ray studies of the transformation from high- to low-density amorphous water.

Here we report about the structural evolution during the conversion from high-density amorphous ices at ambient pressure to the low-density state. Using high-energy X-ray diffraction, we have monitored the transformation by following in reciprocal space the structure factor SOO( Q) and derived in real space the pair distribution function gOO( r). Heating equilibrated high-density amorphous ice (eHDA) at a fast rate (4 K min-1), the transition to the low-density form occurs very rapidly, while domains of both high- and low-density coexist. On the other hand, the transition in the case of unannealed HDA (uHDA) and very-high-density amorphous ice is more complex and of continuous nature. The direct comparison of eHDA and uHDA indicates that the molecular structure of uHDA contains a larger amount of tetrahedral motives. The different crystallization behaviour of the derived low-density amorphous states is interpreted as emanating from increased tetrahedral coordination present in uHDA. This article is part of the theme issue 'The physics and chemistry of ice: scaffolding across scales, from the viability of life to the formation of planets'.

Experimental and Theoretical Insights into the Structure of Tellurium Chloride Glasses.

The structure of the binary chalcohalide glasses Te1- xCl x (0.35 ≤ x ≤ 0.65) is considered by combining experimental and theoretical results. The structural network properties are influenced by a competition between ionic and covalent bonding in such glasses. At first, a focus is placed on the detailed information available by using the complementary high-energy X-ray and the neutron diffractions in both the reciprocal and real spaces. The main characteristic suggested by the structure factors S( Q) concerns the presence of three length scales in the intermediate range order. The total correlation function T( r) lets us also suppose that the structure of these glasses is more complicated than Te-chain fragments with terminal Cl as demonstrated in crystalline Te3Cl2. Molecular dynamics simulations were subsequently performed on Te3Cl2 and Te2Cl3, and coupled with the experimental data, a highly reticulated network of chalcogen atoms, with a fair amount of chlorine atoms bonded in a bridging mode, is proposed. The simulations clearly lead to a glass description that differs markedly from the simple structural model based on only Te atom chains and terminal Cl atoms. Solid-state NMR experiments and NMR parameters calculations allowed validation of the presence of Te highly coordinated with chlorine in these glasses.

Structure and Liquid Fragility in Sodium Carbonate.

The relationship between local structure and dynamics is explored for molten sodium carbonate. A flexible fluctuating-charge model, which allows for changes in the shape and charge distribution of the carbonate molecular anion, is developed. The system shows the evolution of highly temperature-dependent complex low-dimensional structures which control the dynamics (and hence the liquid fragility). By varying the molecular anion charge distribution, the key interactions responsible for the formation of these structures can be identified and rationalized. An increase in the mean charge separation within the carbonate ions increases the connectivity of the emerging structures and leads to an increase in the system fragility.

The structure of liquid alkali nitrates and nitrites.

High energy X-ray diffraction has been combined with containerless techniques to determine the structure of a series of alkali and ammonium nitrate and nitrite liquids. The systems have been modelled using molecular dynamics simulation which allows for the flexibility of, and movement of charge within, the molecular anions. The model reproduces the experimentally-determined scattering functions in both the low- and high-Q regimes reflecting the inter- and intra-molecular length-scales. For ammonium nitrate the best fit to the diffraction data is obtained by assuming the NH4+ cation to have a radius closer to that for Cs+ rather than a smaller cation such as Rb+ as often previously assumed. The alkali nitrites show an emergent length scale, attributed to the nitrogen-nitrogen spatial correlations, that depends on both temperature and the identity of the alkali cation. The corresponding nitrates show a more subtle effect in the nitrogen-nitrogen correlations. As a result, the nature of this N-N length-scale appears different for the respective nitrites and nitrates.