Comparison of validation methods of simulations for final state interactions in hadron production experiments
Not Available
SEARCH · Search NASA
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.
Not Available
Background: Proton elastic scattering at intermediate energy is widely employed as a tool for determining the matter radius of atomic nuclei. Here, the sensitivity of the approach relies on high-resolution measurements at small scattering angles and low-momentum transfer. Under these conditions, the Glauber multiple scattering theory accurately describes the proton-nucleus elastic cross section. Purpose: Investigate the sensitivity of the Glauber multiple scattering theory to uncertainties associated with input parameters such as the nuclear-matter density distribution and nucleon-nucleon data. Method: A joint Bayesian inference was performed using 12 angular distributions of elastic scattering at different energies on 58 Ni, 90 Zr, and 208 Pb targets. A Metropolis-Hastings algorithm was implemented to make an uncertainty quantification analysis for the input parameters used in the Glauber multiple scattering theory. Results: The experimental cross sections were fitted simultaneously using a joint Bayesian inference approach. Posterior probability density distributions of 42 input parameters were obtained from the analysis. A moderate correlation between the nuclear density parameters and the nucleon-nucleon cross sections was found. This correlation impacts the extraction of the nuclear-matter radius. Conclusions: The present analysis provided a consistent method for extracting the nuclear-matter density distribution of 58 Ni, 90 Zr, and 208 Pb from data across different incident energies. Due to the correlation of the nucleon-nucleon cross sections with the other input parameters, a constrained Bayesian inference using free nucleon-nucleon cross section data was performed. The nuclear-matter radii obtained from the analysis are in good agreement with multiple results reported in the literature.
A number of accelerator-based isotope production facilities utilize $100-200$ MeV proton beams due to the high production rates enabled by high-intensity beam capabilities and the greater diversity of isotope production brought on by the long range of high-energy protons. However, nuclear reaction modeling at these energies can be challenging because of the interplay between different reaction modes and a lack of existing guiding cross section data. A Tri-lab collaboration has been formed between the Lawrence Berkeley, Los Alamos, and Brookhaven National Laboratories to address these complexities by characterizing charged-particle nuclear reactions relevant to the production of established and novel radioisotopes. In the inaugural collaboration experiments, stacked-targets of niobium foils were irradiated at the Brookhaven Linac Isotope Producer ($E_p=200$ MeV) and the Los Alamos Isotope Production Facility ($E_p=100$ MeV) to measure $^{93}$Nb(p,x) cross sections between $50-200$ MeV. The results were compared with literature data as well as the default calculations of the nuclear model codes TALYS, CoH, EMPIRE, and ALICE. The default code predictions largely failed to reproduce the measurements. Therefore, we developed a standardized procedure, which determines the reaction model parameters that best reproduce the most prominent reaction channels in a physically justifiable manner. Overall, the primary focus of the procedure was to determine the best parameterization for the pre-equilibrium two-component exciton model. This modeling study revealed a trend towards a relative decrease for internal transition rates at intermediate proton energies ($E_p=20-60$ MeV) in the current exciton model as compared to the default values. The results of this work are instrumental for the planning, execution, and analysis essential to isotope production.
The evaluated $^{16}O(n,α)$ cross sections in the ENDF/B-VIII.0 nuclear data library remain uncertain because of systematic discrepancies in the measured data. In the energy region below the first excited state, R-matrix analyses rely heavily on the inverse reaction, and in particular, they rely on the measured $^{13}C(α,n)^{16}O$ cross sections reported by Bair et al. in 1973 and Harissopulos et al. in 2005. The Harissopulos cross section values are systematically lower than those previously reported by Bair. Here, drawing on the available experimental information, this paper briefly describes and demonstrates that the two sets of measured cross sections are consistent.
The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the accurate calculation of cross sections for the different reaction channels. Many common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account in these models: LD at neutron separation energy coming from average spacings of s- and p-wave resonances ( D 0 and D 1 , respectively) whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only in very specific regions of excitation energy, and for specific spins and parities. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for fitting and the determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD as a starting point as the HFB LD provide a more realistic spin and parity distributions than phenomenological models such as Gilbert-Cameron (GC). We then associate variations predicted by the HFB model with the structure observed in double-differential cross sections at low outgoing neutron energy, a region that is dominated by the LD input. We also use (n, p) on 56 Fe , as an example case where angle-integrated cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or values of D 0 or D 1 are unknown. Finally, it also predicts neutron-induced inelastic γ cross sections that in some cases can differ significantly from more standard phenomenological LD models such as GC.
We report the Hauser-Feshbach fission-fragment decay (HF 3 D) model is extended to calculate the prompt fission neutron spectrum (PFNS) for the thermal-neutron-induced fission on 235 U, where the evaporated neutrons from all possible fission-fragment pairs are aggregated. By studying model parameter sensitivities on the calculated PFNS as well as nonstatistical behavior of low-lying discrete level spin distribution, we conclude that discrepancies between the aggregation calculation and the experimental PFNS seen at higher neutron emission energies can be attributed to both the primary fission-fragment yield distribution and the possible high spin states that are not predicted by the statistical theory of nuclear structure.
The 10 B + n system has been studied by measuring charged particles from neutron-induced reactions in the 1 to 20 MeV energy range. Protons, deuterons, tritons, and α particles were measured at four angles using the pulsed white neutron spectrum at the Los Alamos Neutron Science Center Weapons Neutron Research facility. Differential cross sections for each species are reported. These new data are combined with literature data in an R-matrix analysis of the 11 B compound system. Furthermore, a comparison is made to the molecular and cluster states candidates observed in these reactions.
Radioactive decays from 42 Ar and its progeny 42 K are potential background sources in large-scale liquid-argon-based neutrino and dark matter experiments. In the atmosphere, 42 Ar is produced primarily by cosmogenic activation on 40 Ar. The use of low radioactivity argon from cosmogenically shielded underground sources can expand the reach and sensitivity of liquid-argon-based rare event searches. We estimate 42 Ar production underground by nuclear reactions induced by natural radioactivity and cosmic-ray muon-induced interactions. At 3,000 mwe, 42 Ar production rate is 1.8×10 -3 atoms per ton of crust per year, 7 orders of magnitude smaller than the 39 Ar production rate at a similar depth in the crust. In conclusion, by comparing the calculated production rate of 42 Ar to that of 39 Ar for which the concentration has been measured in an underground gas sample, we estimate the activity of 42 Ar in gas extracted from 3,000 mwe depth to be less than two decays per ton of argon per year.
A single crystal chemical vapor deposited (sCVD) diamond detector is used as an active target to measure neutron-induced reactions on natural carbon using the neutrons produced by spallation, with a broad energy spectrum at LANSCE. Additionally, the neutron-induced reactions are detected in the diamond as low as E n = 400 keV and up to approximately 100 MeV. Relative cross sections for C 12 ( n , α 0 ) , C 12 ( n , p 0 ) , C 12 ( n , d 0 + p 1 ) , and C 13 ( n , α 0 ) are reported up to E n = 22 MeV and comparisons on detected pulse-height spectra and detector response of scattering reactions are made with GEANT4 simulations using the ENDF/B-VIII.0 evaluated nuclear data library up to 20 MeV. The results are compared with past experimental data, including other works that incorporate diamond detectors as an active carbon target. In addition, R-matrix calculations for the C 13 + n system are presented.
We report nuclear reaction data for neutron induced reactions on unstable nuclei are critical for a wide range of applications spanning studies of nuclear astrophysics, nuclear reactor designs, and radiochemistry diagnostics. However, nuclear data evaluations of the reaction cross sections are largely based on calculations due to the difficulty in performing this class of measurements and the resulting lack of experimental data. For neutron induced charged particle reactions at fast neutron energies, at the MeV scale, these cross section predictions are predominately driven by statistical Hauser-Feshbach calculations. In this work, we present partial and total 59 Ni(n, p) and 59 Ni(n, α) cross sections, measured directly with a radioactive 59 Ni target, and compare the results to the present nuclear data evaluations. In addition, the results from this work are compared to a recent study of the 59 Ni(n, xp) reaction cross section that was performed via an indirect surrogate ratio method. The expected energy trend of the cross section, based on the current work, is inconsistent with that of the surrogate work. This calls into question the reliability of that application of the surrogate ratio method and highlights the need for direct measurements on unstable nuclei, when feasible.
Herein we report the measurement of the total kinetic energy (TKE) release in the fast neutron induced fission of 240 Pu and 242 Pu. The results are compared to the predictions of the GEF model, the CGMF model, and the model of Denisov and Sedykh as well as previous exptl. work on these reactions. Our absolute measurements of the TKE release are in good agreement with the previous measurements of Nethaway et al. for the interaction of 14.8 MeV neutrons with 240 Pu [Phys. Rev. C16, 1907 (1977)] and of Winkelmann and Aumann for the interaction of 15 MeV neutrons with 242 Pu [Phys. Rev. C30, 934 (1984)]. The general trends of the measured TKE values agree with phenomenol. models but the variances of the TKE distributions are significantly less than predicted by various models. The mean postneutron emission TKE release decreases nonlinearly with increasing neutron energy and can be represented as TKE(MeV) = 175.8 ± 0.3 - (2.4 ± 0.8) lo g 10 E n - (1.4 ± 0.4) lo g 10 $E^{2}_{n}$ for 240 Pu and TKE(MeV) = 177.1 ± 0.3 - (1.2 ± 0.9) lo g 10 E n - (1.8 ± 0.5)lo g 10 $E_{n}^{2}$ for 242 Pu.
The characterization of fission-driven nuclear systems primarily relies on calculations of neutron-induced chain reactions, and these calculations require evaluated nuclear data as input. Calculation accuracy heavily depends on input nuclear data evaluation accuracy, and thus high precision on the experimental input to the nuclear data evaluation is essential for fundamental quantities like the energy spectrum of neutrons emitted from neutron-induced fission (i.e., the prompt fission neutron spectrum, PFNS). Despite decades of measurement efforts, prior to the measurements described in this work there were only three literature data sets for the 235 U(n,f) PFNS at incident neutron energies above 1.0 MeV considered reliable for inclusion in nuclear data evaluations and no reliable data sets above 3.0 MeV incident neutron energy. In this work we report on new measurements of the 235 U(n,f) PFNS spanning a grid of 1.0–20.0 MeV in incident neutron energy and 0.01–10.0 MeV in outgoing (PFNS) neutron energy. These measurements were carried out at the Weapons Neutron Research facility at the Los Alamos Neutron Science Center and used a multifoil parallel-plate avalanche counter target with both a Li-glass and a liquid scintillator detector array in separate experiments to span the quoted outgoing neutron energy ranges. The PFNS results are shown in terms of the energy spectra themselves as well as the average PFNS energy $(\langle{E}\rangle)$ and ratios of $\langle{E}\rangle$ at forward and backward angles. Here, the results are compared with literature data and selected nuclear data evaluations. Generally, the data agree with the ENDF/B-VIII.0 evaluation below 5.0-MeV incident neutron energy and more closely with the JEFF-3.3 evaluation above 5.0 MeV, though no evaluations considered for comparison in this work agree with the data across all of the incident and outgoing neutron energies shown, especially in regions where the third-chance fission process becomes available. Additionally, we show a ratio of the present PFNS results for 235 U(n, f) with a recent and highly correlated experiment to measure the 239 Pu(n, f) PFNS at the same experimental facility and with nearly identical equipment and analysis procedures. Many observations reported in this work are the first of their kind and represent significant advancements for knowledge of the 235 U(n, f) PFNS.
We estimate the amount of 37 Ar produced in natural xenon via cosmic-ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth’s surface. We then calculate the resulting 37 Ar concentration in a 10-tonne payload (similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage, and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea-level production rate of 37 Ar in natural xenon is estimated to be 0.024 atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1 tonne/month, the average 37 Ar activity after 10 tons are purified and transported underground is 0.058 − 0.090 𝜇Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic 37 Ar will appear as a noticeable background in the early science data, while decaying with a 35-day half-life. This newly noticed production mechanism of 37 Ar should be considered when planning for future liquid-xenon-based experiments.
Measurements of the 182 W(n, 2n) 181 W cross section have been performed in the neutron energy range between 8 and 15 MeV using the activation technique. Such data are needed to help interpret results of laser shots at the National Ignition Facility using a new DT capsule design, featuring a high-Z inner shell, with tungsten as the favored material, and an outer shell made of a low-Z material. Our data are in very good agreement with the previous data of Frehaut et al., which are based on a different technique, and in fair agreement with the ENDF/B-VIII.0 and the JEFF-3.3 evaluations.
Not Available
Not Available
The radionuclide 22 Na is a target of γ-ray astronomy searches, predicted to be produced during thermonuclear runaways driving classical novae. The 22 Na (p, γ) 23 Mg reaction is the main destruction channel of 22 Na during a nova, hence, its rate is needed to accurately predict the 22 Na yield. However, experimental determinations of the resonance strengths have led to inconsistent results. In this Rapid Communication, we report a measurement of the branching ratios of the 23 Al β-delayed protons as a probe of the key 204-keV (center-of-mass) 22 Na (p, γ) 23 Mg resonance strength. We report a factor of 5 lower branching ratio compared to the most recent literature value. The variation in 22 Na yield due to nuclear data inconsistencies was assessed using a series of hydrodynamic nova outburst simulations and has increased to a factor of 3.8, corresponding to a factor of ≈ 2 uncertainty in the maximum detectability distance. Finally, this is the first reported scientific measurement using the Gaseous Detector with Germanium Tagging system.
We formulate microscopic neutron-nucleus optical potentials from many-body perturbation theory based on chiral two- and three-body forces. The neutron self-energy is first calculated in homogeneous matter to second order in perturbation theory, which gives the central real and imaginary terms of the optical potential. Here, the real spin-orbit term is calculated separately from the density matrix expansion using the same chiral interaction as in the self-energy. Finally, the full neutron-nucleus optical potential is derived within the improved local density approximation utilizing mean-field models consistent with the chiral nuclear force employed. We compare the results of the microscopic calculations to phenomenological models and experimental data up to projectile energies of E = 200 MeV. Experimental elastic differential scattering cross sections and vector analyzing powers are generally well reproduced by the chiral optical potential, but we find that total cross sections are overestimated at high energies.