Effective Field Theory in The Study of Long Range Nuclear Parity Violation on Lattice

A non-zero signal $A_\gamma^\mathrm{np}=(-3.0\pm1.4\pm0.2)\times 10^{-8}$ of the gamma-ray asymmetry in the neutron-proton capture was recently reported by the NPDGamma Collaboration which provides the first determination of the $\Delta I=1$ parity-odd pion-nucleon coupling constant $h_\pi^1=(2.6\pm 1.2\pm 0.2)\times 10^{-7}$. The ability to reproduce this value from first principles serves as a direct test of our current understanding of the interplay between the strong and weak interaction at low energy. To motivate new lattice studies of $h_\pi^1$, we review the current status of the theoretical understanding of this coupling, which includes our recent work that relates it to a nucleon mass-splitting by a soft-pion theorem. We further investigate the possibility of calculating the mass-splitting on the lattice by providing effective field theory parameterizations of all the involved quark contraction diagrams. We show that the lattice calculations of the easier connected diagrams will provide information of the chiral logarithms in the much harder quark loop diagrams and thus help in the chiral extrapolation of the latter.


I. INTRODUCTION
The study of parity (P)-violation in nuclear and atomic systems has continued to be a central topic in the low-energy community despite that the P-violation in Standard Model (SM) electroweak (EW) sector is well-established and all the EW parameters are already quite precisely measured. The reason is that we are really using the hadronic weak interaction (HWI) as a tool to understand the peculiarities in the strong interaction dynamics.
The non-perturbative nature of Quantum Chromodynamics (QCD) in the confinement region resembles a black box that asserts non-predictable dressings to the confined quarks in a hadron. Therefore, in order to examine its properties, the "bare" weak interaction which is well-understood serves as a probe inserted into the black box which then returns the HWI that is experimentally measured. The role of P-violation in this procedure is also obvious: as the effective strong interaction coupling is 10 6 times larger than the weak coupling, one relies entirely on a symmetry-violating signal to disentangle the HWI from the huge strong interaction background. It is therefore not the discovery of a non-zero P-violation signal in HWI, but the precise measurement of its value that will provide us with the opportunity of testing the SM with the interplay between weak and strong interactions.
Effects of the hadronic parity violation (HPV) are usually classified according to their isospin, and among all others the ∆I = 1 HPV possesses a special role as a unique probe of the hadronic weak neutral current. Moreover, it is the only channel that allows for a single pion-exchange, and hence plays a dominant role in the long-range HPV. Also, the ∆I = 1 P-odd pion-nucleon coupling plays a nontrivial role in the pp scattering through two-pion exchange as discussed in Refs. [1,2]. The recent observation of a P-violating 2.2 MeV gamma-ray asymmetry in the polarized neutron capture on hydrogen by the NPDGamma Collaboration [3] provides the first solid experimental confirmation of the isovector HPV, and is promised to create a new stir to the field that has been suffering from a "slow pace of (experimental) results since 1980" [4]. It is therefore timely to review our current knowledge of HPV and discuss how it could be improved by making the fullest use of the new experimental result.
Early attempts to describe HPV at the phenomenological level are based on isospin symmetry and perturbative expansions of small interaction energies, a strategy that is now inherited by the effective field theory (EFT) approach. A well-known example of such a kind is the work by Danilov [5] in the 60s that parameterized the P-odd nucleon-nucleon interaction at very low energy in terms of five S-P transition amplitudes with ∆I = 0, 1, 2. The ground-breaking work by Desplanques, Donoghue and Holstein (DDH) [6] in the early 80s adopted a very different starting point, namely to describe HPV through single exchange of light mesons π, ρ and ω with seven independent nucleon-meson coupling constants. Despite being a model, its succinctness has attracted much attentions and has become the basis of many experimental analysis. The development of the EFT description of HPV [7][8][9][10][11][12][13][14] signifies a switch to a model-independent framework that features pion-exchanges and contact terms, where a systematic power expansion with respect to a typical small momentum scale p ensures the finiteness of the number of operators needed in any given order. Translation tables, sometimes known informally as the "Rosetta stone" [4,15], are available to connect these many different effective descriptions of the same physics [2] (where the cutoff dependence is also discussed for the translation). Finally, nuclear model calculations have been carried out to connect the HPV coupling strengths to the experimental observables in nuclear or atomic systems; examples in the ∆I = 1 channel include Refs. [12,16,17], It is apparent that none of the frameworks above allows by itself a quantitative connection to the SM EW sector. The latter requires predictions of the theory parameters (such as the Danilov amplitudes, the nucleon-meson couplings in the DDH formalism and the low-energy constants (LECs) in the EFTs) in term of the SM EW parameters, which in turn require a precise control of strong dynamics in the confinement regime. In the original DDH paper, the SU(6) quark model was used to predict a "reasonable range" (0−11)×10 −7 and the "best guess" around 4.6×10 −7 for the ∆I = 1 P-odd pion-nucleon coupling h 1 π . Subsequent efforts include the use of quark models [18][19][20], Skyrme models [21][22][23] and QCD sum rules [24,25].
Their predictions of h 1 π , together with the NPDGamma outcome, which could be improved over using the analysis of Ref. [26], are summarized in Table I, and one can see that there is in general no agreement between different model predictions. Recent analyses based on large-N c [27][28][29] suggest a suppression of h 1 π from the naïve dimensional analysis result, and Ref. [15] uses this idea to identify a hierarchical structure of the five Danilov amplitudes.
Lattice QCD is currently the only available approach to compute low-energy hadronic observables from the first principle with a controlled error. Unfortunately, in contrast to the steady progress made in the lattice calculation of ∆I = 2 P-odd amplitudes [30,31], there is so-far only one very preliminary study of h 1 π by Wasem in Ref. [32] with no follow- unquantified assumptions as pointed out in Ref. [31]: (1) the three-quark representation of the N π interpolator; (2) the negligence of the so-called "quark loop diagrams"; (3) the calculation was done with only a single choice of volume, lattice spacing and pion mass; and (4) the lattice renormalization was not performed. We find the current situation not totally satisfactory because although the lattice calculation in the ∆I = 2 channel is technically simpler, there is no existing HPV experiment to our knowledge that depends only on the ∆I = 2 couplings (see, e.g. Ref. [4] for a summary) so that its comparison with experiments will not be straightforward. In contrast, a successful calculation of ∆I = 1 HPV can be directly confronted to the NPDGamma result. Therefore, despite all the technical difficulties, we believe a renewed lattice study of h 1 π is extremely worthwhile, and in this work we discuss how the proper application of a chiral EFT in the continuum space may help in alleviating part, if not all, of such difficulties.
The contents of this paper are as follows. We first introduce the theoretical basis of the ∆I = 1 HPV, including the underlying four-quark operators, their Wilson coefficients and the rigorous definition of the coupling h 1 π as a soft-pion matrix element. Next, we review the soft-pion theorem derived in our previous work [33] and present some of the technical details not included in that Letter. Then, we begin the analysis of contraction diagrams by rigorously defining them in terms of three-point correlation functions. With the aid of the partially-quenched chiral perturbation theory (PQChPT), we derive the theoretical expression for each contraction diagram that contributes to h 1 π as a function of the pion mass; such expressions are useful in performing chiral extrapolations from unphysical light quark masses to the physical ones. We point out that there are only a small number of LECs needed to fix the matrix elements, and provide approximate relations between different LECs that may facilitate their global fit. Finally, we briefly discuss the four-quark operators with strange quark fields and draw our conclusions.

