Pair Distribution Function Research Articles

Correlated Anion Disorder in Heteroanionic Cubic TiOF2.

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO3-type) TiOF2 using X-ray pair distribution function (PDF), 19F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF2, characterized by predominant cis-[O2F4] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF4. We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

Comprehensive Investigation of the Crystal Structure of Cation-Disordered Li3VO4 as a High-Rate Anode Material: Unveiling the Dichotomy between Order and Disorder.

This study investigates mechanochemical synthesis and cation-disordering mechanism of wurtzite-type Li3VO4 (LVO), highlighting its promise as a high-performance anode material for lithium-ion batteries and hybrid supercapacitors. Mechanochemical treatment of pristine LVO using a high-energy ball mill results in a "pure cation-disordered" LVO phase, allowing for meticulous analysis of cation arrangement. The X-ray and neutron diffraction study demonstrates progressive loss of order in LVO crystal with increasing milling duration. High-resolution transmission electron microscopy reveals disrupted lattice fringes, indicating cationic misalignment. Pair-distribution function analysis confirms loss of cation arrangements and the presence of short-range order. Combination of these multiple analytical techniques achieves a comprehensive understanding of cation regularity and clearly demonstrates order/disorder dichotomy in cation-disordered materials, ranging from short (<8 Å) to middle-long range (8-30 Å), using an integrated superstructure model of the cation-disordered LVO crystals. Electrochemical testing reveals that mechanochemically treated LVO exhibits superior rate capability, with a 70% capacity retention at a high current density of 50C-rate. Lithium diffusion coefficient measurements demonstrate enhanced lithium-ion mobility in the mechanochemically treated LVO, attributed to cation-disordering effect. These findings provide valuable insights into mechanochemical cation-disordering in LVO, presenting its potential as an efficient anode material for lithium-ion-based electrochemical energy storage.

Structure evolution of the silica glass under the vibration field through molecular dynamics simulations

AbstractSilica glass has found wide applications, including space shuttle windows, optical fibers, and deep ultraviolet light lens elements, due to its excellent optical, electrical, and mechanical properties. However, understandings of its structure change under the vibration field are limited in literature. In this work, the structure of silica glass and its behavior under the vibration field have been investigated through molecular dynamics simulations. The short and medium‐range structures, including coordination number, Qn species distribution, pair distribution function, and bond angle distribution, were analyzed. It was found that the overall and local structures of silica glass showed periodic changes under the vibration field. The root‐mean‐square displacements as well as the kinetic energies of the atoms in the system were further analyzed. The results show that there is a periodic transformation of kinetic and potential energies. In addition, the kinetic energy decreases and is converted into the thermal energy during the vibration, which is confirmed by the increase of the system temperature. This can explain the thermoelastic loss mechanism in energy dissipation from literature. This work reveals the structure evolution and energy change of glass under the vibration field from atomistic level, which provides important understandings for related researches.

Instrumental broadening and the radial pair distribution function with 2D detectors.

The atomic pair distribution function (PDF) is a real-space representation of the structure of a material. Experimental PDFs are obtained using a Fourier transform from total scattering data which may or may not have Bragg diffraction peaks. The determination of Bragg peak resolution in scattering data from the fundamental physical parameters of the diffractometer used is well established, but after the Fourier transform from reciprocal to direct space, these contributions are harder to identify. Starting from an existing definition of the resolution function of large-area detectors for X-ray diffraction, this approach is expanded into direct space. The effect of instrumental parameters on PDF peak resolution is developed mathematically, then studied with modelling and comparison with experimental PDFs of LaB6 from measurements made in different-sized capillaries.

Large dynamic scissoring mode displacements coupled to band gap opening in the cubic phase of the methylammonium lead halide perovskites

Hybrid perovskites are a rapidly growing research area, having reached photovoltaic power conversion efficiencies of over 25%. There is a increasing consensus that the structures of these materials, and hence their electronic structures, cannot be understood purely from the time and space averaged crystal structures observable by conventional methods. We apply a symmetry-motivated analysis method to analyse x-ray pair distribution function data of the cubic phases of the hybrid perovskites MAPbX 3 (X = I, Br, Cl). We demonstrate that, even in the cubic phase, the local structure of the inorganic components of MAPbX 3 (X = I, Br, Cl), are dominated by scissoring type deformations of the PbX 6 octahedra. We find these modes to have a larger amplitude than equivalent distortions in the A-site deficient perovskite ScF3 and demonstrate that they show a significant departure from the harmonic approximation. Calculations performed on an inorganic perovskite analogue, FrPbBr3, show that the large amplitudes of the scissoring modes are coupled to a dynamic opening of the electronic band gap. Finally, we use density functional theory calculations to show that the organic MA cations reorientate to accommodate the large amplitude scissoring modes.

