arxiv: v3 [hep-ph] 8 Oct 2015

Transcription

1 Diperive approache for three-particle final tate interaction Peng Guo,,,, I. V. Danilkin,,, and Adam P. Szczepaniak Phyic Department, Indiana Univerity, Bloomington, IN 4745, USA Center For Exploration of Energy and Matter, Indiana Univerity, Bloomington, IN 4748, USA. Thoma Jefferon National Accelerator Facility, Newport New, VA 66, USA (Dated: October 9, 5) In thi work, we preent different repreentation of Khuri-Treiman equation and dicu advantage and diadvantage of each repreentation. In particular we focu on the inverion technique propoed by Paquier, which even though developed a long time ago, ha not been ued in modern analye of data on three particle decay. We apply the method to a toy model and compare the enitivity of thi and alternative olution method to the left-hand cut contribution. We alo dicu the meaning and applicability of the Waton theorem when three particle in final tate are involved. JLAB-THY-4-95 arxiv:49.865v [hep-ph] 8 Oct 5 PACS number: Keyword: I. INTRODUCTION Low energy hadronic reaction play an important role in contraining parameter of effective theorie of QCD. Furthermore preciion in the determination of the underlying reaction amplitude i needed when earching for phyic beyond the Standard Model. There are model independent retriction that at low energie help to determine hadronic amplitude. Thee include, for example, chiral ymmetry, partial wave and effective range expanion and unitarity. Thee contraint, together with the requirement that reaction amplitude repect croing relation and are analytical function of the kinematical variable are often ued to formulate diperion relation for the amplitude. A particular implementation lead to the o-called Khuri-Treiman (KT) equation [ 7]. Thee equation were originally written for amplitude involving four external particle, i.e amplitude decribing cattering or decay. Specifically, the KT equation impoe two-body unitarity on a truncated et of partial wave in the three channel of a fourpoint amplitude. In practical application, for example to the analyi of η π decay, the KT approach helped reduce the dicrepancy between the meaured decay width and the NLO chiral perturbation theory prediction [8 5]. With the availability of high-preciion data on production and decay of light hadron, there ha recently been a renewed interet in application of the KT framework. Different repreentation of KT equation and variou method for olving KT equation have been propoed in the pat [ 7]. To the bet of our knowledge, the comparion of thoe different variant of KT equation ha never been demontrated before. In thi paper, we begin with the dicuion of the variou repreentation and method of obtaining their olution. In particular, Electronic addre: pguo@jlab.org we focu on the Paquier inverion technique that wa introduced in [6, 7]. The olution of the KT equation require knowledge of two-body partial wave amplitude. Once thee amplitude are known, different method for olving the equation yield the ame anwer. However, when the Paquier inverion i applied the input twobody amplitude i probed in a different energy region than in the cae of a direct olution of the KT equation. Therefore, for phyical application (where the two-body partial wave amplitude are only known in a limited energy range), variou olution trategy are complimentary to each other and allow to tudy ytematic uncertaintie. Since the two-body partial wave amplitude are only known in a limited energy range, in the cae of the Paquier inverion method approximation are often needed when dealing with the left-hand cut. In [6] the author howed that baed on a realitic et of ππ partial wave and a ingle real parameter repreenting the left-hand cut, one could effectively reproduce the Dalitz plot ditribution in η π + π π. Given the number of unknown parameter in thi approximation, it could be beneficial to tudy the decay of higher ma particle, e.g. η, ω, φ, D etc. where there are no pecific low energy contraint. However, from the other ide it i important to illutrate the Paquier inverion technique and the correponding approximation in a chematic model where the exact olution i known. Thi will hed further light on it uefulne. We conider a toy model of calar particle interacting in a ingle partial wave. There are known cae, e.g. ω/φ π where a ingle partial wave dominate. Alo in our toy model, the left-hand cut i choen to be cloe to the phyical region. That i imilar to what happen with the ππ amplitude, where the left-hand tart at = and the threhold i located at 4 m π. The paper i organized a follow. In Section II we firt introduce kinematic and Khuri-Treiman (KT) equation needed for a diperion decription of the three-body decay. Then we dicu different olution trategie and compare numerical reult on the hand of the chematic

