arxiv:hep-ph/ v1 4 Jul 2005

Transcription

1 Freiburg-THEP 05/06 hep-ph/ arxiv:hep-ph/ v 4 Jul 005 Two-Loop Bhabha Scattering in QED R. Bonciani and A. Ferroglia Fakultät für Mathematik und Phyik, Albert-Ludwig-Univerität Freiburg, D-7904 Freiburg, Germany Abtract In the context of pure QED, we obtain analytic expreion for the contribution to the Bhabha cattering differential cro ection at order α 4 which originate from the interference of two-loop photonic vertice with tree-level diagram and from the interference of one-loop photonic diagram amongt themelve. The ultraviolet renormalization i carried out. The IR-divergent oft-photon emiion correction are evaluated and added to the virtual cro ection. The cro ection obtained in thi manner i valid for on-hell electron and poitron of finite ma, and for arbitrary value of the center of ma energy and momentum tranfer. We provide the expanion of our reult in power of the electron ma, and we compare them with the correponding expanion of the complete order α 4 photonic cro ection, recently obtained in 0. A a by-product, we obtain the contribution to the Bhabha cattering differential cro ection of the interference of the two-loop photonic boxe with the tree-level diagram, up to term uppreed by poitive power of the electron ma. We evaluate numerically the variou contribution to the cro ection, paying particular attention to the comparion between exact and expanded reult. Key word: Feynman diagram, Multi-loop calculation, Bhabha cattering PACS:.5.Bt,.0.D Roberto.Bonciani@phyik.uni-freiburg.de Andrea.Ferroglia@phyik.uni-freiburg.de

2 Introduction The Bhabha cattering proce play a crucial role in the tudy of elementary particle phenomenology, ince it i the proce employed in the luminoity meaurement at e e collider. The mall-angle Bhabha cattering at high 00 GeV center of ma energy and the large-angle Bhabha cattering at intermediate 0 GeV energie have cro ection that are large and QED dominated; thee two characteritic allow for precie experimental meaurement, a well a for a detailed theoretical evaluation of the cro ection. The radiative correction to the Bhabha cattering in pure QED have been extenively tudied ee and reference therein. The Oα 3 correction have been known for a long time, even in the full Electroweak Standard Model. The econd order QED correction, Oα 4, were the ubject of renewed interet in the lat few year, and everal work have been devoted to the tudy of econd order radiative correction, both virtual and real, enhanced by factor of ln n /m with n =,, the c. m. energy, and m the ma of the electron 3, 4, 5. The complete et of thee correction wa finally obtained in 6 by employing the QED virtual correction for male electron and poitron of 7, the reult of 8, and by uing the known tructure of the IR pole in dimenional regularization 9. Very recently, the complete et of photonic Oα 4 correction to the cro ection that are not uppreed by poitive power of the ratio m / were obtained in 0. The ubet of virtual correction of Oα 4 involving a cloed fermion loop, together with the correponding oft-photon emiion correction, were obtained in an analytic and non approximated form in,, 3. The reult contained in thee paper do not rely on any ma expanion, and they are valid for arbitrary value of the c. m. energy, momentum tranfer t, and of the electron ma m. In,, 3, the calculation of the relevant loop diagram wa performed by employing the Laporta-Remiddi algorithm 4 which take advantage of the integration by part 5 and Lorentz-invariance 6 identitie in order to reduce the problem to the calculation of a mall et of mater integral. The mater integral are calculated uing the differential equation method 7; their expreion i given in term of harmonic polylogarithm 8. Both IR and UV divergencie were regularized in the dimenional regularization cheme 9 and they appear, in the intermediate reult, a ingularitie in D 4, where D i the dimenion of the pace-time. It i natural to apply the ame approach to the calculation of the complete et of Oα 4 virtual correction. At preent, the lit of the mater integral required to complete the calculation i available in 0,. However, only very few of the mater integral related to the two-loop photonic box diagram have thu far been calculated, 0,. The ubet of econd order radiative correction due to the interference of one-loop diagram wa tudied in 3. In uing the reult of 4, 5, we calculate in thi paper the following Oα 4 correction to the differential cro ection in QED: correction due to the interference of the two-loop photonic vertex diagram with the tree-level amplitude; correction due to the interference of one-loop photonic diagram amongt themelve;

3 correction due to the emiion of two real oft photon from a tree-level diagram and of one oft photon from a one-loop photonic diagram. All of the contribution mentioned above are calculated for finite electron ma and for arbitrary value of the c. m. energy and momentum tranfer t. We employ dimenional regularization in order to regularize both UV and IR divergencie. The UV renormalization i carried out in the on-hell cheme. By following the ame technique of 3, we pair virtual and oft-photon emiion correction in order to check the cancellation of the IR ingularitie. We alo expand our reult, that retain the full dependence on the electron ma, in the limit in which the electron ma i negligible with repect to the Mandeltam invariant. In thi way, it i poible: to prove that our reult reproduce the correct mall-angle Bhabha cattering cro ection at Oα 4, which i determined by the Dirac vertex form factor only 4; to provide trong cro-check of a large part of the reult of 0; to find, by ubtracting our reult from the cro ection of 0, the contribution to the cro ection of the interference between two-loop photonic boxe and the tree-level amplitude; to invetigate the numerical relevance of the term uppreed by poitive power of the electron ma. Thi paper i tructured a follow: after a brief ummary of our notation in Section, in Section 3, we dicu the irreducible two-loop vertex photonic correction, providing an expreion for their contribution to the virtual differential cro ection. In Section 4, 5, and 6, we calculate the interference of the reducible vertex diagram with the tree-level amplitude, the interference of the one-loop vertex diagram with themelve and with the one-loop box diagram, repectively, and we obtain the correponding contribution to the virtual differential cro ection. In Section 7, we complete the analyi of the interference among one-loop diagram by conidering the interference of one-loop boxe. In Section 8, we dicu the oft-photon emiion at Oα 4 and, in Section 9, we explicitly how how the cancellation of the IR divergencie work between virtual and oft correction. In Section 0, we analyze the expanion of our reult in the limit m / 0; we dicu the behavior of the cro ection at mall angle, compare our calculation with the reult preent in the literature, and dicu the numerical accuracy of the expanion. Section contain our concluion. In Appendix A, we collect the definition of ome function introduced in the paper. Finally, in Appendix B, we provide the expreion of the contribution of the two-loop photonic boxe to the differential cro ection at order α 4 in the limit m / 0; all of the function of the Mandeltam variable introduced and employed throughout the paper are available in electronic format in 6.

4 Kinematic, Notation and Convention In thi paper, we employ the notation and convention adopted in,, 3, which are ummarized in the preent ection. We conider the Bhabha cattering proce: e p e p e p 3 e p 4, where p, p, p 3, and p 4 are the momenta of the incoming electron, incoming poitron, outgoing electron, and outgoing poitron, repectively. All of the external particle are on their ma-hell, i. e. p i = m, i =,...,4, where m i the electron ma. The Mandeltam invariant, t, and u are related to the beam energy E and cattering angle in the center of ma frame of reference θ by the relation P p p = 4E, t Q p p 3 = 4 E m in θ, 3 u V p p 4 = 4 E m co θ. 4 Moreover, the Mandeltam invariant atify the relation t u = 4m. The analytic expreion for the the vertex and box form factor that we employ in the ret of the paper are calculated in the non phyical kinematic region < 0, and are then analytically continued to the phyical region > 4m. In order to expre the reult for < 0 in a compact form, it i convenient to introduce the dimenionle variable x, y, and z, defined through the following relation: m x, x = x t m y, y = y u m z, z = z 4m 4m, 0 x, 5 4m t t 4m t t, 0 y, 6 4m u u 4m u u, 0 z. 7 The analytic continuation of all the form factor to the phyical region > 4m can be readily obtained by replacing x x iǫ, where ǫ i an infiniteimal poitive quantity and where x 4m = 4m, 8 ee for a more detailed dicuion. The Bhabha cattering differential cro ection, calculated by umming over the pin of the final tate and averaged over the pin of the initial one, can be expanded in power of the fine tructure contant α a follow:, t, m = 0, t, m i=v,s α π i, t, m α π i, t, m, 9 3

