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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Power density in direct nuclear-pumped He-3 lasers

The interaction of neutron beams with He-3 gas is of interest for nuclear pumped lasers. The effects of spectral dependence of the neutron beam, neutron attenuation in the gas-filled laser tube, and transport of the charged-particle He-3(n, p)He-3 reaction products are treated in detail. An expression for the energy density as a function of position within the tube, tube radius, operating pressure, and neutron fluence is given. The maximum energy density within the optical cavity is achieved when the tube radius is 3.26/P where P (atm) is the operating pressure. The variation of radius by 50% above and below optimum will change the energy density at most by 10%, although performance degrades quickly for radii outside this range. If the optimum tube radius is used for each operating pressure, then the power density on the centerline (kW/cu cm) is given as 9.3 x 10 to the -18th power times the operating pressure times flux in a thermal neutron environment of flux (n/sq cm sec).

Wilson, J. W.↗
NASA NTRS1987

Simplified curve fits for the transport properties of equilibrium air

New, improved curve fits for the transport properties of equilibruim air have been developed. The curve fits are for viscosity and Prandtl number as functions of temperature and density, and viscosity and thermal conductivity as functions of internal energy and density. The curve fits were constructed using grabau-type transition functions to model the tranport properties of Peng and Pindroh. The resulting curve fits are sufficiently accurate and self-contained so that they can be readily incorporated into new or existing computational fluid dynamics codes. The range of validity of the new curve fits are temperatures up to 15,000 K densities from 10 to the -5 to 10 amagats (rho/rho sub o).

Srinivasan, S.↗
DOE OSTI2023

System-size dependence of the charged-particle pseudorapidity density at $\sqrt{s_{NN}}$ = 5.02 TeV for pp, p–Pb, and Pb–Pb collisions

We present the first systematic comparison of the charged-particle pseudorapidity densities for three widely different collision systems, pp, p–Pb, and Pb–Pb, at the top energy of the Large Hadron Collider ($\sqrt{s_{NN}}$ = 5.02 TeV) measured over a wide pseudorapidity range (-3.5 < η < 5), the widest possible among the four experiments at that facility. The systematic uncertainties are minimised since the measurements are recorded by the same experimental apparatus (ALICE). The distributions for p–Pb and Pb–Pb collisions are determined as a function of the centrality of the collisions, while results from pp collisions are reported for inelastic events with at least one charged particle at midrapidity. The charged-particle pseudorapidity densities are, under simple and robust assumptions, transformed to charged-particle rapidity densities. This allows for the calculation and the presentation of the evolution of the width of the rapidity distributions and of a lower bound on the Bjorken energy density, as a function of the number of participants in all three collision systems. We find a decreasing width of the particle production, and roughly a smooth ten fold increase in the energy density, as the system size grows, which is consistent with a gradually higher dense phase of matter

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Auxiliary function approach for determining symmetry energy at suprasaturation densities

Nuclear symmetry energy E sym (ρ) at density ρ is normally expanded or simply parameterized as a function of χ = (ρ – ρ 0 )/3ρ 0 in the form of Esym(ρ) ≈ S + Lχ + 2 –1 K sym χ 2 + 6 –1 J sym χ 3 + … using its magnitude S, slope L, curvature K sym and skewness J sym at the saturation density ρ0 of nuclear matter. Much progress has been made in recent years in constraining especially the S and L parameters using various terrestrial experiments and astrophysical observations. Furthermore, such kind of expansions/parameterizations do not converge at supra-saturation densities where χ is not small enough, hindering an accurate determination of high-density E sym (ρ) even if its characteristic parameters at ρ0 are all well determined by experiments/observations. By expanding the E sym (ρ) in terms of a properly chosen auxiliary function Π sym (χ,Θ sym ) with a parameter Θsym fixed accurately by an experimental E sym (ρ r ) value at a reference density ρr, we show that the shortcomings of the χ-expansion can be completely removed or significantly reduced in determining the high-density behavior of E sym (ρ). In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher χ-order contributions and converges to the same E sym (ρ) much faster than the conventional χ-expansion at densities ≲3ρ 0 . Several quantitative demonstrations using Monte Carlo simulations are given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Highly accurate and constrained density functional obtained with differentiable programming

Using an end-to-end differentiable implementation of the Kohn-Sham self-consistent field equations, we obtain a highly accurate neural network–based exchange and correlation (XC) functional of the electronic density. The functional is optimized using information on both energy and density while exact constraints are enforced through an appropriate neural network architecture. Here we evaluate our model against different families of XC approximations and show that at the meta-GGA level our functional exhibits unprecedented accuracy for both energy and density predictions. For nonempirical functionals, there is a strong linear correlation between energy and density errors. We use this correlation to define an XC functional quality metric that includes both energy and density errors, leading to an improved way to rank different approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Density functional study of atomic arrangements in CrMnFeCoNi high-entropy alloy and their impact on vacancy formation energy and segregation

