Microscopic optical potentials derived from ab initio translationally invariant nonlocal one-body densities
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
We formulate microscopic optical potentials for nucleon-nucleus scattering from chiral two- and three-nucleon forces. The real and imaginary central terms of the optical potentials are obtained from the nucleon self-energy in infinite nuclear matter at a given density and isospin asymmetry, calculated self-consistently to second order in many-body perturbation theory. Furthermore, the real spin-orbit term is extracted from the same chiral potential using an improved density matrix expansion. The density-dependent optical potential is then folded with the nuclear density distributions of 40,42,44,48 Ca from which we study proton-nucleus elastic scattering and total reaction cross sections using the reaction code talys. We compare the results of the microscopic calculations to those of phenomenological models and experimental data up to projectile energies of E = 180 MeV. While overall satisfactory agreement with the available experimental data is obtained, we find that the elastic scattering and total reaction cross sections can be significantly improved with a weaker imaginary optical potential, particularly for larger projectile energies.
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.
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.
Here, we develop a method to construct a microscopic optical potential from chiral interactions for nucleon-nucleus scattering. The optical potential is constructed by combining the Green’s function approach with the coupled-cluster method. To deal with the poles of the Green’s function along the real energy axis we employ a Berggren basis in the complex energy plane combined with the Lanczos method. Using this approach, we perform a proof-of-principle calculation of the optical potential for the elastic neutron scattering on 16 O. For the computation of the ground-state of 16 O, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A ±1 nuclei we use particle-attached/removed equation-of-motion method truncated at two-particle-one-hole and one-particle-two-hole excitations, respectively. We verify the convergence of the optical potential and scattering phase shifts with respect to the model-space size and the number of discretized complex continuum states. We also investigate the absorptive component of the optical potential (which reflects the opening of inelastic channels) by computing its imaginary volume integral and find an almost negligible absorptive component at low-energies. To shed light on this result, we computed excited states of 16 O using equation-of-motion coupled-cluster method with singles-and- doubles excitations and we found no low-lying excited states below 10 MeV. Furthermore, most excited states have a dominant two-particle-two-hole component, making higher-order particle-hole excitations necessary to achieve a precise description of these core-excited states. We conclude that the reduced absorption at low-energies can be attributed to the lack of correlations coming from the low-order cluster truncation in the employed coupled-cluster method.
Abstract We review recent progress and motivate the need for further developments in nuclear optical potentials that are widely used in the theoretical analysis of nucleon elastic scattering and reaction cross sections. In regions of the nuclear chart away from stability, which represent a frontier in nuclear science over the coming decade and which will be probed at new rare-isotope beam facilities worldwide, there is a targeted need to quantify and reduce theoretical reaction model uncertainties, especially with respect to nuclear optical potentials. We first describe the primary physics motivations for an improved description of nuclear reactions involving short-lived isotopes, focusing on its benefits for fundamental science discoveries and applications to medicine, energy, and security. We then outline the various methods in use today to build optical potentials starting from phenomenological, microscopic, and ab initio methods, highlighting in particular, the strengths and weaknesses of each approach. We then discuss publicly-available tools and resources facilitating the propagation of recent progresses in the field to practitioners. Finally, we provide a set of open challenges and recommendations for the field to advance the fundamental science goals of nuclear reaction studies in the rare-isotope beam era. This paper is the outcome of the Facility for Rare Isotope Beams Theory Alliance (FRIB-TA) topical program ‘Optical Potentials in Nuclear Physics’ held in March 2022 at FRIB. Its content is non-exhaustive, was chosen by the participants and reflects their efforts related to optical potentials.
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.
Recently, a variety of studies have shown the importance of including nonlocality in the description of reactions. The goal of this work is to revisit the phenomenological approach to determining nonlocal optical potentials from elastic scattering. We perform a χ 2 analysis of neutron elastic scattering data off 40 Ca, 90 Zr, and 208 Pb at energies E ≈ 5–40 MeV, assuming a Perey and Buck [Nucl. Phys. 32, 353 (1962)] or Tian et al. [Int. J. Mod. Phys. E 24, 1550006 (2015)] nonlocal form for the optical potential. We introduce energy and asymmetry dependencies in the imaginary part of the potential and refit the data to obtain a global parametrization. Independently of the starting point in the minimization procedure, an energy dependence in the imaginary depth is required for a good description of the data across the included energy range. We then present two parametrizations, both of which represent an improvement over the original potentials for the fitted nuclei as well as for other nuclei not included in our fit. Our results show that, even when including the standard Gaussian nonlocality in optical potential, a significant energy dependence is required to describe elastic-scattering data.
We perform a first-principles calculation of optical potentials for nucleon elastic scattering off medium-mass isotopes. Fully based on a saturating chiral Hamiltonian, the optical potentials are derived by folding nuclear density distributions computed with self-consistent Green's function theory with a nucleon-nucleon t matrix computed with a consistent chiral interaction. The dependence on the folding interaction as well as the convergence of the target densities are investigated. Numerical results are presented and discussed for differential cross sections and analyzing powers, with focus on elastic proton scattering off calcium and nickel isotopes. Our optical potentials generally show a remarkable agreement with the available experimental data for laboratory energies in the range 65–200 MeV. We study the evolution of the scattering observables with increasing proton-neutron asymmetry by computing theoretical predictions of the cross section and analyzing power over the calcium and nickel isotopic chains. Published by the American Physical Society 2024
The inelastic (n,n’) reactions on actinide nuclei cannot easily be measured, especially those which populate states of low excitation energies. Model calculations are therefore need to interpolate between (and extrapolate from) the data for those discrete and continuum states that can be measured. Such models are critically dependent on several ingredients: the choice of neutron optical potentials, the collective deformations of these deformed actinides, and the level densities. I here focus on uncertainties arising from the optical potentials, especially when different optical potentials may be fitted almost as well to the same data, yet produce different predictions for what must be inferred indirectly.
The description of final-state interactions (FSI) in the large phase space probed in neutrino experiments poses a great challenge. In neutrino experiments, which operate under semi-inclusive conditions, cascade models are commonly used for this task, while under exclusive conditions FSI can be treated with relativistic optical potentials (ROP). We formulate conditions under which the ROP approach and cascade model can be directly compared. We feed the NEUT cascade with events from a relativistic distorted-wave impulse approximation calculation that uses the real part of an optical potential. Cuts on the missing energy of the resulting events are applied to define a set of events that can be directly compared to RDWIA calculations with the full optical potential. The NEUT cascade and ROP agree for proton kinetic energies $T_p > 150$ MeV for carbon, oxygen and calcium nuclei when a realistic nuclear density is used to introduce events in the cascade. For $T_p < 100$ MeV the ROP and NEUT cross sections differ in shape and differences in magnitude are larger than 50 %. Single transverse variables allow to distinguish different approaches to FSI, but due to a large non-QE contribution the comparison to T2K data does not give an unambiguous view of FSI. We discuss electron scattering and argue that with a cut in missing energy FSI can be studied with minimal confounding factors in e.g. $e4\nu$. The agreement of the ROP and NEUT for T2K conditions lends confidence to these models as a tool in oscillation analyses for sufficiently large nucleon kinetic energies. These results urge for caution when a cascade model is applied for small nucleon energies. The assessment of model assumptions relevant to this region are strongly encouraged. This paper provides novel constraints on cascade models from proton-nucleus scattering that can be easily applied to other neutrino event generators.
Precise modeling of neutrino interactions on nuclear targets is essential for neutrino oscillations experiments. The modeling of the energy of final state particles in quasielastic (QE) scattering and resonance production on bound nucleons requires knowledge of both the removal energy of the initial state bound nucleon as well as the average Coulomb and nuclear optical potentials for final state leptons and hadrons. We extract the average values of the real part of the nuclear optical potential for final state nucleons ( $U_{opt}^{QE}$ U opt QE ) as a function of the nucleon kinetic energy from inclusive electron scattering data on nuclear targets ( $_\mathbf{6 }^\mathbf{12 }{} \mathbf{C} $ 6 12 C + $_\mathbf{8 }^\mathbf{16 }{} \mathbf{O} $ 8 16 O , $_\mathbf{20 }^\mathbf{40 }{} \mathbf{Ca} $ 20 40 Ca + $_\mathbf{18 }^\mathbf{40 }{} \mathbf{Ar} $ 18 40 Ar , $_\mathbf{3 }^\mathbf{6 }{} \mathbf{Li} $ 3 6 Li , $_\mathbf{18 }^\mathbf{27 }{} \mathbf{Al} $ 18 27 Al , $_\mathbf{26 }^\mathbf{56 }{} \mathbf{Fe} $ 26 56 Fe , $_\mathbf{82 }^\mathbf{208 }{} \mathbf{Pb} $ 82 208 Pb ) in the QE region and compare to calculations. We also extract values of the average of the real part of the nuclear optical potential for a $\varDelta (1232)$ Δ ( 1232 ) resonance in the final state ( $U^\varDelta _{opt}$ U opt Δ ) within the impulse approximation. We find that $U^\varDelta _{opt}$ U opt Δ is more negative than $U_{opt}^{QE}$ U opt QE with $U^\varDelta _{opt}\approx $ U opt Δ ≈ 1.5 $U_{opt}^{QE}$ U opt QE for $_\mathbf{6 }^\mathbf{12 }{} \mathbf{C} $ 6 12 C .
Uncertainty quantification for nuclear theories has gained a more prominent role in the field, with more groups attempting to understand the uncertainties in their calculations. However, recent studies have shown that the uncertainties on the optical potentials are too large for the theory to be useful. The goal of this work is to explore possible experimental conditions that may reduce the uncertainties on elastic scattering and single-nucleon transfer cross sections that come from the fitting of the optical model parameters to experimental data. Using Bayesian methods, we explore the effect of the uncertainties of optical model parameters on the angular grid of the differential cross section, including cross-section data at nearby energies, and changes in the experimental error bars. We also study the effect on the resulting uncertainty when other observables are included in the fitting procedure, particularly the total (reaction) cross sections. We study proton and neutron elastic scattering on 48 Ca and 208 Pb. We explore the parameter space with the Markov-chain Monte Carlo method, produce posterior distributions for the optical model parameters, and construct the corresponding 95% confidence intervals on the elastic-scattering cross sections. We also propagate the uncertainties on the optical potentials to the 48 Ca(d,p) 49 Ca(g.s.) and 208 Pb(d,p) 209 Pb(g.s.) cross sections. We find little sensitivity to the angular grid and an improvement of up to a factor of 2 on the uncertainties by including data at a nearby energy. Although reducing the error bars on the data does reduce the uncertainty, the gain is often considerably smaller than one would naively expect. We also find that the inclusion of total reaction cross section can improve the uncertainty although the magnitude of the effect depends strongly on the cases considered.
We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for 40 Ca and 48 Ca. For the computation of the ground states of 40 Ca and 48 Ca, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = ±1 nuclei we use particle-attached and particle-removed equation-of-motion method truncated at two-particle–one-hole and one-particle–two-hole excitations, respectively. Our calculations are based on the chiral nucleon-nucleon and three-nucleon interaction NNLO sat , which reproduces the charge radii of 40 Ca and 48 Ca, and the chiral nucleon-nucleon interaction NNLO opt . In all cases considered here, we observe that the overall form of the neutron scattering cross section is reproduced for both interactions, but the imaginary part of the potential, which reflects the loss of flux in the elastic channel, is negligible. Furthermore, the latter points to neglected many-body correlations that would appear beyond the coupled-cluster truncation level considered in this work. We show that, by artificially increasing the parameter η in the Green's function, practical results can be further improved.
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.
We combine the ab initio symmetry-adapted no-core shell model (SA-NCSM) with the single-particle Green's function approach to construct optical potentials rooted in first principles. Specifically, we show that total cross sections and phase shifts for neutron elastic scattering from a 4 He target with projectile energies between 0.5 and 10 MeV closely reproduce the experiment. In addition, we discuss an important new development that resolves a long-standing issue with spurious center-of-mass motion in the Green's function formalism for many-body approaches. Furthermore, the new development opens a path for first-principle predictions of cross sections for elastic scattering of single-nucleon projectiles, nucleon capture, and deuteron breakup reactions, feasible for a broad range of open-shell spherical and deformed nuclei in the SA-NCSM approach.
Nanodisc technology was implemented as a platform for voltage nanosensors. A fluorescence (Förster) resonance energy transfer (FRET)-based voltage-sensing scheme employing fluorescent nanodiscs and the hydrophobic ion dipicrylamine was developed and utilized to optically record membrane potentials on the single-nanodisc level. Ensemble and single-nanosensor recordings were demonstrated for HEK293 cells and primary cortical neuron cells. Conjugation of nanodiscs to anti-GABA A antibodies allowed for site-specific membrane potential measurements from postsynaptic sites.