5 p p 3 Q a b c d Figure : Two-loop irreducible diagram contributing to the photonic vertex correction of the electron current in the t-channel. The two diagram, analogou to Fig c and d, which have an inertion of a one-loop vertex and elf-energy correction on the outgoing electron line, are not hown. where the upercript V and S indicate virtual and oft-photon emiion contribution, repectively, while the ubcript 0,, and label the tree-level, Oα 3 and Oα 4 correction, repectively. The tree-level Oα and Oα 3 correction are well known their explicit expreion are collected, for example, in, 3. The contribution of Oα 4 to the virtual cro ection can be further plit a follow: V, t, m = V,phBoxe, t, m V,phVertice, t, m, t, m, 0 where the upercript N F = indicate the UV-renormalized diagram including a cloed fermion loop calculated, together with the correponding oft radiation diagram, in, 3. The upercript ph Vertice indicate the contribution of all of the UVrenormalized photonic correction which include at leat one vertex diagram, while the upercript ph Boxe indicate the contribution to the cro ection of the photonic correction which include box diagram only. V,ph Vertice Thi paper i dedicated to the calculation of /, a well a to the V,ph Boxe calculation of the contribution to / deriving from the interference of oneloop box diagram amongt themelve; the oft-photon emiion correction that cancel the reidual IR pole in the quantitie mentioned above i alo dicued. By employing thee reult and the finding of 0, we alo obtain the contribution of the interference of the two-loop photonic boxe to the cro ection at order α 4, up to term of order m / excluded. V,N F = 3 Irreducible Two-Loop Vertex Correction In thi ection, we obtain the contribution of the irreducible photonic vertex correction to the Bhabha cattering differential cro ection at order α 4. Four of the two-loop irreducible vertex graph, which correct the electron current in the t-channel photon exchange contribution to the Bhabha cattering, are hown in Fig.. There are two other graph contributing to the proce: thee are the pecular image of diagram -c and -d with the fermionic arrow revered, and their contribution to the differential cro ection i identical to that of the diagram in -c and -d. The um of the graph -a and -b, and of the graph -c and -d are eparately gauge independent. 4

6 a b c d Figure : Two-loop photonic vertex diagram contributing to Bhabha cattering. The explicit expreion of diagram -a -d were calculated in 5. The electron current which include the two-loop photonic correction can be written a α Γ µ p, p 3 = F π l,ph tγ µ m F l,ph tσ µν p p 3, where σ µν = i/γ µ, γ ν. The electron pinor and the dependence of the form factor F l,ph i t i =, on the electron ma m are omitted in Eq.. The expreion of the contribution of the ingle diagram to the UV-renormalized form factor F l,ph i t can be found in 5, and in 6. The form factor hown in Eq. till include IR ingularitie, which are regularized within the dimenional regularization cheme. The Laurent expanion of the form factor i t = F l,ph, t F l,ph, t D 4 D 4 F l,ph t = F l,ph, t D 4 F l,ph F l,ph,0 t OD 4, F l,ph,0 t OD 4, where D i the dimenional regulator. The four photonic two-loop vertex correction diagram contributing to Bhabha cattering are hown in Fig., where the haded circle repreent the um of the vertex graph hown in Fig.. The interference of diagram -b and -d with the tree-level amplitude provide a contribution to the differential cro ection that i identical to the one of diagram -a and -c, repectively. Thi can eaily be proved by oberving that diagram -b -d can be obtained from -a -c by applying the tranformation p p 3 and p 4 p, and that thee tranformation leave the Mandeltam invariant and t unchanged. The contribution of the diagram in Fig. -a to the Bhabha cattering differential cro ection at order α 4 can be written a V,ph Vertice = α phirr.ver., a t V l,irr.ver., t l,irr.ver. tv, t. 3 The function V l,irr.ver. i, t i =, have the Laurent expanion V l,irr.ver. i, t = V l,irr.ver., i, t V l,irr.ver., i, t D 4 D 4 V l,irr.ver.,0 i, t OD 4, 4 The diagram in Fig. -c and -d enter the um with a multiplicity factor of. 5

7 with V l,irr.ver., i, t = c i, tref l,ph, t, 5 V l,irr.ver., i, t = c i, tref l,ph, t c i, tref l,ph, t c i3, tref l,ph, t, 6 V l,irr.ver.,0 i, t = c i, tref l,ph,0 t c i, tref l,ph, t c i4, tref l,ph, t c i3, tref l,ph,0 t c i6, tref l,ph, t. 7 The function c ij are polynomial of the Mandeltam variable; their explicit expreion are collected in Appendix A, while the Laurent coefficient of the form factor F l,ph i i =, were introduced in Eq.. Even if the form factor are real for a phyical pace-like t, we write ReF l,ph i t i =, in the equation above for convenience of later ue. The contribution to the differential cro ection of the interference of the diagram in Fig. -c with the tree-level amplitude can be obtained from the contribution of diagram -a by replacing p p 3. Thi i equivalent to exchange t, o that one find V,ph Vertice = α phirr.ver., c l,irr.ver. V t, t V l,irr.ver. t,. 8 According to the definition in Eq. 5,6,7, the r. h.. of Eq. 8 involve the function F l,ph i, which develop then an imaginary part above threhold > 4m. Thee imaginary part do not contribute to the differential cro ection at thi order. Finally, the total contribution of the four diagram in Fig. to the Bhabha cattering differential cro ection at order α 4 i V,ph Vertice = α l,irr.ver. phirr.ver. V t, l,irr.ver. tv, t V l,irr.ver., t V l,irr.ver. t,. 9 t The firt, econd, and third term within quared bracket in the r. h.. of Eq. 9 are the -, t-t, and -t channel interference amplitude, repectively. 4 Reducible Two-Loop Vertex Correction In thi ection, we conider the interference of the two-loop reducible diagram, hown in Fig. 3, with the tree-level Bhabha cattering amplitude. The two reducible diagram contain a one-loop vertex correction in both fermionic current. Their contribution to the differential cro ection at order α 4 can be written in term of the one-loop UVrenormalized vertex form factor F l i i =, ee, 6 and of the function c ij i =,,j =, 6 introduced in the previou ection. It mut alo be oberved that the term linear in D 4 i needed in the Laurent expanion of the one-loop vertex form factor. 6