II. THEORETICAL BASIS
We start by reviewing the electroweak interaction Lagrangian of the first two generations of quarks in the SM, where the electromagnetic, charged weak and neutral weak currents are defined as Here, ψ = (c u d s) T is the quark fields while the matrices {Q, C + , C 3 } are defined as where θ C is the Cabibbo angle. A single exchange of a W or Z boson leads to a P-odd interaction between a pair of quarks. At the energy scale E m W , m Z , the W or Z propagator shrinks to a point, so we obtain effective four-quark interactions involving the product of two weak currents.
In this work we focus on the ∆I = 1 P-violation in nucleon-nucleon interactions, and one may deduce from Eqs. (2) and (3) that they are dominated by neutral current interactions.
An easy way to understand this is to realize that in the θ C → 0 limit, the first and the second generations of quarks completely decouple in the current level, and the charged weak current involving light quarks then reads J µ W =ūγ µ (1 − γ 5 )d which is purely an isovector. Therefore, the symmetric combination J µ † W J W,µ can only form ∆I = 0, 2 objects but not ∆I = 1. In reality, the Cabibbo angle is not zero but the charged weak current contribution is suppressed by sin 2 θ C ≈ 0.05 so the neutral current contribution is still dominant. This is an important observation as it identifies the ∆I = 1 HWI as one of the very few direct experimental probes of the quark-quark neutral current effects at low energy.
Perturbative QCD modifies the structure of quark-quark weak interactions and introduces operators that do not appear in the original current-current product. Such an effect can be implemented by the QCD renormalization group (RG) running of the Wilson coefficients of the four-quark operators from the EW scale to the hadronic scale. At low energy, the the lightest hadronic DOFs following the spirit of EFT. The longest-range P-odd nuclear potential always consists of the pion exchange which, according to the Barton's theorem [36], only survives in the ∆I = 1 channel. Thus, the same Lagrangian L w PV can be expressed at low energy as where N = (p n) T is nucleon isospin doublet and the ellipses denote the remaining HPV interactions of shorter range. Eq. (7) may serve as a definition of the P-odd pion-nucleon coupling constant h 1 π , but we could equivalently express the latter in terms of a soft-pion matrix element, where m N is the averaged nucleon mass. The relation above can be obtained by taking the nπ + | ... |p matrix element at both sides of Eq. (7) and approximating the nucleon spinor product byū n u p ≈ 2m N , neglecting the small neutron-proton mass splitting and any small momentum transfer. Eq. (8) serves as the starting point for any first-principle or model-based calculation of h 1 π .

III. FROM P-ODD TO P-EVEN MATRIX ELEMENT
This section mainly serves as a review of the results in our previous work [33] with some more technical details added.