Ultrafast photochemistry and electron diffraction for cyclobutanone in the S2 state: Surface hopping with time-dependent density functional theory.

We simulate the photodynamics of gas-phase cyclobutanone excited to the S2 state using fewest switches surface hopping (FSSH) dynamics powered by time-dependent density functional theory (TDDFT). We predict a total photoproduct yield of 8%, with a C3:C2 product ratio of 0 trajectories to 8 trajectories. One primary S2 → S1 conical intersection is identified involving the compression of an α-carbon-carbon-hydrogen bond angle. Excited state lifetimes computed with respect to electronic state populations were found to be 3.96ps (S2 → S1) and 498fs (S1 → S0). We also generate time-resolved difference pair distribution functions (ΔPDFs) from our TDDFT-FSSH dynamics results in order to generate direct comparisons with ultrafast electron diffraction experiment observables. Global and target analysis of time-resolved ΔPDFs produced a distinct set of lifetimes: (i) a 0.548ps decay and (ii) a 1.69ps decay, both resembling the S2 minimum, as well as (iii) a long decay that resembles the S1 minimum geometry and the fully separated C2 products. Finally, we contextualize our results by considering the impact of the most likely sources of significant errors.

Collective dynamics and pair-distribution function of active Brownian ellipsoids in two spatial dimensions

While the collective dynamics of spherical active Brownian particles is relatively well understood by now, the much more complex dynamics of nonspherical active particles still raises interesting open questions. Previous work has shown that the dynamics of rod-like or ellipsoidal active particles can differ significantly from that of spherical ones. Here, we obtain the full state diagram of active Brownian ellipsoids in two spatial dimensions without hydrodynamic interactions depending on the Péclet number and packing density via computer simulations. The system is found to exhibit a rich state behavior that includes cluster formation, local polar order, polar flocks, and disordered states. Moreover, we obtain numerical results and an analytical representation for the pair-distribution function of active ellipsoids. This function provides useful quantitative insights into the collective behavior of active particles with lower symmetry and has potential applications in the development of predictive theoretical models.

Different or the same? exploring the physicochemical properties and molecular mobility of celecoxib amorphous forms

The influence of different preparation methods on the physicochemical properties of amorphous solid forms have gained considerable attention, especially with recent publications on pharmaceutical polyamorphism. In the present study, we have investigated the possible occurrence of polyamorphism in the drug celecoxib (CEL) by investigating the thermal behavior, morphology, structure, molecular mobility and physical stability of amorphous CEL obtained by quench-cooling (QC), ball milling (BM) and spray drying (SD). Similar glass transition temperatures but different recrystallization behaviors were observed for CEL-QC, CEL-BM and CEL-SD using modulated differential scanning calorimetry analysis. A correlation between the different recrystallization behaviors of the three CEL amorphous forms and the respective distinct powder morphologies, was also found. Molecular dynamics simulations however, reveal that CEL presents similar molecular conformational distributions when subjected to QC and SD. Moreover, the obtained molecular conformational distributions of CEL are different from the ones found in its crystal structure and also from the ones found in the lowest-energy structure obtained by quantum mechanical calculations. The type and strength of CEL hydrogen bond interactions found in CEL-QC and CEL-SD systems are almost identical, though different from the ones presented in the crystal structure. Pair distribution function analyses and isothermal microcalorimetry show similar local structures and structural relaxation times, respectively, for CEL-QC, CEL-BM and CEL-SD. The present work shows that not only similar physicochemical properties (glass transition temperature, and structural relaxation time), but also similar molecular conformational distributions were observed for all prepared CEL amorphous systems. Hence, despite their different recrystallization behaviors, the three amorphous forms of CEL did not show any signs of polyamorphism.

π‐d Conjugated Copper Chloranilate with Distorted Cu‐O4 Site for Efficient Electrocatalytic Ammonia Production

