Search NASA⌕ Search

SEARCH · Search NASA

Results for “optical, coupled-channel and distorted wave models”

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

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↗

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↗

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗

Near-Threshold Dipole Strength in 10 Be with Isoscalar Character

Isoscalar dipole transitions are a distinctive fingerprint of cluster structures. A 1 − resonance at 7.27(10) MeV, located just below the 𝛼-emission threshold, has been observed in the deuteron inelastic scattering reactions off 10 Be. The deformation lengths of the excited states in 10 Be below 9 MeV have been inferred from the differential cross sections using coupled channel calculations. This observed 1 − resonance has isoscalar characteristics and exhausts approximately 5%–15% of the isoscalar dipole energy-weighted sum rule, providing evidence for pronounced 𝛼 cluster structure in 10 Be. The Gamow coupled channel approach supports this interpretation and suggests the near-threshold effect might be playing an important role in this excitation energy domain. Here, the 𝛼 +𝛼 +𝑛 +𝑛 four-body calculation reproduces the observed enhanced dipole strength, implying that the four-body cluster structure is essential to describe the 1 − states in 10 Be.

Cluster models↗

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-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↗

Global Microscopic Description of Nucleon-Nucleus Scattering with Quantified Uncertainties

We develop for the first time a microscopic global nucleon-nucleus optical potential with quantified uncertainties suitable for analyzing nuclear reaction experiments at next-generation rare-isotope beam facilities. Within the improved local density approximation and without any adjustable parameters, we begin by computing proton-nucleus and neutron-nucleus optical potentials from a set of five nuclear forces from chiral effective field theory for 1800 target nuclei in the mass range 12 ≤ A ≤ 242 for energies between 0 MeV < E ≲ 150 MeV. We then parameterize a global optical potential for each chiral force that depends smoothly on the projectile energy as well as the target nucleus mass number and isospin asymmetry. Uncertainty bands for elastic scattering observables are generated from a full covariance analysis of the parameters entering in the description of our global optical potential and benchmarked against existing experimental data for stable target nuclei. Since our approach is purely microscopic, we anticipate a similar quality of the model for nucleon scattering on unstable isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Reduction of the neutron imaginary potential off the stability line and its possible impact on neutron capture rates

The effect of the isovector imaginary optical potential has been studied experimentally by using the particle-evaporation technique for the 11 B + 48 Ca reaction with a 21.8 MeV 11 B beam. Spectra of neutron, proton, and α particles emitted from the neutron-rich compound nucleus 59 Mn have been measured and analyzed with traditional optical model potentials with their original parametrizations as well as with adjusted isovector imaginary components. It is shown that the isovector component of the imaginary potential is indispensable in the reproduction of proton and α-particle yields from this reaction and even needs to be enhanced compared with the suggestions of the original model parametrizations. Furthermore, this can lead to important consequences for astrophysical reaction-rate calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Properties of isoscalar giant multipole resonances in medium-heavy closed-shell nuclei: A semimicroscopic description

The semimicroscopic particle-hole dispersive optical model (PHDOM) is implemented to describe main properties of Isoscalar Giant Multipole Resonances (up to L=3) in medium-heavy closed-shell nuclei. Here, the main properties are characterized by the strength distribution, transition density, partial and total probabilities of direct one-nucleon decay. Calculation results obtained for the 208 Pb nucleus are compared with available experimental data.

190 ≤ A ≤ 219↗

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗

Uncertainty quantification of optical models in fission fragment deexcitation

Here, we take the first step towards incorporating compound nuclear observables at astrophysically relevant energies into the experimental evidence used to constrain optical models, by propagating the uncertainty in two global optical potentials, one phenomenological and one microscopic, to correlated fission observables using the Monte Carlo Hauser-Feshbach formalism. We compare to a wide range of historic and recent experimental fission measurements, and discuss in detail regions of disagreement. We find that the parametric optical model uncertainty in neutron-fragment correlated observables involving neutron energy is significant. On the other hand, we observe that other experimental features, particularly neutron-fragment correlations near the 132 Sn shell closure and the high energy component of neutron spectra, are unlikely to be explained by the optical potential, and will require further experimental and theoretical effort to explain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Elastic scattering in the 12 N + 197 Au system at 𝐸 lab = 70 MeV

