Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
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 OSTI2025

Nuclear β − -decay with statistical de-excitation

he accurate description of nuclear β − -decay has far-reaching consequences for applications spanning nuclear reactors to the creation of heavy elements in astrophysical environments. We present the nuclear particle spectra associated with the β -decay of neutron-rich nuclei calculated with the well benchmarked coupled Quasi-particle Random Phase Approximation and Hauser–Feshbach (QRPA+HF) model. This approach begins with the population of the daughter nucleus via semi-microscopic Gamow-Teller or First-Forbidden strength distributions (QRPA) and follows the statistical de-excitation (HF) until the initial available excitation energy is exhausted. At each stage of de-excitation the emission by neutrons and $γ$-rays is considered obeying quantum mechanical selection rules. For completeness we also provide parsed Auger and Internal Conversion (IC) electron spectra from Evaluated Nuclear Data Files (ENDF). Our results are tabulated and provided in parsable ASCII formatted tables that are suitable for inclusion in various applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Investigation of the isoscalar monopole response in the proton-rich nucleus 14 O

Deuteron inelastic scattering on 14 O was measured in inverse kinematics using an active-target time projection chamber and a magnetic spectrograph. The experimental technique enabled precise measurements of deuteron recoiling particles in coincidence with beam-like fragments detected in the spectrograph focal plane. The double differential cross section was reconstructed for scattering angles of 3–6 degrees and excitation energies up to 26 MeV. The monopole strength distribution was obtained from the data using a multipole decomposition analysis. The results were compared to quasiparticle random-phase approximation (QRPA) and generator coordinate method (GCM) calculations. The QRPA calculation accurately describes experimental data in the energy range of 13 to 26 MeV. GCM calculations assuming a 12 C (g . s .) + p + p cluster configuration were used to determine the 0 + strength in 14 O below 13 MeV. The monopole transition strength of these cluster states provides a good description of the experimental distribution in the 9–11 MeV region, while the 0$^{+}_{2}$ state accounts for only a small fraction of the experimental strength around 6 MeV.

Active target↗
DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

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

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

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

The Influence of Charge Correlation and Ion Solvation on the Phase Behavior of Single-Ion Conducting Polymer Blend Electrolytes Using SAXS/SANS

Single-ion conducting polymer blends (SICPBs) have demonstrated exceptional electrochemical performance as solid-state battery electrolytes; however, their nanoscale morphology and thermodynamic behavior remain unexplored. In this work, we investigate blends composed of deuterated poly(ethylene oxide) and poly[lithium sulfonyl(trifluoromethane sulfonyl)imide methacrylate], dPEO/P(LiMTFSI), and report the first experimental study of the nanostructures of charge-neutral polymer blends using small-angle neutron scattering (SANS) and small-angle X-ray scattering (SAXS). Despite the macroscopic miscibility indicated by a single glass-transition temperature, SANS and SAXS results reveal disordered, charge-correlated nanostructures that are strongly influenced by blend composition and temperature. At low concentrations of charge polymer, the scattering is dominated by concentration fluctuations, and the random phase approximation is applied to extract values of the Flory–Huggins interaction parameter, χ SC . At higher charged polymer content, concentration fluctuations are suppressed, and a correlation model is used to characterize the nanostructures of the charge correlations. We find that the structures of the charge correlations are highly dependent on blend composition─consistent with predictions from Sing’s self-consistent field theory-liquid state models. Understanding these features is essential for uncovering the ion transport mechanism that leads to improved electrochemical performance previously reported in SICPB systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
DOE OSTI2026

Mechanisms of Metal Additive-Induced Ordering During SNIPS Membrane Formation

