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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Covariant density functional theory input for r-process simulations in actinides and superheavy nuclei: The ground state and fission properties

A systematic investigation of the ground-state and fission properties of even-even actinides and superheavy nuclei with Z = 90–120 from the two-proton up to two-neutron drip lines with proper assessment of systematic theoretical uncertainties has been performed for the first time in the framework of covariant density functional theory (CDFT). Furthermore, these results provide a necessary theoretical input for the r-process modeling in heavy nuclei and, in particular, for the study of fission cycling. Four state-of-the-art globally tested covariant energy density functionals (CEDFs), namely, DD-PC1, DD-ME2, NL3*, and PC-PK1, representing the major classes of the CDFT models are employed in the present paper. Ground-state deformations, binding energies, two-neutron separation energies, α-decay Q α values and half-lives, and the heights of fission barriers have been calculated for all these nuclei. Theoretical uncertainties in these physical observables and their evolution as a function of proton and neutron numbers have been quantified and their major sources have been identified. Spherical shell closures at Z = 120, N = 184, and N = 258 and the structure of the single-particle (especially, high-j) states in their vicinities as well as nuclear matter properties of employed CEDFs are two major factors contributing to theoretical uncertainties. However, different physical observables are affected in a different way by these two factors. For example, theoretical uncertainties in calculated ground-state deformations are affected mostly by the former factor, while theoretical uncertainties in fission barriers depend on both of these factors.

190 ≤ A ≤ 219↗
DOE OSTI2025

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗
DOE OSTI2026

Isotopic discrepancy in microscopic fusion of 16 O + 40,48 Ca

We present a detailed comparison of fusion in 16 O + 40,48 Ca within a fully microscopic time-dependent Hartree-Fock framework over sub-, near-, and above-barrier energies. While 16 O + 40 Ca fusion cross sections are reproduced within the available experimental energy range, calculations for the neutron-rich 16 O + 48 Ca system systematically overestimate the measured near-barrier cross sections. The discrepancy persists across the tested energy density functionals, indicating that it is not tied to a particular parametrization. Analysis of the dynamical interaction potentials reveals an outward shift and softening of the barrier region at the lowest collision energies. These results are then compared to a frozen density description where the barrier modifications are absent and the predicted cross sections are brought closer to experiment. Here, this demonstrates that the observed discrepancy is closely linked to dynamical surface rearrangement and isovector effects, and highlights the sensitivity of low-energy fusion to both the structure and dynamics of neutron-rich nuclei. This suggests that limitations in the isovector sector of conventional Skyrme EDFs in describing the structure of 48 Ca, as highlighted by CREX, are also manifested in the fusion dynamics of the 16 O + 48 Ca system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electronic structure of van der Waals ferromagnet CrI 3 from self-consistent vertex corrected GW approaches

Electronic structure of the layered van der Waals ferromagnet CrI 3 is studied with self consistent diagrammatic approaches beyond the GW approximation. Considerable improvement in the calculated band gap as compared to the non-self-consistent G0W0 results has been found. Certain spectral features in the valence bands discovered recently by the angle resolved photoemission spectroscopy, are reproduced better when we use full frequency dependent self energy. Density functional theory and the quasiparticle self-consistent GW method which are based on the frequency-independent self energy are unable to resolve these features. The non-locality effects in the diagrams beyond the GW approximation are large for both polarizability and self energy. This finding can potentially have an impact on the development of methods like the GW+DMFT (dynamical mean field theory).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Weak decays in superheavy nuclei

Superheavy nuclei represent the extreme atoms and nuclides known at the limit of mass and charge. The observed superheavy nuclei are all proton-rich; they decay primarily by emitting 𝛼 particles and by fission with a possible small electron capture (EC) branch. Here, due to the huge atomic numbers and associated relativistic effects, EC decays of superheavy systems are expected to differ from what is known in lighter nuclei. In this letter, using the quantified relativistic nuclear density functional theory and the quasiparticle random-phase approximation with the interaction optimized to experimental EC/𝛽 ± -decay half-lives, and Gamow-Teller resonance energies, we study the EC/𝛽 ± -decays in 𝑍=101–118 nuclei. Both allowed (1 + ) and first-forbidden (0 − ,1 − and 2 − ) transitions are considered. We show that the first-forbidden 1 − transitions dominate the decay rates in almost all studied nuclei. For proton-rich nuclei, EC dominates over 𝛽 + decay. Based on calculations with two relativistic energy density functionals, we identify 45 candidate nuclei in which a competition between weak decays and 𝛼 decay and spontaneous fission is expected.

