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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1975

Photoionization cross sections for atomic chlorine using an open-shell random phase approximation

The use of the Random Phase Approximation with Exchange (RPAE) for calculating partial and total photoionization cross sections and photoelectron angular distributions for open shell atoms is examined for atomic chlorine. Whereas the RPAE corrections in argon (Z=18) are large, it is found that those in chlorine (Z=17) are much smaller due to geometric factors. Hartree-Fock calculations with and without core relaxation are also presented. Sizable deviations from the close coupling results of Conneely are also found.

Starace, A. F.↗
DOE OSTI2025

Embedded random phase approximation for magnetic systems: H 2 dissociative adsorption on Fe(110)

The random phase approximation (RPA), a method for treating electron correlation, has been shown to be superior to standard density functional theory (DFT) approximations in numerous cases. However, the RPA’s computational cost is substantially higher than that of DFT, particularly restricting its application to extended surfaces. The recently introduced embedded RPA (emb-RPA) approach [Wei et al., J. Chem. Phys. 159(19), 194108 (2023)] reduces this computational cost by approximately two orders of magnitude. While previous applications of emb-RPA focused on non-spin-polarized systems, here we extend the approach to ferromagnetic ones. Unlike other embedded correlated wavefunction methods, such as embedded complete active space self-consistent field theory, emb-RPA is advantageous for spin-polarized systems because the RPA is compatible with unrestricted DFT solutions, which are eigenfunctions of the spin angular momentum operator S z but not the total spin-squared operator S 2 . By applying emb-RPA with specific magnetization constraints, we achieved a speedup of two to three orders of magnitude (one order when accounting for the one-time embedding potential optimization cost) with only small errors (∼50 meV) compared to full periodic RPA. Moreover, emb-RPA significantly reduces the over-binding errors of DFT approximations. In conclusion, we anticipate that the acceleration enabled by the spin-polarized emb-RPA approach will broaden the applicability of RPA to magnetic materials.

Density functional theory↗
NASA NTRS1974

Partition function for a two dimensional plasma in the random phase approximation

The partition function for a two-dimensional plasma is evaluated within the random phase approximation. The periodic boundary conditions are fully taken into account by including the periodic image interactions. In the guiding-center limit, the negative temperature threshold energy is evaluated, and a value different from previous calculations results. When an identical random phase evaluated, and a value different from previous calculations results. When an identical random phase evaluation is applied to the finite gyroradius plasma, the Salzberg-Prager-May equation of state is recovered.

Seyler, C. E., Jr.↗
DOE OSTI2026

Optimization of Random Phase Approximation Calculations for Improved Energies of Molecules, Solids, and Surfaces

We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.

Adsorption↗
NASA NTRS1974

Partition function for a two-dimensional plasma in the random-phase approximation

The partition function for a two-dimensional plasma is evaluated within the random-phase approximation. The periodic boundary conditions are taken fully into account by including the periodic image interactions. In the guiding-center limit, the 'negative temperature' threshold energy is evaluated, and a value different from previous calculations results. When an identical random-phase evaluation is applied to the finite-gyroradius plasma, the Salzberg-Prager-May equation of state is recovered.

Seyler, C. E., Jr.↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

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

Analysis of Small-Angle Neutron Scattering from Blends of Charged and Neutral Polymers Based on Rod–Coil Random Phase Approximation

Blends of charged and neutral polymers are of interest due to potential applications in rechargeable batteries. In this study, concentration fluctuations in blends of charged poly[lithium 3-(methacryloyloxy)propylsulfonyl-1-(trifluoromethanesulfonyl)imide] (PLiMTFSI) and neutral poly(ethylene oxide) (PEO) were investigated by small-angle neutron scattering (SANS). The scattering data were analyzed in the framework of the random phase approximation (RPA). Since ion dissociation can lead to stiffening, the charged polymers were approximated as rods, while the neutral polymers were assumed to be random coils. This approach works reasonably well at low weight fractions of charged polymers. For blends with higher weight fractions of the charged polymer, concentration fluctuations were highly suppressed, resulting in q-independent coherent structure factors that are inconsistent with the rod-coil RPA.

Lee, Jaeyong↗
DOE OSTI2025

General trends of superconducting pairing and magnetic correlations in the Ruddlesden-Popper nickelate 𝑚-layered superconductors La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1