Here, using the density functional theory-coupled Monte Carlo approach, we explored the chemical short-range order (SRO) and element segregation in equimolar CrMnFeCoNi alloy. We found that state-of-the-art approximation of random element distribution is only applicable at > 1100 K close to the melting temperature, while the Cr-Cr repulsion driving the system stabilization and accompanying the formation of cubic Cr sublattice, and mild Ni-Ni attraction are the most prominent pair interactions at lower temperatures. Chemical potential and vacancy formation energy calculations indicate that Cr is most sensitive to the local chemical environment, making Cr atoms most stabilized when the preferred SRO is introduced. While the vacancy formation is predicted equally probable among five constituting elements in the random solid solution, Cr and Ni atoms show the lowest vacancy formation energies in the structure with SRO. Furthermore, distinct element segregation was predicted in the vicinity of planar defects, including symmetric tilt grain boundary and stacking fault, which we correlated to the site- and chemistry-dependent atomic volume and bond lengths. It suggests that the local mechanical strain and bond energy induce the SRO development and element segregation: Namely, the system takes advantage of segregation of Ni atoms having large atomic volume or Cr-Cr pairs having elongated bond lengths to fill in the excess volume at defects that relaxes the mechanical strain field and optimizes bond energy distribution. The correlation between the SRO and properties of CrMnFeCoNi alloy needs further investigations, which is expected to greatly help understand and control the properties of high-entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
DOE OSTI2022

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

New Engineering Concepts to High Energy Density Li-S Batteries

Li-S batteries (LSBs) with energy densities > 500 Wh/kg @ C/5 and >1000 cycles, cycle life operating over wide temperatures has been identified to be critical for meeting the global energy demands. Progress is however, very much limited due to several major hurdles related to the following: (a) poor Li' reaction kinetics (conductivity/diffusivity) enabling full reversible S conversion to Li2S, (b) inferior electron (e) conductivity, specific capacity and cycle life due to electrolyte soluble species of polysulfides (PSs) formed with the reaction of Li ions with sulfur migrating to anode, (c) unstable solid electrolyte interphase (SEI) formation and finally, (d) dendrites that are formed on the Li metal anode severely limiting the cycle life due to polarization at the solid/liquid interfaces causing safety, inferior capacity, energy density, rate, and cycle life issues. Creation of a fully functional economical ≤$80/kWh, high energy density LSBs is vital for enabling its widespread commercial deployment and use in current and next generation electric vehicles.

25 ENERGY STORAGE↗
DOE OSTI2021

A framework for quantifying uncertainty in DFT energy corrections

In this work, we demonstrate a method to quantify uncertainty in corrections to density functional theory (DFT) energies based on empirical results. Such corrections are commonly used to improve the accuracy of computational enthalpies of formation, phase stability predictions, and other energy-derived properties, for example. We incorporate this method into a new DFT energy correction scheme comprising a mixture of oxidation-state and composition-dependent corrections and show that many chemical systems contain unstable polymorphs that may actually be predicted stable when uncertainty is taken into account. We then illustrate how these uncertainties can be used to estimate the probability that a compound is stable on a compositional phase diagram, thus enabling better-informed assessments of compound stability.

42 ENGINEERING↗
DOE OSTI2020

A step in the direction of resolving the paradox of Perdew–Zunger self-interaction correction. II. Gauge consistency of the energy density at three levels of approximation

The Perdew-Zunger(PZ) self-interaction correction (SIC) was designed to correct the one-electron limit of any approximate density functional for the exchange-correlation (xc) energy, while yielding no correction to the exact functional. Unfortunately, it spoils the slowly-varying-in-space limits of the uncorrected approximate functionals, where those functionals are right by construction. The right limits can be restored by locally scaling down the energy density of the PZ SIC in many-electron regions, but then a spurious correction to the exact functional would be found unless the self-Hartree and exact self-xc terms of the PZ SIC energy density were expressed in the same gauge. Only the local density approximation satisfies the same-gauge condition for the energy density, which explains why the recent local-scaling SIC (LSIC) is found here to work excellently for atoms and molecules only with this basic approximation, and not with the more advanced generalized gradient approximations (GGAs) and meta-GGAs, which lose the Hartree gauge via simplifying integrations by parts. The transformation of energy density that achieves the Hartree gauge for the exact xc functional can also be applied to approximate functionals. Doing so leads to a simple scaled-down self-interaction (sdSIC) correction that is typically much more accurate than PZ SIC in tests for many molecular properties (including equilibrium bond lengths). The present work shows unambiguously that the largest errors of PZ SIC applied to standard functionals at three levels of approximation can be removed by restoring their correct slowlyvarying- density limits. Here, it also confirms the relevance of these limits to atoms and molecules.

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

Invariant surface elastic properties in FCC metals and their correlation to bulk properties revealed by machine learning methods