Isoporous membranes can be fabricated by combining self-assembly with nonsolvent induced phase separation (SNIPS) using an amphiphilic block copolymer like polystyrene-b-poly(4-vinylpyridine) (SV). Poly(4-vinylpyridine) (V) is known to complex with metal salts, which are hypothesized to stabilize solution ordering and preserve structure during casting. We explored how the molar ratio of metal additive to the poly(4-vinylpyridine) block affected the final membrane morphology via scanning electron microscopy (SEM). Dynamic light scattering (DLS), small-angle X-ray scattering (SAXS), and in situ grazing-incidence SAXS were used to track changes in solution ordering and chain conformation as a function of the molar ratio of the additive to the V block. Additives induced aggregation, promoted the formation of more compact conformations in solution, and facilitated micelle ordering onto lattices at optimal ratios. Furthermore, these experimental results were supported by random phase approximation calculations, which helped explain how the thermodynamic order–disorder transition shifts with additive binding strength. Stronger additive–polymer interactions reduced the block copolymer volume fraction required for ordering in solution, allowing ordered domains to form at lower polymer concentrations.

Additives↗
DOE OSTI2025

Microscopic calculations with noniterative finite amplitude methods and the application to neutron radiative captures and inelastic scatterings

We derive the fully self-consistent quasiparticle random-phase approximation (QRPA) equations with noniterative finite amplitude methods and calculate the transition strengths of giant resonances. Then, we apply the QRPA results to both neutron radiative capture calculations based on the statistical Hauser-Feshbach theory and inelastic scattering calculations based on distorted-wave Born approximation (DWBA). We compare the calculated results with available experimental data and demonstrate how our approach can reproduce giant resonances and various nuclear reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electric dipole excitations near the neutron separation energies in 96 Mo

Electric dipole strength near the neutron separation energy significantly impacts nuclear structure properties and astrophysical scenarios. These excitations are complex in nature and may involve the so-called pygmy dipole resonance (PDR). Transition densities play a crucial role in understanding the nature of nuclear excited states, including collective excitations, as well as in constructing transition potentials in DWBA or coupled-channels equations. In this work, we focus on electric dipole excitations in spherical molybdenum isotopes, particularly 96 Mo, employing fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) methods. We analyze the dipole strength near the neutron separation energy, which represents the threshold for neutron capture processes, and examine the isospin characteristics of PDR states through transition density calculations. Examination of proton and neutron transition densities reveals distinctive features of each dipole state, indicating their isoscalar and isovector nature. We observe that the primary component in the enhanced low-energy region exhibits isovector character. The PDR displays a mixture of isoscalar and isovector nature, distinguishing it from the isovector giant dipole resonance (IVGDR). These findings lay the groundwork for future investigations into the role of transition densities in reaction models and for their application to inelastic scattering calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electronic structure, and magnetic and superconducting pairing tendencies of the alternating single layer–bilayer stacking nickelate La 5 ⁢Ni 3 ⁢O 11 under pressure

Nickelates have continued to surprise since their unconventional superconductivity was discovered. Recently, the layered nickelate La 5 ⁢Ni 3 ⁢O 11 with hybrid single-layer and bilayer stacking showed superconductivity under high pressure. This compound combines features of single-layer La 2 ⁢NiO 4 and bilayer La 3 ⁢Ni 2 ⁢O 7 , but its pairing mechanism remains to be understood. Motivated by this finding, here we report a comprehensive theoretical study of this system. Our density functional theory calculations reveal that the undistorted P4/mmm phase without pressure is unstable due to three distortion modes. As pressure increases, these modes are suppressed, leading to the high-symmetry P4/mmm phase without NiO 6 octahedron tilting. Moreover, the “charge transfer” between the single-layer and bilayer sublattices was observed, leading to hole doping in the single-layer blocks. Our random-phase approximation calculations indicate a leading 𝑑 𝑥 2 −𝑦 2 -wave pairing state that arises from spin-fluctuation scattering between Fermi surface states mainly originating from the single-layer blocks and additional weaker contributions from the bilayer blocks. These spin fluctuations could be detected by inelastic neutron scattering as a strong peak at 𝐪 = (𝜋,𝜋). Our findings distinguish La 5 ⁢Ni 3 ⁢O 11 from other nickelate superconductors discovered so far and the high-𝑇 𝑐 cuprates. We also discuss both similarities and differences between La 5⁢ Ni 3 ⁢O 11 and other hybrid stacking nickelates.

