Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical model”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Avoiding biases in binned fits

Abstract Binned maximum likelihood fits are an attractive option when analysing large datasets, but require care when computing likelihoods of continuous PDFs in bins.For many years the widely used statistical modelling package evaluated probabilities at the bin centre, leading to significant biases for strongly curved probability density functions.We demonstrate the biases with real-world examples, and introduce a PDF class to that removes these biases.The physics and computation performance of this new class are discussed.

Instruments & Instrumentation↗
DOE OSTI2025

Improving tropical cyclone rapid intensification forecasts with satellite measurements of sea surface salinity and calibrated machine learning

Forecasting rapid intensification (RI) of tropical cyclones (TC) is a mission known for large errors. One under-researched factor that affects TC intensification is salinity, which is important for density stratification in certain ocean regions and can affect the surface enthalpy flux under a strengthening hurricane. To investigate the impact and efficacy of using salinity information in state-of-the-art forecasting, we use a statistical model consisting of a variety of machine learning (ML) methods. For salinity data, we use satellite measurements of pre-storm sea surface salinity (SSS) as a proxy for the salinity stratification. We train and test the model on various ocean basins, including the Atlantic, eastern North Pacific and western North Pacific. A calibrator is trained on top of the ML models to correct and enhance probability forecasts. The calibrator significantly improves probability forecasts relative to recent works. The ML model performance is improved with the addition of SSS in the Eastern North Pacific, western North Pacific, and the Caribbean subregion of the North Atlantic, and the overall model performance is better than previous studies. SSS decreases model skill for a model trained on the full Atlantic basin. In the Indian Ocean, SSS is also notably correlated with RI occurrence, but the TC samples are not sufficient to train ML models.

hurricane↗
DOE OSTI2025

Predicting weather impacts on corn production in a data-limited region using a transfer learning approach

The stability of food supply and prices may depend more on annual changes in yields from year-to-year variability in weather than on longer-term average changes from changing climatic conditions. However, the absence of high-quality data on crop yields at fine spatial resolutions in many regions of the world makes it challenging to statistically model their response to interannual variability in weather patterns. Therefore, there is a need for empirical methods that can project annual crop yield changes even in limited data regions. Here, we propose a transfer learning algorithm that uses high spatial resolution data from one region to project yields in another region with more limited data. The goal of our work is to understand what data types can be beneficial for transferring learning from a source region to a very different target region with more limited data. We utilize Long Short-Term Memory to develop a transfer learning model that is trained on historical county-level corn yield in the United States and predicts district-level corn yield variations in India. Even using smaller amounts of data in India, simulating a data-scarce region, we achieve an average root mean square error of 0.48 bu acre−1 in predicting interannual yield variations. Using Shapley values to interpret results, we explore the contribution of the different weather parameters to interannual yield variability and find a larger influence of precipitation-related variables. Our study demonstrates the usefulness of this method for transferring models of weather impacts on crop yields trained on a data-rich country to one with more limited data. It suggests the potential of applying the transfer learning model to mitigate the need for extensive raw data globally.