In this work, we present a combination of machine-learned models that predicts the surface elastic properties of general free surfaces in face-centered cubic (FCC) metals. These models are built by combining a semi-analytical method based on atomistic simulations to calculate surface properties with the artificial neural network (ANN) method or the boosted regression tree (BRT) method. The latter is also used to link bulk properties and surface orientation to surface properties. The surface elastic properties are represented by their invariants considering plane elasticity within a polar method. The resulting models are shown to accurately predict the surface elastic properties of seven pure FCC metals (Cu, Ni, Ag, Au, Al, Pd, Pt). The BRT model reveals the correlations between bulk and corresponding surface properties in terms of invariants, which can be used to guide the design of complex nano-sized particles, wires and films. Finally, by expressing the surface excess energy density as a function of surface elastic invariants, fast predictions of surface energy as a function of in-plane deformations can be made from these model constructs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simplified tuning of long-range corrected density functionals for use in symmetry-adapted perturbation theory

Long considered a failure, second-order symmetry-adapted perturbation theory (SAPT) based on Kohn–Sham orbitals, or SAPT0(KS), can be resurrected for semiquantitative purposes using long-range corrected density functionals whose asymptotic behavior is adjusted separately for each monomer. Here, as in other contexts, correct asymptotic behavior can be enforced via “optimal tuning” based on the ionization energy theorem of density functional theory, but the tuning procedure is tedious, expensive for large systems, and comes with a troubling dependence on system size. Here, we show that essentially identical results are obtained using a fast, convenient, and automated tuning procedure based on the size of the exchange hole. In conjunction with “extended” (X)SAPT methods that improve the description of dispersion, this procedure achieves benchmark-quality interaction energies, along with the usual SAPT energy decomposition, without the hassle of system-specific tuning.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

The Fermi–Löwdin self-interaction correction for ionization energies of organic molecules

(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. Here, the inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction correction on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator zs , only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.

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

Machine-learning-aided density functional theory calculations of stacking fault energies in steel

A combined large-scale first principles approach with machine learning and materials informatics is proposed to quickly sweep the chemistry-composition space of advanced high strength steels (AHSS). AHSS are composed of iron and key alloying elements such as aluminum and manganese. A systematic exploration of the distribution of aluminum and manganese atoms in iron is used to investigate low stacking fault energies configurations using first principles calculations. To overcome the computational cost of exploring the composition space, this process is sped up using an automated machine learning tool: DeepHyper. Here our results predict that it is energetically favorable for Al to stay away from a stacking fault, but Mn atoms do not affect the stacking fault energy and can stay in the vicinity of the fault. The distribution of Al and Mn atoms in systems containing stacking faults and the effects of their interactions on the equilibrium distribution are systematically analyzed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluating Material Design Principles for Calcium-Ion Mobility in Intercalation Cathodes

Multivalent-ion batteries offer an alternative to Li-based technologies, with the potential for greater sustainability, improved safety, and higher energy density, primarily due to their rechargeable system featuring a passivating metal anode. Although a system based on the Ca 2+ /Ca couple is particularly attractive given the low electrochemical plating potential of Ca 2+ , the remaining challenge for a viable rechargeable Ca battery is to identify Ca cathodes with fast ion transport. In this work, a high-throughput computational pipeline is adapted to (1) discover novel Ca cathodes in a largely unexplored space of empty intercalation hosts and (2) develop material design rules for Ca-ion mobility. One candidate from the screening, W 2 O 3 (PO 4 ) 2 , is confirmed to have a low Nudged Elastic Band (NEB) barrier of 168 meV within a one-dimensional (1D) ion percolation topology. This candidate is subsequently synthesized and electrochemically tested, achieving reversible Ca cycling with a capacity of 25 mA h/g. To further accelerate the screening for promising Ca intercalation electrodes, machine learning (ML) Random Forest (RF) and Extreme Gradient Boosting (XGB) classification models are created with local environment descriptors based on a large, structurally and chemically diverse dataset of minimum energy pathways, spanning over 5,000 density functional theory (DFT) site energy calculations. Accuracies of 92% are achieved, material design metrics are quantified, ML force-fields are leveraged in an accelerated iteration of the screening, and a total of 27 novel Ca cathode materials are highlighted for further investigation.

25 ENERGY STORAGE↗
DOE OSTI2021

Quantum chemical calculations of lithium-ion battery electrolyte and interphase species

Abstract Lithium-ion batteries (LIBs) represent the state of the art in high-density energy storage. To further advance LIB technology, a fundamental understanding of the underlying chemical processes is required. In particular, the decomposition of electrolyte species and associated formation of the solid electrolyte interphase (SEI) is critical for LIB performance. However, SEI formation is poorly understood, in part due to insufficient exploration of the vast reactive space. The Lithium-Ion Battery Electrolyte (LIBE) dataset reported here aims to provide accurate first-principles data to improve the understanding of SEI species and associated reactions. The dataset was generated by fragmenting a set of principal molecules, including solvents, salts, and SEI products, and then selectively recombining a subset of the fragments. All candidate molecules were analyzed at the ω B97X-V/def2-TZVPPD/SMD level of theory at various charges and spin multiplicities. In total, LIBE contains structural, thermodynamic, and vibrational information on over 17,000 unique species. In addition to studies of reactivity in LIBs, this dataset may prove useful for machine learning of molecular and reaction properties.

25 ENERGY STORAGE↗