arxiv: v1 [hep-lat] 22 Dec 2017 Neutron Electric Dipole Moment on the Lattice

Transcription

1 arxiv: v1 [hep-lat] 22 Dec 217 Neutron Electric Dipole Moment on the Lattice Boram Yoon 1,, Tanmoy Bhattacharya 1, an Rajan Gupta 1 1 Los Alamos National Laboratory, Los Alamos, NM 87545, USA 1 Introuction Abstract. For the neutron to have an electric ipole moment (EDM), the theory of nature must have T, or equivalently CP, violation. Neutron EDM is a very goo probe of novel CP violation in beyon the stanar moel physics. To leverage the connection between measure neutron EDM an novel mechanism of CP violation, one requires the calculation of matrix elements for CP violating operators, for which lattice QCD provies a first principle metho. In this paper, we review the status of recent lattice QCD calculations of the contributions of the QCD Θ-term, the quark EDM term, an the quark chromo-edm term to the neutron EDM. Electric ipole moment (EDM) of the neutron measures the separation of positive an negative charge in the neutron an is necessarily aligne along the spin axis. For a neutron to have an EDM, the theory of elementary particles must violate parity (P) an time reversal (T) invariance, an charge conjugation an parity (CP) symmetry if CPT is conserve. Uner a parity transformation, the EDM of any particle changes its sign while the irection of the spin remains the same. Uner time-reversal transformation, the spin flips its irection while the EDM stays unchange. So far, a nonvanishing neutron EDM (nedm) has not been observe, an the current experimental upper boun is e cm [1]. In recent experiments, the nedm is measure by the change in the spin precession frequency of ultracol-neutrons aligne in a magnetic fiel uner a flip in the irection of a strong backgroun electric fiel. Figure 1 shows the history of the experimental upper boun on the neutron EDM, as well as the target precision, e cm [2 6], of propose experiments. In the stanar moel, the leaing contribution to nedm ue to the CP violating phase in the CKM matrix comes from three-loop an higher orer iagrams, an the expecte size is more than 5 orers of magnitue below the current experimental boun [7]. In extensions of the stanar moels, however, nedm can appear at one-loop ue to novel CP violating interactions. In some of the most popular moels, such as supersymmetric (SUSY) moels, the expecte size of the nedm is between e cm as shown in Figure 1. Therefore, a nedm smaller that 1 28 e cm will put a serious constraint on those moels [8 1]. There are two main outcomes of neutron EDM experiments for BSM physics. First, a non-zero measurement of nedm woul establish new sources of CP violation, which is one of Sakharov s three conitions for weak-scale baryogenesis [11]. The CP violation in the stanar moel is not sufficient to explain observe the baryon asymmetry of the universe, an it requires new sources of Speaker, boram@lanl.gov

2 Neutron EDM Upper Limit (e cm) Supersymmetry Preictions Stanar Moel Preictions Previous Experiments Future Experiments Year of Publication Figure 1. The timeline of the upper boun on the neutron EDM from previous an future experiments. The yellow squares are the previous bouns, an the blue ot is the sensitivity target of next generation experiments. In the stanar moel, the leaing contribution to nedm from the CP violating phase in the CKM matrix arises at three-loops, an the expecte size of nedm is more than 5 orers of magnitue below the current experimental boun, as shown by the green shae area in the plot [7]. In extensions of the stanar moels, however, nedm can appear at one-loop with the new CP violating interactions. In some of the most popular moels, such as SUSY, the expecte size of nedm covers the regions of the current an future experimental boun, as shown by the re shae area [8 1]. CP violation from beyon the stanar moel (BSM) [12]. Neutron EDM is a goo probe for such novel CP violation. Secon, most BSM theories have aitional sources of CP violation. There are many BSM scenarios preicting a nedm between e cm, an upcoming experiments will put constraints on them [2 6]. This requires both the measurement of, or a boun on, the neutron EDM, an the calculations of the matrix elements of novel CP violating operators within the neutron states. Current non-lattice estimates of the matrix elements have large uncertainties [1]. Lattice QCD can provie a first principle, non-perturbative metho for calculating these matrix elements. To analyze novel CP violation, we work within the framework of effective fiel theories an classify operators by their canonical imension [13]. At haronic scale, the effective Lagrangian for the CP violating interactions at imension 4, 5 an 6 can be written as L =4,5,6 CPV = g2 32π 2 θg G i 2 q=u,,s + w g 6 f abc G a µν G νρ,b G µ,c ρ + q q(σ µν F µν )γ 5 q i 2 i q=u,,s q q(σ µν G µν )γ 5 q C (4q) i O (4q) i, (1) where g is the QCD coupling constant, an G µν,b = ε µναβ G b αβ /2. The first term on the right han sie (r.h.s) is the imension-4 Θ-term alreay allowe in QCD, an the next two terms are the imension-5 quark EDM (qedm) an quark chromo-edm (cedm) terms. There are two types of imension-6 terms in the secon line in the above equation: the Weinberg three-gluon operator an various fourquark operators.

