RMCProfile7: reverse Monte Carlo for multiphase systems.

This work introduces a completely rewritten version of the program RMCProfile (version 7), big-box, reverse Monte Carlo modelling software for analysis of total scattering data. The major new feature of RMCProfile7 is the ability to refine multiple phases simultaneously, which is relevant for many current research areas such as energy materials, catalysis and engineering. Other new features include improved support for molecular potentials and rigid-body refinements, as well as multiple different data sets. An empirical resolution correction and calculation of the pair distribution function as a back-Fourier transform are now also available. RMCProfile7 is freely available for download at https://rmcprofile.ornl.gov/.

Structure of Polaronic Centers in Proton-Intercalated AWO4 Scheelite-Type Tungstates.

The studies of polaronic centers in a homologous series of scheelite-type compounds AWO4 (A = Ca, Sr, Ba) were performed using the W L3-edge and Sr K-edge X-ray absorption spectroscopy combined with the reverse Monte Carlo simulations, X-ray photoelectron spectroscopy (XPS), and first-principles calculations. Protonated scheelites HxAWO4 were produced using acid electrolytes in a one-step route at ambient conditions. The underlying mechanism behind this phenomenon can be ascribed to the intercalation of H+ into the crystal structure of tungstate, effectively resulting in the reduction of W6+ to W5+, i.e., the formation of polaronic centers, and giving rise to a characteristic dark blue-purple color. The emergence of the W5+ was confirmed by XPS experiments. The relaxation of the local atomic structure around the W5+ polaronic center was determined from the analysis of the extended X-ray absorption fine structures using the reverse Monte Carlo method. The results obtained suggest the displacement of the W5+ ions from the center of [W5+O4] tetrahedra in the structure of AWO4 scheelite-type tungstates. This finding was also supported by the results of the first-principles calculations.

The ordered lattice host framework induced guest Li+ disordering in high performance cobalt-free Ni-rich cathode materials.

Recently, due to high price, resource shortage and unstable supply of cobalt, the development of low-cost cobalt-free Ni-rich cathodes has attracted extensive attention with the ever-increasing lithium-ion batteries (LIBs) industry. Selecting cost-effective elements to replace cobalt in Ni-rich cathodes is urgent. However, the principle of structural design of Ni-rich cathode remains unclear, hampering the selection of alternative elements. Herein, the cobalt-free cathodes of LiNi0.95Mg0.05O2 (NiMg) and LiNi0.95Mn0.05O2 (NiMn) are designed as alternatives to LiNi0.96Co0.04O2 (NiCo). NiMg has comparable cycle stability with NiCo, while NiMn has inferior cycle performance. Reverse Monte Carlo modelling was used to generate structural model and uncover local structure by fitting pair distribution function. It reveals Mn causes more severe Jahn-Teller distortions and disordered lattice host framework (Ni0.95M0.05O2, M = Co/Mn/Mg) than Co and Mg due to the strong size effect and coulomb interactions of Mn in Ni0.95Mn0.05O2 layer. The outstanding cycle stability of NiMg and NiCo originates from the ordered lattice host frameworks, which relieve stress and inhibit particle breakage during cycle. Meanwhile, the ordered lattice host framework induced guest Li+ disordering reduces Li+ diffusion energy barrier, improving the rate capability. This study provides a new perspective for the structural design of cobalt-free Ni-rich cathodes.

Evolution of Short-Range Order of Amorphous GeTe Upon Structural Relaxation Obtained by TEM Diffractometry and RMC Methods.

Glasses frequently reveal structural relaxation that leads to changes in their physical properties including enthalpy, specific volume, and resistivity. Analyzing the short-range order (SRO) obtained from electron diffraction by transmission electron microscopy (TEM) in combination with Reverse-Monte-Carlo (RMC) simulations is shown to provide information on the atomic arrangement. The technique elaborated here features several benefits including reliability, accessibility, and allows for obtaining detailed structural data quickly. This is demonstrated with a detailed view of the structural changes in the as-deposited amorphous phase change material (PCM) GeTe. The data show a significant increase in the average bond angle upon thermal treatment. At the same time the fraction of tetrahedrally coordinated Ge atoms decreases due to an increase in octahedrally distorted and pyramidal motifs. This finding provides further evidence for the atomic processes that govern structural relaxation in amorphous GeTe and other PCMs. A thorough literature review finally unveils possible origins of the large discrepancies reported on the structure of amorphous GeTe.