A. PCAC relation
The matrix element in Eq. (8) involves a soft pion in the final state that greatly complicates its analysis. We shall illustrate this point by considering a possible lattice QCD calculation of such a matrix element. First, one needs to choose a form of the interpolator for the nπ + state. The most natural choice with the largest overlap with the physical state is obviously a five-quark interpolator, for example, ε abc d a (u bT Cγ 5 d c )d e γ 5 u e . Such a choice will however lead to many contraction diagrams, some of which involving up quark propagators between n and π + are noisy and expensive. Another possible choice is a three-quark interpolator with negative parity, as adopted in Ref. [32], ε abc γ 5 u a (d bT Cγ 5 u c ). This choice avoids the calculations of the nπ + contraction diagrams, but cannot avoid a large overlap with single-nucleon excited states, e.g., the N (1535). Thus, the use of a three-quark interpolator would be unjustified without properly taking into account the excited-state contaminations.
Next, the rescattering effect between the final-state nπ + modifies the finite-volume correction on lattice from an exponentially-suppressed effect to a power-suppressed effect. Finally, while we want the final-state pion to have a vanishing momentum squared, lattice QCD only computes matrix elements of on-shell states. As a result, the lattice calculation returns not just h 1 π but its linear combination with the LECs of total-derivative operators that must be introduced to compensate the energy difference between the initial p and the final nπ + [37]. Although the leading effect can be canceled by considering the difference between p → nπ + and n → pπ − , but the m q -suppressed terms still retain, and they are in principle indistinguishable from the m q -dependent terms of h 1 π and lead to a sizable systematic error. The situation above can be greatly improved by a simple observation that the pion in the external state plays a special role in QCD: it is the pseudo-Nambu-Goldstone (pNG) boson that arises due to the spontaneously-broken chiral symmetry SU(2) R × SU(2) L → SU(2) V . That is, the axial charge operatorsQ i A do not annihilate the vacuum but produce soft pion states. Consequently, any matrix element involving an external soft pion can be replaced by another matrix element without the soft pion via the partially-conserved axial current (PCAC) relation, where F π = 92.1 MeV is the pion decay constant, and a, b are hadrons. This applies exactly to our case: Instead of computing nπ + | L w PV |p , one may compute n| [L w PV ,Q − A ] |p which is much simpler. To that end, it is beneficial to introduce the following four-quark operators, Explicitly evaluating the commutators in the PCAC relation yields where {θ q(j) , θ (s) q(j) } are defined as the P-even four-quark operators {θ q , θ (s) q } in Eq. (10) with the replacement τ 3 → τ j . Thus, the P-odd hadronic matrix elements nπ + | θ q |p and nπ + | θ (s) q |p can be mapped to the P-even hadronic matrix elements n| θ q(j) |p and n| θ (s) q(j) |p by PCAC, respectively. We may go one step further by transforming the latter into flavor-diagonal hadronic matrix elements through isospin rotation. By doing so, the operators in the matrix elements turn into those in Eq. (10). The final result is where L w PC is an auxiliary P-even Lagrangian, The reader should be alerted that L w PC is not the actual P-conserving weak four-quark interaction in SM, but simply an auxiliary Lagrangian introduced to facilitate the calculation of h 1 π , and the Wilson coefficients have to be identical with those in L w PV . At this level we have successfully mapped a P-odd N → N π matrix element to a flavordiagonal, P-even N → N matrix element. The right hand side (RHS) of Eq. (12) can be rewritten in terms of the neutron-pion mass splitting (δm N ) 4q ≡ (m n − m p ) 4q induced by the P-even Lagrangian L w PC using the Feynman-Hellmann theorem, Thus, combining Eqs. (8), (12) and (14), we obtain an approximate relation between h 1 π and (δm N ) 4q , which is one of the central results in Ref. [33].
We would like to point out that the idea above is not at all new. To our knowledge, the first application of PCAC in the study of the ∆S = 0 weak pion-baryon vertex appeared in Ref. [38] in the late 60s; it was also adopted in the DDH paper [6] as well as Ref. [21] as a starting point of their model-based estimation of h 1 π . The originality of Ref. [33] is really not in its application of PCAC, but rather in its quantitative analysis of the higher-order corrections which determines the degree of accuracy of the PCAC result, as we shall describe later.