3 If the matrix elements of all the operators are O(1) an the anomalous imensions are small, then the O(1/Λ D 4 scale ) suppression implies that the lower mass imension operators are more important. However, it turns out that the current boun on the nedm alreay makes the imension-4 QCD Θ- term unnaturally tiny; θ O( ) [13 16], where θ is the total effective CP violating angle θ = θ + arg et M CKM with the quark mixing matrix M CKM. The unexpecte smallness of this number is known as the strong CP problem. One of the popular moels explaining the strong CP problem is the Peccei-Quinn mechanism [17, 18], which requires the existence of the, as yet unobserve, axion. The two imension-5 operators typically arise at the TeV scale as imension-6 operators involving a Higgs. At the haronic scale, 1 GeV, the Higgs fiels are replace by its vacuum expectation value, v EM, an the operator becomes imension-5 with the coefficients suppresse by v EM /MBSM 2. How well these imensional arguments work an which terms in Eq. (1) make the largest contribution to nedm epens on the etails of the BSM theory. But, while the couplings epen on the starting BSM moel, the matrix elements of these operators within the neutron states are moel inepenent results. Thus, lattice QCD calculation of the matrix elements of these operators will play an important role in connecting the measure neutron EDM an novel CP violation in BSM scenarios, i.e., knowing the matrix elements along with the boun on (or the value of) the nedm, one can put bouns on the couplings an thus on the parameter space of allowe BSM theories. In this paper, we review recent progress in Lattice QCD calculations of the matrix elements of these CP violating operators. We start by iscussing the contribution of the QCD Θ-term to the nedm an use it to illustrate the various methos use in the calculations. 2 QCD Θ-term There are three strategies for the calculation of the contribution of CP violating operators to the neutron EDM on the lattice. We illustrate these three approaches for the case of the CP-violating Θ-term whose volume integral is the topological charge Q top. They are the external electric fiel metho [19 22], expansion in θ [23 27], an simulation with imaginary θ [28 3]. 2.1 External Electric Fiel Metho The first metho uses an external electric fiel [19 22], in which the neutron EDM, N, is extracte from the energy ifference of the neutron states with spin S aligne [anti]parallel to the external electric fiel E: E θ S Eθ S 2 Nθ S E. (2) For a given value of θ, the effect of the Θ-term is inclue by reweighting the nucleon correlation function with the topological charge: N N θ ( E, t) = N(t) N()e iθq top E, (3) where the expectation value on the r.h.s is evaluate on lattice configurations generate without the Θ-term. This metho extracts the neutron EDM from a spectral quantity, not the form factor F 3 as iscusse next, so the results are not affecte by the mixing uner parity violation problem iscusse in Section 3.

4 2.2 Expansion in θ Base on the current experimental boun of the neutron EDM an moel stuies, we know that the coupling θ is tiny. Hence it can be treate as a small expansion parameter. Any expectation value in presence of a small QCD Θ-term can then be written as follows [23 27]: O(x) θ = 1 [U, q, q]o(x)e S QCD+iθQ top Z θ = O(x) θ= + iθ O(x)Q top θ= + O(θ2 ). (4) In the secon line of the above equation, the QCD Θ-term effect is inclue at the lowest orer by measuring the correlation of the topological charge with the observable O, using configurations generate without the Θ-term. To calculate the nedm, the operator O is chosen as the matrix element of the vector current within the neutron states. Assuming PT invariance, the matrix element, calculate as a function of q µ, the momentum transfere by the vector current, can be ecompose in terms of Lorentz covariant form factors: [ NV µ N θ = ū F 1 (q 2 )γ µ + i F 2(q 2 ) σ µν q ν F 3(q 2 ] ) σ µν q ν γ 5 u. (5) 2m N 2m N Of these, the form factor of interest is F 3 at zero-momentum transfer. This is extracte by extrapolating F 3 (q 2 ), measure at finite q 2, because the term containing F 3 only contributes at q ν as can be seen from Eq. (5). Then the contribution to N is given by F 3 (q 2 ) N = lim. (6) q 2 2m N This approach, which relies on calculating the form-factor F 3 from nucleon matrix elements, requires a careful hanling of the neutron spinors in a theory with P violation. Most previous calculations i not correctly account for this subtlety [31]. The problem can be correcte retroactively, an in Figure 3 we show the resulting reuction in the value reporte in recent calculations [26, 27]). 2.3 Imaginary θ Simulation The Θ-term in Eucliean space-time is purely imaginary, an Monte Carlo simulations incluing the Θ-term face the sign problem. One way to circumvent this issue is to make the action real by taking θ to be purely imaginary [28 3]. Assuming that the theory is analytic aroun θ =, the results can then be continue to real θ for small θ. Furthermore, this calculation can be carrie out using the axial anomaly whereby the QCD Θ-term is chirally rotate to the fermionic term S q θ as θ i θ, S q θ = m l m s θ q(x)γ 5 q(x). (7) 2m s + m l Lattice simulation can then be performe by aing S q θ to the original QCD action. In this approach, the EDM is again extracte from the form factor analysis, so the results are affecte by the mixing inuce by parity violation iscusse in Section 3. Again, in Figure 3 we show one of the recent results [3] with an without the correction ue to this mixing problem. x

