Coarse graining $\pi\pi$ scattering

We carry out an analysis of $\pi\pi$ scattering in the $IJ=00$, $11$ and $20$ channels in configuration space up to a maximal center-of-mass energy $\sqrt{s}=1.4$ GeV. We separate the interaction into two regions marked by an elementarity radius of the system; namely, a long distance region above which pions can be assumed to interact as elementary particles and a short distance region where many physical effects cannot be disentangled. The long distance interaction is described by chiral dynamics, where a two-pion-exchange potential is identified, computed and compared to lattice calculations. The short distance piece corresponds to a coarse grained description exemplified by a superposition of delta-shell potentials sampling the interaction with the minimal wavelength. We show how the so constructed non-perturbative scattering amplitude complies with the proper analytic structure, allowing for an explicit N/D type decomposition in terms of the corresponding Jost functions and fulfilling dispersion relations without subtractions. We also address renormalization issues in coordinate space and investigate the role of crossing when fitting the scattering amplitudes above and below threshold to Roy-equation results. At higher energies, we show how inelasticities can be described by one single complex and energy dependent parameter. A successful description of the data can be achieved with a minimal number of fitting parameters, suggesting that coarse graining is a viable approach to analyze hadronic processes.


I. INTRODUCTION
Hadronic interactions at low and intermediate energies are typically characterized by a combination of elementary and composite particle features. While at long distances hadrons behave as elementary particles and their interactions can be described in terms of purely color singlet degrees of freedom, at short distances their composite character becomes manifest in terms of quark and gluon fields in the fundamental and adjoint representations of the color group, respectively. The relevant scale separating between this dual description marks the onset of a confinement scale and we expect it to be of the order of the hadron size, which generally is found to be about 1 fm. While the hadronic dynamics can be organized quite often as a long distance perturbative hierarchy with an increasing number of exchanged particles, it is by itself incomplete; some further either ab initio or phenomenological information reflecting the underlying quark-gluon structure is needed to provide a full description of the scattering process.
The way how this separation is visualized in the complex energy plane is not completeley straightforward. Traditionally, and within a genuinely hadronic picture, one appeals to Mandelstam analyticity [1], i.e. the assumption that a scattering amplitude can be expressed by double dispersive integrals in terms of double-spectral density functions, where the integration ranges extend over those regions in the Mandelstam plane where the corresponding double-spectral functions have non-vanishing support [2]. This viewpoint is ulti-mately grounded in the Mandelstam conjecture, which holds in lowest order in the coupling constant in quantum field theory [1,3] or to all orders within a non-relativistic context in potential scattering [4], and, which, in the ππ scattering case, has been rigorously proved in a finite domain [5,6]. It is noteworthy that under this same assumption an equivalent local and energy dependent optical potential of non-relativistic form was derived many years ago by Cornwall and Ruderman [7,8]. For a balanced review on these issues at the textbook level see, for e.g., [9,10]. The existence of a finite analyticity domain suggests in turn the very existence of a finite cut-off on a purely hadronic basis but without an explicit reference to the underlying quark-gluon dynamics and in particular to the confinement scale, so that the cut-off may be determined phenomenologically from data.
Pion-pion scattering is the simplest reaction in QCD mediated by strong interactions involving the lightest hadrons. Tight theoretical constraints based on analyticity, crossing, unitarity, chiral symmetry and Regge behavior can be imposed (see for e.g. [11] for an early review). The machinery of effective field theories (EFT) [12] and in particular its implementation in Chiral Perturbation Theory (χPT) [13] has enabled as a consequence, the most precise theoretical extraction of the ππ S-wave scattering lengths to date with about an order of magnitude more precision than the experiment [14][15][16][17][18][19][20][21], an unprecedented case in strong interactions, where invariably just the opposite situation happens. A historic overview is given in [22]. Along these lines, the most precise ππ-scattering analyses to date have been obtained in [16,21]. The latter corresponds to a ππ description up to √ s = 1.42 GeV, obtained by fitting the available experimental data from πN → ππN and K e4 decays while imposing as further constraints Roy and arXiv:1807.10837v1 [hep-ph] 27 Jul 2018 Roy-like equations, and with statistical uncertainties satisfying the necessary normality requirements of the residual distributions [23], (see for e.g. [24,25] for reviews). We stress that despite all these tight mathematical constraints, most of its non-perturbative setup rests upon the validity of the Mandelstam conjecture [1,3], a result which, as already mentioned, has not yet been rigorously proven since it was first proposed in 1958. This tacit assumption will also be made throughout our work.
In the present paper we invoke the equivalent local and energy dependent optical potential approach suggested long ago in [7,8] to describe ππ scattering in coordinate space. In order to do so, we consider a relativistic Schrödinger equation and define a potential to describe the ππ interaction by matching the field theoretical result to an equivalent quantum mechanical problem in perturbation theory. Phenomenological precursors of ππ scattering analyses in coordinate space were prompted in [26,27] within the boundary condition model of strong interactions [28]. Equivalent coordinate space potentials using the Mandelstam representation or the Bethe-Salpeter equation as a starting point were also proposed to all orders in [29][30][31]. As it will become clear below, it is remarkable if not surprising that so little work on ππ scattering has been conducted within this approach as compared to more popular momentum space methods. Our work fills this gap by implementing Wilsonian ideas inspired by recent developments in the NN case [32][33][34]. These NN investigations had as a consequence a selection of the largest np and pp database up to energies about pion production threshold of 3σ mutually consistent data. Our present investigation within ππ is in a sense of exploratory character and it pretends also to provide some training playground with an eye put on the more compelling πN case, where the selection of the currently existing database is largely needed (see for e.g. [35][36][37] and references therein).
At short distances, where the interaction is non perturbative, we will assume a complete ignorance of the strong interaction behavior and consider a coarse graining of the interaction instead, very much in the spirit of the work done in [32][33][34] for the NN case. The basic idea is to separate the ππ interaction into an inner and outer region at a given separation distance, r c , located at about some elementarity radius. This radius is defined so that at larger distances pions behave effectively as point-like particles. We will assume that in this long distance regime their interactions are ruled by chiral symmetry, and hence they become calculable within χPT. Thus, for r > r c , we will construct a chiral potential with the correct low-energy analytic properties by matching both quantum mechanical and field theoretical scattering amplitudes in perturbation theory 1 . On the contrary, the inner region, r < r c , is regarded as unknown and sampled with the minimal de Broglie wave length determined by the maximum energy we want to describe. This corresponds to a coarse graining of the short range piece and, in its simplest realization, the inner potential will be written as a superposition of equidistant delta-shell interactions. A key issue is to confidently determine the numerical value of the separation scale r c , since, as noted in [40,41] and we will see below, the combination pr c will fix the total number of independent fitting parameters. The longest range interaction corresponds to a 2π-exchange which is O(e −2m π r ), so that a naive estimate suggests r c ∼ 2/(2 m π ) ∼ 1.4 fm 2 , a number which will be corroborated by our numerical analysis.
While the potential approach has been explained in great detail in previous works within the NN context (see for instance [41]), it is unconventional within the ππ scattering folklore. Thus, we will assume no previous knowledge from the side of the reader and for the sake of completeness we will briefly go through all the important issues along the paper. Moreover, ππ scattering is characterized because at resonance energies relativistic effects cannot be ignored. For instance, for the prominent case of the ρ-meson √ s = m ρ 2m π . Unlike the NN case, a new important aspect in the discussion is related to crossing symmetry, which actually intertwines the s, t and u channels 3 . In addition, the current extraordinary precision achieved theoretically in extracting the Swaves scattering lengths or the lightest ππ resonance pole parameters [42][43][44][45] provides a great confidence on the theoretical ideas supporting these benchmarking extractions. The fact that the coarse graining approach works for NN scattering in a regime where relativistic and inelastic effects become important, such as pp scattering up to √ s ∼ 2 GeV [41], suggests extending the method to other hadronic reactions under similar operating conditions 4 .
Finally, for the sake of completeness let us mention that lattice calculations are naturally formulated in coordinate space. These calculations attack the problem on the finite lattice spacing and the finite volume in two different fashions: either an (energy dependent) potential is determined and the Schrödinger equation is solved subsequently in the continuum or, alternatively, the energy level shifts are determined on the lattice and converted into phase-shifts by means of Luscher's formula [46]. Actually, in a pioneering work [47] , the I = 2 ππ S-wave scattering phase shifts from Lattice QCD have been determined. Later on, ππ scattering has been studied from N f = 2 + 1 and N f = 2 + 1 + 1 flavors in [48] and [49], respectively. Connected and Disconnected Contractions have also been analyzed in [50]. In addition, the ππ I = 2 channel has also been studied within the potential approach [51]. A comparison between potential and Luscher's approaches has been undertaken in [52] for the I = 2 case, with rather similar results. We remind that both methods have potential draw-backs. On the one hand, the potential method uses interpolating fields which may distort the physics at short distances, and we will explicitly show that in a chiral expansion such potential presents a short distance singularity, which evades the conventional solutions of the Schrödinger equation. On the other hand, the current applicability of this Luscher's method [46] requires the interaction to sharply vanish at the edges of the volume (in the relative coordinate), a fact that has been often ignored in momentum space treatments (see for e.g. [53,54]) but needs to be established for the ππ case. Our analysis below supports this assumption.
The paper is organized as follows. In Section II, we provide a general and brief field theoretical overview of ππ scattering to fix our notation in a way that our problem can be easily formulated. In Section III, we show our choice for a quantum mechanical description in terms of a complex local and energy dependent optical potential. We analyze the long-range contributions within χPT in Section IV, where an expression for the potential is obtained from the discontinuities of the ππ scattering amplitude in the t-channel. This requires introducing a short distance cut-off to handle the strong short distance power divergences of the chiral potential, an issue which we discuss at length in Section V. In Section VI, we analyze the concept of effective elementarity in order to display in two examples how the elementarity radius depends on the particular process. In Section VII, we address the problem of coarse graining ππ interactions with and without the longrange contributions. The analytical properties of the scattering amplitude and the relation of our approach with the N/D method is discussed in Section VIII. The implementation of inelasticities within a coarse grained perspective is explained in Section IX. We also analyze some aspects concerning low energy constants and the number of parameters in Section X. Finally, in Section XI we summarize our main results and provide some outlook for future work.