B. Chiral perturbation theory analysis
The PCAC relation in Eq. (9) holds rigorously only in the exact chiral limit, i.e., when m π = 0. For example, it predicts that the matrix element at the left hand side should vanish ifÔ is chirally-invariant, which is obviously incorrect. Here we shall provide an immediate counter-example in a closely-related problem, namely the study of the P, T-odd pion-nucleon couplingḡ i π induced by higher-dimensional operators. Eq. (9) suggests that chirally-invariant operators such as the Weinberg three-gluon operator f ABCGA µν G Bν ρ G Cρµ would not contribute toḡ i π ; however, we know in reality that the contribution of such an operator is non-zero, but just suppressed by powers of m π [39]. Therefore, a truly practical application of the PCAC relation will need to take into account all the m π -related corrections to the level of desired precision. The above-mentioned task is made possible by recasting the PCAC statement in the language of chiral perturbation theory (ChPT), where Eq. (9) then becomes a simple consequence of two observables sharing the same LEC at the tree level. Higher-order corrections such as loop diagrams and counterterms to the left and the right sides of the equation can be computed order-by-order; any mismatch will then signifies a quantifiable violation of the tree-level matching. This idea was born of the in-depth studies of the P, T-odd pion-nucleon couplingḡ i π [40][41][42][43][44][45], and Ref. [33] constitutes its first implementation in HPV. Below we shall describe the method in detail.
Let us start by introducing the basic ingredients in a two-flavor ChPT with nucleons and pions. For that purpose it is instructive to first look at the QCD Lagrangian with two quark flavors, where q R,L = (1/2)(1 ± γ 5 )q is the right/left-handed component of the quark field. The quark mass matrix is given by M q = diag(m u , m d ). Direct inspection of Eq. (16) shows that the Lagrangian is invariant under the following SU(2) R × SU(2) L chiral rotation, in the limit of vanishing quark masses. Furthermore, if one would assume a transformation rule M q → RM q L † of the quark mass matrix, then L QCD would be chirally-invariant even with the existence of the quark masses. This is the so-called spurion trick to take into account symmetry breaking terms.
ChPT involves writing down all possible operators with hadronic DOFs that are consistent with the symmetry of L QCD under chiral rotation. There are infinitely many terms of such a kind. Thus, they have to be arranged according to a power counting scheme such that in any given order there are only a finite number of terms. The pions are contained in the matrix U defined as where F 0 is the pion decay constant in the chiral limit. It transforms as U → RU L † under the chiral rotation. The chiral Lagrangian of pions at LO consists of only two terms, In particular, the second term gives rise to the pion mass at LO, m 2 π = B 0 (m u + m d ). In the baryon sector, the nucleon doublet N appears as a matter field and can be chosen The chiral Lagrangian of nucleon at LO reads where m 0 and g 0 are the nucleon mass and axial coupling constant in the chiral limit, the chiral covariant derivative is defined as D µ = ∂ µ + Γ µ , and are the vector connection and axial vector, respectively, in the absence of external fields.
The naïve application of the nucleon Lagrangian in Eq. (21) will cause problems as it contains a large bare nucleon mass term that needs to be treated as a large energy scale. There are different schemes introduced to tackle this issue, for example the heavy baryon (HB) expansion [46][47][48], the infrared regularization [49,50] and the extended onmass-shall scheme (EOMS) [51,52]. While the first approach is technically simplest, the other two approaches possess extra advantages in preserving the analytic structures of the amplitude. In this work we are only interested in static matrix elements that are insensitive to the analytic behaviors around the threshold, so we shall just adopt the simplest HB approach. Below we shall briefly summarize its most important results, and interested readers may refer to standard textbooks such as Ref. [53] for details. In this approach, a redefinition of the nucleon field is performed to remove the large bare mass term in the Lagrangian, Eq. (21). As a consequence, the nucleon field N is effectively replaced by its "light" component N v which appears as a massless excitation, and the Dirac structures are effectively reduced as γ µ → v µ and γ µ γ 5 → 2S µ where v µ is a constant four-velocity vector and S µ = iγ 5 σ µν v ν /2 is the nucleon spin-matrix satisfying S · v = 0. With this, the LO nucleon Lagrangian becomes and the bare nucleon propagator reads i/ The effects of the ∆I = 1 four-quark operators can be most easily implemented to the chiral Lagrangian by adding the P-odd and P-even Lagrangian to obtain whereθ 1 = 2q a γ µ q aqbR γ µ τ 3 q bR ,θ 2 = 2q a γ µ q bqbR γ µ τ 3 q aR , One immediately observes that all the operators listed above would be invariant under SU(2) chiral rotation if the matrix τ 3 would transform as τ 3 → Rτ 3 R † . Therefore, the effect of the combined operators can be implemented to the chiral Lagrangian through a single Hermitian, traceless spurion X R = u † τ 3 u that transforms as X R → KX R K † . There is only one available operator at LO, where the quantity α is an unknown LEC that describes the strength of the dressed weak interaction. Upon expanding the Lagrangian with respect to the pion fields, we find that the first and the second terms correspond exactly to the neutron-proton mass splitting and the P-odd pion-nucleon coupling, respectively. Since they share the same unknown coefficient α, we immediately obtain Eq. (15) which is the PCAC prediction.
Our conclusion above is based on a tree-level analysis of ChPT, so it is legitimate to ask whether it still holds under higher-order corrections, including both the long-range (one loop) and short-range (counterterm) corrections. The way to proceed is to compute explicitly the higher-order corrections to both the left and the right sides of Eq. (15). First, the oneparticle-irreducible (1PI) diagram corrections to h 1 π and (δm N ) 4q at one loop are depicted by the first four diagrams and the last two diagrams in Fig. 1, respectively. They give where I a and I e are two standard loop functions (see, e.g., Ref. [43]) that often occur in one-loop calculations with HBChPT, with λ = 2/(4 − d) − γ E + ln 4π + 1. Notice that dimensional regularization is used to regularize the ultraviolet (UV) divergences in the loops.
Next, we consider the corrections from the LECs of order O(m 2 π /Λ 2 χ ). There are only two independent operators one could write down at this order, and their contributions read wherem = (m u + m d )/2 is the average light quark mass. 2 Finally, one also needs to include the standard pion [54] and nucleon [47,48] wavefunction renormalization as well as the higher-order correction to the pion decay constant F π , 3 where l 4 and B 20 are LECs in the mesonic ChPT and HB ChPT introduced in Ref. [54] and Ref. [55], respectively.
After grouping everything together one obtains [33] δ F π h 1 Namely, the higher-order corrections to F π h 1 π and (δm N ) 4q (up to O(m 2 π /Λ 2 χ )) satisfies exactly the same matching relation in Eq. (15). Hence, the accuracy of such a relation is better than m 2 π /Λ 2 χ ∼ 1% when the pion mass takes its physical value. This is the most important observation in Ref. [33]. It shows that the PCAC relation does not simply serve for an illustration purpose but is also quantitatively accurate. Hence, for all practical purposes, it is sufficient to calculate (δm N ) 4q or equivalently p| L w PC |p on lattice which is much easier 2 Notice that no isospin symmetry is assumed here. 3 We take this opportunity to correct a typo in Eq. (13) of Ref. [33]: The denominator in the first term of than h 1 π . We shall also contrast the result here with similar studies ofḡ i π . For the latter, the degree of accuracy of the tree-level matching depends critically on the involved P, T-odd operators: For the QCD θ-term and the quark chromo-electric dipole moments, the relation is preserved by one-loop corrections but violated by counterterms; for P, T-odd four-quark operators, the relation is violated by both one-loop corrections and counterterms [43,44]. A reason of such worse behavior is that a P, T-odd source may introduce a linear term in the pion field which has to be rotated away. Such a rotation gives rise to an extra piece in the tree-level matching relation in addition to the PCAC prediction, and the matching of this extra piece is usually not preserved at higher orders.

