Search NASA⌕ Search

SEARCH · Search NASA

Results for “structure refinement”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Integrating machine learning interatomic potentials with hybrid reverse Monte Carlo structure refinements in RMCProfile

Structure refinement with reverse Monte Carlo (RMC) is a powerful tool for interpreting experimental diffraction data. To ensure that the under-constrained RMC algorithm yields reasonable results, the hybrid RMC approach applies interatomic potentials to obtain solutions that are both physically sensible and in agreement with experiment. To expand the range of materials that can be studied with hybrid RMC, we have implemented a new interatomic potential constraint in RMCProfile that grants flexibility to apply potentials supported by the Large-scale Atomic/Molecular Massively Parallel Simulator ( LAMMPS ) molecular dynamics code. This includes machine learning interatomic potentials, which provide a pathway to applying hybrid RMC to materials without currently available interatomic potentials. To this end, we present a methodology to use RMC to train machine learning interatomic potentials for hybrid RMC applications.

Cuillier, Paul↗
DOE OSTI2023

Atomic protein structure refinement using all-atom graph representations and SE(3)-equivariant graph transformer

Abstract Motivation The state-of-art protein structure prediction methods such as AlphaFold are being widely used to predict structures of uncharacterized proteins in biomedical research. There is a significant need to further improve the quality and nativeness of the predicted structures to enhance their usability. In this work, we develop ATOMRefine, a deep learning-based, end-to-end, all-atom protein structural model refinement method. It uses a SE(3)-equivariant graph transformer network to directly refine protein atomic coordinates in a predicted tertiary structure represented as a molecular graph. Results The method is first trained and tested on the structural models in AlphaFoldDB whose experimental structures are known, and then blindly tested on 69 CASP14 regular targets and 7 CASP14 refinement targets. ATOMRefine improves the quality of both backbone atoms and all-atom conformation of the initial structural models generated by AlphaFold. It also performs better than two state-of-the-art refinement methods in multiple evaluation metrics including an all-atom model quality score—the MolProbity score based on the analysis of all-atom contacts, bond length, atom clashes, torsion angles, and side-chain rotamers. As ATOMRefine can refine a protein structure quickly, it provides a viable, fast solution for improving protein geometry and fixing structural errors of predicted structures through direct coordinate refinement. Availability and implementation The source code of ATOMRefine is available in the GitHub repository (https://github.com/BioinfoMachineLearning/ATOMRefine). All the required data for training and testing are available at https://doi.org/10.5281/zenodo.6944368.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Structure refinement and anisotropic atomic displacement parameters of 1M Illite: Rietveld and pair distribution function analysis using synchrotron X-ray radiation

Illite, a widespread clay mineral, plays a pivotal role in geological processes, notably as an indicator in diagenetic and hydrothermal alteration environments, and possesses significant industrial relevance in applications including ceramics, construction and catalysis. However, challenges including its nanoscale crystallinity, structural disorder and frequent interstratification with other clay minerals have hindered detailed structural characterization using conventional X-ray diffraction (XRD) techniques. This study employs integrated synchrotron XRD and pair distribution function (PDF) analysis to elucidate the crystal structure of the 1M illite polytype, yielding the first determination of its anisotropic atomic displacement parameters (U aniso ). TheseU aniso parameters provide critical insights into atomic dynamics and static disorder within the structure, enabling a more refined understanding of structure–property relationships. This integrated approach, combining synchrotron XRD, Rietveld refinement and PDF analysis, yields a comprehensive structural characterization, capturing both average crystallographic and local atomic arrangements. Considering illite's widespread geological occurrence and industrial importance, this high-precision structural dataset, especially the determinedU aniso values, provides a crucial benchmark for future modeling and simulation efforts targeting accurate prediction of its physicochemical behavior.

Chemistry↗
DOE OSTI2020

On single-crystal neutron-diffraction in DACs: quantitative structure refinement of light elements on SNAP and TOPAZ

Quantitative single crystal neutron-diffraction in diamond anvil cells has so far been limited by the neutron flux available at the various neutron sources. As a result, highly precise measurements of the exact position of light elements have not been possible preventing, for example, structural studies of hydrogen and hydrogen bonds under pressure. Here we report experiments carried out on SNAP at the Spallation Neutron Source (ORNL, TN, USA) to explore the possibility and current limits of such studies. Furthermore, we benchmarked the obtained data quality with reference experiments carried out on TOPAZ, a dedicated single-crystal instrument.We show that measuring single-crystal diffraction intensities on SNAP is possible to such a precision that we are able to resolve the hydrogen bonds in potassium dideuterium phosphate (DKDP) as well as in ice VI.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗
DOE OSTI2023

Symmetry-mode analysis for local structure investigations using pair distribution function data

Symmetry-adapted distortion modes provide a natural way of describing distorted structures derived from higher-symmetry parent phases. Structural refinements using symmetry-mode amplitudes as fit variables have been used for at least ten years in Rietveld refinements of the average crystal structure from diffraction data; more recently, this approach has also been used for investigations of the local structure using real-space pair distribution function (PDF) data. Here, the value of performing symmetry-mode fits to PDF data is further demonstrated through the successful application of this method to two topical materials: TiSe2, where a subtle but long-range structural distortion driven by the formation of a charge-density wave is detected, and MnTe, where a large but highly localized structural distortion is characterized in terms of symmetry-lowering displacements of the Te atoms. Here, the analysis is performed using fully open-source code within the DiffPy framework via two packages developed for this work: isopydistort, which provides a scriptable interface to the ISODISTORT web application for group theoretical calculations, and isopytools, which converts the ISODISTORT output into a DiffPy-compatible format for subsequent fitting and analysis. These developments expand the potential impact of symmetry-adapted PDF analysis by enabling high-throughput analysis and removing the need for any commercial software.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Refining perovskite structures to pair distribution function data using collective Glazer modes as a basis

Structural modelling of octahedral tilts in perovskites is typically carried out using the symmetry constraints of the resulting space group. In most cases, this introduces more degrees of freedom than those strictly necessary to describe only the octahedral tilts. It can therefore be a challenge to disentangle the octahedral tilts from other structural distortions such as cation displacements and octahedral distortions. This paper reports the development of constraints for modelling pure octahedral tilts and implementation of the constraints in diffpy-CMI, a powerful package to analyse pair distribution function (PDF) data. The model in the program allows features in the PDF that come from rigid tilts to be separated from non-rigid relaxations, providing an intuitive picture of the tilting. The model has many fewer refinable variables than the unconstrained space group fits and provides robust and stable refinements of the tilt components. It further demonstrates the use of the model on the canonical tilted perovskite CaTiO 3 which has the known Glazer tilt system α + β – β – . The Glazer model fits comparably to the corresponding space-group model Pnma below r = 14 Å and becomes progressively worse than the space-group model at higher r due to non-rigid distortions in the real material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structure Analysis of Natural Wangdaodeite—LiNbO3-Type FeTiO3

This paper reports the first structure refinement of natural wangdaodeite, LiNbO3-type FeTiO3 from the Ries impact structure. Wangdaodeite occurs together with recrystallized ilmenite clasts in shock melt veins which have experienced peak shock pressures of between 17 and 22 GPa. Comparison of natural and synthetic wangdaodeite points toward a correlation between the distortion of ferrate- and titanate-polyhedra and the c/a ratio of the unit cell. The Raman spectrum of wangdaodeite is calculated based on the refined structure. Comparison to the reported spectrum of the type-material shows that the Raman peak at 738–740 cm−1 is indicative for this phase, whereas other features in type-wangdaodeite are tentatively assigned to disordered ilmenite.

58 GEOSCIENCES↗
DOE OSTI2024

Transferable Hirshfeld atom model for rapid evaluation of aspherical atomic form factors

Form factors based on aspherical models of atomic electron density have brought great improvement in the accuracies of hydrogen atom parameters derived from X-ray crystal structure refinement. Today, two main groups of such models are available, the banks of transferable atomic densities parametrized using the Hansen–Coppens multipole model which allows for rapid evaluation of atomic form factors and Hirshfeld atom refinement (HAR)-related methods which are usually more accurate but also slower. In this work, a model that combines the ideas utilized in the two approaches is tested. It uses atomic electron densities based on Hirshfeld partitions of electron densities, which are precalculated and stored in a databank. This model was also applied during the refinement of the structures of five small molecules. A comparison of the resulting hydrogen atom parameters with those derived from neutron diffraction data indicates that they are more accurate than those obtained with the Hansen–Coppens based databank, and only slightly less accurate than those obtained with a version of HAR that neglects the crystal environment. The advantage of using HAR becomes more noticeable when the effects of the environment are included. To speed up calculations, atomic densities were represented by multipole expansion with spherical harmonics up to l = 7, which used numerical radial functions (a different approach to that applied in the Hansen–Coppens model). Calculations of atomic form factors for the small protein crambin (at 0.73 Å resolution) took only 68 s using 12 CPU cores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Structure and properties of the Sr 2 In 1-x Sn x SbO 6 double perovskite

A series of n -type oxide double perovskite semiconductors, Sr 2 In 1-x Sn x SbO 6 (0 ≤ x ≤ 0.3) has been synthesized; Sn 4+ partially substitutes for In 3+ . 121 Sb and 119 Sn Mössbauer spectroscopy are employed to investigate the B-site cation ordering because this issue cannot be resolved by conventional diffraction techniques alone. Rigid ordering between In 3+ /Sn 4+ and Sb 5+ sites is revealed by the spectroscopic method, and hence in combination with the structural parameters extracted from the XRD structural refinements, the crystallographic structure of this series of compounds is depicted. Furthermore, the temperature dependent magnetic susceptibilities, band gaps, and carrier type are characterized, and the calculated band structure is presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Multilevel atomic structural model for interstratified opal materials

The structure of opal has long fascinated scientists. It occurs in a number of structural states, ranging from amorphous to exhibiting features of stacking disorder. Opal-CT, where C and T signify cristobalite- and tridymite-like interstratification, represents an important link in the length scales between amorphous and crystalline states. However, details about local atomic (dis)order and arrangements extending to long-range stacking faults in opal polymorphs remain incompletely understood. Here, a multilevel modeling approach is reported that considers stacking states in correlation with the abundance of C and T segments as a high-level structural parameter (i.e. not each atom). Optimization accounting for inter-tetrahedral bond lengths and angles and the regularity of the silicate tetrahedra is included as lower levels of structural parameters. Together, a set of parameters with both coarse-grained and atomistic features for different levels of structural details is refined. Structural disorder at the ~10–100 Å distance scale is evaluated using experimental pair distribution function and diffraction datasets, comparing peak intensities, widths and asymmetry. Here this work presents a complete multilevel structural description of natural opal-CT and explains many of the unusual features observed in X-ray powder diffraction patterns. This modeling approach can be adopted generally for analyzing layered materials and their assembly into 3D structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Using Complementary Methods of Synchrotron Radiation Powder Diffraction and Pair Distribution Function to Refine Crystal Structures with High Quality Parameters—A Review

Determination of the atomic-scale structures of certain fine-grained minerals using single-crystal X-ray diffraction (XRD) has been challenging because they commonly occur as submicron and nanocrystals in the geological environment. Synchrotron powder diffraction and scattering techniques are useful complementary methods for studying this type of minerals. In this review, we discussed three example studies investigated by combined methods of synchrotron radiation XRD and pair distribution function (PDF) techniques: (1) low-temperature cristobalite; (2) kaolinite; and (3) vernadite. Powder XRD is useful to determine the average structure including unit-cell parameters, fractional atomic coordinates, occupancies and isotropic atomic displacement parameters. X-ray/Neutron PDF methods are sensitive to study the local structure with anisotropic atomic displacement parameters (ADP). The results and case studies suggest that the crystal structure and high-quality ADP values can be obtained using the combined methods. The method can be useful to characterize crystals and minerals that are not suitable for single-crystal XRD.

58 GEOSCIENCES↗
DOE OSTI2021

Experimental Verification of Ir 5d Orbital States and Atomic Structures in Highly Active Amorphous Iridium Oxide Catalysts

In iridium oxide catalysts, the electronic states whose energies are in the range of energetics and charge transfer kinetics of the oxygen evolution reaction (OER) originate from the Ir 5d orbital states. However, the understanding of the atomic structures and orbital states underlying catalytic reactivity in amorphous iridium oxide oxygen evolving catalysts (Ir-OECs) is incomplete compared to that of crystalline oxides, owing to a lack of direct experimental verification. Here, we present experimental approaches using resonant inelastic X-ray scattering (RIXS) to directly access Ir 5d orbital excitations at the Ir L 3 edge and atomic pair distribution function (PDF) measurements to characterize electronic and coordination structures at the atomic scale. The so-called iridium blue layer (IrBL) and IrO x were formed from the organometallic precursor complex [Cp*Ir(H 2 O) 3 ]SO 4 and the inorganic precursor IrCl 3 , respectively. Ex situ IrBL and IrO x films for RIXS and PDF measurements were prepared by conditioning electrodeposited films at a low voltage. The incident energy RIXS profile of IrO x exhibited extra weak resonantly enhanced excitation below 2 eV energy loss. The feature was clearly different from a single high-energy excitation above 3 eV of IrBL related to the interband transition between π- and σ-antibonding states. The atomic structure refinement based on PDF measurements revealed the atomic structure domains to have edge- and corner-shared IrO 6 octahedra with trigonal-type distortion. Furthermore, density functional theory calculations guided by the refined atomic structures shed light on the electronic structure corresponding to experimental results, including insulating and metallic phases in ex situ IrBL and IrO x films, respectively. Our study establishes different Ir 5d orbital states and atomic structures in two amorphous Ir oxide OER catalysts in their reduction states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Octahedral oxide glass network in ambient pressure neodymium titanate

Abstract Rare-earth titanates form very fragile liquids that can be made into glasses with useful optical properties. We investigate the atomic structure of 83TiO 2 -17Nd 2 O 3 glass using pair distribution function (PDF) analysis of X-ray and neutron diffraction with double isotope substitutions for both Ti and Nd. Six total structure factors are analyzed (5 neutron + 1 X-ray) to obtain complementary sensitivities to O and Ti/Nd scattering, and an empirical potential structure refinement (EPSR) provides a structural model consistent with the experimental measurements. Glass density is estimated as 4.72(13) g cm −3 , consistent with direct measurements. The EPSR model indicates nearest neighbor interactions for Ti-O at $$\overline{r}_{TiO}$$ r ¯ TiO = 1.984(11) Å with coordination of $$n_{TiO}$$ n TiO = 5.72(6) and for Nd-O at $$\overline{r}_{NdO}$$ r ¯ NdO = 2.598(22) Å with coordination of $$n_{NdO}$$ n NdO = 7.70(26), in reasonable agreement with neutron first order difference functions for Ti and Nd. The titanate glass network comprises a mixture of distorted Ti-O 5 and Ti-O 6 polyhedra connected via 71% corner-sharing and 23% edge-sharing. The O-Ti coordination environments include 15% nonbridging O-Ti 1 , 51% bridging O-Ti 2 , and 32% tricluster O-Ti 3 . This structure is highly unusual for oxide glasses melt-quenched at ambient pressure, as it consists of Ti-O x predominantly in octahedral (with nearly no tetrahedral) coordination.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Operando single crystal neutron diffraction reveals insight into the field response mechanisms in the hydrogen-bonded KH 2 PO 4 ferroelectric

The mechanism that facilitates polarization reorientation in KH 2 PO 4 (KDP) was investigated using operando single-crystal neutron diffraction. Diffraction data were measured from a KDP single crystal during the application of alternating electric fields and were then binned into 40 increments to enable field-dependent single-crystal structure refinements. The field-dependent structures are compared with an as-grown crystal to determine how the lattice and atomic sites evolve in response to the applied electric fields. These analyses provide evidence that the reorientation of the macroscopic polarization is facilitated through a cooperative change in hydrogen bonding, which results in the reversal of the spontaneous dipole. In addition, a decrease in secondary extinction near the coercive field indicates that the inversion of the macroscopic polarization is achieved through the nucleation and subsequent growth of new domains.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Structure of Liquid and Glassy Carbamazepine

To enhance the solubility of orally administered pharmaceuticals, liquid capsules or amorphous tablets are often preferred over crystalline drug products. However, little is known regarding the variation in bonding mechanisms between pharmaceutical molecules in their different disordered forms. In this study, liquid and melt-quenched glassy carbamazepine have been studied using high energy X-ray diffraction and modeled using Empirical Potential Structure Refinement. The results show significant structural differences between the liquid and glassy states. The liquid shows a wide range of structures; from isolated molecules, to aromatic ring correlations and NH-O hydrogen bonding. Upon quenching from the liquid to the glass the number of hydrogen bonds per molecule increases by ~50% at the expense of a ~30% decrease in the close contact (non-bonded) carbon-carbon interactions between aromatic rings. During the cooling process, there is an increase in both singly and doubly hydrogen-bonded adjacent molecules. Although hydrogen-bonded dimers found in the crystalline states persist in the glassy state, the absence of a crystalline lattice also allows small, hydrogen-bonded NH-O trimers and tetramers to form. This proposed model for the structure of glassy carbamazepine is consistent with the results from vibrational spectroscopy and nuclear magnetic resonance.

Benmore, Chris J. (ORCID:0000000170077749)↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