The Contribution of Novel CP Violating Operators to the nEDM using Lattice QCD

In this talk, we motivate the calculation of the matrix elements of novel CP violating operators, the quark EDM and the quark chromo EDM operators, within the nucleon state using lattice QCD. These matrix elements, combined with the bound on the neutron EDM, would provide stringent constraints on beyond the standard model physics, especially as the next generation of neutron EDM experiments reduce the current bound. We then present our lattice strategy for the calculation of these matrix elements, in particular we describe the use of the Schr\"odinger source method to reduce the calculation of the 4-point to 3-point functions needed to evaluate the quark chromo EDM contribution. We end with a status report on the quality of the signal obtained in the lattice calculations of the connected contributions to the quark chromo EDM operator and the pseudoscalar operator it mixes with under renormalization.


Introduction: Baryogenesis and the need for novel CP violation
The observed universe has 6.1 +0. 3 −0.2 × 10 −10 baryons for every black body photon [1], whereas in a baryon symmetric universe, we expect no more that about 10 −20 baryons for every photon [2]. It is difficult to include such a large excess of baryons as an initial condition in an inflationary cosmological scenario [3]. The way out of the impasse lies in generating the baryon excess dynamically during the evolution of the universe.
In the early history of the universe, if the matter-antimatter symmetry was broken post inflation and reheating, then one is faced with Sakharov's three necessary conditions [4]: the process has to violate baryon number, evolution has to occur out of equilibrium, and charge-conjugation and T (or equivalently CP if CPT remains a good symmetry) invariance has to be violated.
CP violation (CP) exists in the standard model (SM) of particle interactions due to a phase in the Cabibo-Kobayashi-Maskawa (CKM) quark mixing matrix [5], and possibly by a similar phase in the leptonic sector, given that the neutrinos are not massless [6].CP arising from the CKM matrix in the SM contributes at O(10 −32 ) e-cm, much smaller than the current experimental bound d n < 2.9 × 10 −26 e-cm. Thus the nEDM puts no constrain on the SM and the strength of theCP in the CKM matrix is much too small to explain Baryogenesis. In principle, the SM has an additional source ofCP arising from the effect of QCD instantons. The presence of these finite action non-perturbative configurations in a non-Abelian theory often leads to inequivalent quantum theories defined over various 'Θ'-vacua [7]. Because of asymptotic freedom, a Speaker. e-mail: rajan@lanl.gov. Los Alamos Report LA-UR-16-29579 arXiv:1701.04132v1 [hep-lat] 16 Jan 2017 EPJ Web of Conferences all non-perturbative configurations including instantons are strongly suppressed at high temperatures where baryon number violating processes occur. Because of this,CP due to such vacuum effects do not lead to appreciable baryon number production. In short, additional much largerCP is needed from physics beyond the SM (BSM).
To determine whether such additionalCP exists, a very promising approach is to measure the static electric dipole moments of elementary particles, atoms and molecules, which are necessarily proportional to their spin. Since under time-reversal the direction of spin reverses but the electric dipole moment does not, a non-zero measurement confirms T violation or equivalentlyCP. Of the elementary particles, atoms and nuclei that are being investigated, the electric dipole moment of the neutron (nEDM) is the laboratory where lattice QCD can provide the theoretical part of the calculation needed to "connect" the experimental bound (value) on the nEDM to the strength ofCP in a given BSM theory.
Most extensions of the SM have new sources ofCP. Each of these contributes to the nEDM and for some models it can be as large as 10 −26 e-cm. Planned experiments are aiming to reduce the bound on d n < 2.9 × 10 −26 e-cm to d n 3 × 10 −28 e-cm. The strategy for finding out which class of BSM theories are viable candidates is as follows: As the bound on the nEDM is lowered in current and planned experiments, BSM theories withCP giving rise to a nEDM larger than this bound get ruled out provided the "connection" between the bound and the couplings is known sufficiently accurately. This "connection" is the matrix elements of the novel CP violating interactions within the neutron state that simulations of lattice QCD are gearing up to provide.