2 4 = a() (a) f * () 4 + FIG. : A diagrammatic repreentation of dicontinuity relation in Eq.(6). model. A a further analyi, we how how the olution from variou repreentation of the KT equation change depending on the left-hand cut approximation. In Section III we give a brief dicuion of the role of the Waton final tate interaction theorem. Summary and concluion are preented in Section IV. II. THE KHURI-TREIMAN MODEL A. Definition The reaction amplitude decribing interaction of four calar particle i = 4 i given by a ingle complex function, A(, t, u), of the three Mandeltam variable atifying + t + u = M + m. In the following i =,, 4 label identical particle with ma m and i = refer to a particle with ma M. We are pecifically intereted in the decay of particle i.e. the cae M > m. The amplitude A(, t, u) decribe four ditinct procee. Thee are (with the bar denoting an antiparticle) the -channel cattering, + + 4, the t-channel cattering + + 4, the u-channel cattering, and the particle decay channel, In term of particle momenta the three Mandeltam variable are = (p + p ), t = (p + p ) = (p p ), u = (p + p 4 ) = (p p 4 ). In the -channel, the partial wave (p.w.) expanion, A(, t, u) = 6π a(t) 4 (b) f * () (l + ) A l () P l (z ), () l= define the -channel partial wave, A l (), with z = co θ being the coine of the -channel center of ma cattering angle, z = (t u) λ / (, M, m )λ / (, m, m ). () The triangle function λ i given by λ(a, b, c) = a + b + c a b b c a c. Similar p.w. expanion can be written in the t- and the u-channel with the center of ma cattering angle given by t ( u) z t = λ / (t, M, m )λ / (t, m, m ), u (t ) z u = λ / (u, M, m )λ / (u, m, m ), () 4 repectively. A truncated p.w. erie define an amplitude which i a regular function of the cattering angle and it only ingularitie are with repect to the channel energy variable, i.e. the -channel p.w. expanion ha ingularitie in but not in z or equivalently t or u. Since A(, t, u) ha ingularitie in all three variable, it implie that p.w. erie in any channel diverge outide the phyical region of that channel e.g. the -channel p.w. expanion diverge outide the > (M + m) > 4m and z <. In other word, p.w. expanion in a pecific channel need to be analytically continued outide the phyical region where the erie i defined to obtain the phyical amplitude repreenting reaction in the other channel. In the Khuri-Treiman model ingularitie of A(, t, u) in all three variable are recovered by approximating the amplitude a the um of p.w. erie in the three channel imultaneouly [ 7], with each um truncated at ome finite value, l = L, A(, t, u) = 6 π L (l + ) l= [a l () P l (z ) + a l (t) P l (z t ) + a l (u) P l (z u )]. (4) The truncation on the number of partial wave alleviate the problem of double counting, ince the removed, high partial wave in one channel are being replaced by the low partial wave in the croed channel. In Eq. (4), a l denote the partial wave amplitude of the KT model and hould be ditinguihed from the A l in Eq.(). The relation between the two i obtained by projecting Eq. (4) on to the -channel partial wave, A l () = a l ()+ dz P l (z ) L (l + ) l = [a l (t(, z )) P l (z t (, z )) + (t u)]. (5) Each of the three term on the right hand ide of Eq.(4) ha ingularitie in one variable only i.e., t and u, repectively. Thee ingularitie are aumed to originate from unitarity, o that even though a l () ha only the right hand, unitary cut, A l () alo ha the left hand cut due to the exchange term. Truncation of the partial wave erie lead to an incorrect aymptotic behavior of A(, t, u) at large value of the Mandeltam variable and the KT model i intrinically limited to low-energie. For implicity, in the following analyi we truncate the partial wave to include S-wave only i.e. et L = and denote a() a (). At low energie, below the firt inelatic threhold, unitarity i aturated by two-particle intermediate tate, the dicontinuity relation of a() by croing unitarity cut in i thu given by, a() = (a( + iɛ) a( iɛ)) i [ = f () ρ() a() + K() t+() t () dt a(t) ]. (6)

