the modern effective field theory ( eft ) analysis of the nn interaction using chiral symmetry as a constraint has a recent but prolific history  @xcite ( for comprehensive reviews see e.g. ref .
most theoretical setups are invariably based on a perturbative determination of the chiral potential  @xcite and the subsequent solution of the scattering problem  @xcite .
actually , the theory encounters many problems in the low partial waves and in particular in the @xmath2waves ( see however ref .
@xcite for a more optimistic view ) .
indeed , there is at present an ongoing debate on how an eft program should be sensibly implemented within the nn context and so far no consensus has been achieved ( see e.g. @xcite ) .
the discussion is concerned with the issue of renormalization vs.  finite cutoffs , a priori ( power counting ) vs.  a posteriori error estimates or the applicability of perturbation theory both on a purely short distance theory or around some non - perturbative distorted waves . at the moment
, it seems fairly clear that an eft scheme with a cutoff - independent and systematic perturbative power counting for the s - matrix ( the so - called ksw counting ) fails  @xcite . on the other hand , the original eft inspired scheme  @xcite ( the so - called weinberg counting )
has recently been shown to produce many results which turn out to be strongly cutoff dependent  @xcite and hence to be incompatible with renormalizability  .
thus , some acceptable compromise must be made .
actually , the chosen approach between this dichotomy depends strongly on the pursued goals and it is fair to say that any choice has both advantages and disadvantages . in any case , the reason for both failures can be traced back to the nature of the long - distance chiral potentials ; while pion - exchange potentials fall off exponentially at long distances they include strong power - law singularities at short distances .
those singularities become significant already at distances comparable with the smallest de broglie wavelength probed in nn scattering below pion production threshold .
obviously , any development of nn interactions based on chiral dynamics will presumably require a deeper understanding and proper interpretation of the peculiarities of these highly singular chiral potentials .
although singular potentials where first analyzed many years ago  @xcite ( for an early review see e.g.  @xcite and for a more updated view within an eft context see  @xcite ) , their short distance singular character within the nn interaction has seriously been faced more recently within a renormalization context for the one - pion exchange ( ope ) potential  @xcite and the two - pion exchange ( tpe ) potential  @xcite .    in the np scattering problem
the @xmath0 channel is very special since the scattering length is unnaturally large as compared to the range of the strong interaction , @xmath3 .
in fact , even at zero energy the wave function probes relatively short distance components of the chiral potential  @xcite  \ ] ] where @xmath4 is the zero energy wave function , fulfilling the asymptotic condition @xmath5 . most of the integrand is located in the region around @xmath6 which is in between ope and tpe ranges . moreover , the low energy theorem of ref .
@xcite allows to write @xmath7 which in the extreme limit @xmath8 yields @xmath9 .
numerically it is found that @xmath10 is far more dependent than @xmath11 and @xmath12 when being evaluated at lo , nlo and nnlo . ] .
higher energies become even more sensitive to short distance interactions . consequently , this channel looks like an ideal place to learn about the size of the most relevant short range corrections to the nn force in the elastic scattering region . actually , in the @xmath0 channel , most eft inspired schemes yield at leading order ( lo ) ( which consists of ope plus a nonderivative counterterm ) an almost constant phase of about @xmath13 around @xmath14 , and an effective range of @xmath15 . on the other hand
, most determinations from partial wave analysis  @xcite and high quality potential models  @xcite yield an almost vanishing phase at this center - of - mass ( cm ) momentum while the experimental effective range is about twice the ope value , @xmath16 .  .
] it is quite unbelievable that such large changes can be reliably accommodated by perturbation theory starting from this lo result despite previous unsuccessful attempts treating ope and tpe perturbatively  @xcite .
actually , for the singlet channel case , short distance components are enhanced due to the large value of the scattering length and the weakness of the ope interaction in this channel .
this is why tpe contributions have been treated with more success in a non - perturbative fashion  @xcite . despite phenomenological agreement with the data
, the inclusion of finite cutoffs suggests that there might be some regulator dependence in those calculations .    in a series of recent papers  @xcite , two of us ( m.p.v . and e.r.a .
) have proposed not only to iterate but also to renormalize to all orders the nn chiral potential within a long distance expansion . by allowing the minimal number of counterterms to yield a finite result , long - distance regulator - independent correlations
are established . in practice ,
the potential must be computed within some power counting scheme .
while the potential is used within weinberg s power counting to lo , next - to - leading order ( nlo ) , and next - to - next - to - leading order ( n2lo ) , we only allow for those counterterms which yield a finite and unique scattering amplitude . in the @xmath0 np channel a single energy and momentum independent counterterm @xmath17
is considered which is determined by adjusting the physical scattering length . in the derivation of this result the mathematical requirements of completeness and self - adjointness for the renormalized quantum mechanical problem for a local chiral potential play a decisive role .
this surprising result is in contrast to the standard weinberg counting where an additional counterterm @xmath18 is included already at nlo .
this @xmath18 counterterm could , in fact , be determined by fitting the experimental value of the effective range ; the physics of @xmath18 is to provide a short distance contribution to the effective range in addition to the contribution from the known long distance chiral potential . within this context
it is remarkable that according to ref .
@xcite , where such a short distance contribution vanishes ( or equivalently @xmath19 when the cut - off is removed ) , rather accurate values are _ predicted _ yielding @xmath20 and @xmath21 _
after renormalization_. this latter value is less than @xmath22 larger than the experimentally accepted value and it suggests that most of the effective range is saturated by n2lo tpe contributions and calls for pinning down the remaining discrepancy .
this trend to convergence and agreement is also shared by higher order slope parameters in the effective range expansion without strong need of specific counterterms although there is still room for improvement .
motivated by this encouraging result , one goal of the present paper is to analyze the size of the next - to - next - to - next - to - leading order ( n3lo ) corrections to the results found in  @xcite .    the calculations in ref .
@xcite exploit explicitly the local character of the chiral potential by conducting the calculations in coordinate space which makes the analysis more transparent .
many results , in particular the conditions under which a renormalized limit exists , can be established a priori analytically .
moreover , the highly oscillatory character of wave functions at short distances is treated numerically using efficient adaptive - step differential equations techniques .
this situation contrasts with momentum space calculations where , with the exception of the pion - less theory , there is a paucity of analytical results , and one must mostly rely on numerical methods .
moreover , the existence of a renormalized limit is not obvious a priori and one may have to resort to some trial and error to search for counterterms .
finally , renormalization conditions are most naturally formulated at zero energy for which the momentum space treatment may be challenging , at times .
of course , besides these technical issues , there is no fundamental difference between proceeding in momentum or coordinate space , particularly _ after _
renormalization , provided the same renormalization conditions are specified , since disparate regulators stemming from either space are effectively removed . indeed , we will check agreement for the phase - shifts determined in different spaces whenever such a comparison becomes possible .
this equivalence is in itself a good motivation for renormalization .
however , at n3lo some unavoidable non - localities appear in the chiral long distance potential .
although they could be treated in configuration space , we adopt here a momentum space treatment .
this will also allow us to answer an intriguing question which was left open in the coordinate space analysis of previous works  @xcite , namely , the role played by the conventional momentum - polynomial representations of the short distance interaction used in most calculations  @xcite in the renormalization problem . more specifically , ref .
@xcite showed that taking @xmath19 was consistent , and a regularization scheme exists where a fixed @xmath18 was irrelevant , but could not discriminate whether @xmath23 was inconsistent as far as it was readjusted to the effective range parameter for any cut - off value .
, becomes irrelevant in the limit @xmath24 ( see below ) .
] in this regard , the present paper yields a definite answer making the surprising agreement of the effective range found in ref .
@xcite an inevitable mathematical consequence of renormalization .
the paper is organized as follows .
the problem is stated in sec .
[ sec : renorm - prob ] where a general overview of the renormalization problem is given both in momentum as well as coordinate space . in sec .
[ sec : mom - space ] , we particularize the momentum space formulation of the scattering problem with counterterms for the @xmath0 channel within a sharp three - momentum cutoff scheme .
likewise , in sec .
[ sec : coor - space ] we proceed similarly in the coordinate space formulation within a boundary condition regularization with a short distance cutoff @xmath25 . in sec .
[ sec : coor - mom ] , we discuss some features of both coordinate and momentum space formulations in the pion - less theory and try to connect the high momentum cutoff @xmath26 with the short distance radial cutoff @xmath25 .
this allows a one - to - one mapping of counterterms in both spaces which will prove useful later on in the pion - full theory .
the identification between the sharp momentum cutoff and the short distance radius found in the pion - less theory is discussed further in appendix  [ sec : nyquist ] in the presence of a long distance potential in the light of the nyquist theorem . in sec .
[ sec : n3lo ] , we come to the central discussion on n3lo corrections to the phase shifts when the scattering amplitude is renormalized with only one short distance counterterm . a wider perspective is achieved by further considering the role of explicit @xmath1-excitations in intermediate states and the subsequent one - counterterm renormalization of the scattering amplitude .
we also discuss the role of three pion exchange as well as how the results depend on the renormalization scheme used to compute the potential based on cut - off independent counterterms .
the mathematical justification for using just one counterterm is provided in sec .
[ sec : weinberg ] where the standard weinberg scheme is pursued both in momentum and coordinate space at nlo and n2lo and shown to potentially have some problems . finally , in sec .
[ sec : conc ] , we summarize our main points .
let us define the scope and goals of the present work .
the standard non - perturbative formulation of the renormalization problem starts with an effective lagrangian or hamiltonian  ( see e.g. @xcite and @xcite and references therein ) , from which a certain set of irreducible feynman diagrams ( usually up to a certain order ) is calculated .
these irreducible diagrams are defined to represent a potential @xmath27 .
the potential is then inserted into a scattering equation where it is iterated infinitely many times or , in other words , re - summed non - perturbatively . in the cm frame , where the np kinetic energy is given by @xmath28 , with @xmath29 ,
the scattering process is governed by the lippmann - schwinger equation @xmath30 with @xmath27 the potential operator and @xmath31 the resolvent of the free hamiltonian .
the outgoing boundary condition corresponds to @xmath32 .
using the normalization @xmath33 one has @xmath34 here @xmath26 means a generic regulator and represents the scale below which all physical effects are taken into account _
explicitly_. the degrees of freedom which are above @xmath26 are taken into account _ implicitly _ by including a suitable cutoff dependence in the potential . the precise equation governing this cutoff dependence
was described and studied in some detail in ref .
@xcite with particular emphasis on infrared fixed points .
we will analyze the cutoff dependence below focusing on the ultraviolet aspects of the interaction .
motivated by the low energy nature of the effective theory , the potential is usually separated into short and long distance components in an _ additive _ form @xmath35 where the long distance contribution is usually given by successive pion exchanges @xmath36 and the short distance component is characterized by a power series expansion in momentum @xmath37 where for simplicity we assume a spin singlet channel  . note that @xmath17 and @xmath18 contribute to s - waves while @xmath38 contributes to p - waves , and so on .
one should face the fact that , although the decomposition given by eq .
( [ eq : long - short ] ) is perturbatively motivated and seems quite natural , the additivity between short and long range forces is actually an _ assumption _ which has important consequences , as we will see .  )
does not foresee for instance terms of the form @xmath39 , i.e. terms which are not polynomial but influence the renormalization process .
of course , once additivity is relaxed there are many possible representations , in particular for the short distance components . ]    in principle , for a given regularization scheme characterized by a cutoff @xmath26 , the counterterms @xmath17 , @xmath38 , @xmath18 and so on are determined by fixing some observables .
one naturally expects that the number of renormalization conditions coincides with the number of counterterms in a way that all renormalization conditions are fully uncorrelated .
the statement of uv - renormalizability is that such a procedure becomes always possible when the cutoff @xmath26 is removed by taking the limit @xmath24
. this may not be the case as there may appear redundant contributions ( see the discussion below in sec .
[ sec : coor - mom ] ) meaning that one counterterm or counterterm combination can take any value .
another possible situation is just the opposite ; one may want to impose more renormalization conditions than possible .
in this case some counterterms or counterterm combinations are forbidden . rather than being intricate mathematical pastimes
, these features have already been investigated recently  @xcite for large cutoffs casting some doubt on the regulator independence of the original proposal  @xcite .    even if one admits generically eq .
( [ eq : long - short ] ) as well as eq .
( [ eq : long ] ) and eq .
( [ eq : short ] ) , it is not obvious how many terms should be considered and whether there is a clear way of defining a convergence criterium or identifying a convergence pattern .
it is fairly clear that , on physical grounds , one should consider an expansion of the potential that starts at long distance and decreases in range as the number of exchanged particles increases .
note , however , that , while ope is well - defined , the general form of the tpe potential is not uniquely determined . within this context ,
one of the main attractive features of the eft approach has been the definition of a power counting scheme which provides a hierarchy and an _ a priori _ correlation between long and short range physics . in the present paper
, we will assume weinberg power counting  @xcite ( see below for a more precise definition ) , where the long distance potential is determined in a dimensional power expansion .
the long distance component @xmath40 is obtained by particle exchanges and , in some simple cases , depends only on the momentum transfer   channel . ] .
in such a case , if we formally take a fourier transformation of the long distance potential @xmath41 and take the limit @xmath42 one has the standard schrdinger equation in coordinate space , @xmath43 where the coordinate space potential is @xmath44 +    \dots \ , \ , .
\label{eq : pot - coor}\end{aligned}\ ] ] the whole discussion , which has been carried on for years now , concerns the precise meaning of these delta and derivatives of delta interactions , particularly when a long distance potential is added to the short distance one .
a crucial finding of the present paper based on a direct analysis in momentum space is that _ non - perturbative renormalization imposes restrictions _ on the number of terms and form of the short distance potential which depend also on the particular long distance potential .
some of these restrictions were discussed in previous works  @xcite .
remarkably these new renormalizability restrictions were conjectured in coordinate space in ref .
@xcite for the @xmath0 channel based on self - adjointness and completeness of states and apply to the tpe potential ; a single @xmath45 counterterm is allowed while two counterterms , @xmath45 and @xmath46 , are forbidden  .      in the @xmath0
channel the scattering process is governed by the lippmann - schwinger equation @xmath47 where @xmath48 and @xmath49 are the scattering amplitude and the potential matrix elements , respectively , between off - shell momentum states @xmath50 and @xmath51 in that channel and the sharp three - momentum cutoff @xmath26 represents the scale below which all physical effects are taken into account _
explicitly_. from the on - shell scattering amplitude the phase shift can be readily obtained @xmath52 the short range character of the nuclear force implies that at low energies one has the effective range expansion ( ere ) @xmath53
\nonumber \\   & = & -\frac1{\alpha_0 }    + \frac12 r_0 p^2 + v_2 p^4 + v_3 p^6 + \dots \label{eq : ere}\end{aligned}\ ] ] where @xmath54 is the scattering length , @xmath55 the effective range and @xmath56 , @xmath57 etc .
are slope parameters .    in the @xmath0 channel ,
the potential is decomposed as the sum of short and long range pieces @xmath58 in the standard weinberg counting , the short distance contribution is written as follows @xmath59 where the counting is related to the order of the momentum which appears explicitly .
the long distance component of the potential is taken to be the sum of explicit pion exchanges @xmath60 where  @xcite @xmath61 using dimensional power counting .
ideally , one should determine the physically relevant long range regulator independent correlations , i.e. , long distance effects of similar range .
this would amount to consider _ all _ @xmath62 exchange effects on the same footing , since they yield a long distance suppression @xmath63 modulo power corrections . at present , the only way how these long distance potentials can be systematically computed is by dimensional power counting in perturbation theory , as represented schematically in eqs .
( [ eq : potpi ] ) and ( [ eq:1pi2pi3pi ] ) .    in the standard weinberg
counting one has @xmath64 note that this counting involves both unknown short - distance physics and chiral long - distance physics in a _ uncorrelated way_. note also that there is no first order contribution and that there is no third order contribution to the short distance potential .
regarding eq .
( [ eq : long - short-1s0 ] ) one should stress that the separation between long and short range contributions to the potential is not unique .
in fact , there is a polynomial ambiguity in the long range part which can freely be transferred to the short distance contribution .
however , the non - polynomial part is unambiguous as it is directly related to the left hand cut of the partial wave amplitude which for @xmath65 exchange is located at @xmath66 but presumably becomes incomplete for @xmath67 .   for any non - vanishing radius
( see e.g. the discussion in ref .  @xcite . ) ]      in coordinate space , the problem in the @xmath0 channel is formulated as follows  @xcite .
assuming a _ local _ long distance potential @xmath68 one has to solve the schrdinger equation @xmath69 where @xmath70 is the reduced potential ( in fact , the fourier transformation of @xmath71 ) and @xmath72 the reduced wave function for an s - wave state . here
@xmath25 is the short distance cutoff and the reduced wave function is subject to the boundary condition at @xmath73 and the standard long distance free particle behaviour @xmath74 where @xmath75 is the short distance phase - shift encoding the physics for @xmath76 . in the case of a vanishing long range potential @xmath77
the phase shift is given by @xmath78 .
on the other hand , if we take @xmath79 we get a standard problem with a hard core boundary condition , @xmath80 which for @xmath81 becomes the standard regular solution at the origin . at low energies both
the _ full _ phase - shift @xmath82 and the _ short distance _ phase - shift @xmath83 can be described by some low energy approximation , like e.g. , an effective range expansion , @xmath84 where @xmath85 is the short range scattering length , @xmath86 the short range effective range , and @xmath54 and @xmath55 the full ones   for a further discussion on this . ] .
if we also make an expansion at low energies of the reduced wave function @xmath87 we get the hierarchy of equations @xmath88 and @xmath89 and so on .
the standard way to proceed would be to integrate the equations for @xmath90 , @xmath91 , etc . from infinity downwards , with a known value of @xmath92 , using eq .
( [ eq : u0 ] ) to obtain @xmath93 and then one can use eq.([eq : sch_k ] ) together with eq .
( [ eq : bc_sch ] ) and eq .
( [ eq : ere_short ] ) to compute @xmath94 for any energy with a given truncated boundary condition .
this procedure provides by construction the low energy parameters we started with and takes into account that the long range potential determines the form of the wave function at long distances .
the only parameter in the procedure is the short distance radius @xmath25 , which is eventually removed by taking the limit @xmath95 .
one should mention at this point that the coordinate space is particularly suited for the case of local long distance potentials , but the renormalization with an arbitrary number of countertersm requires an _ energy dependent _ but real boundary condition at short distances which eventually violates self - adjointness . on the other hand ,
the momentum space formulation allows the discussion of nonlocal long distance potentials and the renormalization is done in terms of a _ momentum dependent _ short distance polynomial potential .
although this looks like a self - adjoint problem , we will see that in this formulation the counterterms may in fact become complex .
the renormalization of the pion - less theory , i.e. , a set of pure contact interactions , has been treated with great detail in the literature  @xcite although without much consideration on how this problem might be embedded into the wider and certainly more realistic situation where the finite range and short distance singular chiral nn potentials are present .
in fact , much of the understanding of non - perturbative renormalization within the modern nn context has been tailored after those and further studies based on the non - singular ope singlet @xmath0 potential  @xcite plus the standard perturbative experience . in previous  @xcite and in the present work ,
we pursue exactly the opposite goal : we will only consider renormalization procedures which can directly be implemented in the _ presence _ of long distance potentials since , after all , contact nn interactions are always assumed to approximate truly finite range interactions in the long wavelength limit .
thus , it is useful to review here those developments with an eye on the new ingredients which appear in the non - perturbative renormalization of singular pion exchange potentials as analyzed in later sections .
in addition , the deduced running of the counterterms in the contact theory in the infrared domain serves as a useful starting point when the long distance pion exchange potential is switched on .
finally , we will also discuss the size of finite cutoff corrections to the renormalized result depending both on the particular regularization as well as the corresponding representation of the short distance physics .
although the previously described momentum space framework has extensively been used in the past to describe successfully the data  @xcite with a _ finite cutoff _
@xmath96 it is worth emphasizing some puzzling features regarding the off - shell ambiguities of the short distance potential when finite range corrections , encoded in the @xmath18 , @xmath97 etc .
counterterms , are included .    in momentum space
, the pion - less theory corresponds to taking @xmath98 . in such a case
the lippmann - schwinger equation reduces to a simple algebraic equation  @xcite .
at very small values of the cutoff @xmath99 , the long range part of the potential may be neglected since they scale with powers of momentum and a simple contact theory of the form of eq .
( [ eq : vs(k , k ) ] ) may be used .
for instance , when @xmath100 , the lippmann - schwinger equation ( lse ) may be directly solved . using the basic integral   corresponds to making the replacement @xmath101 . in dimensional regularization ( minimal subtraction scheme )
the integral is just the unitarity piece , @xmath102 . ]
@xmath103 for a sharp momentum cutoff @xmath26 and @xmath104 , one obtains for the phase shifts @xmath105 at zero energy , @xmath106 and , thus , the running of @xmath17 is given by @xmath107 in this case , the phase shift is given by @xmath108 which corresponds to an ere with @xmath109 in the limit @xmath24 .
note that finite cutoff corrections scale as @xmath110 .
this indicates a relatively slow convergence towards the infinite cutoff limit and hence that finite cutoff effects are quantitatively important and might even become a parameter of the theory .
actually , one might determine @xmath26 by fixing the effective range from the first of eq .
( [ eq : pcotdmom0 ] ) , @xmath111 , yielding the accurate numerical value @xmath112 . in this case , this particular three - momentum regularization method becomes itself a model , since we have no control on the remainder . in any case
, it is straightforward to check that for any finite cutoff there is no off - shellness : @xmath113 .
the running given by eq .
( [ eq : c0-short ] ) must be used for _ any _ cutoff @xmath26 if we want to renormalize in the end . however , thinking of the more general case where finite range corrections are relevant such a running is only reliable for very small cutoffs @xmath114 .
if we consider also a @xmath115 coefficient in the potential , the corresponding lse can be solved with the ansatz @xmath116 which yields a set of three linear equations for @xmath117 , @xmath118 and @xmath119 .
after some algebraic manipulation , the final result for the phase shift can then be written in the form @xmath120 matching at low energies to the ere , eq .
( [ eq : ere ] ) , we get the running of @xmath17 and @xmath18 @xmath121 the first equation allows to eliminate uniquely @xmath17 in favour of @xmath54 and @xmath18 , but as we see there are two branches for the solutions . however , we choose the branch for which @xmath18 decouples in the infrared domain , i.e. fulfills @xmath122 for @xmath123 .
in fact , at small cutoffs , one gets for this branch @xmath124 the factor @xmath125 appearing in the small cutoff expansion for @xmath17 differs already from the coefficient in the case @xmath19 .
eliminating @xmath17 and @xmath18 in favour of @xmath54 and @xmath55 , the phase shift becomes @xmath126 note that the finite cutoff corrections are , after fixing @xmath55 , again @xmath127 .
so , fixing more low energy constants in the contact theory does not necessarily imply a stronger short distance insensitivity , as one might have naturally expected  .
in other words , the inclusion of a higher dimensional operator such as @xmath18 does not improve the ultraviolet limit , at least in the polynomial representation given by eq .
( [ eq : vs(k , k ) ] ) . in sec .
[ sec : n3lo ] we will show , however , that with just one counterterm @xmath17 the inclusion of pion exchange long distance contributions generates a much faster convergence towards the renormalized limit as anticipated in refs .  @xcite ( see also ref .
@xcite for a quantitative estimate ) . in sec .
[ sec : weinberg ] we will also show that when a @xmath128 counterterm is added this scaling behaviour is not only broken but also the phase shift fails to converge in the limit @xmath129 .
thus , we see that one can establish a one - to - one mapping between the counterterms @xmath17 , @xmath18 and the threshold parameters @xmath54 and @xmath55 . nevertheless , this is done at the expense of operator mixing , i.e. , both @xmath17 and @xmath18 are intertwined to determine both the scattering length and the effective range . in other words
the cutoff dependence of @xmath130 is different depending on the presence of @xmath18 . as we have seen
this is not a problem since for small cutoffs we expect the running of @xmath17 to be fully independent of @xmath18 and hence on @xmath55 .
however , unlike the one counterterm case , where @xmath19 , the solutions of eq .
( [ eq : c0-c2 ] ) may become complex when @xmath131 for the physical @xmath0 threshold parameters this happens already for @xmath132 ( the other two roots are complex ) . above this critical value
the potential violates self - adjointness . for @xmath133 one
has @xmath134 .
thus , the cutoff can only be fully removed with a self - adjoint short distance potential if @xmath135 .
this is consistent with the violation of the wigner causality condition reported in @xcite .
note that the violation of self - adjointness is very peculiar since once @xmath17 and @xmath18 have been eliminated the phase - shift ( [ eq : ps - lambda ] ) remains real   between off - shell momentum states , is violated , since the schwartz s reflection principle fails @xmath136^\dagger \neq t ( e-{\rm i } 0^+ ) $ ] .
this would also have far reaching consequences for the three body problem , since three body unitarity rests on two - body off - shell unitarity and self - adjointness of three body forces . ] .
one feature in the theory with two counterterms @xmath17 and @xmath18 is that the off - shell @xmath137matrix becomes on - shell only in the infinite cutoff limit , @xmath138 this is unlike the theory with one counterterm @xmath17 where there is no off - shellness at any cutoff .
thus , finite cutoff effects are also a measure of the off - shellness in this particular problem .
this will have important consequences in sec .
[ sec : weinberg ] when attempting to extend the theory with two counterterms in the presence of the long distance pion exchange potentials since the off - shellness of the short distance contribution of the potential becomes an issue in the limit @xmath24 .
the situation changes qualitatively when the fourth order corrections depending on two counterterms @xmath139 and @xmath140 are considered .
obviously , we can not fix both @xmath139 and @xmath140 _ simultaneously _ by fixing the slope parameter @xmath56 of the effective range expansion , eq .
( [ eq : ere ] ) . clearly , one expects some parameter redundancy between @xmath139 and @xmath140 or else an inconsistency would arise since a sixth order parameter in the effective range expansion @xmath57 should be fixed .
the situation worsens if higher orders in the momentum expansion are considered due to a rapid proliferation of counterterms while there is only one more threshold parameter for each additional order in the expansion .
this required parameter redundancy is actually a necessary condition for consistency which is manifestly fulfilled within dimensional regularization but not in the three - momentum cutoff method  .
moreover , it was realized some time ago  @xcite that the finite cutoff regularization and dimensional regularizarion in the minimal subtraction scheme yielded _ different _ renormalized amplitudes for a _ truncated _ potential . this non - uniqueness in the result due to
a different regularization happens when a non - vanishing @xmath18 counterterm is considered . in any case
, the dimensional regularization scheme has never been extended to include the long range part of the tpe potential which usually appear in the present nn context .
thus , for the momentum space cutoffs which have been implemented in practice the short distance representation is somewhat inconsistent at least for a finite value of the cutoff @xmath26 .
alternatively , one may choose an energy dependent representation of the short distance physics as @xmath141 in this case the correspondence between counterterms and threshold parameters @xmath54,@xmath55 , @xmath56 , etc .
is exactly one - to - one , and the parameter redundancy is manifest , since the on shell @xmath142-matrix depends only on the on - shell potential . actually , under dimensional regularization the representations of the potential eq .
( [ eq : vs(k , k ) ] ) and eq .
( [ eq : vs(p ) ] ) yield the same scattering amplitude .
although this on - shell equivalence is certainly desirable it is also unnatural , if the long distance potential is energy independent .
we will nevertheless analyze such a situation in the next subsection in coordinate space .
the previous discussion highlights the kind of undesirable but inherent off - shell ambiguities which arise when finite range corrections are included in the short distance potential   is considered .
it vanishes on the mass shell @xmath143 but nonetheless generates non trivial on shell scattering for the three - momentum cutoff . ] . in our view
these are unphysical ambiguities which have nothing to do with the unambiguous off - shell dependence of the long distance potential .
of course , one way to get rid of the ambiguities is to take the limit @xmath129 which corresponds to the case where a truly zero range theory is approached .
however , even for a finite cutoff there is a case where one is free from the ambiguities , namely when the short distance potential is _ both _ energy and momentum independent for s - wave scattering @xmath144 the key point is that we allow only this counterterm to be cutoff dependent and real , as required by self - adjointness .
of course , the discussion above for the contact theory suggests the benefits of using just one @xmath17 counterterm but does not exactly provide a proof that one _ must _ take further counterterms such as @xmath18 to zero .
the extension of this analysis to the case of singular chiral potentials in sec .
[ sec : weinberg ] will yield the definite conclusion that renormalizability is indeed equivalent to take @xmath19 .
the previous renormalization scheme is the momentum space version corresponding to the coordinate space renormalization adopted in a previous work by two of us ( mpv and era )  @xcite . actually , in the pure contact theory
, we can relate the renormalization constant with the momentum space wave function explicitly . at large values of the short distance cutoff @xmath25 ,
the zero energy wave function reads , @xmath145 thus , the following relation holds @xmath146 comparing with eq .
( [ eq : c0-short ] ) , we get @xmath147 where the momentum cutoff @xmath26 and the short distance cutoff @xmath25 are related by the equation @xmath148 which is nothing but an uncertainty principle relation between cutoffs  . ] .
note that for the standard regular solution @xmath149 one has a vanishing counterterm @xmath150 .
in contrast , @xmath151 for the irregular solution . in the case of the singular attractive potentials
the solution is regular but highly oscillatory and the @xmath130 takes all possible values for @xmath152 .
a more detailed discussion on these issues can be seen in refs .
of course , strictly speaking both eq .
( [ eq : c0-equiv ] ) and eq .
( [ eq : lambda - rc ] ) are based on a zero energy state , and in the finite energy case we will assume these relations having the limit @xmath152 or @xmath42 in mind .
let us now deal with finite energy scattering states .
since there is no potential , @xmath153 , for @xmath154 we have the free wave solution @xmath155 in the theory with one counterterm , we fix the scattering length @xmath54 by using the zero energy wave function and matching at @xmath73 so we get @xmath156 yielding @xmath157 and thus @xmath158 this is in qualitative agreement with the momentum space result when the cutoff is being removed , eq .
( [ eq : pcotdmom0 ] ) and , as we can see , the approach to the renormalized value is similar if the identification @xmath159 is made .
note further that since the boundary condition is energy independent the problem is self - adjoint and hence orthogonality between different energy states is guaranteed .
the theory with two counterterms where both @xmath54 and @xmath55 are fixed to their experimental values opens up a new possibility , already envisaged in ref .
@xcite , related to the non - uniqueness of the result both for a finite cutoff as well as for the renormalized phase - shift .
as pointed out above , this non - uniqueness was noted first in momentum space ref .
@xcite when using a finite three dimensional cutoff or dimensional regularizarion ( minimal subtraction ) .
remarkably , within the boundary condition regularization we will be able to identify both cases as different short distance representations .
actually , when fixing @xmath54 and @xmath55 we are led to @xmath160 where @xmath161 and @xmath162 are defined in sec .
[ sec : coor - space ] . note
that now self - adjointness is violated from the beginning due to the energy dependence of the boudary condition . within the second order approximation in the energy the neglected terms are @xmath163 , so _
any _ representation compatible to this order might in principle be considered as equally suitable .
the close similarity to a pad approximant suggests to compare the following three possibilities for illustration purposes @xmath164 \ , , \nonumber \\
d_p^{\rm iii } & = & \frac{u_0 ' ( r_c)^2}{u_0 ( r_c ) u_0'(r_c ) + p^2 \left [ u_2 ( r_c)u_0'(r_c ) - u_0 ( r_c ) u_2 ' ( r_c ) \right ] } \ , , \nonumber \\   \label{eq : d_p's}\end{aligned}\ ] ] and study what happens as the cutoff is removed , @xmath152 .
note that all three cases possess _ by construction _ the same scattering length @xmath54 and effective range @xmath55 and no potential for @xmath165 .
straightforward calculation yields @xmath166 as we see , all three representations provide the _ same _ threshold parameters , but _ do not _ yield identical renormalized amplitude for finite energy .
actually , cases i and ii coincide with the three - dimensional cutoff regularization method ( see sec .
[ sec : mom - pion - less ] ) , whereas case iii corresponds to dimensional regularization ( ms ) .
moreover , the finite cutoff corrections to the renormalized result are , generally , @xmath167 while the rational representation yields corrections @xmath168 .
these observations survive at higher orders when @xmath56 , @xmath57 , etc .
threshold parameters are further taken into account .
this indicates that not all short distance representations are equally `` soft '' in the uv - cutoff .
the generalization of these results to the case of singular tpe chiral potentials was studied in ref .
@xcite and will be also re - analyzed in sec .
[ sec : weinberg ] while discussing the consistency of the standard weinberg s power counting .    in any case
, when finite range corrections are considered within the boundary condition regularization , there are two possible renormalized solutions depending on the particular parameterization of short distance physics .
a nice feature of this regularization is that they can be identified with similar results found already in the momentum space analysis of ref .
@xcite when confronting three - momentum cut - off and dimensional regularization .
we note also here that no ambiguity arises when the boundary condition is assumed to be energy independent , in which case self - adjointness is guaranteed .
the study of the contact theory in sec .
[ sec : coor - mom ] provides suggestive arguments why it is highly desirable to carry out a regularization with a single counterterm in the @xmath0 channel by adjusting it to the physical scattering length for any cutoff value . in this section
we want to extend that study when the long distance chiral potential organized according to weinberg power counting enters the game and the cutoff is removed . by taking the cutoff to infinity
, we are actually assuming that all degrees of freedom not included in the present calculation become infinitely heavy . this way we expect to learn about missing physics in a _ model _ and _ regularization _ independent fashion .
the traditional strategy of adjusting an increasing number of counterterms may obscure the analysis .
in other words , by using this minimal number of counterterms , we try not to mock up what might be still missing in the long range description .    in this context
there is of course the question of convergence or cutoff insensitivity of the phase shift , when @xmath42 , provided we keep at any rate the scattering length @xmath54 to its physical value by suitably adjusting the unique counterterm @xmath169 . in coordinate space this is a rather trivial matter if the long distance potential is local  @xcite , as it happens in the lo , nlo , and n2lo weinberg counting .
the analysis in momentum space involves detailed large momenta behaviour of the lippmann - schwinger equation and one must resort to trial and error .
indeed , the n3lo case analyzed below includes nonlocalities and it _ turns out _ to provide convergent results .    for numerical calculations , we take the values for the @xmath170 and @xmath171 parameters appearing in the pion exchange potential @xmath172 used in ref .
@xcite , which do a good job for peripheral waves , where re - scattering effects are suppressed and where one is , thus , rather insensitive to cutoff effects
. we will only consider tpe contributions to the n3lo potential .      to determine the running of the counterterm we start from low cutoffs @xmath173 since the long range part of the potential
is suppressed and eq .
( [ eq : c0-short ] ) may be used .
actually , the analytical result is well reproduced by the numerical method used to solve the lippmann - schwinger equation in the pure contact theory with @xmath17 .
once this identification has been done , the value of the cutoff is increased steadily so that the scattering length is always fixed to the experimental value , @xmath174 .
this adiabatic switching on of the long range physics guarantees that we are always sitting on the correct branch which smoothly goes into the contact theory at low cutoffs   fitting @xmath54 .
they are physically unacceptable unless they behave as @xmath175 for @xmath176 since they do not evolve into the theory where the long range components are decoupled . ] .
the running of the counterterm @xmath177 is depicted for the lo , nlo , n2lo , and n3lo potentials in fig .
[ fig : c0 - 1s0.eps ] in the low cutoff region . as expected ,
the deviations from the simple result of the pure contact theory , eq .
( [ eq : c0-short ] ) , start at @xmath178 for lo because of the @xmath179 exchange potential . for higher cutoffs
the nlo , n2lo , and n3lo counterterm displays a cycle structure very similar to what has been observed in coordinate space  @xcite .
the convergence of the @xmath180 phase shift for fixed values of the lab energy , @xmath181 is displayed in fig .
[ fig:1s0-cutoff.eps ] .
of course , one observes a faster convergence for small energies . for the maximal value of @xmath182 , cutoff values @xmath183 are needed to change the phase shift by less than @xmath184 .    in agreement with the analytical estimates of ref .
@xcite , the convergence of the regulated phase - shifts towards their renormalized values follows a computable power like pattern , @xmath185 .
the more singular the potential at large momenta the faster the convergence .
thus , the expected increased insensitivity at short distances is indeed confirmed .
finally , the renormalized @xmath0 phase shift is presented in fig .
[ fig : ps1s0.eps ] for lo , nlo , n2lo , and n3lo . as a check of the present calculation in momentum space ,
let us mention that we reproduce the coordinate space renormalized phase shifts of @xcite at lo , nlo , and n2lo .
for instance , at the maximal cm momentum of @xmath186 , the maximal discrepancy between the coordinate space and momentum space phase shifts is less than half a degree when @xmath187 and @xmath188 , respectively .
the clear converging pattern can be observed all over the elastic scattering region , actually at @xmath189 , one has @xmath190 , @xmath191 , @xmath192 , and @xmath193 .
however , there is still a discrepancy with the nijmegen pwa result which at this energy is @xmath194 for np scattering .
to have an idea on the uncertainty of the calculated phase shift , we vary the scattering length @xmath195 , the value of , @xmath196 and @xmath197 .
actually the error in the latter is correlated through the goldberger - treiman relation , so we will quote both as an error in @xmath198 only .
in addition , for the chiral constants we take the central values used in @xcite which provided a good description for the peripheral waves and for the uncertainties we assume as an educated guess those from the @xmath199 study  @xcite . the only exception is @xmath200 , for which the error from @xmath199  @xcite is much smaller than the systematic discrepancy with the nn determination  @xcite .
so it is more realistic to take the systematic discrepancy as the error .
thus we take @xmath201 , @xmath202 , @xmath203 , @xmath204 , @xmath205,@xmath206 , @xmath207 .
the results for the particular variations are presented in table   [ tab : table1 ] and as we clearly see , the uncertainty in @xmath208 dominates the total error . at @xmath189 ,
one has @xmath209 , @xmath210 , @xmath211 , @xmath212 . as we see the statistical uncertainty stemming from the input parameters increases with the order . if we take the difference @xmath213 as an estimate of the systematic error , adding in quadrature statistical and systematic errors we have , @xmath214 , still a smaller quantity than the discrepancy to the nijmegen phase shift .
this may suggest that _
after renormalization _ there are some physical effects missing even at the n3lo level beyond pure tpe .    [ cols="^,^,^,^,^,^,^,^,^,^,^,^",options="header " , ]      in the previous section we have seen that all tpe effects included to n3lo display a convergent pattern after renormalization , but there is still some missing physics .
note that while at lo and nlo the only parameters are @xmath198 , @xmath215 , @xmath216 and @xmath217 , at n2lo and n3lo there appear new low energy constants ( the @xmath170 and @xmath171 , respectively ) which can be related to @xmath199 scattering and encode short range physics not considered explicitly .
the @xmath1 resonance is an outstanding feature of @xmath218 scattering and explains a great deal of the low energy constants .
thus , it is interesting to analyze the role of explicit @xmath1 excitations as intermediate states in the @xmath219 potential .
the importance of explicit @xmath1 degrees of freedom has been emphasized on power counting grounds in several previous works with finite cutoffs where the @xmath220 splitting is regarded as a small parameter @xmath221  @xcite .
the crucial role played in the renormalization problem has been stressed in ref .  @xcite . in this section
, we analyze the nn potential with the nlo terms of @xcite together with the @xmath222 and @xmath223 in the box diagrams as computed in ref .
one advantage of such an approach is that , as compared to the standard @xmath1-less theory , there only appears the @xmath220 splitting as a parameter .
the renormalization of the @xmath0 channel proceeds along the lines discussed in the @xmath1-less theory .
the results for lo , nlo , nlo+@xmath222 , and nlo+@xmath224 renormalized phase shifts in the @xmath0 channel are plotted and compared to the nijmegen partial wave analysis  @xcite as a function of the cm momentum in fig .
[ fig : ps1s0-del.eps ] .
the most striking result is the very strong resemblance between the n2lo @xmath1-less vs nlo+@xmath222 and the n3lo @xmath1-less vs nlo+@xmath224 .
these results sustain the treatment of the @xmath225 splitting as a small parameter @xmath226 corresponding to the @xmath1-counting nlo@xmath227 = nlo + @xmath222 + @xmath224 .
the fact that all the contributions fall off at large distances as @xmath228 suggests that despite the ability to mimic higher order corrections in the weinberg counting for the long distance potential in the @xmath0 channel there is still missing shorter range physics beyond tpe .      at n3lo
, three - pion exchange occurs for the first time .
these contributions have been calculated by kaiser  @xcite and found to be small , which is why present n3lo nn potentials omit these contributions when renormalization is not implemented .
however , it should be noted that , for small distances , the @xmath229 diagrams are proportional to @xmath230 and , thus , will ultimately dominate at short distances .
a rough estimate of the results published in refs .
@xcite suggests that the sum of all @xmath229 graphs is attractive .
thus , we infer from there that the scattering length will still be a free input parameter . on the other hand
the @xmath231-exchange contribution falls off as @xmath232 at long distances so it becomes active at rather short distances , and so the effect is expected to be small because of the short distance suppression ( modulo oscillations ) of the wave function @xmath233 typical of potentials with a short distance power like singularity  @xcite .
this agrees with the rule that the more singular the potential the more convergent is the calculation , as we have extensively discussed above .
therefore , in a complete and renormalized n3lo calculation of the @xmath0 phase shifts , we expect the @xmath229 effects not to be large , although the predictions may be closer to the empirical values than in fig .
[ fig : ps1s0.eps ]
. an accurate investigation of the impact of @xmath229 exchange at n3lo on nn phase shifts represents an interesting and challenging project for the future .
as we have mentioned , the very definition of the potential is ambiguous as it requires fixing the polynomial terms in the momentum . in the calculations presented above
, we have taken the renormalization scheme for the potential where the fixed and cut - off independent choice @xmath19 for the potential is made .
we have also analyzed the situation when a different renormalization scheme is taken , namely a non - vanishing arbitrary @xmath18 coefficient which does not run with @xmath26 .
however , we allow @xmath177 to run in a way that the scattering length @xmath92 is reproduced .
this generates a different renormalization trajectory for @xmath234 as compared to the case @xmath19 .
we find by actual calculations that this fixed @xmath18 coefficient is irrelevant , i.e. the renormalized phase shift does not depend on this fixed value in the limit @xmath129 . roughly speaking ,
the reason is that while the polynomial combination @xmath235 is large at large momenta , the pion exchange part behaves as @xmath236 with an additional logarithm and thus dominates for fixed @xmath18 .
this irrelevance of @xmath18 was highlighted in the coordinate space analysis of ref .
@xcite where the regularization based on a radial cut - off @xmath25 would provide a compact support not sensing the details of the distributional contributions in eq .
( [ eq : pot - coor ] ) , regardless on how many derivatives of delta s are included .
a quite different situation arises when @xmath18 depends on the cut - off in a way that the effective range is fixed as we discuss in the next section  [ sec : weinberg ] .
there , it will be shown that if @xmath237 is relevant then the phase shift is not convergent .
in previous sections , we have seen that the lo , nlo , n2lo , and n3lo chiral potentials can be renormalized when one counterterm @xmath169 is determined for any value of the cutoff @xmath26 by fixing the scattering length to its experimental value , @xmath238 and the limit @xmath42 is subsequently taken .
this agrees with the observation of refs .
@xcite that attractive singular potentials can be renormalized with a single counterterm .
however , weinberg counting requires further counterterms in the short distance potential , see eq .
( [ eq : vs(k , k ) ] ) . for instance , both at nlo and n2lo a @xmath128 counterterm should be included .    in this section
, we discuss whether the @xmath0 scattering amplitude is renormalizable , i.e. whether the scattering amplitude has a well defined limit , when both the @xmath177 and @xmath239 counterterms are included and fixed by fitting the scattering length @xmath240 as well as the effective range @xmath241 for any value of the cutoff @xmath26 and the limit @xmath24 is pursued .
actually , when the nlo ( n2lo ) long distance potential , eq .
( [ eq:1pi2pi3pi ] ) is used , this way of proceeding corresponds to renormalizing the nlo ( n2lo ) approximation of the @xmath0 channel in the standard weinberg counting .
the method we follow in practice is a straightforward extension of the case with only just one counterterm @xmath17 . in the case
when @xmath242 , we get a finite ( renormalized ) effective range @xmath243 , which is close but still differs significantly from the experimental value . at n2lo with @xmath19
one gets @xmath21 .
this last value is so close to the experimental one that one would not expect big changes when a @xmath18 is added to exactly fit the experimental @xmath55 .
thus , it makes sense to investigate what would happen if one uses a non - vanishing @xmath18 to account for the missing , and in principle tiny , contribution to the effective range .
the surprissing result , to be discussed below in detail , is that trying to fit the discrepancy in the effective range with a @xmath18 counterterm is incompatible with renormalizability .    in the contact theory we found that for @xmath244 the counterterms @xmath17 and @xmath18 diverge before becoming complex numbers .
thus , the corresponding hamiltonian becomes non - self - adjoint , in harmony with the wigner causality bound violations unveiled in ref .
when the nlo contribution to the potential is included , the critical value of the cutoff is slightly shifted towards the higher values @xmath245 ; above those values the short distance contribution to the potential becomes non - self - adjoint and causality bounds are violated .
nevertheless , and similarly to the pion - less theory , the phase shifts remain real beyond this critical cutoff . in fig .
[ fig:1s0-ps - c0-c2-mom.eps ] the value of the phase shift for both nlo and n2lo approximations is depicted for several cutoff values . as we see , at nlo one observes large variations when the cutoff is changed from @xmath246 to @xmath247 . in fig .
[ fig:1s0-ps - c0-c2-mom.300.eps ] , we plot the value at a fixed lab energies as a function of the cutoff for both nlo and n2lo . as we see
the phase shift does not seem to converge to any particular value when @xmath26 is increased ; in fact large variations can be clearly seen at moderate cutoff values , @xmath248 for @xmath249 and @xmath250 for @xmath251 . on the other hand ,
let us note that a plateau region not always appears , and when it does the residual contribution from pion - exchange effects is less important and the counterterms may dominate the calculation .
in other words , if the cutoff was too small , one would be driven back to an effective range expansion with no visible contribution from chiral potentials whatsoever .
thus , any finite cutoff calculation within a higher cutoff regime turns out to be strongly cutoff dependent , or else the cutoff must be fine tuned to intermediate energy data , hence becoming an esssential , and not an auxiliary , parameter of the theory .
note that all the problems are triggered by _ insisting _ on fitting the effective range parameter to the experimental value by introducing the @xmath18 counterterm required by weinberg power counting on the short distance interaction .
in contrast , if @xmath19 , as advocated in refs .
@xcite and sec .
[ sec : n3lo ] above , not only is the phase shift convergent in the limit @xmath42 ( in practice @xmath252 ) but also most of the effective range is saturated by the chiral potential   and under the assumption of van ver waals dominance one gets at n2lo the analytical result ( see ref .
@xcite for details ) , @xmath253^\frac14\ ] ] which yields @xmath254 , a surprissingly good approximation .
further studies exploiting the chiral van der waals correlations can be found in ref .
@xcite . ] .
although we have checked that these results are numerically robust within the lse by increasing the number of grid points , there is always a reasonable doubt , since these are demanding calculations .
thus , it would be nice to understand the lack of convergence from a different perspective , as we do next .
let us compare the previous findings in momentum space with related investigations in coordinate space  @xcite .
a thorough study in coordinate space has revealed that the very existence of the @xmath42 limit may actually depend on the specific representation of the short distance physics .
this analytical result has been verified by numerical calculations and will become extremely helpful in analyzing the momentum space calculations presented above .
thus , for the sake of completeness it is worth reviewing the emerging pattern from ref .
@xcite .    as we have discussed in sec .
[ sec : coor - space ] , within the boundary condition regularization the _ unknown _ short distance interaction is represented by the logarithmic derivative of the wave function at the short distance cutoff radius @xmath25 , @xmath255 . in an energy expansion of the wave function at short distances ,
one has @xmath256 and its logarithmic derivative for which a continuity condition is required .
this introduces an energy dependence which eludes the wigner causality bound discussed in refs .
@xcite , since self - adjointness is broken from the start . as we discussed previously in the contact theory described in sec .
[ sec : coor - pion - less ] , at second order in the energy the neglected terms are @xmath163 and _ any _ of the three representations displayed by eq .
( [ eq : d_p s ] ) might be equally acceptable .
actually , we found that renormalized amplitudes fall into two classes and that not all of them approach the renormalized limit in the same way , see eq .
( [ eq : coor - pion - renor ] ) . here
, we will extend that study by inquiring what happens when the same three representations displayed by eq .
( [ eq : d_p s ] ) are used and the long distance potential @xmath257 is taken to be the nlo and n2lo for @xmath154 as the cutoff is removed , @xmath152 .
note that all three cases possess _ by construction _ the same scattering length @xmath54 and effective range @xmath55 and the same long distance potential for @xmath154 .
thus , any difference is clearly attributable to the different short distance representation .
the results are displayed in fig .
[ fig:1s0-ps - c0-c2-coor.eps ] for fixed short distance cutoff values @xmath25 as a function of the cm momentum @xmath258 . for finite values of @xmath25
we see a difference which can naturally be explained by the different off - shell behaviour of the short distance physics . as we see
, the difference persists as the cutoff is being removed and in fact is magnified .
as was already pointed out in a previous work  @xcite , only the case @xmath259 representing a rational function turns out to yield a unique and well defined finite value for the phase shift  ) . ] .
this fact is clearly seen from inspection of fig .
[ fig:1s0-ps - c0-c2-coor.300.eps ] where we plot the phase shifts for a fixed value of the cm momentum @xmath260 as a function of the equivalent sharp momentum cutoff @xmath261 derived for the pion - less theory in sec .
[ sec : coor - mom ] .
the striking similarity between the momentum space calculation presented in fig .
[ fig:1s0-ps - c0-c2-mom.300.eps ] and the coordinate space calculation displayed in fig .
[ fig:1s0-ps - c0-c2-coor.300.eps ] for _ finite cutoffs _ is noteworthy although not completely surprising in the light of the analysis of appendix  [ sec : nyquist ] .
there , the finite cutoff equivalence , @xmath262 , deduced in the contact theory is shown to hold also in the presence of a local potential .    of course
, all these features depend on the singular character of the interaction at short distances and do not depend on the specific form of the potential , so we expect them to hold also when @xmath1 degrees of freedom are explicitly taken into account .
although we can not prove it analytically in momentum space , the remarkable similarity of the coordinate space analysis with the present momentum space calculations strongly suggests that the standard polynomial representation of the short distance potential is _ incompatible _ with renormalization .
of course , this does not preclude the possible existence of a suitable potential representation of the short distance interaction , most likely with complex counterterms , such that both low energy parameters can be fixed and the cutoff can at the same time be removed .
we leave such an interesting investigation for the future .
in the present paper , we have analyzed the renormalization of the singlet np phase shift in the @xmath0 channel incorporating one- and two - pion exchange effects . for the long distance potential , the standard weinberg scheme based on dimensional power counting
is adopted to n3lo .
however , the short distance physics is parameterized in terms of one unique energy and momentum independent counterterm whose cutoff dependence is determined by adjusting the scattering length to its physical value .
the present analysis is carried out in momentum space in a somewhat complementary manner as previously done in coordinate space  @xcite .
actually , we have analyzed and reproduced those results directly in the more popular momentum space after the cutoff has been effectively removed . in order to stress the equivalence of both approaches , we have displayed many results in parallel .
this is not just a stylistic matter of presentation ; besides numerical simplicity , much understanding of the renormalization problem has been achieved by a direct analysis in coordinate space at least for local potentials .
the present work provides a further example in this respect .
the momentum space framework allows a direct extension to include n3lo contributions which include non - local pieces in the long distance chiral potential .
the main outcome of such a calculation is that the n3lo chiral potential induces rather small corrections as compared to the n2lo results in all of the elastic scattering region .
as anticipated previously  @xcite , this happens to be so even when the n3lo potential is more singular than the n2lo one at short distances .
moreover , the analytical scaling behaviour for large cutoffs @xmath26 of the scattering amplitude predicted in a previous coordinate space analysis  @xcite is confirmed qualitatively by the present momentum space calculations ; increasing the order in the expansion of the long distance potential suppresses more the finite cutoff dependence . in summary ,
the scheme is convergent , but differs from the expected phase shifts obtained from partial wave analyses .
an error analysis of the results shows that although there are some uncertainties induced by the input parameters , the corresponding error bands are not large enough to be compatible with the partial wave analysis of the nijmegen group .    the n3lo approximation to the long distance chiral potential
is computed within a heavy baryon expansion and assuming only explicit nucleon degrees of freedom in the nn potential .
this implies in particular that the @xmath220 splitting is considered a non - small parameter .
however , this number is about twice the pion mass , so it is not clear whether such an assumption is fully justified .
therefore , we have analyzed the nn scattering problem when explicit @xmath1 intermediate state excitations are included in the potential .
the net result is that , after renormalization , the @xmath222 and @xmath224 box diagrams contributions mimic extremely accurately the @xmath1-less n2lo and n3lo potentials , respectively , and thus fail to describe the higher energy region of the @xmath0 phase shift .
given the fact that these contributions also fall off at large distances as @xmath263 , this result reinforces the conclusion that there is some shorter range missing physics beyond that provided by tpe .
further , it is interesting to compare the present renormalized results with those obtained for a two pion exchange potential computed within a relativistic baryon framework which sums up all heavy nucleon components  @xcite thus contains a full determination of tpe contributions assuming all other degrees of freedom ( including explicit @xmath1 s ) are infinitely heavy .
there , all np partial waves with @xmath264 are analyzed and for the particular case of the @xmath0 channel the results are rather similar to those found here .
this , again , suppports the present conclusion that , presumably , the components of the nn potential with a shorter range than tpe might finally provide the missing repulsion needed to reduce the @xmath22 overshooting of the @xmath0 effective range as well as the too low value of the phase shift in the region with cm momenta @xmath265 .
we have also investigated the convergence of the standard weinberg counting for both the short distance as well as the long distance potential . at nlo and
n2lo this corresponds to including a polynomial momentum dependence in the short distance interaction by means of two counterterms which may then be fixed by adjusting the scattering length and the effective range to their experimental values .
we find that , with these two renormalization conditions , the counterterms turn out to be complex for not too large cutoffs ( around @xmath266 ) , signaling the breakdown of self - adjointness of the potential , and more generally suggesting a violation of wigner s causality condition . nonetheless , despite phase shifts remaining real , a unique renormalized limit does not exist .
this lack of convergence agrees with similar findings in coordinate space anticipated in ref .
the physical understanding of the situation is as follows . on the one hand the long range physics
is fixed , thus , the off - shell propagation within the long range region is unambiguous . on the other hand
, the short range potential must be adjusted to provide the threshold parameters such as the scattering length , the effective range etc .
if the short range potential is just determined from fixing the scattering length only , this is at zero energy and becomes zero range when the cutoff is removed , so there is no off - shell ambiguity .
in contrast , fixing further the effective range requires non - zero energy and there are in fact infinitely many ways of parameterizing this , with different off - shell behaviour even when the cutoff is removed .
in other words , fixing a finite range and removing the cutoff depends on details of the method .
these conflicts between off - shellness and finiteness are not new in field theory .
green s functions which are renormalized on - shell do not necessarily provide finite off - shell amplitudes .
the standard weinberg parameterization as a polynomial in momenta is one choice which may not turn out to be consistent with renormalizability . in this regard ,
it is interesting to mention that coordinate space detailed studies based on renormalization group properties  @xcite suggest the existence of a suitable short distance representation yielding convergent results in theories with more than one counterterm .
the generalization of those results to momentum space would require an in - depth study of the renormalization group in the presence of eventually nonlocal but singular potentials .
we note that the energy dependence of the coordinate space solution violates self - adjointness explicitly and we expect that most likely an energy independent non - polynomial momentum space solution would spoil self - adjointness as well .
it is important to note that with a finite cutoff and four counterterms the @xmath0 phase shift has been successfully described within the standard weinberg counting  @xcite . in this regard ,
it is natural to question the usefulness of taking the infinite cutoff limit and to carry out a renormalization process . from a physical point of view
the limit @xmath42 corresponds to consider other degrees of freedom than those not considered explicitly to be infinitely heavy .
thus , the aim of the renormalization program is far more stringent than previous finite cutoff calculations .
within such a framework any failure can be unambiguously attributed to missing physical information on the long distance potential , and the renormalization process highlights this in a rather vivid manner as one clearly sees when going from ope to tpe potentials . from a mathematical perspective the renormalizability requirement imposes tight constraints on the admissible forms of the short distance physics for some given long distance interactions .
these are powerful conditions which have traditionally been the great strength of the renormalization ideas to avoid unnecessary proliferation of interactions .
since the power counting of long distance potentials is not uniquely determined yet ( one may e.g. include or not explicit @xmath1 s ) , it would be very helpful to see what kind of short distance constraints on those potentials might be imposed as well on the basis of renormalizability or other principles .
the problems with the weinberg counting for the short distance polynomial form in the momenta of the interaction with renormalization found in this paper for nlo and n2lo is to be added to the other ones noted in previous works @xcite . on the mathematical side ,
it is noticeable that the coordinate space analysis of ref .
@xcite conjectured this result by exploiting the compelling requirement of completeness and self - adjointness for the renormalized quantum mechanical problem for a local chiral potential which at first sight may seem completely germane concepts to the eft machinery . on the phenomenological side
, it should be noted that this inconsistency result does not explain why the renormalized phase shift with just one counterterm comes out reasonably close to the accepted ones , but certainly makes this unique and finite prediction more credible and inevitable from a theoretical perspective if self - adjointness is maintained , and provides further confidence on the virtues of a renormalization principle within the chiral approach to the nn interaction .
one of us ( e.r.a . ) thanks diego pablo ruiz padillo ( departamento de fsica aplicada , universidad de granada ) for patient instruction on the virtues of nyquist theorem some years ago .
we also thank evgeny epelpaum for a critical reading of the manuscript and daniel phillips for remarks on the paper .
the work of d. r. e. has been partially funded by the ministerio
de ciencia y tecnologa under contract no .
fpa2004 - 05616 , the junta de castilla y len under contract no .
the work of e.r.a . is supported in part by funds provided by the spanish dgi and feder funds with grant no .
fis2005 - 00810 , junta de andaluca grants no .
fqm225 - 05 , eu integrated infrastructure initiative hadron physics project contract no .
rii3-ct-2004 - 506078 .
m. p. v. has been funded by the deutsche forschungsgemeinschaft ( sfb / tr 16 ) , helmholtz association ( contract number vh - ng-222 ) .
r. m. has been supported in part by the u. s. national science foundation under grant no .
in this appendix we discuss further the relation between momentum and coordinate space for _ finite cutoffs _ and local potentials . we will show that solving the lippmann - schwinger equation for a long range potential with a sharp cutoff @xmath26 is equivalent to the solution of the schrdinger equation for the same potential in a discretized grid with @xmath267 .
this equivalence reminds the previous identification between sharp momentum cutoff @xmath26 and the short distance cutoff @xmath25 found in the contact theory in sec .
[ sec : coor - mom ] , @xmath268 .
we can choose the initial grid point at @xmath269 , thus recovering the previously mentioned equivalence . actually , by invoking nyquist theorem on the cut - off ls equation , we deduce a discretized version of the variable phase equation  @xcite which allows to discuss both the renormalization as well as the decimation problem of the nn - force based on chiral potentials .
let us consider as a starting point the lippmann - schwinger equation for @xmath270 wave scattering , written in the form @xmath271 where @xmath272 are the matrix elements of the @xmath137 matrix between initial and final cm momentum states @xmath50 and @xmath51 respectively and the corresponding potential matrix element is given by @xmath273 for a local potential .
since we are integrating the intermediate momentum up to infinity , we are implicitly assuming that @xmath49 is a regular potential .
the singular potential case will be discussed later .
now , if we cutoff the potential in momentum space  ) . ] , we get the regularized potential @xmath274 where @xmath275 is the heaviside step function . in
coordinate space the cutoff potential becomes @xmath276 which is obviously nonlocal and becomes local only when @xmath129 ( second line ) .
the cutoff lse , where @xmath277 , can be solved by standard means . in the spirit of an eft , based on the idea that low energy dynamics does not depend on the details at short distances , it may perhaps be appropriate to proceed a bit differently .
actually , if high momentum states are cutoff from the theory the same idea should apply to small _ resolution _ scales , i.e. regardless whether they are short or long .
that means that using _ too much _ information on the potential @xmath278 even if it is exactly known point - wise may be illusory at wave - lengths longer than a given resolution , @xmath279 ,   is supplemented by highly oscillatory ripples at short resolution scales , @xmath280 the phase shifts should be insensitive to them at @xmath281 .
] so that one can sample @xmath278 with _ some _ resolution @xmath282 .
this obviously reduces the number of coordinate mesh points in the integration .
in the present context , nyquist theorem  @xcite is remarkably useful because it provides an optimal way of sampling signals which have a bandwidth in fourier space , i.e. , functions for which @xmath283 then , for the original function we have @xmath284 if we define the sampling function @xmath285 of the function @xmath286 at the equidistant points @xmath287 , with the optimal @xmath288 separation , @xmath289 and compute its fourier transform we get @xmath290 note that for the optimal sampling @xmath291   is bandwidth limited to the @xmath292 $ ] interval , it can be expressed as a fourier sum @xmath293 where @xmath294 .
thus it is trivial to see that @xmath295 when the sampling is done with @xmath296 . ] . inverting the fourier
transform we get @xmath297}{{\pi } ( x - x_n ) } \delta x \ , .\end{aligned}\ ] ] thus , if we sample the function according to eq .
( [ eq : sampl ] ) , the following identity holds at the sampling points @xmath298 hence , there is no loss of information on the sampling points if the sampling is done equidistantly with the optimal nyquist frequency , @xmath299 . in particular
, it does not really make sense to sample the function with smaller @xmath300 . in the next section
we apply this sampling principle to the potential .
nyquist theorem also applies for the special case of the lippmann - schwinger equation in s - wave scattering , although the derivation is different in some details to the one presented in the previous section . for clarity , we present here these details
.    we will consider first the general case of a non - local potential , for which the schrdinger equation in s - wave reads @xmath301 we can sample this non - local potential as @xmath302 for the momentum space representation of the sampled potential we get @xmath303 where @xmath304 . as a consequence of the cut - off @xmath26
, the potential @xmath305 can be expressed as a sum of spherical bessel functions @xmath306 where @xmath307 . then by taking @xmath267 , the sampled potential recovers the original one , i.e. @xmath308 , for @xmath309 . by fourier - transforming @xmath310 back to coordinate space
, it can be checked that it reproduces the original sampling points , i.e. @xmath311 for the nyquist sampling frequency .    in the case of a local potential , the one which interests us most , we sample the following way @xmath312 after double fourier - transforming , we get @xmath313 which makes the potential local for the grid points .
it is in fact a finite cut - off version of @xmath314 once we notice that @xmath315 @xmath316 .      for a lippmann schwinger equation with a finite cut - off @xmath26
all the matrix elements become a linear combination of separable terms . thus the ls equation becomes a linear matrix equation which can be solved by standard techniques by writing @xmath317 and defining propagator matrix elements @xmath318 therefore
we get @xmath319 this equation can be reduced to a finite - dimensional @xmath320 linear algebra problem by cutting the sums to @xmath321 .
the effect of this simplification can be seen by sampling @xmath322 in coordinate space ( see eqs .
( [ eq : sampl ] ) and ( [ eq : sampl - local ] ) ) @xmath323 from which we can check that cutting the sum is equivalent to introducing the ( harmless ) infrared cut - off @xmath324 .
note that for local potentials , the sampled potential matrix elements at the grid points are diagonal @xmath325 thus eq .  ( [ eq : ls_grid ] ) becomes @xmath326 which looks like a multiple scattering equation , with on - shell propagation between delta - shell scatterers . after matrix inversion , the on - shell solution is then given by @xmath327      in the previous section we have seen that for the lippmann - schwinger equation with a cut - off @xmath26 , we can either use the original momentum space potential @xmath305 and solve by standard means , or expand this potential , i.e. use the sampled potential @xmath310 , and solve as a linear algebra problem ( as they are both the same potential for momenta below the cut - off ) .
alternatively we can directly solve the schrdinger equation for the sampling potential @xmath328 .
this procedure will give an excellent approximation to the solution of the ls equation ( although not the exact solution , as explained at the end of this section ) but at a much smaller computational cost .
for this purpose we make the replacement @xmath329 which is a superposition of equally spaced delta - shell potentials .
it should be noted that the point @xmath330 does not contribute , as it lies at the integration boundary of the schrdinger equation .
then , we have @xmath331 thus , for a potential of size @xmath332 where we do not want to describe energies higher than @xmath26 we should do with an infrared cut - off larger than the potential s length scale , @xmath333 , or equivalently @xmath334
. we can solve the schrdinger equation piecewise , @xmath335 where @xmath336 is the amplitude and @xmath337 can be understood as the accumulated phase - shift due to adding a new delta shell at @xmath338 chosen to be located at the midpoint @xmath339 ( for reasons to become clear soon ) .
note that with this choice the lowest possible location of the phase - shift , @xmath340 corresponds to the point @xmath341 which we may identify with a short distance ( ultraviolet ) cutoff .
this identification between the momentum - space cutoff agrees with the one obtained for the pure short range theory by solving the lse without any discretization ( see sec .
[ sec : coor - mom ] ) .
it should be noted that the previous method for solving the schrdinger equation does not really generate the exact solution of the ls equation with a cut - off , but only a close approximation . although @xmath328 is the best sampling function for @xmath342 , it is not the same quantum mechanical potential .
the solution to the ls equation with @xmath342 can be exactly reproduced by solving the non - local schrdinger equation , eq .
( [ eq : schroedinger_nonlocal ] ) , with the potential @xmath343 which comes from inverse fourier transforming @xmath310 .
the schrdinger equation should be solved with trivial initial conditions , @xmath344 , as all the physically relevant information is included in @xmath343 .    on the contrary , if we solve @xmath328 , a sum of delta shells , we must include a non trivial initial condition   sampling point .
] , even in the absence of any short range physics . if we solve the discretized schrdinger equation with a trivial initial boundary condition , @xmath345 , then a certain error will be included in the final solution .
in the worst case we can expect an error of order @xmath279 . by means of the variable phase equation
@xcite we can perform a better assessment of the error  , the variable phase equation will yield @xmath346 as a curiosity , we can see that @xmath340 will scale as an inverse power of @xmath279 for a singular potential , thus signaling the need of a counterterm .
] , yielding , for example , an @xmath347 error for a yukawa potential or @xmath348 for a square well .    from the previous discussion
, it is apparent how to include a counterterm in the computation .
it enters throught the initial sampling point of @xmath349 , which maps onto a constant term in @xmath310   counterterm when projected to the s - wave takes the form @xmath350 thus , if we write the sampled potential @xmath351 , and take into account eq .
( [ eq : v_ij_nyquist ] ) , we see that @xmath17 maps only onto @xmath352 , giving a constant contribution in momentum space , as expected . ] .
thus it is ignored when solving the schrdinger equation for @xmath349 , and it enters through the initial condition @xmath340 .
similar remarks can be made for singular potentials , in which a counterterm must be included in order to obtain a stable result for @xmath353 .
the previous discussion can be elaborated further to reach interesting results . matching the wave functions at the points where the delta shells are located , @xmath354
, we simply get @xmath355 where @xmath356 is the reduced potential .
this is a recursion relation for the phase - shift at the interval midpoint @xmath357 . unlike the matrix equation which has traditional storage limitations for large number of grid points @xmath358
, this equation does not posses this shortcoming allowing for rather large @xmath358 values .
similar equations were deduced many years ago  @xcite as a practical tool to attack the inverse scattering problem and to determine the nn potential on the grid points .
defining the discretized effective range function , @xmath359 we get @xmath360 which can be rewritten as a continuous fraction .
note that for the cutoff theory both eqs .
( [ eq : vf_dis ] ) and ( [ eq : m_dis ] ) are almost exact .
the only approximation comes from the finiteness of the cutoff @xmath26 .
an important property which will be used later on is the reflection property , namely the symmetry under the replacement @xmath361 which means that on the grid running the relation upwards or downwards are inverse operations of each other .
obviously this property may fail in practice only due to accumulation of computer round - off errors over large evolution distances .
this is also the reason why we choose the midpoint @xmath362 for the accumulated phase shift : in any other case we would lose the reflection property .
the regular solution at the origin reads @xmath363 if we take the limit @xmath42 we can define @xmath364 , to get @xmath365 which is the variable phase equation  @xcite up to finite grid corrections and can be interpreted as the change in the accumulated phase when a truncated potential of the parametric form @xmath366 is steadily switched on as a function of the variable @xmath367 .
( [ eq : vf_dis ] ) is thus a discretized variable phase equation , corresponding to a discretized potential sampled at the optimal nyquist frequency .
this equation and its generalization to coupled channels has extensively been used to treat the renormalization problem in nn scattering in refs .
@xcite .        and similar equations for higher low energy parameters , @xmath370 .
note that the potential only enters explicitly in the equation involving the discretized scattering length .
one of the nice features of this equation is the way how it handles the case of singular points , when @xmath371 or other parameters diverge , since standard integration methods for the corresponding differential equation are based on the smoothness of the solution and hence fail .
the way to proceed in practice is quite simple . for a given number of grid points @xmath358
we take the scattering length , @xmath372 , the effective range , @xmath373 and run eq .
( [ eq : alpha_dis ] ) and eq .
( [ eq : r_dis ] ) to determine @xmath374 and @xmath375 .
then , we use , e.g. @xmath376 and run eq .
( [ eq : m_dis ] ) upwards to obtain @xmath377 . due to the reversibility of the algorithm one
_ exactly _ has @xmath378 in the limit @xmath379 for any finite grid size @xmath380 ( modulo computer arithmetic round - off errors ) . in this way we can fix the initial boundary conditions to exactly reproduce any given scattering length , effective range , etc .
this method was used successfully  @xcite to extract low energy threshold parameters in all partial waves with @xmath381 from high quality potentials  @xcite . on that case
, we ran eq .
( [ eq : alpha_dis ] ) and eq .
( [ eq : r_dis ] ) with trivial initial boundary conditions @xmath382
, @xmath383 , etc , and obtain the potential s threshold parameters as @xmath384 , @xmath385 , etc .
as we have said the discretized running parameters located at the lowest possible radius @xmath386 correspond to the short distance interactions .
nowhere in the equations does the potential at the origin appears explicitly , since according to eq .
( [ eq : m_dis ] ) , one starts with @xmath387 . from this point of view
the treatment of regular and singular potentials at the origin is on equal footing .
the previous equations provide the accumulated phase shifts due to the addition of equidistant delta - shells sampling the original potential in coordinate space .
we want to show that they actually provide a solution of the decimation problem of the lse where the low energy states are cutoff . assuming that we have the lse with a given cutoff @xmath26 related to potential , we ask how does the physical phase - shift change when we make the transformation @xmath388 , in the potential fixed . applying
the process explained above based on nyquist theorem we see that this corresponds to double the grid resolution @xmath389 .
obviously , by repeating the process we may effectively have @xmath390 ( being @xmath332 the range of the potential ) and hence the potential never contributes and @xmath391 .
thus , for @xmath392 , the short range theory is recovered . in the opposite limit @xmath393
we get the full short range plus long range physics .
the discretized representation of the lse is algebraically closed , since it corresponds to a separable interaction , but any term contains _ all _ powers in momentum .
moreover , they also vanish for large @xmath50 , unlike any polynomial approximation to it .
the numerical procedure explained above of computing the phase - shift from some low energy parameters like @xmath397,@xmath55 , etc . and
the long distance potential @xmath278 can be also carried out by assuming that the potential at the lowest grid point , @xmath398 , is energy dependent yielding @xmath399 expanding in powers of momentum we can identify the terms @xmath400\end{aligned}\ ] ] with @xmath17 , @xmath18 , etc .
note that the radius of convergence of the expansion in @xmath50 and @xmath401 is the whole complex plane , due to the meromorphic character of the spherical bessel functions .