# Convergent perturbative nuclear effective field theory

###### Abstract

We consider the nuclear effective field theory including pions in the two-nucleon sector in the S waves up to including the next-to-next-to-leading order (NNLO) terms according to the power counting suggested by the Wilsonian renormalization group analysis done in a previous paper. We treat only the leading contact interaction nonperturbatively, and the rest, including the long-distance part of pion exchange, are treated as perturbations. To define the long-distance part, it is important to introduce a separation scale, or a cutoff. We employ a hybrid regularization, in which the loops with only contact interactions are regularized with Power Divergence Subtraction (PDS), while the loops with (long-distance part of) pion exchange are regularized with a Gaussian damping factor (GDF), to simplify the (nonperturbative) leading-order amplitudes. The scale introduced by PDS is identified with the cutoff of GDF up to a numerical factor. We emphasize that the introduction of the GDF requires a careful definition of the coupling constant for the pion exchange. We obtain the analytic expressions for the phase shifts for the and - channels. By fitting them to the Nijmegen partial wave analysis data, it is shown that the effective theory expansion with perturbative long-distance part of pion exchange is converging.

^{†}

^{†}preprint: KYUSHU-HET-138

## I Introduction

Since the seminal papers by Weinberg Weinberg (1990, 1991), nuclear effective field theory (NEFT), an effective field theory describing systems with more than one nucleons at low energies, has attracted much attention. See Ref. Epelbaum et al. (2009) for a recent review. With a vast number of papers written over twenty years, one might think that it has become a matured discipline. In fact, there are NLO calculations Epelbaum et al. (2005); Entem and Machleidt (2003) in the literature.

However, a very fundamental issue is still under discussions: the power counting and renormalization. Since an effective field theory contains an infinitely many operators, a power counting rule is necessary to organize calculations to achieve a certain order of accuracy. Cutoff dependence of physical quantities must be absorbed in the coupling constants order by order, that is, consistently with the power counting. The original Weinberg’s power counting, which is widely used in most of numerical calculations, is known to be inconsistent, i.e., cutoff dependence arising at a certain order can be absorbed only by terms of higher orders Kaplan (1997). Several authors (e.g. Nogga et al. (2005)) consider variants of the Weinberg’s scheme and discuss the nonperturbative renormalization of them.

An alternative, consistent power counting scheme (known as “KSW power counting”) was proposed by Kaplan, Savage and Wise Kaplan et al. (1998a, b), and independently by van Kolck van Kolck (1999). In their scheme, pion exchanges are treated as perturbation. Fleming, Mehen, and Stewart Fleming et al. (2000) however showed that the effective theory expansion fails to converge at NNLO, due to the strong tensor force of pion exchange at short distances. Beane, Bedaque, Savage, and van Kolck Beane et al. (2002) proposed a remedy, in which the part of pion exchange is treated nonperturbatively.

In a series of papers Harada et al. (2006); Harada and Kubo (2006); Harada et al. (2009, 2011), we advocate that the power counting should be based on the scaling dimensions obtained by Wilsonian renormalization group (RG) analysis. (See Refs. Birse et al. (1999); Birse (2006); Nakamura (2005) for related works.) The idea is that power counting is an order of magnitude estimate based on the dimensional analysis, and that the quantum notion of the dimension of an operator is the scaling dimension. A nonperturbative RG analysis is necessary because in the S waves the physical two-nucleon system is governed by a nontrivial fixed point which is inaccessible in perturbation theory.

In a previous paper Harada et al. (2011), we perform the Wilsonian RG analysis for the nucleon-nucleon scattering in the S waves in the NEFT including pions. We emphasize that it is important to divide pion exchange into its short-distance part (S-OPE) and the long-distance part (L-OPE) separated by a cutoff scale, because they play different roles. The S-OPE is represented as contact terms. A part of the S-OPE is relevant in the RG sense and thus should be treated nonperturbatively while the L-OPE is treated as perturbation. It turns out that the power counting for the contact interactions and the L-OPE is very similar to the KSW power counting.

