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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Group-Additivity–Embedded Multiscale Modeling for Electric Field-Enhanced Nanocatalysis

Elucidating structure-performance relationships remains a central challenge in field-enhanced catalysis, where nanoparticles exhibit nonuniform surface sites with site-dependent responses to electric fields. Low-coordination sites (edges, corners, and tips) are particularly electric field-sensitive (EF), leading to nonuniform charge distribution, adsorption energies, and catalytic activity. Here, using ammonia decomposition on a ruthenium cluster as a model system, we develop a transferable multiscale framework integrating density functional theory, group additivity (GA), Brønsted-Evans-Polanyi scaling, and microkinetic modeling to predict EF-dependent activity across nonuniform cluster sites. Across sites and fields, the nitrogen adsorption energy (E N ) emerges as the governing descriptor, yielding robust volcano relationships whose optimum shifts systematically with field: negative fields strengthen N binding via electron accumulation, while positive fields weaken N binding via charge depletion, moving the optimal E N toward weaker binding. Microkinetic analysis shows that N≡N bond formation remains the key kinetic bottleneck over most conditions; positive fields lower the effective barrier and, critically, increase the fraction of near-optimal active sites, leading to a net enhancement in overall activity relative to zero-field and negative-field cases. By capturing EF- and site-dependent energetics with high accuracy and low computational cost, this GA-embedded multi-scale simulation workflow provides a physically interpretable route to predict and design field-enhanced nanocatalysis.

ammonia decomposition↗
DOE OSTI2023

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Ternary PtIrNi Catalysts for Efficient Electrochemical Ammonia Oxidation

Due to high energy density, carbon-free feature, and easiness to be liquefied, ammonia (NH 3 ) has proved to be an effective alternative to hydrogen in low temperature fuel cells via its direct the ammonia oxidation reaction (AOR) for electricity generation. However, the kinetically sluggish AOR has prohibitively hindered the attractive direct ammonia fuel cell (DAFC) applications. Here we report an efficient AOR catalyst, in which ternary PtIrNi alloy nanoparticles well dispersed on a binary composite support consisting of porous silicon dioxide (SiO 2 ) and carboxyl-functionalized carbon nanotube (PtIrNi/SiO 2 -CNT-COOH) through a sonochemical-assisted synthesis strategy. The PtIrNi alloy nanoparticles, with the aid of abundant OHad provided by porous SiO 2 , and the improved electrical conductivity by CNTs, exhibit remarkable catalytic activity for the AOR in alkaline media. It is evidenced by a lower onset potential (~0.40 V vs. RHE) at room temperature, than that of commercial PtIr/C (ca. 0.43 V vs. RHE). Increasing NH 3 concentrations and operation temperatures significantly improve AOR performance of this catalyst. Specifically, AOR activity of the optimal PtIrNi nanoparticle catalyst can be significantly enhanced by elevating the temperature to 80ºC, with a much lower onset potential (~0.32 V vs. RHE), indicating that DAFC can be operated at higher temperature for increased performance. Constant-potential density functional theory (DFT) calculations showed that the Pt-Ir ensembles on {100}-terminated surfaces serve as the active site. Importantly, the introduction of Ni raises the center energy of the density of states projected onto the group d-orbitals of surface sites and thus lowers the theoretical onset potential for *NH 2 dehydrogenation to *NH when compared to Pt and Pt 3 Ir alloy.

25 ENERGY STORAGE↗
DOE OSTI2025

Xenon–metal pair formation in UO 2 investigated using DFT + U