3 Here ρ() = 4 m / i the two body phae pace factor and f() i the S-wave cattering amplitude decribing two-body interaction between pair of the particle, and. The firt term on the right hand ide of Eq.(6) correpond to the cae when particle + 4 are produced from particle + in -channel, S-wave projection of the unitarity relation. It i illutrated in the diagram in Fig. (a). One half of the econd term, illutrated in Fig. (b), give the contribution from the by t-channel exchange, and it i obtained by projecting the t-channel two-particle intermediate tate expreed through the t-channel partial wave erie (i.e. the econd term on r.h. of Eq.(4) truncated to include S-wave only) onto the -channel S-wave. The factor of in front of the integral take into account the contribution for the u-channel (lat term on the r.h. in Eq.(4)). The -channel partial wave projection of the cro channel partial wave erie i obtained by integrating the t-channel amplitude over z, which, uing Eq.() i expreed a an integral over t at fixed. The integration limit t ± () correpond to z = ±, t ± () = M + m ± K(), (7) where K() i given by [8], +κ(), 4 m (M m), K() = i κ(), (M m) (M + m), κ(), (M + m) < +, κ() = λ(, M, m ) λ(, m, m ) /. (8) For real and > 4m, the integration of t-argument in Eq.(6) follow the path hown in Fig.. In particular when i decreaed toward (M m) i.e. it approache the phyical boundary of the decay region + + 4, a expected, the integration path followed by dt approache the poitive real axi. Even though, however, at = (M m), t + () = t () the integration path remain finite a it run below the unitary cut, from point c on Fig. though point d to point g above the t- channel unitarity cut of a(t). Thi wa the key obervation made in [] where it wa alo hown to be conitent with perturbation theory. That i to ay, a() in the decay region i obtained by an analytical continuation from the cattering region. Thi i the cae becaue the integration path hown in Fig. avoid any ingularitie a i decreaed from the cattering, > (M + m) to the domain of the decay region 4m < < (M m). In ummary, the KT model i a low-energy approximation to the amplitude A(, t, u) in which the exchange force are approximated by a finite number of partial wave. Given the dicontinuity relation Eq.(6), a() can be contructed by uing Cauchy theorem (auming no ubtraction are needed), a() = π d a( ) 4m (9) i h b t + () 4m a d t () g c e t plane f (M m) FIG. : Motion of the upper, t +() and lower, t () limit of integration in Eq.(6) in the complex plane a a function of (real). The arrow indicate the direction of increaing in the interval from 4 m to. The zigzag line repreent the location of ingularitie (right hand cut) of a(t). The point labeled a through i correpond to pecific value of, with (a) t ( ) =, (b) t ((M + m) ) = m (m M), (c) t ((M m) ) = m (M + m), (d) t ( M m ) = 4 m, (e) t ±(4 m ) = M m, (f) t +(m (M + m)) = (M m), (g) t +((M m) ) = m (M +m), (h) t +((M +m) ) = m (m M), and (i) t +( ) =, repectively. with input to thi equation provided by the two-body partial wave, f(). B. Method for olving Eq.(9) If the exchange term in Eq.(6) i ignored, the reulting dicontinuity relation for a() i of the Mukhelihvili- Omnè type [7, 8] and can be olved uing the tandard N/D method [9, ]. Elatic unitarity determine the dicontinuity of the S-wave two-body cattering amplitude, f() = i (f( + iɛ) f( iɛ)) = ρ() f(). () The algebraic olution of Eq.() i f() = ei δ() in δ(), () ρ() where δ() i the S-wave phae hift. Writing the partial wave amplitude a f() = N()/D(), with D() containing the elatic cut and all other cut aborbed into the function N(), linearize Eq.() and yield an analytical parametrization for f(), with D() = exp ( π d δ( ) ) 4m (. () ) By convention we normalized D() o that D() =. The function /D() i referred to a the Mukhelihvili- Omnè (MO) function [7, 8]. Conitency with Eq.(-) require that N() i a regular function above

4 4 elatic threhold. From Eq.() it follow that the aymptotic behavior of D() for i given by δ( )/π. Writing a() = G() D(), () and ignoring the contribution to a() from the exchange channel, one find G() =, i.e. G() i an analytical function in the entire complex plane. Thu if Eq.(9) i conidered a an integral equation for the KT amplitude, in general, the olution i not unique ince it i parametrized by the cla of entire function G(). The only retriction on G() i that it i bounded by δ( )/π o that the integral in Eq.(9) converge. However, ince in the KT model the high energy behavior i not contrained Eq.(9) can be ubtracted arbitrary number of time, in which cae there i no retriction on the large- behavior of G(). It i therefore more appropriate to regard the diperion relation, like the one in Eq.(9) a a contraint on the amplitude rather than a a dynamical equation for the amplitude. Keeping the exchange term, with a() given by the full expreion of Eq.(6) and with the parametrization a() = G()/D(), the dicontinuity of the function G() i given by the exchange term [, ], and G() = ρ() N() K() t+() = ρ() f () D () K() G() = π dt G(t) D(t) t () t+() t () dt G(t) D(t), (4) 4m d G( ). (5) The main advantage of uing Eq.() i that the integral in Eq.(5) depend on the two-body partial wave amplitude f() in the phyical region, i.e. on the real axi above the elatic threhold where f() i entirely determined by the phae hift. The diadvantage i that computation involve a double integral, one in Eq.(4) to obtain the dicontinuity, G() and the other in Eq.(5). An alternative repreentation for a() i obtained by writing it in term of f() intead of the denominator, D() function, a() = f() g(). (6) If the left hand cut in f() i ignored, i.e. N() a approximated by an entire function then the dicontinuity of g() become proportional to that of G() ince, comparing Eq.() and Eq.(6) one find G() = N()g(). In general, however, g() mut aborb the cut of N() o that the latter are abent in a(), which by contruction ha only the unitary cut. One therefore find, g() = g L () + g R (), (7) where g R () = ρ() K() g L () = N() N () t+() t () dt g(t)f(t), (8) g(), (9) and the diperion relation for g() follow directly from Eq.(8),(9), g() = g L () + g R () = π L d g L( ) + π 4m d g R( ). () Here L mark the beginning of the left hand cut branch point of f(). At firt ight, it eem that repreentation given by Eq.(6) i more complicated compared to that of Eq.() ince the former require knowledge of the left-hand cut contribution to the two-body cattering amplitude. Even though precie knowledge of the left-hand cut i neceary for determining the low-energy behavior of the calar-iocalar ππ amplitude and, for example, propertie of the σ meon [, 4], it doe not need to be the cae when the left-hand cut i far away from the phyical region. In thi cae, the left-hand contribution to the phyical region could be well approximated by a uitably choen conformal expanion [5 7]. In the following ection, we tet the enitivity of the KT equation olution to the particular form of the left-hand cut. We alo remark that Eq.(6) idetep computation of the MO function, which would require knowledge of amplitude phae at high energie. Neverthele, it till require an analytical repreentation for the S-wave twobody cattering amplitude f() in the complex plane, which appear under the integral for g R () in Eq.(8). In both cae, either uing Eq.() or Eq.(6), evaluation of a() require computation of a double integral and the equation for G() or g() can be olved only by iteration. The double integral i of the type I() = 4m d A( ) t+() t () dtb(t), () where B(t), which i proportional to the unknown function G(t) or g(t), ha a cut on the real axi for t > 4m. The integral over run over the real axi and the t integral run over the complex contour Γ hown in Fig.. A hown in [6, 8] the two contour can be deformed imultaneouly in uch a way that the real integration i deformed onto a contour C in the complex plane while the t complex contour i brought onto the contour Γ along the real axi with the location of t along Γ depending on the location of on C. More detail of the deformation of contour are ummarized in Appendix A. Thi procedure allow to interchange the order of and t integral, and it i referred to a the Paquier inverion. It reult in the following expreion of I(), I() = dt B(t) Γ C (t) d A( ). ()