5 R Figure 2. CP violating phase α (enote as α 1 in the figure) as a function of the cutoff raius R presente in Ref. [32]. The central value of the α appears to saturate after R is larger than about 15 lattice spacings, which correspons to 1.71 fm in physical unit, while the statistical error grows as R is increase. 2.4 Variance Reuction using Cluster Decomposition In the calculation of neutron EDM inuce by the QCD Θ-term using θ-expansion or external electric fiel metho, one nees to calculate correlators reweighte by the topological charge. Recently, the authors of Ref. [32] reporte an error reuction technique using cluster ecomposition. The main iea is to calculate the topological charge only in the vicinity of the sink in the momentum projection summation. Let us consier a neutron two-point correlator reweighte with the topological charge Q: C R (t) = N( x, t) N(G, )Q R (x), Q R (x) = x x y <R q(y), (8) where G is the source gri. In conventional calculations, Q is calculate by summing the topological charge ensity q(y) from all lattice points. In the new formula, however, the topological charge is calculate only up to a certain raius R from the sink point uner the assumption that correlations ue to long-istance points are negligible. Figure 2 presents the CP violating phase, labele α 1 an calculate using the Eq. (8), as a function of the cutoff raius R. Data show that the mean value of α 1 saturates at some reasonably small value of R. Incluing the contributions of the topological ensity from larger R only increases the statistical noise. Consiering that the mean value still exhibits significant fluctuations, a more etaile costbenefit analysis of the approximation is require. 3 Extracting F 3 in a Theory with Parity Violation In a theory with P an CP symmetry, the spinor u of the neutron state satisfies the following Dirac equation: ( ipµ γ µ + m ) u =, (9)

6 with γ 4 the parity operator: u p γ 4 u p uner P. When CP, but not PT, is violate, the Dirac equation is only moifie by a phase factor to (ipµ γ µ + me 2iαγ5 ) ũ =, (1) where the spinor that solves the moifie equation is efine to be ũ. Furthermore, for this equation, γ 4 is no longer the parity operator, rather uner parity ũ p e 2iαγ 5 γ 4 ũ p. The two solutions u an ũ are relate as ũ = e iαγ 5 u, ũ = ūe iαγ 5. (11) When calculating the form factors F 2 an F 3, it is important to enforce that F 2 is the parity-even magnetic form factor while F 3 is the P an CP o electric form factor. There are two ways to enforce these symmetry requirements. The first is to properly inclue the phase efine in Eq. (11) in the efinition of the matrix element. In that case the ecomposition into the form factors is the same as in a theory with CP symmetry with unmoifie spinors u. Alternatively, one can calculate the stanar matrix element, an then uno the mixing between F 2 an F 3 ue to the non-stanar parity transformation with phase α. Then, the two form factors F 2 an F 3 are given in terms of F 2 an F 3 as in Eq. (14). In both cases, the phase α has to be extracte from the nucleon 2-point function. In this extraction, it is important to note that the phase α is state epenent. In Section 5, we show that the α corresponing to the groun state is given by the long-time behavior of the nucleon 2-point function, i.e., when the groun state ominates. In the first approach, one inclues the phase in the calculation of the n-point functions: e iαγ 5 ÑV µ Ñ e iαγ 5 = NV µ N, e iαγ 5 Ñ Ñ e iαγ 5 = N N ( = ū F 1 (q 2 )γ µ + i [γ µ, γ ν ] F 2 (q 2 ) q ν [γ µ, γ ν ] F 3 (q 2 ) ) q ν γ 5 u. (12) 2 2m N 2 2m N The resulting F 3 ()/2M N is the esire contribution to the nedm. In the secon approach, one calculates ( ÑV µ N, Ñ N = ũ F 1 (q 2 )γ µ + i [γ µ, γ ν ] F 2 (q 2 ) q ν [γ µ, γ ν ] F 3 (q 2 ) ) q ν γ 5 ũ. (13) 2 2m N 2 2m N an one has to extract F 3 from the following mixing structure: F 2 = cos(2α) F 2 sin(2α) F 3, F 3 = sin(2α) F 2 + cos(2α) F 3. (14) The two approaches are equivalent. Unfortunately, as pointe out in Ref. [31], all calculations base on the F 3 extraction prior to their work, starting with the 25 work in Ref. [23], use the secon approach but i not carry out the rotation efine in Eq. (14) to get the F 3. Their quote results are, therefore, for F 3. The authors of Ref. [31] show that correcting for this omission reuces the value of an alreay har to measure F 3 by about a factor of ten. Thus, all previous estimates of the contribution of the Θ-term an the cedm base on evaluating F 3 nee to be revise. Figure 3 shows some of the recent lattice results of F 3 inuce by the QCD Θ-term, before an after the parity mixing correction. The interesting point is that all correcte numbers are close to zero. Previous phenomenological estimates were an orer of magnitue smaller than lattice QCD results. This tension may isappear after a proper analysis of the phase α introuce in the neutron states ue to the breaking of parity. Further iscussion of the mixing between F 2 an F 3 when parity is violate can be foun in Ref. [31].