ANALYSIS
The objective of this paper is to investigate the possibility of performing a high-precision calculation of the coupling h 1 π on the lattice. To that end, it is instructive to go through the existing limitations of the original Wasem calculation (as pointed out in Ref. [31]) and ask ourselves how many of them can be properly taken into account from a theoretical point of view. Clearly, our strategy in Ref. [33] avoids the unjustified application of the N π interpolator and alleviates the effect of the finite-volume correction, but does not resolve the other issues and therefore is not the end of the story. We shall devote the rest of this paper to the discussion of the remaining problems that could be at least addressed partially in a continuous field theory. In particular, we shall discuss the properties of different Wick contraction diagrams that occur in the the calculation of (δm N ) 4q on lattice and the chiral extrapolation of the lattice result. We shall assume isospin symmetry throughout the rest of the paper.
We start by reviewing the standard procedure to compute a hadronic matrix element of a four-quark operator on lattice. Consider the following ∆I = 1, P-even four-quark operator without strange quark fields:Ô where Γ µ is a general Dirac structure. To be concrete we choose the color contraction to be of the type aabb but the discussion of the operators with color contraction of the type abba proceeds exactly the same way. We are interested in the matrix element ofÔ with respect to the static proton state. Following the spirit of lattice QCD, we define a two-point and a three-point correlation functions (with τ > τ 0 > τ ): whereÔ p is the so-called proton interpolator which is an operator with the same quantum number as the proton state, x, x 0 , x are some fixed spatial points on the lattice and τ, τ 0 , τ are Euclidean time: τ = it. One may insert a compete set of states between any two operators, with L the lattice size. Here, the summation also runs over all possible quantized momentum modes. 4 When τ − τ 0 → +∞ and τ 0 − τ → +∞, only the static proton contribution that scales as exp{−m N (τ − τ )} survives. Therefore one obtains In actual lattice calculations, one may choose to sum over a part of the spatial points { x, x 0 , x } in order to improve the signal; this is however not directly relevant to our study of contraction diagrams so we shall not discuss it any further.
Using the Wick theorem, the three-point correlation function F 3 (τ, τ 0 , τ ) can be separated into different terms according to the types of quark contractions. To rigorously define a specific contraction as a mathematical quantity, we need to first choose the three-quark interpolating operator of the ground-state spin-1/2 baryons. Here we find the most convenient choice for theoretical manipulation to be where C is the charge conjugation matrix. It satisfies the following symmetry relations under the exchange of quark flavors, χ(q 1 , q 3 , q 2 ; x, τ ) = χ(q 1 , q 2 , q 3 ; x, τ ), A contraction function X(τ, τ 0 , τ ) may now be defined as the contribution to the following general three-point function, where {i, i , j, j } are flavor indices, from a definite contraction diagram X. With the exchange symmetry of χ as shown in Eq. (38), it is straightforward to demonstrate that there are altogether eight types of independent contraction diagrams as depicted in Fig. 2. In each diagram, the ellipse at the left representsχ(q 1 , q 2 , q 3 ; x , τ ), that at the right represents χ(q 1 , q 2 , q 3 ; x, τ ), whereas a solid line represents the contraction between a pair of quark fields and the short dashed line denotes the four-quark vertex. We shall name the first two diagrams as "connected", the next four as "quark loop" and the last two as "vacuum" diagrams, respectively. One may read off the explicit form of each contraction function X(τ, τ 0 , τ ) from Fig. 2. For instance, where we have suppressed the spacetime coordinates: ( x, τ ) for the first three quark fields, ( x 0 , τ 0 ) for the next four inside the square brackets and ( x , τ ) for the last three. Notice also that in such a definition the contraction function X(τ, τ 0 , τ ) includes all possible minus signs due to the switching of positions between quark fields. For the vacuum contractions V 1 (τ, τ 0 , τ ) and V 2 (τ, τ 0 , τ ), we choose to normalize them according to the two-point function F 2 (τ, τ ) which is independent of the choice of baryon in the flavor degenerate limit. With these, we may now define the quantity M X as We will show that the limits exist in a short while. The matrix element of any four-quark operator with respect to any external spin-1/2 static baryon states can now be expressed as a linear combination of {M X }. For instance, This is how one could systematically express a four-quark matrix element in terms of contributions from different contraction diagrams on lattice.
In practice, the lattice computations of different types of X(τ, τ 0 , τ ) contractions will involve very different techniques. While the C i (τ, τ 0 , τ ) can usually be computed quite economically with high precision, the quark loop contractions D ai (τ, τ 0 , τ ) and D bi (τ, τ 0 , τ ) contain a quark propagator that starts and ends at the same spacetime point and is extremely noisy. For such propagators, one needs to average its value over all lattice points in order to improve its signal, but then it requires the use of all-to-all propagators that are computationally expensive [57][58][59][60][61][62][63][64][65][66][67][68][69][70][71]. As a consequence, with a given computational power, The net effect of this extension is that one introduces quark flavors that can only appear in external states but not in loops, and the strong interaction theory of such system is known as the partially-quenched QCD (PQQCD) [72,73]. Within this framework, any contraction diagram of interest can be constructed by appropriately choosing the quark contents in either the external states or the operators [74]; such an idea has been previously applied in studies of the hadronic vacuum polarization [75], the pion scalar form factor [76] and the ππ scattering amplitudes [77]. In our case, one could easily demonstrate that each quantity M X can be written as a linear combination of four-quark matrix elements in PQQCD; an explicit example is given in Appendix B. With that we also show that each contraction function X(τ, τ 0 , τ ) has the correct asymptotic exponential behavior of exp{−m N (τ − τ )} that guarantees the existence of the limits in Eq. (41).

