Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “models and 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 19 records

Constraints from isoscaling on the source size in energetic heavy ion collisions

In the framework of the statistical multifragmentation model, the nuclear isoscaling analysis is extended to constrain the ratio between the sizes of the decaying sources formed in a collision between two heavy ions. It is found that the ratio between the probabilities of observing n fragments in each event, for each of the sources, follows a scaling law, similar to the traditional nuclear isoscaling. However, the corresponding slope is sensitive to the source sizes. This property is explained analytically using the grand-canonical ensemble. Furthermore, the extent to which our findings are affected by finite size effects and by the deexcitation of the hot primary fragments is investigated. The scaling turns out to be robust and weakly affected by effects implied by these two aspects. We also find that the Poisson distribution is a fairly good approximation to the above mentioned probabilities, associated with both the primordial fragments, produced at the breakup stage, and the final ones, found at the end of the fragment deexcitation process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Considering nonlocality in the optical potentials within eikonal models

Background: For its simplicity, the eikonal method is the tool of choice to analyze nuclear reactions at high energies (E > 100 MeV/nucleon), including knockout reactions. However, so far, the effective interactions used in this method are assumed to be fully local. Purpose: Given the recent studies on nonlocal optical potentials, in this work we assess whether nonlocality in the optical potentials is expected to impact reactions at high energies and then explore different avenues for extending the eikonal method to include nonlocal interactions. Method: We compare angular distributions obtained for nonlocal interactions (using the exact R-matrix approach for elastic scattering and the adiabatic distorted wave approximation for transfer) with those obtained using their local-equivalent interactions. Results: Our results show that transfer observables are significantly impacted by nonlocality in the high-energy regime. Because knockout reactions are dominated by stripping (transfer to inelastic channels), nonlocality is expected to have a large effect on knockout observables too. Three approaches are explored for extending the eikonal method to nonlocal interactions, including an iterative method and a perturbation theory. Conclusions: None of the derived extensions of the eikonal model provide a good description of elastic scattering. Here, this paper suggests that nonlocality removes the formal simplicity associated with the eikonal model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quenching of single-particle strengths in direct reactions

A discrepancy in the asymmetry dependence of spectroscopic factors extracted with different reaction probes calls into question whether the corresponding reaction models are properly understood. In this work, we present extracted spectroscopic factors from the Ar 46 , 34 ( p , d ) Ar 45 , 33 transfer reactions in inverse kinematics at a beam energy of 70 MeV/nucleon. The results are consistent with previous measurements of these reactions at a lower beam energy [Lee et al. , Phys. Rev. Lett. 104 , 112701 (2010) ], indicating that the transfer reaction is a reliable probe for the nuclear structure of exotic nuclei across a wide energy range. Results from a large body of transfer reaction measurements, ( p , p N ) measurements, and theoretical nuclear structure studies make a compelling case for much weaker asymmetry dependence than what is observed with single-nucleon knockout reactions on beryllium or carbon targets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On modelling nuclear reactions in meteorites

An original method of experimental modeling depth distribution of radionuclides in sphere thick targets irradiated by protons in 4(pi)-geometry on JINR (Dubna) synchrocyclotron is described. Some results are presented.

Ustinova, G. K.↗

Investigating the time dependence of neutron-proton equilibration using molecular dynamics simulations

Neutron-proton equilibration has previously been studied experimentally in dynamically deformed nuclei in heavy ion collisions [Phys. Rev. Lett. 118, 062501 (2017); Phys. Rev. C 95, 044604 (2017)] using the NIMROD detector array at the Cyclotron Institute at Texas A&M University. Results indicated the composition of the two heaviest fragments originating from the excited projectile-like fragment evolved exponential with respect to their orientation angle to be more similar. Constrained Molecular Dynamics and Antisymmetrized Molecular Dynamics simulations were performed for different formulations of the density dependence of the asymmetry energy term of the nuclear equation of state. The simulations are compared to experimental results, which indicate a better agreement with a softer interaction in the nuclear equation of state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comparison of heavy-ion transport simulations: Mean-field dynamics in a box

Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared among each other and to exact calculations. For BUU-like codes, employing the test particle method, the results depend on the combination of the number of test particles and the spread of the profile functions that weight integration over space. These parameters can be properly adapted to give a good reproduction of the analytical zero-sound features. QMD-like codes, using molecular dynamics methods, are characterized by large damping effects, attributable to the fluctuations inherent in their phase-space representation. Moreover, for a given nuclear effective interaction, they generally lead to slower density oscillations, as compared to BUU-like codes. The latter problem is mitigated in the more recent lattice formulation of some of the QMD codes. Furthermore, the significance of these results for the description of real heavy-ion collisions is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Systematic study of the validity of the eikonal model including uncertainties

Nuclear reactions at intermediate beam energies are often interpreted using the eikonal model. In the analysis of complex reaction probes, where few-body reaction methods are needed, the eikonal method may be used as an efficient way for describing the fragment-target reaction process. In this work, we perform a systematic study to test the validity of the eikonal approximation for nucleon-nucleus reactions. We also quantify uncertainties due to the nucleon optical potential on reaction observables. We inspect the validity of the eikonal model and its semiclassical correction by comparing it to exact solutions (obtained from solving the optical-model equation with a finite-differences method) for a wide range of reactions. We also study the effect of relativistic corrections, both kinematic and dynamic, by effectively incorporating the relativistic effects at intermediate energies. The uncertainties from a Bayesian global optical potential (KDUQ) are propagated to the observables of interest. Our study includes neutron and proton reactions on 27 Al , 40 Ca , 90 Zr , and 208 Pb , for a wide range of energies 𝐸 lab = 0–400 MeV. We calculate neutron-total cross sections (elastic and reactions) as well as proton-absorption cross sections as a function of beam energy, using the eikonal model, the eikonal model with a semiclassical correction, and the exact solution. Here, we also compute angular distributions for the methods above. Our results show that for the proton-absorption cross section, the eikonal model can be used down to around 60 MeV and the semiclassical correction extends its use to 30 MeV. However, the validity of the eikonal model for the neutron-total cross section only goes down to ≈120 MeV, a range extended to ≈ 50 MeV when using the semiclassical correction. We find the semiclassical correction to the eikonal model to be less effective in describing the angular distributions. The 1⁢𝜎 uncertainty intervals on the observables we studied is less than 5% for most of the energies considered, but increases rapidly for higher energies, namely energies outside the range of KDUQ (𝐸 lab > 200MeV).

Cluster models↗

First simultaneous measurement of the γ-ray and neutron emission probabilities in inverse kinematics at a heavy-ion storage ring

The probabilities for γ-ray and particle emission as a function of the excitation energy of a decaying nucleus are valuable observables for constraining the ingredients of the models that describe the deexcitation of nuclei near the particle emission threshold. These models are essential in nuclear astrophysics and applications. In this paper, we have for the first time simultaneously measured the γ-ray and neutron emission probabilities of 208 Pb. The measurement was performed in inverse kinematics at the Experimental Storage Ring (ESR) of the GSI/FAIR facility, where a 208 Pb beam interacted through the 208 Pb (p,p') reaction with a hydrogen gas jet target. Instead of detecting the γ rays and neutrons emitted by 208 Pb, we detected the heavy beamlike residues produced after γ and neutron emission. These heavy residues were fully separated by a dipole magnet of the ESR and were detected with outstanding efficiencies. The comparison of the measured probabilities with model calculations has allowed us to test and select different descriptions of the γ-ray strength function and the nuclear level density available in the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