Zhang, Yang [Univ. of Tennessee, Knoxville, TN (Un↗
DOE OSTI2026

Beta-Decay Half-Lives beyond 54 Ca: A Systematic Survey of Decay Properties Approaching the Neutron Dripline

In an experiment performed at the Facility for Rare Isotope Beams (FRIB) using the FRIB Decay Station initiator, 15 new half-lives of isotopes near 54 Ca were measured. A new method of extracting lifetimes from experimental data, taking into account the unknown 𝛽-delayed neutron emission branches of very neutron-rich nuclei, was developed to enable systematic uncertainty analysis. The experiment observed a dramatic change in the half-life systematics for the isotopes with neutron number 𝑁 = 34. Beyond 𝑁 = 34, the decline of nuclear lifetime is much slower, leading to longer than anticipated lifetimes for near-dripline nuclei. State-of-the-art shell-model calculations can explain the experimental results, revealing the imprint of shell effects and the need for modification of single-particle neutron states. In conclusion, the results from a newly developed quasiparticle random phase approximation model with potential for making global predictions were also tested against the experimental results and good agreement was found.

Physics - Nuclear physics and radiation physics↗
DOE OSTI2026

Compressive strain turns 𝑠 ± - into 𝑑-wave pairing in a one-unit-cell La 3 ⁢Ni 2 ⁢O 7 thin film via substrate-induced hole doping

Motivated by recent reports of ambient-pressure superconductivity in La 3 ⁢Ni 2⁢ O 7 films grown on LaSrAlO 4 , we investigate the superconducting instability in a one-unit-cell (1UC) thin film using ab initio and random-phase approximation techniques. Compared to the high-pressure bulk system, the ratio of interlayer 𝑑 3⁢𝑧 2 −𝑟 2 hopping to intralayer 𝑑 𝑥 2 −𝑦 2 hopping is suppressed in the 1UC thin film, and the crystal-field splitting of the 𝑒 𝑔 orbitals is increased. Here, our calculation indicates that spin-fluctuation-driven pairing correlations are weak for the stoichiometric case at ambient pressure, but increase significantly under hole doping. The leading pairing symmetry is also found to change by hole doping. Specifically, we obtain a leading 𝑑 𝑥 2 −𝑦 2 pairing state at moderate hole doping, followed by a 𝑑 𝑥⁢𝑦 state at higher doping. These states are driven by intraband spin-fluctuation scattering within the 𝛾 hole pocket centered around the 𝑀 point, and arise primarily from states in the Ni layer farther from the substrate. These results strongly suggest that the thin-film superconducting samples are hole-doped and that pairing in this system predominantly arises in the layer, as opposed to the interlayer pairing in the pressurized bulk system.

Zhang, Yang [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2025

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

We compute moments of the isoscalar monopole response of 𝑁 = 𝑍 closed-shell nuclei based on chiral nucleon-nucleon plus three-nucleon interactions. We employ the random-phase approximation (RPA) and two ab initio many-body approaches, the in-medium similarity renormalization group (IMSRG) and coupled-cluster theory (CC). In the IMSRG framework, the moments are obtained as ground-state expectation values, whereas in the CC approach, they are evaluated through excited-state calculations. We find good agreement between the IMSRG and CC results across all nuclei studied. RPA provides a reasonable approximation to the correlated methods if the interaction is soft. From the calculated moments, we extract average energies of the monopole response, compute finite-nucleus incompressibilities, and estimate the incompressibility of symmetric nuclear matter by a fit to a leptodermous expansion. Our extrapolated values are lower than those obtained in nuclear-matter calculations with the same interactions, but the values are consistent with phenomenological ranges.

Bonaiti, Francesca [Michigan State Univ., East Lan↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