II. FORMALISM FOR ππ SCATTERING
We start by summarizing the relevant formulae for ππ scattering to fix our notation and to provide a proper perspective of our subsequent analysis. A comprehensive presentation at the textbook level can be seen in [11] and also in the lecture [55]. More recent upgrades can be consulted in [24,25] A. Kinematics For a pion state ϕ α with α = {±, 0}, the π α (p 1 ) + π β (p 2 ) → π γ (p 1 ) + π δ (p 2 ) relativistically invariant scattering amplitude can be written as with s = (p 1 + p 2 ) 2 , t = (p 1 − p 1 ) 2 and u = (p 1 − p 2 ) 2 the standard choice of Mandelstam variables. If we take ϕ ± = (φ 1 ±iφ 2 )/ √ 2 and ϕ 0 = φ 3 , with φ a ·φ b = δ ab , in the Cartesian basis we obtain T ab;cd = A(s,t, u)δ ab δ cd + B(s,t, u)δ ac δ bd +C(s,t, u)δ ad δ bc , where A(s,t, u) stands for the π + π − → π 0 π 0 amplitude. This amplitude is the the only independent one thanks to isospin, crossing and Bose-Einstein symmetries, B(s,t, u) = A(t, s, u) and C(s,t, u) = A(u,t, s). Denoting T I (s,t, u) as the isospin combination with well defined isospin I in the s-channel, one has For the normalization, we will use here the conventions in [56][57][58]. The partial-wave decomposition in the s-channel becomes where z = 1 + 2t/(4m 2 π − s) is the s-channel scattering angle, m π = 139.57 MeV the pion mass, P J (cos θ ) the Legendre polynomials and t IJ (s) is the partial-wave projection of the ππ scattering amplitude with isospin I and total angular momentum J. Thus, for waves fulfilling the relation (−1) J+I = 1 one has with the ππ phase factor and δ IJ the scattering phase shift. The inelasticity η IJ (s) = 1 for s < 16m 2 π and the unitarity condition for the partial wave amplitude reads in the elastic region Imt IJ (s) = σ (s)|t IJ (s)| 2 for 4m 2 π ≤ s ≤ 16m 2 π . (6) Of course, for s > 16m 2 π one has absorption η IJ (s) < 1 and inelastic processes such as 2π → nπ take place at √ s = 0.56, 0.84, 1.12 and 1.40 GeV for n = 4, 6, 8 and 10, respectively, as well as KK and ηη at √ s ∼ 1 GeV, etc. In our discussion we will also use the quantum mechanical amplitude f IJ (p) defined by with p the CM momentum. For elastic scattering one has f IJ (s) −1 = p cot δ IJ − ip, so that at low energies one has the threshold expansion with a IJ and b IJ the lowest threshold parameters. An equivalent way of representing the low energy behavior is tan δ IJ (s) p 2J+1 = a IJ + b IJ p 2 + . . . , (9) or by an effective range expansion where α IJ = −a IJ and r IJ /2 = −b IJ /a 2 IJ is the effective range, which is generally positive (see below). Usually, the expansion (9) works for small scattering lengths, such as ππ whereas (10) works for large scattering lengths, such as NN (see, e.g. , [59,60] for a discussion)

B. Anatomy of the ππ interaction
The purpose of the present paper is to coarse grain the unknown pieces of the ππ interaction in configuration space. It is thus important to gather some features emerging from comprehensive studies over the last decades [14][15][16][17][18][19][20][21]. According to these findings the partial wave expansion in (3) is decomposed into two contributions: the low energy contribution described by means of a partial pave (PW) expansion to finite order and the high energy contribution assumed to be given by the leading Regge trajectories, which accounts for the long and short distance behavior of the scattering amplitude respectively. A standard quantum mechanical argument based on the impact parameter provides in the semi-classical limit and for an interaction of finite range r c , the number of necessary partial waves 5 . The impact parameter is defined as b = L/p with p the CM momentum and L the orbital angular momentum, which in our case equals the total angular momentum J. The quantization condition for the angular moment yields L ≈ J(J + 1) ∼ (J + 1/2) for J 1. For a finite range, the maximal impact parameter where scattering happens is b max ∼ r c . Thus, for a maximum CM momentum p max , the maximum angular momentum J max for which the phase shift is compatible with zero within uncertainties is For a maximum energy s max = 2 GeV 2 , corresponding to p max ∼ 0.7 GeV, it was found in [21] that waves beyond J max = 4 are vanishingly small for ππ scattering. Therefore, one obtains from (12) a range r c ∼ 1.3 fm. This simple estimate will be explicitly exploited below as an educated guess. Low energies close to threshold are encoded by the threshold parameters, see (8) and (10). The S-wave scattering 5 These arguments provide in addition a justification for analyticity [61,62]. lengths are α 00 = −0.3 fm and α 20 = 0.03 fm whereas for the P-wave we have α 11 = −(0.48 fm) 3 [21]. These are unnaturally small numbers compared with our above estimate of the range of the interaction, r c ∼ 1.3 fm and the elementarity radius, r e ∼ 1.2 fm (see the discussion in Section VI). While the behavior of the isotensor S-wave resembles a repulsive core, with positive effective range r 20 = 131.4 fm, the effective range in the isoscalar S-wave and isovector P-wave are negative, r 00 = −8.08 fm and r 11 = −5.25 fm, respectively. For S waves the Wigner causality bound [63] (see also [64]) restricts the maximum value of the effective range by the inequality which for I = 2 implies r c ≥ 0.95 fm. The positivity of the effective range is not implied by this condition, and is usually violated in the presence of resonances. This requires some unconventional shape for the S-wave potential as we will see.

C. Chiral Perturbation theory
The scattering amplitude can be computed perturbatively in Quantum Field Theory and in particular in χPT as a sum of Feynmann diagrams in an expansion in 1/ f , with f ∼ 86 MeV the pion weak decay constant in the chiral limit. In the partial waves basis the expansion can schematically be written as where t (n) . To one loop order, they were first computed in [12,13] and the relevant non-polynomial contributions are reproduced for completeness in appendix A. Explicit analytical expressions for the corresponding partial wave amplitudes are displayed in [39]. They obey the perturbative unitarity relation Imt (4) IJ (s) = σ (s)|t (2) IJ (s)| 2 , 4m 2 π ≤ s ≤ 16m 2 π . (15) At lowest order (LO) in the chiral expansion the threshold parameters are unnaturally small, a fact naturally accommodated by χPT with pions coupled derivatively.

D. Unitarization vs Crossing
The requirement of crossing is a fundamental one which stems from the local character of Quantum Field Theories. Chiral Perturbation Theory implements this symmetry at any order in the chiral expansion. The problems with perturbation theory, however, are on the one hand the lack of exact unitarity given by (6) and on the other hand the impossibility of describing outstanding non-perturbative features such as the generation of resonances, which emerge as poles of the scattering amplitude on unphysical Riemann sheets. Within a χPT framework, many methods have been proposed (see for instance [39,65,66] and references therein) based on imposing exact unitarity while matching perturbation theory at low energies. Most of them are nothing but algebraic tricks or functional solutions to a set of a priori conditions. As such, unitarization methods are not unique but strongly driven by experimental information, which explains partly their success. The Bethe-Salpeter method discussed at length in [39] preserves an identification of Feynman diagrams but it is not free from field reparameterizations or off-shell ambiguities. In addition, they violate crossing symmetry, which, in general, is only fulfilled order by order, although these violations can be statistically not-significant [67].
In Sections IV and VII, we will propose yet a new method based on first defining an equivalent quantum mechanical problem and, more importantly, on coarse graining the interaction. Of course, above the inelastic threshold s ≥ 16m 2 π one may wonder what condition should be imposed instead of just (6) 6 . We will extend the coarse graining idea to the case with inelasticities.

