Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models & methods for nuclear reactions”

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

Comparing pion production in transport simulations of heavy-ion collisions at 270⁢𝐴 MeV under controlled conditions

Within the TMEP, we present a detailed study of the performance of different transport models in Sn+Sn collisions at 270A MeV, and put particular emphasis on the production of pions and Delta resonances, which have been used as probes of the nuclear symmetry energy. We prescribe a common and rather simple physics model, and follow in detail the results of 4 BUU models and 6 QMD models. The nucleonic evolution of the collision and the nucleonic observables in these codes do not completely converge, but the differences among the codes can be understood as being due to several reasons: the basic differences between BUU and QMD models in the representation of the phase-space distributions, computational differences in the mean-field evaluation, and differences in the adopted strategies for the Pauli blocking in the collision integrals. For pionic observables, we find that a higher maximum density leads to an enhanced pion yield and a reduced pi- pi+ yield ratio, while a more effective Pauli blocking generally leads to a slightly suppressed pion yield and an enhanced pi- pi+ yield ratio. We specifically investigate the effect of the Coulomb force, and find that it increases the total pi- pi+ yield ratio but reduces the ratio at high pion energies, although differences in its implementations do not have a dominating role in the differences among the codes. Taking into account only the results of codes that strictly follow the homework specifications, we find a convergence of the codes in the final charged pion yield ratio to a 1 sigma deviation of about 5%. Furthermore, the uncertainty is expected to be reduced to about 1.6% if the same or similar strategies and ingredients, i.e., an improved Pauli blocking and calculation of the non-linear term in the mean-field potential, are similarly used in all codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Proton and neutron density distributions at supranormal density in low- and medium-energy heavy-ion collisions. II. Central Pb + Pb collisions

This paper represents a continuation of our investigation into the limitations on total, proton, and neutron particle number densities, as well as the asymmetry of proton and neutron density distributions achievable in central heavy-ion collisions. We explore these aspects at low and medium energies within the framework of the Boltzmann-Uhlenbeck-Uehling (pBUU) transport and time-dependent Hartree-Fock (TDHF) models. Previous studies, focusing on symmetric and asymmetric collisions of Ca and Sn nuclei [1], and initial results on Pb-nuclei collisions [2], emphasized the role of the Coulomb interaction in these events. Our findings indicated that: (i) the highest total densities predicted at 𝐸 beam = 800 MeV/nucleon were on the order of ≈ 2.5⁢𝜌 0 (𝜌 0 = 0.16 fm −3 ), (ii) the proton-neutron asymmetry for maximal densities, 𝛿=(𝜌$^{max}_{n}$− 𝜌$^{max}_{p}$)/(𝜌$^{max}_{n}$+𝜌$^{max}_{p}$), did not generally exceed the asymmetry in the initial state of the collision at all beam energies and tended to decrease during the reaction, and (iii) a significant portion of this asymmetry had its microscopic origin in Coulomb forces, masking the pure nuclear contribution. These new findings, particularly relevant in the astrophysical context, are further examined in this work, focusing on the heaviest target-projectile combination 212,208 Pb accessible in an experiment. We introduce the SkT3 Skyrme force model, not previously used for the Pb system, and compare it to the SV-bas and SV-sym34 models to explore the symmetry-energy dependence of the results. Contour plots of nucleonic densities are presented, contrasting the time evolution of the density distributions in low (TDHF) and high (pBUU) models. We also present the evolution of normalized maximal proton, neutron, and total nucleon number density with increasing beam energy in the full pBUU model and the Vlasov approximation, aiming to explore the impact of correlations in the reaction. The time evolution of the proton and neutron density distributions in the plane transverse to the beam direction is illustrated at both low and high beam energy. In conclusion, our detailed examination of the Pb system in this work provides further essential evidence that the aforementioned findings (i)–(iii) are only weakly dependent on system size and a symmetry-energy model, and thus, they are of more general importance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Statistical tools for a better optical model

