Search NASA⌕ Search

SEARCH · Search NASA

Results for “relativistic density functional theory”

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 OSTI2025

Green’s Function Perspective on the Nonlinear Density Response of Quantum Many-Body Systems

We derive equations of motion for higher order density response functions using the theory of thermodynamic Green’s functions. We also derive expressions for the higher order generalized dielectric functions and polarization functions. Moreover, we relate higher order response functions and higher order collision integrals within the Martin–Schwinger hierarchy. We expect our results to be highly relevant to the study of a variety of quantum many-body systems such as matter under extreme temperatures, densities, and pressures.

Density Functional Theory↗
DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Emerging edge states on the surface of the epitaxial semimetal CuMnAs thin film

Epitaxial thin films of CuMnAs have recently attracted attention due to their potential to host relativistic antiferromagnetic spintronics and exotic topological physics. Here, we report on the structural and electronic properties of a tetragonal CuMnAs thin film studied using scanning tunneling microscopy (STM) and density functional theory (DFT). STM reveals a surface terminated by As atoms, with the expected semi-metallic behavior. An unexpected zigzag step edge surface reconstruction is observed with emerging electronic states below the Fermi energy. DFT calculations indicate that the step edge reconstruction can be attributed to an As deficiency that results in changes in the density of states of the remaining As atoms at the step edge. Furthermore, this understanding of the surface structure and step edges on the CuMnAs thin film will enable in-depth studies of its topological properties and magnetism.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Neutron skins probed in proton knockout from neutron-rich nuclei

Proton-induced quasifree knockout reactions provide a powerful probe of nuclear single-particle structure and reaction dynamics in both stable and neutron-rich nuclei. Here, in this work, we develop a unified theoretical framework for the calculation of inclusive (𝑝, 2⁢𝑝) and sequential (𝑝, 3⁢𝑝) reaction cross sections and fragment momentum distributions at intermediate and relativistic energies. The approach is based on a probabilistic extension of Glauber multiple-scattering theory combined with microscopic nuclear densities obtained from Hartree-Fock-Bogoliubov calculations using Skyrme energy-density functionals. We focus in particular on the sensitivity of total cross sections and longitudinal momentum dispersions to neutron-skin thickness along isotopic chains. Our results indicate that both (𝑝, 2⁢𝑝) and (𝑝, 3⁢𝑝) reactions exhibit a systematic decrease of cross section and momentum width with increasing neutron excess, reflecting enhanced attenuation and surface bias induced by neutron skins. The effect is significantly stronger for two-proton removal, suggesting that (𝑝, 3⁢𝑝) reactions may offer enhanced sensitivity to isovector nuclear structure. These findings establish proton-induced knockout reactions as complementary hadronic probes of neutron skins and the density dependence of the nuclear symmetry energy.

direct reactions↗
DOE OSTI2024

Global properties of nuclei at finite-temperature within the covariant energy density functional theory

In stellar environments nuclei appear at finite temperatures, becoming extremely hot in core-collapse supernovae and neutron-star mergers. However, due to theoretical and computational complexity, most model calculations of nuclear properties are performed at zero temperature, while those existing at finite temperatures are limited only to selected regions of the nuclide chart. In this study we perform the global calculation of nuclear properties for even-even 8 ≤ Z ≤ 104 nuclei at temperatures in range 0 ≤ T ≤ 2 MeV. Calculations are based on the finite-temperature relativistic Hartree-Bogoliubov model supplemented by the Bonche-Levit-Vautherin vapor subtraction procedure. We find that near the neutron-drip line the continuum states have significant contribution already at moderate temperature T ≈ 1 MeV, thus emphasizing the necessity of the vapor subtraction procedure. Results include neutron emission lifetimes, quadrupole deformations, neutron-skin thickness, proton and neutron pairing gaps, entropy and excitation energy. Up to the temperature T ≈ 1 MeV, the nuclear landscape is influenced only moderately by the finite-temperature effects, mainly by reducing the pairing correlations. Here, as the temperature increases further, the effects on nuclear structures become pronounced, reducing both the deformations and the shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Dirac's equation and its implications for density functional theory based calculations of materials containing heavy elements