QCD Lagrangian and the EM current in the presence of CP Violation
Many candidate BSM theories have been proposed by theorists. While the true BSM theory is not known, one can write down all possibleCP interactions at the energy scale of hadronic matter (∼ few GeV) in terms of quark and gluon fields based on symmetry and organized by their dimension; operators with higher dimension being, in general, suppressed. At the lowest dimension five, there are two leading operators called the quark EDM (qEDM) and the quark chromo EDM (CEDM). The QCD Lagrangian in the presence of the Θ-term, the qEDM and the CEDM operators is where the sum is over all the quark flavors, Σ µν = [γ µ , γ ν ]/2, and F µν and G µν are the electromagnetic and chromo field strength tensors. The couplings d γ q are the qEDMs and the d G q are the CEDMs. They parameterize newCP in BSM theories. The goal is to constrain BSM models by bounding these couplings using the experimental bound on the neutron EDM and lattice QCD calculations of the matrix elements of the corresponding operators. In other words, the matrix elements of the electromagnetic current J EM µ between neutron states in the presence of CP violation provides the "connection" between the BSM couplings and the nEDM.
In the presence ofCP, the electromagnetic current, defined as δL/δA µ , gets an additional term: The matrix elements of this leading qEDM (second) term are the flavor diagonal tensor charges g q T : The connected and disconnected (with u, d, s and c quark loops) Feynman diagrams for the qEDM contribution are shown in Fig. 1 (left). The first lattice results for the qEDM were given in [8] and phenomenological consequences for a particular BSM theory, Split SUSY, were presented in Ref. [9].
The key computational challenge to evaluating Eq. (3) is the calculation of the disconnected contributions to g q T . These were shown to be small and noisy, and decrease with the quark mass [8]. The errors in the strange disconnected contribution were small enough to yield the continuum limit estimate g s T = 0.008(9) [8]. In most promising BSM theories, the newCP interactions are Yukawa like, i.e., the couplings are proportional to the quark mass. Consequently, the impact of the small value of g s T with O(1) error, gets magnified by the quark mass ratio 2m s /(m u + m d ) = 27. Thus, it is important to improve the estimate of g s T and, for the same reason, calculate g c T .
Contributions from the CEDM and the next order in quark EDM arise due to the change in the action, In an ideal world, the best way to calculate these effects would be to simulate the lattice theory using a discretized version of L & & CP QCD and compute the matrix element This ideal approach is not practical becauseCP interactions are complex and lattice simulations of theories with a complex action are not yet realistic.
The method of choice is to treat theseCP interactions as perturbations in the small couplings d γ q and d G q , and expand the theory about L QCD (see Sec. 3). Then, the leading terms that contribute to the nEDM are the matrix elements of the product of the electromagnetic current J EM µ and the 4-volume integral of the operators. For the CEDM operator, the matrix elements that have to be calculated are as discussed in Sec. 3. Similarly for the next order qEDM This second qEDM term has two electromagnetic interactions and is therefore expected to be smaller than the leading term given in Eq. (3). For this reason we neglect its contribution. Finally, the contribution of the Θ-term requires calculating the matrix element: This correlation of J µ with topological charge Q gives rise to the 2 diagrams shown in Fig. 1 (right).
. CEDM contribution is the reweighting factor times the sum of connected and disconnected diagrams.