Background: Modern statistical tools provide the ability to compare the information content of observables and provide a path to explore which experiments would be most useful to give insight into and constrain theoretical models. Purpose: Here we study three such tools in the context of nuclear reactions with the goal of constraining the optical potential. Method: The three statistical tools examined are (i) the principal component analysis, (ii) the sensitivity analysis based on derivatives, and (iii) the Bayesian evidence. We first apply these tools to a toy-model case, comparing the form of the imaginary part of the optical potential. Then we consider two different reaction observables, elastic angular distributions and polarization data for reactions on 48 Ca and 208 Pb at two different beam energies. Results: For the toy-model case, we find significant discrimination power in the sensitivities and the Bayesian evidence, showing clearly that the volume imaginary term is more useful to describe scattering at higher energies. When comparing between elastic cross sections and polarization data using realistic optical models, sensitivity studies indicate that both observables are roughly equally sensitive but the variability of the optical model parameters is strongly angle dependent. The Bayesian evidence shows some variability between the two observables, but the Bayes factor obtained is not sufficient to discriminate between angular distributions and polarization. Conclusions: From the cases considered, we conclude that, in general, elastic scattering angular distributions have similar impact in constraining the optical potential parameters compared with the polarization data. The angular ranges for the optimum experimental constraints can vary significantly with the observable considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Low-energy 17 O(𝑛,𝛾)⁢ 18 O reaction within the microscopic potential model and its role for the weak 𝑟 process

The neutron radiative capture reaction 17 O ⁡(𝑛,𝛾) ⁢18 O plays a pivotal role in both nuclear structure studies and astrophysical nucleosynthesis, particularly in the formation of elements during hydrostatic and explosive stellar environments. We calculated the 17 O ⁡(𝑛,𝛾) ⁢18 O cross section within the Skyrme Hartree-Fock potential model and analyzed electric dipole 𝐸⁢1 transitions to both positive- and negative-parity states below the α-decay threshold in 18 O. Our cross sections are significantly different from the data available in commonly used libraries. We further investigate the impact of the new calculated cross section on weak 𝑟-process nucleosynthesis using large-scale reaction network calculations across a wide range of electron fractions and entropies. Our results show that the 17 O ⁡(𝑛,𝛾) ⁢18 O reaction rate significantly influences the production of first 𝑟-process peak elements, such as strontium, under specific astrophysical conditions. This study highlights the importance of accurate nuclear data for light isotopes in modeling heavy-element synthesis and provides updated reaction rates for future nucleosynthesis simulations.

6 ≤ A ≤ 19↗
DOE OSTI2024

Negative components in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ of heavy-ion fusion cross sections

Negative components are observed in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ in heavy-ion fusions related to recent studies of noticeable oscillations in the high-energy region of the fusion excitation functions. A modified multiGaussian model has been developed, which includes contributions for both positive and negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ representation. The model reproduces the heavy-ion fusion excitation functions very well, including oscillations. The property of these negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ spectrum and the reaction mechanism of fusion at energies above the Coulomb barrier are investigated. Finally, the model also reproduces well light heavy-ion fusion reactions between 12 C and 16 O at high energy; results imply that the compound nucleus channel effect may have an important contribution at the high energy region for all fusion systems.

low & intermediate energy heavy-ion reactions↗
DOE OSTI2024

Universal function for heavy-ion fusion cross sections

A universal function for heavy-ion fusion cross sections, Y = $\sqrt{π}$ Xerfc (- X ) + exp(- X 2 ) is proposed. By scaling both the cross section σ(E) and the energy E, heavy-ion fusion cross section data are found to follow closely a universal function over the whole energy range. The scaling is developed from either a simple, empirical single-Gaussian barrier distribution model for the representation d 2 (σE)/dE 2 , or the modified Siwek-Wilczynski model. The cross section expressions of these models are analytical, and can be easily used for all heavy-ion fusion excitation functions. Thus a bench marking of heavy-ion fusion excitation functions has been achieved. Finally, a general discussion regarding the universal function is given.

low & intermediate energy heavy-ion reactions↗
DOE OSTI2021