AbstractUnderstanding the relationship between electrocatalytic performance and local structure at the molecular level is of great significance. Herein, a bifunctional electrocatalyst CuCA (CA = chloranilate) is constructed for both nitrogen reduction reaction (NRR) and nitrate reduction reaction (NO3RR). Combined structural analyses using Rietveld refinement, extended X‐ray adsorption fine structure (EXAFS), and pair distribution function (PDF) revealed a significant distortion of the Cu‐O4 structure. Benefitting from the unique local structure, Cu‐CA shows an impressive NH3 yield rate of 286.00 ug h−1 mg−1 (FE = 18.25%, ‐0.85 V vs RHE), 3180.00 ug h−1 mg−1 (FE = 90.3%, ‐0.9 V vs RHE) for NRR and NO3RR, respectively. In contrast, the pyrazine (Pyz) decorated compound Cu‐CA‐Pyz with a less distorted Cu‐O4 structure and fewer active sites show much lower activity. Density functional theory (DFT) calculations shed light on that the distorted nature can effectively regulate the electron density distribution, which can lower the energy barrier of adsorption and activation of the intermediate species, leading to the enhanced activity. These findings may give new insight into the structural‐property relationship and open up opportunities for the exploration of efficient electrocatalysts.

RZera: A customized and user-friendly software for the total scattering data reduction at CSNS

Reduction Zone Era (RZera), a customized and user-friendly software designed for the reduction of time-of-flight (TOF) data from the Multi-physics Instrument (MPI) at China Spallation Neutron Source (CSNS), is developed and presented in this paper. After inputting the required information, the software can transform and apply a series of corrections from raw data to user data that only reflects the intrinsic properties of the sample. The software has two modes: diffraction and pair distribution function reduction. Each mode is capable of processing batch data quickly. The nickel has been used to test the diffraction and the pair distribution function reduction mode separately. The diffraction result has been fitted by GSAS, and the pair distribution function result has been fitted by the DiffPy-CMI suite. All the fittings demonstrate the effectiveness of RZera. The future versions of RZera will integrate additional instrument data reduction capabilities, and incorporate advanced algorithms and features that address emerging challenges in neutron scattering research.

Influence of local structures on amorphous alumina exhibiting resistance random-access memory function

Amorphous alumina resistance random-access memory is a promising candidate as a next-generation nonvolatile memory. It is intriguing that the nonvolatile memory function emerges in only amorphous samples, unlike crystalline samples. We studied local structures of amorphous alumina samples and Al2O3 polycrystalline using atomic pair distribution function measurements. We derived the Al–Al, O–O, and Al–O atomic distances for each sample. By comparing them, we revealed that the subtle difference in the local structure significantly influences the performance of a nonvolatile memory function.

Improving Hydrogen Release From Oxygen‐Functionalized LOHC Molecules by Ru Addition to Pt/C Catalysts

AbstractBimetallic PtRu/C catalysts have been prepared by depositing Ru onto a commercial Pt/C catalyst and subsequent thermal annealing. Different Pt : Ru ratios of 1 : 1, 2 : 1, 4 : 1, 8 : 1, and 16 : 1 have been investigated with annealing temperatures of 350 and 500 °C. Using X‐ray diffraction (XRD), synchrotron powder X‐ray diffraction (SPXRD) combined with pair distribution function (PDF) analysis and scanning transmission electron microscopy with energy dispersive X‐ray spectrum imaging (STEM‐EDXS), we show that diffusion of Ru into the fcc crystal structure of the Pt nanoparticles takes place during thermal annealing. Partial segregation and the formation of Pt nanoclusters on the bimetallic surface was observed depending on the Pt : Ru ratio. The annealing process does not only cause a high Pt dispersion but also affects the electronic structure of the Pt surface by broadening the d‐band thus facilitating the product's desorption from the surface of the bimetallic particle. These beneficial effects have been exemplified for the catalytic dehydrogenation of dicyclohexylmethanol, a promising hydrogen carrier compound with 7.2 wt% hydrogen capacity. The most active catalyst material among the tested supported alloys and monometallic reference materials was Pt4Ru1/C after annealing at 500 °C, which increased the hydrogen release productivity by 65 % compared to the unmodified Pt/C catalyst.

Multiscale investigation on the formation path of the apatite phase in bioactive glasses