Halo structure is an interesting exotic configuration developed in some light weakly bound nuclei, where a valence particle orbits a nuclear core. Signatures of halo structure can be observed in the angular distributions of the elastic scattering induced by these nuclei at energies around the Coulomb barrier. There are some well-studied reactions with neutron-rich halo nuclei, such as 6 He and 11 Li . However, the information is scarce on the proton-rich side. Recent works confirm the halo structure in the 8B nuclei but still lack more experimental studies for other proton-halo candidates, such as 12 N and 17 F. In this work, we report experimental data for the elastic scattering of 12 N on 197 Au target at 𝐸 lab = 70 MeV. The 12 N is a proton-rich nucleus with proton separation energy 𝑆 𝑝 = 600 keV, which is higher than the 8 B (𝑆 𝑝 = 137 keV) and almost the same as the 17 F (𝑆 𝑝 = 601 keV). Data were obtained at the Cyclotron Institute of Texas A&M University where the 12 N radioactive beam was produced by the momentum achromatic recoil spectrometer. The calculation of the optical model was used to fit the measured angular distribution for the elastic scattering and to obtain the reaction cross section 𝜎 𝑅 . We also performed continuum discretized coupled-channel calculations to compare with the experimental data. The angular distribution of the elastic cross sections exhibits a suppression of the Fresnel peak. From the fitting of the optical model, we obtain the total reaction cross section, 𝜎 𝑅 = 1269 ± 41 mb. The agreement between the CDCC calculation and the experimental elastic scattering data is limited and the breakdown does not exhaust the measured 𝜎 𝑅 . The resulting reduced reaction cross section 𝜎 Red for 12 N + 197 Au is large and comparable to the one obtained for the 8 B + 208 Pb system. This suggests a strong decoupling of the valence proton from the core because of the low binding energy and a dynamic polarization effect. Here, more research is required to estimate the contribution of core excitation.

Elastic scattering reactions↗

Dynamical model of φ meson photoproduction on the nucleon and 4 He

We investigate φ meson photoproduction on the nucleon and the 4 He targets within a dynamical model approach based on a Hamiltonian which describes the production mechanisms by the Pomeron-exchange, meson-exchanges, φ radiations, and nucleon resonance excitations mechanisms. The final φN interactions are included being described by the gluon-exchange, direct φN couplings, and the box-diagrams arising from the couplings with πN, ρN, KΛ, and KΣ channels. The parameters of the Hamiltonian are determined by the experimental data of γp → φp from the CLAS Collaboration. The resulting Hamiltonian is then used to predict the coherent φ-meson production on the 4 He targets by using the distorted-wave impulse approximation. For the proton target, the final φN rescattering effects, as required by the unitarity condition, are found to be very weak, which supports the earlier calculations in the literature. For the 4 He targets, the predicted differential cross sections are in good agreement with the data obtained by the LEPS Collaboration. Furthermore, the role of each mechanism in this reaction is discussed and predictions for a wide range of scattering angles are presented, which can be tested in future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Microscopic calculation of fission product yields for odd-mass nuclei

Fission data are essential inputs to reaction networks involved in nucleosynthesis simulations and nuclear forensics. In such applications as well as in the description of multichance fission, the characteristics of fission for odd-mass nuclei are just as important as those for even-even nuclei. The fission theories that aim at explicitly describing fission dynamics are typically based on some variant of the nuclear mean-field theory. In such cases, the treatment of systems with an odd number of particles is markedly more involved, both formally and computationally. In this article, we use the blocking prescription of the Hartree-Fock-Bogoliubov theory with Skyrme energy functionals to compute the deformation properties of odd-mass uranium isotopes. We show that the resulting fission fragment distributions depend quite significantly on the spin of the odd neutron. By direct calculation of the spin distribution of the fissioning nucleus, we propose a methodology to rigorously predict the charge and mass distributions in odd-mass nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Unitary coupled-channel three-body amplitude with pions and kaons

