arxiv: v3 [hep-ph] 28 Nov 2012

Transcription

1 PREPARED FOR SUBMISSION TO JHEP Baryon electric dipole moments from strong CP violation arxiv:.5887v3 hep-ph] 8 Nov Feng-Kun Guo a and Ulf-G. Meißner a,b a Helmholtz-Institut für Strahlen- und Kernphysik and Bethe Center for Theoretical Physics, Universität Bonn, D-535 Bonn, Germany b Institute for Advanced Simulation, Institut für Kernphysik and Jülich Center for Hadron Physics, JARA- FAME and JARA-HPC, Forschungszentrum Jülich, D-545 Jülich, Germany fkguo@hiskp.uni-bonn.de, meissner@hiskp.uni-bonn.de ABSTRACT: The electric dipole form factors and moments of the ground state baryons are calculated in chiral perturbation theory at next-to-leading order. We show that the baryon electric dipole form factors at this order depend only on two combinations of low-energy constants. We also derive various relations that are free of unknown low-energy constants. We use recent lattice QCD data to calculate all baryon EDMs. In particular, we find d n =.9 ±.9 and d p =. ±. in units of 6 e θ cm. Finite volume corrections to the electric dipole moments are also worked out. We show that for a precision extraction from lattice QCD data, the next-to-leading order terms have to be accounted for. KEYWORDS: Chiral Lagrangians, CP violation, Lattice QCD

2 Contents Introduction Effective Lagrangian for strong CP violation 3 3 Baryon electric dipole form factors in the infinite volume 5 3. Baryon electric dipole form factors up to NLO 6 3. Numerical results for the loop contributions of the baryon EDMs A first determination of the LECs from lattice results 3.4 Counterterm-free relations 3 4 Finite-volume corrections in the p-regime 5 5 Summary 8 A Baryon masses up to NLO B Loop integrals in infrared regularization C Expressions for the baryon EDFFs up to NLO D Finite volume corrections to loops 4 E Asymptotic expansion of finite volume corrections 5 Introduction The neutron electric dipole moment EDM) is a sensitive probe of CP violation in the Standard Model SM) and beyond. The current experimental limit d n.9 6 e cm ] is still orders of magnitude larger than the SM prediction due to weak interactions. Furthermore, in quantum chromodynamics QCD) the breaking of the U) A anomaly allows for strong CP violation, which is parameterized through the vacuum angle θ. Therefore, an upper bound on d n allows to constrain the magnitude of θ. Furthermore, such electric dipole moments are very sensitive to physics beyond the SM, see e.g. ]. Many extensions of the SM in fact lead to larger EDMs than the tiny SM predictions, so that any limit of d n leads to bounds on the scale of the new physics. In this paper, we concentrate on the CP violation generated by the θ-term of QCD. New and on-going experiments with ultracold neutrons strive to improve the aforementioned bounds even further, see e.g. 3] for a recent review. Furthermore, there are new experimental proposals to measure the EDM of the proton and the deuteron in storage rings, which in principle allow for an even higher sensitivity than obtained with the instable neutron, see e.g. refs. 4 7].

3 Besides these challenging experimental activities, first full lattice QCD calculations of the neutron and the proton electric dipole moment are becoming available. There exist three different methods of calculating the nucleon EDM on the lattice. The EDM can be related to the energy difference of the nucleon with different spin alignments in the presence of an external electric field 8 ]. It can also be obtained by calculating the electric dipole form factor EDFF) at finite momentum transfer q, see the definition in eq. 3.), and extrapolating to the point with q =, 3], d N = F 3,N) m N..) In addition, the nucleon EDM can also be calculated by analytically continuing θ to a purely imaginary quantity 4, 5], as the QCD action in the presence of the θ-term becomes real in Euclidean space. For a brief review of these methods, see ref. 6]. We do not discuss here another method 7] that relates certain Fock state components of the EDM to the Fock state expansion of the magnetic moment derived in light-front QCD) 8] since it is not clear what these relations imply for the observable quantities 9, ]. These lattice studies require a careful study of the quark mass dependence of the nucleon EDM to connect to the physical light quark masses. In addition, CP-violating atomic effects can be sensitive to the nuclear Schiff moment, which receives a contribution from the radius of the nucleon electric dipole form factor EDFF), see e.g. ]. It is thus of paramount interest to improve the existing calculations of these fundamental quantities in the framework of chiral perturbation theory CHPT). In ], the electric dipole moments of the neutron and the Λ were calculated within the framework of U3) L U3) R heavy-baryon chiral perturbation theory and an estimate for θ was given for earlier works utilizing chiral Lagrangians, see 3 5]). In 6], the electric dipole form factor of the nucleon was analyzed to leading one-loop accuracy in chiral SU). In that calculation, the form factor originates entirely from the pion cloud. The strength of the form factor was shown to be proportional to a non-derivative and CP-violating pion nucleon coupling ḡ πnn that was estimated from dimensional analysis in 6]. In ref. 7], the results of refs., 6] were extended to higher order based on a covariant version of U3) L U3) R baryon CHPT, and an expression for ḡ πnn was given in terms of the measurable quantities. For other recent work on these issues, see, 8]. Here, we extend the studies of ref. 7] to the baryon octet and make contact to recent lattice QCD studies. Furthermore, the leading contributions to the neutron EDM at finite volume and in partiallyquenched calculations were considered in 9], and in 3] the leading order extrapolation formula using a mixed action chiral Lagrangian is given. Here, we work out the finite volume expression for the whole baryon octet at next-to-leading order NLO). In particular, we show that the formally suppressed NLO corrections are sizeable, and even dominate for some baryons due to large cancellations at leading order. The manuscript is organized as follows. The underlying effective chiral U3) L U3) R Lagrangian is given in section. In section 3 we work out the baryon EDFFs and EDMs in the infinite volume for varying quark masses. We derive new relations for all baryon EDMs and show that these only depend on two combinations of unknown low-energy constants LECs). Using recent lattice data, we can give predictions for all baryon EDMs. Then, in section 4, we calculate the finite volume corrections for the baryons at NLO and show that these NLO corrections are substan-