7 F Phenomenology Shintani 216 Alexanrou 216 Guo 215 Shintani Blue: Correcte F 3 Re : Ol F ~ m π (MeV) Figure 3. Recent lattice results of F 3 inuce by the QCD Θ-term calculate with M π < 55 MeV. Re points are the original results, which give F 3, an the blue points are the correcte results. Original results are obtaine from Refs. [27] (labele Shintani 216), [26] (Alexanrou 216), [3] (Guo 215), an [23] (Shintani 25), an the correcte results are obtaine from Refs. [31, 33]. As note in Ref. [31], some assumptions are mae to calculate the correcte F 3 numbers as some of the original papers o not contain full information neee for the parity inuce rotation Eq. (14), so the F 3 results presente in this plot may not be precise. 4 quark EDM The electromagnetic current, V µ δl/δa µ, gets an aitional contribution when the quark EDM operator is ae to L SM. Thus, while the stanar conserve vector current interacts with a quark with charge q f, this aitional contribution, which is is just the flavor iagonal tensor operator for each quark flavor, couples with strength γ f g f T, where γ f are the BSM moel epenent CP violating couplings an g u T, gs T, gs T,... are the tensor charges, i.e., matrix elements of the quark bilinear tensor operator within the nucleon states. Furthermore, the effects analogous to the mixing iscusse in Section 3 are suppresse by powers of the electromagnetic coupling α EM 1/137. Thus, the leaing contribution of the EDM of the quarks to the nedm is given by the flavor iagonal tensor charges: N qσµν q N = g q T ūnσ µν u N, N = u g u T + g T + sg s T. (15) In many moels, such as the supersymmetric moels, the quark EDMs are proportional to the corresponing quark masses ( q m q ), because in these theories all connections between left an righthane quarks are meiate by a common set of Yukawa interactions. Since the strange quark is much heavier than the u an the quarks, the strange quark contribution gets enhance, by m s /m 2. On the other han gt s gl T. Since the contribution is the prouct of the two, it is, therefore, important to etermine the strange quark, an perhaps even the charm quark, tensor charge precisely. There are two classes of iagrams in the calculation of flavor iagonal tensor charges: quarkline connecte an the quark-line isconnecte iagrams as illustrate in Figure 4. In case of the strange (charm) quark tensor charge, only the isconnecte iagram with strange (charm) quark loop contributes. This is the basis of the expectation that gt s gl T. In Refs. [34, 35], the tensor charges were calculate using clover fermions on the HISQ lattices, an extrapolate to continuum an physical pion mass limit, as illustrate in Figure 5 (left). The is-

8 Figure 4. Quark-line connecte (left) an isconnecte (right) iagrams contribute to the quark EDM. Black ellipses are the neutrons, an the cross is the quark bilinear operator with tensor structure. - - ( ) u 1 25 e cm n e cm 9 CL g T s e cm Figure 5. (Left) Extrapolation of g s T to the continuum limit. The fit is performe simultaneously in a an M π but only the projection to the a-plane with ata extrapolate to the physical M π, is presente. (Right) Bouns on u an are obtaine using the formula given in Eq. (15), the tensor charges given in Eq. (16) an the experimental bouns N < e cm. Since the strange quark tensor charge is consistent with zero within statistics, it imposes no constraint on s. Both figures are taken from Ref. [34]. connecte iagrams were estimate using a stochastic metho, an it turne out that the isconnecte iagram contribution to the tensor charges is very small. Their results are g u T =.23(3), g T =.79(7), gs T =.8(9). (16) There are BSM moels in which the quark EDMs are the ominant sources of CP violation at low energy. In these scenarios, constraints can be place on the couplings of the EDM of each quark flavor using the experimental boun on the neutron EDM an the lattice results for the tensor charges. Neglecting all other sources of CP violation, figure 5 (right) shows the allowe region of u an using the lattice results, Eq. (16), escribe in Ref. [34, 35]. Of all the CP violating operators liste in Eq. (1), the calculation of the quark EDM, incluing renormalization, is theoretically wellestablishe an the results in Eq. (16) for g u T an g T inicate that these are alreay available at the 15% level. Further improvement in precision will occur as higher precision ata are generate in future calculations.

