SEARCH · Engineering Papers

Results for “random phase approximation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Quasiparticle self-consistent GW band structures and phase transitions of LiAlO 2 in tetrahedrally and octahedrally coordinated structures

In this work, a first-principles computational study is presented of various phases of LiAlO 2 . The relative total energies and equations of state of the α, β, and γ phases are determined after structural relaxation of each phase. The β and γ tetrahedral phases are found to be very close in energy and lattice volume with the γ phase having the lowest energy. The octahedral α phase is a high-pressure phase and the transition pressure from the γ and β phases to α is determined to be about 1 GPa. The electronic band structures of each phase at their own equilibrium volume are determined using the quasiparticle self-consistent (QS) GW method as well as using the 0.8 Σ approach in which the QS GW self-energy is reduced by a factor of 0.8 to correct for the underscreening of W in QS GW. The effective masses of the band edges and the nature of the band gaps are presented. The lowest energy γ phase is found to have a pseudodirect gap of 7.69 eV. The gap is direct at Γ but corresponds to a dipole forbidden transition. The imaginary part of the dielectric function and the absorption coefficient are calculated in the long-wavelength limit and the random phase approximation, without local field or electron-hole interaction effects for each phase, and their anisotropies are discussed. Si doping on the Al site is investigated as a possible n-type dopant in γ-LiAlO 2 using a 128-atom supercell corresponding to 3.125% Si on the Al sublattice in the generalized gradient approximation and a smaller 16-atom cell with 25% Si in the QS GW approximation. The Si is found to significantly perturb the conduction band and lower the gap but a clearly separated deep donor defect level is not found. However, the donor binding energy is still expected to be relatively deep, on the order of a few tenths eV in the hydrogenic effective mass approximation.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Large-scale calculations of 𝛽-decay rates and implications for 𝑟-process nucleosynthesis

Nuclear 𝛽 decay is a key element of the astrophysical rapid neutron capture process (𝑟 process). In this work, we present state-of-the-art global 𝛽-decay calculations based on the quantified relativistic nuclear energy density functional theory and the deformed proton-neutron quasiparticle random-phase approximation. Our analysis considers contributions from allowed and first-forbidden transitions. We used two point-coupling functionals with carefully calibrated time-odd terms and isoscalar pairing strength. The new calculations display consistent results for both employed functionals, especially near the neutron drip line, suggesting slower 𝛽 decays past the 𝑁=126 neutron shell closure than in commonly used 𝛽-decay models. The new rates, along with the existing rates based on the recent nonrelativistic global calculations, are found to slow down the synthesis of heavy elements in the 𝑟 process and significantly reduce the contribution of neutron-induced fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Cooperative spin coupling in PrFeO3

A transition from crystal-field to spin wave excitations has been experimentally observed in PrFeO3 using inelastic neutron scattering. Theoretical modeling with density functional theory and the mean field random phase approximation successfully reproduces the experimental observations and allows extraction of the values for the various exchange interactions. The study reveals that the strong coupling between the two magnetic sublattices of Pr and Fe is the fundamental mechanism underlying the low-temperature Pr magnetism. This coupling, analogous to the Dicke cooperative interaction between atoms and photons, polarizes the Pr moments, enhances the Pr-Pr interactions, and contributes to the formation of a quasidoublet, leading to the transition from single-particle to collective excitations within the Pr sublattice. Furthermore, this coupling induces spin reorientation on the Fe sublattice. The manipulation of magnetic moments via spin waves is achievable through the coupling between different magnetic sublattices in a bulk material containing multiple magnetic ions.

Garlea, Vasile [ORNL] (ORCID:0000000253227271)↗

DOE OSTI2026

Experimental confirmation of the magnetic ordering transition induced by an electronic structure change in the metallic triangular antiferromagnet Co 1/3⁢ TaS 2

We report angle-resolved photoemission spectroscopy (ARPES) studies combined with DFT+DMFT calculations to confirm that the magnetic ordering vector transition from Q=(1/2,0,0) to Q=(1/3,0,0) in the metallic triangular antiferromagnets Co 1/3±𝜀⁢ TaS 2 (𝜀 ≈ 0.007) is induced by the electronic structure change in the system. The ARPES-measured Fermi surface (FS) maps of Co 0.325 ⁢TaS 2 show two hexagonal and one circular holelike FSs around Γ, which matches well with the triple- Q state by taking into account the contribution of nesting vectors occurring between Co 3⁢𝑑 and Ta 5⁢𝑑 orbitals. In the case of Co 0.340 ⁢TaS 2 , a new electron pocket around K appears, and the FS geometry changes as a result of the correlation effect of Co 4 ⁢S 18 tripods forming in the system. The magnetic susceptibility calculations based on the charge-self-consistent DFT+DMFT band structures and the random phase approximation indicate that the most stable magnetic ordering vector (1/2,0,0) split into (1/6,0,0) and (1/2,0,0), which is consistent with the magnetic phase transition around 𝑥=1/3 in Co 𝑥 ⁢TaS 2 .

