Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio phasing”

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 OSTI2020

Ab Initio Phase Diagram of Tungsten

The phase diagram of tungsten (W) to a pressure (P) of 2500 GPa is investigated using a comprehensive ab initio approach that includes (i) the calculation of the zero temperature (T) free energies (enthalpies) of different solid structures, (ii) the quantum molecular dynamics simulation of the melting curves of different solid structures, (iii) the derivation of the analytic form for the solid-solid phase transition boundary, and (iv) the simulations of the solidification of liquid W into the final solid states on both sides of the solid-solid phase transition boundary, in order to confirm the corresponding analytic form. There are two solid structures confirmed to be present on the phase diagram of W, the ambient body-centered cubic (bcc) and the high-pressure double hexagonal close-packed (dhcp). At T = 0, the bcc-dhcp transition occurs at 1060 GPa, and the transition boundary has a positive slope dT/dP : the bcc-dhcp-liquid triple point is at (P, T) = (1675 GPa, 23680 K).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab Initio Phase Diagram of Chromium to 2.5 TPa

Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.

equation of state↗
DOE OSTI2020

Structure of disordered TiO 2 phases from ab initio based deep neural network simulations

Amorphous TiO 2 (a-TiO 2 ) is widely used in many fields, ranging from photo-electrochemistry to bio-engineering, hence detailed knowledge of its atomic structure is of scientific and technological interest. Here we use an ab initio-based deep neural network potential (DP) to simulate large scale atomic models of crystalline and disordered TiO 2 with molecular dynamics. Our DP reproduces the structural properties of all (11) TiO 2 crystalline phases, predicts the densities and structure factors of molten and amorphous TiO 2 with only a few percent deviation from experiments, and describes the pressure dependence of the amorphous structure in agreement with recent observations. Furthermore, it can be extended to model additional structures and compositions and can be thus of great value in the study of TiO 2 -based (nano-)materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Macromolecular phasing using diffraction from multiple crystal forms

A phasing algorithm for macromolecular crystallography is proposed that utilizes diffraction data from multiple crystal forms – crystals of the same molecule with different unit-cell packings (different unit-cell parameters or space-group symmetries). The approach is based on the method of iterated projections, starting with no initial phase information. The practicality of the method is demonstrated by simulation using known structures that exist in multiple crystal forms, assuming some information on the molecular envelope and positional relationships between the molecules in the different unit cells. With incorporation of new or existing methods for determination of these parameters, the approach has potential as a method for ab initio phasing.

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

Physics in the Machine: Integrating Physical Knowledge in Autonomous Phase-Mapping

Application of artificial intelligence (AI), and more specifically machine learning, to the physical sciences has expanded significantly over the past decades. In particular, science-informed AI, also known as scientific AI or inductive bias AI, has grown from a focus on data analysis to now controlling experiment design, simulation, execution and analysis in closed-loop autonomous systems. The CAMEO (closed-loop autonomous materials exploration and optimization) algorithm employs scientific AI to address two tasks: learning a material system’s composition-structure relationship and identifying materials compositions with optimal functional properties. By integrating these, accelerated materials screening across compositional phase diagrams was demonstrated, resulting in the discovery of a best-in-class phase change memory material. Key to this success is the ability to guide subsequent measurements to maximize knowledge of the composition-structure relationship, or phase map. In this work we investigate the benefits of incorporating varying levels of prior physical knowledge into CAMEO’s autonomous phase-mapping. This includes the use of ab-initio phase boundary data from the AFLOW repositories, which has been shown to optimize CAMEO’s search when used as a prior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stromataxic Stabilization of a Metastable Layered ScFeO3 Polymorph

Metastable polymorphs--materials with the same stoichiometry as the ground state but a different crystal structure--enable many critical technologies. This work describes the development of a stabilization approach for metastable polymorphs that are difficult to achieve through other stabilization techniques (such as epitaxy or quenching) called stromataxy. Stromataxy is a method based on controlling the precursor structure during the initial stages of material growth to dictate phase formation. To illustrate this approach, we controlled the atomic layering of the precursors of ScFeO3 and stabilized the metastable P63cm phase, under conditions that previously led to the ground-state Ia3¯ bixbyite phase. Ab initio mechanistic calculations highlight the importance of the variable oxidation state of Fe and the layer stability during layer-by-layer growth. The broad applicability of a stromataxy approach was demonstrated by stabilizing this metastable phase on substrates that have previously been shown to stabilize other polymorphs under continuous growth. Stromataxy is shown as a viable option for accessing polymorphs that are close in energy, difficult to differentiate by strain, or that lack a well epitaxially matched substrate.