A. Relativistic equation
At the maximum CM energy we will be considering in this work s max = 2 GeV 2 , relativity and inelasticities are crucial physical ingredients since firstly √ s max m π and secondly we can produce up to n = ( √ s max − 2m π )/m π ∼ 8 pions as well as one KK and ηη pair in the final state. From a field theoretical point of view, this could be solved by using a multichannel Bethe-Salpeter equation for the several 2π, 4π, 6π, 8π, KK and ηη coupled channels, but it would be an extremely difficult task, which has never been accomplished to our knowledge. Even in the simplest elastic case the offshell ambiguities are present for calculations with a truncated kernel [38,39]. In order to grasp the nature of the ambiguities, consider for instance the case of π 0 (p 1 )π 0 (p 2 ) → π 0 (k 1 )π 0 (k 2 ) scattering, in the elastic regime. The Bethe-Salpeter (BS) equation reads, where p = p 1 −p 2 2 , k = k 1 −k 2 2 , P = p 1 + p 2 and q ± = P/2 ± q. ∆(q ± ) = 1/(q 2 ± − m 2 π + i0 + ) is the free pion propagator and V P (p, k) and T P (p, k) stand for the two-particle irreducible kernel or potential and the scattering amplitude, respectively. The factor 1/2 comes from the scattering of identical particles.
While the BS equation has been the subject of extensive research for a given potential, the main point of [38,39] was the flexible interpretation of the BS equation within χPT or more generally within EFT. Indeed, while the potential V P (p, k) 6 Usually the coupled channel unitarity condition is implemented instead.
Typically analyses within such a setup leave out the "small" multiple production channels, 2π → nπ, see e.g. [68] and works cited therein.
Thus, there is an inherent ambiguity in the definition and form of the potential, which has no consequences perturbatively but become relevant in the solution of the BS equation (16) where the off-shellness enters explicitly. This was mended in [38,39] by invoking an on-shell scheme, namely considering only on shell intermediate states, i.e. q 2 = s/4 − m 2 π and P · q = 0, so that the on-shell amplitude T (s,t) depends only on the on-shell potential V (s,t). Unfortunately, it also gives rise to pathologies in the coupled channel case producing spurious singularities due to an improper treatment of the crossed-channel exchanges [68]. The present paper pretends to address crossed-channel exchanges without invoking the on-shell scheme.

B. Invariant mass and equivalent Schrödinger equation
We will follow here the invariant mass formulation [69] 7 , already used for NN scattering with an optical potential [41,70]. This is the simplest way of retaining relativity without solving a BS equation but with a phenomenological optical potential that we review here for completeness. The idea is to write the total squared mass operator as where W represents the (invariant) interaction, which can be determined in the CM frame by matching in the nonrelativistic limit to a non-relativistic potential V ( x). This yields for ππ scattering after quantizationM 2 = 4(p 2 + m 2 π ) + 4m π V , withp = −i∇. Thus, the relativistic equation can be written asM 2 Ψ = 4(p 2 + m 2 π )Ψ, with p the CM momentum, i.e. as a non-relativistic Schrödinger equation This corresponds to the simple rule that one may effectively implement relativity by just promoting the non-relativistic CM momentum to the relativistic CM momentum. This minimal relativity ansatz is as good as the more fundamental one based on the Bethe-Salpeter equation as long as we use scattering data to determine the corresponding potential rather than an ab initio determination (see Ref. [39] for an in-depth discussion). 7 These authors wondered if there was a way to promote non-relativistic fits of NN scattering to a relativistic formulation without refitting parameters. The answer is in the affirmative by just reinterpreting the CM momentum by its relativistic counterpart.
To take into account the inelasticity within the masssquared construction, we assume a local and energydependent phenomenological potential, V ( r, s) = ReV ( r, s) + i ImV ( r, s), which could be obtained by fitting inelastic scattering data. Due to causality, the optical potential in the schannel satisfies a dispersion relation for each CM radial distance r of the form [7] ReV (r, s) = V (r) + 1 where √ s 0 = 4m π is the first 4π inelastic threshold and V (r) is an energy independent component. The complete potential includes also the crossed u channel component. The simple looking equation (19), together with the fixed-r dispersion relation (20), incorporates the necessary physical ingredients present in any theoretical approach: relativity and inelasticity consistent with analyticity.

C. Isospin and exchange potential
The incorporation of isospin into the game is straightforward. Rotational, isospin and particle exchange invariance requires the representation of the potential to be given by where V D and V X stand for the direct and exchange potential pieces, respectively. P 12 is the particle exchange operator, which implements the Bose-Einstein symmetry and that can be factorized as P 12 = P x P I . Moreover, for states with a well defined total isospin I = I 1 + I 2 , we can use the relation I 1 · I 2 = I(I + 1)/2 − 2 with I = 0, 1, 2, so that P I = (−1) I . In addition, for angular momentum eigenstates P x = (−1) J , so that P 12 = (−1) I+J . Therefore, in the isospin basis the potential can be decomposed as where we have introduced the projection operators fulfilling the orthogonality relations P I P I = δ II P I . In the partial wave representation, the exchange symmetry of the potential is preserved by just solving the Schrödinger equation for the direct potential for the allowed IJ channels with (−1) J+I = 1. In addition, for a spherically symmetric potential we have the usual factorization of the wave function [71] where Y lm l (x) are the spherical harmonics and u l (r) is the reduced wave function, fulfilling the radial Schrödinger equation where U I (r) = U I ( r) is the central potential with isospin I. This equation is indeed regular at the origin 8 u l (r) → r l+1 (26) and it satisfies the asymptotic scattering condition at infinity Thus, the partial wave expansion for the quantum mechanical scattering amplitude with isospin I in the CM system is defined by:

D. Inverse scattering problem
Although we will be determining the potentials from fits to phase shifts, it is worth reminding that the inverse scattering problem allows one to determine a local and continuous potential directly from scattering data by solving for each partial wave either the Gelfand-Levitan or Marchenko equations (see for e.g. [72] for a review). It can be shown that for holomorphic S-matrix functions both methods yield the same local potential. While usually the discussion is conducted within a non-relativistic setup, according to our discussion above, the analysis can directly be overtaken and interpreted at the relativistic level.
This inverse scattering approach was adopted in [73], where a holomorphic S-matrix was used to parameterize the scattering data. In that work it was found that the S and P-wave potentials have a range around 0.25 fm with strengths between 100 − 200 GeV. Quite remarkably they also found a barrier in the isoscalar S-wave and a repulsive core in the isotensor S-wave. While this is a very insightful and mathematically rigorous approach, this method requires exact knowledge of the phase shifts at all energies. In practice, a meromorphic function is fitted up to a maximum energy corresponding to a maximum momentum p max . As we will see below, this puts in practice a limitation to the resolution ∆r ∼ 1/p max with which the potential V (r) may be determined, so that a suitable coarse graining makes sense.

IV. THE CHIRAL ππ LOCAL POTENTIAL
In this section we outline the perturbative matching procedure between quantum mechanic (QM) and quantum field theory (QFT) calculations in order to determine the local and energy dependent chiral potential. The connection between the QFT and QM scattering amplitudes is given by The potential appearing in this equation will be determined in perturbation theory. For the quantum mechanical problem we have the Born series where r 12 = | r 1 − r 2 |, p and p are the initial and final CM momenta, respectively, and q = p − p = 2 p sin (θ /2) is the momentum transfer. The potential U(r) is directly defined from the two-particle irreducible states included in the scattering amplitude. We will define the potential through the tchannel exchanges of the amplitude, so that crossing symmetry will be incorporated exactly when simmetryzing the partial wave expansion 9 . Moreover, in a coordinate space description, contact terms are irrelevant as long as the field theoretical potential is not extended to the origin r = 0 since withP(q 2 ) a generic polynomial in q 2 . Thus, any polynomial part of the scattering amplitude gives a vanishing contribution to the long-range piece of the potential. Therefore, we will analyze only the effect of pion loop contributions on the tchannel.

A. Leading order
We will first discuss the lowest non-trivial order since it provides just contact terms. In the Born approximation, i.e. just taking the first term in (30), the scattering amplitude just becomes the Fourier transform of the potential [71] 9 We have checked that one can either work either in the particle or the isospin basis, the resulting potential is the same.
where q = 2p sin (θ /2) is the momentum transfer. This equation can be inverted to give so in the Born approximation, the scattering amplitudes can be related to the potential by: where T I (s,t) B denotes the disconnected part of the amplitude, i.e. contact terms and t-channel exchange. In the same way, the potential (defined in spatial coordinates) is defined from the disconnected part of the amplitude by: Using the χPT lowest order amplitudes [13] we get These algebraic manipulations are purely formal, and in fact is unspecified what is the meaning of solving the wave equation with these highly singular potentials. Already at this level, we can see the need of introducing a regularization 10 .