In this paper, we consider the NEFT including pions in the two-nucleon sector in the S waves with the power counting suggested by the Wilsonian RG analysis mentioned above. In order to separate the pion exchange into two parts, we introduce an explicit separation scale, or a cutoff, , and propose a new hybrid regularization in which the diagrams with only contact interactions are regularized with Power Divergence Subtraction (PDS)Kaplan et al. (1998a, b), but those with potential pion exchange (L-OPE) are regularized with a Gaussian damping factor (GDF). We obtain the analytic expressions for the phase shifts for and at NNLO and fit them to Nijmegen partial wave analysis (PWA) data to determine the low-energy constants (LECs). All the technical details together with the lengthy analytic expressions will be given elsewhere.

Our approach is similar to that by Beane, Kaplan, and Vuorinen Beane et al. (2009) (BKV) in the respect that a separation scale is introduced. There are however important differences: (i) We use the same regularization both for the and - channels, though BKV introduce the separation scale only for the channel. (ii) We use a GDF to regularize the pion potential, while BKV use a Pauli-Villars type regulator, which we find insufficient to render several diagrams convergent. (iii) We interpret the separation scale as an analog of the floating cutoff in the Wilsonian RG analysis so that it should not exceed the physical cutoff 350 MeV above which the effective field theory description does not hold, while BKV consider a rather large value in the range 600 MeV 1000 MeV, although they consider it as a low-momentum scale of . (iv) We interpret the “renormalization scale” appeared in the PDS as the separation scale too, and take . BKV consider that is independent of and choose . (v) In our formulation, the separation scale is smaller than or equal to the physical cutoff , but otherwise arbitrary. On the other hand, BKV tune the value of to optimize the perturbation expansion.

## Ii Hybrid regularization

Since we consider that the introduction of the separation scale is essential, we could work only with a GDF even for the contact interactions, but such a scheme is very complicated because the leading order (relevant) operator becomes a linear combination of several operators when higher order operators are included. On the other hand, in the PDS regularization the leading operator stays the same even when the higher order operators are included. The hybrid regularization takes advantage of both regularizations.

Our interpretation of the “renormalization scale” introduced in the PDS as a separation scale equivalent to might look strange. It however comes from the comparison of power-divergent loop integrations calculated in the PDS and those with an explicit cutoff. The power divergences of the integral are represented as polynomials of in the PDS. They are typically related as , which is inferred by calculating the simplest one-loop diagram with two regularizations. Note that the interpretation is different from the usual one for the dimensional regularization such as scheme in the relativistic field theory, where only the logarithmic divergences are explicitly treated. There the scale is arbitrary, but to avoid the large logarithms, it is taken to be the typical scale of the process in question. In the present case, however, is not a low-momentum scale, and a typical momentum scale satisfies .

We introduce a GDF for the pion exchange so that they represent the L-OPE. For the scalar part, we can explicitly extract the S-OPE from the decomposition

(1) |

where the first term may be considered as a contact interaction and the loops containing this term are regularized with the PDS. The second term is the L-OPE, for which we introduce a GDF . For the tensor part, on the other hand, because of the tensor structure, it is impossible to extract local operators in a similar procedure. We assume that the effects of the S-OPE is already encoded in the coupling constants of the contact interactions, and consider only the L-OPE.

A diagram with pion exchange may be expanded in powers of , where is an external momentum or the pion mass. Thus a single diagram produces a series of contributions of different orders. For example, a diagram which appears at NLO may also contain NNLO contributions. This is a new feature of the hybrid regularization.

A very nontrivial point with the GDF regularization comes from the requirement that the definition of the pion exchange coupling constant should be independent of the cutoff. In order to satisfy the requirement, one needs to introduce an extra factor . Including the coupling constant, the scalar part of the pion exchange may be written as

(2) |

where stands for a contact operator that should be treated with PDS. One can see that this definition of the coupling constant has several favorable features: (i) The coupling constant defined by the residue of the pole of the Yukawa potential becomes independent of the separation scale . (ii) The results of loop integrals including pion propagators with GDFs contain the factor , which is canceled by the extra factor. If the extra factor is not included, there would be (disastrous) nonlocal contributions to higher orders when expanded in powers of .

The calculated amplitudes depend on (and ). All the dependence is actually eliminated by making the LECs -dependent, i.e., there is no nonlocal -dependence.

## Iii NNLO Fitting to the Nijmegen data