V. PQCHPT ANALYSIS
In the previous section we have successfully separated each contraction diagram into well-defined matrix elements in PQQCD, and here we shall proceed to study the low-energy behavior of each individual contraction that contributes to the four-quark matrix element p|Ô |p of our interest. This involves the application of the low-energy EFT of PQQCD as follows.
In the massless limit, PQQCD with four fermionic and two ghost quarks has a "graded" SU(4|2) chiral symmetry, namely the Lagrangian is invariant under the transformation where R ∈ SU(4|2) R , L ∈ SU(4|2) L are elements of a special unitary (4|2) graded symmetry group. At low energy, this graded chiral symmetry is spontaneously broken as SU(4|2) L × SU(4|2) R → SU(4|2) V which generates 35 pNG particles in complete analogy to the ordinary chiral symmetry breaking of two-flavor QCD: SU(2) L ×SU(2) R → SU(2) V . One may proceed to write down the PQChPT [74,[78][79][80] which is the low-energy EFT pf PQQCD that incorporates the interactions between the pNG particles and other matter fields.
The essential information of heavy baryon PQChPT can be found in Appendix C and the references listed above, so here we shall concentrate on the implementation of the weak interaction in the chiral Lagrangian. Let us consider a generic four-quark operator in SU(4|2), where τ A , τ B are SU(4|2) generators; we shall assumeÔ 4q to be Hermitian for simplicity. One may decompose the quark fields inÔ 4q into left-and right-handed components and obtain It is now obvious that the implementation ofÔ 4q in the chiral Lagrangian will involve which is a straightforward generalization of the spurion introduced in Section III. With this, there are eight independent operators in the baryon sector one could write down at LO, where The corresponding Lagrangian can be written as Here we choose to label a given operator according to the name of the contraction diagram in Fig. 2 that contracts the quark indices in the same way as the operator. We observe from Eqs. (38) and (C5) Therefore, a combination of Eqs. (42) and (49)  and M D b1 + M D b2 in terms of the diagonal matrix element ofÔ with respect to the baryon states p,Σ 0 andΛ. One may then compute the latter to one loop using HB PQChPT to obtain the former. Below we summarize the main results. First, at tree-level, we obtain where we have defined for compactness.
Next we consider the one-loop contributions which can be further sub-divided into 1PI contributions and the wavefunction renormalization. There are two types of 1PI diagrams with topologies identical to those of Fig. 1 (e) and (f), and their contributions are labeled as 1PI(e) and 1PI(f), respectively. The results are and Finally, the wavefunction renormalization Z N is identical to that in the ordinary two-flavor ChPT result in Eq. (31) and it applies to all fermionic baryons due to the exact SU(4|2) V symmetry.
The one-loop calculations above contain UV divergences which need to be canceled by counterterms of order O(m 2 π /Λ 2 χ ). Here we shall not bother to write down the full counterterm Lagrangian because there are too many available terms, and eventually what we only need to know is that M C 1 , M D a1 + M D b1 and M D a2 + M D b2 will acquire independent combinations of counterterms. Therefore, by collecting tree-level, one-loop and counterterm contributions, we obtain Here m N is the physical nucleon mass that has the full quark-mass dependence, and {δα C 1 , δα D 1 , δα D 2 } represent the total counterterm contribution to each respective quantity.
Let us try to understand the results above. The first important observation is that the two LECs {α C 1 , α C 2 }, which contribute only to connected diagrams at tree level, enter the quark loop diagrams in the form of chiral logarithms; this feature can be understood diagrammatically as depicted in Fig. 3. On the other hand, there is no way that the LECs {α D i } can induce connected diagrams through loop corrections, and therefore we observe that both the tree-level and chiral logarithms of M C 1 depend only on {α C i }. Since connected diagrams can be readily computed on lattice, they may be computed with several values of m π which then, using our derived formula, allow for a determination of the LECs {α C i }. By doing so, we do not just acquire the full information of the connected diagrams, but also fix a part of the chiral logarithms in the quark loop diagrams. Even though the leading terms and the remaining chiral logarithms of the latter can only be fixed by direct lattice calculations of such diagrams, now they depend on a smaller amount of unknown LECs (i.e., just α D 1 and α D 2 ), making their chiral extrapolation much easier.
We should end this section by mentioning another technical detail. From Eq. (54) we find that the calculation of M C 1 can only fix the combination 4α C 1 + α C 2 but we do need to know the two LECs separately. Therefore, we shall supply also the theoretical formula for the other contraction diagram M C 2 . The easiest way to calculate it is to realize that if we and thus we only need to compute one more matrix element Σ 0 |Ô |Σ 0 . We simply quote the final result, Thus, by computing M C 1 and M C 2 on lattice one is able to fix α C 1 and α C 2 simultaneously.

VI. SPIN-FLAVOR SYMMETRY
In Section V, we show that in order to determine the tree-level and chiral logarithmic terms in p|Ô |p it is necessary to perform a global fit of the pion mass dependence using four LECs: It would be beneficial to obtain approximate relations among the LECs, especially between the more difficult α D 1 and α D 2 in order to provide a guidance to the starting point of the global fit. To facilitate such a discussion, we start by providing another matrix representation of the combination M D a1 + M D b1 and M D a2 + M D b2 .
If we choose the four-quark operator, then we have Approximate relations between the matrix elements above can then be obtained by considering the spin-flavor symmetry [81] among baryons. In a QCD with N fermionic quark flavors, such a symmetry means that the quarks {q i } of a definite flavor q and spin i form the fundamental representation of a SU(2N ) group and all hadrons can be grouped into irreducible representations of that symmetry group. For instance, the baryon octet and decuplet collectively form an irreducible 56-plet of the spin-flavor SU (6). It is well-known that the spin-flavor symmetry is a direct consequence of the large-N c limit (with N c the number of colors), and for most of the practical purposes it is simply equivalent to the non-relativistic quark model [82]. Therefore, for the discussion here let us consider the spin-flavor wave function of a spin-upΣ 0 state in the quark-model (QM) representation, where the quark creation and annihilation operators satisfy the anti-commutation relation and the baryon state is normalized as QM Σ 0 |Σ 0 QM = 1. With this, we can compute the matrix element ofÔ uu andÔ kk with respect to the spin-upΣ 0 state, The last equality is due to rotational symmetry. Similarly, one also obtains Therefore, combining Eqs. (61), (62) and the u ↔ k flavor symmetry, we arrive at QM Σ 0 ↑ |Ô uk |Σ 0 ↑ QM = 0 which implies the following approximate relation according to Eq. (58), as a consequence of the spin-flavor symmetry. At LO in HB PQChPT, this in turns implies an approximate relation to the LECs, The equation above may now serve as a starting point for the global fit of the four LECs i } in Eq. (10). They give rise to contraction diagrams of type D b , except that now the quark in the loop has a heavier mass. It is unavoidable that these quark loop diagrams must be calculated directly on lattice if we are to study their contributions to h 1 π . Nevertheless, the EFT analysis is still beneficial as it provides an extrapolation formula of the matrix element with respect to the light quark mass, which can be read off directly from Eq. (32): Notice that both the LO term and the chiral logarithms depend only on one single LEC α (s) i . One could therefore compute this matrix element on lattice with unphysical pion mass, which is presumably easier, and then extrapolate the result to the physical region using the formula above.