8 a b Figure 3: Two-loop reducible vertex diagram contributing to Bhabha cattering. The contribution of diagram 3-a to the cro ection i given by V,ph Vertice = α phred.ver.,3 a t V l,red.ver., t l,red.ver. tv, t, 0 where the Laurent expanion of V l,red.ver. i i V l,red.ver. i = V l,red.ver., i V l,red.ver., i D 4 D 4 V l,red.ver.,0 i OD 4, with V l,red.ver., i, t = c i, t ReF l, t ImF l, t, V l,red.ver., i, t = c i, t ReF l, t ImF l, t c i, t ReF l, tref l,0 t ImF l, timf l,0 t c i3, t ReF l, tref l,0 t ImF l, timf l,0 t, V l,red.ver.,0 i, t = c i4, t ReF l, t c i, t ReF l, c i6, t ReF l, c i, t ReF l, c i3, t ReF l, c i, t ReF l,0 t c i3, t c i5, t ImF l, tref l,0 t ImF l, timf l,0 tref l,0 t ImF l, timf l,0 tref l, t ImF l, timf l, tref l, t ImF l, timf l, ImF l,0 ReF l,0 tref l,0 t ImF l,0 timf l,0 ReF l,0 t 3 t t t t t t t t. 4 ImF l,0 Similarly, the contribution of diagram 3-b to the cro ection at order α 4 i given 7

9 by V,ph Vertice phred.ver.,3 b = α t V l,red.ver. t, l,red.ver. V t,. 5 It i then poible to conclude that the interference of the diagram in Fig. 3 with the tree-level amplitude give the following contribution to the cro ection V,ph Vertice phred.ver. = α l,red.ver. V t, l,red.ver. tv, t V l,red.ver., t V l,red.ver. t,. 6 t The firt, econd, and third term within quared bracket in the r. h.. of Eq. 6 are the -, t-t, and -t channel interference amplitude, repectively. 5 Interference of One-Loop Vertex Diagram In thi ection, we obtain the Bhabha cattering cro ection at order α 4 deriving from the interference of the diagram in Fig. 4. We begin by conidering the contribution of the amplitude of diagram 4-a quared; one find that V,ph Vertice phver.ver.,4 a = α t I, t. 7 In Eq. 7, we introduce the function I that ha the following Laurent expanion in power of D 4: I, t = I, t D 4 I, t D 4 I0, t OD 4, 8 with, t = c, t ReF l, t ImF l, t, 9, t = c, t ReF l, t ImF l, t c, t ReF l, tref l,0 t ImF l, timf l,0 t c 3, t ReF l, tref l,0 t ImF l, timf l,0 t, 30 I 0, t = c, t ReF l, tref l,0 t ImF l, timf l,0 t c, t ReF l, tref l, t ImF l, timf l, t c, t ReF l, tref l,0 t ImF l, timf l,0 t c 3, t ReF l, tref l, t ImF l, timf l, t I I 8

10 a b c d Figure 4: One-loop vertex diagram contributing to Bhabha cattering. c, t ReF l,0 t ImF l,0 t c 3, t ReF l,0 tref l,0 t ImF l,0 timf l,0 t c 7, t ReF l,0 t ImF l,0 t. 3 The expreion of the function c ij in term of the Mandeltam invariant can be found in Appendix A, while the expreion of the one-loop vertex form factor appearing in the function I can be found in, 5 and are collected in 6. Alo the contribution of the quare of the diagram in Fig. 4-c can be eaily expreed by employing the function I : V,ph Vertice phver.ver.,4 c = α I t,. 3 The interference between the diagram in Fig. 4-a and Fig. 4-b generate a term in the cro ection that can be written a V,ph Vertice phver.ver.,4 a b where the function I i given by = α t I, t, 33 with I, t = I, t D 4 I, t D 4 I0, t OD 4, 34 I, t = I, t, 35 I, t = I, t, 36 I 0, t = I0, t c 5, t c 7, t ReF l,0 t ImF l,0 t. 37 Similarly, the interference between the diagram in Fig. 4-c and Fig. 4-d generate the following contribution to the differential cro ection V,ph Vertice phver.ver.,4 c d 9 = α I t,. 38

11 The interference of the diagram in Fig. 4-a and Fig. 4-c can be expreed in term of a third function, I 3 : V,ph Vertice phver.ver.,4 a c = α t I 3, t. 39 Alo in thi cae, it i convenient to explicitly write the Laurent expanion of the function I 3 : with I 3, t = I 3, t D 4 I 3, t D 4 I0 3, t OD 4, 40 I 3, t = c t, F l, tref l,, 4 I 3, t = c t, F l, tref l, I 0 c t, c 3 t, F l, tref l,0 F l,0 tref l, F l, tref l,0 F l,0 tref l, 3, t = c 4 t, F l, tref l, c t, F l, tref l,0 F l,0 tref l, c t, F l, tref l, F l, tref l,, 4 c 6 t, F l, tref l,0 c 6, tf l,0 tref l, c 3 t, F l, tref l, c 3, tf l, tref l, c 3, tf l,0 tref l,0 c 3 t, F l,0 tref l,0 c 7 t, F l,0 tref l,0 c t, F l,0 tref l,0. 43 Finally, it i eay to verify that all of the interference between - and t-channel diagram in Fig. 4 give a contribution to the cro ection identical to the one in Eq. 39, while the quared amplitude of the diagram 4-b 4-d coincide with the quared amplitude of diagram 4-a 4-c. We can then conclude that the total contribution of the interference of the diagram in Fig. 4 to the Bhabha cattering cro ection i given by V,ph Vertice phver.ver. = α t I, t I, t 4 t I 3, t I t, I t, Interference of One-Loop Vertex and One-Loop Box Diagram In the preent ection, we conider the interference of the diagram in Fig. 4 with the diagram in Fig. 5. 0

12 a b c d Figure 5: One-loop box diagram contributing to Bhabha cattering. In order to write down the contribution that thee interference provide to the Bhabha cattering cro ection at order α 4, it i convenient to introduce the function I 4,i x, x, x 3, where i =,..., 3 and where x j j =,...,3 repreent one of the Mandeltam invariant,t, and u; thee function have the following Laurent expanion in D 4: I 4,i x, x, x 3 = I 4,i x, x, x 3 I 4,i x, x, x 3 D 4 D 4 I 0 4,i x, x, x 3 OD The coefficient of the Laurent expanion are given by I 4,i x, x, x 3 = x 3 Re F l, l, x B i x, x 3, 46 I 4,i x, x, x 3 = x 3 Re F l, l,0 x B i x, x 3 F l,0 l, x B i x, x 3 F l,0 l, x B i3 x, x 3, 47 I 0 4,i x, x, x 3 = x 3 Re F l, l, x B i x, x 3 F l,0 l,0 x B i x, x 3 F l, l, x B i x, x 3 F l,0 l,0 x B i3 x, x 3 F l, l, x B i3 x, x Beide for the one-loop vertex form factor already employed in the previou ection, one encounter in the equation above the function B l,j i i =,...,3, j =, 0, introduced in to decribe the contribution of the diagram in Fig. 5 to the Bhabha cattering cro ection at Oα 3. The auxiliary function B l,j i i = 4,..., 6, j =, 0 are introduced for the firt time here, and their expreion in term of HPL and the dimenionle variable x, y, and z can be found in 6. The interference of diagram 4-a with diagram 5-a give the following contribution While the analytic expreion of the function B l, i expreion of the function B l, i of the function B l i of HPL and dimenionle variable can be found in 6. and B l,0 i are explicitly given in the are not. They are the coefficient of D 4 in the Laurent expanion Eq. 45 in, that were not needed in that context. Their expreion in term