A ≥ 220↗
DOE OSTI2020

Electric and magnetic dipole strength in 112,114,116,118,120,124 Sn

Background: There is renewed interest in electric dipole strength distributions for a variety of reasons including the extraction of the dipole polarizability related to properties of the symmetry energy and a measure for the neutron skin thickness, understanding the structure of low-energy $\textit{E}$1 strength in nuclei with neutron excess, and establishing the systematics of the isovector giant dipole resonance (IVGDR). Inelastic proton scattering at energies of a few hundred MeV and very forward angles including 0° has been established as a tool for the study of electric and magnetic dipole strength distributions in nuclei. Purpose: The present work aims at a systematic investigation of the electric and magnetic dipole strength distributions in the chain of stable even-mass tin isotopes. Methods: Inelastic proton scattering experiments were performed at the Research Center for Nuclear Physics, Osaka, with a 295-MeV beam covering laboratory angles 0°–6° and excitation energies 6–22 MeV. Cross sections due to $\textit{E}$1 and $\textit{M}$1 excitations were extracted with a multipole decomposition analysis (MDA) and then converted to reduced transition probabilities with the “virtual photon method” for $\textit{E}$1 and the “unit cross section method” for $\textit{M}$1 excitations, respectively. Including a theory-aided correction for the high-excitation-energy region not covered experimentally, the electric dipole polarizability was determined from the $\textit{E}$1 strength distributions. Results: In this work, total photoabsorption cross sections derived from the $\textit{E}$1 and $\textit{M}$1 strength distributions show significant differences compared to those from previous ($\textit{γ, xn}$) experiments in the energy region of the IVGDR. The widths of the IVGDR deduced from the present data with a Lorentz parametrization show an approximately constant value of about 4.5 MeV in contrast to the large variations between isotopes observed in previous work. The IVGDR centroid energies are in good correspondence to expectations from empirical systematics of their mass dependence. Furthermore, a study of the dependence of the IVGDR energies on bulk matter properties is presented. The $\textit{E}$1 strengths below neutron threshold show fair agreement in this work with results from ($\textit{γ, γ'}$) experiments on 112,116,120,124 Sn in the energy region between 6 and 7 MeV, where also isoscalar $\textit{E}$1 strength was found for 124 Sn. At higher excitation energies, large differences are observed, pointing to a different nature of the excited states with small ground-state branching ratios. The isovector spin-$\textit{M}$1 strengths exhibit a broad distribution between 6 and 12 MeV in all studied nuclei. Conclusions: The present results contribute to the solution of a variety of nuclear structure problems including the systematics of the energy and width of the IVGDR, the structure of low-energy $\textit{E}$1 strength in nuclei, new constraints to energy density functionals (EDFs) aiming at a systematic description of the dipole polarizability across the nuclear chart, from which properties of the symmetry energy can be derived, and the systematics of the isovector spin-$\textit{M}$1 strength in heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Measurement of the $\textit{B}(\textit{E}2 ↑)$ strengths of 36 Ca and 38 Ca

In this work, the $B(E2, 0^+_1 → 2^+_1)$ strengths of 36 Ca and 38 Ca are measured to be 131(20) $e^2\text{fm}^4$ and 101(11) $e^2\text{fm}^4$, respectively. The $\textit{B}(\textit{E}2)$ value for 36 Ca required a measurement of the p/γ branching ratio because the 2 + state is proton unbound. This branching ratio is $B_p$=0.087(8). These $\textit{B}(\textit{E}2)$ and branching-ratio values can be reproduced in the shell-model with the ZMB2 interaction, an interaction that predicts the $\textit{Z} = 20 \textit{sd}$-shell closure is incomplete with large proton pf-shell occupancies in the ground state. These occupancies are at odds with other shell-model and energy-density-functional calculations of 36 Ca. New data are used to provide an update on constraints of the density dependence of the symmetry energy through mirror charge-radii differences as well as to help reduce uncertainties of the astrophysical important 35 K(p,γ) reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Pauli energy contribution to the nucleus-nucleus interaction

Background: The Pauli exclusion principle plays a crucial role as a building block of many-body quantal systems comprised of fermions. It also induces a “Pauli repulsion” in the interaction between di-nuclear systems. It has been shown in [Phys. Rev. C 95, 031601(R) (2017)] that the Pauli repulsion widens the nucleus-nucleus potential barrier, thus hindering sub-barrier fusion. Purpose: To investigate the proton and neutron contributions to the Pauli repulsion, both in the bare potential neglecting shape polarization and transfer between the reactants, as well as in the dynamical potential obtained by accounting for such dynamical rearrangements. Methods: As the basis of our study we utilize the Pauli kinetic energy (PKE) obtained by studying the nuclear localization function (NLF). Recently this approach has been generalized to incorporate all of the dynamical and time-odd terms present in the nuclear energy density functional. This approach is employed in the density constrained frozen Hartree-Fock (DCFHF) and in the density constrained time-dependent Hartree-Fock (DC- TDHF) microscopic methods. Results: The PKE spatial distribution shows that a repulsion occurs in the neck between the nuclei when they first touch. Inside the barrier, neutrons can contribute significantly more to the Pauli repulsion in neutron-rich systems. Dynamical effects tend to lower the Pauli repulsion near the barrier. Proton and neutron dynamical contributions to the PKE significantly differ inside the barrier for asymmetric collisions, which is interpreted as an effect of multinucleon transfer. Conclusions: The PKE is shown to make a significant contribution to nuclear interaction potentials. Protons and neutrons can play very different roles in both the bare potential and in the dynamical rearrangement. Further microscopic studies are required to better understand the role of transfer and to investigate the effect of pairing and deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

Nonlinear viscoelastic response of highly filled elastomers under multiaxial finite deformation

A biaxial tester was used to obtain precise biaxial stress responses of highly filled, high strain capability elastomers. Stress-relaxation experiments show that the time-dependent part of the relaxation response can be reasonably approximated by a function which is strain and biaxiality independent. Thus, isochronal data from the stress-relaxation curves can be used to determine the stored energy density function. The complex behavior of the elastomers under biaxial deformation may be caused by dewetting.

Peng, Steven T. J.↗
DOE OSTI2021

Microscopic origin of reflection-asymmetric nuclear shapes

The presence of nuclear ground states with stable reflection-asymmetric shapes is supported by rich experimental evidence. Theoretical surveys of odd-multipolarity deformations predict the existence of pear-shaped isotopes in several fairly localized regions of the nuclear landscape in the vicinity of near-lying single-particle shells with Δℓ = Δj = 3. We analyze the role of isoscalar, isovector, neutron-proton, neutron-neutron, and proton-proton multipole interaction energies in inducing the onset of reflection-asymmetric ground-state deformations. The calculations are performed in the framework of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory using two Skyrme energy-density functionals and density-dependent pairing force. We show that reflection-asymmetric ground-state shapes of atomic nuclei are driven by the odd-multipolarity neutron-proton (or isoscalar) part of the nuclear interaction energy. This result is consistent with the particle-vibration picture, in which the main driver of octupole instability is the isoscalar octupole-octupole interaction giving rise to large E3 polarizability. The necessary condition for the appearance of localized regions of pear-shaped nuclei in the nuclear landscape is the presence of parity doublets involving Δℓ = Δj = 3 proton or neutron single-particle shells. This condition alone is, however, not sufficient to determine whether pear shapes actually appear, and, if so, what the corresponding reflection-asymmetric deformation energies are. The predicted small reflection-asymmetric deformation energies result from dramatic cancellations between even- and odd-multipolarity components of the nuclear binding energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Linearized Pair-Density Functional Theory for Vertical Excitation Energies

Multiconfiguration pair-density functional theory (MC-PDFT) is a computationally efficient method that computes the energies of electronic states in a state specific or state average framework via an on-top functional. However, MC-PDFT does not include state interaction among these states since the final energies do not come from the diagonalization of an effective model-space Hamiltonian. Recently, multistate extensions such as linearized PDFT (L-PDFT) have been developed to accurately model the potentials near conical intersections and avoided crossings. However, there has not been any systematic study evaluating their performance for predicting vertical excitations at the equilibrium geometry of a molecule, when the excited states are generally well separated. In this paper, we report the performance of L-PDFT on the extensive QUESTDB data set of vertical excitations using a database of automatically selected active spaces. We show that L-PDFT performs well on all these excitations and successfully reproduces the performance of MC-PDFT. These results further demonstrate the potential of L-PDFT, as its scaling is constant with the number of states included in the state-average manifold, whereas MC-PDFT scales linearly in this regard.

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

GW190814: Impact of a 2.6 solar mass neutron star on the nucleonic equations of state

Is the secondary component of GW190814 the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system? This is the central question animating this paper. Covariant density functional theory provides a unique framework to investigate both the properties of finite nuclei and neutron stars, while enforcing causality at all densities. By tuning existing energy density functionals we were able to: (i) account for a 2.6 M ⊙ neutron star, (ii) satisfy the original constraint on the tidal deformability of a 1.4 M ⊙ neutron star, and (iii) reproduce ground-state properties of finite nuclei. Yet, for the class of models explored in this work, we find that the stiffening of the equation of state required to support supermassive neutron stars is inconsistent with either constraints obtained from energetic heavy-ion collisions or from the low deformability of medium-mass stars. Furthermore, we speculate that the maximum neutron star mass can not be significantly higher than the existing observational limit and that the 2.6 M ⊙ compact object is likely to be the lightest black hole ever discovered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Electron capture of superheavy nuclei with realistic lepton wave functions

The superheavy nuclei push the periodic table of the elements and the chart of the nuclides to their limits, providing a unique laboratory for studies of the electron-nucleus interactions. The most important weak decay mode in known superheavy nuclei is electron capture (EC). In the standard calculations of EC, the lepton wave functions are usually considered in the lowest-order approximation. In this work, we investigate the sensitivity of EC rates on the choice of the electron wave functions by (i) assuming the single-particle approximation for the electron wave functions, and (ii) carrying out Dirac-Hartree-Fock (DHF) calculations. The nuclear response is generated based on the state-of-the-art quasiparticle random phase approximation employing relativistic nuclear energy density functional theory. Here, we show that using the improved lepton wave functions reduces the EC rates up to 40% in the superheavy nucleus oganesson (𝑍=118). Interestingly, because of screening effects, the difference between the EC rates obtained with the DHF and single-particle calculations is fairly small.

Atomic orbital↗
DOE OSTI2021

Implications of PREX-2 on the Equation of State of Neutron-Rich Matter

Laboratory experiments sensitive to the equation of state of neutron rich matter in the vicinity of nuclear saturation density provide the first rung in a “density ladder” that connects terrestrial experiments to astronomical observations. In this context, the neutron skin thickness of 208 Pb ($R$ $^{208}_{skin}$) provides a stringent laboratory constraint on the density dependence of the symmetry energy. In turn, an improved value of $R$ $^{208}_{skin}$ has been reported recently by the PREX collaboration. Exploiting the strong correlation between $R$ $^{208}_{skin}$ and the slope of the symmetry energy L within a specific class of relativistic energy density functionals, we report a value of L = (106 ± 37) MeV—which systematically overestimates current limits based on both theoretical approaches and experimental measurements. Here, the impact of such a stiff symmetry energy on some critical neutron-star observables is also examined.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Landau–Devonshire analysis of strain effects on ferroelectric Al 1-x Sc x N

We present a thermodynamic analysis of the recently discovered nitride ferroelectric materials using the classic Landau–Devonshire approach. Electrostrictive and dielectric stiffness coefficients of Al 1-x Sc x N with a wurtzite structure (6 mm) are determined using a free energy density function assuming a hexagonal parent phase (6/mmm), with the first-order phase transition based on the dielectric stiffness relationships. The results of this analysis show that the strain sensitivity of the energy barrier is one order of magnitude larger than that of the spontaneous polarization in these wurtzite ferroelectrics, yet both are less sensitive to strain compared to classic perovskite ferroelectrics. These analysis results reported here explain experimentally reported sensitivity of the coercive field to elastic strain/stress in Al 1-x Sc x N films and would enable further thermodynamic analysis via phase field simulation and related methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Landau-Devonshire Approach of Elastic Strain Effects on Ferroelectric Al1-xScxN

Landau-Devonshire thermodynamic modeling enables the phenomenological description of ferroelectrics. 1 The approach has been validated across a wide range of ferroelectric materials from the first ferroelectric discovered, Rochelle Salt, to recent HfO2-based materials. 2-4 Extension to what is now referred to as the Landau-Ginsburg-Devonshire model enables rigorous investigation of polarization boundaries and reorientation, phase transformations, and domain formation/evolution. We present a thermodynamic analysis of the recently discovered nitride ferroelectric materials 5-7 using the classic Landau-Devonshire approach. The electrostrictive and dielectric stiffness coefficients of Al1-xScxN with wurtzite structure (6mm) are determined using a free energy density function assuming a hexagonal parent phase (6/mmm), with the first order phase transition based on the dielectric stiffness relationships (Fig. 1). The results of this analysis show that the strain sensitivity of the energy barrier is one order of magnitude larger than that of the spontaneous polarization in these wurtzite ferroelectrics, yet both are less sensitive to strain compared to classic perovskite ferroelectrics (Fig. 2). These analysis results reported here explain experimentally reported sensitivity of coercive field to elastic strain/stress in Al1-xScxN films, and would enable further thermodynamic analysis via phase field simulation and related methods.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2021

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. In this study, we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