Measurement of U 238 ( n , n ' γ ) cross section data and their impact on reaction models

A better knowledge of (n,xn) reaction cross sections is important for both reaction modeling and energy applications. This article focuses on inelastic scattering of neutrons off 238 U for which improvements are needed to better constrain evaluations and solve inconsistencies in nuclear power reactor calculations. A new precise measurement of (n,xnγ) reaction cross sections on 238 U has been performed at the GELINA (Geel Electron LINear Accelerator) neutron facility operated by EC-JRC-Geel (Belgium) with the GRAPhEME (GeRmanium array for Actinides PrEcise MEasurements) setup. The prompt γ-ray spectroscopy method coupled to time-of-flight measurements is used to extract (n,xnγ) cross section values which can be further combined to infer the total neutron inelastic scattering cross section. Cross section data for 18 γ transitions (five never measured before) are presented and compared to the data in the literature. Here, emphasis is especially given to the uncertainty determination to produce partial cross section data as accurate as possible. Due to intrinsic limitations of the experimental method, the use of additional nuclear structure information coupled with theoretical modeling is required to determine the total (n, n') cross section over the whole neutron energy range. We have investigated modeling aspects of the 238 U(n,n'γ) cross sections related to the description of compound nucleus and preequilibirum mechanisms as well as the discrete part of nuclear structure. Through comparison between experimental and calculated (n,n'γ) cross sections, we pinpoint inaccuracies in the description of specific reaction mechanisms and challenge recently implemented models. This helps improving the whole modeling of the (n, n') reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Uncertainty quantification in breakup reactions

Breakup reactions are one of the favored probes to study loosely bound nuclei, particularly those in the limit of stability forming a halo. In order to interpret such breakup experiments, the continuum discretized coupled channel method is typically used. In this study, the first Bayesian analysis of a breakup reaction model is performed. Herein we use a combination of statistical methods together with a three-body reaction model (the continuum discretized coupled channel method) to quantify the uncertainties on the breakup observables due to the parameters in the effective potential describing the loosely bound projectile of interest. The combination of tools we develop opens the path for a Bayesian analysis of not only breakup processes, but also a wide array of complex processes that require computationally intensive reaction models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Cross sections for the formation of Rb84m,g, Rb83, and Rb82m in Sr86(d,x) reactions up to deuteron energies of 49 MeV: Competition between α-particle and multinucleon emission processes

Cross sections of Sr86(d,x) reactions leading to the products Rb84m,g, Rb83, and Rb82m were measured by the stacked-sample activation technique up to deuteron energies of 49 MeV. Nuclear model calculations were performed using the codes talys and empire, which combine the statistical, precompound, and direct interaction components. In all cases, the empire results were much higher than the talys calculation. Fairly good agreement was obtained between measured data and the talys calculation after some optimization of the input model parameters. Insight into competition between α-particle and multinucleon emission in the Y88 compound-nucleus system was also gained.

59 ≤ A ≤ 89↗
DOE OSTI2023

Quantifying uncertainties due to optical potentials in one-neutron knockout reactions

One-neutron knockout reactions have been widely used to extract information about the single-particle structure of nuclei from the valley of stability to the drip lines. The interpretation of knockout data relies on reaction models, where the uncertainties are typically not accounted for. In this work, we quantify uncertainties of optical potentials used in these reaction models and propagate them, for the first time, to knockout observables using a Bayesian analysis. We study two reactions in the present paper, the first of which involves a loosely bound halo projectile, 11 Be, and the second a tightly bound projectile, 12 C. We first quantify the parametric uncertainties associated with phenomenological optical potentials. Complementing this approach, we also quantify the model uncertainties associated with the chiral forces that can be used to construct microscopic optical potentials. For the phenomenological study, we investigate the impact of the imaginary terms of the optical potential on the breakup and stripping components of the knockout cross sections as well as the impact of the angular range. For the 11 Be case, the theoretical uncertainty from the phenomenological method is on the order of the experimental uncertainty of the knockout observables; however, for the 12 C case, the theoretical uncertainty is significantly larger. Additionally, the widths of the uncertainty bands for the knockout observables obtained for the microscopic study and the phenomenological approach are of similar orders of magnitude. Based on this work we conclude that structure information inferred from the ratio of the knockout cross sections will carry a theoretical uncertainty of at least 20% for halo nuclei and at least 40% for tightly bound nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Halo effective field theory analysis of one-neutron knockout reactions of Be 11 and C 15

Background: One-nucleon knockout reactions provide insightful information on the single-particle structure of nuclei. When applied to one-neutron halo nuclei, they are purely peripheral, suggesting that they could be properly modeled by describing the projectile within a halo effective field theory (halo-EFT). Purpose: We reanalyze the one-neutron knockout measurements of 11 Be and 15 C —both one-neutron halo nuclei—on beryllium at about 60 MeV/nucleon. We consider halo-EFT descriptions of these nuclei which already provide excellent agreement with breakup and transfer data. Method: Here we include a halo-EFT description of the projectile within an eikonal-based model of the reaction and compare its outcome to existing data. Results: Excellent agreement with experiment is found for both nuclei. The asymptotic normalization coefficients inferred from this comparison confirm predictions from ab initio nuclear-structure calculations and values deduced from transfer data. Conclusions: Halo-EFT can be reliably used to analyze one-neutron knockout reactions measured for halo nuclei and test predictions from state-of-the-art nuclear structure models on these experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Coupled fission fragment angular momenta

Here, nuclear fission produces fragments whose spins are coupled to the relative angular motion via angular momentum conservation. It is shown how ensembles of such spins can readily be obtained by either direct microcanonical sampling or by sampling of the associated normal modes of rotation. The resulting distribution of the spin-spin opening angle is illustrated in various three- and two-dimensional scenarios and it is demonstrated how recent mutually conflicting model calculations can be well reproduced with different assumptions about the scission geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗
DOE OSTI2023

Uncertainty-quantified phenomenological optical potentials for single-nucleon scattering

Optical-model potentials (OMPs) continue to play a key role in nuclear reaction calculations. However, the uncertainty of phenomenological OMPs in widespread use—inherent to any parametric model trained on data—has not been fully characterized, and its impact on downstream users of OMPs remains unclear. Here we assign well-calibrated uncertainties for two representative global OMPs, those of Koning-Delaroche and Chapel Hill '89, using Markov-chain Monte Carlo for parameter inference. By comparing the canonical versions of these OMPs against the experimental data originally used to constrain them, we show how a lack of outlier rejection and a systematic underestimation of experimental uncertainties contributes to bias of, and overconfidence in, best-fit parameter values. Our updated, uncertainty-quantified versions of these OMPs address these issues and yield complete covariance information for potential parameters. Scattering predictions generated from our ensembles show improved performance both against the original training corpora of experimental data and against a new “test” corpus comprising many of the experimental single-nucleon scattering data collected over the last twenty years. Finally, we apply our uncertainty-quantified OMPs to two case studies of application-relevant cross sections. We conclude that, for many common applications of OMPs, including OMP uncertainty should become standard practice. Furthermore, to facilitate their immediate use, digital versions of our updated OMPs and related tools for forward uncertainty propagation are included as Supplemental Material.

150 ≤ A ≤ 189↗
DOE OSTI2023

Collective enhancement in the exciton model

The preequilibrium reaction mechanism is considered in the context of the exciton model. A modification to the one-particle–one-hole state density is studied which can be interpreted as a collective enhancement. The magnitude of the collective enhancement is estimated by simulating the Lawrence Livermore National Laboratory (LLNL) pulsed-spheres neutron-leakage spectra. Here the impact of the collective enhancement is explored in the context of the highly deformed actinide, 239 Pu. A consequence of this enhancement is the removal of fictitious levels in the distorted-wave Born approximation often used in modern nuclear reaction codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

α -cluster microscopic study of C 12 + C 12 fusion toward the zero energy limit

The carbon burning process is a fundamental step of stellar evolution and governs the synthesis of chemical elements important for the formation of life. In this work, we utilize the microscopic hybrid α cluster (HαC) model and an analytical approach, both in the framework of the Imaginary Time Method (ITM), to study the carbon fusion reaction towards zero energy. We obtain the values of the cross sections, astrophysical factors and correlate our results to collective motion. We also include a calculation for the 2 + carbon fusion and discuss a possible experimental investigation. Our results confirm direct experimental and theoretical results close to the barrier, while suggest possible 2 + mixtures in the indirect experimental data. Furthermore, our study offers an accurate view of the burning process in the somewhat unexplored low energy region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum physics of stars

Stars are slowly developing objects; the lifetimes of the different burning phases are determined by the strength of nuclear reactions, which in turn are defined by the quantum structure of the associated nuclei at the threshold and the respective reaction mechanisms. Stars, from the nuclear physics perspective, are cold environments where only a few of the key nuclear reactions have been measured at the actual stellar plasma temperatures. This is also the case for more dynamic astrophysical phenomena from the big bang to stellar explosions. Most of the nuclear reaction rates are therefore based on theoretical extrapolations. A number of discrepancies between these predictions and the associated stellar signatures have been observed, and many may be due to low-energy or near-threshold quantum effects. These effects need to be understood in order to reliably model nuclear reaction processes, not only for stars but also for low-temperature plasma environments such as controlled magnetic or inertial confinement fusion systems, which operate in similar temperature regimes. This review summarizes the various theoretical techniques presently used for deriving reaction rates and discusses possible quantum effects that may impact the reaction cross section near the reaction threshold. These resemble enhanced single-particle and cluster structures near threshold and associated interference effects. New experimental techniques such as deep-underground accelerators or the study of transfer reactions to mimic the quantum-mechanical transition strength, the so-called Trojan horse method, provide ways to directly or indirectly probe the reaction features that determine the reaction rates at stellar energies. Furthermore, this is demonstrated on a number of key nuclear reactions for different nucleosynthesis environments. Finally, current inconsistencies between experimental predictions and observations are discussed.

Models & methods for nuclear reactions↗