A recent experimental study on a spent uranium dioxide (UO 2 ) fuel sample from Belgium Reactor 3 identified a unique pair structure formed by the noble metal phase (NMP) and fission gas [xenon (Xe)] precipitate. However, the fundamental mechanism behind this structure remains unclear. The present study aims to provide an understanding of the interaction between five different metal precipitates [molybdenum (Mo), ruthenium (Ru), palladium (Pd), technetium (Tc), and rhodium (Rh)] and the Xe fission gas atoms in UO 2 , by using density functional theory (DFT) in combination with the Hubbard U correction to compute the formation energies involved. All DFT + U calculations were performed with occupation matrix control to ensure antiferromagnetic ordering of UO 2 . The calculated formation and binding energies of the Xe and solid fission products in the NMP reveal that these metal precipitates form stable pair structures with Xe. Notably, the formation energy of Xe–metal pairs is lower than that of the isolated single defects in all instances, with Pd and Mo showing the most favorable binding energy, likely accounting for the observed pair structure formation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Computationally Predicted High-Throughput Free-Energy Phase Diagrams for the Discovery of Solid-State Hydrogen Storage Reactions

The design of multinary solid-state material systems that undergo reversible phase changes via changes in temperature and pressure provides a potential means of safely storing hydrogen. However, fully mapping the stabilities of known or newly targeted compounds relative to competing phases at reaction conditions has previously required many stringent experiments or computationally demanding calculations of each compound’s change in Gibbs energy with respect to temperature, G(T). Here, we have extended the approach of constructing chemical potential phase diagrams based on ΔG f (T) to enable the analysis of phase stability at non-zero temperatures. We first performed density functional theory calculations to compute the formation enthalpies of binary, ternary, and quaternary compounds within several compositional spaces of current interest for solid-state hydrogen storage. Temperature effects on solid compound stability were then accounted for using our recently introduced machine learned descriptor for the temperature-dependent contribution G δ (T) to the Gibbs energy G(T). From these Gibbs energies, we evaluated each compound’s stability relative to competing compounds over a wide range of conditions and show using chemical potential and composition phase diagrams that the predicted stable phases and H2 release reactions are consistent with experimental observations. This demonstrates that our approach rapidly computes the thermochemistry of hydrogen release reactions for compounds at sufficiently high accuracy relative to experiment to provide a powerful framework for analyzing hydrogen storage materials. This framework based on G(T) enables the accelerated discovery of active materials for a variety of technologies that rely on solid-state reactions involving these materials.

08 HYDROGEN↗
DOE OSTI2021

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Atomic-scale vacancy engineering unlocks basal-plane catalytic activity in metallic WSe 2 for reversible oxygen electrocatalysis

Two-dimensional metallic transition metal dichalcogenides offer high electrical conductivity and large surface areas for electrocatalysis, yet their inherent basal planes are catalytically inert. Here, we present an atomic-scale vacancy engineering strategy to activate the basal surfaces of metallic WSe 2 for reversible oxygen electrocatalysis. This approach, based on intentionally designed substitutional metal doping, promotes the spontaneous formation of selenium vacancies while preserving the metallic 1 T′ phase, thereby creating highly reactive and oxygen-affinitive sites. Density functional theory calculations reveal that these vacancy-mediated metal complexes dramatically lower the energy barriers for initial oxygen adsorption, enabling dissociative oxygen adsorption. Operando and ex-situ spectroscopic analyses confirm that vacancy-mediated metal complexes transform into dynamic Se/W-oxide intermediates under operating conditions. Se/W-oxides on the surface experimentally and theoretically prove electrocatalytic activity and reversibility. Applying this strategy in lithium–oxygen batteries, the basal-plane activated WSe 2 shows high discharge capacities (9868 mA h g −1 , corresponding to 3947 mA h g$^{-1}_{cathode}$), impressive cycle retention over 550 cycles at 1000 mA h g −1 , and outstanding rate–capability over a wide current–density range (100–3000 mA g −1 ) during 256 cycles.

2D materials↗
DOE OSTI2023

Material Effects on Electron-Capture Decay in Cryogenic Sensors