calculations↗
DOE OSTI2024

Smaller Is Better: The Case for Lower-Order Iodoplumbate Species Dominating MAPbI 3 /Dimethylformamide Solutions

Here, using complementary experimental measurements and computational predictions of spectroscopic measurements (EXAFS, XANES, and UV–vis), we have determined the identity of the most stable iodoplumbate species in dilute lead halide perovskite precursor solutions. We have determined which species are most likely to be thermodynamically stable compared to others that are unstable or metastable. Condensed phase ab initio models were constructed, and the resulting ensembles were used to directly compare the computed signals to the experimental results of the EXAFS, XANES, and UV–vis spectra of PbI 2 :MAI in DMF. The results of this study suggest that only Pb 2+ , PbI + , and PbI 2 are dominant in the dilute lead perovskite precursor solutions as thermodynamically stable entities. Our interpretation of the relative stability of iodoplumbate species in solution, based on an analysis of EXAFS and XANES spectra, provides critically important new insight into the species most likely to be responsible for crystal nucleation and growth in these materials. This insight will have a significant consequence on the broad scientific community and will necessitate the reinterpretation of peaks in the UV–vis spectra of lead halide perovskite precursor solutions.

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

Determining Catalytically Relevant Surfaces through Coverage-Dependent Lattice Gas Models: Carbon Adsorption on Fe(100)

Here, we have quantified the C–C lateral interactions on Fe(100) using a density functional theory (DFT)-parameterized lattice gas cluster expansion (LG CE) model trained using 265 unique configurations spanning a C coverage from 0 to 1 monolayer (ML). Our LG CE model shows high predictive accuracy with a leave-multiple-out cross-validation score of 10.2 and 16.6 meV/site for systems with and without the top two layers of Fe atoms fixed, respectively. Electronic ground-state structures identified from the lattice gas model (including the structures at 0 and 1 monolayers) were further used to generate ab initio phase diagrams under a range of temperatures and pressures. At low temperatures (<400 K), we found that the 1.0 monolayer structure is dominant, whereas at higher temperatures (>500 K), the 0.88 ML structure is most likely to form on the Fe surface. Interestingly, our model identified a c (2 × 2) ordered structure at 1/2 ML, which correlates well with previous DFT studies for carbon adsorption on iron surfaces and matches with the experimentally observed low-energy electron diffraction structure. Overall, the DFT-parameterized energies for the C/Fe system including effects of coverage and configurational space can further help in developing multiscale models for various heterogeneous reactions involving C–C and C–Fe interactions.

08 HYDROGEN↗
DOE OSTI2025

Multifunctional Catalysts for the Tandem Reactions of Oxygenates

Industrially-relevant catalytic reactions rarely consist of a simple sequence of elementary steps. Moreover, kinetic coupling of multiple reactions on a catalyst surface is highly desired for process intensification and improved energy efficiency for large scale chemical transformations. The shifting landscape of hydrocarbon chemical feedstocks in the US also motivates research on the selective conversion of more complex molecules. One desirable type of catalytic reaction is the reduction of carboxylic acids that are produced from biomass feedstocks to their corresponding alcohols. The proposed research explores the fundamental importance of hydrogen spillover on a multifunctional catalyst for carboxylic acid reduction with H 2 composed of metal particles coupled to metal oxide particles. Recent work has demonstrated the excellent performance of supported tungsten oxide clusters for carboxylic acid reduction, but only after they are promoted with a late transition metal such as palladium. Elucidating the active state of the catalyst and the associated reaction mechanism for acid reduction on that active state are the overall goals of the proposed project and successful completion will enable future design of efficient multifunctional catalysts. The critically important role of the metal promoter is hypothesized to be its ability to dissociate H2 and spillover atomic H to the support. Although hydrogen spillover is a well-recognized phenomenon in catalysis, its role in both catalyst activation and catalytic turnover are still unresolved. The study combined materials synthesis, characterization, reactivity testing, and molecular simulations, to explore the effect of hydrogen chemical potential on the formation of the active catalytic sites and on the steady state catalytic reduction of carboxylic acid. Varying the hydrogen chemical potential through modification of the gas conditions, support composition, and metal loading to modulated the structure and catalytic performance of the reducible metal oxide. Dual function catalysts containing supported Pd and WO x species co-located on a non-reducible carrier (silica) and a reducible carrier (titania) were synthesized and characterized by electron microscopy, temperature-programmed reduction, and chemisorption. Spectroscopic methods such as X-ray absorption and UV-vis were also used to evaluate the catalysts, which were used in the reduction of propionic acid to aldehyde and alcohol. Quantum chemical calculations, including ab initio phase diagrams provided molecular insights into the H spillover phenomenon.