B. Next-to-leading order
The NLO contribution becomes more cumbersome. Firstly, we take the potential to be expanded as so we get the matching condition Due to the Dirac delta functions in the LO potential (36), we have a divergence for the real part of e ipr 12 /r 12 , albeit it can be absorbed in the real part of the NLO potential U (4) (r, s). Besides, the non-polynomial pieces in T (4) amplitude corresponding to the t-channel exchange can generally be written as where P(s,t) is a polynomial in both t and s, which analytical expression can be read from (A3), and J(t) denotes the oneloop 2π function. In order to integrate this amplitude, we will take advantage of the analytic structure of the loop function J(t), which is analytic in the whole complex plane but for a cut above 4m 2 π with a discontinuity, disc J(t) = 2 i Im J(t) = 2π iσ (t)/16π, with σ (t) the phase-space factor defined in (5). Thus, up to subtractions one finds the dispersion relation where C.T. is a subtraction constant that can be fixed by setting the value of J(4m 2 π ) = 1/8π 2 . Likewise we have Thus, taking into account the Yukawa integral and the inversion of (38), the NLO potential becomes where t = µ 2 and the spectral function ρ I (µ, s) is defined as with P I (s, µ 2 ) polynomials in s and µ 2 of fourth degree (see Appendix A) and C.T. map into contact terms, which are distributions at the origin and hence vanish elsewhere. All necessary integrals can be obtained from the general integral valid for r > 0, with K n (x) the modified Bessel function of order n and Γ(z) the Euler's Gamma. Polynomials in µ can be generated from derivation with respect to r. The chiral potentials obtained directly from the spectral representation (43), read then From a more general point of view, these potentials play for the ππ-system the role of relativistic van der Waals interactions (see [74] for a review in the atomic case) and hence display their characteristic features: they are attractive and diverge at short distances as and have the expected exponentially suppressed long distance behavior ∼ e −2m π r , namely U 0 (r, s) = − 23m 9/2 π e −2m π r 256π 5/2 f 4 r 5/2 √ s + . . .
In Fig. 1 we show the threshold combination V I (r, 4m 2 π ) ≡ m π U I (r, 4m 2 π ) and, as we can see, they are attractive at all distances. The energy dependence generates a repulsive effect for increasing values of s, i.e.
On the lattice, the energy dependence of the potential is generated from the Nambu-Bethe wave function [52]. As already stated in the introduction, the ππ potential in the I = 2 channel has been computed on the lattice by the HAL QCD collaboration [52,75] for a ≈ 0.12 fm on a 163 × 32 lattice and with a pion mass m π ≈ 870 MeV. For these pion masses the value the chiral potentials in (46) become smaller than for the physical case depicted in Fig. 1. In addition, the HAL QCD lattice potential presents a repulsive core below 0.5 fm. This is a feature one can not obtain using the chiral potentials in (46) which display strong short distance singularities. This fact already suggests that they can not be used at arbitrary short distances. At this point it is worth stressing that both the lattice as well as the present approach based on chiral perturbation theory assume point-like sources, a unrealistic feature.
In the next sections we analyze this topic in more detail.

V. RENORMALIZATION
The renormalization of non-perturbative amplitudes is a tricky matter, particularly with the highly power-divergent kernels deduced from χPT (see for e.g. [38,39] for a discussion within the Bethe-Salpeter framework in momentum space). The chiral potential deduced in coordinate space by a perturbative matching procedure presents an energy dependence. In this section we show how the scattering amplitude stemming from the iteration of the two-pion exchange (TPE) chiral potentials in (46) can be renormalized from a coordinate space point of view if the energy dependence is ignored by taking, say, the threshold value √ s = 2m π . Hence, we will implement as renormalization conditions the scattering amplitude at threshold. We will see that, while this is a mathematically viable approach, it fails phenomenologically. Furthermore, the consideration of energy dependence will prevent a sensible non-perturbative renormalization procedure. For large values of the coordinate space cut-off r c 1.2fm, the results will not be affected by taking either U I (s, r) or U I (4m 2 π , r).

A. Discussion
One of the advantages of the energy dependent coordinate space representation of the potential is that off-shell and field reparameterization ambiguities manifest as contact interactions at the origin. Thus, they reflect the cut-structure of the amplitude, which is hence unambiguously defined. This is unlike their momentum space counterpart, where both poly-nomial and cut contributions are treated on equal footing [39].
On the other hand, a difficulty with the chiral ππ potentials in the previous section is that they become singular at short distances. Hence, the solution of the Schrödinger equation is not well defined in a conventional sense, since the short distance behavior is not dominated by the centrifugal barrier and the regular solution given in (26) is not suitable. An early review on the subject can be found in [76]. Singular potentials are commonplace within EFT and finite solutions exist in a renormalization sense within well specified conditions, as discussed at length in the NN scattering case [77][78][79]. Applications for αα-scattering [80,81] and atom-atom scattering [82,83] are well documented by now (see for e.g. [84] for a sucint and pedagogical presentation). While these renormalized solutions represent theoretically a viable solution to the problem, we will consider here a more phenomenological interpretation by introducing a short distance cut-off r c , which value reflects short distance effects not taken into account in the derivation of the chiral ππ potential 11 . This leaves undefined the short distance dynamics.
The energy dependence of the potential takes into account retardation effects. This can be seen as follows; if the potential is given as a function of the difference of two space-time explaining the 1/ √ s factor in (46). In the case of "heavy pions" or very low energies p m π , we can take √ s ∼ 2m π and work within a non-relativistic approximation. A compelling consequence of causality is the verification of dispersion relations in the complex energy plane. We will dedicate Section VIII to prove that the right analytical properties hold for the partial wave amplitudes.
The spectral representation (40) suggests the use of a spectral regularization consisting of introducing a cut-off Λ at a given value of µ, so that the potential at short distances becomes regular. We find that for Λ > 5m π the regularization quenches the potential for r < 1/m π . We explore this issue in more detail in Sect. VI C.

B. The short distance cut-off and boundary condition regularization
One way to analyze the range of validity of the chiral potentials is to discuss the zero momentum scattering, p = 0. Thus, we solve the Schrödinger equation using the asymptotic solutions at zero momentum (we discard the isospin label here), Color online: zero momentum integrated-in wave functions (solid) from large distances using the scattering length as input and the chiral 2π exchange potentials V ππ at threshold √ s = 2m π as a function of the distance for I = 0 (left) , I = 1 (middle) and I = 2 (right). We also draw the asymptotic zero energy wave function (dashed) for comparison.
with α l ≡ − lim p→0 δ l (p)/p 2l+1 and integrating inward. The result is illustrated in Fig. 2 where we see that the zero momentum wave function presents oscillations at short distances. This behavior can qualitatively be understood if we write the potential at short distances in the form where R ∝ 1/ f is the van der Waals scale, which in our case and from (47) takes the values R = 0.94, 0.68, 0.79 fm for I = 0, 1, 2, respectively. Actually, for r R the centrifugal barrier, l(l + 1)/r 2 , can be neglected and a semi-classical approximation holds since λ (r) 1, where λ (r) is the local wavelength, λ (r) ≡ 1/k(r) = 1/ −U(r), with k(r) the local wavenumber. Thus, one finds λ (r) = 7/2 (r/R) 5/2 1 for (51). In this case, the WKB wave function reads [71] u(r) ≈ 1 [k(r)] where ϕ is an arbitrary phase, which is fixed by the long distance solution.
According to the oscillation theorem [71], the number of nodes of the wave function at zero momentum corresponds to the number of bound states. However, as we know, there are no bound states in the ππ scattering case. Thus, if we want to avoid these spurious solutions, we cannot remove the cut-off completely, but it should be larger than the outmost right node of the wave functions depicted in Fig. 2, which are located at r ∼ 0.5, 0.4, 0.4 fm for I = 0, 1, 2, respectively. In practice, we will take a cut-off slightly above the van der Waals scale. Note that there is a priori no other reason to discard the chiral potential down to these scales. In fact, phase shifts in any partial wave δ l (p) are convergent when the short distance cut-off goes to zero, provided the scattering length α l is fixed. This corresponds to a renormalization program already developed in previous studies [77][78][79][80][81][82][83][84] that will not pursued any further here. The upshot of these considerations in the ππ case is that generally one needs a cut-off r c > 0.5 fm to prevent the appearance of unphysical bound states generated by the TPE potential. However, for these kind of short distance attractive singularity U(r) ∝ −1/r 7 finite r c effects are minor in physical observables if the energy dependence of the potential is neglected. Take for instance the effective range defined in (10) and depicted in Fig. 3. As one can see, the chiral potential and the scattering length lead to a finite result at short distances, r < R. Up to minor oscillations, it provides values which differ from the experimental ones. In Section VI, we will see that this turns out to be much smaller than the elementarity radius, r e ∼ 1.2 fm. Moreover, there is no finite cut-off r c which reproduces the experimental values r 00 = −8.08 fm and r 20 = 131.4 fm.
In the previous discussion the energy dependence of the potential was neglected. On the one hand, if we take the energy dependence into account, we see in Fig. 3 a quite different trend at short distances, namely the effective range is not convergent when the cut-off is removed, i.e. for r c → 0. On the other hand, we also see that for r c 1.2fm this energy dependence in the chiral potential becomes irrelevant.

VI. EFFECTIVE ELEMENTARITY
As we have seen, even at the perturbative level we must introduce a short distance cut-off, r c . In this section we elaborate on sensible choices of this cut-off on the light of the onset of effective elementarity and its corresponding radius r e . This is the scale above which particles interact as if they were pointlike. Thus, they can be taken as elementary, so that a hadronic field can be attached to the particle.