13 to the differential cro ection: V,ph Vertice phver.box,4 a 5 a = α 4 t I 4,t,, t. 49 The notation ha been choen to clarify that Eq. 49 i related to the interference of two t-channel diagram. Similarly, it i poible to write all the interference of diagram 4-a and 4-c with the diagram in Fig. 5 in term of the function I 4,i : V,ph Vertice phver.box,4 a 5 b V,ph Vertice phver.box,4 a 5 c V,ph Vertice phver.box,4 a 5 d V,ph Vertice phver.box,4 c 5 a V,ph Vertice phver.box,4 c 5 b V,ph Vertice phver.box,4 c 5 c V,ph Vertice phver.box,4 c 5 c = α 4 t I 4,t, u, t, 50 = α 4 t I 4,t, t,, 5 = α 4 t I 4,3t, u,, 5 = α 4 t I 4,,, t, 53 = α 4 t I 4,3, u, t, 54 = α 4 I 4,, t,, 55 = α 4 I 4,, u,. 56 The interference of diagram 4-b 4-d with a diagram in Fig. 5 i identical to the interference of diagram 4-a 4-c with the ame diagram in Fig. 5. Therefore, the total contribution of the interference between one-loop vertex diagram and one-loop box diagram to the Bhabha cattering cro ection i twice the um of Eq ; namely V,ph Vertice = α phver.box t I 4,t,, t I 4, t, u, t t I 4,t, t, I 4,3 t, u, I 4,,, t I 4,3, u, t I 4,, t, I 4,, u, Interference of One-Loop Box Diagram With reference to Eq. 0, it i poible to further plit σ V,phBoxe in the um of two contribution; the firt originate from the interference of two-loop photonic box diagram and tree-level diagram, the econd originate from the interference of two one-loop box

14 diagram: V,ph Boxe = V,ph Boxe V,ph Boxe. 58 phboxbox LBox The two-loop photonic box diagram in the m 0 cae are at the moment till unknown, and the calculation of the econd term in the r. h.. of Eq. 58 currently remain a an open problem ee 0,. On the contrary, the calculation of the firt term in Eq. 58 i, in principle, traightforward. The interference of every pair of the one-loop box diagram hown in Fig. 5 provide a contribution to the differential cro ection of the form V,ph Boxe = α 4 C ijx, y, z, 59 phboxbox,ij where the indice i, j run over the diagram label i, j = a, b, c, d and where we have introduced new function C ij. The explicit expreion of the latter, already continued to the phyical region > 4m, i particularly long and can be found in 6. 8 Soft-Photon Emiion at Order α 4 All of the two-loop photonic correction dicued in the previou ection are UV renormalized, but they till include double and ingle pole in D 4. Thee ingularitie have an IR origin, and they can be eliminated by adding the contribution of the real oft-photon emiion diagram at order α 4 to the virtual cro ection. Before dicuing the oft correction to the Bhabha cattering differential cro ection of order α 4, the reader i reminded that the oft correction at order α 3, dicued in detail in 3, can be written in the factorized form α π S, t, m α D = 0, t, m S IR, 60 π where σ D 0, t, m i the tree-level cro ection obtained by calculating the trace of Dirac matrice in D dimenion and where S IR i defined a S IR 4 4 J j, J j = ǫ j p p j I j, 6 j= with ǫ = ǫ 4 = and ǫ = ǫ 3 =, and I j = Γ 3 D π D 4/ m D 4 4π ω d D k k 0 p k p j k. 6 The integral in Eq. 6 can be found in 7 ee alo Appendix A of 3; the integration over the momentum of the oft photon k i retricted to the region k = k 0 < ω, where 3

15 a b c d e f g h Figure 6: Example of double photon emiion from tree-level diagram and ingle photon emiion from one-loop diagram. ω i the cut-off on the energy of the unoberved oft photon. The expanion of σ D 0, t, m in power of D 4 i D 0, t, m = 0, t, m D 4 0, t, m D 4 0, t, m O D 4 3, 63 where σ 0, t, m i the well known tree-level cro ection, and 0, t, m 0, t, m = α = α { t t 4 4 t t } t m t, 64 4 t t. 65 The contribution of the - and t- channel diagram, and of their interference to σ 0, t, m and σ 0, t, m, i evident in Eq. 64,65. There are two different kind of oft-photon emiion diagram contributing to the real correction to the cro ection at order α 4 : i the tree-level diagram with the emiion of two oft photon ome diagram belonging to thi cla are hown in Fig. 6-a d, and ii the diagram which include a one-loop correction and the emiion of a oft photon from one of the external leg ee the example in Fig. 6-e h. Since the oft-photon correction in QED exponentiate, the contribution of the double emiion diagram to the Bhabha cattering differential cro ection i given by α π S,double, t, m = α 0 D, t, m S IR. 66 π The emiion of a ingle photon from one-loop diagram, interfered with the tree-level ingle photon emiion graph, amount to α π S,ingle, t, m α = π V,D, t, m S IR, 67 4

16 where σ V,D i the UV-renormalized virtual cro ection at order α 3. The upercript D wa introduced a a reminder that the Laurent expanion of σ V,D mut be known up to and including term linear in D 4; the reaon for thi i that thee linear term give rie to a finite contribution when multiplied by the ingle pole preent in S IR. The real correction of order α 4 which originate from ingle-photon emiion diagram and which include a fermionic loop a, for example, diagram 6-h, were calculated in 3; they include ingle IR pole that cancel againt the virtual correction of order α 4 that alo include a photon elf-energy inertion. Thi et doe not play any role to the preent dicuion and i thu ytematically ignored. With the above in mind, the one-loop virtual cro ection appearing in Eq. 67 can be written a V,D, t, m = V,D, t, m V,D, t, m l,v, 68 l,b where the ubcript V indicate the contribution of the interference between vertex graph and tree-level amplitude; B tand for the cro ection generated by the interference between one-loop boxe and tree-level diagram. Furthermore, we plit the term on the r. h.. of Eq. 68 a follow:, t, m = V, t, m l,j D 4 V,, t, m l,j V,D l,j O D 4, 69 with j = V, B. The ingular and finite part of the one-loop virtual cro ection, correponding to the firt term in the r. h.. of Eq. 68, are well known. Their expreion in term of Mandeltam invariant and HPL can be found in Eq. 43 and Eq. 49 of. The term proportional to D 4 which arie from box and vertex one-loop correction are given by, t, m V, l,b = α m ReB l,, t ReB l, t, B l, 3 u, t ReB l, u, m t ReB l,, t ReB l, t, B l, u, t ReB l, 3 u,, 70 and V,, t, m = α l, l,v V t, l, tv, t V l, t, V l,, t, t 7 repectively, with V l, i, t = c i, tref l, t c i, tref l,0 t c i4, tref l, t c i6, tref l,0 t c i3, tref l, t, i =,. 7 5

17 All of the function appearing in the r. h.. of Eq were introduced in the previou ection, and their expreion in term of Mandeltam invariant and HPL have been collected in 6. A i mentioned above, the calculation of the integral I j which appear in Eq. 6 ha been carried out in 3, up to term linear in D 4 excluded. At firt glance, it appear that the calculation of uch term i needed, ince they provide, in the limit D 4, a non vanihing contribution to both Eq. 66 and Eq. 67. However, it i poible to prove that thi i not the cae. In order to proceed with our proof, one need to plit the one-loop UV-renormalized virtual correction in an IR-divergent part and a finite reminder, a i exemplified in the following: V, t, m = l,j D 4, t, m V,0, t, m l,j. 73 l,j V, The Laurent expanion of S IR ha the form S IR = S IR D 4 S0 IR D 4S IR O D The cancellation of IR divergencie in the order α 3 cro ection guarantee that 3, t, m V,, t, m l,b 0, t, m S IR = l,v V, By employing Eq. 66,67 in combination with Eq. 63,73,74, one can prove that the non-vanihing term proportional to S IR appearing in the double emiion cro ection Eq. 66 i S,double, t, m 0, t, m S IR S IR, 76 while the non vanihing term proportional to S IR cro ection at order α 4 i S,ingle, t, m V, appearing in the ingle-photon emiion, t, m V,, t, m l,b l,v S IR. 77 Therefore, we can conclude that, due to Eq. 75, the non vanihing term proportional to S IR cancel out in the total real emiion cro ection at order α4, given by the um of Eq. 66,67. 9 Cancellation of the IR Singularitie The IR divergencie in the real correction at order α 4 Eq. 66,67 hould cancel the IR ingularitie preent in the virtual correction dicued in Section 3-7 and the one ariing from the interference of the yet unknown two-loop photonic box graph. 3 We remind the reader that the photon elf-energy diagram are IR-finite ee 3. 6