Electronic structure calculations based on density functional theory (DFT) give quantitatively accurate predictions of properties of most materials containing light elements. For heavy materials, and in particular for f -electron systems, DFT based methods can fail both qualitatively and quantitatively for two distinct reasons: their failure to describe confinement effects arising from localized f -electron behavior and their incomplete or approximate treatment of relativity. In addition, different methods for incorporating relativistic effects, which give identical results in most light materials, can give different predictions in heavy elements. In order to develop a quantitative capability for calculating the properties of these materials, it is essential to separate the predictions of the underlying equations from the uncertainty introduced in approximations used in computation. Working toward that goal, here we have developed a code, called dirac-fp , which is based directly on solving the Dirac-Kohn-Sham equations and uses the full-potential linear muffin-tin orbital (FP-LMTO) approach to electronic structure. In order to assess the performance of dirac-fp , we perform calculations on three different face-centered cubic materials using different approximate treatments of relativity: the scalar relativistic (SR) approach commonly used in most solid-state DFT codes, the scalar relativistic plus spin-orbit coupling corrections (SR+SO) approach which includes spin-orbit coupling self-consistently using the SR states inside the muffin tins, and the Dirac-Kohn-Sham (Dirac) approach implemented in dirac-fp . Performing calculations on thorium, in which relativistic effects should be strong, aluminum, in which relativistic effects should be negligible, and gold, in which relativistic effects play an intermediate role, we find that the Dirac approach is able to provide theoretically consistent results in the electronic structure and ground-state properties across all three materials.

