Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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

Machine Learning‐Guided Discovery of High‐Entropy Perovskite Oxide Electrocatalysts via Oxygen Vacancy Engineering

Abstract High‐entropy perovskite oxides (HEPOs) have recently emerged as multifunctional catalysts. However, the HEPOs’ structural and compositional complexity hinders the easy and accurate extrapolation of activity indicators, which are essential for establishing structure‐property correlations. Here, OxiGraphX, is introduced as a novel graph neural network (GNN) model designed to capture the complex relationships among structure, composition, and atomic chemical environments for accurate prediction of oxygen vacancy formation energies (OVFEs) in HEPOs. By integrating machine learning (ML), density functional theory (DFT), and experimental validation, this work demonstrates an efficient framework for rapidly and accurately screening HEPO electrocatalysts for oxygen evolution reaction (OER). The OxiGraphX predicts OVFEs with a precision exceeding existing data, enabling the identification of compositions of higher oxygen vacancy content (OVC) and, thus, higher catalytic activity. Furthermore, the model explores latent spaces that translate effectively into experimental domains, bridging computational predictions with real‐world applications. This approach accelerates the discovery of high‐performance HEPO catalysts while providing deeper insights into their catalytic mechanisms.

Chemistry↗
DOE OSTI2026

Maximizing machine learning interatomic potential transferability for the discovery of the novel stellated octadecagon Bi18-Pt24 cage structure

Achieving true transferability remains the central challenge for Machine Learning Interatomic Potentials (ML-IAPs) in modeling complex bimetallic nanoclusters across their vast potential energy surfaces. We systematically investigate data selection strategies to optimize the Chebyshev Interaction Model for Efficient Simulation (ChIMES) potential for the Bi-Pt nanoclusters by comparing three innovative sampling methods: Principal Component Analysis (PCA)/k-means (structural diversity), t-distributedStochasticNeighborEmbedding (t-SNE)/k-means (force-space diversity), and hierarchical clustering. Quantitatively, the PCA/k-means strategy proved most effective for global accuracy, yielding the lowest force errors and achieving energy root mean square errors (RMSE) values competitive with Density Functional Theory (DFT), demonstrating excellent accuracy (19.16meV/atom). Structural validation on 34 unique DFT-optimized isomers further confirmed the potential’s high fidelity, with the best model PCA/k-means reproducing structures with an average root mean square deviation (RMSD) of 0.10 Å. However, the t-SNE methods, by maximizing diversity in the force space, demonstrated superior extrapolative power, leading to the more precise prediction of a novel stellated octadecagon Bi18⁢Pt24 cage structure, demonstrating the potential for exploring previously unseen morphologies. Our results establish a clear methodology for strategic data sampling that successfully maximizes ML-IAP transferability, providing an accurate and computationally efficient tool that accelerates the theoretical discovery of complex bimetallic architectures.