18 J J 3 J J 3 = IR finite 4J 3 J 3 = IR finite Figure 7: Example of the cancellation of the IR divergencie in the two-loop irreducible vertex correction firt line and in the interference of the one-loop direct box in the -channel with itelf econd line. J J J J = double IR pole ζ J J J J = double IR pole ζ J J J J = double IR pole ζ J J J J = double IR pole ζ Figure 8: Reidual IR pole proportional to ζ in the -channel cro ection. A in the cae of the cancellation of the IR divergencie of the virtual cro ection at order α 4 N F = dicued in 3, it i poible to organize the contribution to the cro ection in IR-finite block by pairing the virtual correction originating from a certain et of diagram with an appropriate ubet of the oft-photon emiion correction. A an example, in the firt line of Fig. 7 we illutrate the cancellation of the IR pole preent in the interference of the graph in Fig. -a with the t-channel tree-level diagram; in all the term in the l. h.., the product of two graph repreent the contribution of their interference to the Bhabha cattering differential cro ection. Once again, the gray circle repreent the um of the UV-renormalized two-loop photonic vertex correction to the electron current in the t-channel photon-exchange diagram. In the econd line of Fig. 7, we provide another example of the cancellation of the IR divergencie in the box-by-box ector. Similar relation can be found for all of the contribution to the cro ection at order α 4. A pecial cae i repreented by the virtual correction in which the coefficient of the IR pole, when calculated in the non-phyical region < 0, include HPL with two or more zeroe in the rightmot poition in the weight lit. When performing the analytic 7

19 J J double IR J J = pole ζ J J J J = double IR pole ζ J J J J J J = double IR pole ζ Figure 9: Reidual IR pole proportional to ζ in the -t-channel interference contribution to V,ph Vertice σ. continuation to the phyical region > 4m, thee HPL generate real term proportional to ζ, where ζ i the Riemann Zeta function. For example, by replacing the non phyical dimenionle variable x according to x x iǫ, one find that or H0, 0; x H0, 0; x iǫ = H0, 0; x 3ζ iπh0; x, 78 H, 0, 0; x H, 0, 0; x iǫ = H, 0, 0; x 3ζ iπh, 0; x. 79 The part of the IR pole, that i proportional to the ζ factor ariing from analytical continuation doe not cancel in the combination with the oft-radiation contribution that eliminate the other IR ingularitie. Thi kind of behavior, already oberved in 3 in the dicuion of the Oα 4 N F = cro ection, i again encountered in the cae illutrated in Fig. 8. In the total differential cro ection, the reidual pole of Fig. 8 cancel themelve out. Combining the variou contribution hown in Fig. 9, the reidual IR pole do not cancel. Reidual ingle pole proportional to ζ alo arie in the combination hown in Fig. 0 and. Clearly, uch reidual pole cancel once the contribution of the two-loop box graph i included in the cro ection at order α 4. The five two-loop photonic box topologie are hown in Fig.. 0 Logarithmic Expanion In order to check our calculation againt the reult available in the literature, the contribution to the cro ection at order α 4 decribed in the previou ection are expanded in the limit in which the ma of the electron i mall with repect to all of the Mandeltam 8

20 4 J J 3 J J 3 J 3 J J = ingle IR pole ζ 4 J J 4 J J 4 J 4 J J = ingle IR pole ζ Figure 0: Reidual IR ingle pole proportional to ζ in the interference of one-loop vertex by one-loop box diagram. 4J J = double IR pole ζ 4J J 3 J 3 J = ingle IR pole ζ 4J J 4 J 4 J = ingle IR pole ζ Figure : Reidual IR pole proportional to ζ in the interference amongt one-loop box diagram. 9

21 a b d d d Figure : Two-loop photonic box topologie. invariant, t, and u 4. In thi limit, it i cutomary to write the two-loop photonic cro ection a follow ee 6, 0: ph 0 i V,i S,i = δ ln δ 0 m ln δ 0 m O m. 80 In the equation above, i = ph Vertice, ph Boxe. Following the notation adopted in the previou ection, V,ph Vertice Vertice V,ph = V,ph Vertice phirr.ver. phred.ver. Vertice V,ph V,ph Vertice, 8 phver.ver. phver.box V,ph Boxe while σ wa introduced in Section 7. The correponding σ S,i cro ection were obtained by pairing virtual and oft-photon emiion contribution a decribed in Section 9. In everal point of our dicuion, we treed the fact that a non-approximated calculation of the contribution of the two-loop photonic box diagram to the cro ection i till miing. However, all of the coefficient δ i i =,, 0 in the expanion in Eq. 80 are completely known; the firt two can be found in 6, while δ 0 wa recently obtained in 0. Therefore, by employing the m / 0 limit of the reult preented here in combination with 6, 0, it i poible to indirectly obtain the m / 0 limit of the contribution of the yet unknown two-loop photonic boxe and correponding oft-photon emiion correction. In Appendix B, we report the expreion of uch a contribution, both before and after adding the correponding oft-photon correction. Moreover, the expanion in the m / 0 limit of all the contribution to the Bhabha cattering cro ection dicued in the preent paper can be found in 6. It i known that the mall-angle Bhabha cattering cro ection i completely determined by the Dirac vertex form factor 4. In particular, one find that for the virtual cro ection V,ph 0 θ 0 = 6 F l l,ph t 4F t, 8 where F l and F l,ph are the UV-renormalized vertex form factor already employed in thi paper. The IR pole preent in the form factor are eaily removed by adding the 4 Note that thi expanion i not valid for very mall cattering angle, correponding to t < m, and for almot-backward cattering, correponding to u < m. 0

