Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗
DOE OSTI2023

Mass spectrometry and spectroscopic characterization of a tetrameric photosystem I supercomplex from Leptolyngbya ohadii , a desiccation-tolerant cyanobacterium

Cyanobacteria inhabiting desert biological soil crusts face the harsh conditions of the desert. They evolved a suite of strategies toward desiccation-hydration cycles mixed with high light irradiations, etc. In this study we purified and characterized the structure and function of Photosystem I (PSI) from Leptolyngbya ohadii, a desiccation-tolerant desert cyanobacterium. We discovered that PSI forms tetrameric (PSI-Tet) aggregate. We investigated it by using sucrose density gradient centrifugation, clear native PAGE, high performance liquid chromatography, mass spectrometry (MS), time-resolved fluorescence (TRF) and time-resolved transient absorption (TA) spectroscopy. MS analysis identified the presence of two PsaB and two PsaL proteins in PSI-Tet and uniquely revealed that PsaLs are N-terminally acetylated in contrast to non-modified PsaL in the trimeric PSI from Synechocystis sp. PCC 6803. Chlorophyll (Chl) a fluorescence decay profiles of the PSI-Tet performed at 77 K revealed two emission bands at ∼690 nm and 725 nm with the former appearing only at early delay time. The main fluorescence emission peak, associated with emission from the low energy Chls a, decays within a few nanoseconds. TA studies demonstrated that the 725 nm emission band is associated with low energy Chls a with absorption band clearly resolved at ∼710 nm at 77 K. As a result, our work suggests that the heterogenous composition of PsaBs and PsaL in PSI-Tet is related with the adaptation mechanisms needed to cope with stressful conditions under which this bacterium naturally grows.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Topological Dirac Semimetal Phase in Bismuth Based Anode Materials for Sodium-Ion Batteries

Bismuth has recently attracted interest in connection with Na-ion battery anodes due to its high volumetric capacity. It reacts with Na to form Na 3 Bi which is a prototypical Dirac semimetal with a nontrivial electronic structure. Density-functional-theory based first-principles calculations are playing a key role in understanding the fascinating electronic structure of Na 3 Bi and other topological materials. In particular, the strongly-constrained-and-appropriately-normed (SCAN) meta-generalized-gradient-approximation (meta-GGA) has shown significant improvement over the widely used generalized-gradient-approximation (GGA) scheme in capturing energetic, structural, and electronic properties of many classes of materials. Here, we discuss the electronic structure of Na 3 Bi within the SCAN framework and show that the resulting Fermi velocities and s-band shift around the Γ point are in better agreement with experiments than the corresponding GGA predictions. SCAN yields a purely spin-orbit-coupling (SOC) driven Dirac semimetal state in Na 3 Bi in contrast with the earlier GGA results. Our analysis reveals the presence of a topological phase transition from the Dirac semimetal to a trivial band insulator phase in Na 3 Bi x Sb 1-x alloys as the strength of the SOC varies with Sb content, and gives insight into the role of the SOC in modulating conduction properties of Na 3 Bi.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2020

Characterizing Peatland Microtopography Using Gradient and Microform-Based Approaches

Abstract Peatlands represent an important component of the global carbon cycle, storing 180–621 Gt of carbon (C). Small-scale spatial variations in elevation, frequently referred to as microtopography, influence ecological processes associated with the peatland C cycle, including Sphagnum photosynthesis and methane flux. Microtopography can be characterized with measures of topographic variability and by using conceptual classes (microforms) linked to function: most commonly hummocks and hollows. However, the criteria used to define these conceptual classes are often poorly described, if at all, and vary between studies. Such inconsistencies compel development of explicit quantitative methods to classify microforms. Furthermore, gradient-based characterizations that describe spatial variability without the use of microforms are lacking in the literature. Therefore, the objectives of this study were to (1) calculate peatland microtopographical elevation gradients and measures of spatial variability, (2) develop three microform classification methods intended for specific purposes, and (3) evaluate and contrast classification methods. Our results suggest that at spatial scales much larger than microforms, elevation distributions are unimodal and are well approximated with parametric probability density functions. Results from classifications were variable between methods and years and exhibited significant differences in mean hollow areal coverages of a raised ombrotrophic bog. Our results suggest that the conceptualization and classification of microforms can significantly influence microtopographic structural metrics. The three explicit methods for microform classification described here may be used and built upon for future applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

From Molecules to Solids: A vdW-DF-C09 Case Study of the Mercury Dihalides

The mercury dihalides show a remarkable diversity in the structural preferences in their minimum energy structure types, spanning molecular to strongly bound ionic solids. A challenge in the development of density functional methods for extended systems is to arrive at strategies that serve equally well such a broad range of bonding modes or structural preferences. The chemical bonding and the stabilities of mercury dihalides and the general utility and reliability of the van der Waals density functional with C09 exchange (vdW-DF-C09) in predicting or describing the energetics and structural preferences in these metal dihalides is examined. Here, we show that, in contrast with the uncorrected generalized gradient approximation of the Perdew-Burke-Erzenhoff (PBE) exchange-correlation functional, qualitative and quantitative patterns in the bonding of the mercury dihalide solids are well reproduced with vdW-DF-C09 for the full series of HgX 2 systems for X = F, Cl, Br, and I. The possible existence of a low-temperature cotunnite polymorph for HgF 2 and PbF 2 is posited.

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

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I↗
DOE OSTI2024

Training quantum neural networks using the quantum information bottleneck method

