Search

We present a novel approach to calculate theoretical uncertainties in few-nucleon calculations that makes use of automatic differentiation. We demonstrate this method in deuteron bound state and nucleon - nucleon scattering calculations. Backpropagation, implemented in the Python pytorch library, is used to calculate the gradients with respect to model parameters and propagate errors from these parameters to the deuteron binding energy and selected phase-shift parameters. The uncertainty values obtained using this approach are validated by directly sampling from the potential parameters. We find very good agreement between two ways of estimating that uncertainty., Comment: 20 pages, 7 figures, 5 tables

We discuss new spin observables presently accessible to measurement, namely polarization transfer coefficients from doubly spin-polarized initial state to the outgoing proton in the elastic proton-deuteron (pd) scattering and in the deuteron breakup reactions. The sensitivity of these observables to three-nucleon force (3NF) effects is investigated and compared to sensitivities of the constituent standard single polarization transfer coefficients. $K_{y,y}^{y'}$ in elastic pd scattering, for which large 3NF effects, up to 40\%, have been found at higher energies, seems the most promising observable to measure., Comment: 23 pages, 7 figures

A measurement of the double-spin asymmetry from electron-$^{3}$He scattering in the threshold region of two- and three-body breakup of $^{3}$He was performed at Jefferson Lab, for Q$^{2}$ values of 0.1 and 0.2 (GeV/$c$)$^{2}$. The results of this measurement serve as a stringent test of our understanding of few-body systems. When compared with calculations from plane wave impulse approximation and Faddeev theory, we found that the Faddeev calculations, which use modern nuclear potentials and prescriptions for meson-exchange currents, demonstrate an overall good agreement with data.

The differential cross sections for the $^2$H(p,pp)n reaction have been measured for 84 angular configurations of the outgoing protons in the range of polar angles from 13 to 33 degrees with a proton beam of 108 MeV. Data have been collected in the Cyclotron Center Bronowice (CCB) at the Institute of Nuclear Physics PAS in Cracow, using the BINA detector setup. Analysis leading to determination of the breakup cross section values is described. Absolute normalization is obtained by normalization to the simultaneously measured $^2$H(p,d)p scattering events.Experimental results are compared to the state-of-the-art theoretical calculations. Global analysis shows significant influence of the Coulomb interaction and small effects of three nucleon force in the studied phase space region.

We discuss two approaches which, by applying the screening method, permit one to include the long range proton-proton (pp) Coulomb force in proton-deuteron (pd) momentum-space scattering calculations. In the first one, based on Alt-Grassberger-Sandhas (AGS) equation, presented in Phys. Rev. C{\bf{71}}, 054005 (2005) and {\bf{73}}, 057001 (2006), one needs to renormalize elastic scattering amplitude before calculating observables. In the second treatment, proposed by us in Eur. Phys. Journal A {\bf{41}}, 369 (2009), {\bf{41}}, 385 (2009), and arXiv:2310.03433 [nucl.th], this renormalization is avoided. For the proton induced deuteron breakup reaction both approaches require renormalization of the corresponding transition amplitudes. We derive the basic equations underlying both methods under the assumption that all contributing partial wave states are included and explain why in our approach renormalization of the elastic scattering amplitude is superfluous. We show that in order to take into account in the screening limit all partial waves it is required that four additional terms, based on the 3-dimensional and partial-wave projected pp Coulomb t-matrices, identical for both approaches, must appear in transition amplitudes. We investigate importance of these terms for elastic pd scattering below the breakup threshold., Comment: 17 pages, 6 figures

We propose a simplified approach to incorporate the long-range proton-proton (pp) Coulomb force in the three-nucleon (3N) scattering calculations, based on exact formulation presented in Eur. Phys. Journal A {\bf{41}}, 369 (2009) and {\bf{41}}, 385 (2009). It permits us to get elastic proton-deuteron (pd) scattering and breakup observables relatively simply by performing standard Faddeev calculations as known for the neutron-deuteron (nd) system. The basic ingredient in that approach is a 3-dimensional screened pp Coulomb t-matrix obtained by numerical solution of the 3-dimensional Lippmann-Schwinger (LS) equation. Based on this t-matrix pure Coulomb transition terms contributing to elastic scattering and breakup are calculated without any need for partial wave decomposition. For elastic scattering such a term removes the Rutherford amplitude for point deuteron proton-deuteron (pd) scattering. For breakup it has never been applied in spite of the fact that its contributions could become important in some regions of the breakup phase space. We demonstrate numerically that the pd elastic observables can be determined directly from the resulting 3N amplitudes without any renormalization, simply by increasing the screening radius in order to reach the existing screening limit. However,for pd breakup the renormalization of the contributing on-shell amplitudes is required. We apply our approach in a wide energy range of the incoming proton for pd elastic scattering as well as for pd breakup reaction., Comment: 50 pages, 21 figures

We build a relativistic model to perform calculations of exclusive, semi-exclusive and inclusive unpolarized cross sections and various polarization observables in electron and neutrino scattering experiments with deuteron targets. The strong interaction dynamics is defined by an explicit dynamical unitary representation of the Poincar\'e group, where representations of space translations and rotations in the interacting and non-interacting representations are the same. The Argonne V18 potential is used to construct a relativistic nucleon-nucleon interaction reproducing the experimental deuteron binding energy and nucleon-nucleon scattering observables. Our formalism does not include the pion production channel and neglects two-body contributions in the electromagnetic as well as in the weak nuclear current operator. We show that it is applicable to processes at kinematics, where the internal two-nucleon energy remains below the pion production threshold but the magnitude of the three-momentum transfer extends at least to several GeV., Comment: 59 pages, 29 figures

The pi- + 2H -> n + n, pi- + 3H -> n + n + n, pi- + 3He -> n + d and pi- + 3He -> p + n + n capture reactions from the lowest atomic orbitals are studied under full inclusion of final state interactions. Our results are obtained with the single-nucleon and two-nucleon transition operators derived at leading order in chiral effective field theory. The initial and final three-nucleon states are calculated with the chiral nucleon-nucleon SMS potential up to N4LO+ augmented by the consistently regularized chiral N2LO three-nucleon potential. We found that absorption rates depend strongly on the nuclear pion absorption operator used, and its two-body parts change the rates by a few orders of magnitude. The final state interactions between nucleons generated by the two-nucleon forces are also important, while the three-nucleon interaction plays a visible role only in the pi- + 3He -> n + d reaction. Our absorption rate for the pi- + 2H -> n + n process is in good agreement with the experimental data from the hadronic ground-state broadening in pionic deuterium. The capture rates on 3He are also generally consistent with the spectroscopic data within error bars, though our central values are found to be systematically below the data. We show that for the three-body breakup processes the dominant contributions to the absorption rates arise from the quasi-free scattering and final-state interaction kinematical configurations., Comment: 30 pages, 11 figures; authors' list corrected

We present a comprehensive investigation of few-nucleon systems as well as light and medium-mass nuclei up to $A=48$ using the current Low Energy Nuclear Physics International Collaboration two-nucleon interactions in combination with the third-order (N$^2$LO) three-nucleon forces. To address the systematic overbinding of nuclei starting from $A \sim 10$ found in our earlier study utilizing the N$^2$LO two- and three-nucleon forces, we take into account higher-order corrections to the two-nucleon potentials up through fifth order in chiral effective field theory. The resulting Hamiltonian can be completely determined using the $A=3$ binding energies and selected nucleon-deuteron cross sections as input. It is then shown to predict other nucleon-deuteron scattering observables and spectra of light $p$-shell nuclei, for which a detailed correlated truncation error analysis is performed, in agreement with experimental data. Moreover, the predicted ground state energies of nuclei in the oxygen isotopic chain from $^{14}$O to $^{26}$O as well as $^{40}$Ca and $^{48}$Ca show a remarkably good agreement with experimental values, given that the Hamiltonian is fixed completely from the $A \leq 3$ data, once the fourth-order (N$^3$LO) corrections to the two-nucleon interactions are taken into account. On the other hand, the charge radii are found to be underpredicted by $\sim 10\%$ for the oxygen isotopes and by almost $20\%$ for $^{40}$Ca and $^{48}$Ca., Comment: 20 pages, 14 figures, 8 tables

We employ two models of the nucleon-nucleon force: the OPE-Gaussian as well as the chiral N4LO and N4LO+ interactions with semilocal regularization in momentum space to study correlations among two-nucleon and three-nucleon elastic scattering observables. These models contain a number of free parameters whose values and covariance matrices are evaluated from a fit to the two-nucleon data. Such detailed knowledge of parameters allows us to create, using various sets of statistically generated parameters, numerous versions of these potentials and next apply them to two- and three-nucleon scattering to make predictions of various observables at the reaction energies up to 200 MeV. This permits a systematic analysis of correlations among two-nucleon and three-nucleon observables, basing on a relatively big sample of predictions. We found that most observables in neutron-proton and neutron-deuteron systems are uncorrelated, but there are exceptions revealing strong correlations, which depend on the reaction energy and scattering angle. This information may be useful for precise fixing free parameters of two-nucleon and three-nucleon forces and for understanding dependencies and correlations between potential parameters and observables., Comment: 20 pages, 8 pdf figures

We investigate the importance of the three-nucleon (3N) force contact terms in elastic nucleon-deuteron (Nd) scattering by applying the N$^4$LO$^+$ chiral semi-local momentum space (SMS) regularized nucleon-nucleon (NN) chiral potential supplemented by N$^2$LO and all subleading N$^4$LO three-nucleon force (3NF) contact terms. Strength parameters of the contact terms were obtained by least squares fitting of theoretical predictions to cross section and analyzing powers data at three energies of the impinging nucleon. Although the N$^3$LO contributions to the 3N force were completely neglected, the results calculated with the contact terms multiplied by the fitted strength parameters yield an improved description of the elastic Nd scattering observables in a wide range of incoming nucleon energies below the pion production threshold., Comment: 43 pages, 13 figures