Fermi surface↗

DOE OSTI2020

First-principles study of magnetism and electric field effects in 2D systems

This review article provides a bird's-eye view of what first-principles based methods can contribute to next-generation device design and simulation. After a brief overview of methods and capabilities in the area, the authors focus on published work by their group since 2015 and current work on CrI3. The authors introduce both single- and dual-gate models in the framework of density functional theory and the constrained random phase approximation in estimating the Hubbard U for 2D systems vs their 3D counterparts. A wide range of systems, including graphene-based heterogeneous systems, transition metal dichalcogenides, and topological insulators, and a rich array of physical phenomena, including the macroscopic origin of polarization, field effects on magnetic order, interface state resonance induced peak in transmission coefficients, spin filtration, etc., are covered. For CrI3, the authors present their new results on bilayer systems such as the interplay between stacking and magnetic order, pressure dependence, and electric field induced magnetic phase transitions. The authors find that a bare bilayer CrI3, graphene|bilayer CrI3|graphene, hexagonal boron nitride (h-BN)|bilayer CrI3|h-BN, and h-BN|bilayer CrI3|graphene all have a different response at high field, while at small field, the difference is small except for graphene|bilayer CrI3|graphene. The authors conclude with discussion of some ongoing work and work planned in the near future, with the inclusion of further method development and applications.

Cheng, Hai-Ping (ORCID:0000000159901725)↗

DOE OSTI2020

Isoscalar and Isovector Giant Resonances in Closed Shells Nuclei and Bulk Properties of Nuclear Matter

Centroid energies, E CEN , of the isoscalar (T = 0) and isovector (T = 1) giant resonances of multipolarities L = 0 – 3 in 40,48 Ca, 68 Ni, 90 Zr, 116 Sn, 144 Sm and 208 Pb, were calculated within the fully self-consistent spherical Hartree-Fock (HF)-based random phase approximation (RPA) theory, using 33 different energy density functionals associated with Skyrme-type effective nucleon-nucleon interactions of the standard form commonly employed in the literature. Here, we also calculate the Pearson linear correlation coefficients between each E CEN and each bulk property of nuclear matter (NM), associated with the Skyrme interactions used in the calculations, and determine the sensitivity of E CEN to bulk properties of NM. By comparing the calculated values of E CEN to the experimental data, we constrain the values of the bulk NM properties. We find that interactions associated with the values of the NM effective mass, m* / m = 0.70 to 0.90, incompressibility coefficient, K NM = 210 to 240 MeV, and the enhancement coefficient of the energy weighted sum rule of the isovector giant dipole resonance, κ = 0.25 to 0.70, best reproduce the experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Nucleon-induced inelastic scattering with statistical strength functions and the ECIS direct reaction code

Modern theoretical descriptions of inelastic scattering make use of multi-step direct reaction approaches together with transition potentials obtained from sophisticated nuclear structure models. Here we demonstrate how the complexity of such calculations can be reduced to permit simpler ones, also using the ECIS code, but providing an almost equally precise alternative to a much more detailed calculation. We have studied the transition form factors within the random phase approximation (RPA), where these are obtained as linear combinations of particle–hole states. At moderate to high excitation energies, where interference effects tend to disappear, we have proposed an independent particle–hole formalism in which particle–hole states are spread in energy with an appropriate strength function obtained from the RPA. The effects of more complex modes, such as 2p–2h ones, are simulated with widths calculated in a semi-classical context. In this work, we verify the validity of our approximations for pre-equilibrium proton-induced reactions on 90 Zr target. Our calculations provide a good description of the reaction data and point toward a simplification of the description of nucleon-induced reactions based on averages of microscopic details of the projectile–target interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Impact of HF-CRPA CCQE Model on the Latest NOvA Results

NOvA is a long-baseline neutrino experiment based at Fermilab, dedicated to measuring various neutrino oscillation parameters. Recently, the NOvA collaboration presented new results at NEUTRINO 2024. A significant enhancement in the modeling of charged-current quasi elastic (CCQE) interactions has been achieved with the implementation of the Hartree-Fock (HF) mean-field model, incorporating continuum random phase approximation (CRPA) corrections. This HF-CRPA model offers substantial improvements in the low-energy region. In this presentation, I will discuss the impact of the HF-CRPA model on the latest three-flavor oscillation results of the NOvA experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Modeling 2p-2h Interactions in Neutrino-Nucleus Scattering: A Comparative Overview of the Valencia and SuSAv2-MEC Models