9 5 Quark Chromo-EDM The last imension-5 term is the quark chromo-electric ipole moment (cedm) operator. The moification to the action in presence of the CP violating cedm term is: S = S QCD i 4 x q q(σ G)γ 5 q, (17) 2 where q are again BSM couplings evolve to the haronic scale for each quark flavor. The lattice calculation consists of evaluating the matrix element of the insertion of the prouct of the electromagnetic current V µ an the q(σ G)γ 5 q operator within the neutron states. Then, the contribution of cedm operators to the nedm, ignoring the complicate operator mixing problem uner renormalization iscusse in Section 5.4 below, is obtaine from q q N V µ q(σ G)γ 5 q N. Lattice QCD stuies of the cedm operator have starte only recently, an three methos are being explore: an expansion in q [31], external electric fiel metho [31], an the Schwinger source metho [36, 37]. 5.1 Expansion in q The expansion in q metho procees analogously to the expansion in θ in the QCD Θ-term calculation. For small BSM cedm couplings, q, the neutron matrix element of the electromagnetic current can be written in terms of the expectation values evaluate on the stanar CP even lattices, NVµ N CPV = NV µ N + q NVµ N O cedm (x) + O( q) 2, (18) x O cedm = i 2 q(σ G)γ 5q. (19) Their contribution to the neutron EDM is extracte from the CP violating form factor, F 3 as efine in Eq. (5). The calculation is challenging because it requires calculating the expectation value of fourpoint functions because one has to insert both the vector current an the cedm operator between the neutron source an sink, with the cedm operator efine as a sum over all space-time points. Figure 6 (top) shows all the require quark-line connecte iagrams. These four-point correlators are constructe using five types of propagators use as builing blocks that are shown in the bottom panel of Figure 6 (bottom). F an B are the usual forwar an backwar propagators. C is the cedmsequential propagator. It is constructe by starting with the regular quark propagator an inserting the chromo-edm operator at all lattice sinks points. This is then use as the source for calculating the chromo-edm inserte propagator. Results from Ref. [31] for the signal in F 3 an the source-sink separation epenence of the signal are shown in Figure 7 (left). F 3 has larger value when the cedm is inserte in the -quarks of the neutron. In Figure 7 (right), the results are extrapolate to zero-momentum to obtain the contribution F 3 () as neee for the neutron EDM. 5.2 External Electric Fiel Metho The chromo-edm contribution to the neutron EDM also can be calculate using the external electric fiel. Similarly to the QCD Θ-term, the neutron EDM is extracte from the energy ifference of nucleon states uner spin flip in presence of a uniform external electric fiel ie: E cedm S E cedm S 2 N S i E, (2)

10 z z y y y x z z x x z Figure 6. (Top) Connecte iagrams neee for the calculation of the cedm operator. Crosses are the cedm operator insertion, an ots are the vector current insertion. (Bottom) Propagators neee for the quark-line connecte four-point correlators. In this figure, the electromagnetic current is enote as J. Figures are taken from Ref. [31]. where the neutron correlators in the presence of the cedm term are obtaine by reweighting, N N cedm (i E, t) = N(t) N()O cedm i E. (21) This metho nees only two- an three-point functions of the neutron, the latter with the insertion of the cedm bilinear operator. In Ref. [31], F 3 is extracte using the external electric fiel metho on the same ensemble of lattices use for the calculation using the q -expansion metho, an the results are presente in Figure 7 (right). The signal is poorer in the external electric fiel metho compare to those in expansion in q metho, but the results of the two methos are consistent. Note that the external electric fiel metho is not affecte by the parity mixing problem escribe in Section 3, while the q -expansion metho is. Therefore, the consistency of the results from the two methos, after taking care of the parity mixing rotation given in (14), suggests that both methos are yieling a reliable signal. 5.3 Schwinger Source Metho Another metho to calculate the contribution of the cedm operator to the neutron EDM is the Schwinger source metho [36, 37]. Noting that the cedm operator is a quark bilinear, i q(σ G)γ 5 q, one can a it to the QCD fermion action: D clov D clov + iεσ µν γ 5 G µν, (22) where D clov is the Dirac operator for the Wilson-clover action, an ε is a small control parameter the epenence on which will be remove by taking a erivative with respect to it at the en of the calculation. Because cedm is a quark bilinear, the integration over the fermion egrees of freeom in the path integral can still be carrie out as before. In the case of the clover fermions, the aition of the cedm operator is equivalent to shifting the coefficient of the imension-5 clover term by iεγ 5 : c sw σ µν G µν σ µν (c sw + iεγ 5 )G µν. (23)