We calculate the LO, NLO, and NNLO amplitudes for the nucleon-nucleon scattering in the and - channels analytically and fit the phase shifts obtained from the amplitudes to the Nijmegen PWA data. The Lagrangian is the same as that of Fleming et al. Fleming et al. (2000). We include diagrams containing an - mixing contact term together with the diagrams given in Ref. Fleming et al. (2000), and calculate with the hybrid regularization explained above. The so-called “radiation pion” and “soft pion” contributions, and the isospin breaking terms are not included. As input parameters, we set the nucleon mass MeV, the pion mass MeV, and the pion exchange coupling constant MeV. We choose MeV, which is about the physical cutoff. We have however checked that the quality of the fitting does not depend on the value of , though the fitted values of the LECs are different for different values of . A detailed account based on RG analysis will be given in a separated paper.

In order to fit the calculated phase shifts to the Nijmegen PWA data, we use the Mathematica command NonlinearModelFit with a weight function , emphasizing the low-momentum region. The range of the data varies with the order: the ranges of center-of-mass momentum – MeV, – MeV, and – MeV are used for LO, NLO, and NNLO respectively for both the and phase shifts. If one dares to try to fit to an unreasonably wider range, the quality of the fit in the low-momentum region becomes worse. We should not do so because NEFT must describe the lower-momentum region better than the higher-momentum region. Given the order of the expansion, the range in which the fitting is successful may be viewed as the validity region.

We show the calculated scattering phase shifts by using the fitted values of the LECs in Fig. 1 for the channel, and in Fig. 2 for the channel. Most importantly, no breakdown of the effective theory expansion is observed. By going to higher orders, the effective field theory fits better and with wider validity regions. The NNLO results are excellently fitted up to MeV.

There are one, three and six (seven) coupling constants to be fitted in the LO, NLO and NNLO calculations respectively in the (-) channel. We renormalize only the -dependence before fitting the phase shifts to the data, in order to avoid the cancellations among large numbers arising from the terms with positive powers of .

Tables 1 and 2 show how fitted values of the (-renormalized) dimensionless coupling constants change as we go higher orders. The dimensionless coupling constants for with , which are used in Ref. Fleming et al. (2000), are defined as , where stands for the channel. It is important to note that most of these coupling constants are of order one. It implies that our estimate of the magnitude of each term is correct. There are however large ambiguities in determining the values of coupling constants at NNLO. There seem to be “flat directions.”; sets of considerably different values of the LECs give almost the same phase shifts.

LO | -0.958 | — | — | — | — | — |
---|---|---|---|---|---|---|

NLO | -0.964 | 0.517 | 0.479 | — | — | — |

NNLO | -0.967 | 0.205 | 0.099 | -0.055 | 1.180 | 1.203 |

LO | -1.22 | — | — | — | — | — | — |
---|---|---|---|---|---|---|---|

NLO | -1.31 | 1.47 | 0.68 | — | — | — | — |

NNLO | -1.42 | 0.89 | -0.11 | -11.72 | 4.28 | 31.1 | -4.80 |

## Iv Summary and Discussions

We perform NNLO calculations for the scattering amplitudes for the nucleon-nucleon scattering in the S waves with the power counting suggested by the Wilsonian RG analysis done in a previous paper, which is very similar to the KSW power counting. A novel hybrid regularization is employed to introduce the scale , which separates the pion exchange into the S-OPE and the L-OPE. We fit the calculated phase shifts to the Nijmegen PWA data. The fitted values of most of the coupling constants are of the natural size and the effective field theory expansion seems converging up to including the NNLO.

In our approach, the “renormalization scale” appearing in the PDS is identified with the scale up to a numerical constant. The scale plays an analogous role to that of the floating cutoff in the Wilsonian RG analysis. The low-momentum physical quantities should not depend on the values of . This requirement leads to a set of renormalization group equations (RGEs) for the coupling constants. The requirement is satisfied order by order. With the RGEs and their solutions, we have a complete control over the dependence. The RGEs and their solutions will be given elsewhere.

We find that the introduction of the GDF requires a careful definition of the pion exchange coupling constant and an additional factor is necessary. This finding is possible only with analytic expressions and the RG analysis. We think that it is also important for numerical studies with similar damping factors.