2p-2h interactions are crucial for describing ν-A scattering, especially in the "dip region" between the Quasi-Elastic (QE) peak and Δ-resonance production. They involve the ejection of two nucleons, leading to two holes in the nuclear ground state. The underlying nuclear dynamics involve short-range correlations (SRC), long-range correlations (Random Phase Approximation, RPA), and the interplay of one- and two-body currents. Several models attempt to describe 2p-2h interactions, including Microscopic Models (e.g., Valencia), Scaling-Based Models (e.g., SuSAv2-MEC), and Transport Models (e.g., GiBUU). Despite progress, discrepancies exist between models and with data. I Will give an overview of 2p2h from the perspective of two Models Valencia and SuSAv2-MEC.

Hassinin, Karim [Houston U.]↗

DOE OSTI2025

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB+QRPA calculations with the Gogny interaction

The Zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98Zr. This region includes a benchmarking case and offers insights into nuclear deformation phenomena. To investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially-symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the Quasiparticle Random Phase Approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the K π = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the K π = 0 − and K π = 1 − components. Intrinsic transition densities and associ ated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e. we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

MEC modeling

2p-2h interactions are crucial for describing ν-Ar scattering, especially in the "dip region" between the Quasi-Elastic (QE) peak and Δ-resonance production. They involve the ejection of two nucleons, leading to two holes in the nuclear ground state. The underlying nuclear dynamics involve short-range correlations (SRC), long-range correlations (Random Phase Approximation, RPA), and the interplay of one- and two-body currents. Several models attempt to describe 2p-2h interactions like Valencia and SuSAv2-MEC. Despite progress, discrepancies exist between models and with data. I Will give an overview of 2p2h from the perspective of two Models Valencia and SuSAv2-MEC. I will introduce my current work to apply reweighting from SuSAv2 to Valencia on argon which will help improving the systematic uncertainty on MEC for SBN and DUNE.