11 F 3n, (cedm) U T = 8a T = 1a T = 12a Q 2 =.21 GeV 2 F 3n, (cedm) D Q 2 =.21 GeV 2 Q 2 =.41 GeV 2 Q 2 =.41 GeV 2 rotation givenq in 2 =.59 (14). GeV 2 Same lattices use Q 2 =.59 in GeV (left) 2 plot are use for this calculation. E =.39 GeV 2. Note that 1 2 these 1results are not renormalize. Figures are taken from Ref. [31] t op t sep /2 t op t sep /2 F3n, (cedm) U F3n, (cedm) D E/E = ±1 E/E = ±2 NEW F3(T = 8) NEW F3(T = 1) OLD F3(T = 8) OLD F3(T = 1) Q 2 [GeV 2 ] Figure 7. (Left) Plateau plots for F 3 inuce by the cedm operator calculate by expaning in q. Simulations were performe on , a =.11 fm an M π = 34 MeV DWF lattices. (Right) Results for cedm contribution to F 3 extrapolate to Q 2 = from the q -expansion results compare with the F 3 obtaine by using the external electric fiel metho. OLD an NEW F 3 label the results before an after fixing the parity mixing In this approach, the insertion of the cedm term is carrie out by calculating valence quark propagators incluing this moifie clover term, an using this propagator in all iagrams requiring a cedm insertion. Having taken care of the cedm term, the calculation of the original 4-point functions reuces to three-point functions, with an insertion of the vector current within the neutron state, mae up of propagators with an without cedm insertion. To perform this calculation on ensembles generate with just the stanar QCD action, however, requires taking into account the change in the Boltzmann factor uner the aition of the cedm term to the action as shown in Eq. (22). This, a priori, non-unitary formulation between sea an valence quarks can be accounte for at leaing orer in ε by reweighting each configuration by ratio of the fermion eterminant with an without the cedm term: et ( ) D clov + iεσ µν γ 5 G µν et ( ) exp [ iεtr ( σ µν γ 5 G µν Dclov)] 1. (24) D clov Diagrammatically, all the terms that nee to be calculate in this Schwinger source metho are shown in Figure 8. There are three classes of iagrams at leaing orer in ε: the reweighting factor for each configuration, an the quark-line connecte an isconnecte iagrams on that configuration. Current calculations are at the stage of emonstrating a signal in the full evaluation of Figure 8 with both the cedm an the γ 5 operator with which it mixes uner renormalization as iscusse in Section 5.4. First, in Figure 9, we show the success at extraction of the phase α that is inuce in the groun state nucleon spinor by the CP violating cedm an γ 5 operators. In the right panel of Figure 9, we also show that this α is linear in ε, which allows us to tune the value of ε. We want to use a large ε to increase the signal but remain in the linear response regime. First examples of the quality of the signal in the contribution of the connecte iagrams to the form factor F 3 inuce by the cedm an γ 5 operators are shown in Figures 1 for two ifferent source an sink separations. At present, the signal is consistent with zero for both operators.

12 Figure 8. Quark line iagrams neee for the Swinger source metho. The reweighting factor (left), the quarkline isconnecte (mile four), an the quark-line connecte (right four) iagrams. Re crosse box is the cedm operator insertion, P is the regular quark propagator without the cedm operator insertion, an P ε is the quark propagator with the cedm operator insertion. An ientical calculation has to be one with the cedm operator replace by the γ 5 operator in orer to efine a finite renormalize quantity..5 a9m a9m α cedm.4.35 ε = t α γ ε γ5 = t α N a12m31 cedm a9m31 cedm a12m31 γ 5 a9m31 γ ε Figure 9. The left panel shows the extraction of the CP violating phase α inuce by the cedm operator from the long-time behavior of the 2-point correlator. The mile panel shows the extraction of the phase for the γ 5 operator insertion. The right panel shows that both kins of phases are linear in ε when ε is small. These ata are from the Clover-on-HISQ analysis on the M π 31 MeV ensembles with lattice spacing a =.12 (a12m31) an.9 (a9m31) fm. 5.4 Renormalization of cedm Operator The renormalization of the cedm operator is stuie both in one-loop perturbation theory with Twiste-mass fermions [38], an in nonperturbative RI-SMOM scheme [39]. The most challenging outcome is the ivergent mixing with lower-imensional operators: in particular the 1/a 2 mixing with the pseuoscalar quark bilinear operator: O cedm = a 2 qσ µν γ 5 G µν q, O P = qγ 5 q. (25) Because the mixing is 1/a 2 ivergent, it nees to be calculate precisely so that a finite operator can be constructe, ensemble by ensemble, by subtracting the two terms. In Ref. [38], the authors presente the results of the 1/a 2 mixing coefficients for the Twiste-mass fermions calculate nonperturbatively an using the one-loop perturbation theory. In Ref. [39], authors efine a momentum subtraction scheme, RI- SMOM, for the non-perturbative renormalization of the cedm operator on the lattice, an provie one-loop matching coefficients to the MS scheme in continuum limit.