09 BIOMASS FUELS↗
DOE OSTI2021

Tailored computational approaches to interrogate heavy element chemistry and structure in condensed phase

In this chapter we are presenting a brief review of the challenges encountered in the study of 4f and 5f block elements in the condensed phase. Their recovery, use in molten salt reactors and other interesting applications necessitate the use of molecular dynamics and large-scale models that take into account both the electronic structure and relativistic corrections. Sampling the multitude of electronic and atomic configurational states is at the heart of reliable predictions of structure, reactivity, dynamics and transport of heavy metals in complex environments. We present three examples that combine lanthanide elements with large scale models: i) our recent developments of a versatile adaptive learning method that enables global optimization in high dimensional spaces, ii) results of computed pKa values of lanthanide aqua complexes, and iii) structure and computed EXAFS of heavy elements in molten salts. The latter two employ our recently-optimized lanthanide pseudo-potentials and companion basis sets for condensed phase ab initio molecular dynamics.

Nguyen, Manh Thuong↗
DOE OSTI2021

Ab-initio predictions of phase stability, electronic structure, and optical properties of (0001)-MAX surfaces in M 2 AC (M = Cr, Zr, Hf; A = Al, Ga)

In this work, we report MAX phases' surface properties, which are essential for thin-film technology due to their excellent resistance to high-temperature oxidation, corrosion, and wear. The surface stability, electronic, and optical properties of 0001-surfaces in M 2 AC (M = Zr, Hf, Cr; A = Al, Ga) are investigated and compared with their bulk counterparts. The interplay between chemical bonding and charge distribution is discussed from electronic structure, including the Fermi surfaces. Four possible (0001)-terminated surfaces are considered by breaking M - C and M-A bonds in which cleavage energy of M - C is higher than M-A. The Cr–Al bond in Cr 2 AlC is stronger than other M-A bonds. The charge density of valance A-p electrons redistributes in the surface area, distinct from that of the bulk. The A- and M(C)-terminated (0001)-surfaces are the most stable and energetically favorable terminations due to lower surface energies. The optical properties of the most stable (0001)-surfaces were also investigated to understand the dielectric and photoconductive behavior in the (0001)-terminated surfaces of M 2 AC.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crystal structure of the $τ_{11}$ Al 4 Fe 1.7 Si phase from neutron diffraction and ab initio calculations

The intermetallic τ 11 Al 4 Fe 1.7 Si phase is of interest for high-temperature structural application due to its combination of low density and high strength. We determine the crystal structure of the τ 11 phase through a combination of powder neutron diffraction and density functional theory calculations. Using Pawley and Rietveld refinements of the neutron diffraction data provides an initial crystal structure model. Since Al and Si have nearly identical neutron scattering lengths, we use density-functional calculations to determine their preferred site occupations. The τ 11 phase exhibits a hexagonal crystal structure with space group P6 3 /mmc and lattice parameters of a = 7.478 Å and c = 7.472 Å. The structure comprises five Wyckoff positions; Al occupies the 6h and 12k sites, Fe the 2a and 6h sites, and Si the 2a sites. Here we observe site disorder and partial occupancies on all sites with a large fraction of 80% Fe vacancies on the 2d sites, indicating an entropic stabilization of the τ 11 phase at high temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab initio Geochemistry of Hydrous Phases (Final Report)

Hydrous phases are among the most important Earth components. They are important for a broad suite of Earth processes, including the origin of life. From both thermodynamic and structural perspectives, however, they represent some of the most complex naturally occurring materials: their bonding often includes a combination of covalent solid, weak ionic, van der Waals, and hydrogen bonding, all within large unit cells. Most are solid solutions, and many are prone to variations in layer packing. Thermodynamic modeling of these materials is fundamental for understanding present and past natural processes, including those that shaped—and continue to shape—the structure and evolution of our planet. Yet, the thermodynamic properties of these materials at appropriate conditions are challenging to measure. To make significant progress and attain a deep understanding of these materials requires an atomistic theoretical approach.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