Three-body dynamics above threshold is required for the reliable extraction of many amplitudes and resonances from experiment and lattice QCD. The S-matrix principle of unitarity can be used to construct dynamical coupled-channel approaches in which three particles scatter off each other, rearranging two-body subsystems by particle exchange. This paper reports the development of a three-body coupled-channel, amplitude including pions and kaons. The unequal-mass amplitude contains two-body S- and P-wave subsystems (“isobars”) of all isospins, 𝐼 = 0,1/2,1,3/2,2 , and it also allows for transitions within a given isobar. The 𝑓 0 ⁡(500)⁢(𝜎),𝑓 0 ⁡(980),𝜌⁡(700),𝐾$^{*}_{0}$⁡(700)⁢(𝜅), and 𝐾*⁡(892) resonances are included, apart from repulsive isobars. Different methods to evaluate the amplitude for physical momenta are discussed. Production amplitudes for 𝑎 1 quantum numbers are shown as a proof of principle for the numerical implementation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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↗

Semimicroscopic description of isoscalar giant multipole resonances in medium-mass closed-shell nuclei

Main properties of isoscalar giant multipole resonances (up to L=3), including L=0,and 2 overtones, in medium-mass closed-shell nuclei are described within the semimicroscopic particle-hole dispersive optical model. Further, main properties are characterized by the energy-averaged strength distribution, projected (one-body) transition density, and probabilities of direct one-nucleon decay. Calculation results obtained for characteristics of the mentioned resonances in the 48 Ca, 90 Zr, and 132 Sn nuclei are compared with available experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Role of the likelihood for elastic scattering uncertainty quantification

In the last decade, uncertainty quantification (UQ) for optical model potentials (OMPs) has become a focal point for nuclear reaction theory, and several competing approaches for OMP UQ have recently been developed. Here, we clarify recent efforts to compare frequentist and Bayesian approaches in the context of OMP UQ [G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019)]. We replicate a portion of that OMP UQ study but use independent statistical tools. Specifically, we compare two methods for OMP parameter inference from elastic scattering data: the Levenberg-Marquardt algorithm for χ 2 minimization on one hand and Markov chain Monte Carlo (MCMC) sampling on the other. Separately, we assess the common practice of using a renormalized likelihood (χ 2 /N), N being the number of data points, instead of the canonical weighted-least-squares likelihood (χ 2 ), as a way of accounting for unknown data correlations. Here, we show that for a generic linear model and for a five-parameter OMP analysis, frequentist and uniform-prior Bayesian approaches recover the same optimum and uncertainty estimates—not systematically larger uncertainties for the Bayesian approach, as was concluded in G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019). Further, we show that if an additional, near-degenerate parameter is introduced into the same OMP analysis such that the parameter posterior becomes non-Gaussian, then covariance-based estimates of uncertainty become unreliable. Finally, we show that regardless of optimization approach, if χ 2 /N is used for the likelihood, the resulting parametric uncertainties increase by $\sqrt{N}$, and that this is responsible for the conclusions drawn in the revisited study. Based on our replication results, we find that a fortuitous cancellation of unreported errors and the renormalization factor can lead to improvement in empirical coverages, as was the case in the original comparative study. We emphasize that developing and applying a realistic likelihood function is an essential task in a UQ analysis, and that several recent UQ studies that employed a renormalized likelihood (i.e., including a 1/N factor) may have yielded unrealistically large uncertainties for elastic-scattering observables. If the parameter posterior deviates from multivariate-normal, a sampling-based approach like MCMC has a clear advantage over methods that assume the Laplace approximation holds. We note that empirical coverage can serve as an important internal check for the analyst whose model or data may have additional, unaccounted-for uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