Simulating QCD Including Complex¨CP CEDM and γ 5 operators
The QCD Lagrangian becomes complex with the addition of any CP violating interaction. A robust cost-effective method for simulating theories with complex actions does not exist. To calculate the contribution of theCP interactions to the nEDM, we work within the framework of the Schwinger source method. We start by noting that the fermion part of the Lagrangian, L & & CP QCD in Eq. (1), involves only quark bilinear operators in presence of the qEDM and CEDM interactions. The fermion fields in the path integral can, therefore, still be integrated out. There are, however, two modifications to the standard calculations of 3-point functions due to the presence ofCP interactions treated as sources. First, the Dirac operator, in the presence of say the CEDM term, becomes Quark propagators calculated with this modified Dirac operator include the insertion of the CEDM operator at all space-time points. Second, the Boltzmann weight of all gauge configurations generated without the CEDM term in the action need to be modified by reweighting them by the ratio of the determinants of the Dirac operators for the two theories-with and withoutCP terms: Since the coupling for the CEDM term for all quark flavors is small, we expect the leading order reweighting factor to be a good approximation. This 'reweighting factor' for each gauge configuration is, to leading order, the integral over all spacetime points of the value of a closed quark loop with the insertion of the CEDM operator.
There is an additional challenge: The CEDM operator has an UV divergent mixing with the "γ 5operator, qγ 5 q/a 2 [10]. Thus, one has to (i) calculate the matrix element of the "γ 5 -operator", which we do in the same way as for the CEDM operator, (ii) calculate a similar reweighting factor, and (iii) handle the 1/a 2 divergent mixing. Since the reweighting is an overall phase, and configuration to configuration fluctuations can be large, it is therefore important to demonstrate the quality of the signal in F 3 (0) after reweighting for both the CEDM and γ 5 terms. In Sec. 5, we describe the progress CONF12 we have made towards addressing these challenges. There is also a mixing of the CEDM operator with the G µνGµν operator. Addressing this mixing is part of the calculation of the Θ−term per Eq. (6) [10].
The Schwinger source method has allowed us to recast the challenging calculation of 4-point functions given in Eq. (4) to the difficult calculation of 3-point functions [11]. The Feynman diagrams needed to calculate N|J EM µ (q)|N in the presence of CEDM and γ 5 interactions are illustrated in Fig. 2. Starting with an ensemble of gauge configurations generated with a standard lattice action, for example the Wilson-clover action, the steps in the calculation are: • Calculate propagators, labeled P in Fig. 2, using the standard Wilson-clover Dirac operator. We assume isospin symmetry so the propagators for u and d quarks are numerically the same.
• Calculate a second set of propagators with the CEDM term included with coefficient in the Wilsonclover matrix: These propagators, labeled P , include the full effect of inserting the CEDM operator at all possible points along them. The cost of inversion increases by about 7% with respect to P, however, using P as the starting guess in the inversion for P reduces the number of iterations required by 20-40% depending on the quark mass. The overall average cost of P is found to be about 80% of P.
• Using P and P , construct 4 kinds of sequential sources, labeled P • If using the coherent sequential source method, then N sources at different time-slices for the N different calculations being done simultaneously on that configuration, are added together. Sequential propagators with these (coherent) sequential sources are calculated by inverting the Dirac operator with and without the CEDM term as appropriate. The 4 types of (coherent) sequential propagators are shown by the block of 4 correlation functions in the right half of Fig. 2.
• Using the two original and the four sequential propagators, calculate the connected 3-point function. These contractions give the four 3-point functions shown on the right in Fig. 2. Note that a separate insertion of J EM µ is done on the u and d quark lines. • Calculate the disconnected quark loop with the insertion of electromagnetic current at zero momentum for each of the quark flavors, u, d, s, c (and b if needed) using the Dirac propagator with and without the CEDM term.
• Calculate the correlation of these quark loops with the appropriate nucleon 2-point correlation functions. The 4 types of disconnected Feynman diagrams required are illustrated in the left of Fig. 2.
• To correct for the omission of the chromo EDM operator in the action during the generation of the gauge configurations, calculate the volume integral of the quark loop with the CEDM insertion for each flavor. Multiply the sum of the connected and disconnected contributions by the reweighting factor-exponential of the estimate of the "CEDM loop" for the configuration with the appropriate coupling factor i . This is shown by the overall factor in Fig. 2.
• Repeat the calculation for different values of that bracket the expected numerical values of the couplings d G q for the various quark flavors. • Repeat the whole calculation for the γ 5 -operator instead of the CEDM operator.