5 5 The integral over C i now independent of the function B(t) and for pecific A( ) the integration can be done analytically or numerically reulting in a kernel function K A (, t), giving I() = dt B(t) K A (, t). () Γ Uing thi method computation of G() or g() reduce to a ingle integral equation which can be olved either by iteration or by matrix inverion. The Paquier inverion eem to favor the repreentation of Eq.(6) over that given in Eq.(). Thi i becaue the diperion relation for g R (), with B(t) = g(t)f(t) and A( ) = ρ( )/K( ) reult in a univeral kernel, K g i.e. one that doe not depend on the input two-body amplitude, K g (, t) = θ(t) g (, t) θ( t) Σ g (, t) (4) and the function g and Σ g are given by Eq.(A7) and (A8) repectively. The integral equation for the function g() i then given by, g() = π + π L d N( ) N ( ) g( ) dt g(t) f(t) K g (, t). (5) On the other hand, the function A( ) entering the diperive repreentation for G(), doe depend on the input f() amplitude, A( ) = f( ) ρ( ) D( )/K( ) = N( ) ρ( )/K( ), leading to where G() = π dt G(t) D(t) K G(, t), (6) K G (, t) = θ(t) G (, t) θ( t) Σ G (, t) (7) and the function G and Σ G, given by Eq.(A4) and (A5), depend on the input two-body amplitude. In thi cae the kennel K G (, t) i model dependent and in general ha to be computed numerically. C. Model analyi In the previou ection we dicued two common parametrization of the KT partial wave amplitude. Thee parametrization are ditinguihed by how the elatic cut i implemented; either through the MO function, a in Eq.(), or uing the partial wave two body cattering amplitude, a in Eq.(6). It i alo poible to compute a() directly from Eq.(6) and (9), i.e. without factoring out the elatic cattering amplitude contribution to the KT amplitude. Given input partial wave amplitude, f() and boundary condition (cf. ee dicuion in Sec. II B) the olution i obtained by numerically olving an integral equation. The method baed on the Paquier inverion i potentially bet uited for analyi of large data et, ince it reduce the problem to a one-dimenional integral equation that can be olved by matrix inverion. Furthermore for thi method the repreentation of Eq.(6) i mot natural ince the kernel function are univeral in thi cae. However, the input involve p.w. two-body amplitude outide the phyical region. The amplitude in the phyical region contribute only to the econd integral on the r.h. of Eq.(5) in the interval between 4m and (M m). The olution baed on Eq.(6) alo depend on the extrapolation outide the phyical region. The integral over t involve the MO function below threhold and the kernel function K G (, t) depend on the twobody amplitude evaluated along the complex contour C (cf. Eq.(A4) and (A5)). Since the extrapolation of the two-body amplitude i largely-model dependent it i of interet to tudy enitivity of the variou repreentation dicued above to model of the left-hand cut that determine the input amplitude f() outide the phyical region. In thi ection we ue a imple analytical model for f() to analyze thi enitivity. The model amplitude we ue i dominated by a ingle reonance and incorporate the left hand cut in an analytical form. Specifically we chooe, f() to be given by, f() = α m R α I R() β I L (), (8) where I R,L () = I R,L () ReI R,L (m R ), o that the real part of the denominator vanihe at the reonance ma, = m R. The function I R i given by the integral over the two-body phae pace I R () = d 4m / π 4m (. (9) ) I R () i reponible for the right hand cut and guarantee unitarity. The contribution from the left hand cut i modeled after a Born-term exchange contribution to the invere amplitude with a threhold at 4m in the croed channel, which reult in the left hand cut in the direct channel partial wave tarting at L = I L () = π 4m ln 4m. () For large-, f() = O(/) and the phae hift reulting from thi model, δ = tan (Imf/Ref) approache π. For a particular choice of parameter, α =., β =., m R =.8 GeV, and m =.4 GeV and M =.4 GeV, the phae hift, and the function N() and D () are hown in Fig.. We numerically olve Eq.(9), (5) and (6). The boundary condition on the olution i impoed by demanding a() = f() which we obtain by ubtracting the

