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 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↗
DOE OSTI2024

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗
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↗
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↗
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↗
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↗