13 ε F 3 D, cedm, -u Q 2 =.18 GeV 2 ε =.3 a9m31 t sep =1 t sep = τ - t sep /2 ε F 3 D, γ 5, -u Q 2 =.18 GeV 2 ε =.3 a9m31 t sep =1 t sep = τ - t sep /2 Figure 1. (Left) Signal of F 3 form-factor for two ifferent separation of source an sink on a9m31 lattices. cedm is inserte on -quarks. Only the connecte iagrams are consiere. (Right) Same plot as the (mile), but γ 5 is inserte instea of the cedm operator. Results shown are for the bare operators. 6 Summary In this talk, we reviewe the recent lattice QCD calculations of the contribution of imension four an five CP violating operators to the EDM of the neutron. We reviewe the three approaches, external electric fiel metho, expansion in θ, an simulation with imaginary θ, use for the calculation of the QCD Θ-term in Section 2. We also summarize the observation mae in Ref. [31] concerning the subtlety of properly efining the neutron spinor in the presence of parity violating interactions an its consequences vis-a-vis mixing between the form factors F 2 an F 3 in Section 3. All calculations base on extracting F 3, previous to Ref. [31], ha misse this mixing. After correction, current lattice QCD results for the QCD Θ-term, are reuce by about a factor of ten an o not show a statistically significant non-zero signal. The quark EDM contribution to the neutron EDM can be written in terms of the neutron tensor charges g T q. Currently, the g T u, are etermine to within 15% uncertainty incluing the quark-line isconnecte iagrams, an g T s is known to be very small. The results are summarize in Section 4. Lattice QCD calculation for the quark chromo-edm operator have just starte. Three methos for the calculation of the contribution of cedm an the construction of a finite renormalize operator are escribe in Section 5. The methoology for calculation of the 3-point (or 4-point) correlation functions is uner control but a non-zero signal in all the iagrams has yet to be emonstrate. Once a signal is achieve, we will still be a long way away from a full calculation ue to the ivergent mixing between the cedm an γ 5 operators. A few brief wors on other lattice QCD calculations relate to the neutron EDM that have not been covere in this review. First, preliminary results for the Weinberg three-gluon operator, one of the imension-six operators, were presente at the Lattice 217 conference [4]. Calculations of the Weinberg operator are at the stage of emonstrating a signal an have not even begun to aress the issue of renormalization that is potentially even more complicate than that for the cedm operator. Calculations exploring the four-fermion imension six operators are yet to begin. Secon, there are efforts to stuy the neutron EDM using lattice spectroscopy, an other matrix elements combine with chiral perturbation theory [41], in aition to the irect calculation of the matrix elements with the CP violating operators. To summarize, the theoretical calculations of matrix elements neee to exten the importance an reach of nedm experiments to constrain BSM physics have begun in earnest but are very challenging an will require many new innovations in both theory an computations.

