Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nucleon–Nucleon potentials”

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 37 records · Page 2

Order-by-order uncertainties of nucleon-nucleon Wolfenstein amplitudes in chiral effective field theory

Quantum mechanical invariance principles dictate the most general operator structure that can be present in the nucleon-nucleon (NN) interaction. Five independent operators appear in the on-shell NN amplitude together with five corresponding coefficient functions. The usual choice for these coefficient functions is known as the NN Wolfenstein amplitudes. We analyze the order-by-order convergence of each of the five NN Wolfenstein amplitudes predicted by a semilocal coordinate space potential implementation of chiral effective field theory (𝜒⁢EFT). We do this at laboratory kinetic energies between 25 and 200 MeV for both neutron-proton and proton-proton scattering. Our analysis uses the Gaussian-process methods developed by the BUQEYE Collaboration to describe the contributions of each 𝜒⁢EFT order, and so yields truncation uncertainties for each Wolfenstein amplitude that are correlated across scattering angles. We combine information on the size of different orders in the EFT to infer the 𝜒⁢EFT breakdown scale for each amplitude, finding, on average, Λ 𝑏 between 750 and 800 MeV. Furthermore, with this choice of Λ 𝑏 , the EFT truncation uncertainties cover both higher-order results and empirical Wolfenstein amplitudes well for all orders other than the leading order.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗

Probing axionlike particles with multimessenger observations of neutron star mergers

Axion-like particles (ALPs) can be copiously produced in binary neutron star (BNS) mergers through nucleon-nucleon bremsstrahlung if the ALP-nucleon couplings 𝑔 𝑎⁢𝑁 are sizable. The ALP-photon coupling 𝑔 𝑎⁢𝛾 may trigger conversions of ultralight ALPs into photons in the magnetic fields of the merger remnant and of the Milky Way. This effect would lead to a potentially observable short gamma-ray signal, in coincidence with the gravitational-wave signal produced during the merging process. This event could be detected through multimessenger observation of BNS mergers employing the synergy between gravitational-wave detectors and gamma-ray telescopes. Here, we study the sensitivity of current and proposed MeV gamma-ray experiments to detect such a signal. As an explicit example, we consider ALP couplings related as in the Kim-Shifman-Vainshtein-Zakharov axion model, and show that in this case the proposed instruments can reach a sensitivity down to 𝑔 𝑎⁢𝛾 ≳ few ×10 −13 GeV −1 for 𝑚 𝑎 ≲ 10 −9 eV, comparable with the SN 1987A limit.

Axion-like particles↗

Neutron-neutron quasifree scattering in neutron-deuteron breakup at 10 MeV

New measurements of the neutron-neutron quasifree scattering cross section in neutron-deuteron breakup at an incident neutron energy of 10.0 MeV are reported. The experiment setup was optimized to evaluate the technique for determining the integrated beam-target luminosity in neutron-neutron coincidence cross-section measurements in neutron-deuteron breakup. The measurements were carried out with a systematic uncertainty of ± 5.6%. Our data are in agreement with theoretical calculations performed using the CD-Bonn nucleon-nucleon potential in the Faddeev formalism. The measured integrated cross section over the quasifree peak is 20.5 ± 0.5(stat) ± 1.1 (sys) mb/sr 2 in comparison with the theory prediction of 20.1 mb/sr 2 . Overall, these results validate our technique for determining the beam-target luminosity in neutron-deuteron breakup measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Benchmark calculations of pure neutron matter with realistic nucleon-nucleon interactions

Here, we report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner–Bethe–Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne v 18 interaction and two of its simplified versions, and four of the new Norfolk Δ -full chiral effective field theory potentials. Primarily because of the advancements in the auxiliary-field diffusion Monte Carlo, we observe good agreement among the three many-body techniques up to nuclear saturation density—the maximum difference in the energy per particle is within 1.5 MeV for all the potentials we consider. At higher densities, the divergences become more important, and are mainly driven by the Fermi hypernetted chain/single-operator calculations. We also study the connection between nucleon-nucleon scattering data and the energy per particle of pure neutron matter. Our results suggest that fitting to higher-energy nucleon-nucleon scattering helps reduce the spread of energies among the models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantifying Uncertainties In Neutron-Alpha Scattering With Chiral Nucleon-Nucleon And Three-Nucleon Forces