Vangheluwe, Raphaël [Université Paris-Saclay, CNRS↗
DOE OSTI2020

Extend NdJ relationship with the size, multiple exchanges and Dzyaloshinskii-Moriya interaction for Néel skyrmions in hexagonal magnetic interfaces

For the integration with other parts of spintronic devices, it is desired to minimize the size of skyrmions to a nanometer scale. To provide a theoretical guidance for this endeavour, we established an extended NdJ relationship between the skyrmion size and fundamental magnetic properties of two-dimensional hexagonal ferromagnetic lattices, including multiple exchanges between different neighbors (J 1 -J 4 ), Dzyaloshinskii-Moriya interactions (d 1 and d 2 ) and on-site magnetic anisotropy energy. We demonstrated the applicability of this approach through density functional theory and spin spiral calculations for several interface, e.g., Co/Gra, Rh/Gra, Co/Pt 3 and MgO/Co/Pt 3 , and showed that the size of Néel skyrmions can be effectively reduced by either enhancing d or more conveniently reducing J through cancelation of J i s.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modified microstructures in proton irradiated dual phase 308L weldment filler material

In this study, the effect of proton irradiation on the microstructure of δ ferrite—γ austenite mixed phase 308L filler material in a 508–304 dissimilar metal weldment was investigated over a depth of 0 to 10 µm. Ni–Si–Mn G-phase precipitates were observed with SEM and TEM in δ ferrite but not in γ austenite. Our density functional theory based calculations show that the G/Fe interface energy in δ-Fe is significantly lower than that in γ-Fe (0.35 versus 1.25 J/m 2 ), which provides a thermodynamics-based explanation for our experimental observations of preferential formation of G-phase in δ ferrite. STEM-EDS, TEM dark field imaging, and diffraction patterns confirmed the Ni–Si–Mn enriched precipitates were G-phase precipitates with a stoichiometry of Mn 6 Ni 16 Si 7 . Intragranular voids and Ni–Si enriched clusters were observed in irradiated γ austenite. Additionally, Ni and Si segregation was observed along the void interfaces. In both cases, Ni–Si clusters and segregation to voids, selected area diffraction patterns did not reveal the existence of a second phase. Proton irradiation induced Cr depletion and Si and Ni enrichment at γ-γ austenite grain boundaries that was characterized with STEM/EDS.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Level density within a micro-macroscopic approach

Statistical level density $ρ$($E, A$) is derived for nucleonic system with a given energy $E$, particle number $A$ and other integrals of motion in the micro-macroscopic approximation beyond the standard saddle-point method of the Fermi gas model. This level density reaches the two limits; the well-known Fermi gas grand-canonical ensemble limit for a large entropy $S$ related to large excitation energies, and the finite micro-canonical limit for a small combinatorical entropy $S$ at low excitation energies. In conclusion, the inverse level density parameter $K$ as function of the particle number $A$ in the semiclassical periodic orbit theory, taking into account the extended Thomas-Fermi and Strutinsky shell corrections, is calculated and compared with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Determining the Structure-Property Relationships of Quasi-Two-Dimensional Semiconductor Nanoplatelets.

We report a theoretical study of CdSe nanoplatelets aimed at identifying the main factors determining their photophysical properties. Using atomic configurations optimized with density functional theory calculations, we computed quasiparticle and exciton binding energies of nanoplatelets with two to seven monolayers. We employed many body perturbation theory at the GW level and solved the Bethe-Salpeter equation to obtain absorption spectra and excitonic properties. Our results, which agree well with recent experiments, were then used to design a model that allows us to disentangle the effects of quantum confinement, strain induced by passivating ligands, and dielectric environment on the electronic properties of nanoplatelets. We found that, for the model to accurately reproduce our first principle results, it is critical to account for surface stress and consider a finite potential barrier and energy-dependent effective masses when describing quantum confinement. Our findings call into question previous assumptions on the validity of an infinite barrier to describe carrier confinement in nanoplatelets, suggesting that it may be possible to optimize interfacial charge transfer and extraction by appropriately choosing passivating ligands. The model developed here is generalizable to core-shell platelets and enables the description of system sizes not yet directly treatable by first-principles calculations.

Greenwood, Arin R.↗
DOE OSTI2023

Impact of Disorder and Defects on Bond Energies in Highly Fluorinated Graphene

Fluorinated graphite materials are of interest for an assortment of applications and can be synthesized under a variety of synthetic conditions from many different types of carbon. Due to such variations, structural disorders in the form of defects and polymorphism are often present. Here, we investigate the impact of local structural variations on the C–F bond dissociation energies (BDEs) in carbon-based fluoride materials using density functional theory (DFT) computational methods. Employing fluorographene (FG) cluster models, we determine the impact of different C–F bonding configurations in the core of each platelet on the equilibrium BDEs for each C–F bond. The introduction of structural disorder decreases the first C–F BDE by approximately 1 eV compared to the canonical arrangement of axial C–F bonds ordered as in a network of cyclohexane “chairs”. Variability of calculated BDEs among the different polymorphs decreases upon subsequent F removal. Common structural tendencies of the adiabatic defluorination pathways for each polymorph are identified. In conclusion, our analysis suggests that at F/C ratios near 1.0, disorder in the local structure can play a significant role in the energetics of the initial carbon fluoride defluorination and that the influence of this configurational disorder diminishes with decreasing F/C ratios.

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

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Testing the predictive power of theory for Pd x Ir (100-x) alloy nanoparticles for the oxygen reduction reaction

Here, density functional theory (DFT) calculations of O and OH binding energies on triatomic surface ensembles of Pd x Ir (100-x) nanoalloys successfully predicted the overall trend in experimental oxygen reduction reaction (ORR) activity as a function of nanoparticle (NP) composition. Specifically, triatomic Pd 3 ensembles were found to possess optimal O and OH binding energies and were predicted to be highly active sites for the ORR, rivaling that of Pt(111). However, DFT calculations suggest that the O binding energy increases at active sites containing Ir, thereby decreasing ORR activity. Pd x Ir (100-x) nanoalloys were synthesized using a microwave-assisted method and their activity towards the ORR was tested using rotating disk voltammetry (RDV). As predicted, the bimetallic electrocatalysts exhibited worse catalytic activity than the Pd-only NPs. The strong qualitative correlation between the theoretical and experimental results demonstrates that the activity of individual active sites on the surface of NPs can serve as a proxy for overall activity. This is a particularly useful strategy for applying DFT calculations to electrocatalysts that are too large for true first-principle analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

[Fe 4 S 4 ] cubane in sulfite reductases: new insights into bonding properties and reactivity

The dissimilatory sulfite reductase enzyme has very characteristic active site where the substrate binds to an iron site, ligated by a siroheme macrocycle and a thiol directly connected to a [Fe 4 S 4 ] cluster. This arrangement gives the enzyme remarkable efficiency in reducing sulfite and nitrite all the way to hydrogen sulfide and ammonia. For the first time we present a theoretical study where substrate binding modalities and activation are elucidated using active site models containing proton supply side chains and the [Fe 4 S 4 ] cluster. Density functional theory (DFT) was deployed in conjunction with the energy decomposition scheme (as implemented in AMS), the quantum theory of atoms in molecules (QTAIM), and conceptual DFT (cDFT) descriptors. We quantified the role of the electrostatic interactions inside the active site created by the side chains as well as the influence of the [Fe4S4] cluster on the substrate binding. Furthermore, using conceptual DFT results we shed light of the activation process, thus, laying foundation for further mechanistic studies. We found that the bonding of the ligands to the iron complex is dominated by electrostatic interactions, but the presence of the [Fe 4 S 4 ] cubane leads to substantial changes in electronic interaction. The spin state of the cubane, however, affects the binding energy only marginally. The conceptual DFT results show that the presence of the [Fe 4 S 4 ] cubane affects the reactivity of the active site as it is involved in electron transfer. This is corroborated by an increase in the electrophilicity index, thus making the active site more prone to react with the ligands. The interaction energies between the ligand and the siroheme group are also increased upon the presence of the cubane group, thus, suggesting that the siroheme group is not an innocent spectator but plays an active role in the reactivity of the dSIR active site.

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

GW with hybrid functionals for large molecular systems

A low-cost approach for stochastically sampling static exchange during time-dependent Hartree–Fock-type propagation is presented. This enables the use of an excellent hybrid density functional theory (DFT) starting point for stochastic GW quasiparticle energy calculations. Generalized Kohn–Sham molecular orbitals and energies, rather than those of a local-DFT calculation, are used for building the Green function and effective Coulomb interaction. The use of an optimally tuned hybrid diminishes the starting point dependency in one-shot stochastic GW, effectively avoiding the need for self-consistent GW iterations.

Chemistry↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Segregation of Al and its effect on coercivity in Nd-Fe-B

Doping a small amount of Al can effectively enhance coercivity in Nd-Fe-B magnets. Here we investigated the partitioning behavior of Al and its effect on coercivity in Nd-Fe-B using first principles DFT (density functional theory) calculation and micromagnetic simulation. The calculated substitution energies of Fe by Al are negative at the crystallographic sites of 4c and 8j 2 while they are positive values at the other sites in Nd 2 Fe 14 B (2:14:1), implying a small solubility of Al in 2:14:1. Further, Al prefers to segregate at grain boundary (GB) and stabilize the Nd-Fe-Al phase with a Nd 6 Fe 13 Si-type tetragonal structure (6:13:1). The formation of the antiferromagnetic or weak ferrimagnetic 6:13:1-like phase depletes Fe and reduces the amount of ferromagnetic Nd–Fe type grain boundary phase (GBP), which weakens the inter-grain magnetic interaction. Micromagnetic simulations indicate that the 6:13:1-like GBP increases the pinning field of magnetic domain wall at GB and suppresses the nucleation of reversal magnetic domain on the grain surface of 2:14:1 during demagnetization process. The formation of Al-rich shell on 2:14:1 grain surface can further moderately increase the domain pinning field at GB and the nucleation field of 2:14:1 grain. Developing novel processing method to tailor Al segregation and promote formation of 6:13:1-like phase at GB can be a promising approach to improve coercivity in Nd-Fe-B magnet.

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

Many-fermion simulation from the contracted quantum eigensolver without fermionic encoding of the wave function

Quantum computers potentially have an exponential advantage over classical computers for the quantum simulation of many-fermion quantum systems. Nonetheless, fermions are more expensive to simulate than bosons due to the fermionic encoding—a mapping by which the qubits are encoded with fermion statistics. Here we generalize the contracted quantum eigensolver (CQE) to avoid fermionic encoding of the wave function. In contrast to the variational quantum eigensolver, the CQE solves for a many-fermion stationary state by minimizing the contraction (projection) of the Schrödinger equation onto two fermions. We avoid fermionic encoding of the wave function by contracting the Schrödinger equation onto an unencoded pair of particles. Solution of the resulting contracted equation by a series of unencoded two-body exponential transformations generates an unencoded wave function from which the energy and two-fermion reduced density matrix (2-RDM) can be computed. We apply the unencoded and the encoded CQE algorithms to the hydrogen fluoride molecule, the dissociation of oxygen O 2 , and a series of hydrogen chains. Both algorithms show comparable convergence towards the exact ground-state energies and 2-RDMs, but the unencoded algorithm has computational advantages in terms of state preparation and tomography.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Variation of carrier density in semimetals via short-range correlation: A case study with nickelate NdNiO 2

Carrier density is one of the key controlling factors of material properties, particularly in controlling the essential correlations in strongly correlated materials. Typically, carrier density is externally tuned by doping or gating and remains fixed below room temperature. Strangely, the carrier density in correlated semimetals is often found to vary sensitively against weak external controls such as temperature, magnetic field, and pressure. Here, we develop a realistic simulation scheme that incorporates interatomic noncollinear magnetic correlation without a long-range order. Using the recently discovered nickelate superconductor as an example, we demonstrate a rather generic low-energy mechanism that in semimetals short-range correlation can reversely modulate the carrier density as well. Such a mutual influence between correlation and carrier density provides an extra ingredient for sensitive bifurcating behavior. This special feature of correlated semimetals explains their versatile carrier density at low energy and opens up new possibilities of functionalizing these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electronic structure and small-hole polarons in YTiO 3

As a prototypical Mott insulator with ferromagnetic ordering, YTiO 3 (YTO) is of great interest in the study of strong electron correlation effects and orbital ordering. Here we report the first molecular beam epitaxy (MBE) growth of YTO films, combined with theoretical and experimental characterizations of the electronic structure and charge transport properties. The obstacles of YTO MBE growth are discussed and potential routes to overcome them are proposed. DC transport and Seebeck measurements on thin films and bulk single crystals identify p-type Arrhenius transport behavior with an activation energy of ~0.17 eV in thin films, consistent with the energy barrier for small hole polaron migration from hybrid density functional theory calculations. Hard x-ray photoelectron spectroscopy measurements show the lower Hubbard band at 1.1 eV below the Fermi level, whereas a Mott-Hubbard band gap of ~1.5 eV is determined from photoluminescence measurements. These findings provide critical insight into the electronic band structure of YTO and related materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evidence of antiferromagnetism in ultrathin metallic (111)-oriented LaNiO 3 films

Antiferromagnets with exotic spin textures promise low-power spintronic devices with extremely high operating frequencies and resistance to external perturbations. In particular, the combination of highly tunable correlated electron physics, as in complex oxides, with metallicity and antiferromagnetism is desirable but exceedingly rare. LaNiO 3 , the lone example of a perovskite nickelate which is metallic across all temperatures, has long been a promising candidate, but the antiferromagnetic metallic state has remained elusive. Here, we demonstrate the emergence of this state in ultrathin films of (111)-oriented LaNiO 3 using a combination of polarized neutron reflectometry, low-energy muon spectroscopy and anomalous Hall effect measurements, and density functional theory calculations with a Hubbard U term. Further, we find a highly strained symmetry-breaking interfacial region which may support canting of the AFM moments leading to a vanishingly small net magnetization at the film-substrate interface, providing a convenient route toward control of the Néel order. Evidence of antiferromagnetic metallic behavior in (111)LNO films highlights the role of crystal symmetry in tuning the novel quantum states in complex oxides.

36 MATERIALS SCIENCE↗