Here, we report a comprehensive theoretical analysis of the Ruddlesden-Popper layered nickelates La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1 (𝑚 = 1 to 6) under pressure. These materials have recently received significant attention due to the discovery of superconductivity in some nickelates under pressure. Our results suggest that, while these Ruddlesden-Popper layered nickelates display many similarities, they also show noticeable differences. One of the common features of La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1 is that the electronic states near the Fermi level are mainly contributed by Ni 3⁢𝑑 orbitals, slightly hybridized with O 2⁢𝑝 orbitals. The Ni 𝑑 3⁢𝑧 2 −𝑟 2 orbitals display bonding-antibonding, or bonding-antibonding-nonbonding, characteristic splittings, depending on the even or odd number of stacking layers 𝑚. In addition, the ratio of the in-plane interorbital hopping between 𝑑 3⁢𝑧 2 −𝑟 2 and 𝑑 𝑥 2 −𝑦 2 orbitals and in-plane intraorbital hopping between 𝑑 𝑥 2 −𝑦 2 orbitals was found to be large in La 𝑚+1 ⁢Ni 𝑚 ⁢O 3⁢𝑚+1 (𝑚 = 1 to 6), and this ratio increases from 𝑚 = 1 to 𝑚 = 6, suggesting that the in-plane hybridization will increase as the layer number 𝑚 increases. In contrast to the dominant 𝑠 ± -wave state driven by spin fluctuations in the bilayer La 3 ⁢Ni 2⁢ O 7 and trilayer La 4 ⁢Ni 3 ⁢O 10 , two nearly degenerate 𝑑 𝑥 2 −𝑦 2 -wave and 𝑠 ± -wave leading states were obtained in the four-layer stacking La 5⁢ Ni 4 ⁢O 13 and five-layer stacking La 6 ⁢Ni 5 ⁢O 16 . The leading 𝑠 ± -wave state was recovered in the six-layer material La 7 ⁢Ni 6 ⁢O 19 with slightly higher calculated pairing strength 𝜆 than that of the 𝑑 𝑥 2 −𝑦 2 -wave state. All this evidence suggests that both 𝑠 ± -wave and 𝑑 𝑥 2 −𝑦 2 -wave channels are strongly competing in the high-order niceklates based on our random-phase approximation calculations. In general, at the level of the random-phase approximation treatment, the superconducting transition temperature 𝑇 𝑐 decreases in stoichiometric bulk systems from the bilayer La 3 ⁢Ni 2 ⁢O 7 to the six-layer La 7 ⁢Ni 6 ⁢O 19 , despite the 𝑚-dependent dominant pairing. Both in-plane and out-of-plane magnetic correlations are found to be quite complex. Within the in-plane direction, we obtained the peak of the magnetic susceptibility at 𝐪 = (0.6⁢𝜋, 0.6⁢𝜋) for La 5 ⁢Ni 4 ⁢O 13 (𝑚 = 4) and La 7 ⁢Ni 6 ⁢O 19 (𝑚 = 6) and at 𝐪 = (0.7⁢𝜋, 0.7⁢𝜋) for La 6⁢ Ni 5⁢ O 16 (𝑚 = 5). Along the out-of-plane direction, four layers are coupled as ↓−↑−↑−↓ in La 5 ⁢Ni 4⁢ O 13 , five layers are coupled as ↑−↑−↓−↑−↑ in La 6 ⁢Ni 5 ⁢O 16 , and six layers are coupled as ↑−↓−↓−↑−↑−↓ in La 7 ⁢Ni 6 ⁢O 19 .

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

Magnetic correlations and pairing tendencies of the hybrid stacking nickelate superlattice La 7 Ni 5 O 17 (La 3 ⁢Ni 2 ⁢O 7 /La 4 Ni 3 ⁢O 10 ) under pressure

Motivated by the recent rapid progress in high-𝑇 𝑐 nickelate superconductors, we comprehensively study the physical properties of the alternating bilayer trilayer stacking nickelate La 7 ⁢Ni 5 ⁢O 17 . The high-symmetry phase of this material, without the tilting of oxygen octahedra, is not stable at ambient conditions but becomes stable under high pressure, where a small hole pocket 𝛾 0 , composed of the 𝑑 3⁢𝑧 2 −𝑟 2 states in the trilayer sublattice, appears. Here, this pocket was identified in our previous work for trilayer La 4 ⁢Ni 3 ⁢O 10 as important to develop superconductivity. Moreover, using random-phase approximation calculations, we find a leading 𝑠 ± pairing state for the high-symmetry phase under pressure with similar pairing strength as that obtained previously for the bilayer La 3⁢ Ni 2 ⁢O 7 compound, suggesting a similar or higher superconducting transition temperature 𝑇 𝑐 , at the random-phase approximation level. In addition, we find that the dominant magnetic fluctuations in the system driving this pairing state have antiferromagnetic structure both in-plane and between the planes of the top and bottom trilayer and bilayer sublattices, while the middle trilayer is magnetically decoupled.

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

Magnetic phase diagram of a two-orbital model for bilayer nickelates with varying doping

Motivated by the recently discovered high-T c bilayer nickelate superconductor La 3⁢ Ni 2 ⁢O 7 , we comprehensively research a bilayer 2×2×2 cluster for different electronic densities n by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying n. Based on the spin structure factor S(q), we have obtained a rich magnetic phase diagram in the plane defined by n and U/W, at fixed Hund coupling, where U is the Hubbard strength and W the bandwidth. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. At half-filling, n=2 (two electrons per Ni site, corresponding to N=16 electrons), the canonical superexchange interaction leads to a robust G-AFM state (π,π,π) with antiferromagnetic couplings both in-plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the “half-empty” and “half-full” mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At n=1.5 (or N=12), density corresponding to La 3 ⁢Ni 2 ⁢O 7 , we obtained the “Stripe 2” ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the z axis) in the 2×2×2 cluster. In addition, we obtained a much stronger AFM coupling along the z axis than the magnetic coupling in the xy plane. The random-phase approximation calculations with varying n give very similar results as Lanczos, even though both techniques are based on quite different procedures. Additionally, a state with q/π=(0.6,0.6,1) close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to n=1.25, possibly responsible for the E-phase SDW recently observed in experiments. In conclusion, our predictions can be tested by chemically doping La 3 ⁢Ni 2 ⁢O 7 .

36 MATERIALS SCIENCE↗
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↗