22 oft emiion contribution. By introducing the IR-finite form factor one find that F l t = F l t J J 3, F l,ph t = F l,ph t J J 3, 83 ph 0 θ 0 l l,ph = 6 F t 4 F t, { = ln 9 4ω 4ω ξ ξ m ln 6 ln ln 6ζ 3 3ζ 93 4ω m 8 9 lnξ 4 ln ln ξ 4ω ln 7 6 ln ξ 9ζ ζ ζ ln 4 5 ζ ζ 3ln ξ 3ζ ln ξ 8 ln ξ 9 ln ξ 4ω ln 4 lnξ ln ξ 4ω 8 ln 4 lnξ 6 ln ξ O ξ}. 84 The variable ξ i defined a ξ = coθ, 85 with θ the cattering angle in the c. m. frame. By expanding the exact photonic correction involving vertex diagram in the θ 0 limit and by neglecting term proportional to the electron ma, we recover the expreion in Eq. 84, which agree with 0 and the theorem in 4. Conequently, thi repreent a non-trivial tet of our calculation. A expected, we oberve that the interference between one-loop box diagram and one-loop vertex correction doe not contribute to the mall-angle cro ection. The interference of one-loop box diagram amongt themelve ha a non-zero mall angle limit, and all of the reidual term are proportional to ζ. Thee reidual term cancel out once the contribution of the two-loop photonic boxe i added; for thi reaon they are excluded from the preent dicuion. The agreement with the reult of 0 i clarified in Fig. 3, where we plot a a function of the cattering angle θ the reult of 0 and the correction of order α 4 originating from vertex graph Section 3 6, plu correponding oft emiion contribution. It i eaily een that, at mall angle, the vertex correction completely determine the cro ection. In addition, the expanion of the interference of one-loop boxe provide another trong tet of both our calculation and the one dicued in 0. The interference of ome pair of the diagram in Fig. 5 give origin, in the ratio / 0, to term proportional to ln n ξ, n =,, 4 ; m =,. 86 ξ m

23 reult of 0 ph Vertice / Figure 3: Order α 4 Bhabha cattering differential cro ection divided by the cro ection in Born approximation, a a function of the cattering angle θ. The continuou line repreent the reult of 0, while the dahed line repreent the correction involving at leat a vertex graph. The beam energy i choen equal to 0.5 GeV and the oft-photon energy cut-off ω i et equal to E. θ It i poible to oberve that, in the um of all the one-loop box interference, uch term cancel out, and that they do not appear in the complete photonic cro ection at order α 4 0. The photonic correction to the Bhabha cattering cro ection are now known up to term of Om / excluded; it i poible to ue the correction calculated exactly in thi work to etimate the relevance of the Om / term. We define i,ph = V,i S,i = i,ph m O L, m t, m u, 87 where the index i = Vertice, Box Box repreent the contribution dicued in Section 3 6, and 7, repectively. In Fig. 4 and 5, we plot, a a function of the beam energy, the quantitie α i,ph D i = π i,ph L 0 α π. 88 It can be een from thee plot that the term proportional to the electron ma become negligible for value of the beam energy that are very mall with repect to the one encountered in practically all of the e e experiment. It i alo reaonable to expect that the term proportional to the electron ma are negligible in the correction due to the two-loop photonic boxe. In thi ene, the approximated cro ection obtained in 6, 0,, 3 hould be ufficient for all phenomenological tudie.

24 e-05.5e-05 D Vertice θ, E e-05 5e E GeV Figure 4: D Vertice a a function of the beam energy, for θ = 0 olid line and θ = 90 dahed line. The oft-photon energy cut-off i et equal to E. 4e-06 3e-06 D Box Box θ, E e-06 e e E GeV Figure 5: D Box Box a a function of the energy, for θ = 0 olid line and θ = 90 dahed line. The oft-photon energy cut-off i et equal to E. 3

25 reult of 0 ph Vertice / 0 ph Vertice d σ / Figure 6: Comparion of the vertex correction with the complete photonic cro ection at order α 4. E and ω a in Fig. 3. θ Finally, in Fig. 6 and 7, it i poible to compare the variou contribution to the Bhabha cattering differential cro ection known at preent with the complete photonic cro ection in the m / 0 limit 0. The dahed line in Fig. 6 correpond to the contribution to the cro ection defined in Eq. 87 for i = Vertice, plotted a a function of the cattering angle. The local minimum at θ 80 and the maximum in the backward direction are due to puriou term proportional to the monomial ζ ln, ζ ln m m ln ξ, and ζ ln m ln ξ. 89 Thee term are not preent in the complete cro ection they cancel out againt analogou contribution deriving from the interference of two-loop box diagram with the tree-level amplitude and from the interference one-loop box diagram amongt themelve. Removing them from the vertex correction, one obtain the dahed-dotted curve in Fig. 6, which i moother than the dahed one. The olid curve repreent the complete photonic cro ection. A i already oberved above, the vertex contribution reproduce the full reult in the mall-angle region. The dahed line in Fig. 7 repreent the quantity ph Boxe 0 = ph 0 ph Vertice 0 = phboxbox 0 phlbox In the equation above, σ phboxbox wa introduced in Eq. 87, while σ phlbox i the contribution to the cro ection of the two-loop photonic boxe interfered with the treelevel amplitude. The expreion of the latter i given in Appendix B. Alo in thi cae, the behavior at θ 80 and at θ 80 i dominated by the puriou term of the kind of Eq. 89; by removing them, we obtain the dahed-dotted curve. 4

26 reult of 0 ph Boxe / 0 ph Boxe d σ / Figure 7: Comparion of the box correction with the complete photonic cro ection at order α 4. E and ω a in Fig. 3. θ For completene, in Fig. 8 we plot the photonic 0 and N F = 3 contribution to the Bhabha cattering cro ection at order α 4. The dotted line repreent the photonic correction. The correction of Oα 4 N F = dahed line have, for thi choice of ω ω = E, an oppoite ign with repect to the photonic correction. However, it i neceary to ay that the Oα 4 N F = of 3 include large contribution proportional to ln 3 /m that cancel out once the contribution of the oft pair production i included. Of the latter, only the term proportional to ln n /m n =,, 3 are known ee 5; we checked that they cancel the ln 3 /m term of the Oα 4 N F = cro ection. The dahed-dotted line repreent the um of the Oα 4 N F = cro ection with the known term of the pair production correction 5. The olid line i the complete order α 4 QED Bhabha cattering cro ection, including photonic, N F = and pair production contribution. Concluion In the preent paper, we obtained analytic non approximated expreion for all of the photonic correction to the QED Bhabha cattering differential cro ection at econd order Oα 4, except for the one deriving from the interference of two-loop photonic box diagram with the tree-level amplitude; at preent, the integral neceary to the 5 The pair production, calculated in 5 up to term enhanced by ln n /m included, depend upon a cut-off on the energy of the oft electron-poitron pair. In the numerical evaluation of Fig. 8 we et lnd = 4Ω ln, 9 with lnd defined in 5 and Ω numerically equal to the oft-photon cut-off: Ω = ω. 5

27 reult of 0 N F = / 0 N F =P / 0 tot/ Figure 8: Photonic, N F =, and total contribution to the cro ection at order α 4. E and ω a in Fig. 3. The pair production cut-off i et equal to ω. θ calculation of the latter are not known. The calculation were carried out by retaining the full dependence on the electron ma m. The reult are valid for arbitrary value of the c. m. energy and momentum tranfer t. We included a dicuion of the oft-photon emiion at order α 4, employing dimenional regularization to handle the IR-divergent term. We proved that the term proportional to D 4 in the Laurent expanion of the photon phae-pace integral I j in Eq. 6 doe not contribute to the Bhabha cattering cro ection. After ubtracting the IR ingularitie by adding the contribution of the oft-photon emiion graph, we expanded the contribution to the cro ection dicued in Section 3 7 in the m / 0 limit. In thi way, it wa poible to cro-check large part of the reult of 0, a well a to efficiently tet our calculation. By ubtracting the contribution to the cro ection obtained in thi paper from the cro ection of 0, it wa alo poible to indirectly obtain the contribution due to the interference of the twoloop photonic diagram with the tree-level amplitude, up to term uppreed by poitive power of the electron ma. By comparing the non approximated reult with the correponding m / 0 limit, we explicitly checked that the contribution to the cro ection of the term proportional to poitive power of the ratio m / i negligible at high- and intermediate-energy e e collider. It i reaonable that the ame concluion applie to the contribution to the order α 4 cro ection involving the two-loop photonic box diagram. For what concern phenomenological tudie, the reult of 6, 0,, 3 therefore provide a complete expreion of the virtual and oft-photon emiion correction to the Bhabha cattering cro ection at order α 4. Neverthele, a future calculation of the two-loop photonic box graph for m 0 would repreent a very intereting reult in the field of multi-loop calculation. 6