Examining the Effect of Cu and Mn Dopants on the Structure of Zinc Blende ZnS Nanopowders.

It is known that doping zinc sulfide (ZnS) nanoparticles with Mn or Cu ions significantly affects their luminescent properties. Herein, we investigated how dopant atoms are incorporated into the structure of ZnS using X-ray diffraction and multi-edge X-ray absorption spectroscopy. The observed broadening of the X-ray diffraction patterns indicates an average crystallite size of about 6 nm. By analyzing the Zn, Mn, and Cu K-edge extended X-ray absorption fine structure (EXAFS) spectra using the reverse Monte Carlo method, we were able to determine the relaxations of the local environments around the dopants. Our findings suggested that upon the substitution of Zn by Mn or Cu ions, there is a shortening of the Cu-S bonds by 0.08 Å, whereas the Mn-S bonds exhibited lengthening by 0.07 Å. These experimental results were further confirmed by first-principles density functional theory calculations, which explained the increase in the Mn-S bond lengths due to the high-spin state of Mn2+ ions.

Comparative study of the short-range structure of α-V2O5, α-TeO2 and xV2O5-(100 - x)TeO2 glasses using X-ray diffraction, Rietveld analysis and reverse Monte Carlo simulations.

Vanadium-tellurite glasses, tetragonal TeO2 and orthorhombic V2O5 crystalline samples were characterized for their atomic structure properties by synchrotron X-ray diffraction, pair distribution function analysis, reverse Monte Carlo simulations (RMC) and Rietveld analysis. The pair correlation function, G(r), of V2O5 shows the first peak at 1.61 Å. G(r) of TeO2 shows three peaks at 1.57, 2.13 and 2.88 Šdue to Te-O linkages of three different lengths, whereas the Te-Te atomic pair correlation shows a peak at 3.85 Å. The average coordination number of V with O in crystalline V2O5 is 4.39 while that of Te with O in crystalline TeO2 is 3.71. G(r) of the vanadium tellurite glass shows the first peak at 1.90 Šdue to overlapping Te-O and V-O atomic pair correlations. The RMC analysis on diffraction data of glasses found that the V-O coordination number is in the range 5.27-5.59 and the Te-O coordination number is 5.39-5.67. However, it is found that these coordination numbers cannot be clearly defined due to short-range disorder.

Influence of the Graft Length on Nanocomposite Structure and Interfacial Dynamics.

Both the dispersion state of nanoparticles (NPs) within polymer nanocomposites (PNCs) and the dynamical state of the polymer altered by the presence of the NP/polymer interfaces have a strong impact on the macroscopic properties of PNCs. In particular, mechanical properties are strongly affected by percolation of hard phases, which may be NP networks, dynamically modified polymer regions, or combinations of both. In this article, the impact on dispersion and dynamics of surface modification of the NPs by short monomethoxysilanes with eight carbons in the alkyl part (C8) is studied. As a function of grafting density and particle content, polymer dynamics is followed by broadband dielectric spectroscopy and analyzed by an interfacial layer model, whereas the particle dispersion is investigated by small-angle X-ray scattering and analyzed by reverse Monte Carlo simulations. NP dispersions are found to be destabilized only at the highest grafting. The interfacial layer formalism allows the clear identification of the volume fraction of interfacial polymer, with its characteristic time. The strongest dynamical slow-down in the polymer is found for unmodified NPs, while grafting weakens this effect progressively. The combination of all three techniques enables a unique measurement of the true thickness of the interfacial layer, which is ca. 5 nm. Finally, the comparison between longer (C18) and shorter (C8) grafts provides unprecedented insight into the efficacy and tunability of surface modification. It is shown that C8-grafting allows for a more progressive tuning, which goes beyond a pure mass effect.

Combining reverse Monte Carlo analysis of X-ray scattering and extended X-ray absorption fine structure spectra of very small nanoparticles.