A. Hadron sizes and form factors
Hadrons have a finite size, which is usually characterized by their form factors and corresponding radii. Roughly speaking, we expect that for two hadrons of size r 1 and r 2 they will not overlap at relative separations above the mean average distance (r 1 + r 2 )/2, and their interaction will correspond to that of two elementary particles. However, hadronic sizes obtained from, say, electroweak form factors are specific on that particular process, as we will illustrate below. In fact, in order to describe ππ interactions we are interested on the corresponding effective elementary size, r e , as seen by the interaction. Un- , as a function of the distance for I = 0 (left) and I = 2 (right) It is obtained using the scattering length as input and the chiral 2π exchange potentials V ππ at threshold √ s = 2m π . In addition, we also plot the energy dependence case (blue-dashed), which leads to a divergent value at r → 0.
fortunately, without a microscopic calculation, there is no way to know this size in the case of strong interactions. Nonetheless, for our discussion on the relevant scales it is important to estimate a priori the separation distance r c between the inner and outer pieces of the potential, as this determination has an impact on the minimal number of fitting parameters. Quite generally, this separation distance should be larger than the elementarity radius, r e ≤ r c , and the optimal choice would be to take both distances equal.
Since the pion size is around 1 fm, it is natural to expect that the inelastic non-perturbative description of the interaction will take place in the region below 2 fm.
In order to illustrate our point, let us consider the implications of elementarity for electric and gravitational interactions between pions. The electromagnetic pion form factor reads π + (p )|J µ (0)|π + (p) = (p µ + p µ )F em (q) , whereas the gravitational form factor is given by where q = p − p, P = p + p and Θ 2 (q 2 ) stands for the Heaviside step function. In the Breit frame, where there is no energy transfer, form factors can be interpreted as the Fourier transform of a density [85] F(q) = d 3 r e iq·r ρ(r) , so that the charge form factor determines the charge density and the gravitational form factor the mass density of the pion. Quite generally, form factors are matrix elements of local operators between hadronic states. For the case of operators with well-defined J PC quantum numbers, a generalized meson dominance is expected to work in harmony with the high energy behavior deduced from QCD counting rules (see for e.g. [86] for a thorough discussion and comparison with lattice QCD data). For example, in the electromagnetic case, it can be parameterized according to vector meson dominance as, where we have kept only the ρ meson with m ρ = 0.77 GeV, as we are only interested in the long distance properties.

B. Point-like vs extended particles interactions
According to the previous discussion, the electrostatic potential can be written as which is depicted in the left panel of Fig. 4 and reflects that for r > 1.2 ∼ 1.5 fm, pions start to interact as expected from point-like particles 12 . In the gravitational case, a good description is found with the tensor f 2 (1270) meson [86] with the mass scale given by m f = 1.2 GeV, so that 12 The value at r = 0 corresponds to the electromagnetic mass difference between charged and neutral pions, V el ππ (0) = ∆m π | EM = m π + − m π 0 | EM = e 2 m ρ /2 = 2.8 MeV, which provides a reasonable value. In this case, right panel of Fig. 4, the interaction becomes elementary at r e ∼ 1 fm, a shorter scale. The previous discussion illustrates our point, namely effective elementarity depends on the particular process.
Before proceeding further, it should be noted that within a more microscopic point of view, for instance a cluster quark model, these are not the only possible contributions to the interactions between pions (see e.g. [87][88][89]). Actually, in a Hartree-Fock approximation they would correspond to the direct interaction (Hartree) term. In addition, one also has the exchange (Fock) term where the quarks inside different pions are interchanged. These terms are genuinely non-local at short distances, so that they are exponentially suppressed at long distances. Therefore, they are expected to contribute below the elementarity radius, r e , and hence can be regarded as finite size effects.

C. Spectral regularization
While we might estimate the elementarity radius for the chiral potentials from the corresponding folding of "strong" densities, we prefer to analyze instead the effect of introducing a cut-off Λ in the spectral function in (43), i.e.
which corresponds to the two-pion invariant mass spectrum. In Fig. 5 we show the ratio U I Λ (r)/U I (r) for spectral cutoffs of Λ = 1.4, 1.12, 0.84 GeV at threshold √ s = 2m π as a function of the distance. This estimate yields a larger elementarity radius, which becomes r e = 1.2 fm for the largest spectral cut-off and a quenching of about a half for Λ ∼ m ρ . This numerical exercise shows that the naive estimate r e ∼ 1/Λ is numerically rather inaccurate. Of course, one could consider Λ as a fitting parameter in the analysis of ππ scattering. Nevertheless, in our view this has the disadvantage that a model dependence is introduced by the cut-off procedure above the elementarity radius.
Note also that the larger r e the smaller the TPE potential, since U(r e ) = O(e −2m π r e ), so that one may end up a the situation where a model independent treatment becomes only possible when the chiral contribution is actually vanishingly small. In the next section we will present a different strategy. Overall, our numerical results will confirm this pessimistic expectation.

D. Potential separation
The previous discussion suggests that we should decompose the potential for each isospin channel as U I (r) = U I Short (r)θ (r c − r) +U I Long (r)θ (r − r c ) , (60) where U Short (r) is a short distance and U Long (r) stands for the long distance contribution. The natural choice is to take for the long-range part the unregularized chiral TPE potential i.e. potential computed in χPT to a given order in the chiral expansion. As we have seen in the previous section, the lowest order effect in the chiral potential comes from the two pion exchange, which gives a contribution that at long distances is U χ (r) = O(e −2m π r / f 4 ). For r c ∼ 1/m π this contribution will in principle play a role. Of course, these most peripheral contributions to the interaction will contain perturbative corrections to all orders in m π / f , which are expected to modify slightly the tail of the potential. Therefore the question is what is the relative importance of the unknown short distance piece and the known long distance contribution. Finally, let us mension that while we will keep the energy dependence already present in the chiral potential we will assume for simplicity that the short distance potential is energy independent, as long as the inelasticity is negligible (see also Section IX).

VII. COARSE GRAINING
The basic notion of coarse graining in the scattering problem was outlined by Avilés as early as 1972 in an insightful but forgotten paper [90] within the context of NN interactions. In this article the potential was effectively represented as a sum of delta-shell potentials. This form has important simplifications and a recent comprehensive mathematical study of this specific case has been carried out [91]. Here we extend the approach to account for the χPT potential tail at long distances (see also Refs. [32][33][34] for a paralell treatment of the NN case).

A. Short distance potential
The basic idea is as follows. If we want to describe the two-particle CM wave-functions limited to the range ∆p, only gross information can be determined in an interval ∆r, with ∆r∆p ∼ 1. Thus, for a maximal s max = 2 GeV 2 , we have ∆p ∼ s max /4 − m 2 π ∼ 0.70 GeV and one obtains ∆r ∼ 0.3 fm. This uncertainty suggests that for a limited energy range the potential only needs to be known in a limited number of points. With this in mind, we consider for r < r c the ππ potential as a sum of a subsequent number of δ functions separated by about ∆r. Thus, up to r c , it requires to introduce N δ delta shells 13 Thus, we have two relevant scales in our setup: the separation distance r c and the coarse graining scale ∆r. While one expects the results and main properties of the potential to be independent of the particular choices of r c and ∆r, this can only happen when accurate information on the interaction is available at all energies. In our case ∆r acts as an UV regulator whereas r c works as an IR regulator. We will analyze this problem for four different values of r c . For convenience, we will set r c = 0.9, 1.2, 1.5 and 1.8 fm. The solution of (25) for the delta-shell potential in (62) is straightforward. One has u l,n (r) =ĵ l (pr) − tan δ l,nŷl (pr) , r n < r < r n+1 , whereĵ l (x) = x j l (x) andŷ l (x) = x y l (x) are the reduced Bessel functions of first and second kind, respectively, and δ l,n is the accumulated phase shift. The discontinuity in the logarithmic derivative at r = r n with r + n ≡ r n + 0 + and r − n ≡ r n + 0 − leads after using the unit Wronskian conditionĵ l (x)ŷ l (x) −ŷ l (x)ĵ l (x) = −1 to the bilinear recursion relation for tan δ l,n tan δ l,n+1 = A l,n (p) + B l,n (p) tan δ l,n C l,n (p) + D l,n (p) tan δ l,n , with A l,n (p) = ∆rU(r n ) j l (pr n ) 2 , B l,n (p) = ∆rU(r n ) j l (pr n )y l (pr n ) + k , C l,n (p) = −∆rU(r n ) j l (pr n )y l (pr n ) + k , D l,n (p) = ∆rU(r n )y l (pr n ) 2 , with the initial and final conditions δ l,0 (k) = 0 , δ l (k) = δ l,N (k) .
As shown in [92], these equations can be interpreted as the discrete version of the variable phase equation of Calogero [93], corresponding to the limit ∆r → 0. Of course, the previous equations define an integration method in this limit. We stress that the idea of coarse graining is that if U(r) is determined from data with a maximum CM momentum p max , the natural resolution of the problem is ∆r ∼ 1/p max , and U(r n ) are the natural fitting parameters 14 .

