How to renormalize integral equations with singular potentials in effective field theory

We briefly review general concepts of renormalization in quantum field theory and discuss their application to solutions of integral equations with singular potentials in the few-nucleon sector of the low-energy effective field theory of QCD. We also describe a particular subtractive renormalization scheme and consider a specific application to a toy-model with a singular potential serving as its effective field theoretical leading-order approximation.


I. INTRODUCTION
The issue of renormalization is a subject of long-standing controvercy in the few-nucleon sector of the low-energy effective field theory (EFT) of the strong interactions, pioneered in Ref. [1]. Recent reviews and references addressing this problem can be found, e.g., in Refs. [2][3][4][5][6][7][8][9][10]. The Lagrangian of any consistent EFT includes all local interactions allowed by the underlying symmetries. Therefore, all ultraviolet (UV) divergences appearing in physical quantities are in principle canceled by the corresponding counterterms [11]. However, the practical implementation of such a quantum field theoretical renormalization is a highly non-trivial problem in the few-body sector of the low-energy EFT of the strong interactions. The difficulty is caused by the fact that due to the large absolute values of the scattering lengths, an infinite number of renormalized diagrams has to be summed up at any order. To perform this summation one defines an effective potential as a sum of irreducible diagrams contributing to the process under consideration and then obtains the scattering amplitude by solving the corresponding integral equation [1]. In practice, solving integral equations corresponds to summing up regularized but not renormalized diagrams. Regularized non-perturbative expressions obtained by solving such integral equations, when expanded in powers of (corresponding to the loop expansion), reproduce the regularized perturbative diagrams. While non-perturbative solutions to integral equations may indeed contain pieces which yield vanishing contributions to the perturbative series, the integral equations in the low-energy EFT are used solely as tools for performing partial resummations of diagrams. Therefore, properly renormalized non-perturbative expressions in a self-consistent EFT must reproduce the renormalized perturbative series, when expanded in .
It is well known how to renormalize Feynman diagrams in quantum field theories. However, it turns out to be a challanging problem to properly renormalize the non-perturbative solutions to integral equations since this usually requires taking into account contributions of an infinite number of counterterms. The problem is actually caused by the non-renormalizability of EFT in the traditional sense, manifested by the presence of singular potentials in the few-body sector.
The low-energy EFT of the strong interactions yields a systematic expansion of the effective potentials in terms of small masses and momenta, where small is taken relative to the typical hard (breakdown) scale of about 500 MeV. However, these expressions are to be trusted only at low energies. Translated into coordinate space, the EFT yields an expansion of the potential applicable at large distances, while the (strongly scheme-dependent) short-range part is determined by the desired resolution. Accordingly, extrapolations of the long-range part of the effective potential to short distances do not root in the underlying QCD. Moreover, QCD features only shallow bound states in multi-nucleon systems thus indicating that the singular van der Waals-like behavior of one-and multi-pion exchange potentials is not a valid approximation of the nuclear interaction at short distances.
Reading the literature in the few-body sector of low-energy chiral EFT and references therein, one often gets the impression that there is no well-established concept of renormalization for non-perturbative problems in quantum field theories, and that we are currently seeking for such a self-consistent definition. In this quest for the solution of a technically challenging problem, fundamental concepts that are well established and should be well-known to the practitioners are often ignored or even disputed in the nuclear physics community. Our understanding is that the very definition of the EFT and renormalization as it is understood in quantum field theory (QFT) specifies uniquely what needs to be done. The implementation of a proper quantum field theoretical renormalization is, however, a rather challenging problem and subject to further investigations.
In our recent work on the renormalization in the few-body sector of the EFT of the strong interactions [12] we have put the main emphasis on discussing common misconceptions, i.e. what should not be done when renormalizing the solutions to integral equations with singular EFT potentials. In this paper we again address this issue, but focus on what should be done to achieve the proper renormalization in EFT with singular potentials.
Our paper is organized as follows. In section II we discuss the definition and conceptual issues of renormalization in QFTs. In section III we consider a particular QFT renormalization scheme which can be applied to the solutions of integral equations with singular potentials. An application of the introduced renormalization scheme to a singular potential for a specific toy model is given in section IV. We summarize our work in Sec. V.