6 a() FIG. 4: Real part (olid) and imaginary part (dahed) of the function a() obtained by olving Eq.(5). Alo hown are the olution, red and blue quare, obtained by olving directly the integral Eq.(9) and Eq.(6) repectively. N() D - () FIG. : Phae hift (top pane), real part (olid curve) and imaginary part (dahed curve) of N() (middle pane) and /D() (bottom pane) correponding to the model amplitude of Eq.(8). diperion relation at =. A expected, regardle of repreentation, the three method give the ame reult a illutrated in Fig. 4. We tudy the enitivity to the left hand cut by firt analyzing it effect on the exact olution. We do thi by varying the parameter β in Eq.(8), which control it trength, i.e. β = reult in an amplitude f() without left hand ingularitie. Thi dependence i illutrated in Fig. 5. Even though the left-hand cut contribution ha ome effect on g() = a()/f() in phyical decay region, which i a large a about 5% the phae of a() i barely affected by left-hand cut of f() (only about % change), a een in lower plot in Fig. 5. In the following we how how the olution from different repreentation of the KT equation change depending on left-hand cut approximation. Since the KT amplitude can be repreented in term of either D() or f() we expect a difference in the KT amplitude a() depending on how the left-hand cut i approximating in the two repreentation. In the extreme cae, which we tudy here, in olving the diperion relation for G() or g() we retrict the integral to be a much a poible determined by the phyical region of the two-body amplitude. That i, intead of Eq.(5) we firt ue g() = π dt g(t)f(t)k g (, t). () where, a a imple approximation, we parametrized contribution from g L by a contant. In principle thi contant hould be fitted to the data. However, here we tudy the extreme cae when ubtraction contant i alway fixed by a() = f() that correpond to a cenario

7 g().6. a() FIG. 5: Real part (olid) and imaginary part (dahed) of the function, g() = a()/f() with β varied in the range [,.]. The colored curve correpond to β = (black),. (red),. (blue). The phae of a() (lower pane), φ() = tan (Ima()/Rea()), the color cheme i ame a in upper panel. when the contribution from g L i et to be zero. In econd tudy we further limit the range of integration to only include the phyical region, g() = π 4m dt g(t)f(t)k g (, t). () A far a tudying enitivity of the repreentation in term of D() to the unphyical region in the two-body amplitude in Eq.(6) we retrict the integration range to FIG. 6: Real part (olid) and imaginary part (dahed) of the function a() (upper pane) with different approximate cheme for β =., black, blue, red and green curve repreent full olution, olution of Eq.(), olution of Eq.() and the olution of Eq.() repectively. The phae of a() (lower pane), φ() = tan (Ima()/Rea()), the phae are hifted by a contant, o that φ(m R) = π, the color cheme i ame a in upper pane. the phyical region a well, G() = π 4m dt G(t) D(t) K G(, t). () The reult are hown in Fig. 6. The overall hape of the olution eem to be quite inenitive to the lefthand cut, alo ee the phae of olution in lower plot in Fig. 6. The normalization can vary by a much a 5%, however, thi can be reduced by changing the ubtraction contant. Indeed, the ubtraction contant