3-dimensional systems↗
DOE OSTI2020

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to leading order (–0.221 Z 5/3 and –0.021 Z ln Z respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to E x (N,Z) ≈ –0.354 N 2/3 Z (as known before only for Z >> N >> 1) and E c ≈ –0.02 N ln N . These asymptotic estimates are most correct for atomic ions with large N and Z >> N, but we find that they are qualitatively and semi- quantitatively correct even for small N and for N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. Here, it is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nuclear charge densities in spherical and deformed nuclei: Toward precise calculations of charge radii

Precise measurements of atomic transitions affected by electron-nucleus hyperfine interactions offer sensitivity to explore basic properties of the atomic nucleus and study fundamental symmetries, including the search for new physics beyond the standard model of particle physics. Specifically, such measurements, augmented by atomic and nuclear calculations, will enable extraction of the higher-order radial moments of the charge-density distribution in spherical and deformed nuclei. The new data impose higher precision requirements on a theoretical description. The nuclear charge density is composed of the proton point distribution folded with the nucleonic charge distributions. The latter induce subtle relativistic corrections due to the coupling of nucleon magnetic moments with the nuclear spin-orbit density. Additional corrections come from the effect of center-of-mass projection. We assess the precision of nuclear charge density calculations by studying the behavior of relativistic and center-of-mass motion corrections to the second and fourth charge radial moments. Special attention has been paid to the magnetic spin-orbit density associated with the local variations of the spin-orbit current. The calculations for semimagic and open-shell nuclei are performed in the framework of self-consistent mean-field theory using quantified energy density functionals and density-dependent pairing forces. We used the general expression for the spin-orbit form factor that is valid for spherical and deformed nuclei. We studied the impact of various correction terms on the charge radii, fourth radial moments, diffraction radii, and surface thickness of spherical and deformed nuclei. The spin-orbit corrections to charge radial moments and surface thickness show strong shell fluctuations which can make an appreciable effect when aiming at high-precision predictions of isotopic shifts. The inclusion of relativistic and center-of-mass corrections impacts the quality of energy density functionals optimized to charge radii data. To establish reliable constraints on the existence of new forces from isotope shift measurements, precise calculations of nuclear charge densities of deformed nuclei are needed. The proper inclusion of the spin-orbit charge density and other correction terms is essential when aiming at extraction of subtle effects which become particularly visible in isotopic trends. It is also important when developing high-quality nuclear energy density functionals optimized using heterogeneous datasets involving absolute charge radii, differential charge radii, and charge form factor properties deduced from electron-scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Electronic and vibrational properties of bulk Cr 2 Ge 2 Te 6 from first-principles calculations

Here we use density functional theory based ab initio calculations to investigate the structural, vibrational, magnetic, and electronic properties of the layered ferromagnet Cr 2 Ge 2 Te 6 (CGT) that has attracted attention for potential spintronic applications. We optimize the structure for a fixed c/a ratio. Our results are in excellent agreement with experimental data on structure, phonons, and electronic properties. The use of van der Waals interactions and relativistic spin-orbit coupling yields accurate lattice constants and interlayer distances. Simultaneously, the adopted theoretical methods lead to a rigorous description of the vibrational normal modes, as well as the valence bands, that are in excellent agreement with Raman spectroscopy and angle-resolved photoemission spectroscopy spectra, respectively. The magnetic moment is slightly overestimated, and the magnetic anisotropy has the correct sign but is greater in magnitude than that in experiment. We find that inclusion of the on-site Coulomb repulsion on Cr−3d orbitals worsens agreement with experiment, in contrast to previous studies. Our results provide a path toward ab initio analysis of magnetic heterostructures based on CGT.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Spin alignment of vector mesons by glasma fields

We explain how spin alignment of vector mesons can be induced by background fields, such as electromagnetic fields or soft gluon fields. Our study is based on the quantum kinetic theory of spinning quarks and antiquarks and incorporates the relaxation of the dynamically generated spin polarization. The spin density matrix of vector mesons is obtained by quark coalescence via the Wigner function and kinetic equation. Our approach predicts a local spin correlation that is distinct from the nonlocal expressions previously obtained in phenomenological derivations. We estimate the magnitude of such local correlations in the glasma model of the preequilibrium phase of relativistic heavy ion collisions. It is found that the resulting spin alignment could be greatly enhanced and may be comparable to the experimental measurement in order of magnitude. We further propose new phenomenological scenarios to qualitatively explain the transverse-momentum and centrality dependence of spin alignment in a self-consistent framework.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

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

Comprehensive simulation of heavy-ion collisions at nonzero baryon chemical potential

Here, we present results of hydrodynamic modeling of Au-Au collisions from $\sqrt{s_{NN}}$ = 7.7 to 200 GeV. Our simulations have three novel components. First, we use a Linear EXtrapolation of Ultratrelativistic nucleon-nucleon Scattering to nucleus-nucleus collisions (LEXUS) inspired Monte Carlo initial-state model. Second, we use a crossover equation of state at finite baryon densities without a critical point. Finally, we use departure functions derived from the quasiparticle theory of transport coefficients for hadronic matter at nonzero baryon densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Impact of symmetry breaking and spin-orbit coupling on the band gap of halide perovskites

Halide perovskite (HP) materials have recently emerged as a class of semiconductors with immense promise for various optoelectronic applications, ranging from solar cells to light-emitting diodes. One of the unique attributes of HPs is their tunable band gaps with different factors governing their value. Furthermore, the first factor is related to relativistic corrections [“mass-Darwin,” connected to the ns 2 lone pairs, and spin-orbit coupling (SOC)] that induce an orbital shift or degeneracy splitting, resulting in a band-gap reduction. The second factor involves the structural configuration: in HPs the local symmetry of each Wyckoff position tends to be broken, inducing an opening of the band gap. Based on high-throughput density functional theory calculations, this paper systematically studies a possible self-cancelation on the band-gap correction for HPs when the polymorphous configuration—structural effects—and the SOC—electronic effects—are included. Our results indicate that the nature of interplay between SOC and symmetry breaking (SB) is that they are independent decoupling effects to describe the band-gap magnitude in halide perovskites. As a result of that, we observe a transitivity of the band-gap description; i.e., if we know the band gap of halide perovskites without SB and SOC, we can independently add the effects of band-gap reduction due to SOC and band-gap opening due to SB, regardless of the order in which these effects are considered.

14 SOLAR ENERGY↗
NASA NTRS1994

X-ray and gamma-ray emission of Sagittarius A* as a wind-accreting black hole

If, as many believe, Sgr A* is a massive black hole at the Galactic center, one should expect it to be a source of X-ray and gamma-ray activity, behaving basically as a scaled-down active galactic nucleus. An unavoidable source of accretion is the wind from IRS 16, a nearby group of hot, massive stars. Since the density and velocity of the accreting matter are known from observations, the accretion rate is basically a function of the putative black hole mass, M(sub h), only; this value represents a reliable lower limit to a real rate, given the other possible sources of accreting matter. Based on this and on the theories about shock acceleration in active galactic nuclei, we have estimated the expected production of relativistic particles and their hard radiation. These values turn out to be a function of M(sub h) as well. Comparing our results with available X-ray and gamma-ray observations which show Sgr A* to have a relatively low activity level, we conclude tentatively that the putative black hole in the Galactic center cannot have a mass greater than approximately 6 x 10(exp 3) solar mass. This conclusion is consistent with the upper limits to the black hole mass found by different methods earlier, although much more work is needed to make calculations of shock acceleration around black holes more reliable.

Mastichiadis, A.↗
DOE OSTI2022

Investigating effects of relativistic kinematics, dimensionality, interactions, and short-range correlations on the ratio of quartic over quadratic nuclear symmetry energies

While ample evidence for the so-called empirical parabolic law of the Equation of State (EOS) of isospin asymmetric nuclear matter (ANM) has been obtained in many studies within both non-relativistic and relativistic nuclear many-body theories using various interactions, it has been unclear if there is any fundamental physics reason for the small quartic symmetry energy compared to the quadratic one even as the ANM approaches pure neutron matter. Within both relativistic and non-relativistic Free Fermi Gas (FFG) models in coordinate spaces of arbitrary dimension d with and without considering Short-Range Correlations (SRC) as well as non-linear Relativistic Mean Field (RMF) models, we study effects of relativistic kinematics, dimensionality, interactions and SRC on the ratio Ψ(ρ) of quartic over quadratic symmetry energies in ANM EOSs. We found that the ratio Ψ(ρ) in the FFG model depends strongly on the dimension d. While it is very small already in the normal 3D space, it could be even smaller in spaces with reduced dimensions for sub-systems of particles in heavy-ion reactions and/or whole neutron stars due to constraints, collectivities and/or symmetries. Here, we also found that the ratio Ψ(ρ) could theoretically become very large only at the ultra-relativistic limit far above the density reachable in neutron stars. On the other hand, nuclear interaction directly and/or indirectly through SRC-induced high-momentum nucleons affect significantly the density dependence of Ψ(ρ) compared to the relativistic FFG model prediction. The SRC affects significantly not only the kinetic energy of symmetric nuclear matter but also the ratio Ψ(ρ) while the relativistic corrections are found negligible. The results may help better understand the EOS of dense neutron-rich matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Temperature-dependent lattice dynamics in iridium

The characterization of simple elemental systems is key to benchmarking first-principles modeling of electronic and vibrational behaviors of materials. A large body of literature has been built for most elemental systems in the periodic table. However, surprisingly little neutron work has been performed to probe the vibrational properties of iridium, likely due to its large neutron absorption cross section. Nonetheless, iridium is of significant scientific and technological interest due to large relativistic electron effects and electron-phonon coupling, particularly in strongly correlated iridate compounds. In this paper, we report temperature-dependent inelastic neutron scattering measurements of the vibrational properties of iridium, from which we extract key thermodynamic properties. To overcome the challenge of the large neutron absorption of iridium, we developed a simple post-processing correction procedure. The measured densities of phonon states compare well with quasiharmonic density functional theory calculations, although the obtained experimental phonon Grüneisen parameters are much larger than expected, reaching as high as γ=4.5, indicating substantial anharmonicity.

36 MATERIALS SCIENCE↗