Engineering Papers⌕ Search

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.

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 OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2024

Electronic structure, self-doping, and superconducting instability in the alternating single-layer trilayer stacking nickelates La 3 ⁢Ni 2 ⁢O 7

Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La 3 ⁢Ni 2 ⁢O 7 , we comprehensively study this system using ab initio and random-phase approximation techniques. Here, our analysis unveils similarities between this novel La 3⁢ Ni 2 ⁢O 7 structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap value and orbital-selective behavior of the eg orbitals. Pressure primarily increases the bandwidths of the Ni e g bands, suggesting an enhancement of the itinerant properties of those e g states. By changing the cell volume ratio V/V 0 from 0.9 to 1.10, we found that the bilayer structure in La 3 ⁢Ni 2 ⁢O 7 always has lower energy than the single-layer trilayer stacking La 3 ⁢Ni 2 ⁢O 7 . In addition, we observe a “self-doping” effect (compared to the average 1.5 electrons per eg orbital per site of the entire structure) from the trilayer to the single-layer sublattices and this effect will be enhanced by overall electron doping. Moreover, we find a leading d x 2 -y 2 -wave pairing state that is restricted to the single layer. Because the effective coupling between the single layers is very weak, due to the nonsuperconducting trilayer in-between, this suggests that the superconducting transition temperature T c in this structure should be much lower than in the bilayer structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2024

Prediction of s ± -Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La 4 Ni 3 O 10 under Pressure

Motivated by the recently reported signatures of superconductivity in trilayer La 4 ⁢Ni 3 ⁢O 10 under pressure, here we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d 3z 2 –r 2 orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p z orbitals inducing a “trimer” lattice in La 4⁢ Ni 3 ⁢O 10 , analogous to the dimers of La 3 ⁢Ni 2 ⁢O 7 . The Fermi surface consists of three electron sheets with mixed e g orbitals, and a hole and an electron pocket made up of the d 3⁢z 2 –r 2 orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s ± -wave channel due to partial nesting between the M = (π,π) centered pockets and portions of the Fermi surface centered at the Γ = (0,0) point. With changing electronic density n, the s ± instability remains leading and its pairing strength shows a domelike behavior with a maximum around n = 4.2 ( ~6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La 3 ⁢Ni 2 ⁢O 7 , correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T c superconductivity of La 3⁢ Ni 2 ⁢O 7 does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La 4 ⁢Ni 3 ⁢O 10 with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