28 Acknowledgment The author wih to thank H. Czyz, B, Tauk, and L. Trentadue for ueful dicuion and uggetion, a well a A. Penin for cro-check, dicuion, and for providing u with the expreion of Eq. 3 of 0 in electronic format. We are grateful to J. Vermaeren for hi kind aitance in the ue of the algebra manipulating program FORM 8. A C ij Function In thi appendix, we collect the explicit expreion of the c ij function employed in the main text: c, t = t 4m, 9 c, t = t t m t, 93 c 3, t = t 3 tm 34 t, 94 c 4, t = t, 95 4 c 5, t = t t 4m tm 3 4 t t3 8m, 96 c 6, t = t 4tm t, 97 4 c 7, t = 4 t t m, 98 m c, t = m t t, 99 c, t = t, 00 c 3, t = t t m, 0 c 4, t = 0, 0 c 5, t = 3t, 03 c 6, t = t, 04 c 7, t = t t 4m t 4m 3 4 t. 05 B Photonic Double Boxe We provide here the expanion of the interference of the two-loop photonic boxe with the tree-level amplitude in the limit m / 0. Below one can find the contribution of 7

29 the two-loop photonic boxe to the virtual cro ection a well a the ame contribution after the ubtraction of the correponding oft-photon correction. Without including oft radiation V,phLBox 0 = δ V,LBox,3 log 3 m δ LBox,0 O δ V,LBox, log δ V,LBox, m log m m, 06 where δ V,LBox,3 = δ V,LBox, = δ V,LBox, = { 7 ξ ξ ξ 7ξ 4 3 ξ ξ4 ln ξ 6 ξ 7 ξ 7 ξ3 7 } 3 ξ4 lnξ, 07 { ξ ξ 6 ξ 8ξ ξ 3 6ξ 4 ln ξ 3ξ 9ξ 9ξ 3 D 4 6ξ 4 lnξ ξ ξ 9 ξ63 4 ξ 6ξ 3 ζ 9 4 ξ7 4 ξ 9 4 ξ3 ξ 4 ln ξ 4 ξ3 ξ4 ln ξ 7 4ξ 69 4 ξ 9 ξ3 4ξ 4 lnξ ln ξ ξ 4 ξ 7 4 ξ3 3 ξ4 lnξ 8 ξ 45 8 ξ 33 } 8 ξ3 ξ 4 ln ξ, 08 { ξ ξ 8 6ξ 4ξ 6ξ 3 8ξ 4 ln ξ 4ξ ξ ξ 3 D 4 D 4 8ξ 4 lnξ 4 30ξ 33ξ ξ 3 ζ 0 5ξ ξ 5ξ 3 0ξ 4 ln ξ 4 9ξ3ξ 9ξ 3 4ξ 4 ln ξ ξ ξ 3ξ 3 0ξ 4 lnξ 4ξ3ξ 8ξ 3 8ξ 4 lnξ ln ξ 9 ξ 9 ξ 5 ξ3 4ξ 4 ln ξ 7 ξ 33 ξ 7 ξ3 ζ 8 5 ξ 4ξ 5 ξ3 38ξ 85 ξ

30 δ V,LBox,0 = 89ξ 3 46ξ ζ3ζ 4 ln ξ 3 ξ 4 ξ 6 ξ3 3 ξ4 ln 3 ξ 3 ξ 5 ξ 3 ξ3 ln ξ 3ξ 3ξ 3 ξ 4 ln ξli ξ 67 4 ξ49 ξ 67 4 ξ3 ξ 4 ln ξ ξ 7 ξ 4ξ 3 ξ 4 Li 3 ξ 4ξ 7 ξ ξ Li 3 3 ξ 5 ξ3ξ 5 ξ3 ξ 4 Li 3 ξ ξ 4 57 ξ 7 ξ 3 ξ3 48ξ 4 ζ lnξ 3 5 ξ 8ξ 9 ξ3 ξ 4 lnξ ln ξ ξ 8ξ ξ3 lnξ ln ξ 4 3 ξ 3 ξ 7 5 ξ3 ξ 4 lnξli ξ 4 ξ 49 4 ξ 3 ξ3 ξ 4 lnξ ξξ ξ3 ξ 4 ln ξ ln ξ 8 ξ 5 4 ξ 7 8 ξ3 ln ξ 9 6 ξ 5 6 ξ 3 ξ3 } 3 ξ4 ln 3 ξ, 09 { ξ ξ 4 36ξ 36ξ ξ 3 ζ 8 6ξ 4ξ 6ξ 3 D 4 8ξ 4 ln ξ4 8ξξ 8ξ 3 4ξ 4 ln ξ4ξ ξ ξ 3 8ξ 4 lnξ 8 6ξ 4ξ 6ξ 3 8ξ 4 lnξ ln ξ ξ 6ξ 6ξ 3 4ξ 4 ln ξ ξ 5ξ ξ 3 ζ4 0ξ69ξ 74ξ 3 40ξ 4 ζln ξ D 4 8 ξ6ξ ξ 3 8ξ 4 ln ξ ξ 3ξ ξ 3 ln ξ ξ ξ ξ 3 ln 3 ξ 4 3ξ57ξ 54ξ 3 40ξ 4 ζ lnξ ξ 8ξ 0ξ 3 8ξ 4 lnξ 0 ξ 8ξ lnξ ln ξ 4 9ξ 8ξ 3ξ 3 lnξ ln ξ ξ 3 ξ ξ3 ln ξ 6 7 ξ 7 ξ 9 5 ξ3 ln ξ ln ξ ξ 7 ξ 3 ξ3 ln 3 ξ 5 ξ 4 ξ 7 ξ3 3ξ 4 ξ 3ξ 3 ξ3 ζ ζ3 43 ξ ξ 9ξ 3 ξ3 4ξ 4 ζ3 ln ξ ξ 4 ζ ln ξ 49 5 ξ8ξ 39 5 ξ3 5 ξ4 ζ 4 ξ 57ξ 3 36ξ 4 73 ξ ξ73 48 ξ 5 4 ξ3 7 ξ4 ln 4 ξ 9 ζ ln ξ 4 ξ 7ξ 3 6 ξ 6 ξ

31 3 ξ3 3 ξ4 ln 3 ξ 5 5 ξ5 ξ3 5ξ 4 ln ξli ξ 5 4 ξ 9 4 ξ 5 4 ξ3 ξ4 ln ξ 5 ξ 5 ξ3 ξ 4 ln ξli ξ 6 9ξ 7 ξ 4ξ 3 8ξ 4 ln ξli 3 ξ ξ 3 ln ξli 3 ξ 8 7 4ξ 7 ξ ξ ξ4ξ 3 ξ3 ξ 4 ln ξli 3 ξ 8 9ξ 9 ξ 9ξ 3 8ξ 4 ln ξ ξ7ξ 6ξ 3 38ξ 4 ζ lnξ 6 5 ξ ξ 3 ξ3 4ξ 4 ζ3 lnξ 40 5 ξ 39 ξ 3 ξ3 4ξ ζ 4 lnξ ln ξ 5 5 ξ5 3 ξ 9 ξ3 8 3 ξ4 lnξ ln 3 ξ 3 4 ξ 8 ξ 4 ξ3 ξ 4 lnξ ln ξ 33 ξ 9ξ 3 ξ3 lnξ ln ξli ξ 43 4 ξ0ξ 3 ξ3 ξ 4 lnξ ln ξ 45 4 ξ ξ 0 4 ξ3 ξ 4 ζ ln ξ 49 8 ξ 5 4 ξ ξ3 ln ξ ln ξ 6ξ 33 8 ξ 4 ξ3 ξ 4 ln ξ ln ξ ξ ξ ξ ξ4 ln 3 ξ ln ξ 48 ξ 7 6 ξ ξ3 5 ξ4 ln 4 ξ 6 ξ 5 4 ξ 5 4 ξ3 3 ξ4 ln 3 ξ ξ 3 4 ξ ξ3 3ξ 4 ln ξli ξ 4 ξ 8 ξ 4 ξ3 ξ4 ln ξ 4 4ξ 3 ξ 3 ξ3 ξ 4 lnξli ξ 6 4ξ 3 ξ 9ξ 3 ξ 4 lnξli 3 ξ 4ξ 7 ξ ξ 3 lnξli 3 ξ 6 3 ξ ξ ξ 9 ξ3 4ξ 4 lnξli 3 ξ ξ 9 4 ξ 0ξ 3 8ξ 4 lnξ 3ξ 93 ξ 3 44ξ 3 3ξ ζli 4 ξ ξ 7 ξ 7 ξ3 ξ 4 Li 3 ξ 4ξ 7 ξ ξ 3 Li 3 ξ ξ ξ ξ3 ξ 4 Li 3 ξ 7ξ 9 ξ 4ξ 3 6ξ Li 4 4 ξ 6 4ξ 9 ξ 7ξ Li 3 4 ξ 7 ξ ξ 7 } ξ3 ξ Li 4 4 ξ. 0 30

