Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Polaron-induced metal-to-insulator transition in vanadium oxides from density functional theory calculations

Vanadium oxides have been extensively studied as phase-change memory units in artificial synapses for neuromorphic computing due to their metal-insulator transitions (MIT) at or near room temperature. Recently, injection of charge carriers into vanadium oxides, e.g., via optically via a heterostructure, has been proposed as an alternative switching mechanism and also potentially as a means to tune the MIT temperature. In this study, we explore the formation of small polarons in the low temperature (LT) insulating phases for V 3 O 5 ,VO 2 , and V 2 O 3 , and the barriers to their migration using density functional theory calculations. We find that V 3 O 5 exhibits very low hole and electron polaron migration barriers (<100 meV) compared to V 2 O 3 and VO 2 , leading to much higher estimated polaronic conductivity. We also link the relative migration barriers to the amount of distortion that has to travel when the polaron migrate from one site to another. Polarons in V 3 O 5 also have smaller binding energies to vanadium and oxygen vacancy defects. Furthermore, these results suggest that the triggering of the MIT via injection of charge carriers are due to the formation of small polarons that can migrate rapidly through the crystal.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

HydraGNN_Predictive_GFM_2026 - Ensemble of predictive graph foundation models for atomistic materials modeling

This release contains data and parameters of HydraGNN-based graph foundation models trained as a result of the work published in the pre-print "Exascale Multi-Task Graph Foundation Models for Imbalanced, Multi-Fidelity Atomistic Data" by M. Lupo Pasini et al. (https://arxiv.org/abs/2604.15380). We jointly train on 16 open first-principles datasets (544+ million structures covering 85+ elements) using a multi-task architecture with per-dataset heads and a scalable ADIOS2/DDStore data pipeline. On Frontier, we execute six large-scale DeepHyper hyperparameter optimization campaigns in FP64 and promote the top-performing message-passing models to sustained 2,048-node training, yielding a PaiNN-based lead model. The version of HydraGNN used to generate the outputs provided in this release is HydraGNN v5.0 (https://github.com/ORNL/HydraGNN/releases/tag/v5.0) The list of datasets used for the training of the graph foundation model is the following: 1) Alexandria [1] 2) ANI1x [2] 3) MPTrj [3] 4) Open Catalyst 2020 (OC20) [4] 5) Open Catalyst 2022 (OC22) [5] 6) Open Catalyst 2025 (OC25) [6] 7) Open Direct ir Capture 2023 (ODAC23) [7] 8) Open Materials 2024 (OMat24) [8] 9) Open Molecules 2025 (OMol25) [9] 10) OMol25-neutral (subset of OMol25 that contains only molecules with zero total charge) 11) OMol25-non-neutral (subset of OMol25 that contains only molecules with non-zero total charge) 12) Open Polymers 2026 (OPoly2026) [10] 13) Nabla2DFT [11] 14) QCML [12] 15) QM7X [reference 13] 16) transition1x [14] Dataset references: [1] J. Schmidt et al., “A dataset of 175k stable and metastable materials calculated with the PBEsol and SCAN functionals,” Scientific Data, vol. 9, p. 64, 2022. [2] J. S. Smith et al., “The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules,” Scientific Data, vol. 7, p. 134, 2020. [Online]. Available: https: //www.nature.com/articles/s41597-020-0473-z [3] A. Jain et al., “Commentary: The Materials Project: A materials genome approach to accelerating materials innovation,” APL Materials, vol. 1, no. 1, p. 011002, 07 2013. [Online]. Available: https://doi.org/10.1063/1.4812323 [4] L. Chanussot et al., “Open catalyst 2020 (oc20) dataset and community challenges,” ACS Catalysis, vol. 11, no. 10, pp. 6059–6072, 2021. [Online]. Available: https://doi.org/10.1021/acscatal.0c04525 [5] K. Tran et al., “Open catalyst 2022 (oc22) dataset and challenges for oxidation electrocatalysts,” ACS Catalysis, vol. 13, no. 5, pp. 3066–3084, 2023. [Online]. Available: https://doi.org/10.1021/acscatal.2c05426 [6] S. J. Sahoo et al., “The open catalyst 2025 (oc25) dataset and models for solid-liquid interfaces,” arXiv preprint arXiv:2509.17862, 2025. [Online]. Available: https://arxiv.org/abs/2509.17862 [7] A. Sriram et al., “The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture,” ACS Central Science, vol. 10, no. 5, pp. 923–941, 2024. [8] L. Barroso-Luque et al., “Open materials 2024 (omat24) inorganic materials dataset and models,” 2024. [Online]. Available: https://arxiv.org/abs/2410.12771 [9] D. S. Levine et al., “The open molecules 2025 (OMol25) dataset, evaluations, and models,” 2025. [Online]. Available: https://arxiv.org/abs/2505.08762 [10] D. S. Levine et al., The open polymers 2026 (OPoly26) dataset and evaluations,” arXiv preprint arXiv:2512.23117, 2025. [Online]. Available: https://arxiv.org/abs/2512.23117 [11] K. Khrabrov et al., “Nabla2dft: A universal quantum chemistry dataset of drug-like molecules and a benchmark for neural network potentials,” in NeurIPS 2024 Datasets and Benchmarks Track, 2024. [Online]. Available: https://openreview.net/forum?id=ElUrNM9U8c [12] S. Ganscha et al., “The QCML dataset, quantum chemistry reference data from 33.5M DFT and 14.7B semi-empirical calculations,” Scientific Data, vol. 12, p. 406, 2025. [13] J. Hoja et al., “QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules,” Scientific Data, vol. 8, p. 43, 2021. [Online]. Available: https://www.nature.com/articles/s41597-021-00812-2 [14] M. Schreiner et al., “Transition1x - a dataset for building generalizable reactive machine learning potentials,” Scientific Data, vol. 9, p. 779, 2022. The folder "datasets_ADIOS2_format" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "datasets_ADIOS2_format" directory contains 2 sub-directories, one for the version "v1" of the datasets and one for the version "v2" of the datasets. The version "v1" of the datasets provides values of the total energy as they are extracted from the original data as it was released by the respective institutions. The version "v2" of the datasets provides values of the energy that have been realigned. The realignment was performed by training a linear regression model that predicts the total energy as a function of the chemical composition of the atomistic structure, and then subtract such prediction from the original value of the total energy. Both folders "v1" and "v2" contain 16 sub-directories, each corresponding to an ADIOS2-formatted dataset The folder "DeepHyper-results" contains the configurational files and model's parameters for all the 186 HPO trials that were successfully completed by the scalable hyperparameter optimization (HPO) runs on Frontier. The content of the folder "DeepHyper-results" I structured as follows: 1) task-list.txt: list of mpnn name, jobid, and deephyper task id 2) gfm_${MPNN}_${JOBID}_0.${TASKID}: run directory with checkpoint files 3) gfm_${MPNN}: deephyper summary directory (*.csv) for each specific MPNN type 4) deephyper-experiment-${JOBID}: output and error logs for each job The file "deephyper-sorted.csv" contains the details of each HydraGNN model built and tested by HPO, obtained by merging the (*.csv) filed from each HPO run executed. Out of all the HPO trials, we selected 10 to continue the training of the respective HydraGNN models. Due to limited computational budget available in the LRN070 allocation we could not complete the training till convergence for all these 10 selected models. The folder "models" contains multiple sub-folders, one per each HydraGNN model trained. Each model sub-folder contains the parameters of each HydraGNN model, with multiple checkpoint-restarts. The list of sub-folders are as follows: 1) multidataset_hpo-BEST1-fp64 2) multidataset_hpo-BEST2-fp64 3) multidataset_hpo-BEST3-fp64 4) multidataset_hpo-BEST4-fp64 5) multidataset_hpo-BEST5-fp64 6) multidataset_hpo-BEST6-fp64 7) multidataset_hpo-BEST7-fp64 8) multidataset_hpo-BEST8-fp64 9) multidataset_hpo-BEST9-fp64 10) multidataset_hpo-BEST10-fp64 Within each one of these folders, additional auxiliary log files are provided with descriptions about how the training proceeded. The lead PaiNN-model is contained inside "multidataset_hpo-BEST6-fp64". The file "mlp_branch_weights" contains the parameters of the multi-layer perceptron (MLP) used to reconcile the predictions of the 16 output decoding heads of the HydragNN architectures. The MLP takes in input the chemical composition of the atomistic structure and predicts averaging weights to linearly mix the predictions of each output decoding head toward consolidating them into a single one. The folder "1.1billion-structure-inference" contains 1.1 billion atomistic structures randomly generated. Each structures is associated with energy and forces predicted with the lead-PaiNN model combined with the MLP model for reconciliation of the multi-branch predictions generated by the 16 output decoding heads. The folder "1.1billion-structure-inference" contains 9,300 (*.tar.gz) subdirectories, one per Frontier compute node used to execute the inference at exascale. Once uncompressed, each (*.tar.gz) subdirectory contains an ADIOS2 (*.bp) file container, where each atomistic structure is stored as a PyTorch-Geometric Data object. The file "export_dataset_environment_variables.sh" contains the environment variables that need to be set before running the HydraGNN code to reproduce the results provided in this dataset release. The code that can be used to load the ADIOS2 files, load HydraGNN models, and run inference is available at: https://github.com/ORNL/HydraGNN/releases/tag/v5.0

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio calculations of third-order elastic coefficients

Third-order elasticity (TOE) theory predicts strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the natural, or zero pressure, reference state. This definition is inconvenient for predictions of SOECs under significant initial stresses. Therefore, when TOE theory is necessary to study the strain dependence of elasticity, the seismological community has resorted to an empirical version of the theory. This study reviews the thermodynamic definition of the third-order elastic tensor and proposes using an “effective” third-order elastic tensor. An explicit expression for the effective third-order elastic tensor is given and verified. Additionally, we extend the ab initio approach to calculate third-order elastic tensors under finite pressure and apply it to two cubic systems, namely, NaCl and MgO. As applications and validations, we evaluate (a) strain-induced changes in SOECs and (b) pressure derivatives of SOECs based on ab initio calculations. Good agreement between third-order elasticity-based predictions and numerically calculated values confirms the validity of our theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Vacancy surface migration mechanisms in dilute nickel-chromium alloys

We investigate the unitary mechanisms related to the surface migration of vacancies in dilute Ni-Cr alloys via first-principle calculations. We survey a complete set of surface and sub-surface migration paths for vacancies near the (100) free surface and calculate the corresponding migration barriers. Our results show that a vacancy migrating towards the free surface will face lower energy barriers to migrate via atomic exchange with a neighboring Cr atom rather than with a Ni atom. Once a vacancy reaches the free surface, it will be trapped there. Our results also reveal that, when a Cr atom sits in the atomic plane just below the free surface, any in-plane vacancy hopping jump that would result in the vacancy sitting on top a subsurface Cr atom is energetically unfavorable. Taken together, these fundamental unitary surface migration mechanisms offer insights into the complex interactions between surface segregation and vacancy migration phenomena in Ni-Cr-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Topological antiferromagnetic semimetal for spintronics: A case study of a layered square-net system EuZnSb 2

Here, we use the first-principles and effective Hamiltonian methods to study the electronic structure and magnetic properties of a recently synthesized layered antiferromagnetic square-net topological semimetal EuZnSb2 [Phys. Rev. Res. 2, 033462 (2020)]. The main message of the paper is that effects of small changes in the band structure produced by the magnetic ordering and changes in the orientation of the Néel vector are amplified in such transport properties as the spin Hall conductivity. We predict that the effects of the broken symmetry introduced by the ordering of the Néel vector, being very weak in the bulk, are pronounced in the surface electronic dispersion, suggesting that surface probes may be more suited to measure them. The coexistence of the magnetism with many other competing phases make this material interesting and possibly useful for quantum spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Li interaction-induced phase transition from black to blue phosphorene

A comprehensive first-principle calculation has been carried out and revealed that at sufficiently high Li concentration and certain well-defined configurations a phase transition from black to blue phosphorene can take place. Blue phosphorene, a newly predicted allotrope of phosphorus, possesses unique crystalline and electronic structure and is a promising candidate, not only for fundamental research but also for electronic and optoelectronic applications. Methods to growth high quality blue phosphorene layers are highly desirable but challenging. Here, a novel kinetic pathway to grow blue phosphorene layers from black phosphorene layers via Li intercalation is proposed based on first principle study. This study pointed out that Li atoms intercalated in black phosphorene could act as ‘catalyst’ in the ‘reactive region’ of the lone pair of P atoms, leading to a P-P bond breaking and subsequently, a local structural transformation from orthorhombic lattice to an assembly of parallel narrow nanoribbons with rhombohedra-like symmetry. During Li deintercalation, these nanoribbons are self-mended and form blue phosphorene layers. The interlayer distance was found 4.60 Å for double layer with AA stacking and 4.13 Å for multilayer with ABC stacking, respectively, indicating a monolayer blue phosphorene can be mechanically exfoliated. Furthermore, this study also points out the possibility of new phases in other systems, where intercalation can lead to an unexpected structural phase transition and even a discovery of novel materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine learning guided discovery of superconducting calcium borocarbides

Pursuit of superconductivity in light-element systems at ambient pressure is of great experimental and theoretical interest. In this work, we combine a machine learning (ML) method with first-principles calculations to efficiently search for the energetically favorable ternary Ca-B-C compounds. Three layered borocarbides (stable CaBC 5 and metastable Ca 2 BC 11 and CaB 3 C 3 ) are predicted to be phonon-mediated superconductors at ambient pressure. The stable CaBC 5 and the low-energy metastable Ca 2 BC 11 (with formation energy only 9.5 meV/atom above the convex hull) have a superconducting T c of 5.2 and 8.9 K, respectively. While the hexagonal CaB 3 C 3 possesses a T c of 26.1 K, it is metastable with formation energy of 153 meV/atom above the convex hull. The ML-guided approach opens up a way for greatly accelerating the discovery of new high-T c superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Split Ga vacancies and the unusually strong anisotropy of positron annihilation spectra in β–Ga 2 O 3

We report a systematic first-principles study on positron annihilation parameters in the β–Ga 2 O 3 lattice and Ga monovacancy defects complemented with orientation-dependent experiments of the Doppler broadening of the positron-electron annihilation. We find that both the β–Ga 2 O 3 lattice and the considered defects exhibit unusually strong anisotropy in their Doppler broadening signals. This anisotropy is associated with low symmetry of the β–Ga 2 O 3 crystal structure that leads to unusual kind of one-dimensional confinement of positrons even in the delocalized state in the lattice. In particular, the split Ga vacancies recently observed by scanning transmission electron microscopy produce unusually anisotropic positron annihilation signals. Furthermore, we show that in experiments, the positron annihilation signals in β–Ga 2 O 3 samples seem to be often dominated by split Ga vacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Equation of state, phase transitions, and band-gap closure in Pb Cl 2 and Sn Cl 2

The equations of state and band-gap closures for Pb Cl 2 and Sn Cl 2 were studied using both experimental and theoretical methods. Here we measured the volume of both materials to a maximum pressure of 70 GPa using synchrotron-based angle-dispersive powder x-ray diffraction. The lattice parameters for both compounds showed anomalous changes between 16–32 GPa, providing evidence of a phase transition from the cotunnite structure to the related Co 2 Si structure, in contrast to the postcotunnite structure as previously suggested. First-principles calculations confirm this finding and predict a second phase transition to a Co 2 Si-like structure between 75– 110 GPa in Pb Cl 2 and 60–75 GPa in Sn Cl 2 . Band gaps were measured under compression to ~70 GPa for Pb Cl 2 and ~66 GPa for Sn Cl 2 and calculated up to 200 GPa for Pb CI 2 and 120 GPa for Sn Cl 2 . We find an excellent agreement between our experimental and theoretical results when using the Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional, which suggests that this functional could reliably be used to calculate the band gap of similar AX 2 compounds. Experimental and calculated band-gap results show discontinuous decreases in the band gap corresponding to phase changes to higher-coordinated crystal structures, giving insight into the relationship between interatomic geometry and metallicity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydrogen Dissociation Controls 1-Hexyne Selective Hydrogenation on Dilute Pd-in-Au Catalysts

Increasing the selectivity of the catalytic hydrogenation of alkynes to alkenes is of major importance for the processing of petrochemicals and the production of fine chemicals. Achieving high selectivity for alkene formation at high conversions, however, remains a long-standing challenge in heterogeneous catalysis. Here, the mechanism and origin of the high selectivity of dilute Pd-in-Au catalysts has been studied by a combination of first-principles calculations, microkinetic simulations, and isotopic exchange hydrogenation experiments using Pd 0.04 Au 0.96 nanoparticles embedded in raspberry colloid-templated silica. The Pd is predominantly in the form of isolated atoms, only surrounded by Au atoms, based on prior studies. The simulations indicate that the rate-limiting process for 1-hexyne hydrogenation on Pd monomers in Au(111) is H 2 dissociation, which has a large free energy barrier of 0.86 eV at 363 K and 0.2 bar of H 2 . The C–H bond formation steps, on the other hand, proceed with lower barriers, which contrasts with previous studies of extended Pd catalysts. The microkinetic simulations identify the sizable H 2 dissociation barrier and the small barrier for the hydrogenation of 1-hexyne as key factors that lead to a high selectivity for the production of 1-hexene from 1-hexyne, even at high conversion. The unconventional H 2 dissociation limiting process in combination with the low coverage of weakly bound hydrocarbon intermediates explains the near-zero order of 1-hexyne found experimentally. Furthermore, the partial hydrogenation of 1-hexyne to form 1-hexene is shown to be an irreversible process from our isotopic exchange hydrogenation experiments and is explained by the strongly exothermic nature of the reaction. Diluting active species, like Pd, in a less active host metal, like Au, hence appears promising as a means of tuning the binding energy of reactants and altering reaction profiles, leading to distinct kinetic behavior for an optimal catalytic activity and selectivity. Furthermore, the combination of microkinetic modeling, density functional theory calculations, and isotopic exchange experiments is thus demonstrated to be an effective approach to modeling important catalytic phenomena.

1-hexyne↗
DOE OSTI2025

HfNb_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 hafnium-niobium (Hf-Nb). 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,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetoelastic coupling and spin contributions to entropy and thermal transport in biferroic yttrium orthochromite

Direct engineering of material properties through exploitation of spin, phonon, and charge-coupled degrees of freedom is an active area of development in materials science. However, the relative contribution of the competing orders to controlling the desired behavior is challenging to decipher. In particular, the independent role of phonons, magnons, and electrons, quasiparticle coupling, and relative contributions to the phase transition free energy largely remain unexplored, especially for magnetic phase transitions. Here, we study the lattice and magnetic dynamics of biferroic yttrium orthochromite using Raman, infrared, and inelastic neutron spectroscopy techniques, supporting our experimental results with first-principles lattice dynamics and spin-wave simulations across the antiferromagnetic transition at T N ~ 138 K. Spectroscopy data and simulations together with the heat capacity (C p ) measurements, allow us to quantify individual entropic contributions from phonons (0.01 ± 0.01k B atom -1 ), dilational (0.03 ± 0.01k B atom -1 ), and magnons (0.11 ± 0.01k B atom -1 ) across TN. High-resolution phonon measurements conducted in a magnetic field show that anomalous T-dependence of phonon energies across TN originates from magnetoelastic coupling. Additionally, phonon scattering is primarily governed by the phonon–phonon coupling, with little contribution from magnon–phonon coupling, short-range spin correlations, or magnetostriction effects; a conclusion further supported by our thermal conductivity measurements conducted up to 14 T, and phenomenological modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Local positional and spin symmetry breaking as a source of magnetism and insulation in paramagnetic EuTi O 3

Here we consider theoretically the paramagnetic phases of EuTiO 3 that represent configurations created by two sets of microscopic degrees of freedom (m-DOFs): positional symmetry breaking due to octahedral rotations and magnetic symmetry breaking due to spin disorder. The effect of these sets of m-DOFs on the electronic structure and properties of the para phases is assessed by considering sufficiently large (super) cells with the required nominal global average symmetry, allowing, however, the $\mathcal{local}$ positional and magnetic symmetries to be lowered. We find that tendencies for local symmetry breaking can be monitored by following total energy lowering in mean-fieldlike density-functional theory, without recourse for strong correlation effects. While most nominally cubic ABO 3 perovskites are known for their symmetry breaking due to the B-atom sublattice, the case of f-electron magnetism in EuTiO 3 is associated with A-sublattice symmetry breaking and its coupling to structural distortions. We find that (i) paramagnetic $\mathcal{cubic}$ EuTiO 3 has an intrinsic tendency for both magnetic and positional symmetry breaking, while paramagnetic $\mathcal{tetragonal}$ EuTiO 3 has only magnetic symmetry lowering and no noticeable positional symmetry lowering with respect to low-temperature antiferromagnetic tetragonal phase. (ii) Properly modeled paramagnetic tetragonal and cubic EuTiO 3 have a nonzero local magnetic moment on each Eu ion, consistent with the experimental observations of local magnetism in the para phases of EuTiO 3 significantly above the Néel temperature. Interestingly, (iii) the local positional distortion modes in the short-range ordered para phases are inherited from the long-range ordered low-temperature antiferromagnetic ground-state phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting the Activity and Selectivity of Bimetallic Metal Catalysts for Ethanol Reforming using Machine Learning

Machine learning is ideally suited for the pattern detection in large uniform datasets, but consistent experimental datasets on catalyst studies are often small. Here we demonstrate how a combination of machine learning and first-principles calculations can be used to extract knowledge from a relatively small set of experimental data. The approach is based on combining a complex machine-learning model trained on an extensive computational library of transition-state energies with simple linear regression models of experimental catalytic activities and selectivities from the literature. Using the combined model, we identify the key C–C bond scission reactions involved in ethanol reforming and perform a computational screening for ethanol reforming on monolayer bimetallic catalysts with architectures TM-Pt-Pt(111) and Pt-TM-Pt(111) (TM = 3d transition metals). The model also predicts four promising catalyst compositions for future experimental studies. In conclusion, the approach is not limited to ethanol reforming but is of general use for the interpretation of experimental observations as well as for the computational discovery of novel catalytic materials.

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

First-principles thermodynamic assessment of Sr-containing secondary phase formation in strontium-substituted lanthanum manganites for solid oxide cell applications

Sr-secondary phase formation is a potentially significant degradation mode with direct impact upon solid-oxide cell (SOC) commercial viability. A first-principles based thermodynamic study was performed for La 1−x Sr x MnO 3±δ (LSM) perovskites to assess their stability against formation of different Sr-secondary phases, including SrO, SrCrO 4 , SrSO 4 , SrCO 3 , and Sr(OH) 2 , for SOC applications. The Sr-secondary phase formation reaction free energies were determined via a thermodynamic model by combining ab initio lattice dynamics calculations for the solid phases and ab initio thermodynamic data for the gas phases. The current approach expands the previously reported thermodynamic modeling studies by integrating first-principles based point defect equilibria into the thermodynamic analysis. The modeling results obtained using this new approach indicate an increased tendency to form the SrO oxide upon decreasing the oxygen partial pressure. Additionally, the enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with the contaminant gas species as a function of temperature and pressure.

defect thermodynamics modeling↗