14 Acknowlegments We thank Keh-Fei Liu, Gerrit Schierholz, Eigo Shintani, an Sergey Syritsyn for valuable iscussions an for sharing results. The work was supporte by the U.S. DoE HEP Office of Science contract number DE-KA-1412 an the LANL LDRD Program. The simulations for the qedm contribution an the cedm contribution using Swinger source metho are carrie out on computer facilities at (i) the USQCD Collaboration, which are fune by the Office of Science of the U.S. Department of Energy, (ii) the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supporte by the Office of Science of the U.S. Department of Energy uner Contract No. DE-AC2-5CH11231, (iii) Oak Rige Leaership Computing Facility at the Oak Rige National Laboratory, which is supporte by the Office of Science of the U.S. Department of Energy uner Contract No. DE-AC5- OR22725, an (iv) Institutional Computing at Los Alamos National Laboratory. References [1] J.M. Penlebury et al., Phys. Rev. D92, 923 (215), [2] Search for the neutron electric ipole moment at PSI, [3] R. Picker, JPS Conf. Proc. 13, 15 (217), [4] A next generation measurement of the electric ipole moment of the neutron at the FRM-II, [5] Search for a neutron electric ipole moment at the SNS, positionpapers/fsn_position_nedm_214_5.pf [6] A cryogenic experiment to search for the EDM of the neutron, inex_files/cryoedm.htm [7] S. Dar (2), hep-ph/8248 [8] M. Pospelov, A. Ritz, Annals Phys. 318, 119 (25), hep-ph/54231 [9] M.J. Ramsey-Musolf, S. Su, Phys. Rept. 456, 1 (28), hep-ph/61257 [1] J. Engel, M.J. Ramsey-Musolf, U. van Kolck, Prog.Part.Nucl.Phys. 71, 21 (213), [11] A. Sakharov, Pisma Zh.Eksp.Teor.Fiz. 5, 32 (1967) [12] W. Bernreuther, Lect. Notes Phys. 591, 237 (22), [,237(22)], hep-ph/25279 [13] M. Pospelov, A. Ritz, Nucl. Phys. B573, 177 (2), hep-ph/99858 [14] R.J. Crewther, P. Di Vecchia, G. Veneziano, E. Witten, Phys. Lett. 88B, 123 (1979), [Erratum: Phys. Lett.91B,487(198)] [15] A. Abaa, J. Galan, A. Le Yaouanc, L. Oliver, O. Pene, J.C. Raynal, Phys. Lett. B256, 58 (1991) [16] E. Mereghetti, J. e Vries, W.H. Hockings, C.M. Maekawa, U. van Kolck, Phys. Lett. B696, 97 (211), [17] R.D. Peccei, H.R. Quinn, Phys. Rev. Lett. 38, 144 (1977) [18] R.D. Peccei, H.R. Quinn, Phys. Rev. D16, 1791 (1977) [19] S. Aoki, A. Gocksch, Phys. Rev. Lett. 63, 1125 (1989), [Erratum: Phys. Rev. Lett.65,1172(199)] [2] S. Aoki, A. Gocksch, A.V. Manohar, S.R. Sharpe, Phys. Rev. Lett. 65, 192 (199) [21] E. Shintani, S. Aoki, N. Ishizuka, K. Kanaya, Y. Kikukawa, Y. Kuramashi, M. Okawa, A. Ukawa, T. Yoshie, Phys. Rev. D75, 3457 (27), hep-lat/61132 [22] E. Shintani, S. Aoki, Y. Kuramashi, Phys. Rev. D78, 1453 (28),

15 [23] E. Shintani, S. Aoki, N. Ishizuka, K. Kanaya, Y. Kikukawa, Y. Kuramashi, M. Okawa, Y. Tanigchi, A. Ukawa, T. Yoshie, Phys. Rev. D72, 1454 (25), hep-lat/5522 [24] F. Berruto, T. Blum, K. Orginos, A. Soni, Phys. Rev. D73, 5459 (26), hep-lat/5124 [25] A. Shinler, T. Luu, J. e Vries, Phys. Rev. D92, (215), [26] C. Alexanrou, A. Athenoorou, M. Constantinou, K. Hajiyiannakou, K. Jansen, G. Koutsou, K. Ottna, M. Petschlies, Phys. Rev. D93, 7453 (216), [27] E. Shintani, T. Blum, T. Izubuchi, A. Soni, Phys. Rev. D93, 9453 (216), [28] T. Izubuchi, S. Aoki, K. Hashimoto, Y. Nakamura, T. Sekio, G. Schierholz, PoS LAT27, 16 (27), [29] S. Aoki, R. Horsley, T. Izubuchi, Y. Nakamura, D. Pleiter, P.E.L. Rakow, G. Schierholz, J. Zanotti (28), [3] F.K. Guo, R. Horsley, U.G. Meissner, Y. Nakamura, H. Perlt, P.E.L. Rakow, G. Schierholz, A. Schiller, J.M. Zanotti, Phys. Rev. Lett. 115, 621 (215), [31] M. Abramczyk, S. Aoki, T. Blum, T. Izubuchi, H. Ohki, S. Syritsyn, Phys. Rev. D96, 1451 (217), [32] K.F. Liu, J. Liang, Y.B. Yang (217), [33] Private communication with Eigo Shintani. [34] T. Bhattacharya, V. Cirigliano, R. Gupta, H.W. Lin, B. Yoon, Phys. Rev. Lett. 115, 2122 (215), [35] T. Bhattacharya, V. Cirigliano, S. Cohen, R. Gupta, A. Joseph, H.W. Lin, B. Yoon (PNDME), Phys. Rev. D92, (215), [36] T. Bhattacharya, V. Cirigliano, R. Gupta, E. Mereghetti, B. Yoon, PoS LATTICE215, 238 (216), [37] T. Bhattacharya, V. Cirigliano, R. Gupta, B. Yoon, PoS LATTICE216, 225 (216), [38] M. Constantinou, M. Costa, R. Frezzotti, V. Lubicz, G. Martinelli, D. Meloni, H. Panagopoulos, S. Simula, Phys. Rev. D92, 3455 (215), [39] T. Bhattacharya, V. Cirigliano, R. Gupta, E. Mereghetti, B. Yoon, Phys. Rev. D92, (215), [4] J. Dragos, A. Shinler, T. Luu, J. e Vries (217), Lattice 217 conference, ugr.es/event//session/95/contribution/249 [41] J. e Vries, E. Mereghetti, C.Y. Seng, A. Walker-Lou, Phys. Lett. B766, 254 (217),