Matrix Elements and Form Factors from 3-point Functions
The matrix elements, defined in Eqs. (3), (4) and (6), are extracted from the vacuum to vacuum 3-point correlation function: insertion of the electromagnetic current at times t between the neutron source and sink operators at time 0 and T : where in the right hand side we have used the expansion in a complete basis of states |N , |N , . . . that couple to the neutron interpolating operator N. We use N ≡ abc [d aT Cγ 5 (1 + γ 4 )u b ]d c where C = γ 0 γ 2 is the charge conjugation matrix, a, b, c are the color indices, u, d are the quark flavors, u N is the free neutron spinor and M N is the neutron mass. In the presence ofCP, the nucleon 2-point function is where the phase angle α N arises due to theCP coupling and depends on it. To determine the desired region of linearity of α versus we show the imaginary part of the neutron 2-point function in Fig. 3.
It demonstrates the quality of the signal for α for appropriately small values of on two different ensembles. Fig. 4 shows the expected linear behavior of α versus over a range of small . Once the matrix element of the electromagnetic current, J EM µ (q) defined in Eq. 2, within the nucleon state is extracted using Eq. (11), it is parameterized in terms of four Lorentz covariant form CONF12 Figure 4. Linear behavior of the phase α versus with the insertion of the CEDM and the γ 5 operators.
factors F 1 , F 2 , F A and F 3 : F 1 and F 2 are the standard Dirac and Pauli form factors. The anapole form factor F A violates parity P and the electric dipole form factor F 3 violates P and CP. F 3 is extracted from the different matrix elements by using different combinations of the momentum transfer q µ and spin projections. The zero momentum limit of these form factors gives the charges and dipole moments: the electric charge is F 1 (0) = 1 and the anomalous magnetic dipole moment is F 2 (0)/2M N . The contribution of the matrix element of eachCP interaction defined in Eqs. (3), (4), and (6) to the electric dipole moment of the neutron is given by d n = lim q 2 →0 F 3 (q 2 )/2M n .

Status of Numerical Calculations
So far, we have performed numerical calculations on two 2+1+1-flavor HISQ ensembles generated by the MILC collaboration [14]. The first, labeled a12m310, has a = 0.12 fm and M π = 310 MeV and the second, labeled a09m310, has a = 0.09 fm and M π = 310 MeV. The correlation functions are constructed using Wilson-clover fermions. On the a12m310 (a09m310) ensemble we have analyzed 400 (270) configurations with 64 measurements on each configuration. The ensembles used and our strategy for the calculation of 2-and 3-point functions with this clover-on-HISQ approach are described in Ref. [12]. In Fig. 5, we show the data for F 3 from the connected diagrams in the presence of the CEDM term, and in Fig. 6, the data for F 3 in the presence of the γ 5 term. The data are presented for three values of the source-sink separation t sep = 8, 10 and 12. The excited-state contamination in the matrix elements, and thus in F 3 , is removed by taking the t sep → ∞ limit. These figures show that an acceptable signal-to-noise ratio can be obtained to estimate F 3 with O(1) errors, our first goal for the CEDM calculations.
There are theoretical reasons to expect that the connected contribution of the γ 5 operator is, with small corrections, proportional to that of the CEDM operator [13]. In Fig. 7, we show the ratio of the two contributions and find that this expectation is actually realized. EPJ Web of Conferences  Figure 5. Illustration of the signal in F 3 with the inclusion of the CEDM term. The data in the top row are for p = (1, 0, 0) × 2π/La in the following order: insertion on u and d quarks for the a12m310 ensemble followed by insertion on u and d quarks for the a09m310 ensemble. The plots in the second row are in the same order except with p = (2, 0, 0) × 2π/La.  To conclude, the data so far show that an acceptable signal-to-noise ratio can be obtained in the connected contributions to F 3 for both the CEDM and γ 5 operator insertion. We also find that The ratio of the connected contribution of the γ 5 to CEDM term is a constant. This implies that the mixing of the γ 5 term with the CEDM term can be cast as part of the multiplicative renormalization of the CEDM operator. Controlling its O(1/a 2 ) UV divergent coefficient is a hurdle for lattice theories that respect chiral symmetry (domain wall or overlap fermions) and those that don't such as Wilsonclover fermions. Work to address the full mixing of the CEDM operator and the calculation of the disconnected diagrams and the reweighting factors is under progress.  Figure 7. The ratio of the connected contribution of the γ 5 term to the CEDM term in the F 3 form factor. The rest is the same as in Fig. 5.