Finite size effects in partial pair distribution functions generate artefacts in the scattering structure factor and scattering intensity. It is shown how they can be overcome using a binned version of the Debye scattering equation. Accordingly, reverse Monte Carlo simulations are used for very small nanoparticles of LaFeO3 with diameters below 10 nm to simultaneously analyse X-ray scattering data and extended X-ray absorption fine structure spectra at the La K and Fe K edges. The structural information obtained is consistent regarding local structure and long-range order.

A reverse Monte Carlo algorithm to simulate two-dimensional small-angle scattering intensities.

Small-angle scattering (SAS) experiments are a powerful method for studying self-assembly phenomena in nanoscopic materials because of the sensitivity of the technique to structures formed by interactions on the nanoscale. Numerous out-of-the-box options exist for analysing structures measured by SAS but many of these are underpinned by assumptions about the underlying interactions that are not always relevant for a given system. Here, a numerical algorithm based on reverse Monte Carlo simulations is described to model the intensity observed on a SAS detector as a function of the scattering vector. The model simulates a two-dimensional detector image, accounting for magnetic scattering, instrument resolution, particle polydispersity and particle collisions, while making no further assumptions about the underlying particle interactions. By simulating a two-dimensional image that can be potentially anisotropic, the algorithm is particularly useful for studying systems driven by anisotropic interactions. The final output of the algorithm is a relative particle distribution, allowing visualization of particle structures that form over long-range length scales (i.e. several hundred nanometres), along with an orientational distribution of magnetic moments. The effectiveness of the algorithm is shown by modelling a SAS experimental data set studying finite-length chains consisting of magnetic nanoparticles, which assembled in the presence of a strong magnetic field due to dipole interactions.

Pressure cells for in situ neutron total scattering: time and real-space resolution during deuterium absorption.

In situ gas-loading sample holders for two-dimensionally arranged detectors in time-of-flight neutron total scattering experiments have been developed to investigate atomic arrangements during deuterium absorption using time and real-space resolution. A single-crystal sapphire container was developed that allows conditions of 473 K and 10 MPa hydrogen gas pressure. High-resolution transient measurements detected deuterium absorption by palladium that proceeded within a few seconds. A double-layered container with thick- and thin-walled vanadium allowed conditions of 423 K and 10 MPa hydrogen gas pressure. The deuterium occupation sites of a lanthanum-nickel-aluminium alloy are discussed in detail on the basis of real-space high-resolution data obtained from in situ neutron scattering measurements and reverse Monte Carlo structural modeling.

Ferroelectric Ordering in Nanosized PbTiO3.

The insight into the three-dimensional configuration of ferroelectric ordering in ferroelectric nanomaterials is motivated by the application of the development of functional nanodevices and the structural designing. However, the atomic deciphering of the spatial distribution of ordered structure remains challenging for the limitation of dimension and probing techniques. In this paper, a neutron pair distribution function (nPDF) was utilized to analyze the spontaneous polarization distribution of zero-dimensional PbTiO3 nanoparticles in three dimensions, via the application of reverse Monte Carlo (RMC) modeling. The comprehensive identification with transmission electron microscopy verified the linear characteristics of polarization along the c-axis in the main body, while electric polarization distribution on the surface was enhanced abnormally. In addition, the correlation of dipole vectors extending to three unit cells below the surface is retained. This work shows an application of the micro/macroscale information to effectively decode the polarization structure of nanoferroelectrics, providing new views of designing nanoferroelectric devices.

Octahedral Tilting in Homologous Perovskite Series CaMoO3-SrMoO3-BaMoO3 Probed by Temperature-Dependent EXAFS Spectroscopy.

