Abstract
We discuss the current status of chiral effective field theory in the three-nucleon sector and present selected results for nucleon–deuteron scattering observables based on semilocal momentum-space-regularized chiral two-nucleon potentials together with consistently regularized three-nucleon forces up to third chiral order. Using a Bayesian model for estimating truncation errors, the obtained results are found to provide a good description of the experimental data. We confirm our earlier findings that a high-precision description of nucleon–deuteron scattering data below pion production threshold will require the theory to be pushed to fifth chiral order. This conclusion is substantiated by an exploratory study of selected short-range contributions to the three-nucleon force at that order, which, as expected, are found to have significant effects on polarization observables at intermediate and high energies. We also outline the challenges that will need to be addressed in order to push the chiral expansion of three-nucleon scattering observables to higher orders.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The past few years have seen remarkable advances towards pushing the precision frontier of chiral effective field theory (EFT) in the two-nucleon sector; see Refs. [1,2,3] for review articles. The fifth-order (\(\hbox {N}^4\hbox {LO}\)) contributions to the nucleon-nucleon (NN) force have been worked out in Ref. [4] and a new generation of accurate and precise chiral EFT NN potentials up through \(\hbox {N}^4\hbox {LO}\) has been developed in Refs. [5,6,7,8], see also Refs. [9,10,11] for related developments. The \(\hbox {N}^4\hbox {LO}\) interactions of Refs. [7, 8] utilize the values of the pion–nucleon (\(\pi \hbox {N}\)) low-energy constants (LECs) determined from matching chiral perturbation theory to the Roy–Steiner-equation analysis of \(\pi \hbox {N}\) scattering at the subthreshold point [12], but differ substantially by the regularization procedure. The potentials developed by our group in Ref. [7] build upon our earlier studies [5, 6] and employ a local regulator for pion-exchange contributions which, per construction, maintains the analytic structure of the long-range interaction. Differently to the nonlocally regularized potentials of Ref. [8], the interactions constructed in Ref. [7] do not produce long-range artifacts at any finite-order of expansion in inverse powers of the momentum-space cutoff \(\varLambda \). The resulting semilocalFootnote 1 momentum-space regularized (SMS) potentials of Ref. [7] are currently the most precise chiral EFT interactions and provide, at the highest order \(\hbox {N}^4\hbox {LO}^+\)Footnote 2, a nearly perfect and \(\varLambda \)-independent (within the considered cutoff range) description of neutron–proton and proton–proton scattering data below pion production threshold from the self-consistent 2013 Granada data base [13]. The achieved description of the scattering data is comparable to or even better than that based on the so-called high-precision semi-phenomenological potentials like the AV18 [14], CD Bonn [15] and Nijm I, II [16] models.
In spite of this exciting progress in the NN sector, applications of chiral EFT to three- and more-nucleon systems and in processes involving external sources are, with very few exceptions, still limited to the next-to-next-to-leading order (\(\hbox {N}^2\hbox {LO}\)) in the chiral expansion. What makes high-accuracy calculations of such systems/reactions beyond \(\hbox {N}^2\hbox {LO}\) so difficult? The main technical and conceptual difficulties are related to the treatment of many-body forces and exchange current operators. Three- (3NF) and four-nucleon forces (4NF) start to contribute at third (\(\hbox {N}^2\hbox {LO}\)) and fourth (\(\hbox {N}^3\hbox {LO}\)) orders of the chiral expansion, respectively, while the first contributions to the exchange electroweak currents appear already at second order (NLO) relative to the dominant single-nucleon terms. Over the past decade, we have worked out off-shell-consistent expressions for the 3NF [17, 18], 4NF [19, 20] and electroweak charge and current operators [21,22,23,24] completely up through \(\hbox {N}^3\hbox {LO}\) using dimensional regularization (DR) to compute pion loop contributions; see also Refs. [25,26,27,28] for a similar work on electroweak currents by the JLab-Pisa group and Ref. [29] for a related discussion. Furthermore, selected \(\hbox {N}^4\hbox {LO}\) contributions to the 3NF have been worked out in Refs. [30,31,32], and the longest-range part of the 3NF was also analyzed in the framework involving \(\varDelta (1232)\) degrees of freedom. On the technical side, the implementation of the 3NFs and exchange currents in few-body calculations requires their partial-wave decomposition. This nontrivial task can nowadays be accomplished numerically for a general 3NF specified in momentum space as described in Refs. [33, 34]. The main conceptual challenge that still needs to be addressed concerns a consistent regularization of many-body forces and exchange operators. While this issue is irrelevant in the NN sector, a naive regularization of many-body interactions and exchange currents based on the expressions derived using DR violates chiral symmetry and leads to inconsistent/wrong results at the one-loop level (i.e. at \(\hbox {N}^3\hbox {LO}\)) and beyond [35, 36]. A possible solution could be provided by using higher-derivative regularization instead of DR, which has to be chosen in a way compatible with the regularization scheme of Ref. [37] in the NN sector. Work along these lines is in progress.
In this paper we update our recent study [38] based on the semilocal coordinate-space regularized (SCS) NN forces of Refs. [5, 6] and analyze nucleon–deuteron (Nd) scattering observables using the new SMS chiral NN potentials of Ref. [7] in combination with the SMS 3NF at \(\hbox {N}^2\hbox {LO}\). We also refine our previous estimations of truncation uncertainties by employing a Bayesian approach instead of the algorithm proposed in Ref. [5] and used in our earlier studies [38,39,40] in the three-nucleon sector. Last but not least, we explore the role of selected \(\hbox {N}^4\hbox {LO}\) short-range 3NF terms in Nd scattering.
Our paper is organized as follows. In Sect. 2, we specify our Bayesian model for truncation errors by following the approach proposed in Ref. [41]. Our results for Nd scattering observables at \(\hbox {N}^2\hbox {LO}\) are presented in Sect. 3, while the role of selected \(\hbox {N}^4\hbox {LO}\) 3NF operators is discussed in Sect. 4. The main results of our study are summarized in Sect. 5.
2 Bayesian model for truncation errors
A reliable estimation of theoretical uncertainties is an essential ingredient of any systematic approach such as chiral EFT. Cutoff variation offers one possibility to quantify the impact of contributions beyond the truncation level. However, in the few- and many-nucleon sectors, the available cutoff range is often rather limited. Furthermore, cutoff variation does not allow one to probe the impact of neglected long-range interactions. In Ref. [5], a more reliable, universally applicable algorithm for estimating truncation errors using the available information on the chiral expansion for any observable of interest without relying on cutoff variation was proposed. Here and in what follows, this algorithm will be referred to as the EKM approach. For applications of the EKM method to a broad range of low-energy reactions in the single-baryon and few-/many-nucleon sectors, see Refs. [42,43,44,45] and [6, 38,39,40, 46,47,48], respectively. Being very simple and easy to implement, the EKM approach does, however, not directly provide a statistical interpretation of the estimated uncertainties. In Refs. [41, 49], a general Bayesian approach to calculating the posterior probability distribution for truncation errors in chiral EFT was developed. The EKM approach was then shown to essentially correspond to a particular choice of prior probability distribution for the coefficients in the chiral expansion of observables. Using the chiral NN potentials from Ref. [5], the EKM error estimations for the neutron–proton (np) total cross section at selected energies were found in the Bayesian approach of Ref. [49] to be consistent with \(68\%\) degree-of-belief (DoB) intervals.
Throughout this paper, we employ a slightly modified version of the Bayesian model from Ref. [49]. Specifically, consider a two-nucleon scattering observable X(p) with p referring to the center-of-mass (CM) momentum. Calculating X(p) using chiral EFT potentials at various orders \(Q^i\), \(i=0, 2, 3, 4, \ldots \), (but for a fixed cutoff value) yields the corresponding predictions \(X^{(i)} (p)\), and the chiral expansion of X(p) can be written in the form
with \(\varDelta X^{(2)} := X^{(2)} - X^{(0)}\) and \(\varDelta X^{(i)} := X^{(i)} - X^{(i-1)}\) for \(i > 2\). The second equality serves as a definition of the dimensionless expansion coefficients \(c_i\). The reference value \(X_{\mathrm{ref}}\), which sets the overall scale, will be defined below. Here and in what follows, the expansion parameter of chiral EFT is assumed to have the form [5, 6, 41, 49]
where \(\varLambda _\mathrm{b}\) is the breakdown scale of the chiral expansion. The quantity \(M_\pi ^{\mathrm{eff}}\) serves to model the expansion of NN observables around the chiral limit. In Refs. [5, 6, 41, 49], this scale was set to the physical pion mass \(M_\pi \). However, as pointed out in Ref. [50], the error plots in Ref. [5] indicate that the transition between the two expansion regimes in the NN sector actually appears at a scale \(M_\pi ^{\mathrm{eff}}\) higher than \(M_\pi \). On the other hand, both Bayesian model parameters \(\varLambda _\mathrm{b}\) and \(M_\pi ^{\mathrm{eff}}\) can be determined/tuned empirically by calculating the success rates for a given set of observables and/or energies. In particular, Ref. [49] confirmed the EKM estimation \(\varLambda _\mathrm{b} \sim 600~\hbox {MeV}\) based on the results for the total np cross section and using the SCS potentials of Ref. [5], but it also found somewhat larger values of \(\varLambda _\mathrm{b}\) to be statistically consistent; see also a related discussion in [51]. A similar empirical analysis was performed in Ref. [50] for both \(\varLambda _\mathrm{b}\) and \(M_\pi ^{\mathrm{eff}}\) yielding the values of \(\varLambda _\mathrm{b} \sim 650 \ldots 700~\hbox {MeV}\) and \(M_\pi ^{\mathrm{eff}} \sim 200~\hbox {MeV}\).
Suppose the results for the observable X(p) are available up through the order \(X^{(k)}\), \(k \ge 2\). The goal is then to estimate the truncation error \(\delta X^{(k)} \equiv \sum _{i > k} \varDelta X^{(i)}\) resulting from neglecting the unknown higher-order contributions, i.e. to compute the posterior probability distribution function (pdf) for \(\delta X^{(k)}\) given the explicit knowledge of \(\{ X^{(0)}, \, X^{(2)}, \, \ldots , X^{(k)} \}\). The Bayesian model of Ref. [41] uses the leading-order (LO) result for X(p) to set the overall scale
in order to define the dimensionless expansion coefficients \(c_i\) with \(c_0=1\) in Eq. (1). As we will argue below, this approach may, for certain choices of the prior pdf, be too restrictive in the kinematical regions near the points where the LO contribution vanishes e.g. by changing the sign. Such situations are not uncommon if one looks at observables which depend on continuously varying parameters such as energy or scattering angle, see also a discussion in Ref. [51]. In such circumstances it is advantageous to set the overall scale from the next-to-leading order (NLO) contribution \(\varDelta X^{(2)}\) via \(X_{\mathrm{ref}} = \varDelta X^{(2)}/Q^2\) in order to avoid underestimating \(X_{\mathrm{ref}}\). To have an approach applicable when both \(X^{(0)}\) and/or \(\varDelta X^{(2)}\) are accidentally small, we set the scale \(X_{\mathrm{ref}}\) via
for \(k = 2\) and
for \(k \ge 3\). Let \(c_m = 1\), \(m \in \{0, 2, 3\}\), be the coefficient used to define the overall scale. Assuming that the remaining coefficients \(c_i\) are distributed according to some common pdf \(\mathrm{pr} (c_i | {\bar{c}} )\) with a hyperparameter \({\bar{c}}\) and performing marginalization over h chiral orders \(k+1, \ldots , k+h\), which are assumed to dominate the truncation error, the probability distribution for the dimensionless residual \(\varDelta _k \equiv \sum _{n=k+1}^\infty c_n Q^n \simeq \sum _{n=k+1}^{k+h} c_n Q^n\) to take a value \(\varDelta _k = \varDelta \), given the knowledge of \(\{c_{i \le k} \}\), is given by [41]
where the set A is defined as \(A = \{n \in {\mathbb {N}}_0 \, | \, n \le k \, \wedge \, n \ne 1 \, \wedge \, n \ne m \}\) and
Here and in what follows, we employ the Gaussian prior of “set C” from Ref. [41], namely
and assume a log-uniform probability distribution, [52]
The nice feature of this prior is that all integrations in Eq. (6) can be carried out analytically. For the sake of completeness, we give in Appendix A the corresponding expressions for \(\mathrm{pr}_h ( \varDelta | \{ c_{i \le k} \})\) from Ref. [41]. The posterior pdf is an even function of \(\varDelta \), and for any given DoB interval, the corresponding value of the residual \(\varDelta _k\) and the resulting truncation error \(\delta X^{(k)} = X_{\mathrm{ref}} \varDelta _k\) can be readily obtained by numerically integrating \(\mathrm{pr}_h ( \varDelta | \{ c_{i \le k} \})\) over \(\varDelta \).
As a first application, we employ the Bayesian model of Ref. [41] and set the overall scale \(X_{\mathrm{ref}}\) solely from the corresponding LO contribution; see Eq. (3). In the bottom row of Figs. 1 and 2, we show the \(68\%\) and \(95\%\) DoB intervals for selected np scattering observables at \(E_{\mathrm{lab}}= 143~\hbox {MeV}\) for the non-informative prior \(C_\epsilon \) [49] with \({\bar{c}}_< = \epsilon \) and \({\bar{c}}_> = 1/\epsilon \), which makes no assumption as regards either the maximum or the minimum size of \({\bar{c}}\) by taking the limit \(\epsilon \rightarrow 0\); see Eq. (23). Here, the breakdown scale was set to \(\varLambda _\mathrm{b} = 700~\hbox {MeV}\), and the resulting Bayesian model is referred to as \(C_\epsilon ^{700}\).
Estimated theoretical uncertainty for the chiral EFT results for np differential cross section \(\mathrm{d} \sigma /\mathrm{d} \varOmega \) (left panel) and polarization transfer coefficient \(D_t\) at laboratory energy of \(E_{\mathrm{lab}} = 143~\hbox {MeV}\). The light- (dark-) shaded yellow, green, blue and red bands of decreasing width depict \(95\%\) (\(68\%\)) DoB intervals at NLO, \(\hbox {N}^2\hbox {LO}\), \(\hbox {N}^3\hbox {LO}\) and \(\hbox {N}^4\hbox {LO}\), respectively. Dashed lines show the LO predictions. Open circles refer to the results of the Nijmegen partial-wave analysis [53]. Data for the cross section are at \(E_{\mathrm{lab}} = 142.8\,\hbox {MeV}\) and taken from [54]. The first, second, third and fourth rows correspond to the Bayesian models \(C_{0.5-10}^{650}\), \({\tilde{C}}_{0.5-10}^{650}\), \({\bar{C}}_{0.5-10}^{650}\) and \(C_{\epsilon }^{700}\). All results shown are based on the SMS NN potentials using the cutoff of \(\varLambda = 450~\hbox {MeV}\)
Being non-informative, the prior set \(C_\epsilon \) is generally expected to yield conservative estimates for truncation errors. However, setting \({\bar{c}}_< \rightarrow 0\) yields a \(\delta \)-function-like posterior \(\mathrm{pr}_h^{C_\epsilon } ( \varDelta | \{ c_{i \le k}\})\) for \({\varvec{c}}_k^2 \rightarrow 0\) as can be seen from Eq. (23), i.e. this model fails to provide an adequate estimation of the truncation error if the corrections \(\varDelta ^{(2)}, \; \ldots , \; \varDelta ^{(k)}\) happen to be accidentally small. For the examples shown in Figs. 1 and 2, this is the case at NLO for the differential cross section \(d \sigma /d \varOmega \) around \(\theta _{\mathrm{CM}} \sim 85^\circ \), for the polarization transfer \(D_t\) around \(\theta _{\mathrm{CM}} \sim 140^\circ \), for the spin-correlation coefficients \(C_{kp}\) at \(\theta _{\mathrm{CM}} \sim 65^\circ \) and for the coefficient \(C_{kk}\) at \(\theta _{\mathrm{CM}} \sim 25^\circ \), and \(\theta _{\mathrm{CM}} \sim 115^\circ \). In all these cases, the corresponding functions \(\varDelta X^{(2)} (\theta _{\mathrm{CM}} )\) change their sign. To circumvent the problem with the underestimation of the truncation error in such kinematical regions, the authors of Ref. [41] suggested to use a more informative (but not too restrictive) prior set \(C_{0.25-10}\) corresponding to \({\bar{c}}_< = 0.25\) and \({\bar{c}}_> = 10\). Here and in what follows, we make the choice \({\bar{c}}_< = 0.5\), which we found to be more efficient in resolving the above-mentioned issue while still sufficiently general. As shown in the upper row of Figs. 1 and 2, the prior set \(C_{0.5-10}^{650}\) indeed yields reasonable estimates of the truncation errors for \(\mathrm{d} \sigma /\mathrm{d} \varOmega \). However, the more informative prior with \({\bar{c}}_< \ne 0\) suffers from another issue: it yields a vanishingly small truncation error at all orders in the cases when \(X^{(0)}\) happens to be accidentally small. This is the case for \(D_t\) at \(\theta _{\mathrm{CM}} \sim 100^\circ \) and for \(C_{kk}\) at \(\theta _{\mathrm{CM}} \sim 10^\circ \), \(\theta _{\mathrm{CM}} \sim 75^\circ \) and \(\theta _{\mathrm{CM}} \sim 180^\circ \); see the plots in the upper row of Figs. 1 and 2. The most extreme situation is observed for the coefficient \(C_{kp}\), for which the LO contribution appears to be small for all scattering angles. The problem can be traced back to the misidentification of the overall scale by Eq. (3) in such accidental cases. Writing \(X^{(0)}\) as \( X^{(0)}= \alpha {{\tilde{X}}}^{(0)} \) with \(\alpha \rightarrow 0\) being a dimensionless parameter, one finds \(\varDelta _k \sim \alpha ^{-1} \) for the prior set \(C_\epsilon \) while \(\varDelta _k \sim \alpha ^0\) leading to \(\delta X^{(k)} = \alpha {\tilde{X}}^{(0)} \varDelta _k \sim \alpha \) for the prior set \(C_{{\bar{c}}_< - {\bar{c}}_>}\) with \({\bar{c}}_> < \infty \). The problem can be easily fixed by replacing Eq. (3) with Eq. (4) as shown in the second row of Figs. 1 and 2. Here and in what follows, the resulting Bayesian model is referred to as \({\tilde{C}}\). However, while highly unlikely, it is still possible that both the LO contribution and the NLO correction are simultaneously accidentally small. For the considered observables, this happens for \(D_t\) at backward and for \(C_{kk}\) at forward angles. To prevent underestimating the truncation errors in such rare cases, we replace, for \(k \ge 3\), Eq. (4) with Eq. (5) and refer to the resulting Bayesian model as \({\bar{C}}\). This model is found to yield robust results for NN scattering observables in all kinematical regions; see e.g. the third row in Figs. 1 and 2 and will be employed in nucleon–deuteron scattering calculations considered in the next sections.
Same as Fig. 1 but for the spin-correlation parameters \(C_{kp}\) and \(C_{kk}\)
It is important to emphasize that the considered examples have been selected to visualize the possible issues in certain accidental situations, and that the differences between the considered Bayesian models are fairly minor for most observables. This is exemplified in Fig. 3, where the truncation errors for the total cross section are shown for a variety of considered models including the original EKM approach (with \(M_\pi ^{\mathrm{eff}} = M_\pi \) and \(\varLambda _\mathrm{b} = 600~\hbox {MeV}\)).
Neutron–proton cross sections at \(E_{\mathrm{lab}} = 50~\hbox {MeV}\) (left panel) and \(E_{\mathrm{lab}} = 150~\hbox {MeV}\) (right panel) at various orders of the chiral expansion using the SMS potentials of Ref. [7] with the cutoff \(\varLambda = 450~\hbox {MeV}\). The smaller (red) error bars correspond to \(68\%\) DoB intervals while the larger (blue) ones indicate \(95\%\) DoB intervals for a variety of Bayesian models as described in the text. The gray shaded bands show the uncertainty estimates using the original EKM approach. Dashed lines show the result of the next-order calculation (\(\hbox {N}^2\hbox {LO}\) at NLO, \(\hbox {N}^3\hbox {LO}\) at \(\hbox {N}^2\hbox {LO}\), \(\hbox {N}^4\hbox {LO}\) at \(\hbox {N}^3\hbox {LO}\) and \(\hbox {N}^4\hbox {LO}^+\) at \(\hbox {N}^4\hbox {LO}\))
As pointed out in Ref. [49], the dependence on a particular Bayesian model and/or assumed prior set decreases with an increasing order, i.e. with the increasing amount of information about the actual pattern of the chiral expansion. Notice that differently to the EKM approach, the considered Bayesian models exploit only the information up to the order, at which the truncation error is estimated. The more conservative error estimations at \(E_{\mathrm{lab}} = 50~\hbox {MeV}\) with the models C and \({\bar{C}}\) as compared to the original EKM approach are mainly due to the larger value of \(M_\pi ^{\mathrm{eff}}\). For the considered set of the total cross section calculations, counting the success rate for the next-higher order result to lie within the estimated uncertainty as shown in Fig. 3 yields the values of \(62.5\ldots 75\%\) (\(100\%\)), which are statistically consistent with the DoB intervals of \(68\%\) (\(95\%\)).Footnote 3 We have verified that this conclusion also holds true for the larger set of energies considered in Ref. [50]. We refrain from performing similar statistical tests for angular distributions due to their correlated nature. This could be done using the method proposed in Ref. [51], which is based on Gaussian processes and encodes correlation structure of coefficients \(c_i (\theta _\mathrm{CM} )\).
3 Nucleon–deuteron scattering at \(\hbox {N}^2\hbox {LO}\) using SMS nuclear potentials
We now turn to the main topic of our study and consider nucleon–deuteron scattering in chiral EFT. The \(\hbox {N}^2\hbox {LO}\) three-nucleon force is given by [55, 56]
where the subscripts refer to the nucleon labels, \({\mathbf {q}}_{i} = {\mathbf {p}}_i \, ' - {\mathbf {p}}_i\) with \({\mathbf {p}}_i \, '\) and \({\mathbf {p}}_i\) being the final and initial momenta of the nucleon i and \({\mathbf {\sigma }}_i\) (\({\varvec{\tau }}_i\)) are the Pauli spin (isospin) matrices. Further, \(c_i\), D and E denote the corresponding low-energy constants (LECs), while \(g_A\) and \(F_\pi \) refer to the nucleon axial coupling and pion decay constant, respectively. It is customary to express the LECs D and E in terms of the corresponding dimensionless coefficients via
where, following [38, 56], we use \(\varLambda _\chi = 700~\hbox {MeV}\). For the LECs \(c_i\), we employ the central values from the Roy–Steiner-equation analysis of Ref. [12] at the corresponding chiral order, namely \(c_1= -0.74~\hbox {GeV}^{-1}\), \(c_3 = -3.61~\hbox {GeV}^{-1}\) and \(c_4 = 2.44~\hbox {GeV}^{-1}\). The same values are used in the SMS NN potentials of Ref. [7] at \(\hbox {N}^2\hbox {LO}\). Differently to Ref. [38], we employ the same semilocal momentum-space regulator for the 3NF as in the NN potentials of Ref. [7] by replacing the pion propagators via
For the D-term, the contact interaction between nucleons 1 and 2 is, in addition, regularized by multiplying the matrix elements with a nonlocal Gaussian regulator \(\exp (- ({\mathbf {p}}_{12}^2 + {{\mathbf {p}}_{12}'} ^2)/\varLambda ^2 )\), where \({\mathbf {p}}_{12} = ({\mathbf {p}}_1 - {\mathbf {p}}_2 )/2\), \({{\mathbf {p}}}_{12}^{\, \prime } = ({\mathbf {p}}_1^{\, \prime } - {\mathbf {p}}_2^{\, \prime } )/2\). For the contact interaction proportional to the LEC E, we apply a nonlocal Gaussian regulator in momentum space,
where \({\mathbf {k}}_3 = 2 ({\mathbf {p}}_3 - ({\mathbf {p}}_1 + {\mathbf {p}}_2)/2)/3\) and \({\mathbf {k}}_3^{\, \prime } = 2 ({\mathbf {p}}_3^{\, \prime } - ({\mathbf {p}}_1^{\, \prime } + {\mathbf {p}}_2^{\, \prime })/2)/3\) are the corresponding Jacobi momenta.
It is important to emphasize that the SMS NN potentials of Ref. [7] employ additional (local) short-range subtractions to ensure that the coordinate-space expressions of the regularized pion-exchange contributions and derivatives thereof vanish at the origin. This convention ensures that regularized pion-exchange contributions contain only long-range pieces. On the other hand, using the regulator in Eq. (12), the resulting TPE contributions still contain admixtures of short-range terms of the D- and E-types.
Following Ref. [38], we determine the LECs \(c_D\) and \(c_E\) from the \(^3\)H binding energy and the nucleon–deuteron differential cross section minimum at \(E_{\mathrm{lab}}^N = 70~\hbox {MeV}\). Specifically, we fit to the experimental data in the range of \(\theta _{\mathrm{CM}} = 107^\circ \ldots 141^\circ \). The resulting \(c_D\)- and \(c_E\)-values are strongly dependent on the cutoff \(\varLambda \) and change from \(c_D = 8.9\) (\(c_E = 1.15\)) for \(\varLambda = 400~\hbox {MeV}\) to \(c_D = -5.4\) (\(c_E = -0.25\)) for \(\varLambda = 550~\hbox {MeV}\). Such a strong cutoff dependence has to be expected given that these LECs refer to bare quantities. The calculated observables, on the other hand, show only a weak residual cutoff dependence, consistent with the estimated truncation uncertainty.
In Figs. 4, 5, 6 and 7 we show the results for the differential cross section and selected polarization observables in elastic nucleon–deuteron scattering at NLO and \(\hbox {N}^2\hbox {LO}\) at energies of \(E_{\mathrm{lab}} = 10~\hbox {MeV}\), \(E_{\mathrm{lab}} = 70~\hbox {MeV}\) and \(E_{\mathrm{lab}} = 135~\hbox {MeV}\) for the cutoff \(\varLambda = 500~\hbox {MeV}\), along with the truncation errors corresponding to the DoB intervals of \(68\%\) and \(95\%\).
Predictions for the differential cross section, nucleon and deuteron analyzing powers \(A_y^n\) and \(A_y^d\) as well as deuteron tensor analyzing powers \(A_{yy}\), \(A_{xz}\) and \(A_{xx}\) in elastic nucleon–deuteron scattering at laboratory energy of \(E_\mathrm{lab}^N = 10~\hbox {MeV}\) at NLO (yellow bands) and \(\hbox {N}^2\hbox {LO}\) (green bands) based on the SMS NN potentials of Ref. [7] for \(\varLambda = 500~\hbox {MeV}\). The light- (dark-) shaded bands indicate \(95\%\) (\(68\%\)) DoB intervals using the Bayesian model \({\bar{C}}_{0.5-10}^{650}\). The dotted (dashed) lines show the results based on the CD Bonn NN potential [15] (CD Bonn NN potential in combination with the Tucson-Melbourne 3NF [57]). Open circles are neutron–deuteron from Ref. [58] and proton–deuteron data from Ref. [59,60,61], corrected for the Coulomb effects; see Ref. [56] for details
Results for the differential cross section, nucleon and deuteron analyzing powers \(A_y^n\) and \(A_y^d\) as well as deuteron tensor analyzing powers \(A_{yy}\), \(A_{xz}\) and \(A_{xx}\) in elastic nucleon–deuteron scattering at laboratory energy of \(E_\mathrm{lab}^N = 70~\hbox {MeV}\) at NLO (yellow bands) and \(\hbox {N}^2\hbox {LO}\) (green bands) based on the SMS NN potentials of Ref. [7] for \(\varLambda = 500~\hbox {MeV}\). Red dashed lines show the \(\hbox {N}^2\hbox {LO}\) results for the cutoff values of \(\varLambda = 400\), 450, 500 and \(550~\hbox {MeV}\) (the lines with a shorter dash length correspond to smaller cutoff values). Open circles are proton–deuteron data from Ref. [62]. For remaining notation, see Fig. 4
Predictions for the differential cross section, nucleon and deuteron analyzing powers \(A_y^n\) and \(A_y^d\) as well as deuteron tensor analyzing powers \(A_{yy}\), \(A_{xz}\) and \(A_{xx}\) in elastic nucleon–deuteron scattering at laboratory energy of \(E_\mathrm{lab}^N = 135~\hbox {MeV}\) at NLO (yellow bands) and \(\hbox {N}^2\hbox {LO}\) (green bands) based on the SMS NN potentials of Ref. [7] for \(\varLambda = 500~\hbox {MeV}\). Open circles are proton–deuteron data from Ref. [62]. For remaining notation, see Fig. 4
Predictions for polarization transfer coefficients \(K_{xx}^y\), \(K_{y}^y\), \(K_{yy}^y\), \(K_{xz}^y\), \(K_{xx}^y-K_{yy}^y\) and the induced polarization \(P^y\) in elastic nucleon–deuteron scattering at laboratory energy of \(E_{\mathrm{lab}}^N = 135~\hbox {MeV}\) at NLO (yellow bands) and \(\hbox {N}^2\hbox {LO}\) (green bands) based on the SMS NN potentials of Ref. [7] for \(\varLambda = 500~\hbox {MeV}\). Open circles are proton–deuteron data from Ref. [63]. For remaining notation, see Fig. 4
Notice that the truncation errors are symmetric, and the actual results of our calculation lie in the middle of the corresponding error bands.
To estimate the truncation uncertainty we have used the Bayesian model \({\bar{C}}_{0.5-10}^{650}\) introduced in Sect. 2. For scattering reactions involving three and more nucleons, we, however, also need to specify the pertinent momentum p scale that sets the expansion parameter Q in Eq. (2). Consider nucleon-nucleus scattering and let \(E_{\mathrm{lab}}^N\) be nucleon beam energy in the laboratory system. Neglecting the neutron–proton mass difference and the binding energy of the target nucleus, which is assumed to consist of A nucleons, the CM momentum is related to \(E_{\mathrm{lab}}^N\) via
where m is the nucleon mass and “\(\simeq \)” refers to the nonrelativistic approximation. Identifying the scale \(p \equiv | {\mathbf {p}}|\) in Eq. (2) with \(p_{\mathrm{CM}} \equiv | {\mathbf {p}}_{\mathrm{CM}} |\) results in A-dependent values of the expansion parameter Q corresponding to the same excess energy. For example, the pion production threshold in the NN (Nd) system with \(E_\mathrm{lab}^N \sim 290~\hbox {MeV}\) (\(E_{\mathrm{lab}}^N \sim 215~\hbox {MeV}\)) corresponds to \(p_{\mathrm{CM}} \sim 370~\hbox {MeV}\) (\(p_{\mathrm{CM}} \sim 425~\hbox {MeV}\)), leading to the expansion parameter of \(Q = 0.57\) (\(Q = 0.65\)). Alternatively, one can define the momentum scale p in terms of the Lorentz-invariant excess energy \(\sqrt{s} - \sqrt{s_0} =\sqrt{s} - (A + 1) m\) available in the \(A+1\)-nucleon system and define the momentum scale p via the relation
that ensures that p coincides with \(p_{\mathrm{CM}}\) in the NN system. Here, s is the usual Mandelstam variable. One can thus express the scale p in terms of \(E_{\mathrm{lab}}^N\) via
with \(s = m^2 (A + 1)^2 + 2 A m E_{\mathrm{lab}}^N\). In the nonrelativistic approximation, this relation simplifies to
The nonrelativistic approximation holds at a sub-percent level for the energy range considered in this study and we use the relation (17) to define the expansion parameter Q in Eq. (2). The breakdown scale \(\varLambda _\mathrm{b} = 650~\hbox {MeV}\) then corresponds to the excess energy of \(\sim 400~\hbox {MeV}\) independently of the number of nucleons A in the target nucleus. Notice that the employed model leads to less conservative error estimates for \(A>1 \) than the assignment of \(p = p_{\mathrm{CM}}\) in Eq. (2).
We now turn to the results for Nd scattering observables at \(E_\mathrm{lab}^N = 10 \ldots 135~\hbox {MeV}\) shown in Figs. 4, 5, 6 and 7. Except for the differential cross section at \(E_{\mathrm{lab}}^N = 70~\hbox {MeV}\) shown in the upper left panel of Fig. 5, the results at \(\hbox {N}^2\hbox {LO}\) can be regarded as parameter-free predictions. It is reassuring to see that the calculated observables are in a reasonably good agreement with the experimental data, which in most cases lie within the \(95\%\) DoB intervals. One should, however, keep in mind that the estimated truncation errors depend on the Bayesian model, assumed prior sets and the values of parameters \(M_\pi ^{\mathrm{eff}}\) and \(\varLambda _\mathrm{b}\). While model dependence of uncertainty estimates is expected to decrease at high chiral orders, it may still be significant at \(\hbox {N}^2\hbox {LO}\).
While we only show the uncertainty bands for the cutoff value of \(\varLambda = 500~\hbox {MeV}\), the results for different values of \(\varLambda \) are similar. To illustrate this point we plot in Fig. 5 for the intermediate energy of \(E_{\mathrm{lab}}^N = 70~\hbox {MeV}\) the \(\hbox {N}^2\hbox {LO}\) results for all available cutoff values in the range of \(\varLambda = 400 \ldots 550~\hbox {MeV}\). Notice that the residual cutoff dependence of the considered observables at NLO is similar to the one at \(\hbox {N}^2\hbox {LO}\) and is, in most cases, comparable with the \(\hbox {N}^2\hbox {LO}\)\(68\%\) DoB intervals. These results demonstrate that our calculations for different values of \(\varLambda \) are consistent with each other within errors. Notice further that the softest cutoff of \(\varLambda = 400~\hbox {MeV}\) shows the largest deviation from the bulk behavior and from the experimental data, which presumably points to the increasing amount of finite-regulator artifacts.
We also show in Figs. 4, 5, 6 and 7 the results based on the CD Bonn NN potential [15] with and without the Tuscon-Melbourne (TM) 3NF [57]. In particular, for the cases where the TM 3NF is known to provide sizable corrections such as e.g. for the differential cross section around its minimum and for the deuteron analyzing powers at the intermediate energies of \(E_{\mathrm{lab}}^N = 70\) and \(135~\hbox {MeV}\), our \(\hbox {N}^2\hbox {LO}\) results agree well with the CD Bonn NN plus TM 3NF calculations, while the predictions based on the CD Bonn NN force alone are often outside the \(65\%\) and sometimes even \(95\%\) DoB intervals. This should not come as a surprise given the known weak dependence of the Nd elastic scattering observables on the off-shell behavior of the NN potentials [64] and a similar structure of the TM and the leading chiral 3NF at \(\hbox {N}^2\hbox {LO}\) [65, 66] largely driven by the intermediate \(\varDelta \)(1232) excitation.
Similarly to our findings in Ref. [38] based on the SCS interactions of Refs. [5, 6], the nucleon and deuteron vector analyzing powers are also properly described (within errors) at the lowest considered energy of \(E_{\mathrm{lab}}^N = 10~\hbox {MeV}\) showing no evidence of the so-called \(A_y\)-puzzle [64] at this chiral order. The much larger truncation uncertainty for the vector analyzing powers at this low energy as compared with other observables indicates their strongly fine-tuned nature; see also Ref. [38] for a related discussion. We further emphasize that the approximate subtraction of the Coulomb effects from the proton–deuteron data at this energy may lead to sizable uncertainties.
In Refs. [39, 40], we have calculated Nd scattering observables based on the SCS NN chiral potentials of Refs. [5, 6] and estimated the truncation errors using the EKM approach. While these calculations are incomplete starting from \(\hbox {N}^2\hbox {LO}\) due to the missing 3NF, they have demonstrated that the expected accuracy of chiral EFT at high orders such as \(\hbox {N}^4\hbox {LO}\) should be substantially smaller than the observed discrepancies between state-of-the-art calculations and experimental data.
Results for the differential cross section, nucleon and deuteron analyzing powers \(A_y^n\) and \(A_y^d\) as well as deuteron tensor analyzing powers \(A_{yy}\), \(A_{xz}\) and \(A_{xx}\) in elastic nucleon–deuteron scattering at laboratory energy of \(E_\mathrm{lab}^N = 200~\hbox {MeV}\) based on the SMS NN potentials of Ref. [7] at \(\hbox {N}^3\hbox {LO}\) (blue shaded bands) and \(\hbox {N}^4\hbox {LO}^+\) (red shaded bands) combined with the 3NF at \(\hbox {N}^2\hbox {LO}\) using \(\varLambda = 500~\hbox {MeV}\). Blue (red) shaded bands show the expected truncation uncertainty for complete\(\hbox {N}^3\hbox {LO}\) (\(\hbox {N}^4\hbox {LO}\)) calculations and are obtained by multiplying the \(\hbox {N}^2\hbox {LO}\) truncation errors for the model \(\bar{C}_{0.5-10}^{650}\) with the expansion parameter \(Q \simeq 0.55\) (\(Q^2 \simeq 0.3\)). Open circles are proton–deuteron data from Ref. [67]. For remaining notation, see Fig. 4
Figure 8 shows an update of these finding by using the new SMS NN potentials of Ref. [7], including the 3NF at \(\hbox {N}^2\hbox {LO}\) and replacing the EKM approach to estimating truncation errors by the Bayesian model \(\bar{C}_{0.5-10}^{650}\). Specifically, the incomplete \(\hbox {N}^3\hbox {LO}\) and \(\hbox {N}^4\hbox {LO}\) results shown in this figure are based on the \(\hbox {N}^3\hbox {LO}\) and \(\hbox {N}^4\hbox {LO}^+\) NN potentials accompanied by the \(\hbox {N}^2\hbox {LO}\) 3NF with the LECs \(c_D\) and \(c_E\) being readjusted to the \(^3\)H binding energy and the differential cross section at \(E_{\mathrm{lab}}^N = 70~\hbox {MeV}\) in exactly the same way as done at \(\hbox {N}^2\hbox {LO}\). In the 3NF, we have used the values of the LECs \(c_i\) from Ref. [12] consistent with the NN interactions at the corresponding chiral order, namely \(c_1= -1.07\) GeV\(^{-1}\), \(c_3 = -5.32\) GeV\(^{-1}\) and \(c_4 = 3.56\) GeV\(^{-1}\) at \(\hbox {N}^3\hbox {LO}\) and \(c_1= -1.10\) GeV\(^{-1}\), \(c_3 = -5.54\) GeV\(^{-1}\) and \(c_4 = 4.17\) GeV\(^{-1}\) at \(\hbox {N}^4\hbox {LO}\), subject to the additional shifts of
generated by the pion loop contributions to the 3NF at \(\hbox {N}^3\hbox {LO}\) [17]. Since we do not have complete results beyond \(\hbox {N}^2\hbox {LO}\), the error bands in Fig. 8 are obtained by just rescaling the corresponding \(68\%\) and \(95\%\)\(\hbox {N}^2\hbox {LO}\) DoB intervals. The incomplete \(\hbox {N}^3\hbox {LO}\) and \(\hbox {N}^4\hbox {LO}^+\) results may, of course, still be regarded as complete \(\hbox {N}^2\hbox {LO}\) predictions. At \(E_{\mathrm{lab}}^N = 200~\hbox {MeV}\), the \(\hbox {N}^3\hbox {LO}\) uncertainty bands are still quite sizable indicating that the \(\hbox {N}^4\hbox {LO}\) contributions to the 3NF could play a significant role. Thus, fully in line with the findings of Ref. [39, 40], our results suggest that the accurate description of Nd scattering data below pion production threshold will likely require the chiral expansion of the 3NF to be pushed to \(\hbox {N}^4\hbox {LO}\). Notice that the accurate and precise description of neutron–proton and proton–proton data below pion production threshold also required the chiral expansion of the NN force to be pushed to \(\hbox {N}^4\hbox {LO}\) (or even \(\hbox {N}^4\hbox {LO}^+\)) [7].
4 Subleading short-range 3NF: an exploratory study
To include the 3NF contributions beyond \(\hbox {N}^2\hbox {LO}\) one needs to regularize the corresponding pion loop contributions consistently with the NN interactions of [7] in a chirally symmetric manner; see Refs. [35, 36] for discussion. Such consistently regularized pion-exchange contributions to the 3NF are not yet available beyond \(\hbox {N}^2\hbox {LO}\). In addition to long- and intermediate-range interactions generated by pion-exchange diagrams, the chiral 3NF involves ten purely short-range operators at \(\hbox {N}^4\hbox {LO}\), which have been worked out in Ref. [32]. In the exploratory study of Ref. [68], the effects of these subleading short-range terms are investigated in proton–deuteron scattering below \(E_{\mathrm{lab}}^p = 3~\hbox {MeV}\) within the hybrid approach based on phenomenological two- and three-nucleon forces. The authors of Ref. [68] have succeeded to fit the coefficients of the short-range operators to obtain a good description of experimental data at \(E_{\mathrm{lab}}^p = 3~\hbox {MeV}\). However, except for fine-tuned observables like the neutron–deuteron doublet scattering length, which is well known to be correlated with the \(^3\)H binding energy, and the analyzing powers \(A_y\) and \(i T_{11}\); see Fig. 4, the scattering observables at such low energies are dominated by the NN interaction, and the 3NFs are expected to play a minor role [64]. On the other hand, large discrepancies between theory and data are observed in Nd elastic scattering at intermediate and higher energies, where the 3NFs are expected to play a prominent role [69].
To explore the role of the subleading short-range 3NF contributions we choose two out of the ten terms, namely the isoscalar central and spin–orbit interactions
where \({\mathbf {K}}_i = ({\mathbf {p}}_i ' + {\mathbf {p}}_i )/2\), while \(E_1\) and \(E_7\) denote the corresponding LECs. We apply the same nonlocal Gaussian regulator as employed in the \(\hbox {N}^2\hbox {LO}\) short-range part of the 3NF and defined in Eq. (13), and restrict ourselves to the cutoff \(\varLambda = 450~\hbox {MeV}\). In this study we do not attempt to determine the LECs \(E_1\) and \(E_7\) from data but explore effects of these 3NF terms for fixed values of these LECs. Specifically, in a complete analogy with Eq. (11), we express \(E_i\) in terms of dimensionless coefficients \(c_{E_i}\) according to
For the coefficients \(c_{E_i}\), we consider in this study the fixed values of \(c_{E_i} = \pm 2\). Based on the variation of the LEC \(c_E\), which, depending on the cutoff \(\varLambda = 400 \ldots 550~\hbox {MeV}\), takes the values in the ranges of \(c_E = -0.3 \ldots 1.2\), \(c_E = -0.3 \ldots 0.8\), \(c_E = -0.6 \ldots 0.7\) and \(c_E = -0.5 \ldots 0.7\) when using the NN potential at orders \(\hbox {N}^2\hbox {LO}\), \(\hbox {N}^3\hbox {LO}\), \(\hbox {N}^4\hbox {LO}\) and \(\hbox {N}^4\hbox {LO}^+\), respectively, we expect that the actual values of \(c_{E_i}\) should lie well within the range spanned by \(c_{E_i} = \pm 2\). The adopted naturalness estimates are, however, subject to convention-dependent ambiguitiesFootnote 4 and should be taken with care. A more meaningful and reliable assessment of the natural size of the LECs can be carried out in the spectroscopic basis as done in Ref. [7] for the NN potentials. This, however, would require the inclusion of a complete set of independent contact operators in the 3NF at \(\hbox {N}^4\hbox {LO}\), which goes beyond the scope of our study.
For a given observable, the impact of the subleading short-range 3NF terms can be quantified by comparing the results for \(c_{E_i} = \pm 2\) with those for \(c_{E_i} = 0\) after renormalization. Since we are only able to perform implicit renormalization by expressing the bare LECs in terms of low-energy observables, this requires a readjustment of the LECs \(c_D\) and \(c_E\) for every considered set of \(c_{E_1}\), \(c_{E_7}\). To allow for a meaningful interpretation of the obtained results, we follow in each case exactly the same fitting procedure as explained in Sect. 3 by adjusting \(c_D\) and \(c_E\) to the \(^3\)H binding energy and the Nd cross section minimum at \(E_{\mathrm{lab}}^N = 70~\hbox {MeV}\). In Fig. 9 we show the resulting dependence of \(c_D\) and \(c_E\) on \(c_{E_1}\) and \(c_{E_7}\).
In addition to looking at the absolute values of the various dimensionless LECs, it is also instructive to compare the corresponding expectation values in the triton state, which are listed in Table 1. For the considered SMS regulator and \(\varLambda = 450~\hbox {MeV}\), the expectation value of the two-pion exchange 3NF contributions amounts to \(\langle V_{2 \pi } \rangle = -0.19~\hbox {MeV}\).Footnote 5 Notice that the spin–orbit 3NF term \(\propto c_{E_7}\) does not contribute to the S-wave partial waves in the triton state (for the employed angle-independent regulator), and is found to provide a negligible contribution to the \(^3\)H binding energy. The apparent contradiction with the findings of Ref. [68] regarding this term is presumably caused by a different regulator employed in that paper. The expected natural contribution of the 3NF can be estimated based on naive dimensional analysis via \(| \langle V_{3N} \rangle | \sim Q^3 | \langle V_{2N} \rangle | \sim 0.3^3 \times 40 \text{ MeV } \sim 1~\hbox {MeV}\). For \(c_{E_7} = \pm 2\), the individual terms in the 3NF already start exceeding their expected natural size, thus indicating that the considered values for this LEC likely overestimate its natural range. This conclusion is supported by Nd scattering results, as described below.
We are now in a position to discuss the effects of the subleading short-range contributions to the 3NF in selected Nd scattering observables. To that aim, we first perform calculations based on the NN SMS potential of Ref. [7] at \(\hbox {N}^4\hbox {LO}^+\) together with the 3NF at \(\hbox {N}^2\hbox {LO}\).
The resulting predictions for the Nd elastic scattering observables lie in the middle of the blue bands as shown in Figs. 10, 11 for the intermediate energy of \(E_{\mathrm{lab}}^N = 135~\hbox {MeV}\). The light- and dark-shaded blue bands show the \(95\%\) and \(68\%\) DoB intervals for the truncation error at \(\hbox {N}^3\hbox {LO}\). These error bands do not properly reflect the uncertainty of our calculation, which is only complete through \(\hbox {N}^2\hbox {LO}\), but show the expected size of \(\hbox {N}^4\hbox {LO}\) corrections estimated within the employed Bayesian approach. Next, we repeat the calculations by switching on the \(\hbox {N}^4\hbox {LO}\) short-range terms in the 3NF. The resulting predictions for \(c_{E_1} = \pm 2\) and \(c_{E_7} = 0\) (\(c_{E_1} = 0\) and \(c_{E_7} = \pm 2\)) are shown in Figs. 10, 11 by the purple dashed-dotted (blue dashed) lines. As expected from the estimated truncation uncertainty at \(\hbox {N}^3\hbox {LO}\), the considered \(\hbox {N}^4\hbox {LO}\) 3NF terms yield sizable contribution to the Nd scattering observables at this rather high energy, especially in the region of the cross section minimum and at backward angles. In fact, the magnitude of the \(c_{E_1} = \pm 2\) corrections compares well with the width of the \(\hbox {N}^3\hbox {LO}\) error bands, especially with the ones corresponding to \(95\%\) DoB intervals. This finding supports our expectation that the actual value of this LEC should be well within the considered range of \(c_{E_1} = \pm 2\). On the other hand, the contributions of the \(c_{E_7}\)-term lie in most cases outside of the \(\hbox {N}^3\hbox {LO}\) truncation bands, which suggests that the employed values of \(c_{E_7} = \pm 2\) overestimate the natural size of this LEC. This conclusion is in line with the pattern shown in Fig. 9, which indicates a significantly larger shifts in the LECs \(c_D\), \(c_E\) induced by changing \(\delta c_{E_7} = \pm 2\) as compared with the ones induced by \(\delta c_{E_1} = \pm 2\).
Another interesting observation is that both the \(c_{E_1}\) and \(c_{E_7}\) contributions tend to lie well within the truncation bands at forward angles while outside at backward angles. This could indicate a shortcoming of the employed Bayesian model, which relies on Eq. (2) and does not explicitly account for higher momentum scales being probed in backward scattering as compared with forward scattering; see also Ref. [51] for a similar conclusion.
Results for the differential cross section, nucleon and deuteron analyzing powers \(A_y^n\) and \(A_y^d\) as well as deuteron tensor analyzing powers \(A_{yy}\), \(A_{xz}\) and \(A_{xx}\) in elastic nucleon–deuteron scattering at laboratory energy of \(E_\mathrm{lab}^N = 135~\hbox {MeV}\) based on the SMS NN potentials of Ref. [7] at \(\hbox {N}^4\hbox {LO}^+\) in combination with the 3NF at \(\hbox {N}^2\hbox {LO}\) using \(\varLambda = 450~\hbox {MeV}\). Blue light- (dark-) shaded bands show the expected truncation uncertainty for a complete\(\hbox {N}^3\hbox {LO}\) calculation and are obtained by multiplying the \(\hbox {N}^2\hbox {LO}\) truncation error corresponding to \(95\%\) (\(68\%\)) DoB intervals for the model \({\bar{C}}_{0.5-10}^{650}\) with the expansion parameter \(Q \simeq 0.45\). Short-dashed-dotted and long-dashed-dotted red lines show the impact of the \(\hbox {N}^4\hbox {LO}\) central short-range 3NF \(\propto c_{E_1}\) with \(c_{E_1} = -2\) and \(c_{E_1} = 2\), respectively. Similarly, short-dashed and long-dashed blue lines show the impact of the \(\hbox {N}^4\hbox {LO}\) spin–orbit short-range 3NF \(\propto c_{E_7}\) with \(c_{E_7} = -2\) and \(c_{E_7} = 2\), respectively. For remaining notation, see Fig. 6
Same as Fig. 12 but for polarization transfer coefficients \(K_{xx}^y\), \(K_{y}^y\), \(K_{yy}^y\), \(K_{xz}^y\), \(K_{xx}^y-K_{yy}^y\) and the induced polarization \(P^y\)in elastic nucleon–deuteron scattering at laboratory energy of \(E_{\mathrm{lab}}^N = 135~\hbox {MeV}\). For remaining notation, see Figs. 6 and 10
At the lowest considered energy of \(E_{\mathrm{lab}}^N = 10~\hbox {MeV}\), the effects of the considered \(\hbox {N}^4\hbox {LO}\) 3NF terms turn out to be small; see the right panel of Fig. 12 for a representative example, except for the nucleon and deuteron vector analyzing powers \(A_y^N\) and \(A_y^d\). These results provide yet another confirmation of the fine-tuned nature of these observables, see the discussion in Sect. 3, and indicate that the apparent \(A_y\)-puzzle could be naturally resolved at the \(\hbox {N}^4\hbox {LO}\) level by the corresponding short-range contributions to the 3NF [70]; see also Ref. [71] for a related discussion within pionless EFT. While the spin–orbit 3NF is well known to have a strong impact on the vector analyzing power [72], we found, quite surprisingly, that the isoscalar central short-range 3NF term \(\propto c_{E_1}\) also significantly affects \(A_y\). A comparison of effects due to the \(c_{E_i}\)-terms with the estimated truncation error bands at \(E_{\mathrm{lab}}^N = 10~\hbox {MeV}\) leads to the same conclusions as reached at the higher energy of \(E_\mathrm{lab}^N = 135~\hbox {MeV}\).
Last but not least, we have also looked at the differential cross section in the so-called symmetric space star configuration of the Nd breakup reaction at \(E_{\mathrm{lab}}^N = 13~\hbox {MeV}\), which is known to represent another low-energy puzzle in the three-nucleon continuum. Contrary to \(A_y\), this observable is dominated by the S-wave components of the NN force and turns out to be highly insensitive to the 3NFs considered so far; see Ref. [73] for recent results based on the SCS nuclear potentials of Refs. [5, 6, 38]. We found that the inclusion of the \(c_{E_i}\)-terms has a negligible effect on the cross section in this breakup configuration, and the obtained results essentially coinside with the one presented in Ref. [73]. The observed discrepancy between the theoretical calculations and the neutron–deuteron and proton–deuteron data thus indeed appears to be puzzling. It will be interesting to see if this puzzle can be resolved by the inclusion of the \(\hbox {N}^3\hbox {LO}\) and remaining \(\hbox {N}^4\hbox {LO}\) contributions to the 3NF.
5 Summary and conclusions
In this paper we analyzed selected Nd scattering observables at \(\hbox {N}^2\hbox {LO}\) in chiral EFT based on the SMS interactions of Ref. [7]. The main results of our study can be summarized as follows.
Following the approach of Ref. [41], we have explored several pointwise Bayesian models for quantifying truncation uncertainties in chiral EFT and tuned them by calculating angular distributions of neutron–proton scattering.
Using the SMS NN forces of Ref. [7] accompanied with the \(\hbox {N}^2\hbox {LO}\) 3NF regularized in the same way, we have determined the LECs \(c_D\) and \(c_E\) entering the 3NF from the \(^3\)H binding energy and the Nd cross section data of Ref. [62] at \(E_{\mathrm{lab}}^N = 70~\hbox {MeV}\). The resulting \(\hbox {N}^2\hbox {LO}\) results for elastic Nd scattering observables agree within errors with our earlier \(\hbox {N}^2\hbox {LO}\) calculations based on the SCS interactions [38] and with experimental data.
The truncation errors for various Nd scattering observables estimated in [38,39,40] using the approach of Ref. [5] are found to be consistent with \(68\%\) DoB intervals for the employed Bayesian model \({\bar{C}}_{0.5-10}^{650}\). In particular, we confirm our earlier findings in Refs. [39, 40], obtained using the SCS NN forces of Refs. [5, 6], that Nd scattering at intermediate and higher energies provides a “golden window” to study higher-order contributions to the chiral 3NF.
Based on the estimated truncation errors, we argue that a high-precision description of neutron–deuteron scattering below pion production threshold will likely require pushing the chiral expansion of the 3NF at least to \(\hbox {N}^4\hbox {LO}\). This conclusion is in line with the convergence pattern of the chiral expansion in the NN sector [7] and is further substantiated by an exploratory study of the short-range 3NF contributions \(\propto E_{1,7}\), which, as expected from our Bayesian analysis, are found to have significant effect on Nd polarization observables at intermediate energies. Our results show that the contributions to the 3NF beyond \(\hbox {N}^2\hbox {LO}\) are indeed potentially capable of resolving the discrepancies between theory and experiment observed in Nd scattering at intermediate and higher energies [69].
As a next step, this study should be extended to \(\hbox {N}^3\hbox {LO}\), which will require the inclusion of consistently regularized 3NF contributions. Work along these lines is in progress by the LENPIC Collaboration.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: All relevant data are already given in the figures and in the table.]
Notes
The term “semilocal” refers to a local regularization approach for long-range interactions in combination with a nonlocal cutoff for the short-range part of the nuclear force.
The \(\hbox {N}^4\hbox {LO}^+\) potentials of Ref. [7] include four sixth-order (i.e. \(\hbox {N}^5\hbox {LO}\)) short-range operators contributing to F-waves, which are needed to describe certain very precisely measured proton–proton scattering observables at intermediate and high energies. The same operators are included in the \(\hbox {N}^4\hbox {LO}\) potentials of Ref. [8].
For the EKM approach, the success rate equals \(100\%\) per construction.
For example, our notation for the spin–orbit term \(\propto E_7\) in the 3NF differs from the one adopted in Ref. [68] by a factor of 1/2.
We emphasize again that the regularized two-pion exchange contributions contain admixtures of the short-range \(c_D\)- and \(c_E\)-like terms. When adopting the convention of Ref. [7] by explicitly subtracting the short-range pieces of the two-pion exchange, the expectation value changes to \(\langle V_{2 \pi } \rangle = -0.61~\hbox {MeV}\), showing that the \(\hbox {N}^2\hbox {LO}\) 3NF is actually dominated by the long-range pieces.
References
E. Epelbaum, H.W. Hammer, U.-G. Meißner, Rev. Mod. Phys. 81, 1773 (2009)
E. Epelbaum, U.-G. Meißner, Ann. Rev. Nucl. Part. Sci. 62, 159 (2012)
R. Machleidt, D.R. Entem, Phys. Rep. 503, 1 (2011)
D.R. Entem, N. Kaiser, R. Machleidt, Y. Nosyk, Phys. Rev. C 91(1), 014002 (2015)
E. Epelbaum, H. Krebs, U.-G. Meißner, Eur. Phys. J. A 51(5), 53 (2015)
E. Epelbaum, H. Krebs, U.-G. Meißner, Phys. Rev. Lett. 115(12), 122301 (2015)
P. Reinert, H. Krebs, E. Epelbaum, Eur. Phys. J. A 54(5), 86 (2018)
D.R. Entem, R. Machleidt, Y. Nosyk, Phys. Rev. C 96(2), 024004 (2017)
A. Ekström et al., Phys. Rev. Lett. 110(19), 192502 (2013)
A. Gezerlis, I. Tews, E. Epelbaum, S. Gandolfi, K. Hebeler, A. Nogga, A. Schwenk, Phys. Rev. Lett. 111(3), 032501 (2013)
M. Piarulli, L. Girlanda, R. Schiavilla, R. Navarro Pérez, J.E. Amaro, E. Ruiz Arriola, Phys. Rev. C 91(2), 024003 (2015)
M. Hoferichter, J. Ruiz de Elvira, B. Kubis, U.-G. Meißner, Phys. Rev. Lett. 115(19), 192301 (2015)
R. Navarro Pérez, J. E. Amaro, E. Ruiz Arriola, Phys. Rev. C 88(6), 064002 (2013) Erratum: [Phys. Rev. C 91(2), 029901 (2015)]
R.B. Wiringa, V.G.J. Stoks, R. Schiavilla, Phys. Rev. C 51, 38 (1995)
R. Machleidt, Phys. Rev. C 63, 024001 (2001)
V.G.J. Stoks, R.A.M. Klomp, C.P.F. Terheggen, J.J. de Swart, Phys. Rev. C 49, 2950 (1994)
V. Bernard, E. Epelbaum, H. Krebs, U.-G. Meißner, Phys. Rev. C 77, 064004 (2008)
V. Bernard, E. Epelbaum, H. Krebs, U.-G. Meißner, Phys. Rev. C 84, 054001 (2011)
E. Epelbaum, Phys. Lett. B 639, 456 (2006)
E. Epelbaum, Eur. Phys. J. A 34, 197 (2007)
S. Kölling, E. Epelbaum, H. Krebs, U.-G. Meißner, Phys. Rev. C 80, 045502 (2009)
S. Kölling, E. Epelbaum, H. Krebs, U.-G. Meißner, Phys. Rev. C 84, 054008 (2011)
H. Krebs, E. Epelbaum, U.-G. Meißner, Ann. Phys. 378, 317 (2017)
H. Krebs, E. Epelbaum, U.-G. Meißner, Few Body Syst. 60(2), 31 (2019)
S. Pastore, R. Schiavilla, J.L. Goity, Phys. Rev. C 78, 064002 (2008)
S. Pastore, L. Girlanda, R. Schiavilla, M. Viviani, R.B. Wiringa, Phys. Rev. C 80, 034004 (2009)
M. Piarulli, L. Girlanda, L.E. Marcucci, S. Pastore, R. Schiavilla, M. Viviani, Phys. Rev. C 87(1), 014006 (2013)
A. Baroni, L. Girlanda, S. Pastore, R. Schiavilla, M. Viviani, Phys. Rev. C 93(1), 015501 (2016) Erratum: [Phys. Rev. C 93, no. 4, 049902 (2016)] Erratum: [Phys. Rev. C 95,(5), 059901 (2017)]
H. Krebs, E. Epelbaum, U.-G. Meißner, arXiv:2001.03904 [nucl-th]
H. Krebs, A. Gasparyan, E. Epelbaum, Phys. Rev. C 85, 054006 (2012)
H. Krebs, A. Gasparyan, E. Epelbaum, Phys. Rev. C 87(5), 054007 (2013)
L. Girlanda, A. Kievsky, M. Viviani, Phys. Rev. C 84, 014001 (2011)
J. Golak et al., Eur. Phys. J. A 43, 241 (2010)
K. Hebeler, H. Krebs, E. Epelbaum, J. Golak, R. Skibinski, Phys. Rev. C 91(4), 044001 (2015)
H. Krebs, PoS CD2018 (2019) 098
E. Epelbaum, Towards high-precision nuclear forces from chiral effective field theory, in Proceedings of the International Conference ‘Nuclear Theory in the Supercomputing Era—2018’ (NTSE-2018), Daejeon, South Korea, October 29–November 2, ed. by A.M. Shirokov, A.I. Mazur (Pacific National University, Khabarovsk, Russia, 2019), p. 150
A.A. Slavnov, Nucl. Phys. B 31, 301 (1971)
E. Epelbaum et al. [LENPIC Collaboration], Phys. Rev. C 99(2), 024313 (2019)
S. Binder et al. [LENPIC Collaboration], Phys. Rev. C 93(4), 044002 (2016)
S. Binder et al. [LENPIC Collaboration], Phys. Rev. C 98(1), 014002 (2018)
J.A. Melendez, S. Wesolowski, R.J. Furnstahl, Phys. Rev. C 96(2), 024003 (2017)
D. Siemens, V. Bernard, E. Epelbaum, A. Gasparyan, H. Krebs, U.-G. Meißner, Phys. Rev. C 94(1), 014620 (2016)
D.L. Yao, D. Siemens, V. Bernard, E. Epelbaum, A.M. Gasparyan, J. Gegelia, H. Krebs, U.-G. Meißner, JHEP 1605, 038 (2016)
D. Siemens, V. Bernard, E. Epelbaum, A.M. Gasparyan, H. Krebs, U.-G. Meißner, Phys. Rev. C 96(5), 055205 (2017)
A.N. Hiller Blin, Z.F. Sun, M.J. Vicente Vacas, Phys. Rev. D 98(5), 054025 (2018)
J. Hu, Y. Zhang, E. Epelbaum, U.-G. Meißner, J. Meng, Phys. Rev. C 96(3), 034307 (2017)
R. Skibiński et al. Phys. Rev. C 93(6), 064002 (2016)
N. Nevo Dinur, O.J. Hernandez, S. Bacca, N. Barnea, C. Ji, S. Pastore, M. Piarulli, R.B. Wiringa, Phys. Rev. C 99(3), 034004 (2019)
R.J. Furnstahl, N. Klco, D.R. Phillips, S. Wesolowski, Phys. Rev. C 92(2), 024005 (2015)
E. Epelbaum, PoS CD2018 (2019) 006
J.A. Melendez, R.J. Furnstahl, D.R. Phillips, M.T. Pratola, S. Wesolowski, Phys. Rev. C 100(4), 044001 (2019)
M. R. Schindler and D. R. Phillips, Ann. Phys. 324, 682 (2009) Erratum: [Annals Phys. 324, 2051 (2009)]
V.G.J. Stoks, R.A.M. Klomp, M.C.M. Rentmeester, J.J. de Swart, Phys. Rev. C 48, 792 (1993)
A.J. Bersbach et al., Phys. Rev. D 13, 535 (1976)
U. van Kolck, Phys. Rev. C 49, 2932 (1994)
E. Epelbaum, A. Nogga, W. Glöckle, H. Kamada, U.-G. Meißner, H. Witała, Phys. Rev. C 66, 064001 (2002)
S.A. Coon, H.K. Han, Few Body Syst. 30, 131 (2001)
C.R. Howell et al., Few Body Syst. 2, 19 (1987)
K. Sagara et al., Phys. Rev. C 50, 576 (1994)
G. Rauprich et al., Few Body Syst. 5, 67 (1988)
F. Sperisen et al., Nucl. Phys. A 422, 81 (1984)
K. Sekiguchi et al., Phys. Rev. C 65, 034003 (2002)
K. Sekiguchi et al., Phys. Rev. C 70, 014001 (2004)
W. Glöckle, H. Witała, D. Huber, H. Kamada, J. Golak, Phys. Rep. 274, 107 (1996)
E. Epelbaum, H. Krebs, U.-G. Meißner, Nucl. Phys. A 806, 65 (2008)
U.-G. Meißner, A.I.P. Conf. Proc. 1011(1), 49 (2008)
B. von Przewoski et al., Phys. Rev. C 74, 064003 (2006)
L. Girlanda, A. Kievsky, M. Viviani, L.E. Marcucci, Phys. Rev. C 99(5), 054003 (2019)
N. Kalantar-Nayestanaki, E. Epelbaum, J.G. Messchendorp, A. Nogga, Rep. Prog. Phys. 75, 016301 (2012)
L. Girlanda, A. Kievsky, M. Viviani, L.E. Marcucci, EPJ Web Conf. 113, 04009 (2016)
A. Margaryan, R.P. Springer, J. Vanasse, Phys. Rev. C 93(5), 054001 (2016)
A. Kievsky, Phys. Rev. C 60, 034001 (1999)
H. Witała et al., Few Body Syst. 60(1), 19 (2019)
Acknowledgements
Open Access funding provided by Projekt DEAL. This study has been performed within Low Energy Nuclear Physics International Collaboration (LENPIC) project and was supported by BMBF (contract No. 05P18PCFP1), by DFG and NSFC (TRR 110), by DFG (project No. 279384907-CRC 1245), by the Polish National Science Centre under Grants No. 2016/22/M/ST2/00173 and 2016/ 21/D/ST2/01120, by the Chinese Academy of Sciences (CAS) Presidents International Fellowship Initiative (PIFI) (Grant No. 2018DM0034) and by VolkswagenStiftung (Grant No. 93562). Numerical calculations were performed on the supercomputer cluster of the JSC, Jülich.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by V. Somà
Appendices
A Analytic expressions for the posterior pdf
For the Gaussian prior pdf in Eq. (6), the posterior pdf takes the form
where \({\bar{q}}^2 \equiv \sum _{i = k+1}^{k+h} Q^{2i}\), \({\varvec{c}}_k^2 \equiv \sum _{i\in A} c_i^2\) and the incomplete gamma function is defined as
For the non-informative prior \(C_\epsilon \) with \({\bar{c}}_< = \epsilon \) and \({\bar{c}}_> = 1/\epsilon \), the expression for the posterior, after taking the limit \(\epsilon \rightarrow 0\) for \({\varvec{c}}_k^2 \ne 0\) simplifies to
B Partial-wave decomposition of the N4LO 3NF contact terms
For the \(E_1\)-term
we choose the Faddeev component \(V_{3N}^{(1)} = 2\, E_1\, {\mathbf {q}}_1^2\) invariant with respect to the interchange of nucleons 2 and 3 and obtain
with \({\mathbf {p}}\), \({\mathbf {p}} '\), \({\mathbf {k}}\) and \({\mathbf {k}} '\) denoting the corresponding Jacobi momenta.
For the \(E_7\)-term
we choose the Faddeev component
and obtain
with \({{\hat{X}}} \equiv 2 X + 1\). For more details on notation, see Ref. [64].
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Epelbaum, E., Golak, J., Hebeler, K. et al. Towards high-order calculations of three-nucleon scattering in chiral effective field theory. Eur. Phys. J. A 56, 92 (2020). https://doi.org/10.1140/epja/s10050-020-00102-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epja/s10050-020-00102-2