Hassinin, Karim [Houston U.] (ORCID:00090009778353↗

DOE OSTI2022

Transverse Enhancement, Longitudinal Quenching and Coulomb Sum Rule in e-$^{12}$C and e-$^{16}$O Scattering

We present a short summary of a phenomenological analysis of all available electron scattering data on 12 C 12 C (about 6600 differential cross section measurements) and on 16 O 16 O (about 250 measurements) within the framework of the quasielastic superscaling model (including Pauli blocking). All quasielastic and inelastic cross section measurements are included down to the lowest momentum transfer ? q (including photo-production data). We find that there is enhancement of the transverse quasielastic response function and quenching of the longitudinal response function at low ? q (in addition to Pauli blocking). We extract parameterizations of a ? ? ? ? ? ? ? ? ? ? ? ? ? ? multiplicative low ? q "Longitudinal Quenching Factor" and an ? ? ? ? ? ? ? ? additive "Transverse Enhancement" contribution. Additionally, we find that the excitation of nuclear states contribute significantly (up to 30\%) to the Coulomb Sum Rule. We extract the most accurate determination of the sum rule to date and find it to be in disagreement with random phase approximation based calculations but in reasonable agreement with recent theoretical calculations such as "First Principle Green's Function Monte Carlo".

Bodek, Arie↗

DOE OSTI2023

Attosecond Time Delay Trends across the Isoelectronic Noble Gas Sequence

The analysis and measurement of Wigner time delays can provide detailed information about the electronic environment within and around atomic and molecular systems, with one the key differences being the lack of a long-range potential after a halogen ion undergoes photoionization. In this work, we use relativistic random-phase approximation to calculate the average Wigner delay from the highest occupied subshells of the atomic pairings (2p, 2s in Fluorine, Neon), (3p, 3s in Chlorine, Argon), (4p, 4s, 3d, in Bromine, Krypton), and (5p, 5s, 4d in Iodine, Xenon). The qualitative behaviors of the Wigner delays between the isoelectronic pairings were found to be similar in nature, with the only large differences occurring at photoelectron energies less than 20 eV and around Cooper minima. Interestingly, the relative shift in Wigner time delays between negatively charged halogens and noble gases decreases as atomic mass increases. All atomic pairings show large differences at low energies, with noble gas atoms showing large positive Wigner delays, while negatively charged halogen ions show negative delays. The implications for photoionization studies in halide-containing molecules is also discussed.

Grafstrom, Brock↗

DOE OSTI2023

Cross Sections for Coherent Elastic and Inelastic Neutrino-Nucleus Scattering

The prospects of extracting new physics signals in coherent elastic neutrino–nucleus scattering (CEνNS) processes are limited by the precision with which the underlying nuclear structure physics, embedded in the weak nuclear form factor, is known. We present calculations of charge and weak nuclear form factors and CEνNS cross sections on 12C, 16O, 40Ar, 56Fe and 208Pb nuclei. We obtain the proton and neutron densities, and charge and weak form factors by solving Hartree–Fock (HF) equations with a Skyrme (SkE2) nuclear potential. We validate our approach by comparing 208Pb and 40Ar charge form factor predictions with available elastic electron scattering data. Since CEνNS experiments at stopped-pion sources are also well suited to measure inelastic charged–current and neutral–current neutrino–nucleus cross sections, we also present calculations for these processes, incorporating a continuum Random Phase Approximation (CRPA) description on top of the HF–SkE2 picture of the nucleus. Providing both coherent as well as inelastic cross sections in a consistent framework, we aim at obtaining a reliable and detailed comparison of the strength of these processes in the energy region below 100 MeV. Furthermore, we attempt to gauge the level of theoretical uncertainty pertaining to the description of the 40Ar form factor and CEνNS cross sections by comparing relative differences between recent microscopic nuclear theory and widely-used phenomenological form factor predictions. Future precision measurements of CEνNS will potentially help in constraining these nuclear structure details that will in turn improve prospects of extracting new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Comprehensive Analysis of the Neutrino Process in Core-collapsing Supernovae

We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Continuum contribution to charged-current absorption of low-energy $ν_e$ on $^{40}$Ar

Accurate modeling of the absorption of tens-of-MeV $ν_e$ on $^{40}$Ar is needed to enable measurements of astrophysical neutrinos using large liquid argon time projection chamber (LArTPC) detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE). We revisit the MARLEY neutrino interaction model used in present estimates of DUNE sensitivity to supernova and solar neutrino signals. Multiple theoretical refinements are pursued, especially in the unbound continuum region of nuclear excitation energy. Inclusive charged-current neutrino-argon cross sections are calculated using a hybrid strategy. Nuclear transitions to unbound states are treated using a Hartree-Fock Continuum Random Phase Approximation (HF-CRPA) model, including forbidden contributions. Allowed transitions to low-lying discrete levels are also included using indirect measurements and approximate corrections for the momentum transfer dependence. Exclusive predictions are obtained by coupling these calculations with a statistical nuclear de-excitation model. The impact on observables of interest for DUNE and similar experiments is examined in terms of both total and differential cross sections. Our refined calculations predict a lower allowed portion of the cross section relative to the prior MARLEY model. At neutrino energies appreciably below 100 MeV, the inclusion of forbidden transitions does not fully compensate for the loss of allowed strength. For a representative neutrino burst from a galactic core-collapse supernova, our results suggest that MARLEY 1.2.0 overestimates the event yield in a DUNE-like detector by approximately 20%. However, because this overestimation is more severe at backwards angles, use of the charged-current $ν_e$-$^{40}$Ar reaction for supernova pointing may be more feasible than previously expected.

Gardiner, Steven [Fermilab]↗

DOE OSTI2023

Competing magnetic fluctuations and orders in a multiorbital model of doped SrCo 2 ⁢As 2

We revisit the intriguing magnetic behavior of the paradigmatic itinerant frustrated magnet Sr⁢Co 2 ⁢As 2 , which shows strong and competing magnetic fluctuations yet does not develop long-range magnetic order. By calculating the static spin susceptibility Χ⁡(q) within a realistic 16-orbital Hubbard-Hund model, we determine the leading instability to be ferromagnetic (FM). We then explore the effect of doping and calculate the critical Hubbard interaction strength U c that is required for the development of magnetic order. We find that U c decreases under electron doping and with increasing Hund's coupling J, but increases rapidly under hole doping. This suggests that magnetic order could possibly emerge under electron doping but not under hole doping, which agrees with experimental findings. We map out the leading magnetic instability as a function of doping and Hund's coupling and find several antiferromagnetic phases in addition to FM. We also quantify the degree of itinerant frustration in the model and resolve the contributions of different orbitals to the magnetic susceptibility. Lastly, we discuss the dynamic spin susceptibility Χ⁡(q,ω) at finite frequencies, where we recover the anisotropy of the peaks at Q π = (π,0) and (0,π) observed by inelastic neutron scattering that is associated with the phenomenon of itinerant magnetic frustration. By comparing results between theory and experiment, we conclude that the essential experimental features of doped SrCo 2 ⁢As 2 are well captured by an itinerant Hubbard-Hund multiorbital model if one considers a small shift of the chemical potential towards hole doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

How well do one-electron self-interaction-correction methods perform for systems with fractional electrons?

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (z σ ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew–Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of $H$$^{+}_{2}$ and $H$$^{+}_{2}$ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li + to F – at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