Local distortions in perovskites can be induced by cation displacements and/or by the tilting and rotating of cation-anion octahedra. Both phenomena have been subject to intense investigations over many years. However, there are still controversies in the results obtained from experimental techniques that are sensitive to long-range order (X-ray, neutron, or electron diffraction) and those sensitive to short-range order (X-ray absorption spectroscopy). In this study, we probed the details of the local environment in AMoO3 perovskites (A = Ca, Sr, Ba) using extended X-ray absorption fine structure (EXAFS) in a wide temperature range (10-300 K). An advanced analysis of the EXAFS spectra within the multiple-scattering formalism using the reverse Monte Carlo method enhanced by an evolutionary algorithm allowed us (i) to extract detailed information on metal-oxygen and metal-metal radial distribution functions, and metal-oxygen-metal and oxygen-metal-oxygen bond angle distribution functions, and (ii) to perform polyhedral analysis. The obtained results demonstrate the strong sensitivity of the EXAFS spectra to the tilting of [MoO6] octahedra induced by the differences in the sizes of alkaline earth metal cations (Ca2+, Sr2+, and Ba2+).

Zero-Overpotential Redox Reactions of Quinone-Based Molecules Confined in Carbon Micropores.

Quinone-based aromatic compounds have been studied as electrode materials for various energy-storage devices. However, the relatively large activation barrier of the charge-transfer process of these redox-active molecules causes sluggish reactions and a decrease in energy efficiency. To lower the activation barrier, aromatic compounds must be strongly adsorbed on the electrode surface, preferably via π-π stacking interactions. Molecules in slit-shaped micropores strongly adsorb on the graphitic walls, thus experiencing unique micropore-confinement properties. In this study, the micropore-confinement effect is extended to the adsorption of quinone-based redox-active molecules in 0.8 nm slit-shaped micropores of activated carbon, which produces a drastic reduction in the activation barrier of the charge-transfer process and creates a zero-overpotential redox reaction. The property originates from the short distance (approximately 0.3 nm) between the quinone molecules and the graphitic wall due to the strong adsorption of the aromatic compound. Our results provide the first demonstration that the micropore-confinement effect can reduce and nearly eliminate the activation barrier of an electrochemical reaction. We also demonstrate the applicability of this approach via the charge/discharge performance of a two-electrode cell. Cells comprising the aromatic compound/activated carbon material as positive and negative electrodes exhibit a greater retention capacity than those without activated carbon. The technique described herein can guide the development of high-performance, rapid charging/discharging electrodes for energy-storage devices such as batteries, supercapacitors, and hybrid devices using organic materials.

Describing the Influence of Ball-milling on the Amorphization of Flubendazole Using the PDF and RMC Methods with X-ray Powder Diffraction Data.

Flubendazole (FBZ) is a poorly water-soluble drug, and different methodologies have been proposed to improve its oral bioavailability. Obtaining the amorphous drug phase is an alternative to improve its water solubility. Several techniques for drug amorphization, such as spray drying, lyophilization, melt quenching, solvent-evaporation, and ball milling, can yield various types of structural disorder and possibly render variations in physicochemical properties. Herein, we focus on evaluating the influence of the ball-milling process on the amorphization of FBZ. The characterization of the average global and local structures before, during, and after the milling process is described by sequential Rietveld refinements, pair distribution function analysis, and the Reverse Monte Carlo method. We show that preserving the local structure (nearest molecules) can be responsible for avoiding the fast structure recrystallization commonly observed when using the solvent-evaporation process for the studied drug.

Phase behaviour of ammonium bromide-d4under high pressure and low temperature; an average and local structure study.

We revisit the pressure-induced order-disorder transition between phases II and IV in ammonium bromide-d4using neutron diffraction measurements to characterise both the average and local structures. We identify a very sluggish transition that does not proceed to full conversion and local structure correlations indicate a slight preference for ammonium cation ordering along ⟨110⟩ crystallographic directions, as pressure is increased. Simultaneous cooling below ambient temperature appears to facilitate the pressure-induced transition. Variable-temperature, ambient-pressure measurements across the IV → III → II transitions show slower conversion than previously observed, and that phase III exhibits metastability above ambient temperature.

Atomic structure of an FeCrMoCBY metallic glass revealed by high energy x-ray diffraction.

