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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Density dependence of the symmetry energy in the post–PREX-CREX era

The recently published CREX results suggest a rather peculiar picture for the density dependence of the symmetry energy. Whereas PREX favors a large neutron skin thickness in 208 Pb, thereby suggesting a stiff equation of state, CREX suggests instead a much softer equation of state. Here, this discrepancy has caused a large spur in the theoretical community, since no model has been able to simultaneously reproduce within 1σ the PREX and CREX results. Motivated by a novel correlation between a CREX observable and a combination of bulk symmetry energy parameters, we calibrate three new covariant energy density functionals that reproduce binding energies and charge radii of spherical nuclei—and also accommodate the constraints imposed by PREX and CREX. Given that these models suggest a stiff equation of state at high densities, predictions for neutron star properties are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Density dependence of the symmetry energy in the post–PREX-CREX era

The recently published CREX results suggest a rather peculiar picture for the density dependence of the symmetry energy. Whereas PREX favors a large neutron skin thickness in 208 Pb, thereby suggesting a stiff equation of state, CREX suggests instead a much softer equation of state. This discrepancy has caused a large spur in the theoretical community, since no model has been able to simultaneously reproduce within 1σ the PREX and CREX results. Motivated by a novel correlation between a CREX observable and a combination of bulk symmetry energy parameters, we calibrate three new covariant energy density functionals that reproduce binding energies and charge radii of spherical nuclei—and also accommodate the constraints imposed by PREX and CREX. Given that these models suggest a stiff equation of state at high densities, predictions for neutron star properties are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Exploring and enhancing the accuracy of interior-scaled Perdew–Zunger self-interaction correction