B. Fitting procedures
In this section we present our numerical results based on standard χ 2 -fits. Precise ππ-scattering phase shifts have been obtained in [21,23] using Roy equations up to √ s max = 1.42 GeV and we will use their tabulated values to determine our fitting parameters. While the standard strategy in ππ scattering studies has been to use the physical region √ s > 2m π , for reasons to be justified in Section VII F we also include in the S0 fit the subthreshold region, 0 < √ s ≤ 2m π , also deduced in those Roy-equation analyses [21,23].
Anticipating the result, we will carry out the study, first without including the chiral tail, since as we will see in Sect. VII G, it plays a minor role in the resulting fitting parameters λ n .

C. Energy independent coarse graining
To start with, we will assume on purpose sufficiently large separation distances so that the long-distance field theoretical contribution O(e −2m π r ) can be safely neglected.
In order to constrain as much as possible the delta-shell coefficient values λ j , we will introduce each delta-shell adiabatically, one by one. In addition, we will impose at threshold the corresponding ππ scattering lengths and slope parameters, so that the value of the most internal delta-shell parameters, λ 0 and λ 1 , are completely fixed. The ππ S0-, S2and P-wave phase shift obtained in this way for the different values of r c previously chosen are plotted in Fig. 6, whereas the value of the delta-shell parameters are given in Table I. The maximum fitted energy √ s for each partial waves is chosen as the maximum one for which one can find a χ 2 /d.o. f around 1. In addition, as we will discuss in Section VII F, we also include for the S0 fit the subthreshold region m π < s < 2m π .
We can see from these results that with a few free parameters, from 3 to 6 depending on the particular choice of r c , IJ r c (fm) a 0 (m π ) −1 b 0 (m π ) −3 00 0.  this formalism already allows one to describe the ππ scattering in the elastic region. In the case of the S2 wave, where the inelasticities are very small at low energies, it is possible to obtain a perfect description up to the maximum energy at s max = 2 GeV 2 . For the P wave one obtains a good description, up to energies around theKK threshold at √ s = 1 GeV, whereas for the S0-wave the description is limited to the region around 0.9 GeV, where the phase presents a huge rise due to the effect of the inelastic f 0 (980) resonance.
Our results reproduce qualitatively features found in [73] using inverse scattering methods, where by construction the potential is a local and continuous function. In particular, for the S waves we find a short distance barrier in the isoscalar and a repulsive core in the isotensor. Of course, we only get a coarse grained version of those potentials, which befits our idea of a finite resolution sampling.

D. Threshold parameters
Our results for the threshold parameters defined by (8) are presented in Table II. As expected from the quality of the fits, they agree within uncertainties with the results obtained in previous analyses [14][15][16][17][18][19][20][21].

E. Resonance poles
Resonances are determined by looking for poles in the second Riemann sheet, which is defined by the complex wavenumber k R = k r + ik i with k r ≡ Re k R > 0 and k i ≡ Im k R < 0. This can be achieved directly by substituting k → k R in the Schrödinger equation (19), or its coarse-grained delta-shell implementation (65) and (66).
Our numerical results are presented in Table III. The numerical values are slightly different from those quoted in benchmarking studies based on Roy equations and forward disper-  [43]. This is not fully surprising since the analytical properties of our scattering amplitude are not the same as in the analyses based on the Roy equations, but only to leading order in the chiral expansion and hence the extension to the complex plane is not determined from the phase-shift analysis only. Nonetheless, we find that the result is encouraging and we expect to return in the future in order to implement higher orders to analyze the effect.

F. Crossing and comparison with Roy equations
As we have mentioned in the introduction, there are no direct ππ scattering data. The closest thing to it are possibly the outcome of Roy equations, which incorporate by construction crossing, analyticity, unitarity and Regge behavior [16,17,21]. Moreover, a proper identification of the LECs requires by definition the satisfaction of crossing, particularly in the sub-threshold region. For definiteness, we show in Fig. 7 the real part of the amplitude starting at the edge of the left hand cut s = 0, covering the subthreshold region 0 < √ s ≤ 2m π (where the amplitude is purely real) and the physical elastic region, 2m π < √ s < 2m K , where the imaginary part of the amplitude is determined from unitarity and the corresponding real part. Note that we are neglecting possible inelasticities coming from multi-pion channels. In particular, one can clearly see the Adler zeros of the S0 and S2-wave amplitudes, a characteristic features of the subthreshold region and a distinct fingerprint of chiral symmetry [11].
The comparison in Fig. 7 supports the view that fits to phase shifts in the physical region, even to relatively high energies, do not constrain the subthreshold region. This has a large effect on the location of the σ -resonance and the implications will be discussed in more detail elsewhere. In addition, our fits have not imposed crossing correlations and the consideration or not of the subthreshold region in the S0 channel can be grasped by comparing the location of the corresponding resonance. If the subthreshold region m π < s < 2m π is not imposed in the fit, one gets substantially smaller values for the real part. For instance, for r c = 0.9 fm and N δ = 3 delta-shells one gets √ s σ = (324 − i 194) MeV, whereas the pole moves to √ s σ = (453 − i 212) MeV when the subthreshold region is also fitted. This large influence of the subthreshold region in the S0 wave is not surprising, since the pole is rather far from the real axis. In contrast, the location of the ρ-resonance pole proves rather insensitive to the subthreshold region in the P-wave. Although we are not fitting the subthreshold region in this channel, one gets for r c = 0.9 fm and N δ = 3 deltashells √ s ρ = (762 − i 70) MeV, even when the subthreshold extrapolation of the fit, middle panel of Fig. 7, shows a significant discrepancy with the outcome of the Roy equations. Finally, we also report a discrepancy in the isotensor channel, where the phase shifts are compatible up to s max . Actually, in all cases the amplitude has a zero at s = 0, due to the factor √ s in the definition of t IJ (s). Thus, the failure to satisfy the subthreshold behavior in the S2 wave is intrinsic 15 . We have checked that this does not improve by incorporating the TPE tail of the chiral potential (see also next subsection). The fact that the subthreshold region is quantitatively as relevant as the physical region, is another manifestation of the relevance of crossing in ππ scattering, and calls for improvement when all relevant partial waves are considered.

G. Full potential: Inclusion of Two Pion Exchange
Once the chiral tail has been defined, we can finally construct a full potential. At short distances (r < r c ), the nonperturbative regime of QCD is encoded by a sum of deltashells separated by the minimum De Broglie wave-length considered in the analysis. At large energies (r < r c ), the chiral potential constructed in (46) with the correct analytical properties and left-hand cut contribution is included. Higher order 15 This is in common with other unitarization methods, such as IAM where the unitarized amplitude t IJ (s) = t   corrections in the chiral expansion are O( f −6 ) and the determination of the corresponding potential requires going in the quantum mechanical picture to second order in the Born approximation. Their analysis is left for future research.
The resulting fits are reported in Table IV. As one can see, and with the exception of the S2 wave with three delta shells, we do not find a substantial improvement due to the explicit IJ r c (fm) a 0 (m π ) −1 b 0 (m π ) −3 00 0. incorporation of the TPE potential. A similar trend is found for the corresponding low energy threshold and resonance parameters, see Tables V and VI. In any case, chiral effects due to explicit TPE are generally found to play a minor role. In our analysis we have restricted to the lowest Sand Pwaves. The implementation of higher partial waves is straightforward and requires introducing further short range contributions. Due to the centrifugal barrier term, we expect that the number of grid points will be reduced as the angular momentum increases.

VIII. ANALYTIC PROPERTIES AND THE N/D METHOD
Traditionally, much of the discussion of ππ scattering has been marked by analyticity and dispersion relations, which at the partial-wave level and in the unsubtracted case read where 0 < s n < 4m 2 π are the possible bound states, which we keep for generality, and Imt IJ (s) = ρ(s)|t IJ (s)| 2 for elastic scattering. In this section, we will discuss this issue within the context of our coordinate space framework and for the specific chiral potential derived in Section IV, but the results are general. Actually, we will see that the subtractions are not explicitly needed, although they are somewhat encoded into the short distance component of the potential. While many of the issues discussed here have been known for potential scattering since many years [94,95] (see also [9,10]), we feel it is necessary to review the main aspects for completeness.
The analytical properties of the scattering amplitude in the complex energy plane are determined by the long-distance behavior of the interaction. This is the basis for dispersion relations, which own their popularity to their link to axiomatic field theory and its straightforward implementation in terms of leading singularities. A frequent method, which has been used in these regard to implement known analytical properties, is the so-called N/D approach. In the N/D approach the partial wave amplitude is written in the form [96] with N IJ (s) having only left-hand cut singularities and D IJ (s) having only right-hand cut singularities. This method has often been invoked in ππ scattering and implemented in several approximations (see for e.g. [97]).