Amorphous bulk metallic glasses with the composition Fe48Cr15Mo14C15B6Y2have been of interest due to their special mechanical and electronic properties, including corrosion resistance, high yield-strength, large elasticity, catalytic performance, and soft ferromagnetism. Here, we apply a reverse Monte Carlo technique to unravel the atomic structure of these glasses. The pair-distribution functions for various atomic pairs are computed based on the high-energy x-ray diffraction data we have taken from an amorphous sample. Monte Carlo cycles are used to move the atomic positions until the model reproduces the experimental pair-distribution function. The resulting fitted model is consistent with ourab initiosimulations of the metallic glass. Our study contributes to the understanding of functional properties of Fe-based bulk metallic glasses driven by disorder effects.

Orientational disorder in sulfur hexafluoride: a neutron total scattering and reverse Monte Carlo study.

The orientational disorder in crystalline sulfur hexafluoride, SF6, has been studied using a combination of neutron total scattering and the reverse Monte Carlo method. Analysis of the atomic configurations has shown the extent of the disorder through the evaluation of the S-F bond orientational distribution function, consistent with, but improving upon, the results of earlier neutron powder diffraction data. The correlations between orientations of neighbouring molecules have been studied through analysis of the distributions of F-F distances, showing that nearest-neighbour F-F close contacts are avoided, consistent with previous molecular dynamics simulation results. The results present a new case study of the application of neutron total scattering and the reverse Monte Carlo methods for the study of orientational disorder, where in this instance the disorder arises from orientational frustration rather than from a mismatch of molecular and site symmetries.

Orientational order and phase transitions in deuterated methane: a neutron total scattering and reverse Monte Carlo study.

We report a study of the orientational order and phase transitions in crystalline deuterated methane, carried out using neutron total scattering and the reverse Monte Carlo method. The resultant atomic configurations are consistent with the average structures obtained from Rietveld refinement of the powder diffraction data, but additionally enable us to determine the C-D bond orientational distribution functions (ODF) for the disordered molecules in the high-temperature phase, and for both ordered and disordered molecules in the intermediate-temperature phase. We show that this approach gives more accurate information than can been obtained from fitting a bond ODF to diffraction data. Given the resurgence of interest in orientationally-disordered crystals, we argue that the approach of total scattering with the RMC method provides a unique quantification of orientational order and disorder.

Real-space modeling for complex structures based on small-angle X-ray scattering.

A three-dimensional real-space model has been created for hierarchical materials by matching observed and simulated small-angle X-ray scattering patterns. The simulation is performed by arranging the positions of small primary particles and constructing an aggregate structure in a finite-sized cell. In order to avoid the effect of the finite size of the cell, the cell size is extended to infinity by introducing an asymptotic form of the long-range correlations among the primary particles. As a result, simulations for small-angle X-ray scattering patterns can be performed correctly in the low-wavenumber regime (<0.1 nm-1), allowing the model to handle hundred-nanometre-scale structures composed of primary particles of a few nanometres in size. An aerogel structure was determined using this model, resulting in an excellent match with the experimental scattering pattern. The resultant three-dimensional model can generate cross-sectional images similar to those obtained by transmission electron microscopy, and the calculated pore-size distribution is in accord with that derived from the gas adsorption method.

Unravelling Magnetic Nanochain Formation in Dispersion for In Vivo Applications.

Self-assembly of iron oxide nanoparticles (IONPs) into 1D chains is appealing, because of their biocompatibility and higher mobility compared to 2D/3D assemblies while traversing the circulatory passages and blood vessels for in vivo biomedical applications. In this work, parameters such as size, concentration, composition, and magnetic field, responsible for chain formation of IONPs in a dispersion as opposed to spatially confining substrates, are examined. In particular, the monodisperse 27 nm IONPs synthesized by an extended LaMer mechanism are shown to form chains at 4 mT, which are lengthened with applied field reaching 270 nm at 2.2 T. The chain lengths are completely reversible in field. Using a combination of scattering methods and reverse Monte Carlo simulations the formation of chains is directly visualized. The visualization of real-space IONPs assemblies formed in dispersions presents a novel tool for biomedical researchers. This allows for rapid exploration of the behavior of IONPs in solution in a broad parameter space and unambiguous extraction of ​the parameters of the equilibrium structures. Additionally, it can be extended to study novel assemblies formed by more complex geometries of IONPs.