The Perdew–Zunger self-interaction correction (PZ-SIC) improves the performance of density functional approximations for the properties that involve significant self-interaction error (SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by local scaling self-interaction correction (LSIC) of the PZ-SIC to the local spin density approximation (LSDA). Here, we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of the atomic number Z in the asymptotic expansion of the exchange–correlation (xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator z σ , which distinguishes one-electron regions from many-electron regions. LSIC+ applied to the LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA), and almost close to the strongly constrained and appropriately normed (SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled down SIC to functionals beyond the LSDA requires a gauge transformation of the functional’s energy density. Here, the resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds.

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

Global optimization of harmonic oscillator basis in covariant density functional theory

The present investigation focuses on the improvement of the accuracy of the description of binding energies within moderately sized fermionic basis. Using the solutions corresponding to infinite fermionic basis it was shown that in the case of meson exchange (ME) covariant energy density functionals (CEDFs) the global accuracy of the description of binding energies in the finite $N_F$ = 16 - 20 bases can be drastically (by a factor ranging from ~3 up to ~9 dependent on the functional and $N_F$) improved by a global optimization of oscillator frequency of the basis. This is a consequence of the unique feature of the ME functionals in which with increasing fermionic basis size fermionic and mesonic energies approach the exact (infinite basis) solution from above and below, respectively. As a consequence, an optimal oscillator frequency $\hbar\omega_0$ of the basis can be defined which provides an accurate reproduction of exact total binding energies by the ones calculated in truncated basis. This leads to a very high accuracy of the calculations in moderately sized $N_F=20$ basis when mass dependent oscillator frequency is used: global rms differences $\delta B_{rms}$ between the binding energies calculated in infinite and truncated bases are only 0.025 MeV and 0.031 MeV for the NL5(Z) and DD-MEZ functionals, respectively. Optimized values of the oscillator frequency $\hbar\omega_0$ are provided for three major classes of CEDFs, i.e. for density dependent meson exchange functionals, nonlinear meson exchange ones and point coupling functionals.

Binding energy & masses↗
DOE OSTI2021

Semi-classical approximation description of static properties of nuclei

In this report, we present recent development in semi-classical description of static properties of nuclei. By employing the Wigner transform, we derive simple semiclassical approximations for evaluating properties of nuclear system. We present results of calculations of static properties, using the energy density functional associated with Skyrme effective nucleon–nucleon interaction. In particular, we consider properties of the Wigner distribution function, line of beta stability, Coulomb and symmetry energies, deformation energies, equation of state, nuclear radii, and single-particle level density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Pauli potential formalism at finite temperature

At zero temperature, the Pauli potential—the functional derivative of the Pauli kinetic energy density functional—is the key to the accuracy of the orbital-free density functional theory (OFDFT) as it is supposed to capture all the effects associated with the Pauli exclusion principle. Here, we extend this concept to finite temperature by defining Pauli free energy and the modified Pauli free energy, both representing the natural generalizations of the Pauli term from zero- T to finite- T . We discuss their physical interpretation, the mathematical nuances, and the applicability, arguing that the modified Pauli potential should be used as an extension of the zero- T counterpart within the OF-DFT framework. Through analytical and numerical methods, we then analyze some of the exact properties concerning the modified Pauli free- and kinetic-energy terms and examine the temperature dependence of the modified Pauli potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Influence of the symmetry energy on the nuclear binding energies and the neutron drip line position

A clear connection can be established between properties of nuclear matter and finite-nuclei observables, such as the correlation between the slope of the symmetry energy and the dipole polarizability, or between compressibility and the isoscalar monopole giant resonance excitation energy. Establishing a connection between realistic atomic nuclei and an idealized infinite nuclear matter leads to a better understanding of underlying physical mechanisms that govern nuclear dynamics. In this work, we aim to study the dependence of the binding energies and related quantities (e.g., location of drip lines, the total number of bound even-even nuclei) on the symmetry energy S 2 (ρ). The properties of finite nuclei are calculated by employing the relativistic Hartree-Bogoliubov model, assuming even-even axial and reflection symmetric nuclei. Calculations are performed by employing two families of relativistic energy density functionals, based on different effective Lagrangians, constrained to a specific symmetry energy at the saturation density J within the interval of 30–36 MeV. Nuclear binding energies and related quantities of bound nuclei are calculated between 8 ≤ Z ≤ 104 from the two-proton to the two-neutron drip line. As the neutron drip line is approached, the interactions with stiffer J tend to predict more bound nuclei, resulting in a systematic shift of the two-neutron drip line towards more neutron-rich nuclei. Consequentially, a correlation between the number of bound nuclei N nucl and S 2 (ρ) is established for a set of functionals constrained using the similar optimization procedures. Furthermore the direction of the relationship between the number of bound nuclei and the symmetry energy highly depends on the density under consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Free-energy orbital-free density functional theory: recent developments, perspective, and outlook

By summarizing the constraint-based development of orbital-free free-energy density functional approximations, we provide a perspective on progress over the last 15 years, the limitations of existing functionals, and the challenges awaiting resolution. We outline the chronology of the development of noninteracting and exchange-correlation free-energy orbital-free functionals and summarize the theoretical basis of existing local density approximation, second-order approximation, generalized gradient approximation (GGA), and meta-GGAs. We discuss limitations and challenges such as problems with thermodynamic derivatives, free-energy nonadditivity and the closely related issue of all-electron versus valence-only local pseudo-potential performance.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Cluster model of C 12 in the density functional theory framework

We employ the constrained density functional theory to investigate cluster phenomena for the 12 C nucleus. The proton and neutron densities are generated from the placement of three 4 He nuclei (α particles) geometrically. These densities are then used in a density constrained Hartree-Fock calculation that produces an antisymmetrized state with the same densities through energy minimization. In the calculations no a priori analytic form for the single-particle states is assumed and the full energy density functional is utilized. The geometrical scan of the energy landscape provides the ground state of 12 C as an equilateral triangular configuration of three αs with molecular bond like structures. The use of the nucleon localization function provides further insight to these configurations. One can conclude that these configurations are a hybrid between a pure mean-field and a pure α particle condensate. Furthermore, this development could facilitate density functional theory based fusion calculations with a more realistic 12 C ground state.

6 ≤ A ≤ 19↗
DOE OSTI2024

Isospin composition of fission barriers

We employ a microscopic method to study how isospin affects the fission potential of 240 Pu. Our approach uses constrained Hartree-Fock theory (CHF) which allows us to separately investigate the isoscalar and isovector properties of the nuclear energy density functional (EDF). By analyzing the isoscalar and isovector components of the EDF along the fission path we can assess the isovector contribution to fission barriers. Here, we study this effect for the fully adiabatic path to scission. The isovector component of the fission potential is found to increase in magnitude as the nucleus evolves towards scission, exemplifying the importance of stringent constraints on the isovector sector of the nuclear EDF for reliable predictions of fission properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multimodal fission from self-consistent calculations

When multiple fission modes coexist in a given nucleus, distinct fragment yield distributions appear. Multimodal fission has been observed in a number of fissioning nuclei spanning the nuclear chart, and this phenomenon is expected to affect the nuclear abundances synthesized during the rapid neutron-capture process (𝑟-process). In this study, we generalize the previously proposed hybrid model for fission-fragment yield distributions to predict competing fission modes and estimate the resulting yield distributions. Here, our framework allows for a comprehensive large-scale calculation of fission-fragment yields suited for 𝑟-process nuclear network studies. Nuclear density functional theory is employed to obtain the potential energy and collective inertia tensor on a multidimensional collective space defined by mass multipole moments. Fission pathways and their relative probabilities are determined using the nudged elastic band method. Based on this information, mass and charge fission yields are predicted using the recently developed hybrid model. Fission properties of fermium isotopes are calculated in the axial quadrupole-octupole collective space for three energy density functionals (EDFs). Disagreement between the EDFs appears when multiple fission modes are present. Within our framework, the UNEDF⁢1 HFB EDF agrees best with experimental data. Calculations in the axial quadrupole-octupole-hexadecapole collective space improve the agreement with the experiment for SkM*. We also discuss the sensitivity of fission predictions on the choice of EDF for several superheavy nuclei. Fission-fragment yield predictions for nuclei with multiple fission modes are sensitive to the underlying EDF. For large-scale calculations in which a minimal number of collective coordinates is considered, UNEDF⁢1 HFB provides the best description of experimental data, though the sensitivity motivates robust quantification of the uncertainties of the theoretical model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Applications of reduced-basis methods to the nuclear single-particle spectrum

Reduced-basis methods provide a powerful framework for building efficient and accurate emulators. Although widely applied in many fields to simplify complex models, reduced-basis methods have only been recently introduced into nuclear physics. In this Letter we build an emulator to study the single-particle structure of atomic nuclei. By scaling a suitable mean-field Hamiltonian, a “universal” reduced basis is constructed capable of accurately and efficiently reproduce the entire single-particle spectrum of a variety of nuclei. Indeed, the reduced-basis model reproduces both ground- and excited-state energies as well as the associated wave functions with remarkable accuracy. Here our results bode well for more demanding applications that use Bayesian optimization to calibrate nuclear energy density functionals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Information content of the differences in the charge radii of mirror nuclei

Differences in the charge radii of mirror nuclei have been recently suggested to contain information on the slope of the symmetry energy L. To test this hypothesis, we perform statistical correlation analysis using quantified energy density functionals that are consistent with our previous knowledge on global nuclear observables such as binding energies and charge radii. We conclude that the difference in charge radii between a mirror pair, $ΔR$$^{mir}_{ch}$, is an inferior isovector indicator compared to other observables, such at the neutron skin or electric dipole polarizability α D . In particular, this quantity correlates poorly with both the neutron skin and L. Here, we demonstrate that $ΔR$$^{mir}_{ch}$ is influenced by pairing correlations in the presence of low-lying proton continuum in the proton-rich mirror-partner nucleus. Considering the large theoretical uncertainties on $ΔR$$^{mir}_{ch}$, we conclude that the precise data on mirror charge radii cannot provide a stringent constraint on L.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Precision Determination of the Neutral Weak Form Factor of Ca 48

In this work, we report a precise measurement of the parity-violating asymmetry $A_{\rm PV}$ in the elastic scattering of longitudinally polarized electrons from $^{48}{\rm Ca}$. We measure $A_{\rm PV} =2668\pm 106\ {\rm (stat)}\pm 40\ {\rm (syst)}$ parts per billion, leading to an extraction of the neutral weak form factor $F_{\rm W} (q=0.8733$ fm$^{-1}) = 0.1304 \pm 0.0052 \ {\rm (stat)}\pm 0.0020\ {\rm (syst)}$ and the charge minus the weak form factor $F_{\rm ch} - F_{\rm W} = 0.0277\pm 0.0055$. The resulting neutron skin thickness $R_n-R_p=0.121 \pm 0.026\ {\rm (exp)} \pm 0.024\ {\rm (model)}$~fm is relatively thin yet consistent with many model calculations. The combined CREX and PREX results will have implications for future energy density functional calculations and on the density dependence of the symmetry energy of nuclear matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Global study of separable pairing interaction in covariant density functional theory

A systematic global investigation of pairing properties based on all available experimental data on pairing indicators has been performed for the first time in the framework of covariant density functional theory. It is based on the separable pairing interaction. The optimization of the scaling factors of this interaction to experimental data clearly reveals its isospin dependence in the neutron subsystem. However, the situation is less certain in the proton subsystem since similar accuracy of the description of pairing indicators can be achieved both with isospin-dependent and mass-dependent scaling factors. The differences in the functional dependencies of scaling factors lead to uncertainties in the prediction of proton and neutron pairing properties which are especially pronounced at high isospin and could have a significant impact on some physical observables. For a given part of the nuclear chart the scaling factors for spherical nuclei are smaller than those for deformed ones; this feature exists also in nonrelativistic density functional theories. Furthermore, its origin is traced back to particle-vibration coupling in odd-A nuclei which is missing in all existing global studies of pairing. Although the present investigation is based on the NL5(E) covariant energy density functional (CEDF), its general conclusions are expected to be valid also for other CEDFs built at the Hartree level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