Mesoporous bioactive glasses (MBGs) have been reported as promising biomaterials to stimulate and promote the growth of bones. When exposed to Simulated Body Fluid (SBF), MBGs develop an apatite phase called carbonate hydroxyapatite (H(C)A) on their surfaces that closely mimics the mineral phase of bones. In order to gain deeper insights into the key stages of the surface reactions of two different types of MBGs (58S and 70S30C) soaked in SBF, we have used high energy total X-ray scattering coupled to Pair Distribution Function (PDF) analysis, along with Raman spectroscopy and scanning electron microscopy. This type of coupled analytical approach can ultimately help to unravel the details of the reaction kinetics in the complex calcium phosphate environment around the BGs surfaces. The data obtained in our in-vitro study point to the formation of other intermediate phases like amorphous calcium phosphate (ACP) and calcite (CC), that transform after a specific time into H(C)A depending on the composition of the MBGs, their surface properties and their interaction times with SBF.

Approaches and challenges in whole-nanoparticle refinements from neutron total-scattering data

This study considers critical data reduction steps and data analysis approaches required to determine explicitly the atomic arrangements in nanoparticles from time-of-flight neutron total scattering. A practical procedure is described for removing parasitic backgrounds caused by the incoherent scattering of hydrogen inevitably present in most nanoparticle samples and normalizing the recovered coherent scattering intensities onto an absolute scale. A model-free analysis is presented of a pair-distribution function derived from total scattering that can be used to determine thermal expansion coefficients and particle sizes directly from experimental data without knowledge of the material's structure. Finally, atomistic whole-nanoparticle refinements of yttrium-doped ceria nanoparticles from neutron total-scattering data are demonstrated using the reverse Monte Carlo method implemented in the RMCProfile software. These results reveal a strong dependence of the cation–oxygen and oxygen–oxygen distances on the coordination numbers, which leads to gradients of these distances near the particle surface. The details are dependent on the surface coverage by ligands and adsorbates and on the structure of grain boundaries in nanocrystalline agglomerates. The refined models confirm the expectations of more substantial disorder at particle surfaces, with a distorted surface layer extending over several coordination shells. The results highlight the feasibility of whole-nanoparticle refinements from neutron data, calling for further development of data reduction and analysis procedures.

GIWAXS experimental methods at the NFPS-BL17B beamline at Shanghai Synchrotron Radiation Facility.

The BL17B beamline at the Shanghai Synchrotron Radiation Facility was first designed as a versatile high-throughput protein crystallography beamline and one of five beamlines affiliated to the National Facility for Protein Science in Shanghai. It was officially opened to users in July 2015. As a bending magnet beamline, BL17B has the advantages of high photon flux, brightness, energy resolution and continuous adjustable energy between 5 and 23 keV. The experimental station excels in crystal screening and structure determination, providing cost-effective routine experimental services to numerous users. Given the interdisciplinary and green energy research demands, BL17B beamline has undergone optimization, expanded its range of experimental methods and enhanced sample environments for a more user-friendly testing mode. Thesemethods include single-crystal X-ray diffraction, powder crystal X-ray diffraction, wide-angle X-ray scattering, grazing-incidence wide-angle X-ray scattering (GIWAXS), and fully scattered atom pair distribution function analysis, covering structure detection from crystalline to amorphous states. This paper primarily presents the performance of the BL17B beamline and the application of the GIWAXS methodology at the beamline in the field of perovskite materials.

Deciphering diffuse scattering with machine learning and the equivariant foundation model: the case of molten FeO

Bridging the gap between diffuse x-ray or neutron scattering measurements and predicted structures derived from atom–atom pair potentials in disordered materials, has been a longstanding challenge in condensed matter physics. This perspective gives a brief overview of the traditional approaches employed over the past several decades. Namely, the use of approximate interatomic pair potentials that relate three-dimensional structural models to the measured structure factor and its’ associated pair distribution function. The use of machine learned interatomic potentials has grown in the past few years, and has been particularly successful in the cases of ionic and oxide systems. Recent advances in large scale sampling, along with a direct integration of scattering measurements into the model development, has provided improved agreement between experiments and large-scale models calculated with quantum mechanical accuracy. However, details of local polyhedral bonding and connectivity in meta-stable disordered systems still require improvement. Here we leverage MACE-MP-0; a newly introduced equivariant foundation model and validate the results against high-quality experimental scattering data for the case of molten iron(II) oxide (FeO). These preliminary results suggest that the emerging foundation model has the potential to surpass the traditional limitations of classical interatomic potentials.

PyCSH: Automated atomic-level structure generation of bulk C-S-H and investigation of their intrinsic properties