We report that modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Starting from NN potentials up to fifth order (N 4 LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (c E and c D ) contribute to the overall uncertainty budget of many-body calculations of neutron- 4 He (n-α) elastic scattering. We compute structure and reaction observables for three-, four-, and five-nucleon systems within the ab initio frameworks of the no-core shell model and no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte Carlo. We find rapidly converging n-α phase shifts with respect to the chiral order. With the adopted leading-order 3N force, calculations based on the NN interaction at N 4 LO of Entem, Machleidt, and Nosyk are unable to reproduce the experimental phase shifts in the 3/2 - channel within the estimated chiral truncation errors. Closer agreement with empirical data is found when using an older parametrization of the NN interaction at order N 3 LO, and the position and width of the P-wave resonances can be used to reduce the uncertainty of the 3N low-energy constants. The present results point to a lack of spin-orbit strength when the newer parametrization of the chiral NN force up to fifth order is combined with the leading-order 3N force. The inclusion of higher-order 3N-force terms may be required to recover the missing strength. GPMs can act as fast and accurate emulators of ab initio many-body calculations of low-energy scattering and reactions of light nuclei, opening the way to a robust quantification of theoretical uncertainties grounded in the description of the underlying chiral Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Convergence of nuclear effective field theory with perturbative pions

The classic paper by Fleming, Mehen, and Stewart (FMS) cast doubts on the convergence of spin-triplet nucleon-nucleon partial wave scattering amplitudes when following the proposal of Kaplan, Savage, and Wise to construct nuclear effective field theory (EFT) around the unitary fermion limit with perturbative pion exchange. FMS identified the subclass of iterated one-pion exchange potential graphs as the cause of this poor convergence, which they showed persisted in the chiral limit. Theoretical tools are developed in this work to compute these Feynman graphs analytically to high order in all angular momentum channels simultaneously, examining the amplitudes computed to seven loops in the $\textit{L = J}$ channels and three loops in the coupled $\textit{L = J}$ ± 1 channels. One finds that there is nothing pathological about the perturbative expansion of a 1/$r^3$ potential in general and that the expansion converges satisfactorily in all partial waves except those with the lowest angular momentum, particularly the $^3P_0$ and the coupled $^3S_1 - ^3D_1$ channels. The results corroborate work by Birse, which suggests possible avenues to explore for improving the range of validity of the EFT expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ab initio leading order effective potentials for elastic nucleon-nucleus scattering

Calculating microscopic effective interactions (optical potentials) for elastic nucleon-nucleus scattering has already in the past led to a large body of work. For first-order calculations a nucleon-nucleon ( NN ) interaction and a one-body density of the nucleus were taken as input to rigorous calculations of microscopic full-folding calculations. Based on the spectator expansion of the multiple scattering series we employ a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction on the same footing in the structure as well as in the reaction calculation to obtai n an in leading-order consistent effective potential for nucleon-nucleus elastic scattering, which includes the spin of the struck target nucleon. The first order effective folding potential is computed by first deriving a nonlocal scalar density as well as a spin-projected momentum distribution. Those are then integrated with the off-shell Wolfenstein amplitudes $A$, $C$, and $M $. The resulting nonlocal potential serves as input to a momentum-space Lippmann-Schwinger equation, whose solutions are summed to obtain the nucleon-nucleus scattering observables. As a result, we calculate elastic scattering observables for $^4$He, $^6$He, $^8$He, $^{12}$C, and $^{16}$O in the energy r egime between 100 and 200~MeV projectile kinetic energy, and compare to available data. We also explore the extension down to about 70 MeV, and study the effect of ignoring the spin of the struck nucleon in the In our calculations we contrast elastic scattering off closed-shell and open-shell nuclei. We find that fo r closed-shell nuclei the approximation of ignoring the spin of the struck target nucleon is excellent. We only see effects of the spin of the struck target nucleon when considering $^6$He and $^8$He, which are nuclei with a $N/Z$ ratio larger than 1.

6 ≤ A ≤ 19↗

Charge radii of exotic neon and magnesium isotopes

We compute the charge radii and ground-state energies of even-mass neon and magnesium isotopes from neutron number N = 8 to the dripline. Our calculations are based on nucleon-nucleon and three-nucleon potentials from chiral effective field theory that include Δ isobars. These potentials yield an accurate saturation point and symmetry energy of nuclear matter. We use the coupled-cluster method and start from an axially symmetric reference state. Binding energies and two-neutron separation energies largely agree with data, and the dripline in neon is accurate. The computed charge radii are accurate for many isotopes where data exist. Finer details, such as isotope shifts, however, are not accurately reproduced. These chiral potentials indicate a subshell closure at N = 14 for the radii (but not for two-neutron separation energies) and a decrease in charge radii at N = 8 (observed in neon and predicted for magnesium). Furthermore, they yield a continued increase of charge radii as neutrons are added beyond N = 14 yet underestimate the large increase at N = 20 in magnesium.

20 ≤ A ≤ 38↗

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of the neutron matter equation of state and the symmetry energy up to fourth order of chiral effective field theory

We present predictions for the neutron matter equation of state, from leading to fourth order of chiral effective field theory, using recently developed, accurate chiral nucleon-nucleon potentials. For the many-body method, we employ the nonperturbative particle-particle ladder approximation, that is, we solve the G-matrix equation. Furthermore, we find the impact of subleading three-neutron forces to be mild and attractive. We also show order-by-order predictions for the symmetry energy, and discuss its density dependence in relation to empirical constraints. For the nuclear matter equation of state, in this work we adopt an empirical parametrization with good saturation properties. This is to highlight, specifically, the energy and pressure in neutron matter, particularly when comparing with empirical constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleon-nucleon potentials from Δ-full chiral effective-field-theory and implications

We closely investigate NN potentials based upon the Δ-full version of chiral effective-field-theory. We find that recently constructed NN potentials of this kind, which (when applied together with three-nucleon forces) were presented as predicting accurate binding energies and radii for a range of nuclei from A = 16 to A = 132 and providing accurate equations of state for nuclear matter, yield a χ 2 /datum of 60 for the reproduction of the pp data below 100 MeV laboratory energy. This χ 2 is more than three times what the Hamada-Johnston potential of the year of 1962 achieved already some 60 years ago. We perceive this historical fact as concerning in view of the current emphasis on precision. We are able to trace the very large χ 2 as well as the apparent success of the potentials in nuclear structure to unrealistic predictions for P-wave states, in which the Δ-full next-to-next-to-leading order (NNLO) potentials are off by up to 40 times the NNLO truncation errors. In fact, we show that the worse the description of the P-wave states, the better the predictions in nuclear structure. Thus, these potentials cannot be seen as the solution to the outstanding problems in current microscopic nuclear structure physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Overview of symmetric nuclear matter properties from chiral interactions up to fourth order of the chiral expansion

We present and discuss predictions for a cross section of bulk and single-particle properties in symmetric nuclear matter based on recent high-quality nucleon-nucleon potentials at N 3 LO and including all subleading three-nucleon forces. We begin with the equation of state and its saturation properties and proceed to the single-nucleon potential. We also explore short-range correlations as seen through the defect function. The various predictions which we present have a common foundation in an internally consistent ab initio approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Elastic proton scattering off nonzero spin nuclei

In recent years, we constructed a microscopic optical potential (OP) for elastic nucleon-nucleus ( NA ) scattering using modern approaches based on chiral theories for the nucleon-nucleon ( NN ) interaction. The OP was derived at first order of the spectator expansion in Watson multiple scattering theory and its final expression was a folding integral between the NN t matrix and the nuclear density of the target. Two- and three-body forces are consistently included both in the target and in the projectile description. The purpose of this work is to apply our microscopic OP to nuclei characterized by a ground state of spin-parity quantum numbers J π ≠ 0 + . We extended our formalism to include the spin of the target nucleus. The full amplitudes of the NN reaction matrix are retained in the calculations starting from two- and three-body chiral forces. The microscopic OP can be applied in the energy range 100 ≤ E ≤ 350 MeV. Finally, we show a remarkable agreement with experimental data for the available observables and, simultaneously, provide reliable estimates for the theoretical uncertainties. This work paves the way toward a full microscopic approach to inelastic NA scattering, showing that the derivation of optical potentials between states with J π ≠ 0 + is completely under control.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measuring recoiling nucleons from the nucleus with the future Electron Ion Collider

Short range correlated nucleon-nucleon (NN) pairs are an important part of the nuclear ground state. They are typically studied by scattering an electron from one nucleon in the pair and detecting its spectator correlated partner (spectator-nucleon tagging). The Electron Ion Collider (EIC)should be able to detect these nucleons, since they are boosted to high momentum in the lab frame by the momentum of the ion beam. To determine the feasibility of these studies with the planned EIC detector configuration, we have simulated quasi-elastic scattering for two electron and ion beam energy configurations: 5 GeV electrons and 41 GeV/A ions, and 10 GeV electrons and 110 GeV/A ions. We show that the knocked-out and recoiling nucleons can be detected over a wide range of initial nucleon momenta. We also show that these measurements can achieve much larger momentum transfers than current fixed target experiments. By detecting both low and high initial-momentum nucleons,the EIC will provide the data that should allow scientists to definitively show if the EMC effect and short-range correlation are connected, and to improve our understanding of color transparency

Hauenstein, F. (ORCID:0000000212652212)↗

Nuclear states projected from a pair condensate

Atomic nuclei exhibit deformation, pairing correlations, and rotational symmetries. To meet these competing demands in a computationally tractable formalism, we revisit the use of general pair condensates with good particle number as a trial wave function for even-even nuclei. After minimizing the energy of the condensate, allowing for general triaxial deformations, we project out states with good angular momentum with a fast projection technique. To show applicability, we present example calculations from pair condensates in several model spaces and compare against angular-momentum projected Hartree-Fock and full configuration-interaction shell-model calculations. This approach successfully generates spherical, vibrational, and rotational spectra, demonstrating potential for modeling medium- to heavy-mass nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Applicability of semiclassical methods for modeling laser-enhanced fusion rates in a realistic setting

In the context of the potential laser-induced enhancement to the rates of DHe 3 and DT fusion, we discuss the frequently-used Wentzel-Kramers-Brillouin (WKB) method and the imaginary-time method (ITM). For static external electric fields, we find that these methods predict significant enhancement to the fusion cross section for electric-field strengths > 10 14 V/m, especially at low values (≈ keV) for the enter-of-mass (CoM) energy. When considering dynamic electric fields, this enhancement can be amplified by considering increased photon frequencies. However, we also pro- vide a review of the region of laser-parameter phase space where these semiclassical methods are applicable. Here, we conclude that this allowable region decreases for higher photon frequencies in con-junction with lower values for the electric-field strength, motivating the need for future experiments to test the predictions of these methods and their ranges of validity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Renormalization of one-pion exchange in chiral effective field theory for antinucleon-nucleon scattering

The renormalization of iterated one-pion exchange (OPE) is studied in Chiral Effective Field Theory ($\chi$EFT) for the antinucleon-nucleon ($\overline{N}N$) system. The OPE potential is cut off at a certain distance and contact interactions are represented by a complex spherical well with the same radius. We investigate the dependence on the cutoff radius of the phase shifts, inelasticities, and mixing angles for the low partial waves in $\overline{N}N$ scattering. We show that renormalization requires additional contact interactions compared to the expectation based on naive dimensional analysis. Results after renormalization are compared with the state-of-the-art energy-dependent partial-wave analysis of $\overline{N}N$ data. We compare our conclusions with applications of $\chi$EFT to the nucleon-nucleon system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