A. Jost functions
The way of realizing the N/D representation in potential scattering is well-known. We will describe here the formalism within the coarse grained approach for completeness as well as to provide the N l (k) and D l (k) functions explicitly, (we drop the isospin index for simplicity). The discussion is naturally carried in terms of the quantum mechanical scattering amplitude defined in (7), f l (k) = 2t l (s)/ √ s, as a function of the CM momentum, k, with s = 4(k 2 + m 2 π ). Note that s is invariant under k → −k and hence two-valued. The discussion below entitles to take Imk > 0 for the first Riemann sheet and Imk < 0 for the second.
For a regular potential, the Jost functions, F l (k), are defined by the regular solutions of the wave equation at short distances, i.e. u l (r) →ĵ l (kr) and subjected to the asymptotic condition at r → ∞ [71] whereĥ (1,2) l (x) =ĵ l (x) ± iŷ l (x) are the reduced Hankel functions fulfillingĥ (1) l (x) * =ĥ (2) l (x) andĥ l (x). Thus, the S-matrix is then defined as so that the scattering amplitude becomes The Jost functions fulfill the reflection conditions for complex k. Furthermore, due to (74), it can be shown [9,10] that the functions fulfill the relations n l (k) * = n l (−k * ), n l (k) = n l (−k), which means that n l (k) is purely real for Im k = 0. Moreover, Re n l (k) = Re n l (−k * ) and Im n l (k) = −Im n l (−k * ), so that n l (k) is purely imaginary for Re k = 0. In addition, Thus, they have the desired properties for a potential constructed as a superposition of Yukawa potentials. A straightforward consequence of these properties is Levinson's theorem. While the proofs of these statements have been known for a long time [94,95], to our knowledge they have not been considered within the present context. We will review them adapted to our complete potential, which can be decomposed into two pieces: a cut-off potential and a Yukawa superposition. While they are discussed separately in the literature, our case at hand involves both cases at the same time. We will show next the pertinent steps.

B. Analytic properties
The scattering amplitude obtained perturbatively from χPT enjoys analytical properties deduced from the corresponding Feynman diagrams, i.e. particle exchange generated at the partial-wave level a left-hand cut, which discontinuity has been used to reconstruct the chiral potential. A relevant question is whether our resulting full amplitudes obtained by solving the Schrödinger equation, share these properties beyond first order in perturbation theory, i.e. whether they satisfy, as expected, dispersion relations.
In order to check the analytical properties of the quantum mechanical amplitude, we note that our chiral potential is indeed a superposition of Yukawa potentials with the exception of the explicit s-dependence in the spectral function. It is convenient to write the potential as a superposition of exponentials in the form where and ρ I was defined in (44). This is indeed our case for the TPE potential in ππ scattering, as one can see comparing (43) and (78). Being a Laplace transformation, it corresponds to the so-called analytical potential in the complex-r plane for Re r > 0 [98], which fulfills the relation lim ρ→0 U I (ρe iθ , s) = 0 for −π/2 < θ < π/2, with r = ρe iθ . As we will see below, while the exponential falloff of the potential at long distances suffices to prove the analyticity in the strip |Im k| < m π , with k = s/4 − m 2 π , the spectral decomposition is indeed needed to establish the cut along the line Re k = 0 and m π < Im k of the S-matrix. There are several versions of the proof. On the one hand, it can be proved by estimating a bound for the Jost function using directly the spectral representation. On the other hand, by profiting from the analytical character of the potential and deforming the integration in r into the complex plane so that kr > 0. For completeness we will review next and in a sketchy fashion the second method, as it does not require a bounded spectral function or the use of an spectral regularization (see Section VI C).
The determination of the analyticity domain of the quantum mechanical problem is based on the equivalent integral equation for the Jost Functions as follows u k,l (r) =ĵ l (kr) + r 0 dr K k,l (r, r )U(r , s)u k,l (r ), (80) which is a Volterra type integral equation and the kernel is given by Taking the large-r limit and comparing with (71), we get This equation is the basis for the analytical continuation to the complex-k plane. Actually, in the limit r → ∞ one has from (71) that |h (1) l (kr)| = O(e −Im k r ) for Im k > 0. Thus, even when |u k,l (r)| = O(e Im k r ), one has a finite integral provided that the potential goes to zero. As a consequence, F l (k) is analytical for Im k > 0. On the contrary, for −m π < Im k < 0 one has |h (1) l (k r)U(r, s) u k,l (r)| = O(e −2Im k r−2m π r ) which is convergent for Imk > −m π . Therefore, F l (k) is analytical for Imk > −m π and hence F l (−k) is analytical for Im k < m π . In conclusion, S l (k) = F l (−k)/F l (k) is analytical in the strip |Im k| < m π .
For k = |k|e iθ , due to the analytical character of U(r, s), we can deform the contour to r = ρe −iθ , so that, for any term in the spectral integral over µ, one has a pole at k = −iµ/2. In the case of F l (k) and F l (−k), this pole becomes a cut after µ-integration along the lines Re k = 0 and −∞ < Im k < −m π and m π < Im k < ∞, respectively. The fact that F l (k) is analytical for Im k > 0 and lim k→∞ F l (k) → 1 allows one to write a dispersion relation in the upper half-circle where the antisymmetry of ImF l (−k) = −ImF l (−k) from (74) has been used. Thus, passing to the variable s = 4(k 2 + m 2 π ), which is one-valued in Imk > 0, one can define the function 16 so that we can identify D l (s) = F l (k) for Imk > 0 where Im D l (s ) = Im F l (s ) has a right-hand cut in 4m 2 π < s < ∞ and is real for s < 4m 2 π . Furthermore, in the elastic approximation, using (6) we get where is real for s > 4m 2 π . For real s < 0 we get that s + i0 + ↔ 0 + + iκ with κ > m π and hence Thus, the discontinuity is where we have used that for Imk > 0, F l (k) is analytical and hence F l (0 + + iκ) = F l (−0 + + iκ). Finally, note that the s-dependence appearing in the spectral function does not spoil these analytic properties. This completes the proof that the scattering amplitude obtained by solving the Schrödinger equation for the potential in (60) has the correct analytical properties.

C. Coarse graining
In order to calculate the Jost functions, in practice we coarse grain the interaction using the delta-shells and we proceed as before. In the discretized version we define the accumulated Jost functions F l,n (k) and F l,n (−k) as u l,n (r) = 1 2 F l,n (k)ĥ (2) l (kr) + F l,n (−k)ĥ (1) l (kr) , r n < r < r n+1 .
However, we keep track of both the continuity of the wave functions and the discontinuity of the derivative separately, with r + n ≡ r n + 0 + and r − n ≡ r n + 0 − . Again we use the fact that the Wronskian 2i(ĥ (1) l (x)ĥ (2) l (x) −ĥ (1) l (x)ĥ (2) l (x)) = 1 so that F l,n+1 (k) = A l,n (k)F l,n (k) + B l,n (k)F l,n (−k) F l,n+1 (−k) = C l,n (k)F l,n (k) + D l,n (k)F l,n (−k) where we have introduced the coefficients, , , l (kr n ) 2k (94) and the initial conditions are which corresponds to take the normalization u l,0 (r) =ĵ l (k, r).
The final values are For illustration purposes, we plot in Fig. 8 the functions N IJ (s) (which is real) and D IJ (s) (which is complex) separately above the threshold for the five delta-shell case. Of course, these functions reproduce the phase shifts presented before. As expected, Im D IJ (4m 2 π ) = 0. In this representation the Breit-Wigner position of the resonance is given by computing the zeros of Re D IJ (s R ) = 0. While it is not shown in the pictures, the Jost functions display some oscillatory behavior at higher energies due to the explicit delta-shells. Nevertheless, they still go to the expected values D(∞) = 1 and N(∞) = 0.

A. Energy dependent coarse graining
As mentioned in the introduction, at sufficiently high energies particles are produced and elastic scattering happens in the presence of absorption, which one may view as a leak or hole in the probability. The inelastic region is characterized by the recombination of the two-pion internal configuration, which, of course, demands a complex energy-dependent potential. The complexification of the potential can be understood in terms of the loss of probability of the elastic channel. In order to have an idea of the size of the inelastic hole a inel , we show in Fig. 9 the inelastic ππ total cross sections in different isospin channels, defined as If we take σ inel I = 4πa 2 inel , the largest inelastic cross section in Fig. 9 is compatible with an inelastic hole of less than half a fm, i.e. a inel 0.5 fm. Of course, this corresponds according to (97) to the contribution of all partial waves. As we have already seen, this loss of probability at the partial-wave level is parametrized in terms of a momentumdependent inelasticity η J I (s), see (4). The imaginary part of the potential plays exactly the same role. At this point, we will assume that this complexification of the potential can be implemented just in the most inner layer of the potential. This assumption can be justified by analyzing the phenomenological inelasticities as a function of the impact parameter given by the relation b = (l + 1/2)/p, with l the angular momentum quantum number. In order to provide the range of the inelasticity, we show this dependence in Fig. 10. For the S waves it is found that only for the smallest value, b min = 1/(2p max ) ∆r, the inelasticity is η 1. As we see, if we take ∆r ∼ 0.3 we may parameterize the inelasticity by one single energy dependent and complex parameter 17 .
Hence, starting from our elastic previous description, we will assume λ 0 in (62) to be a complex unknown function of the momentum and we will fit again the pseudo data given in [21] for the phase shift and inelasticity in [21] up to the maximum energy value provided √ s max = 1.42 GeV. This procedure allows us to describe each partial wave (phase shift and inelasticities) exactly at each energy point. The results for the four r c -values considered are plotted in Fig. 11, whereas the value of the real and imaginary part of the inner delta-shell layer is depicted in Fig. 12.