Vishwakarma, Srishti [ORNL] (ORCID:000000031674419↗
DOE OSTI2024

Trade can buffer climate-induced risks and volatilities in crop supply

Climate change is intensifying the frequency and severity of extreme events, posing challenges to food security. Corn, a staple crop for billions, is particularly vulnerable to heat stress, a primary driver of yield variability. While many studies have examined climate impact on average corn yields, little attention has been given to the climate impact on production volatility. This study investigates the future volatility and risks associated with global corn supply under climate change, evaluating the potential benefits of two key adaptation strategies: irrigation and market integration. A statistical model is employed to estimate corn yield response to heat stress and utilize NEX-GDDP-CMIP6 climate data to project future production volatility and risks of substantial yield losses. Three metrics are introduced to quantify these risks: Sigma (σ), the standard deviation of year-on-year yield change, which reflects overall yield volatility; Rho (ρ), the risk of substantial loss, defined as the probability of yield falling below a critical threshold; and Beta (β), a relative risk coefficient that captures the volatility of a region's corn production compared to the globally integrated market. The analysis reveals a concerning trend of increasing year-on-year yield volatility (σ) across most regions and climate models. This volatility increase is significant for key corn-producing regions like Brazil and the United States. While irrigated corn production exhibits a smaller rise in volatility, suggesting irrigation as a potential buffer against climate change impacts, it is not a sustainable option as it can cause groundwater depletion. On the other hand, global market integration reduces overall volatility and market risks significantly with less sustainability concerns. Furthermore, these findings highlight the importance of a multidimensional approach to adaptation in the food sector. While irrigation can benefit individual farmers, promoting global market integration offers a broader solution for fostering resilience and sustainability across the entire food system.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

SAIGE-GPU: accelerating genome- and phenome-wide association studies using GPUs

Genome-wide association studies (GWAS) at biobank scale are computationally intensive, especially for admixed populations requiring robust statistical models. SAIGE is a widely used method for generalized linear mixed-model GWAS but is limited by its CPU-based implementation, making phenome-wide association studies impractical for many research groups. We developed SAIGE-GPU, a GPU-accelerated version of SAIGE that replaces CPU-intensive matrix operations with GPU-optimized kernels. The core innovation is distributing genetic relationship matrix calculations across GPUs and communication layers. Applied to 2068 phenotypes from 635 969 participants in the Million Veteran Program, including diverse and admixed populations, SAIGE-GPU achieved a 5-fold speedup in mixed model fitting on supercomputing infrastructure and cloud platforms. We further optimized the variant association testing step through multi-core and multi-trait parallelization. Deployed on Google Cloud Platform and Azure, the method provided substantial cost and time savings. Source code and binaries are available for download at https://github.com/saigegit/SAIGE/tree/SAIGE-GPU-1.3.3. A code snapshot is archived at Zenodo for reproducibility (DOI: [10.5281/zenodo.17642591]). SAIGE-GPU is available in a containerized format for use across HPC and cloud environments and is implemented in R/C++ and runs on Linux systems.

Rodriguez, Alex [Argonne National Laboratory (ANL)↗
DOE OSTI2023

En echelon faults reactivated by wastewater disposal near Musreau Lake, Alberta

We use machine-learning and cross-correlation techniques to enhance earthquake detectability by two magnitude units for the earthquake sequence near Musreau Lake, Alberta, which is induced by wastewater disposal. This deep catalogue reveals a series of en echelon ~N–S oriented strike-slip faults that are favourably oriented for reactivation. These faults require only ~0.6 MPa overpressure for triggering to occur. Earthquake activity occurs in bursts, or episodes; episodes restricted to the largest fault tend to have earthquakes starting near the southern end (distant from injectors) and progressing northwards (towards the injectors). While most events are concentrated along these ~N–S oriented faults, we also delineate smaller faults. Together, these findings suggest pore pressure as the triggering mechanism, where a time-dependent increase in pore pressure likely caused these faults to progressively reawaken. Analysis of the ‘next record-breaking event’, a statistical model that forecasts the sequencing of earthquake magnitudes, suggests that the next largest event would be M L ~4.3. The seismically illuminated length of the largest fault indicates potential magnitudes as large as M w 5.3.

58 GEOSCIENCES↗
DOE OSTI2024

BRAKER3: Fully automated genome annotation using RNA-seq and protein evidence with GeneMark-ETP, AUGUSTUS, and TSEBRA

Gene prediction has remained an active area of bioinformatics research for a long time. Still, gene prediction in large eukaryotic genomes presents a challenge that must be addressed by new algorithms. The amount and significance of the evidence available from transcriptomes and proteomes vary across genomes, between genes, and even along a single gene. User-friendly and accurate annotation pipelines that can cope with such data heterogeneity are needed. The previously developed annotation pipelines BRAKER1 and BRAKER2 use RNA-seq or protein data, respectively, but not both. A further significant performance improvement integrating all three data types was made by the recently released GeneMark-ETP. We here present the BRAKER3 pipeline that builds on GeneMark-ETP and AUGUSTUS, and further improves accuracy using the TSEBRA combiner. BRAKER3 annotates protein-coding genes in eukaryotic genomes using both short-read RNA-seq and a large protein database, along with statistical models learned iteratively and specifically for the target genome. We benchmarked the new pipeline on genomes of 11 species under an assumed level of relatedness of the target species proteome to available proteomes. BRAKER3 outperforms BRAKER1 and BRAKER2. The average transcript-level F1-score is increased by about 20 percentage points on average, whereas the difference is most pronounced for species with large and complex genomes. BRAKER3 also outperforms other existing tools, MAKER2, Funannotate, and FINDER. The code of BRAKER3 is available on GitHub and as a ready-to-run Docker container for execution with Docker or Singularity. Overall, BRAKER3 is an accurate, easy-to-use tool for eukaryotic genome annotation.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2020

Strange hadron production in Au + Au collisions at s N N = 7.7 , 11.5, 19.6, 27, and 39 GeV

In this paper, we present STAR measurements of strange hadron (K$0\atop{s}$, Λ, $\bar{Λ}$, $Ξ$ - , $\bar{Ξ}$ + , Ω - , $\bar{Ω}$ + , and Φ ) production at mid-rapidity (|y|<0.5) in Au+Au collisions at √s NN = 7.7 – 39 GeV from the Beam Energy Scan Program at the Relativistic Heavy Ion Collider (RHIC). Transverse momentum spectra, averaged transverse mass, and the overall integrated yields of these strange hadrons are presented versus the centrality and collision energy. Antibaryon-to-baryon ratios(Λ/$\bar{Λ}$, $Ξ$ - / $\bar{Ξ}$ + , Ω - / $\bar{Ω}$ + ) are presented as well, and used to test a thermal statistical model and to extract the temperature normalized strangeness and baryon chemical potentials at hadronic freeze-out (μ B / T ch μ S / T ch ) in central collisions. Strange baryon-to-pion ratios are compared to various model predictions in central collisions for all energies. The nuclear modification factors (R CP ) and antibaryon-to-meson ratios as a function of transverse momentum are presented for all collision energies. The K$0\atop{s}$ R CP shows no suppression for p T up to 3.5 GeV/c at energies of 7.7 and 11.5 GeV. The $\bar{Λ}$/ K$0\atop{s}$ ratio also shows baryon-to-meson enhancement at intermediate p T (~2.5 GeV/c) in central collisions at energies above 19.6 GeV. In conclusion, both observations suggest that there is likely a change of the underlying strange quark dynamics at collision energies below 19.6 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Proton capture on 34 S in the astrophysical energy regime of O-Ne novae

Nuclear reaction sensitivity studies have shown that the final isotopic abundance of O-Ne nova nucleosynthesis is dependent on the 34 S(p,γ) 35 Cl reaction at astrophysical energies corresponding to peak nova burning temperatures of 0.1–0.4 GK. Isotopic ratios of the S, Cl, and Ar products are all used in various methods of cosmochemical analysis of presolar meteoritic grains. Due to the lack of direct experimental data, the 34 S+p reaction rate has been estimated using statistical modeling or information from indirect nucleon transfer experiments. In order to provide direct reaction information, here the resonance strengths of several low energy resonances, E c.m. = 272–495 keV, in the 34 S(p,γ) 35 Cl reaction were measured for the first time in inverse kinematics using the DRAGON recoil separator located at TRIUMF, Canada’s Particle Accelerator Centre in Vancouver.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

First direct measurement of 59 Cu(p,α) 56 Ni: A step towards constraining the Ni-Cu cycle in the cosmos

Reactions on proton-rich nuclides drive the nucleosynthesis in core collapse supernovae (CCSNe) and in x-ray bursts (XRBs). CCSNe eject the nucleosynthesis products to the interstellar medium and hence are a potential inventory of p nuclei, whereas in XRBs nucleosynthesis powers the light curves. In both astrophysical sites the Ni-Cu cycle, which features a competition between 59 Cu(p,α) 56 Ni and 59 Cu(p,γ) 60 Zn, could potentially halt the production of heavier elements. Here, we report the first direct measurement of 59 Cu(p,α) 56 Ni using a reaccelerated 59Cu beam and a cryogenic solid hydrogen target. Furthermore, our results show that the reaction proceeds predominantly to the ground state of 56 Ni, and the experimental rate has been found to be lower than Hauser-Feshbach based statistical model predictions. New results hints that the νp process could operate at higher temperatures than previously inferred and therefore remains a viable site for synthesizing the heavier elements

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measurement of the Ne 18 ( α , p ) Na 21 reaction with the ANASEN active-target detector system at E c . m . = 2.5 – 4 MeV

The 18 Ne(α,p) 21 Na reaction plays a significant role in Type-I X-ray bursts. It is a major path in the breakout from the hot-CNO cycles to the synthesis of heavier elements in the αp- and rp-processes. An experiment to determine the cross section of this reaction was performed with the ANASEN active-target detector system, determining the cross section at energies between 2.5 and 4 MeV in the center-of-mass frame. The measured cross sections for reactions populating the ground state in 21 Na are consistent with results obtained from the time-inverse reaction, but significantly lower than the previously published experimental data of direct measurements. The total cross sections are also compared with those derived from indirect methods and statistical-model calculations. Furthermore, this experiment establishes a new experimental data set on the excitation function of the 18 Ne(α,p) 21 Na reaction, revealing the significance of the excited states' contributions to the total reaction cross section and allowing us to separate the contribution of the (α,2p) reaction. The impact of the measured cross section on thermal reaction rates is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Cross sections of the 83 Rb⁢(p,γ)⁢ 84 Sr and 84 Kr⁢(p,γ)⁢ 85 Rb reactions at energies characteristic of the astrophysical γ process

We have measured the cross section of the Rb 83 ( p , γ ) Sr 84 radiative capture reaction in inverse kinematics using a radioactive beam of Rb 83 at incident energies of 2.4 and 2.7 A MeV. Prior to the radioactive beam measurement, the Kr 84 ( p , γ ) Rb 85 radiative capture reaction was measured in inverse kinematics using a stable beam of Kr 84 at an incident energy of 2.7 A MeV. The effective relative kinetic energies of these measurements lie within the relevant energy window for the γ process in supernovae. The central values of the measured partial cross sections of both reactions were found to be 0.17 – 0.42 times the predictions of statistical model calculations. Assuming the predicted cross section at other energies is reduced by the same factor leads to a slightly higher calculated abundance of the p nucleus Sr 84 , caused by the reduced rate of the Sr 84 ( γ , p ) Rb 83 reaction derived from the present measurement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Strangeness-correlations on the pseudocritical line in ( 2 + 1 )-flavor QCD

We present some lattice QCD results on first ( χ 1 i ) and second ( χ 2 i ) cumulants of and correlations ( χ 11 i j ) among net baryon-number ( B ), strangeness ( S ) and electric charge ( Q ) along the pseudocritical line [ T p c ( μ B ) ] in the temperature ( T )–baryon chemical potential ( μ B ) phase diagram of ( 2 + 1 )-flavor QCD. We point out that violations of sum rules among second order cumulants, which hold in the isospin symmetric limit of vanishing electric charge chemical potential, are small along the T p c ( μ B ) for the entire range of μ B covered in the RHIC beam energy scan. For the strangeness neutral matter produced in heavy-ion collisions this leads to a close relation between χ 11 B S and χ 11 Q S . We compare lattice QCD results for χ 11 B S / χ 2 S along the T p c ( μ B ) line with preliminary experimental measurements of χ 11 B S / χ 2 S for collision energies 7.7 GeV ≤ s N N ≤ 62.4 GeV . While we find good agreements for s N N ≥ 39 GeV , differences are sizeable at smaller values of s N N . Moreover, we compare lattice QCD results for the ratio of the strangeness ( μ S ) to baryon ( μ B ) chemical potentials, which define a strangeness neutral system with fixed electric charge to baryon number density, with experimental results obtained by the STAR collaboration for μ S / μ B using strange baryon yields on the freeze-out line. Finally, we determine the baryon chemical potential at the freeze-out ( μ B f ) by comparing χ 1 B / χ 2 B along the T p c ( μ B ) with the experimentally measured net-proton cumulants χ 1 p / χ 2 p . We find that { μ B f , T p c ( μ B f ) } are consistent with the freeze-out parameters of the statistical-model fits to experimentally measured hadron yields for s N N ≥ 11.5 GeV . Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

First Measurement of Antideuteron Number Fluctuations at Energies Available at the Large Hadron Collider

The first measurement of event-by-event antideuteron number fluctuations in high energy heavy-ion collisions is presented. The measurements are carried out at midrapidity (|η|<0.8) as a function of collision centrality in Pb-Pb collisions at $\sqrt{^sNN}$=5.02 TeV using the ALICE detector. A significant negative correlation between the produced antiprotons and antideuterons is observed in all collision centralities. The results are compared with a state-of-the-art coalescence calculation. While it describes the ratio of higher order cumulants of the antideuteron multiplicity distribution, it fails to describe quantitatively the magnitude of the correlation between antiproton and antideuteron production. On the other hand, thermal-statistical model calculations describe all the measured observables within uncertainties only for correlation volumes that are different with respect to those describing proton yields and a similar measurement of net-proton number fluctuations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

ψ ( 2 S ) Suppression in Pb-Pb Collisions at the LHC

The production of the ψ ( 2 S ) charmonium state was measured with ALICE in Pb-Pb collisions at s NN = 5.02 TeV , in the dimuon decay channel. A significant signal was observed for the first time at LHC energies down to zero transverse momentum, at forward rapidity ( 2.5 < y < 4 ). The measurement of the ratio of the inclusive production cross sections of the ψ ( 2 S ) and J / ψ resonances is reported as a function of the centrality of the collisions and of transverse momentum, in the region p T < 12 GeV / c . The results are compared with the corresponding measurements in p p collisions, by forming the double ratio [ σ ψ ( 2 S ) / σ J / ψ ] Pb − Pb / [ σ ψ ( 2 S ) / σ J / ψ ] p p . It is found that in Pb-Pb collisions the ψ ( 2 S ) is suppressed by a factor of ∼ 2 with respect to the J / ψ . The ψ ( 2 S ) nuclear modification factor R AA was also obtained as a function of both centrality and p T . The results show that the ψ ( 2 S ) resonance yield is strongly suppressed in Pb-Pb collisions, by a factor of up to ∼ 3 with respect to p p . Comparisons of cross section ratios with previous Super Proton Synchrotron findings by the NA50 experiment and of R AA with higher- p T results at LHC energy are also reported. These results and the corresponding comparisons with calculations of transport and statistical models address questions on the presence and properties of charmonium states in the quark-gluon plasma formed in nuclear collisions at the LHC. © 2024 CERN, for the ALICE Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Data-driven approach to parameterize SCAN + U for an accurate description of 3 d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the DFT + U method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. Here we develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-U correction, and (3) DFT/DFT + U compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for 3d transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve SCAN/SCAN + U compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize DFT + U performance for other target physical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Microscopic theory of angular momentum distributions across the full range of fission fragments

Modern nuclear theory provides qualitative insights into the fundamental mechanisms of nuclear fission and is increasingly capable of making reliable quantitative predictions. Most quantities of interest pertain to the primary fission fragments, whose subsequent decay is typically modeled using statistical reaction models. Consequently, a key objective of fission theory is to inform these models by predicting the initial conditions of the primary fragments. Here, in this work, we employ a framework that combines joint angular momentum and particle number projection with time-dependent configuration mixing to calculate the angular momentum distributions of primary fragments. Focusing on the benchmark cases of neutron-induced fission of 235 U and 239 Pu , we predict—for the first time—microscopic angular momentum distributions for all fragments observed in experiments. Our results reveal a pronounced sawtooth pattern in the average angular momentum as a function of fragment mass, consistent with recent measurements. Additionally, we observe substantial variations in angular momentum distributions along isobaric chains, indicating that commonly used empirical formulas lack sufficient accuracy. We also quantify a strong correlation between the angular momentum and the deformation of the fragments at scission, and a weak correlation in the magnitude of the angular momentum between fragment partners. The generated data will enable estimation of the impact of microscopic distributions on fission spectra, paving the way toward fission modeling based on microscopic inputs.

Physics - Nuclear physics and radiation physics↗