II. RENORMALIZATION
Renormalization is a very natural and familiar procedure which is actually carried out in most theories. Considering a theory depending on the parameters g i , with i = 1, ..., N , one fixes these parameters from N experiments and makes predictions for the results of other experiments. The standard example is QED, where the electron charge and mass need to be determined from experiment and then e.g. the anomalous magentic moment of the electron can be predicted to high precision. In mathematical language this means that we calculate N "reference" physical quantities σ i (E) = f i ( , g 1 , g 2 , . . . , g N , E), for i = 1, ..., N , as functions of the couplings g k , of kinematical variables collectively denoted by E, and , taken to be zero for any classical theory. By taking the expressions of these physical quantities for some fixed kinematics µ j (j = 1, . . . , N ) we express the g i as functions of observables σ i (µ i ) Next, we substitute the g i from Eq. (1) in the calculated expressions of all physical quantities and obtain where the F i are some functions, whose specific form depends on the considered theory. In this way we express all physical quantities in terms of the "reference" physical quantities taken at the normalization points µ i , instead of the original "bare" parameters. Performed in the framework of QFT this procedure is known as renormalization. Notice that renormalization is non-perturbative if calculations are done non-perturbatively and perturbative if we are using perturbation theory. As a rule in QFTs, the coefficients of the expansion of physical quantities in powers of the bare couplings g i are divergent. In some theories divergences disappear after renormalization, that is, the expansion of all physical quantities in terms of the "reference" physical quantities in the limit of a removed regulator turn out to have finite coefficients to all orders. Such theories are called renormalizable (as e.g. the already mentioned case of QED). To render the Green functions finite one often introduces renormalized fields and corresponding field renormalization constants. However, physical quantities do not depend on the choice of fields and, therefore, one can also work with bare fields. To simplify the discussion we formulate our considerations in such a framework.
In practice one uses more conveniently chosen functions g R i (µ) = G i (σ 1 (µ 1 ), σ 2 (µ 2 ), . . . , σ N (µ N ), µ 1 , . . . , µ N ) (j = 1, . . . , N ) as new renormalized parameters, where µ stands collectively for all µ 1 , . . . , µ N . Expressing the original parameters g i in terms of the g R i (µ) and substituting in the Lagrangian of the given QFT, we can express the Lagrangian in terms of the renormalized parameters. Usually, the functions G i are chosen such that leading to where the δg k i are counterterms of k-th order in the loop expansion which themselves can be expanded in powers of renormalized coupling constants. Using the Lagrangian parameterized in terms of renormalized couplings and counterterms, we perform calculations of physical quantities in terms of the renormalized parameters g R i (µ). In renormalizable theories, a perturbative calculation of physical quantities to any finite order in the renormalized coupling constants leads to Taylor series with finite coefficients. All counterterm contributions up to given order in the loop expansion and in the coupling constants have to be taken into account in these calculations. If, for whatever reason, the expansion in any of these renormalized couplings needs to be summed up to an infinite order, then this has to be done systematically by also including the contributions stemming from all counterterms. While such a calculation might be technically very complicated or even unfeasible, the formalism of QFT renormalization defines very precisely and uniquely what has to be done.
From the modern point of view, renormalizable (in the traditional sense) theories are viewed as leading order approximations to effective field theories. The corresponding Lagrangians contain an infinite number of terms compatible with underlying symmetries. A self-consistent EFT is renormalizable in the sense that all divergences appearing in physical quantities can be absorbed in the redefinition of an infinite number of parameters of the effective Lagrangian (i.e. N is equal to infinity in an EFT). Despite the dependence on an infinite number of parameters, EFTs do not loose predictive power if it is possible to find renormalization schemes that ensure that low-energy observables calculated with a specified accuracy depend only on a finite number of renormalized parameters. Notice that physical quantities calculated to some finite order depend on the choice of the renormalization scheme. While this dependence is formally of higher order in the EFT expansion, the higher order corrections to physical quantities are suppressed only for appropriately chosen renormalization conditions. While it is always possible to choose a scheme that destroys the validity of perturbation theory, it is not always possible to find a scheme for which perturbative series make sense. A well-known example is QCD, where, on the one hand, perturbation theory becomes very useful at high energies, due to the asymptotic vanishing of the renormalized coupling, provided the renormalization scale is chosen of the order of the characteristic energy. On the other hand, it is impossible to find a renormalization scheme at low energies which would lead to meaningful perturbative expressions for physical quantities in powers of the renormalized coupling constant.
To be more specific, the Lagrangian of an EFT has the generic form where K is the kinetic part, O i are interaction terms with bare couplings g i and ψ stands collectively for the fields. 1 Because of the UV divergences one needs to introduce some kind of regularization. The final expressions of physical quantities in terms of the renormalized couplings are finite in the limit of removed regularization and the results do not depend on the choice of the specific regularization scheme. As mentioned above, to render the Green functions finite one often introduces renormalized fields and corresponding field renormalization constants, however, physical quantities do not depend on the choice of fields and therefore one can also work with bare ones in the regularized 1 It is understood that the considered EFT leads to a systematic perturbative expansion of physical quantities for energies E Q, where Q is some large scale.
theory. To calculate physical quantities one expresses the bare couplings in terms of the renormalized ones as specified above and rewrites the Lagrangian as where we also indicate the dependence of the counterterms on the regularization scheme by including the parameter Λ in the list of arguments. While one obviously has a huge freedom of choosing the renormalization conditions, performing consistent calculations of physical quantities up to any order (either finite or infinite) in any coupling g R j for any particular choice requires the inclusion of all contributions up to a given order, generated by the Lagrangian of Eq. (5). It is difficult to implement such formally consistent program in the few-body sector of the low-energy EFT of the strong interactions, where one needs to sum up an infinite number of diagrams. The problem is caused by the necessity of taking into account contributions of an infinite number of counterterms with growing complexity of the operator structures O i (ψ) in the (heavy baryon) formulation of the EFT. The best available approach to the non-relativistic formulation of an EFT with pions and nucleons as explicit degrees of freedom is the cutoff EFT which relies on the inclusion of a finite number of counterterms and keeping the cutoff parameter such that the difference between the properly renormalized and actually obtained results is of a higher order than the accuracy of the current calculation [4,[13][14][15][16][17][18].
A convenient way of implementing the quantum field theoretical renormalization is to apply the BPHZ subtractive procedure using Zimmerman's forest formula, see, e.g., Ref. [19]. This subtractive approach is equivalent to the above outlined renormalization by including the contributions of all counterterms appearing at a given order. Within the subtractive approach, instead of explicitly including the counterterm contributions, one renormalizes the loop diagrams by subtracting sub-divergences and overall divergences using some fixed renormalization conditions and treats all coupling constants as renormalized finite parameters. In the next section we consider a particular subtractive quantum field theoretical renormalization for the non-perturbative solutions to integral equations with singular effective potentials which removes all divergences from all diagrams generated by iterations of the integral equations.

III. INTEGRAL EQUATION WITH A SINGULAR POTENTIAL AND RENORMALIZATION
In this section we consider one particular realization of a subtractive EFT renormalization for solutions of integral equations with singular potentials. Our scheme is a generalization of the one employed in Refs. [1,20] to the case of the OPE potential and differs from those of Refs. [21][22][23][24][25]. In particular, Refs. [24,25] apply a subtraction scheme at vanishing energy, whose realization in terms of local counterterms of the effective Lagrangian is unclear to us. On the other hand Refs. [21][22][23] demand that the scattering amplitude calculated at a given order must be independent of the subtraction point. As emphasized in the previous section, our understanding is that perturbative amplitudes calculated up to any finite order are renormalization-scheme dependent. A toy-model example considered in the next section demonstrates how problematic it may turn out to find a renormalization scheme, which is not only conceptually consistent but also useful in describing the data.
We emphasize that there is a conceptual difference between the EFT renormalization of the scattering amplitude obtained by solving the Lippmann-Schwinger (LS) equation with singular potentials, that takes into account contributions of an infinite number of counterterms with operator structures of increasing complexity in each partial wave, and self-adjoint extensions of singular Hamiltonians in the framework of quantum mechanics, that ensure finiteness of the calculated observables in the removed regularization limit [26].
We start by considering the LS integral equation for an off-shell scattering amplitude which we write symbolically as To avoid any non-analytic dependence on the quark masses of the counterterms in the few-body sector of the chiral EFT, we split the potential into a singular V S and a regular V R part, where V S has only a polynomial dependence on the quark masses and V R is such that its overlap in loop diagrams with V S , i.e. V S GV R , is finite. For the potential of Eq. (7) the solution to Eq. (6) can be written in the form where T S and T C satisfy the equations and To carry out subtractive renormalization, we replace T S in Eq. (8) by the renormalized amplitude T r S that satisfies the equation where G s stands for the subtracted Green function. In particular, we expand G ≡ G(E) at a fixed value E = −E µ and subtract first several terms, the number of the subtracted terms depending on the UV behaviour of the singular potential, i.e. we have Solving Eq. (11) corresponds to taking into account contributions of an infinite number of counterterms which take care of all sub-divergences and overall divergences. Let us demonstrate this for the first few iterations of this equation In the first iteration, , and then the overall divergence is subtracted by adding 2 V S G e V S G e V S . Renormalization works analogously for further iterations. The effective potential V r which includes all necessary counterterms subtracting sub-divergences and overall divergences in all iterations of Eq. (11) can be obtained by solving the following equation That is, the subtracted amplitude T r S can also be obtained by solving the equation where V r is obtained by solving Eq. (14).
The final renormalized expression of the amplitude has the form where T r C satisfies the equation Below we show that the final renormalized amplitude T r can be obtained in a single step by solving the integral equation To this end we write T r as T r S + T 1 , substitute in Eq. (11) and obtain the following equation for T 1 Next, we define the amplitude T 2 via T 1 =: T 2 (1 + GT r S ) and obtain from Eq. (19) ( Further, we introduce T x through T 2 =: (1 + T r S G)T x and substitute it in Eq. (20) to obtain Writing the equation (11) in the form and using it in Eq. (21), we obtain Using the definition ofṼ R , Eq. (23) reduces to Comparing Eqs. (24) and (17), we see that T x is identical with T r C and hence indeed T r satisfies Eq. (16). We emphasize that we consider the above subtractive scheme not because it is amongst the best ones for, e.g., nucleon-nucleon scattering in chiral EFT, but rather since it can be easily applied in practice and, therefore, well suited for demonstration purposes. As noticed in Refs. [27,28] in the NN scattering problem, it is preferable to use subtraction schemes with two scales corresponding to the leading-and higher-order potentials.

IV. APPLICATION TO A TOY MODEL
As a specific toy-model example of how to apply the discussed renormalization scheme, we construct an effective potential for the underlying "fundamental" potential considered in Ref. [12] where M is the light mass and the heavy masses m 1 , m 2 represent the large scales. Our choice of parameters is α = −36 GeV −2 , M = 0.1385 GeV, m 1 = 0.75 GeV and m 2 = 1.15 GeV. The strength of the interaction α is taken equal for all terms, so that the potential V (r) vanishes for r → 0 and it behaves as −α e −M r /r 3 for large r. More details on this model can be found in Ref. [12]. We consider the LS equation for the S-wave scattering amplitude in the center-of-mass frame of two particles with unit masses with E = k 2 /1 GeV being the total energy, and p = |p|, p = |p | the relative momenta of the incoming and outgoing particles, respectively. We set to its physical value, i.e. = 1, in these calculations. It is easy to construct the corresponding EFT for the underlying "fundamental" potential in Eq. (25). The LO EFT interaction consists of a delta potential and the long-range part −α e −M r /r 3 , which is singular if extended to the small-r region. Following Ref. [12], we choose the coupling constant α = −36 GeV −2 ≈ −1/(0.167 GeV) 2 such that the full LO potential is non-perturbative for the momenta k ∼ M = 0.1385 GeV. A simple UV analysis shows that the LO potential is perturbatively non-renormalizable, i.e. to remove the divergences from its iterations one needs to introduce counterterms of higher orders (in momenta and energy) which themselves generate new divergences, etc., up to infinity.
We apply the renormalization scheme described in the previous section to the LO effective potential, for which we take where µ 1 is the scale parameterizing the division of the potential into the singular and regular parts. To obtain the singular part we expanded the LO potential in Taylor series of the mass M around M = µ 1 and kept the first two terms in the power series of M − µ 1 . Notice that since our toy model is not generated by a chirally invariant Lagrangian, we obtain a power series expansion in M , not in M 2 . We fit the renormalized coupling c R to the lowenergy phase shifts corresponding to the underlying toy-model potential. Subtractions of loop integrals are performed at E = −µ 2 /(1 GeV). Our results for the phase shift for various choices of the renormalization scales are shown in Fig. 1 together with the phase shift of the underlying toy-model. As it is seen from the figure, one can find such a choice of parameters µ and µ 1 for which the description of the data is good for the LO approximation, however, the results for the momenta M are sensitive to the changes of these values.

V. SUMMARY
In this paper we discussed the problem of renormalization of singular potentials in the low-energy effective field theory of the strong interactions. We outlined the standard procedure of quantum field theoretical renormalization emphasizing its general definition that applies to both non-perturbative as well as perturbative calculations. The main message of our work is that chiral EFT in the few-nucleon sector does not require the invention of any new non-perturbative renormalization that would go beyond what is already well-established in quantum field theory.
To demonstrate that the problem lies not in conceptual issues but rather in the complicated technical implementation of the standard renormalization approach, we proposed a way of implementing the subtractive renormalization for solutions of integral equations with singular potentials. The suggested renormalization scheme corresponds to subtracting all divergences appearing at a given order, and also (over-)subtracting finite pieces with some fixed normalization conditions. While being just one of infinitely many possibilities of fixing the normalization conditions, and very likely not the best one, the suggested scheme has the advantage of being easily applicable in practice and serves well our purposes of demonstrating a self-consistent EFT renormalization.
As an example to demonstrate these rather general statements, we applied the considered subtractive renormalization scheme to the toy model of Ref. [12]. While there exists a region of values of the renormalization scales for which the description of the "data" is good for the LO approximation, the results are sensitive to the changes of these values as one expects on general grounds.