8 8 g ().5 /(m R t) ().5 (a) (b) FIG. 9: (a) A diagrammatic repreentation of Eq.(5), (b) A triangle diagram with a fixed internal ma m R in () ub-channel, thi diagram may be aociated to the function d ρ( ) t+ ( ) 4m K( ) t ( ) dt. m R t FIG. 7: The comparion of real part (olid) and imaginary part (dahed curve ) of a() (black) given by Eq.(6), and ρ()f ()a() (red). a() i the olution of Eq.(9) with β = FIG. : The comparion of real part (black olid curve) and imaginary part (black dahed curve) of g pole given by Eq.(5), and the triangle diagram (red curve), with α =. and m R =.8GeV. The real part of olution are hifted up by a unit for the better diplaying purpoe. FIG. 8: The comparion of phae of f(), φ = tan (Imf/Ref) (dahed black), and phae of a() (olid black) by etting β =.. The φ i normalized to φ(m R) = π for comparion purpoe. a() may be computed approximatively in perturbation theory (uch a χpt) or may be adjuted to empirical data directly. Thi i an important reult a it how that unitarity in the phyical region, where it can be contrained by the data, play a key role in determining the olution of the KT equation. III. WATSON S THEOREM IN A THREE-BODY FINAL STATE Amplitude involving two particle in the final tate atify Waton theorem, which tate that the phae of the correponding partial wave amplitude in the elatic region coincide with the elatic phae hift [9]. Thi i a traightforward conequence of unitarity. The phae relation doe not hold, however, even in the elatic region, when more than two particle are produced and there are interaction among all produced particle [ ]. A a imple illutration, one can conider three particle production in a two-body colliion: a + b c + d + e. Partial wave in the cd two-particle ubytem are labeled by

9 9 the total pin, J and helicity Λ. They depend on the total center of ma energy quared, cd of the cd ubytem, momentum tranfer t be between particle b and e and total invariant ma quared, ab. In the phyical limit, of poitive, and ay large- ab and mall- cd, i.e. in the cd elatic region, Waton theorem would relate the phae of the cd partial wave amplitude with pin J to that of the phae hift of the pin-j, cd elatic cattering amplitude. However, ince ab i alo in the phyical region, phae of the amplitude alo depend on contribution from dicontinuitie in the ab channel, i.e. all intermediate procee a + b X c + d + e, which couple at a given value of ab. In particular, in the limit of large- ab thi phae can be determined from Regge-theory [4]. The nonequivalence between imaginary part of the amplitude and it dicontinuity wa oberved early on when unitary relation for multiparticle amplitude were invetigated, for example in [5, 6]. For the recent dicuion we refer to [7 9]. There i ignificant interet in uing diperion relation to extract two-particle phae-hift from amplitude analyi of three-particle final tate, e.g. produced in pion or photon diffraction or in heavy flavor decay [4 46]. In general, three-body recattering effect between reonance decay product and pectator modify the dicontinuity relation that follow from two-body unitarity and a a conequence lead to violation of Waton theorem. We can ue the pecific model decribed in the previou ection tudy to invetigate the ize of thee effect. Uing the olution from Eq.(9) with β =., in Fig. 7, we compare a() given by Eq.(6) and ρ()f ()a(). A can be clearly een, a() doe have an imaginary part. It arie from the t and u-channel exchange, i.e. recattering between a pion emerging from a decay of a t or u channel reonance an the pectator. We alo plot the phae of a() compared to the phae of f() in Fig. (8). The phae difference of a() and f() i given by the phae of g() = a()/f(). Above the elatic threhold, the phae difference i due to the crochannel exchange. c.f. Eq. 6. Near the reonance, the phae difference between a() and f() i minor, which i conitent with the expectation that three-body effect are at mot a correction when two particle emerge from a reonance. Since a() we oberve that a() acquire a non-zero phae below the phyical threhold 4 m. The non-vanihing phae of a() can be undertood by examining the aymptotic behavior of the reduced amplitude g(). Even in the cae when only the right-hand cut of f() i included, a, we obtain, g() π 4m d ρ( ) K( ) t+( ) t ( ) dt f(t) g(t). (4) In general, f() i a complex function, thu, the phae of g() cannot be zero below 4m a i een in Fig. 8. The cup at threhold originate from the triangle (anomalou) ingularity [47 5] and it i caued by the reduced amplitude g(), which can be een in both Fig. 4 and Fig. 5. To examine thi point, in Eq.(6), we replace f() by a implified expreion, α/(m R ) and obtain, g pole () = π ρ( ) K( ) d 4m t+( ) t ( ) dt α m R t g pole (t). (5) g pole () may be identified with the re-cattering contribution from a perturbative analyi of the triangle diagram a hown in Fig. 9. In perturbation theory, we approximate g pole (t) = under the integral on the righthand ide of Eq.(5) and thi lead to the triangle diagram of Fig. 9(b). When m R > 4 m the contour integration Γ(t) will weep through the real axi and pick up an aborptive part reulting in imaginary part of g pole () below threhold < 4 m generating a cup near threhold. In Fig., we how the olution of Eq.(5), g pole () and compare it to the firt recattering correction (a triangle diagram). IV. SUMMARY We analyzed different repreentation of the Khuri- Treiman equation and demontrated, on the bai of a imple model, that olution from variou repreentation are indeed identical. The main focu of our tudy wa the Paquier inverion, which until recently ha never been ued in the numerical analyi of the three particle decay. Therefore, we tudied the enitivity of the variou repreentation to the left-hand cut contribution. We howed that the Paquier inverion might be indeed a good complementary approach to the direct olution of the KT equation, epecially when the phyical region doe not depend trongly on the accurate form of the left-hand cut. In the end, we dicued the Waton theorem when the interaction among three particle in the final tate i involved. We concluded that the phae of the olution of KT equation did not coincide with the elatic phae hift of cattering amplitude. V. ACKNOWLEDGMENTS We thank Michael R. Pennington for many ueful dicuion, we alo thank Andrew W. Jackura for the final proof reading of thi manucript. Thi reearch wa upported in part by the U.S. Department of Energy under Grant No. DE-FG87ER465, the Indiana Univerity Collaborative Reearch Grant and U.S. National Science Foundation under grant PHY-59. We alo acknowledge upport from U.S. Department of Energy contract DE-AC5-6OR77, under which Jefferon Science Aociate, LLC, manage and operate Jefferon Laboratory.