A new Python code for the automated generation of realistic bulk calcium silicate hydrate (C-S-H) structures is introduced. The code was used to generate 400 structures with Ca/Si of 1.3, 1.5, 1.7, and 1.9. The generated structures are in excellent agreement with experimentally measured C-S-H properties (Ca/Si, 2H/Si, MCL, Si-OH/Si, and Ca-OH/Ca). Molecular dynamics was used to simulate the structures, which were then investigated for their structural features and energetic stability. The results indicate very similar short-range ordering and energetic stability between all 400 structures. Finally, it is shown how computational C-S-H models can be used to understand the experimentally measured pair distribution functions. The code, named pyCSH, is available as open source under GPL-2 license at GitHub https://github.com/hegoimanzano/pyCSH.git and https://github.com/akmiitm/pyCSH.git.

Structural and elastic properties of zirconium silicate glasses: Insight from molecular dynamic simulation

Using Morse potentials, we performed molecular dynamics simulations to study the effect of ZrO2/SiO2 substitution on the structure and mechanical properties of SiO2-ZrO2 glasses with various ZrO2 concentrations (5, 10, 20, 30, 40, and 50 mol%). We computed the pair distribution function (PDF) for all atomic pairs and extracted the cation-oxygen coordination number and bond distance as a function of composition for Zr4+ and other cations. We found that Zr entered the glass network as distorted [ZrO6] octahedra. Our results showed that the ZrO2 content was positively associated with the number of non-bridging oxygens and the emergence of tricluster oxygens, which consequently reduced the glass transition temperature. The elastic constants such as Young’s, bulk and shear modulus, increase with rising zirconium content.

High-throughput and high-resolution powder X-ray diffractometer consisting of six sets of 2D CdTe detectors with variable sample-to-detector distance and innovative automation system.

The demand for powder X-ray diffraction analysis continues to increase in a variety of scientific fields, as the excellent beam quality of high-brightness synchrotron light sources enables the acquisition of high-quality measurement data with high intensity and angular resolution. Synchrotron powder diffraction has enabled the rapid measurement of many samples and various in situ/operando experiments in nonambient sample environments. To meet the demands for even higher throughput measurements using high-energy X-rays at SPring-8, a high-throughput and high-resolution powder diffraction system has been developed. This system is combined with six sets of two-dimensional (2D) CdTe detectors for high-energy X-rays, and various automation systems, including a system for automatic switching among large sample environmental equipment, have been developed in the third experimental hutch of the insertion device beamline BL13XU at SPring-8. In this diffractometer system, high-brilliance and high-energy X-rays ranging from 16 to 72 keV are available. The powder diffraction data measured under ambient and various nonambient conditions can be analysed using Rietveld refinement and the pair distribution function. Using the 2D CdTe detectors with variable sample-to-detector distance, three types of scan modes have been established: standard, single-step and high-resolution. A major feature is the ability to measure a whole powder pattern with millisecond resolution. Equally important, this system can measure powder diffraction data with high Q exceeding 30 Å-1 within several tens of seconds. This capability is expected to contribute significantly to new research avenues using machine learning and artificial intelligence by utilizing the large amount of data obtained from high-throughput measurements.

Tracking copper nanofiller evolution in polysiloxane during processing into SiOC ceramic.

Polymer-derived ceramics (PDCs) remain at the forefront of research for a variety of applications including ultra-high-temperature ceramics, energy storage and functional coatings. Despite their wide use, questions remain about the complex structural transition from polymer to ceramic and how local structure influences the final microstructure and resulting properties. This is further complicated when nanofillers are introduced to tailor structural and functional properties, as nanoparticle surfaces can interact with the matrix and influence the resulting structure. The inclusion of crystalline nanofiller produces a mixed crystalline-amorphous composite, which poses characterization challenges. With this study, we aim to address these challenges with a local-scale structural study that probes changes in a polysiloxane matrix with incorporated copper nanofiller. Composites were processed at three unique temperatures to capture mixing, pyrolysis and initial crystallization stages for the pre-ceramic polymer. We observed the evolution of the nanofiller with electron microscopy and applied synchrotron X-ray diffraction with differential pair distribution function (d-PDF) analysis to monitor changes in the matrix's local structure and interactions with the nanofiller. The application of the d-PDF to PDC materials is novel and informs future studies to understand interfacial interactions between nanofiller and matrix throughout PDC processing.