32 Including oft radiation phlbox = δ LBox, log δ LBox, 0 m log δ LBox,0 m O where δ LBox, = δ LBox, = m, ξ ξ 8ξ 8ξ 6ξ 3 ζ, { ξ ξ D ξ 36ξ ξ 3 ζ 4ω ln ξ3ξ 3ξ 3 ξ 4 lnξ 4ξ6ξ 4ξ 3 ξ 4 ln ξ 4ω ln ξ 3 ξ 3ξ 3 ξ 4 ζ 5 ξ7 ξ 5 ξ3 ξ 4 ln ξ 3 ξ 9ξ ξ3 3ξ 4 ln ξ 5ξ 9ξ 7ξ 3 4ξ 4 Li ξ 3ξ ξ ξ 3 ξ 4 lnξ ln ξ ξ 9 ξ 9 3 ξ3 3ξ 4 lnξ 4 ξ 4 ξ 3 4 ξ3 ξ 4 ln ξ 9ξ 39 4 ξ 5 ξ3 3ξ 4 ζ ξ 4ξ 48ξ 3 4ξ 4 ζ ln ξ ξ 4 ξ 5 4 ξ3 3 3 ξ4 ln ξ 4 ξ 3 4 ξ3 ln ξ 3 5 ξ 7 ξ ξ3 6ξ 4 Li ξ 4 30ξ 36ξ 30ξ 3 4ξ 4 ζ lnξ 3 9 ξ 34 ξ 3ξ 3 3ξ 4 lnξ ln ξ 3 4 ξ lnξ 9 8 ξ 3 8 ξ 9 8 ξ3 3 } ξ4 ln ξ, 3 { δ LBox,0 = ξ ξ D ξ 36ξ ξ 3 ζ 6ξ ξ ξ 3 ζ ξ 4ξ 48ξ 3 4ξ 4 ζ ln ξ D ξ 36ξ 30ξ 3 4ξ 4 ζ lnξ 4ω ln ξ 3ξ ξ 3 ξ 4 ln ξ 4ξ 6ξ 4ξ 3 ξ 4 ln ξ 4ξ 6ξ 4ξ 3 ξ 4 lnξ ln ξ ξ 3ξ 3ξ 3 ξ 4 lnξ ξ 3 ξ 3 ξ3 ξ 4 ln ξ 3

33 4ω ln ξ 3 5 ξ 3ξ 3 ξ 4 ζ ξ 9ξ 3 6ξ 4 ζ ln ξ 5 ξ 7 ξ 5 ξ3 ξ 4 ln 3 ξ ξ 7 ξ ξ 3 ξ 4 ln ξ 4 8ξ ξ 8ξ 3 4ξ 4 ln ξli ξ 4 5 ξ ξ 5 ξ3 4ξ 4 ln ξ 5ξ 9ξ 5 7ξ 3 4ξ 4 Li ξ ξ 9ξ 3 6ξ 4 ζ lnξ 5 ξ ξ ξ3 ξ 4 lnξ ln ξ 5 3 ξ 5ξ ξ 3 ξ 4 lnξ ln ξ 5ξ 9ξ 7ξ 3 4ξ 4 lnξli ξ 3 ξ 6ξ 6ξ 3 4ξ 4 lnξ 5 4 ξ 3 4 ξ 4 ξ3 3ξ 4 ln ξ ln ξ 4 ξ 4 ξ 3 4 ξ3 ξ 4 ln ξ 3 4 ξ 4 ξ 3 4 ξ3 ξ 4 ln 3 ξ 7 ξ 7ξ 4ξ ξ ξ4 ζ 9 ξ 4ξ 3 ξ 4 ζ3 ξ 3 ξ ξ 3 7ξ 4 ζ ξ 37ξ 3 ξ 7ξ 59 ξ3 8ξ 4 ζ ln ξ ζ ln ξ 4 8 ξ 48 ξ 3 ξ3 ln 4 ξ 3 8ξ 3 ξ ln 3 ξ 5 9 ξ7ξ 3 ξ3 ξ 4 ln ξli ξ ξ 4ξ ξ ξ4 ln ξ ξ ξ3 ln ξli ξ 6 0ξ 0ξ 3 6ξ 4 ln ξli 3 ξ 6 6ξξ 4ξ 3 ln ξli 3 ξ 6ζ3ξ ξ ζ3ξ 4ζ3ξ 3 ξ 3 6ζ3 ln ξ 0ξ 5 4 ξ 6ξ 3 8ξ 4 ζ lnξ 8 36ξ 6 ξ 3ξ 3 ξ 4 ζ lnξ ln ξ ξ8 3 ξ 6 ξ3 3 ξ 4 lnξ ln 3 ξ 4ξ9 8 ξ 3 4 ξ3 lnξ ln ξ 8 0ξ ξ 7ξ 3 4ξ 4 lnξ ln ξli ξ ξ ξ 8ξ 3 5ξ 4 lnξ ln ξ 7 ξ 4ξ 5 ξ3 ξ 4 ζ ln ξ 7 4 ξ 5 8 ξ ξ 3 ξ4 ln ξ ln ξ 3 5 ξ9 4 ξ ξ 3 ln ξ ln ξ 3 7 ξ 3 4 ξ 5 ξ3 7 3 ξ4 ln 3 ξ ln ξ 48 ξ 5 48 ξ ξ3 6 ξ4 ln 4 ξ 5 3 ξ 3 ξ 5 4 ξ3 ln 3 ξ 4 ξ 7 4 ξ 4 ξ3 ξ 4 ln ξli ξ 3