10 a c 4m d (t) t plane f (M m) b e (M m) h (M+m) + (t) C plane FIG. : The contour Γ (olid purple line) and C (olid black curve) and in Eq.(A). The arrow indicate the direction that invariant follow along the path of integration. The purple and black wiggle line repreent unitarity cut in t-plane and cut attached to two branch point: (M ± m) in -plane repectively. The point labeled by a f correpond to (a) () =, (b) (4 m ) = M m, (c) ( M m ) = 4m, (d) ±((M m) ) = m (m + M), (e) +(m (M + m)) = (M m), (f) +(4 m ) = M m, (g) +() =, (h) +(m (m M)) = (M + m), and (i) +() =, repectively. Appendix A: Kernel function In a double integral equation of the type given by Eq.(6), I = d N( ) 4m U( ) t+( ) (Γ) dt a(t), t ( ) g i (A) the contour Γ followed by t integration i defined in Fig. and the integration path over i defined on the real axi. The quare root function U(z) i given by U(z) = (z (M m) ) (z (M + m) ) in the complex-z plane. The phae convention for U(z) i choen by U( ± i) = (, i, ±) U() for ([, (M m) ], [(M m), (M + m) ], [(M + m), ]) repectively. In Eq.(A), the phae of U() i choen a the value of U() right below the two cut attached to branch point (M ± m), i.e. K()/ρ() = U( i). A decribed in [6, 8], when the order of the and t integral i revered, Eq.(A) become, [ ] I = dt a(t) (C ) d N( ) Γ Γ (t) U(. (A) ) The contour C and Γ avoid the ingularitie in the integrand, ee Fig.. Whether Γ (t) i + (t) or (t) depend on whether t i above or below the cut in t- plane repectively, and ± (t) are given by the olution of φ( ±, t) =, where φ(, t) = t u m (M m ) and + t + u = M + m. Splitting the integration path, [ = Γ ] (M m) in Eq.(A) (ubcript +/ of integration limit denote the path of integration lying above or below the cut attached to branch point (M m) in t-plane, ee Fig. ), one obtain, I = [ θ(t) +(t) (t) dt a(t) (C ) θ( t) (C ) +(t) ] d N( ) U( ). The kernel function G and Σ G are defined by, G (, t) = Σ G (, t) = +(t) (C ) d (t) (C ) d +(t) and Eq.(A) finally become I = N( ) U( ), N( ) U( ). dt a(t) (θ(t) G (, t) θ( t) Σ G (, t)). (A) (A4) (A5) (A6) In general, kernel function G and Σ G have to be evaluated numerically by contour integration in the complex plane. In particular for the cae N() =, the correponding kernel, which we denote a g and Σ g can be expreed in term of elementary function g (, t) = U() ln R(, t) + U() U(t) R(, t) U() U(t) and θ (φ(, t)) i π U(), (A7) Σ g (, t) = U() ( + (t) ) ( M m + U() ) ln ( + (t) ) ( M m ) + U () U() U( + (t)). where R(, t) = M 4 + ( m )(t m ) + M ( + t). (A8) (A9) For real and t the phyical value of g and Σ g correpond to the limit + i and t + i. [] N. N. Khuri and S. B. Treiman, Phy. Rev. 9, 5 (96). [] J. B. Bronzan and C. Kacer, Phy. Rev., 7