B. Traces of analyticity in the inelastic case
An interesting feature which can be appreciated in Fig 12 is the close resemblance of the real and imaginary parts of the inner delta-shell coefficient with the expected behavior from dispersion relations around a pole or a inelastic threhsold. These two effects reflect in the S0 and P channels higher resonances or inelastic channels not explicitly included in the present optical potential analysis. Interestingly, Cornwall and Ruderman [7] found the fixed-r dispersion relation for the optical potential V (r, s) given in (20). The implementation of such a dispersion relation in our analysis would reduce the number of fitting parameters in the inelastic region but it would also require a clear understanding of the high energy behavior. We leave this interesting investigation for future research.
The results for the real and imaginary part of the energydependet inner delta-shell when the chiral tail is included differ from those without, plotted in Fig. 12, only at low energies. The comparison is depicted in Fig. 13 showing again the rather small effect introduced by chiral corrections in the inelasticity parameter.
Another interesting possibility which deserves some further investigation is the generalization to the coupled channel case, to account explicitly for the opening of the KK and ηη thresholds, while keeping multi-pion channels in the inelasticity factor.

X. SUBTRACTIONS, LOW ENERGY CONSTANTS AND THE NUMBER OF PARAMETERS
In our construction of the chiral ππ potential the low energy constants (LECs) have been discarded as they do not contribute for r = 0. In fact, the spectral representation does provide the same coordinate dependent potential regardless on the number of subtractions.
This raises the problem of where is this counterterm information gone within the present approach. In this last section we want to address the relation among subtraction constants in dispersion relations, low energy constants and the number of independent parameters within our coarse graining approach.
is an analytical function of the energy. (see also the discussion below) We warn the reader that we have not succeeded in finding a unambiguous relation, which may ultimately be traced to two aspects. Firstly, it has to do with known ambiguities in mapping different regularization methods, namely the one used in χPT and the coordinate space regularization used here. Secondly, there is a difficulty in separating the short distance parameters in perturbation theory, particularly if we use a nonperturbative resummation scheme to fit the parameters.
In order to elaborate on this issue and appreciate the difficulties, we proceed in perturbation theory and re-write the problem in terms of a Fredholm integral equation, u k,l (r) =ĵ l (kr) + ∞ 0 dr G k,l (r, r )U l (r ) u k,l (r ) , (98) whereĵ l (x) = x j l (x) is a reduced spherical Bessel function of the first kind and G k,l (r, r ) is the Green's function satisfying An analytic expression for the Green's function G k,l (r, r ) can be written in terms of two function ansatzs u(r) and v(r) as Inserting this in (99), it follows that u(r) and v(r) are solutions of the homogeneous equation with a unit Wronskian, i.e., We choose one of the two solutions to be proportional to the regular one,ĵ l (kr). Then the other linearly independent solution with the desired Wronskian has to be proportional tô y l (kr) = kr y l (kr), i.e. the reduced spherical Bessel function of the second kind. Therefore the Green's function of the ordinary differential equation with the proper normalization can be written as where r < = min{r, r } and r > = max{r, r }. For the normalization of (98) the scattering amplitude can be written as which using the perturbative expansion inferred from reitera- Actually, the singularity structure in coordinate space suggests introducing a short distance cut-off r c and hence a short distance potential U Short (r) 18 . That means that in order to identify numerically the counterterms in perturbation theory we 18 Here we assume for simplicity a local form. More generally we may assume a nonlocal form of the type [Uψ](r) = rc 0 U(r, r )ψ(r )dr with U(r, r ) = 0 for r > r c . The coarse graining interpretation below reduces effectively this non locality to a local form within a sampling distance ∆r. FIG. 12. Real and imaginary part of the inner delta-shell potential for the S0, P and S2 partial waves whereas solid-black, green-dashed, red-dotted and blue dot-dashed lines correspond to the results with r c = 0.9, 1.2, 1.5 and 1.8 fm, respectively. The error bands have been computed using at each energy point a bootstrap with a uniformly distributed sample of 1000 points and taking the 68% of their distribution as the standard deviation. may split the integrals as and apparently the matching to the one-loop result could be undertaken in a straightforward fashion. This is actually not so, since this implies disentangling the fitting parameters in a chiral expansion, namely λ n ≡ U(r n )∆r = λ (2) n + λ where we can arbitrarily shift λ (2) n and λ  Tables I  or IV. This situation is not exclusive to the present approach and in fact is common to all unitarization schemes. For instance in the IAM method the fitted LEC's are different than those of FIG. 13. Real part of the inner delta-shell potential for the S0, P and S2 partial waves when the chiral tail is included or excluded. The only differ at low energies. Gray (light-gray), green (light-green), red (orange) and blue (cyan) bands correspond to the results with (without) the chiral tail for r c = 0.9, 1.2, 1.5 and 1.8 fm, respectively.
χPT or the predicted unitarized amplitudes from χPT develop huge uncertainties [67]. The fact that this perturbative matching can ultimately provide a successful description and still fulfilling the condition λ (2) n λ (4) n is expected. We note, however, that the small changes between the λ n parameters corresponding to the case without TPE potential listed in I and the case with TPE listed in IV suggest this hypothesis.

XI. CONCLUSIONS
The optical potential in ππ scattering is a meaningful object under the most common and general assumption of the validity of the Mandelstam double spectral representation. Therefore, it plays a relevant role in the analysis of such an interaction within a invariant mass formulation of the relativistic twobody problem. Contrary to the more employed Bethe-Salpeter equation, such an approach is free from well-documented spurious singularities, which are triggered by incomplete calculations embodying subsets of Feynman diagrams with particle exchange. In addition, it is a much simpler approach in the CM frame, as it effectively reduces to a Schrödinger equation for equal mass particles.
Within such a framework, in the present paper we have analyzed ππ scattering from a coarse-grained point of view at distances smaller than the elementarity radius of the pion. This means sampling the interaction in coordinate space at a resolution of the order of the shortest de Broglie wavelength. In our case, where we choose a maximal CM energy of s max = 2 GeV 2 , the resolution turns out to be ∆r ∼ 0.3 fm. As a results, we obtain successful fits in the S0, P and S2 partial waves with the expected number of parameters. We have also analyzed the role of inelasticities by an energy dependent coarse grained interaction. The implications from chiral symmetry have also been analyzed in terms of a long-distance potential featuring the two-pion-exchange mechanism. This potential has been determined for the first time.
A non-perturbative renormalization of the amplitude, based on boundary conditions in coordinate space, is precluded by the energy dependence of the chiral potential. For a finite short distance cut-off about 1.2-1.5 fm this energy dependence becomes irrelevant.
A somewhat surprising result of our analysis is that explicit chiral corrections play a minor role, since at the distances above 1.2-1.5 fm chiral potentials are almost negligible. A rewarding consequence in this regard concerns the extraction of phase-shifts from energy shifts calculated in a finite box in the relative distance by means of the Luscher formula [46]; its applicability requires the interaction to vanish above a given size which provides a lower limit to the size of the box, L. While the mere O(e −m π L ) nominal estimate suggests L 2fm, our analysis is compatible with taking L ∼ 1 − 1.5fm, a much smaller value.
We have also shown that the quantum mechanical reinterpretation of the problem does not spoil the proper analytical properties of the partial wave scattering amplitude, which are explicitly fulfilled by the Feynman diagrams. Thus, the conventional dispersion relations with both a left-hand cut due to particle exchange in the crossed channel and the right hand cut due to unitarity are fulfilled. The standard N/D decomposition of the amplitude is explicitly realized, albeit without subtractions. We expect subtractions to be encoded in the short distance components of the potential, but the explicit determination and identification of the subtraction constants in terms of the short-distance parameters remains to be accomplished.
An important advantage of the coarse graining perspective is that the number of independent parameters is determined a priori by the shortest wavelength and the available crossing constraints. These constraints become increasingly large as we increase the angular momentum of the partial wave. This point deserves further investigation and the determination of all partial waves with this minimal number of parameters is left for future research.
A traditional objection to the successful data-driven unitarity methods is based on their lack of a power counting scheme, reflecting field dependence and off-shell ambiguities absent in the conventional bona fide EFT framework. The present coarse graining approach to ππ scattering makes no further assumptions than those usually made. Namely, it implements unitarity in the elastic regime and it matches χPT in perturbation theory above a given separation distance, which is estimated to be about 1.2-1.5 fm. In contrast, it does have the advantage that we can estimate the number of fitting param-eters a priori. More demanding fits due to an increase in the number of (consistent) data should not require more fitting parameters, but rather determining short distance parameters more accurately. Besides, the method is quite simple as it parameterizes the unknown short distance behavior regardless of any power counting. This also allows a discussion a posteriori of the chiral contributions which turn out to be minor.