Several current searches for physics beyond the standard model are based on measuring the electron-capture (EC) decay of radionuclides implanted into cryogenic high-resolution sensors. The sensitivity of these experiments has already reached the level where systematic effects related to atomic state energy changes from the host material are a limiting factor. One example is a neutrino mass study based on the nuclear EC decay of 7 Be to 7 Li inside cryogenic Ta-based sensors. Here, to understand the material effects at the required level, we use density-functional theory to model the electronic structure of lithium atoms in different atomic environments of the polycrystalline Ta absorber film. The calculations reveal that the Li 1s binding energies can vary by more than 2 eV due to insertion at different lattice sites, at grain boundaries, in disordered Ta, and in the vicinity of various impurities. However, the total range of Li 1s shifts does not exceed 4 eV, even for extreme amorphous disorder. Furthermore, when investigating the effects on the Li 2s levels, we find broadening of more than 5 eV due to hybridization with the Ta band structure. Material effects are shown to contribute significantly to peak broadening in Ta-based sensors that are used to search for physics beyond the standard model in the EC decay of 7 Be, but they do not explain the full extent of observed broadening. Understanding these in-medium effects will be required for current- and future-generation experiments that observe low-energy radiation from the EC decay of implanted isotopes to evaluate potential limitations on the measurement sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The effect of elastic strains on the adsorption energy of H, O, and OH in transition metals

The influence of elastic strains on the adsorption of H, O, and OH on the (111) surfaces of 8 fcc (Ni, Cu, Pd, Ag, Pt, Au, Rh, Ir) and on the (0001) surfaces of 3 hcp (Co, Zn, Cd) transition metals was analyzed by means of density functional theory calculations. To this end, surface slabs were subjected to different strain states (uniaxial, biaxial, shear, and a combination of them) up to strains dictated by the mechanical stability limits indicated by phonon calculations. It was found that the adsorption energy followed the predictions of the d-band theory but – surprisingly – the variations in the adsorption energy only depended on the area of the adsorption hole and not on the particular elastic strain tensor applied to achieve this area. The analysis of the electronic structure showed that the applied strains did not modify the shape of the Projected Density of States (PDOS) of the d-orbitals of the transition metals but only led to a shift in the energy levels. Moreover, the presence of the adsorbates on the surfaces led to negligible changes in the PDOS. Thus, the adsorption energies were a function of the Fermi energy which in turn was associated with the change of the area of the adsorption through a general linear law that was valid for all metals. Furthermore, the information in this paper allows the immediate and accurate estimation of the effect of any elastic strain on the adsorption energies of H, O, and OH in 11 transition metals with more than half-filled d-orbitals.

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

Aqueous Diels–Alder reactions for thermochemical storage and heat transfer fluids identified using density functional theory

Thermal storage and transfer fluids have important applications in industrial, transportation, and domestic settings. Current thermal fluids have relatively low specific heats, often significantly below that of water. However, by introducing a thermochemical reaction to a base fluid, it is possible to enhance the fluid's thermal properties. Here, density functional theory (DFT) is used to screen Diels–Alder reactions for use in aqueous thermal fluids. From an initial set of 52 reactions, four are identified with moderate aqueous solubility and predicted turning temperature near the liquid region of water. These reactions are selectively modified through 60 total functional group substitutions to produce novel reactions with improved solubility and thermal properties. Among the reactions generated by functional group substitution, seven have promising predicted thermal properties, significantly improving specific heat (by as much as 30.5%) and energy storage density (by as much as 4.9%) compared to pure water.

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

High-Energy Hybridized States Enable Long-Lived Hot Electrons in Cobaloxime-Silicon Nanocrystal System