11 (96). [] I. J. R. Aitchion, II Nuovo Cimento 5, 44 (965). [4] I. J. R. Aitchion, Phy. Rev. 7, B7 (965); Phy. Rev. 54, 6 (967). [5] I. J. R. Aitchion and R. Paquier, Phy. Rev. 5, 74 (966). [6] R. Paquier and J. Y. Paquier, Phy. Rev. 7, 94 (968). [7] R. Paquier and J. Y. Paquier, Phy. Rev. 77, 48 (969). [8] J. Kambor, C. Wieendanger and D. Wyler, Nucl. Phy. B 465,5 (996). [9] A. V. Aniovich and H. Leutwyler, Phy. Lett. B 75,5 (996). [] S. Lanz, PoS CD, 7 () [arxiv:.78 [hepph]]. [] K. Kampf, M. Knecht, J. Novotny and M. Zdrahal, Phy. Rev. D 84, 45 (). [] S. Decote-Genon and B. Mouallam, Eur. Phy. J. C 74, 946 (4). [] J. Bijnen and J. Gaer, Phy. Scripta T 99 4 (). [4] I. V. Danilkin, C. Fernndez-Ramrez, P. Guo, V. Mathieu, D. Schott, M. Shi and A. P. Szczepaniak, Phy. Rev. D 9, 949 (5). [5] S. P. Schneider, B. Kubi, C. Ditche, JHEP, 8 (). [6] P. Guo, I. V. Danilkin, D. Schott, C. Fernndez-Ramrez, V. Mathieu and A. P. Szczepaniak, Phy. Rev. D 9, 546 (5). [7] N. I. Mukhelihvili, Tr. Tbili. Math Intrum., (958); in Singular Integral Equation, J.Radox, ed. (Noordhoff, Groningen, 985). [8] R. Omnè, Nuovo Cim. 8, 6 (958). [9] G. F. Chew and S. Mandeltam, Phy. Rev. 9, 467 (96). [] G. Frye and R. L. Warnock, Phy. Rev., 478 (96). [] S. Mandeltam, Phy. Rev. Lett. 4, 84 (96). [] W. R. Frazer and J. R. Fulco, Phy. Rev. Lett., 65 (959). [] I. Caprini, G. Colangelo and H. Leutwyler, Phy. Rev. Lett. 96, (6). [4] R. Garcia-Martin, R. Kaminki, J. R. Pelaez and J. Ruiz de Elvira, Phy. Rev. Lett. 7, 7 (). [5] F. J. Yndurain, hep-ph/8. [6] A. Gaparyan and M. F. M. Lutz, Nucl. Phy. A 848, 6 (). [7] I. V. Danilkin, L. I. R. Gil and M. F. M. Lutz, Phy. Lett. B 7, 54 (). [8] I. J. R. Aitchion and J. J. Brehm, Phy. Rev. D 7, 7 (978). [9] K. M. Waton, Phy. Rev. 88, 6 (95). [] I. J. R. Aitchion, II Nuovo Cimento 5, 44 (965). [] I. J. R. Aitchion, II Nuovo Cimento A 5, 49 (967). [] R. D. Amado, Phy. Rev. 58, 44 (967). [] I. J. R. Aitchion and C. Kacer, Phy. Rev. 7, 7 (968). [4] R. C. Brower, C. E. DeTar and J. H. Wei, Phy. Rept. 4, 57 (974). [5] R. Blankenbecler, Phy. Rev., 98 (96). [6] L. F. Cook, Jr. and B. W. Lee, Phy. Rev. 7, 8 (96). [7] S. P. Schneider, B. Kubi and F. Niecknig, Phy. Rev. D 86, 54 (). [8] B. Ananthanarayan, I. Caprini and B. Kubi, Eur. Phy. J. C 74, no., 9 (4). [9] I. Caprini, Phy. Rev. D 9, 44 (5). [4] E. M. Aitala et al. [E79 Collaboration], Phy. Rev. D 7, 4 (6); Erratum-ibid. D 74, 599 (6). [4] G. Bonvicini et al. [CLEO-c Collaboration], Phy. Rev. D 78, 5 (8). [4] I. Caprini, Phy. Lett. B 68, 468 (6). [4] The FOCUS Collaboration: J. M. Link et. al., and M.R. Pennington, Phy. Lett. B 65, (7). [44] P. C. Magalhãe et. al., Phy. Rev. D 84, 94 (). [45] F. Haa on behalf of the COMPASS Collaboration, AIP. Conf. Proc. 74, 7 (). [46] S. Uhl on behalf of the COMPASS Collaboration, arxiv:4.494[hep-ex]. [47] P. V. Landhoff and S. B. Treiman, Phy. Rev. 7, 649 (96). [48] R. J. Eden, P. V. Landhoff, D. I. Olive, and J. C. Polkinghorne, The Analytic S-Matrix (Cambridge Univerity Pre, Cambridge, England, 966). [49] J. -J. Wu, X. -H. Liu, Q. Zhao, and B. -S. Zou, Phy. Rev. Lett. 8, 88 (). [5] Q. Wang, C. Hanhart and Q. Zhao, Phy. Lett. B 75, 6 (). [5] Q. Wang, C. Hanhart and Q. Zhao, Phy. Rev. Lett., ().