VIII. CONCLUSIONS
HPV has been studied for many years. Among others, the ∆I = 1 HPV holds a special role as a unique probe of hadronic neutral weak current as well as one of the main contributors of of long-range nuclear PV. Moreover, with the release of the NPDGamma result, the ∆I = 1 P-odd pion-nucleon coupling h 1 π is now the only DDH coupling with a definite isospin that has been numerically measured through a single experiment. Therefore, the first-principle calculations of h 1 π are highly desirable as they are directly comparable to experimental results.
Despite the above, currently we observe a lack of progress in the lattice study of the ∆I = 1 HPV comparing to its ∆I = 2 counterpart. The latter involves a direct computation of nucleon-nucleon scattering amplitudes with the insertion of ∆I = 2 four-quark operators. It does not require computations of noisy disconnected diagrams but the total amount of contractions is tremendous even in the exact isospin limit. On the other hand, although disconnected diagrams are unavoidable in the ∆I = 1 channel, the total amount of contractions is much less. Furthermore, given its more straightforward relation to experimental data, we believe that the study of the ∆I = 1 HPV on lattice should receive the same amount, if not more, of attention as the ∆I = 2 one. In this paper we investigate in some detail how a continuous EFT may help in the future lattice calculation of h 1 π . In Ref. [33] we show that h 1 π can be recast as a neutron-proton mass splitting induced by a set of ∆I = 1 P-even four-quark operators. Improving from the limitations of PCAC, we show that such a relation holds with a precision better than 1% by considering the longand short-range higher-order corrections. This observation turns the lattice study of h 1 π into computations of P-even three-point correlation functions involving only five sets (three independent combinations) of contractions. Two combinations among them are quark loop contractions which are in principle noisy, but in this work we show that one can obtain partial information of the chiral logarithms in the quark loop diagrams by studying the much easier connected diagrams. We further demonstrate that one only needs to perform a global fit with four independent LECs, two of which can be obtained easily from connected diagrams, in order to completely determine the LO and chiral logarithmic terms in h 1 π induced by nonstrange operators. Approximate relations among LECs based on the spin-flavor symmetry are also derived to facilitate the global fit. For operators with strange quark fields, fitting of one single LEC from the computation of one quark loop contraction is needed to determine the LO and chiral logarithmic terms. We hope that the analysis above will provide extra motivations for the lattice community to perform an up-to-date first-principle computation of h 1 π which will constitute a new breakthrough in our understanding of HWI.
are more commonly used. The problem is that they do not satisfy the exchange symmetry relations in Eq. (38) and are thus not sufficient to specify a spin-1/2 baryon. One may construct the linear combinations (i = 1, 2) that do satisfy Eq. (38), but the price is that now each interpolator contains two terms instead of one. In fact, using Fierz identity one is able to show that so these choices of interpolators are not all independent. Furthermore, it is demonstrated that the interpolator χ 2 has negligible overlap with the ground-state baryon [83]. The explanation is that χ 2 scales as O(p 2 /E 2 ) in the non-relativistic expansion, which implies that it overlaps more with excited states than with the ground state [84]. Therefore, as far as this work is concerned, one could legitimately replace χ(q 1 , q 2 , q 3 ) → χ 1 (q 1 , q 2 , q 3 ) without affecting any of the discussions above.
This also implies that each contraction function X(τ, τ 0 , τ ) is related to a linear combination of the "physical" matrix elements in SU(4|2) and thus possesses the correct asymptotic exponential behavior exp{−m N (τ − τ )}. Therefore, each M X can be obtained on lattice in the same way as how the full matrix element p|Ô |p is obtained.
An immediate application of Eq. (B1) is the calculation of all {M X } at tree level using the LO Lagrangian in Eq. (48). The results read which give clear meaning of each LEC in terms of contraction diagrams.

Appendix C: Essentials of HB PQChPT
In this appendix we summarize the basic results of HB PQChPT that are used in this paper.