4 tial and must be included in any extraction from lattice data. We end with a summary and outlook in section 5. Various technical aspects of our calculations are displayed in the appendices. Effective Lagrangian for strong CP violation The most general gauge-invariant and renormalizable Lagrangian for QCD is L QCD = 4 Ga µνg a,µν + q i D M) q + θ g 3π Ga µν G a,µν a =,..., 8),.) with G a µν the gluon field strength tensor and G a µν = ε µνλσ G a,λσ / its dual, q collects the quark fields of various flavors, D µ is the gauge-covariant derivative, and M is the quark mass matrix. The last term is the so-called θ-term, which breaks the P and CP symmetries. It is a consequence of the U) A anomaly. Because the theta-term is related to chiral U) transformations of the quark fields see, for instance, ref. 3]), only the combination θ = θ + arg det M.) is a measurable quantity. In CHPT, one may treat the θ-term using an external field θx), and the QCD Green functions can be obtained by expanding the generating functional around θx) = θ with real quark masses. Under an axial U) transformation, one has θx) θx) N f α,.3) where N f is the number of flavors, and α = α R α L )/. In the limit of infinitely large number of colors N c, the U) A anomaly is absent. In this case, the spontaneous chiral symmetry breaking of U3) L U3) R, which is a symmetry of the QCD Lagrangian, into U3) V gives nine Goldstone bosons the SU3) flavor octet {π ±, π, K ±, K, K, η 8 } and the flavor singlet η ). Collecting these fields in Ũx), this transforms under the axial U) transformation as Thus, the combination Ũx) e iα R Ũx)e iα L..4) θx) = θx) i ln det Ũx),.5) is invariant under chiral transformations. With θx), the most general chiral effective Lagrangian which is invariant under U3) L U3) R can be constructed. The original construction at order O δ ) can be found in 3], where Oδ) = OM φ, k) with M φ and k denoting the Goldstone boson masses and a small momentum, respectively. Further, we count /N c as O δ ). Here, we adopt the notation used in ], whose formulation is partially based on refs. 33, 34]. The most general chiral effective Lagrangian for mesons to second chiral order, complying with the U3) L U3) R symmetry, reads L = V + V Tr µ Ũ µ Ũ ] + V Tr χ Ũ + χũ ] + iv 3 Tr χ Ũ χũ ] +V 4 Tr Ũ µ Ũ ] Tr Ũ µ Ũ ] + V 5 Tr µ θ µ θ ],.6) where χ = B M q with M q = diagm u, m d, m s ) the real quark mass matrix, µ Ũ = µ Ũ ir µ Ũ + iũl µ, and the V i s are functions of θx). 3

5 In order to use the above Lagrangian, we need to determine the vacuum. Denoting the vacuum expectation value of Ũ by U, we can decompose Ũ as where with φ = Ũ = U U U,.7) ) η U = exp i + i φ..8) 3 F F π π + 6 η 8 π + K + π π + 6 η 8 K K K 6 η 8..9) When θ =, the vacuum is trivial, U =. When the θ angle is switched on, the vacuum is shifted, and U has to be determined by minimizing the zero modes of the Lagrangian. One may parametrize the vacuum as U = diag e iϕu, e iϕ d, e iϕs)..) After the vacuum alignment, one obtains the effective Lagrangian ] ] L φ = V + V Tr µ U µ U + V + BV 3 ) Tr χ +AV 3 Tr U + U ] + V 4 Tr U µ U ] Tr U µ U U + U )] iav Tr U U ] ],.) where χ = B diagm u cos ϕ u, m d cos ϕ d, m s cos ϕ s ), and A, B are complicated functions of the V i s, see e.g. ref. ]. Introducing the notation θ = θ q ϕ q, one may expand the V i s in terms of θ. All of them except for V 3, which is an odd function, are even functions of θ. In the leading approximation, A and B are given by Since A = V ) V ) θ + O δ 4), B = V ) 3 θ V ) + O δ 6)..) ln det U = Tr ln U = i 6 F η, the V i are now functions of θ + 6η /F. The correct normalization of the kinetic terms of the Goldstone boson fields can be obtained by requiring V ) = V ) = F π 4, V 4) = F Fπ )..3) Finally, the quantity θ can be expressed in terms of the measurable quantity θ 7], θ = + 4V ) F π M π 4MK M π M K Mπ)] θ..4) In ref. ], there is one more term iv 3 BV )Trχ U U ) ]. However, this term vanishes exactly because B is defined as V 3/V. 4

6 One sees that θ = O δ ), because V ) is of the zeroth chiral order. Similarly, one can also construct the most general effective Lagrangian in U3) L U3) R CHPT for the baryon octet for a general discussion of effective Lagrangians for the θ-term, see ]) B = Σ + 6 Λ Σ + p Σ Σ + 6 Λ Ξ Ξ n 6 Λ..5) The Lagrangian up to the second chiral order is only the terms relevant to our calculation are displayed; for details, see ref. ]) L φb = i Tr Bγ µ D µ, B] ] m Tr BB] D Tr Bγ µ γ 5 {u µ, B} ] F Tr Bγ µ γ 5 u µ, B] ] + w Tr Bγ µ γ 5 B ] Tru µ ] + b D Tr B { χ+, B} ] + b F Tr B χ+, B] ] + b Tr BB] Tr χ + ] 6 +4A w η Tr F BB] + i w 3 θ 6 ) + w 3 η Tr Bσ µν { γ 5 F + F µν, B }] +i w 4 θ 6 ) + w 4 η Tr Bσ µν γ 5 F + F µν, B ]],.6) where χ + = χ + iau U ), and we use the same notation w = w + 3w / as in ref. 7] with the low-energy constants LECs) w and w defined in ref. ]. We remark here that although there seems to be quite a number of unknown LECs, i.e. w, w, w 3, w 3, w 4 and w 4, only two combinations of these will finally appear in the expressions of the baryon EDFFs. 3 Baryon electric dipole form factors in the infinite volume The electromagnetic form factors of a baryon are defined by Bp ) Jem Bp) ν = ūp ) γ ν F q ) i F q ) σ µν q µ m B +i γ ν q γ 5 m B q ν ) γ 5 FA q ) F 3 q ) ] σ µν q µ γ 5 up), 3.) m B where F q ) and F q ) are the P- and CP-conserving Dirac and Pauli form factors, respectively, while F A q ) is the P-violating anapole form factor, and F 3 q ) if the electric dipole form factor EDFF) which breaks P and CP symmetries. Here, J ν em is the electromagnetic current, the baryon mass is denoted by m B, and q µ = p µ p µ is the four-momentum transfer. In what follows, we will consider the dipole form factor F 3 q ). The electric dipole moment of a baryon is defined as the electric dipole form factor at q = d B = F 3,B) m B. 3.) Note that similar to the case of the neutron electric form factor, we do not include the normalization of the form factor at q = in this definition. 5

7 η a) b) c) d) e) f) Figure. Feynman diagrams contributing to the baryon EDFFs up to next-to-leading order, where denotes a CP violating vertex, black dots represent the second order mesonic vertices, filled squares and diamonds are the first and second order baryonic vertices, respectively. Baryons Tree-level p 4 3 e θ α w 3 + 3w 4 ) + w 3 + 3w 4 ] n 8 3 e θ αw 3 + w 3 ) Λ 4 3 e θ αw 3 + w 3 ) Σ e θ α w 3 + 3w 4 ) + w 3 + 3w 4 ] Σ 4 3 e θ αw 3 + w 3 ) Σ 4 3 e θ α w 3 3w 4 ) + w 3 3w 4 ] Ξ 8 3 e θ αw 3 + w 3 ) Ξ 4 3 e θ α w 3 3w 4 ) + w 3 3w 4 ] Table. Tree-level contribution to the EDFFs of the octet baryons. 3. Baryon electric dipole form factors up to NLO Consider first the tree level contributions to the baryon EDFFs, figure a,b). In particular, these feature the counterterms w 3,4 and w 3,4 at O δ ), 7]. The expressions of the form factor F 3 q )/m) of the ground state octet baryons from the tree-level diagrams are collected in table, where α = 44V ) V ) 3 /F F π M η ). One notices that the tree-level contributions to various baryons depend only on two combinations of the LECs: αw 3 + w 3 and αw 4 + w 4. The dependence on the latter may be written in the general form 4Q B e θ αw 4 + w 4 ), with Q B the baryon electric charge. There are four loop graphs contributing to the baryon EDFFs up to NLO, see figure c-f). Formally, there are more diagrams. However, due to cancellations they only contribute starting at The overall sign of the tree-level expression for the neutron EDM given in ref. 7] should be positive. 6

8 Baryons Loops C cd C ef n {π, p} D + F ) b D + b F ) D + F ) b D + b F ) {K +, Σ } D F ) b D b F ) D F ) b D b F ) p {π, p} D + F ) b D + b F ) {π +, n} D + F ) b D + b F ) {K, Σ + } D F ) b D b F ) {K +, Λ} 3 D + 3F ) b D + 3b F ) {K +, Σ } D F ) b D b F ) {η 8, p} 3 D 3F ) b D 3b F ) {η, p} F π β 3F Table. Possible loops contributing to the EDFFs of the nucleons and the corresponding coefficients C cd and C ef. The intermediate states of the loops are listed in the second column for each nucleon. next-to-next-to-leading order. For details, see refs. 7, 35]. The expression for the sum of loops of a given pair of any meson-baryon intermediate state is F M m 3 q ) m = 8e θ V ) ) ] {C Fπ 4 cd J MM q ) + mm m) + M q J MM m q, m ) } + C cd + C ef ) J M m m ) + O δ 4), 3.3) where m is the mass of the external baryon, M and m are the masses of the meson and baryon in the loops, respectively. We have made use of the fact that the baryon mass difference is O δ ), see appendix A. Keeping only the leading order, the loop expression does not depend on the baryon mass, and reads F M m 3LO q ) m = 8e θ V ) Fπ 4 C cd J MM q ). 3.4) The expressions for the loop functions involved in eqs. 3.3) and 3.4), making use of infrared regularization 36], are collected in appendix B. Different baryons receive contributions from loops with different intermediate states. For instance, the loops for the neutron can be {π, p} and {K +, Σ }, while they are {π +, n}, {π η 8, η ), p}, {K +, Σ Λ)} and {K, Σ + } for the proton. A list of the possible intermediate states and the corresponding coefficients C cd and C ef for the nucleons are given in table, where we have defined β = D 3w ) b D + 3b + 6w ), 3.5) for brevity. The loops and coefficients for the Σ, Λ and Ξ hyperons are collected in table 3. The loops are divergent. Up to O δ 3), the divergences can be absorbed into the renormalization of w 3 7

9 and w 4, w 3 = w r 3µ) + w 4 = w r 4µ) + where the divergence is contained in L = µd 4 4π) 4V ) F 4 π 8V ) 3F 4 π Db F + F b D ) L, 5Db D + 9F b F ) L, 3.6) { d 4 ln4π) + Γ ) + ]}, 3.7) with d the number of space-time dimension, and w 3,4 r µ) are the finite parts of w 3,4. Notice that w 3,4 r µ) depend on the renormalization scale µ through L. The scale dependence cancels with that of the loops, and as a result, the final expressions of the EDFFs are scale-independent. With the coefficients collected in tables and 3, one can easily obtain the explicit expressions of the baryon EDFFs up to NLO. They are listed in appendix C. So far, there is no constraint on the LECs w 3,4 r µ) and w 3,4, which appear at the LO of the baryon EDFFs, except for the N c scaling and the naturalness requirement of the effective field theory. As mentioned in section 3., the tree-level expression of the baryon EDFFs depend only on two combinations w a µ) αw 3 + w 3µ) r 3.8) and αw 4 +w 4 r µ). In principle, these two combinations may be extracted from lattice calculations of the baryon EDMs. On the lattice, the calculations can be performed at different quark masses, or equivalently pion masses. From a fitting to the pion mass dependence of some of the baryon EDMs, one may extract the unknown LECs, and then make a prediction of the other baryon EDMs. Furthermore, up to the order O δ 3), w b µ) 3αw 4 + w 4µ)] r V ) β + 4πF F π M η 3.9) m ave always appear together in the expressions of the baryon EDFFs because of SU3) flavor symmetry. The second term is due to the loops involving the η, and its pion mass dependence starts from O δ 4). Here, all the baryon masses in the second term have been replaced by the average mass of the baryons m ave since the difference is a higher order effect. Consequently, NLO one-loop expressions for all baryon EDFFs are given in terms of just two unknown LECs. 3. Numerical results for the loop contributions of the baryon EDMs Since the values of the LECs w 3,4 r µ) and w 3,4 are not known, we will only focus on the contributions from the loops in this section, and discuss several relations which are free of these parameters in section 3.4. In order to get the numerical results, we use D =.84 and F =.463. From fitting to the baryon mass differences at O δ ), see appendix A, we get b D =.68 GeV and b F =.9 GeV. There are other determinations of b D and b F from higher order analysis of various baryonic properties, see, for instance, refs. 37, 38]. The difference reflects higher order effects, so that they will not be used here. In the large N c -limit, one has F = F π, and 8

10 Baryons Loops C cd C ef Σ + {π, Σ + } 4F b F {π +, Λ} 4 3 Db D {π +, Σ } 4F b F { K, p} D F ) b D b F ) {K +, Ξ } D + F ) b D + b F ) {η 8, Σ + } 4 3 Db D {η, Σ + } F π β 3F Σ {π +, Σ } 4F B F 4F B F {π, Σ + } 4F B F 4F B F {K +, Ξ } D + F ) b D + b F ) D + F ) b D + b F ) {K, p} D F ) b D b F ) D F ) b D b F ) Σ {π, Σ } 4F b F {π 4, Λ} 3 Db D {π, Σ } 4F b F {K, Ξ } D + F ) b D + b F ) {K, n} D F ) b D b F ) {η 8, Σ } 4 {η, Σ } Λ {π +, Σ } 4 3 Db D 3 Db D Fπ β 3F 4 3 Db D {π, Σ + 4 } 3 Db D 4 3 Db D {K +, Ξ } 3 D 3F ) b D 3b F ) 3 D 3F ) b D 3b F ) {K, p} 3 D + 3F ) b D + 3b F ) 3 D + 3F ) b D + 3b F ) Ξ {π +, Ξ } D F ) b D b F ) D F ) b D b F ) {K, Σ + } D + F ) b D + b F ) D + F ) b D + b F ) Ξ {π, Ξ } D F ) b D b F ) {π, Ξ } D F ) b D b F ) { K, Σ } D + F ) b D + b F ) {K, Λ} 3 D 3F ) b D 3b F ) {K, Σ } D + F ) b D + b F ) {η 8, Ξ } 3 D + 3F ) b D + 3b F ) {η, Ξ } β Table 3. Possible loops contributing to the EDFFs of the hyperons and the corresponding coefficients C cd and C ef. The intermediate states of the loops are listed in the second column for each hyperon. F π 3F 9

11 .. loop F 3,n e Θ. loop F 3,n e Θ..3 LO NLO Figure. Loop contribution to the neutron EDFF up to LO left) and NLO right). The solid, dotted, dashed and dot-dashed lines are for the pion mass 38 MeV physical value), MeV, 3 MeV and 4 MeV, respectively.... dn loop e Θ fm..4 d loop eθ fm Π M Π GeV.. loop d e Θ fm.. loop d e Θ fm Π Π Figure 3. Loop contributions to the EDMs of the neutral baryons as a function of the pion mass. The bands, reflecting uncertainties by varying the renormalization scale between µ = M ρ and m Ξ, between solid and dashed boundaries are the NLO and the LO results, respectively. The filled circle and square with error bars are the lattice data from refs. ] and 6], respectively. we take 9. MeV 39] for its value. From an analysis of the η η mixing in the framework of U3) L U3) R CHPT, it was found that V ) = 5 4 GeV 4 and V ) 3 = GeV 4]. With these values, the neutron EDFF calculated at µ = GeV is plotted in figure, where only

12 the loop contribution is taken into account. One sees sizable NLO effects, especially for higher pion masses. One notices that the dependence on q can be well approximated by a linear function, which means that a linear extrapolation from finite to vanishing q can be used on the lattice. The same is true for the other baryons, as none of them shows a strong q -dependence. At the physical pion mass, we get the loop contributions to the baryon EDMs in units of 6 e θ cm, d loop n = 3. ±.8, d loop p = 5.6 ±. 6. β GeV), d loop Λ =.6 ±.4, d loop = 3.8 ±. 4.8 β GeV), Σ + d loop =.8 ±.4, d loop =. ± β GeV), Σ Σ d loop = 3.6 ±.8, d loop = 3.7 ± β GeV), 3.) Ξ Ξ where the uncertainties are estimated by varying the scale µ between the masses of the ρ and Ξ. If we replace the baryon masses in the β-term in the baryon EDFF expressions, eqs. C.), C.4), C.5) and C.8), by the averaged baryon mass m ave = 5 MeV as that in eq. 3.9), this term contributes 5. β GeV) to the proton and Σ + and 5. β GeV) to the Σ and Ξ. The difference reflects part of the higher order uncertainties. In the following, we will use the averaged baryon mass for the β-term and keep different masses for the other terms. In order to compare with the results from lattice simulations, we should study the pion mass dependence of the pertinent quantities. We take the physical values for F π, M η8,η and the baryon masses, since their M π -dependent effects contribute at higher orders. For the kaon mass, we use M K = M K + M π, 3.) where M K = 484 MeV is the kaon mass in the SU) chiral limit with vanishing up and down quark masses. The pion mass dependence of the EDMs of the neutral baryons is shown in figure 3, where the shaded bands between the solid and dashed boundaries are the NLO and LO results, respectively. It is clear that all the EDMs vanishes in the chiral limit. The physical reason is that the θ-term can be rotated away if any of the quarks are massless, see, for instance, ref. 3], which ensures a vanishing value for the θ-term induced EDM. For comparison, the lattice data for the neutron at M π = 53 MeV calculated in refs. ] are shown, which are calculated on lattice with lattice spacing a. fm. We also show the lattice data reported very recently in ref. 6] whose uncertainty is smaller. Both the LO and NLO loop results agree with the lattice data. One notices that the NLO corrections for all the neutral hyperons are dramatic. For the Λ, from table and eq. 3.3), one finds that the loops involving π + and Σ cancel with those of π and Σ + exactly. Thus, the remaining LO contributions come from loops involving a kaon. They are small as can be seen from + ln M π µ =.96, + ln M K µ =.4, with µ = GeV. Therefore, the LO result for the Λ EDM is close to zero as already pointed out in ref. ]). The same happens for the Σ.

13 dp loop eθ fm M Π GeV loop d eθ fm M Π GeV loop d eθ fm M Π GeV loop d eθ fm M Π GeV Figure 4. Loop contributions to the EDMs of the charged baryons as a function of the pion mass. The bands, reflecting uncertainties of the LO contributions by varying the renormalization scale between µ = m ρ and m Ξ. The filled circle and square with error bars are the lattice data from refs. ] and 6], respectively. The solid, dotted and dot-dashed lines represent the NLO results evaluated at µ = GeV with β =, and GeV, respectively. The results for the charged baryons are shown in figure 4, where the bands are the LO loop results. We choose three different values for the unknown combination β to illustrate the NLO effects. The results are quite sensitive to the numerical value of β, especially for the Ξ and Σ. 3.3 A first determination of the LECs from lattice results We may use the neutron and the proton EDMs at unphysical quark masses to determine the two combination of parameters w a µ) and w b µ). Using the lattice data at M π = 53 MeV 6] 3 for the neutron EDM, one may determine the counterterm combination w a µ) w a GeV) =. ±.) GeV, 3.) where the uncertainty merely reflects the uncertainty in the lattice calculations. Using the data for the proton, we get w b GeV) =.4 ±.5) GeV. 3.3) 3 As this pion mass value is at the edge of the range of applicability of our approach, lattice data at lower pion mass are urgently called for for a more reliable determination of the LECs.

14 Baryon p n Σ + Σ Σ Λ Ξ Ξ Combination w a + w b ) w a w a + w b ) w a w b w a w a w a w b w a Table 4. Combinations of the unknown parameters appearing in the EDFFs of the baryons up to NLO. With these determinations, the EDMs at the physical pion mass can be predicted including both the tree and the loop contributions. The ones for the neutron and proton in units of 6 e θ cm are d n =.9 ±.4 ±.8, d p =. ±.5 ±.. 3.4) Here, the first uncertainty reflects the uncertainty in the determination of w a GeV) and w b GeV), and the second one corresponds to varying the scale µ between the ρ-meson mass and m Ξ. Comparing with the loop results given in eq. 3.), one sees that the neutron EDM at the physical pion mass is dominated by the loops. For the proton, the loop contributions have similar size as, albeit slightly larger than, the tree level terms. Combining these two errors, we have d n =.9 ±.9 and d p = d p =. ±. in canonical units). Using the experimental upper limit of the neutron EDM, d n <.9 6 e cm ], the θ angle is constrained to be θ ) Similarly, we can predict the hyperon EDMs again in units of 6 e θ cm) 3.4 Counterterm-free relations d Λ =.5 ±. ±.4, d Σ + =.7 ±.5 ±., d Σ =.7 ±. ±.4, d Σ =. ±.5 ±. d Ξ = 3.4 ±.4 ±.8, d Ξ =.6 ±.5 ±.. 3.6) In this section, we will derive relations that are free of the unknown LECs up to NLO. These relations can e.g. serve as checks of lattice simulations for varying quark masses. In fact, due the SU3) flavor symmetry, there exists an exact relation among the EDFFs of the Σ hyperons. From table, one deduces that the loops involving a pion or an η 8) for the Σ + cancel with those of the Σ. As mentioned before, the pionic loops of the Σ cancel with each other. Furthermore, there are also cancellations for the kaonic loops among the Σ hyperons. These cancellations occur when isospin is a good symmetry. Therefore, in the isospin symmetric case, which is considered throughout, one has F 3,Σ + + F 3,Σ F 3,Σ = O δ 4). 3.7) This relation is exact up to NLO including the tree-level contributions. Similarly, there is another relation which is exact only at LO F 3,Σ + F 3,Λ = O δ 3), 3.8) since the LO loop contributions are independent of the baryon mass, and they cancel each other in the sum. Moreover, the above combination is also independent of β which combines the unknown 3

15 d d 6 e Θ cm dn d 6 e Θ cm d d 6 e Θ cm Π Π Figure 5. Counterterm-free combinations of the baryon EDMs as a function of the pion mass. parameters w and w, so that it can be calculated parameter-free. Taking q =, the O δ 3) expression is rather simple ) 6eV θ d Σ + d Λ = M 3πFπ 4 M K M ) π F b D + Db D b F + 3F b ) F + O δ 4 ). 3.9) K In fact, there are more combinations free of the counterterms. As mentioned before, up to NLO, there are only two combinations of unknown constants, w a µ) and w b µ), in the expressions of the EDFFs of the octet baryons defined in eqs. 3.8) and 3.9). It is easy to find combinations free of any unknown parameters up to O δ 3) utilizing the results collected in table 4. In addition to the two in eqs. 3.7) and 3.9), they include F 3,p F 3,Σ +, F 3,n F 3,Ξ, F 3,Σ F 3,Ξ. 3.) Because the counterterms have been cancelled out, these combinations are also finite and hence independent of the scale µ. However, different from eq. 3.9), their LO loop contributions do not vanish. The expression for the F 3,n F 3,Ξ at q = reads d n d Ξ = ev ) θ Db D + F b F ) ln M K + π M ) π M K π Fπ 4 Mπ + 8π M M K M ) π Db D + F b D b F + Db F K m ave ) ] + O δ 4). 3.) 4

16 The relations for the charged baryons are more complicated, and are not shown here they can be obtained from the expressions given in appendix C). The pion mass dependence of the four combinations in eqs. 3.8) and 3.) are shown in figure 5. Their values at the physical pion mass in units of 6 e θ cm are d Σ + d Λ =.8, d p d Σ + =.8, d n d Ξ =.5, d Σ d Ξ =.6. 3.) Note that although the first combination starts from one order higher than the others, their numerical values are of similar order of magnitude. 4 Finite-volume corrections in the p-regime On the lattice, calculations are performed in a finite volume. As a result, the continuum momentum spectrum becomes quantized. Taking periodic boundary condition for all three spatial dimensions forms a torus. If the volume is L 3, the momentum takes values of π n/l, with n a three-dimensional vector of integers. One easily sees that any integral over the spatial components of momentum in the infinite volume should be replaced by a summation for the momentum modes. For instance, the two-point scalar loop integral needs to be changed as follows, i d 4 k π) 4 k m )k + q) m ] i dk L 3 π n k m )k + q) m 4.) ]. The finite volume corrections to a quantity Q is defined as the difference of Q evaluated in a finite and infinite volume, δ L Q] = QL) Q ). 4.) As we already mentioned, at chiral limit the baryon EDMs should vanish. However, the current lattice data at rather large pion masses do not show a decreasing behavior yet, ]. Thus, results at smaller pion masses are necessary. However, in the section, we will show that the finite volume corrections at small pion masses to most of the baryon EDMs are quite large, and it is necessary to include the NLO contributions. The decomposition of the matrix element of the electromagnetic current in the form of eq. 3.) assumes Lorentz invariance. However, on discretized lattice, Lorentz invariance in infinite volume is lost. Furthermore, the external momentum is also quantized in a finite volume. These effects should be taken into account when calculating the finite volume corrections to the form factors 4, 4]. Since we are not aiming at a very precise calculation of the finite volume effects, we will simply assume Lorentz invariance in the following. As shown in ref. 43], one can use the same LECs as in the infinite volume, and the finite volume corrections come only from loops. This may be understood easily, since the finite volume effects are long-distance physics, while the LECs reflect the short-distance physics. For the same reason, the finite volume corrections should not depend on the choice of method to regularize the ultraviolet divergence. This means that one may change the upper bound of the Feynman parameter integration in eqs. D.5) and D.7) in the infrared regularization from to as already noticed in, for instance, ref. 44]. In the p-regime, where M π L 43, 45], the finite volume corrections are exponentially suppressed. 5

17 . L dn dn loop... 4 M Π 38 MeV M Π MeV M Π 3 MeV M Π 4 MeV L d n NLO L d n LO L fm L fm Figure 6. Left: The ratios of the finite volume corrections to the loop contributions in the infinite volume for the neutron EDM as a function of L. Right: The ratios of the finite volume corrections at NLO LO not included) to those at LO. The solid, dotted, dashed and dot-dashed lines are for the pion mass being physical, MeV, 3 MeV and 4 MeV, respectively. In the left panel, for each value of the pion mass, two lines are plotted with the upper and lower ones representing the NLO and LO results, respectively. In ref. 9], finite volume correction to the neutron EDM was calculated at LO in SU) heavy baryon CHPT. The resulting expression is given by 9] δ L d LO n ] = g A ᾱeθ π F π m u m d K LM π n ), 4.3) m u + m d where ᾱ is a coefficient of the CP-violating NNπ vertex. 4 From eqs. C.) and D.4), our result for the finite volume correction to the neutron EDM at LO is δ L d LO n ] = M πeθ 4π F π n D+F ) bd + b F ) K LM π n ) D F ) b D b F ) K LM K n ) ], n where we have used ] 4.4) θ F π Mπ 8V ) θ, 4.5) which may be obtained from eq..4) considering M π M K. Comparing the SU) part with eq. 4.3), we can identify ᾱ = B b D + b F ). The left panel of figure 6 shows the ratios of the finite volume corrections to the loop contributions in the infinite volume for the neutron EDM as a function of L. The results are obtained with µ = GeV. In order to show the impact of the NLO corrections more clearly, we plot the ratios of the finite volume corrections at NLO LO not included) to the LO in the right panel of figure 6. One sees that at the physical pion mass the NLO contributions increase the LO finite volume corrections by about %. This correction is much larger for a pion mass of 4 MeV. At this point, one might worry about the convergence. In fact, the neutron EDM in the infinite volume 4ᾱ is the α used in ref. 9]. 6

18 . 5 L d d loop.. 4 L d L d LO L fm L fm. loop L d d.. 4 LO L d L d L fm L fm.5 loop L d d... LO L d L d L fm L fm Figure 7. Left: The ratios of the finite volume corrections to the loop contributions in the infinite volume for the EDMs of neutral hyperons as a function of L. Right: The ratios of the finite volume corrections of NLO to those of LO. The solid, dotted, dashed and dot-dashed lines are for the pion mass being 38 Nev physical value), MeV, 3 MeV and 4 MeV, respectively. has a nice convergence property, as can be seen from the left panel of figure 3. Moreover, one can show that the ratio tends to increase with increasing L. Although both the LO and NLO finite volume corrections decrease exponentially, the LO one decreases faster. Using eqs. E.) and E.4) from appendix E, we get in the limit L, ] δ L d NLO n δ L d LO n ] π M π LMπ exp m N LM 3 π 8m N ). 4.6) The finite volume corrections to the neutral hyperons are shown in figure 7. It is obvious that the 7

19 Β GeV LO loop L dp dp L dp dp loop... Β GeV L dp dp loop... Β GeV L fm Figure 8. The ratio of the finite volume corrections to the loop contributions in the infinite volume for the proton EDM as a function of L. The first graph presents the LO results, and the others are the NLO results for different choices of β. The solid, dotted, dashed and dot-dashed lines are for the pion mass of 38 MeV, MeV, 3 MeV and 4 MeV, respectively. effects are very small for the Λ and Σ. The reason is that the pionic loops with positively and negatively charged pions cancel each other exactly for these two baryons. This also makes the NLO corrections much larger than the LO ones at µ = GeV, consistent with the infinite volume results shown in figure 3. The situation for the Ξ is similar to case of the neutron. Although the parameter combination w b µ) can be determined from the lattice data of the proton EDM, β is still unknown since it always appear together with the LECs w 4 and w 4 rµ). However, the finite volume corrections for the charged baryons depend on β, but not on w 4 and µ). It is therefore possible to extract the value of β from the volume effects of the baryon EDMs. In figure 8, we show the ratios of the finite volume corrections to the loop contributions in infinite volume to the proton EDM for different chosen values of β. At this point, we want to emphasize that the EDMs of all the charged baryons depend on the same combination of LECs, i.e. β. Once it is determined from one baryon, it can be used for predicting the finite volume corrections to the others. w r 4 5 Summary In this paper, we extended the previous calculation 7] of the nucleon EDFFs and EDMs to the ground state baryon octet up to NLO in the framework of U3) L U3) R chiral perturbation theory. 8

20 Our main findings can be summarized as follows: ) We have shown that the complete one-loop expressions for the baryon EDFFs and EDMs depend on two combinations of unknown LECs only, w a and w b, cf. eqs.3.8,3.9). In case of the charged baryons, the combination w b µ) combines two LECs from the tree graphs and one LEC that appears only in loops. ) We have shown that the NLO corrections are large for the neutral hyperons Λ and Σ. This is due to a suppression of the LO contributions based on exact cancellations between loops of positively and negatively charged pions. For the charged baryons, we find a strong sensitivity to the LEC combination w b. 3) We have derived a set of relations between various EDMs that are free of unknown LECs. These can be useful for future lattice simulations of baryon EDMs. 4) Based on recent lattice results for the neutron and proton EDMs at M π = 53 MeV, we could pin down the two LEC combinations w a and w b. Based on this, we can predict the baryon EDMs at the physical pion mass. In particular, we find d n =.9 ±.9 and d p =. ±. in units of 6 θ e cm. 5) The finite volume corrections to the baryon EDMs in the p-regime are also studied. Because the loops contribute from LO, it is found that the finite volume corrections are huge for all the baryon EDMs except for the Λ and Σ for the pion mass close to its physical value. For the neutron, the finite volume correction is about % at M π L = 4. Note that the calculation of the finite volume corrections in our paper assumes Lorentz invariance, and additional subtleties in the presence of external electromagnetic fields, see 46], are not taken into account. For a precise calculation of the finite volume corrections of the baryon EDMs, separate analyses need to be done for the lattice calculations using the form factor method, the spectrum method with external field and methods with the twisted boundary conditions. More lattice results of the nucleon EDMs are expected to come out soon. Calculations of the hyperons are welcome based on our analysis. With these upcoming lattice data, one can determine the LECs more reliably, and a next-to-next-to-leading order calculation would be feasible and desirable in view of the bad convergence for some baryons. Furthermore, it would be interesting to see how the LEC-free relations get modified at higher orders. Acknowledgments We are grateful to E. Shintani for providing us with lattice data before publication. We acknowledge discussions with T. Izubichi and G. Schierholz. We also thank B. C. Tiburzi for very useful comments on the finite volume corrections. This work is supported in part by the DFG and the NSFC through funds provided to the Sino-German CRC Symmetries and the Emergence of Structure in QCD and the EU I3HP Study of Strongly Interacting Matter under the Seventh Framework Program of the EU. U.-G. M. also thanks the BMBF for support Grant No. 6BN78). F.-K. G. acknowledges partial support from the NSFC Grant No. 655). 9

21 Baryons N Σ Λ Ξ b π b F b D b D /3 b F b K b D b F 4b D /3 b D + b F Table 5. Coefficients b π,k in the NLO baryon mass formula. A Baryon masses up to NLO From eq..6), it is easy to obtain the baryon masses up to O δ ), which may be written as m B = m b M π + M K) 4 bπ M π + b K M K) + O δ 3 ), A.) where we have expressed the quark masses in terms of the meson masses at the lowest order, and b π,k are listed in table 5. B Loop integrals in infrared regularization Loop integrals in the infinite volume used in this paper are collected in this appendix. Let us denote the masses of the light meson and baryon in the loop by M and m, respectively. The momentum of the external current is q µ, and the momenta of the external baryons are denoted by p µ and p µ = p µ + q µ, respectively. We further define L = µd 4 4π) { d 4 ln4π) + Γ ) + ]}, which contains the divergence at the spacetime dimension d = 4. The -point mesonic-loop function when the two mesons have the same mass reads as J MM q d d k ) = i π) d k M + i ɛ) k + q) M + iɛ] = L + ln M 6π µ σ ln σ ), B.) σ + where σ = 4M /q. When there is a baryon in the loop, the loop evaluated in the normal dimensional regularization spoils power counting 47]. The power counting can be restored in the heavy baryon formalism of baryon CHPT 48, 49]. There are also covariant formalisms, see ref. 5] for a recent review. Here, the infrared regularization 36, 5] will be used. In the infrared regularization, the integral is separated into an infrared regular piece and an infrared singular piece. The infrared regular piece can be expanded analytically in the chiral expansion, and absorbed into the counterterms. Only the infrared singular part of the loop integrals will be shown in the following. The two-point scalar loop integral involving one meson and one baryon is J Mm p d d k ) = i π) d k M + i ɛ) p k) m + iɛ] = ω p L + 3π p ω ln M ) µ + λp, m, M ) arccos ω )] B.), pm

22 where ω = p m + M, p = p, and λx, y, z) = x + y + z xy + yz + xz) is the Källén function. Noting that the external baryons are on shell, p is given by the mass squared of the external baryon. Thus, with the baryon mass expressions in appendix A, one has p = m +O δ ) and ω = O δ ). Expanding the arccos function, we get J Mm p ) = M 6π p + ω p L + 3π + ln M µ )] + O δ 3), B.3) where the first and second term are of Oδ) and O δ ), respectively. The three-point loop integral with two mesons and one baryon can be written as J MMm q, p d d k ) = i π) d k M + i ɛ) k + q) M + iɛ] p k) m + i ɛ] = M dxj Mm p ), B.4) where M = M + xx )q, p µ = p µ + x q µ. In the following, p = p will always be assumed, then one has p = p + xx )q. The leading chiral order of this integral is Oδ ). Keeping terms up to Oδ ), the analytic expression for the infrared singular part is J MMm q, p q ) = 6π p q arctan M L + p + 3π p + ln M ω q + µ q ln σ σ σ + )] + Oδ), B.5) where the first term is of order Oδ ), and the terms in the square brackets are of order Oδ ). For completeness, we also give the expression for the loop with one meson and two baryons, though not used in our calculations. Similar to J MMm q, p ), this loop can also be worked out from the two-point loop with one meson and one baryon, J Mmm p, p d d k ) = i π) d k M + i ɛ) p k) m + i ɛ] p k) m + i ɛ] = m dxj M m p ), where m = m + xx )q. The analytic expression up to NLO is L J Mmm p, p ) = p + 3π p + ln M )] ω µ + 64πp mm + O δ ), B.6) B.7) where the terms in the square brackets are at LO, and the second term is at NLO. Up to NLO, this loop does not depend on the momentum transfer q. C Expressions for the baryon EDFFs up to NLO Up to NLO, which is O δ 3), the baryon EDFFs contain two parts: the tree-level expressions are given in table, and the loop contributions are obtained by taking the LO terms of the baryonic loops J MM m q, m ) and J M m q ). The explicit expressions of the EDFFs of the baryon octet up to NLO are given in the following:

23 For the neutron, F 3,n q ) m N = 8 3 e θ αw3 + w 3µ) r ] { + V ) e θ π Fπ 4 D + F ) b D + b F ) ln M π µ + σ π ln σ π σ π + + π Mπ q ) ] q arctan D F ) b D b F ) ln M K m N q M π µ + σ K ln σ K σ K + + π M K q 8 b D b F ) M q K M ) ) ]} q π arctan, C.) m N M K with α = 44V ) V ) 3 /F F π M η ), and σ πk) = 4MπK) /q. Here, we have used the SU3) mass splitting for the baryon masses m Σ m N = 4 b D b F ) MK M π) + O δ 3 ). Similar relations derived from the Lagrangian eq..6) will be used in the following. For the proton, F 3,p q ) m N = 4 3 e θ α w3 + 3w 4 ) + w 3µ) r + 3w 4µ) r ] { V ) e θ 6π Fπ 4 6D + F ) b D + b F ) ln M π µ + σ π ln σ π σ π + + 3πM π + π M π q ) ] q arctan m N m N q M π +4 Db D + 3F b F ) ln M K µ + σ K ln σ K σ K + + πm ) K m N + 4π q arctan DbD + 3F b F ) M q M K m K q ) N +8 MK Mπ) F ) ] b D + 3b ) F 3 Db D b D 3b F ) + π 6D F ) b D b F ) M K + D 3F ) b D 3b F ) M η8 + F ] } π m N F βm η, C.) where β = D 3w ) b D + 3b + 6w ). For the Σ, F 3,Σ q ) m Σ = 4 3 e θ αw3 + w 3µ) r ] { V ) e θ π Fπ 4 +σ K ln σ ) K + π σ K + q arctan q M K Db F + F b D ) ln M K µ DbF + F b D m Σ M K q ) +8 M K M π) F b D + Db D b F + F b F ) ]}. C.3)

24 For the Σ +, F 3,Σ +q ) m Σ = 4 3 e θ α w3 + 3w 4 ) + w 3µ) r + 3w 4µ) r ] { V ) e θ 3π Fπ 4 Db D + 3F b F ) ln M π µ + σ π ln σ π σ π + + πm ) π m Σ +3D + F ) b D + b F ) ln M K µ + σ K ln σ K σ K + + πm K m Σ + π M K q + 8 b D + b F ) M q K Mπ m Σ ) ) ] q arctan M K DbD + 3F b F ) M m π q ) + 3 Σ 3 Db D M K Mπ ) ] + π q arctan q M π + π 6F b F M π + 3D F ) b D b F ) M K + Db D M η8 + F ] } π m Σ F βm η. C.4) For the Σ, F 3,Σ q ) m Σ = 4 3 e θ α w3 3w 4 ) + w 3µ) r 3w 4µ) r ] { + V ) e θ 3π Fπ 4 Db D + 3F b F ) ln M π µ + σ π ln σ π σ π + + πm π m Σ + π q arctan DbD + 3F b F ) q M π +3D F ) b D b F ) m Σ M π q ) Db D ln M K µ + σ K ln σ K σ K + + πm K m Σ ) M K Mπ ) ] + π M K q + 8 b D b F ) M q K M ) ) q π arctan m Σ M K + π 6F b F M π + 3D + F ) b D + b F ) M K + Db D M η8 + F ] } π m Σ F βm η. C.5) ] For the Λ F 3,Λ q ) = 4 m Λ 3 e θ αw3 + w 3µ) r ] { + V ) e θ π Fπ 4 + π q arctan q M K 8 3 Db F + F b D ) ln M K µ + σ K ln σ ) K σ K + DbF + F b D m Σ M K q ) M K Mπ ) F b D + Db D b F + F b ) ]} F. C.6) 3

25 For the Ξ, F 3,Ξ q ) = 8 m Ξ 3 e θ αw3 + w 3µ) r ] { V ) e θ π Fπ 4 D F ) b D b F ) ln M π µ + σ π ln σ π σ π + + π Mπ q ) ] q arctan D + F ) b D + b F ) ln M K m Ξ q M π µ + σ K ln σ K σ K + + π M K q 8 b D + b F ) M q K M ) ) ]} q π arctan, C.7) m Ξ M K For the Ξ, F 3,Ξ q ) = 4 m Ξ 3 e θ α w3 3w 4 ) + w 3µ) r 3w 4µ) r ] { + V ) e θ 6π Fπ 4 6D F ) b D b F ) ln M π µ + σ π ln σ π σ π + + 3πM π + π M π q ) ] q arctan m Ξ m Ξ q M π +4 Db D + 3F b F ) ln M K µ + σ K ln σ K σ K + + πm ) K m Ξ + 4π q arctan DbD + 3F b F M q M K m K q ) Ξ 8 MK Mπ) F ) ] b D + 3b ) F + 3 Db D b D + 3b F ) + π 6D + F ) b D + b F ) M K + D + 3F ) b D + 3b F ) M η8 + F ] } π m Ξ F βm η. C.8) The expressions for the EDMs can be easily obtained by noticing lim σ ln σ q σ + =, lim q arctan q q M = M. D Finite volume corrections to loops This appendix is dedicated to finite volume corrections to the two- and three-point loop integrals. Let us consider the scalar two-point loop J m m q d 4 k ) = i π) 4 k m + i ɛ) k + q) m + iɛ] d 4 k = dx i π) 4 k m + iɛ), D.) where m = xx )q + x)m + x m. If the loop involves both a meson and a baryon, m should be replaced by the mesonic mass, and the upper bound of the integration over the Feynman 4

26 parameter x should be replaced by in infrared regularization. Performing the contour integral over k, one gets J m m q ) = d 3 k dx 4 π) 3 3/. D.) k + m ) Using the formula L 3 n ) d 3 k π) 3 ) j = j a 3 j K 3/ j L a n ) k + a π 3/ Γj) L a n ) 3/ j D.3) n derived in ref. 5], we get the finite volume correction to the loop J m m q ), δ L J m m q )] = 8π dx n K L m n ), D.4) with K ν z) a modified Bessel function of the second kind. From eq. B.4), it is easy to write down the correction to the three-point loop integral, δ L J MMm q, p )] = M = L dx δ L J Mm p )] 6π n dx n dy y z K z ), D.5) where z = L n yy ) p + ym + y) M ] /. We have made use of the integral representation of K n z) 53] K n z) = Γn + /)z)n π which is valid for Re n > /, arg z < π/. Similarly, one gets δ L J Mmm p, p )] = L 6π n dx n where z = L n yy ) p + y m + y)m ] /. cos t dt. D.6) t + z n+/ ) dy y z K z ), D.7) E Asymptotic expansion of finite volume corrections For simplicity, we will focus on the asymptotic expansion of finite volume corrections with q = and p = m. The general case can be treated similarly. For z, one has )] π K n z) z e z + O. E.) z For L, we can consider the term with n = only in eq. D.4). There are 6 possibilities, so that for L, δ L J MM )] 3 π e LM 4π. E.) LM 5

27 For the three-point loop integral, δ L J MMm, m )] L π 6π dy y) e LCy) LCy)] 3/, E.3) with Cy) = y m + y)m ] /. The leading term of the above integral over the Feynman parameter y can be worked out using Laplace s method. For very large L, the integral receives contributions mostly from the neighborhood of the minimum of Cy), which is at y c = M /m ). Thus, δ L J MMm, m )] 3L π 8π y c) e LCyc) LCy c )] 3/ = 3 8πmM exp LM M 4m The asymptotic expansion for δ L J Mmm m, m )] is the same. dy exp ] y C y c ) ) M + O m )]. E.4) References ] C. A. Baker et al., An improved experimental limit on the electric dipole moment of the neutron, Phys. Rev. Lett. 97 6) 38 hep-ex/6]. ] M. Pospelov and A. Ritz, Electric dipole moments as probes of new physics, Annals Phys. 38 5) 9 hep-ph/543]. 3] S. K. Lamoreaux and R. Golub, Experimental searches for the neutron electric dipole moment, J. Phys. G 36 9) 4. 4] F. J. M. Farley, K. Jungmann, J. P. Miller, W. M. Morse, Y. F. Orlov, B. L. Roberts, Y. K. Semertzidis and A. Silenko et al., A New method of measuring electric dipole moments in storage rings, Phys. Rev. Lett. 93 4) 5 hep-ex/376]. 5] Y. K. Semertzidis Storage Ring EDM Collaboration], A Storage Ring proton Electric Dipole Moment experiment: most sensitive experiment to CP-violation beyond the Standard Model, arxiv:.3378 physics.acc-ph]. 6] F. Rathmann and N. Nikolaev, Precursor experiments to search for permanent electric dipole moments EDMs) of protons and deuterons at COSY, PoS STORI ) 9. 7] A. Lehrach, B. Lorentz, W. Morse, N. Nikolaev and F. Rathmann, Precursor experiments to search for permanent electric dipole moments EDMs) of protons and deuterons at COSY, arxiv:.5773 hep-ex]. 8] S. Aoki and A. Gocksch, The neutron electric dipole moment in lattice QCD, Phys. Rev. Lett ) 5 Erratum ibid 65 99) 7]. 9] S. Aoki, A. Gocksch, A. V. Manohar and S. R. Sharpe, Calculating the neutron electric dipole moment on the lattice, Phys. Rev. Lett ) 9. ] E. Shintani, S. Aoki, N. Ishizuka, K. Kanaya, Y. Kikukawa, Y. Kuramashi, M. Okawa and A. Ukawa et al., Neutron electric dipole moment with external electric field method in lattice QCD, Phys. Rev. D 75 7) 3457 hep-lat/63]. 6