Strong electronic coupling is achieved between the molecular catalyst cobaloxime ([Co]) and silicon nanocrystals (Si NCs) bridged by an ethylenepyridine group derived from vinylpyridine (vpy) covalently bound to the Si NC surface (Si-vpy-[Co]). The ethylenepyridine tether in Si-vpy-[Co] is key to dramatic changes to the system’s physical properties which are not observed in the corresponding formylpyridine (fpy) system (Si-fpy-[Co]) consistent with strong electronic coupling previously observed only in dark electrochemical systems. UV−vis absorption spectroscopy reveals new [Co]-centered electronic states in Si-vpy-[Co], and transient absorption spectroscopy finds a strong absorption feature appearing within 250 fs and persisting for at least 5 ns. Astoundingly, spectroelectrochemical measurements reveal that this absorption feature is consistent with both the singly reduced [Co] − and doubly reduced [Co] 2− complexes, leading to the conclusion that these long-lived charges are derived from high-energy “hot” electrons residing in [Co]-centered states. Detailed analysis using cyclic voltammetry, spectroelectrochemistry, electron paramagnetic resonance spectroscopy, and density functional theory (DFT) calculations provides insight into the unique electronic structure created in Si-vpy-[Co]. DFT reveals that the new electronic states arise from hybridization between deep Si NC band states and high-energy molecular orbitals of the ethylenepyridine tether and the [Co] catalyst and are facilitated by σ-bonding character at the ethylenepyridine linkage. This study demonstrates that strong electronic coupling achieved through precise molecular chemistry can change the paradigm of otherwise fixed energy levels in hybrid photoelectrochemical systems for artificial photosynthesis and related applications.

14 SOLAR ENERGY↗
DOE OSTI2024

The Interplay of Solvation and Polarization Effects on Ion Pairing in Nanoconfined Electrolytes

The nature of ion–ion interactions in electrolytes confined to nanoscale pores has important implications for energy storage and separation technologies. However, the physical effects dictating the structure of nanoconfined electrolytes remain debated. Here we employ machine-learning-based molecular dynamics simulations to investigate ion–ion interactions with density functional theory level accuracy in a prototypical confined electrolyte, aqueous NaCl within graphene slit pores. We find that the free energy of ion pairing in highly confined electrolytes deviates substantially from that in bulk solutions, observing a decrease in contact ion pairing but an increase in solvent-separated ion pairing. These changes arise from an interplay of ion solvation effects and graphene’s electronic structure. Notably, the behavior observed from our first-principles-level simulations is not reproduced even qualitatively with the classical force fields conventionally used to model these systems. The insight provided in this work opens new avenues for predicting and controlling the structure of nanoconfined electrolytes.

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

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗
DOE OSTI2026

Direct cross section measurement of 102 Pd ⁢(𝛾,𝑝) and 102 Pd ⁢(𝛾,𝛼) for the astrophysical 𝑝 process

Background: A handful of neutron-deficient stable nuclei, known as the “p nuclei,” cannot be produced through astrophysical neutron capture processes. Instead, some of these nuclei are proposed to be produced by 𝛾-induced reactions on existing r- and s-process seeds. The specific astrophysical site or sites are not yet identified, however, with uncertainties in the cross sections of these 𝛾-induced reactions playing a role. Databases of reaction rates for astrophysical simulations often rely on theoretical statistical model calculations, such as Hauser-Feshbach, for rates where no experimental information is known. However, reasonable variations in the choice of parametrizations of various nuclear properties can create order-of-magnitude variations in the final predicted cross sections and reaction rates, which are then propagated through the models to the predicted final abundances. Purpose: To better constrain these statistical model calculations and ultimately reduce the uncertainties from the nuclear physics on our understanding of the p nuclei, a measurement of the cross sections of 𝛾-induced reactions on the p-nucleus 102 Pd was undertaken. This work represents the first measurement of its kind, using segmented silicon detectors to measure prompt charged particle emission from 𝛾-induced reactions. Methods: Quasimonoenergetic gamma beams from the High Intensity 𝛾 Source facility bombarded an enriched 102 Pd target. A segmented silicon array was arranged to detect the particles emitted from (𝛾,𝑝) and (𝛾,𝛼) reactions. Results: Reaction cross sections were deduced at multiple 𝛾-beam energies between 10 and 19 MeV, and compared to statistical model calculations using talys-1.96. The 102 Pd ⁢(𝛾,𝑝)⁢ 101 Rh reaction cross section was reasonably well reproduced by a subset of photon strength functions and level densities, though the strength to the ground state of 101 Rh was underestimated at higher incident gamma energies. The 102 Pd ⁢(𝛾,𝛼)⁢ 98 Ru was in general overpredicted by the various alpha-nucleus optical model potentials. Conclusions: While the theoretical cross sections used to model the (𝛾,𝑝) reactions for the p process may be reasonable, a more careful approach is needed in the case of (𝛾,𝛼). Further work to probe gamma-induced reaction cross sections at and near the p nuclei is warranted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Factors Governing Oxygen Vacancy Formation in Oxide Perovskites