Our approach will be applicable to other partial waves. The application to the P waves would be of particular interest, because BKV suspect that NLO calculations would be necessary to have a convergent result. (Gegelia also found serious problems in the P waves in his unpublished work.) Recently we have done a Wilsonian RG analysis for the P waves in a similar manner to that in Ref. Harada et al. (2011), and have argued that the pion exchange in the P waves demotes to higher order so that the counterterms are present to absorb the cutoff dependence arising from the loops containing pion exchangesHarada et al. (2013). The calculations of the scattering phase shifts with hybrid regularization are now in progress.

###### Acknowledgements.

This work is supported by JSPS KAKENHI (Grand-in-Aid for Scientific Research (C)) (22540286).## References

- Weinberg (1990) S. Weinberg, Phys. Lett. B251, 288 (1990).
- Weinberg (1991) S. Weinberg, Nucl. Phys. B363, 3 (1991).
- Epelbaum et al. (2009) E. Epelbaum, H.-W. Hammer, and U.-G. Meissner, Rev. Mod. Phys. 81, 1773 (2009), arXiv:0811.1338 [nucl-th] .
- Epelbaum et al. (2005) E. Epelbaum, W. Glockle, and U.-G. Meissner, Nucl.Phys. A747, 362 (2005), arXiv:nucl-th/0405048 [nucl-th] .
- Entem and Machleidt (2003) D. Entem and R. Machleidt, Phys.Rev. C68, 041001 (2003), arXiv:nucl-th/0304018 [nucl-th] .
- Kaplan (1997) D. B. Kaplan, Nucl. Phys. B494, 471 (1997), arXiv:nucl-th/9610052 .
- Nogga et al. (2005) A. Nogga, R. G. E. Timmermans, and U. van Kolck, Phys. Rev. C72, 054006 (2005), arXiv:nucl-th/0506005 .
- Kaplan et al. (1998a) D. B. Kaplan, M. J. Savage, and M. B. Wise, Phys. Lett. B424, 390 (1998a), arXiv:nucl-th/9801034 .
- Kaplan et al. (1998b) D. B. Kaplan, M. J. Savage, and M. B. Wise, Nucl. Phys. B534, 329 (1998b), arXiv:nucl-th/9802075 .
- van Kolck (1999) U. van Kolck, Nucl. Phys. A645, 273 (1999), arXiv:nucl-th/9808007 .
- Fleming et al. (2000) S. Fleming, T. Mehen, and I. W. Stewart, Nucl. Phys. A677, 313 (2000), arXiv:nucl-th/9911001 .
- Beane et al. (2002) S. R. Beane, P. F. Bedaque, M. J. Savage, and U. van Kolck, Nucl. Phys. A700, 377 (2002), arXiv:nucl-th/0104030 .
- Harada et al. (2006) K. Harada, K. Inoue, and H. Kubo, Phys. Lett. B636, 305 (2006), arXiv:nucl-th/0511020 .
- Harada and Kubo (2006) K. Harada and H. Kubo, Nucl. Phys. B758, 304 (2006), arXiv:nucl-th/0605004 .
- Harada et al. (2009) K. Harada, H. Kubo, and A. Ninomiya, Int. J. Mod. Phys. A24, 3191 (2009), arXiv:nucl-th/0702074 .
- Harada et al. (2011) K. Harada, H. Kubo, and Y. Yamamoto, Phys.Rev. C83, 034002 (2011), arXiv:1012.2716 [nucl-th] .
- Birse et al. (1999) M. C. Birse, J. A. McGovern, and K. G. Richardson, Phys. Lett. B464, 169 (1999), arXiv:hep-ph/9807302 .
- Birse (2006) M. C. Birse, Phys. Rev. C74, 014003 (2006), arXiv:nucl-th/0507077 .
- Nakamura (2005) S. X. Nakamura, Prog. Theor. Phys. 114, 77 (2005), arXiv:nucl-th/0411108 .
- Beane et al. (2009) S. R. Beane, D. B. Kaplan, and A. Vuorinen, Phys. Rev. C80, 011001 (2009), arXiv:0812.3938 [nucl-th] .
- Harada et al. (2013) K. Harada, H. Kubo, T. Sakaeda, and Y. Yamamoto, Few Body Syst. 54, 1483 (2013).