Abstract We provide in this paper a concrete method for training a quantum neural network to maximize the relevant information about a property that is transmitted through the network. This is significant because it gives an operationally well founded quantity to optimize when training autoencoders for problems where the inputs and outputs are fully quantum. We provide a rigorous algorithm for computing the value of the quantum information bottleneck quantity within error ε that requires O ( log 2 ⁡ ( 1 / ϵ ) + 1 / δ 2 ) queries to a purification of the input density operator if its spectrum is supported on { 0 } ⋃ [ δ , 1 − δ ] for δ > 0 and the kernels of the relevant density matrices are disjoint. We further provide algorithms for estimating the derivatives of the QIB function, showing that quantum neural networks can be trained efficiently using the QIB quantity given that the number of gradient steps required is polynomial.

Çatlı, Ahmet Burak (ORCID:0000000152294141)↗
DOE OSTI2020

Meta-GGA performance in solids at almost GGA cost

A recent modification, r 2 SCAN, of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA) exchange-correlation functional mostly eliminates numerical instabilities and attendant integration grid sensitivities exhibited by SCAN. Here, we show that the successful deorbitalization of SCAN to SCAN-L (SCAN with density Laplacian dependence) carries over directly to yield r 2 SCAN-L. A major benefit is that the high iteration counts that hindered the use of SCAN-L are eliminated in r 2 SCAN-L. It therefore is a computationally much faster meta-GGA than its orbital-dependent antecedent. In this work, validation data for molecular heats of formation, bond lengths, and vibration frequencies (G3/99X, T96-R, and T82-F test sets, respectively) and on lattice constants, and cohesive energies (for 55 solids) and bulk moduli (for 40 solids) are provided. In addition, we show that the overmagnetization of bcc Fe, hcp Co, and fcc Ni persists in r 2 SCAN but does not appear in r 2 SCAN-L. Distinct from SCAN, both r 2 SCAN and r 2 SCAN-L give the correct nonmagnetic ground state for bcc V.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Machine-learning and first-principles investigation of lightweight medium-entropy alloys for hydrogen-storage applications

The transition to a low-carbon economy demands efficient and sustainable energy-storage solutions, with hydrogen emerging as a promising clean-energy carrier and with metal hydrides recognized for their hydrogen-storage capacity. Here, we leverage machine learning (ML) to predict hydrogen-to-metal (H/M) ratios and solution energy by incorporating thermodynamic parameters and local lattice distortion (LLD) as key features. Our best-performing ML model provides improvements to H/M ratios and solution energies over a broad class of medium-entripy alloys (easily extendable to multi-principal-element alloys), such as Ti–Nb-X (X = Mo, Cr, Hf, Ta, V, Zr) and Co–Ni-X (X = Al, Mg, V). Ti–Nb–Mo alloys reveal compositional effects in H-storage behavior, in particular Ti, Nb, and V enhance H-storage capacity, while Mo reduces H/M and hydrogen weight percent by 40–50 %. We attributed results in molybdenum-rich alloys to slow hydrogen kinetics, as validated by our pressure-composition-temperature (PCT) isotherm experiments on pure Ti and Ti 5 Mo 95 alloys. Density functional theory (DFT) and molecular dynamics (MD) simulations also confirm that Ti and Nb promote H diffusion, whereas Mo hinders it, highlighting the interplay between electronic structure, lattice distortions, and hydrogen uptake. Notably, our Gradient Boosting Regression model identifies LLD as a critical factor in H/M predictions. Here, to aid material selection, we present two periodic tables illustrating elemental effects on (a) H 2 wt% and (b) solution energy, derived from ML, and provide a reference for identifying alloying elements that enhance hydrogen solubility and storage.

08 HYDROGEN↗
DOE OSTI2023

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

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

Thermodynamics of Uranium Tri-Iodide from Density-Functional Theory

Density-functional theory (DFT) is employed to investigate the thermodynamic and ground-state properties of bulk uranium tri-iodide, UI3. The theory is fully relativistic and electron correlations, beyond the DFT and generalized gradient approximation, are addressed with orbital polarization. The electronic structure indicates anti-ferromagnetism, in agreement with neutron diffraction, with band gaps and a non-metallic system. Furthermore, the formation energy, atomic volume, crystal structure, and heat capacity are calculated in reasonable agreement with experiments, whereas for the elastic constants experimental data are unavailable for comparison. The thermodynamical properties are modeled within a quasi-harmonic approximation and the heat capacity and Gibbs free energy as functions of temperature agree with available calculation of phase diagram (CALPHAD) thermodynamic assessment of the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

TaTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-titanium (Ta-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2024

VZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys vanadium-zirconium (V-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements V and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗
DOE OSTI2023

PBE-GGA predicts the B8↔B2 phase boundary of FeO at Earth’s core conditions

FeO is a crucial component of the Earth’s core, and its thermodynamic properties are essential to developing more accurate core models. It is also a notorious correlated insulator in the NaCl-type (B1) phase at ambient conditions. It undergoes two polymorphic transitions at 300 K before it becomes metallic in the NiAs-type (B8) structure at ~100 GPa. Although its phase diagram is not fully mapped, it is well established that the B8 phase transforms to the CsCl-type (B2) phase at core pressures and temperatures. Here, we report a successful ab initio calculation of the B8↔B2 phase boundary in FeO at Earth’s core pressures. We show that fully anharmonic free energies computed with the Perdew–Burke–Ernzerhof-generalized gradient approximation coupled with thermal electronic excitations reproduce the experimental phase boundary within uncertainties at P > 255 GPa, including the largely negative Clapeyron slope of –52 MPa/K. This study validates the applicability of a standard density functional theory functional to FeO under Earth’s core conditions and demonstrates the theoretical framework that enables complex predictive studies of this region.

58 GEOSCIENCES↗