The control of oxygen vacancy (V O ) formation is critical to advancing multiple metal-oxide-perovskite-based technologies. In this work, we report the construction of a compact linear model for the neutral V O formation energy in ABO 3 perovskites that reproduces, with reasonable fidelity, Hubbard-U-corrected density functional theory calculations based on the state-of-the-art, strongly constrained and appropriately normed exchange-correlation functional. We obtain a mean absolute error of 0.45 eV for perovskites stable at 298 K, an accuracy that holds across a large, electronically diverse set of ABO 3 perovskites. Our model considers perovskites containing alkaline-earth metals (Ca, Sr, and Ba) and lanthanides (La and Ce) on the A-site and 3d transition metals (Ti, V, Cr, Mn, Fe, Co, and Ni) on the B-site in six different crystal systems (cubic, tetragonal, orthorhombic, hexagonal, rhombohedral, and monoclinic) common to perovskites. Physically intuitive metrics easily extracted from existing experimental thermochemical data or via inexpensive quantum mechanical calculations, including crystal bond dissociation energies and (solid phase) reduction potentials, are key components of the model. Beyond validation of the model against known experimental trends in materials used in solid oxide fuel cells, the model yields new candidate perovskites not contained in our training data set, such as (Bi,Y)(Fe,Co)O 3 , which we predict may have favorable thermochemical water-splitting properties. The confluence of sufficient accuracy, efficiency, and interpretability afforded by our model not only facilitates high-throughput computational screening for any application that requires the precise control of V O concentrations but also provides a clear picture of the dominant physics governing V O formation in metal-oxide perovskites.

08 HYDROGEN↗
NASA NTRS1983

Ultrasonic wave propagation in two-phase media: Spherical inclusions

The scattering theory, recently developed via the extended method of equivalent inclusion, is used to study the propagation of time-harmonic waves in two-phase media of elastic matrix with randomly distributed elastic spherical inclusion materials. The elastic moduli and mass density of the composite medium are determined as functions of frequencies when given properties and concentration of the spheres and the matrix. Velocity and attenuation of ultrasonic waves in two-phase media are determined for cases of distributed spheres and localized damage. An averaging theorem that requires the equivalence of the strain energy and the kinetic energy between the effective medium and the original matrix with spherical inhomogeneities is employed to derive the effective moduli and mass density. The functional dependency of these quantities upon frequencies and concentration provides a method of data analysis in ultrasonic evaluation of material properties. Numerical results or moduli, velocity and/or attenuation as functions of concentration of inclusion material, or porosity, are graphically displayed.

Fu, L. S.↗
NASA NTRS1984

Ultrasonic wave propagation in two-phase media - Spherical inclusions

The scattering theory, recently developed via the extended method of equivalent inclusion, is used to study the propagation of time-harmonic waves in two-phase media of elastic matrix with randomly distributed elastic spherical inclusion materials. The elastic moduli and mass density of the composite medium are determined as functions of frequencies when given properties and concentration of the spheres and the matrix. Velocity and attenuation of ultrasonic waves in two-phase media are determined for cases of distributed spheres and localized damage. An averaging theorem that requires the equivalence of the strain energy and the kinetic energy between the effective medium and the original matrix with spherical inhomogeneities is employed to derive the effective moduli and mass density. The functional dependency of these quantities upon frequencies and concentration provides a method of data analysis in ultrasonic evaluation of material properties. Numerical results or moduli, velocity and/or attenuation as functions of concentration of inclusion material, or porosity, are graphically displayed.

Fu, L. S.↗