Boundary-condition constraints for the shell model - a method for nuclear structure and nuclear reactions.
Asymptotic boundary conditions as constraints in shell model for scattering, reactions and nuclear structure calculations
NUCLEAR STRUCTURE↗
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.
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.↗
Impact of composition and symmetry energy on the temperature of quasiprojectiles simulated with antisymmetrized molecular dynamics
The equation of state describes the emergent physical properties of matter. Experimental data is needed to help constrain the equation of state for nuclear matter. These constraints can help distinguish between an “asy-stiff” and an “asy-soft” equation of state, which has astrophysical implications. One path to help constrain the models is to analyze the nuclear caloric curve; some experiments have shown dependence on neutron excess, and may thus be sensitive to the asymmetry. A difference in the caloric curve based on the asymmetry of the reconstructed quasiprojectile (QP) had been observed using 70 Zn on 70 Zn at 35 MeV/nucleon taken with the NIMROD array. Antisymmetrized molecular dynamics calculations were performed for the same system and deexcited with gemini++. Both Gogny (asy-soft) and Gogny-as (asy-stiff) data sets were generated. The particles were then filtered based on detector geometric acceptance and thresholds. From the accepted particles, the excitation energy and temperature were calculated in the same way as for experimental data. Additionally, filter effects on the observed nuclear caloric curves were investigated. A tendency for the asy-stiff nuclear caloric curves to have higher temperatures than their asy-soft counterparts was observed for a number of probes. In addition, some probes may show sensitivity to the reconstructed composition of the QP, but this is inconclusive due to high statistical fluctuations and a large dependence on the exact method of event selection.
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↗
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↗
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↗
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↗
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↗
Impact of tensor forces on quasifission product yield distributions
Quantum shell effects are crucial for the stability and structure of atomic nuclei and play a key role in the discovery of superheavy elements. Furthermore, during nuclear collisions, the dynamical evolution of these shell effects disrupts the equilibration process necessary for forming a compound nucleus, leading to the breakup of the initial composite as a result of quasifission. As such, quasifission reactions hinder the formation of a superheavy element.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Learning nuclear cross sections across the chart of nuclides with graph neural networks
We explore the use of deep learning techniques to learn how nuclear cross sections change as we add or remove protons and neutrons. As a proof of principle, we focus on the neutron-induced reactions in the fast energy regime. Our approach follows a two-stage learning framework. First, we apply representation learning to encode cross section data into a latent space using either variational autoencoders (VAEs) or implicit neural representations (INRs). Then, we train graph neural networks (GNNs) on the resulting embeddings to predict missing values across the nuclear chart by leveraging the topological structure of neighboring isotopes. We demonstrate accurate cross section predictions within a 9 × 9 block of missing nuclei. We also find that the optimal GNN training strategy depends on the type of latent representation used, with VAE embeddings performing best under end-to-end optimization in the original space, while INR embeddings achieve better results when the GNN is trained only in the latent space. Furthermore, using clustering algorithms, we map groups of latent vectors into regions of the nuclear chart and show that VAEs and INRs can discover some of the neutron magic numbers. These findings suggest that deep-learning models based on the representation encoding of cross sections combined with graph neural networks hold significant potential in augmenting nuclear theory models, e.g., by providing reliable estimates of covariances of cross sections, including cross-material covariances.
Machine learning↗
Single-channel and single-energy partial-wave analysis with continuity improved through minimal phase constraints
Single-energy partial-wave analysis has often been applied as a way to fit data with minimal model dependence. However, remaining unconstrained, partial waves at neighboring energies will vary discontinuously because the overall amplitude phase cannot be determined through single-channel measurements. This problem can be mitigated through the use of a constraining penalty function based on an associated energy-dependent fit. However, the weight given to this constraint results in a biased fit to the data. In this paper, for the first time, we explore a constraining function which does not influence the fit to data. The constraint comes from the overall phase found in multichannel fits which, in the present study, are the Bonn-Gatchina and Jülich-Bonn multichannel analyses. The data are well reproduced and weighting of the penalty function does not influence the result. The method is applied to KΛ photoproduction data and all observables can be maximally well reproduced. While the employed multichannel analyses display very different multipole amplitudes, we show that the major difference between two sets of multipoles can be related to the different overall phases.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Probabilistic neural networks for improved analyses with phenomenological R -matrix
Here we present a method for measurement analyses based on probabilistic deep neural networks that provide several advantages over conventional analyses with phenomenological models. These include predicting physical quantities directly from data, the rapid generation of statistically robust uncertainties, and the ability to bypass some parameters that may induce ambiguities and complications in data analysis. As deep learning methods make predictions through “black boxes,” the uncertainty quantification is typically challenging. We use a probabilistic framework that provides thorough uncertainty quantification and is straightforward to follow in practice. With the network architecture based on the Transformer, we demonstrate the current method for predicting nuclear resonance parameters from scattering data using the phenomenological R-matrix model.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Test of the boundary-condition constraint method for nuclear reactions.
Boundary condition constraint method for generalizing R-matrix theory, permitting use of shell model as basis for nuclear-reaction calculations
GREEN THEOREM↗
Development of A Monte Carlo Radiation Transport Code System For HEDS: Status Update
Modifications of the Monte Carlo radiation transport code HETC are underway to extend the code to include transport of energetic heavy ions, such as are found in the galactic cosmic ray spectrum in space. The new HETC code will be available for use in radiation shielding applications associated with missions, such as the proposed manned mission to Mars. In this work the current status of code modification is described. Methods used to develop the required nuclear reaction models, including total, elastic and nuclear breakup processes, and their associated databases are also presented. Finally, plans for future work on the extended HETC code system and for its validation are described.
Townsend, Lawrence W.↗
Gamma Decay of the 154 Sm Isovector Giant Dipole Resonance: Smekal-Raman Scattering as a Novel Probe of Nuclear Ground-State Deformation
𝛾 decays of the isovector giant dipole resonance (IVGDR) of the deformed nucleus 154 Sm were measured using 2$^{+}_{1}$-Smekal-Raman and elastic scattering of linearly polarized, quasimonochromatic photon beams. The two scattering processes were disentangled through their distinct angular distributions. Their branching ratio and cross sections were determined at six excitation energies covering the 154 Sm IVGDR. Both agree with the predictions of the geometrical model for the IVGDR and confirm 𝛾 decay as an observable sensitive to the structure of the resonance. Consequently, the data place strong constraints on the nuclear shape, including the degree of triaxiality. The derived 154 Sm shape parameters 𝛽 = 0.2925(25) and 𝛾 = 5.0(15)° agree well with other measurements and recent Monte Carlo shell-model calculations.
150 ≤ A ≤ 189↗