Grading factor
In PQQCD there are both fermionic and bosonic quarks. To determine whether a pair of quark fields commute or anti-commute, it is convenient to define a quantity η i such that η i = 1(0) when i is a fermionic (bosonic) index. Then, any two quantities A and B are said to have grading factors η A and η B , respectively, if For instance, a quark field q i obviously has a grading factor η i so a pair of quark fields anticommute only when both of them are fermionic. Meanwhile, for a matrix Γ in the flavor space, its matrix element Γ ij in general has a grading factor of η i + η j .

The pNG particles
Similar to the ordinary ChPT, the pNG particles (we refrain from using the word "boson" because they can be either bosonic or fermionic though the spin is always 0) are contained in the matrix U . However, here it is more preferable to parameterize U as and to study the propagator of Φ ij . By doing so we avoid the need to define explicitly all the 9N 2 − 1 generators in a SU(2N |N ) PQChPT. With all quarks degenerate, the pNG propagator takes the following compact form [79], where ε i = (−1) η i +1 , and T ij equals −1 when both i, j are ghost indices and 1 otherwise.
The first term at the RHS of Eq. (C3) is contributed by the neutral particles while the second term is by the charged particles.

Three-index representation of spin-1/2 baryons
The spin-half baryons in PQChPT are usually represented by a three-index form B γ ijk , where {i, j, k} are flavor indices and γ is the Dirac index [85]. It can be most easily understood by comparing with a three-quark representation, Of course in ChPT one does not deal explicitly with quark fields, but Eq. (C4) is still useful in determining the symmetries and transformation rules of B ijk as follows.
• Symmetries under the exchange of two flavor indices: • Transformation of B ijk under the chiral rotation [85]: which can be understood by first going back to Eq. (C4), making the transformation q → Kq , and moving K jj , K kk through the quark fields to the left.
The barred quantity of the baryon field is defined asB kji ≡ (B ijk ). It obviously satisfies the exchange symmetry relations, and the transformation rulē

Constructing chiral invariants
Here we introduce all the remaining ingredients needed to construct the chirally-invariant Lagrangian in the baryon sector. The definitions of the vector connection Γ µ and the axial vector u µ (in the absence of external sources) The appearance of the grading factors can be easily understood as dictated by Eq. (C6).
For instance, (−1) (η j +η j )η i is required for (Γ µ ) jj to pass through the index i in B ij k in order to act on the index j .
As far as the one-loop analysis of contraction diagrams in this work is concerned, the only strong interaction Lagrangian we need is the SU(4|2) HB PQChPT Lagrangian at LO [85], One notices that there are two independent axial couplings ρ and β, while there is only one g 0 in the LO SU(2) HBChPT Lagrangian. However, there are two axial couplings D and F in the SU(3) version (see, e.g., Ref. [86]), where B v denotes the baryon octet in Eq. (C24) and u µ takes the SU (3) which reproduce the relations in Ref. [88] (g 0 is written as g A therein) ρ = 4 3 g 0 + 1 3 once we identify g 0 = D + F and g 1 = 2(F − D).

Independent baryon fields and the baryon propagator
The three-index baryon fields {B ijk } are not all independent due to the exchange symmetry relations in Eq. (C5). Thus, let us denote the independent baryon fields as {B a }, and {B ijk } can be expressed in terms of the independent fields as where {ψ a ijk } are c-numbers, and the complex conjugate is just a convention. In fact, the equation above defines the coefficients {ψ a ijk }, and some useful examples of these coefficients are summarized in Appendix C 7. As shown in Eq. (C11), B ijk is normalized such that B kji iv · DB ijk reproduces the properly-normalized kinetic terms of each independent baryon. This imposes the following orthonormal condition to the coefficients, The inversion of Eq. (C15) is not unique; however, the most convenient form of inversion is which is a direct consequence of the orthonormal condition.

Prescription for the correct usage of Feynman rules
We shall name the Feynman vertices directly extracted from the HB PQChPT Lagrangian as the "naïve vertices", which are in terms of the non-independent baryon fields {B ijk }. For instance, the Lagrangian L = λB kji Φ ii B i jk gives the following naïve vertex: iM(B ijk Φ pq → B lmn ) = iλδ iq δ lp δ jm δ kn . Similarly, the propagator in Eq. (C18) should be known as the "naïve propagator".
A direct application of naïve Feynman propagators and vertices in the computation of physical amplitudes is obviously inappropriate. However, there is a simple prescription that ensures the correctness of the final amplitude. For the simplicity of notation we shall 5 The expression of F 1/2 lmn,ijk in Ref. [85] contains a couple of typos which are corrected here.
where {α i , β j } are summed over but {a i , b j } are not. One can show that the contraction of ψ, ψ * to the external states is equivalent to taking only the independent baryonic DOFs in the loop calculation. A good thing about this relation is that we now need only to know the coefficients {ψ a ijk } for ordinary baryons that exist as external states, i.e., baryons consist of three fermionic (dynamical or valence) quarks. 7. {ψ a ijk } in SU (3) In a theory with only three dynamical quarks (q = u, d, s), the three-index representation B ijk can be directly mapped to the baryon octet as [85,87] where the baryon octet matrix B is defined as With this we can easily obtain the independent non-vanishing coefficients ψ a ijk in SU (3): All other non-vanishing coefficients can be obtained by symmetry: ψ a ijk = ψ a ikj and ψ a ijk = −ψ a jik − ψ a kji (no grading factors because all quarks are fermionic). The results above can also be used to determine the coefficients ψ for the other baryons in PQChPT. For example, the SU(4|2) baryon statesΣ 0 andΛ we introduced in Section V have the following independent non-vanishing coefficients: