Parity- and time-reversal-violating nuclear forces

Parity-violating and time-reversal conserving (PVTC) and parity-violating and time-reversal-violating (PVTV) forces in nuclei form only a tiny component of the total interaction between nucleons. The study of these tiny forces can nevertheless be of extreme interest because they allow to obtain information on fundamental symmetries using nuclear systems. The PVTC interaction derives from the weak interaction between the quarks inside nucleons and nuclei and the study of PVTC effects opens a window on the quark-quark weak interaction. The PVTV interaction is sensitive to more exotic interactions at the fundamental level, in particular to strong CP violation in the Standard Model Lagrangian, or even to exotic phenomena predicted in various beyond-the-Standard-Model scenarios. The presence of these interactions can be revealed either by studying various asymmetries in polarized scattering of nuclear systems, or by measuring the presence of non-vanishing permanent electric dipole moments of nucleons, nuclei and diamagnetic atoms and molecules. In this contribution, we review the derivation of the nuclear PVTC and PVTV interactions within various frameworks. We focus in particular on the application of chiral effective field theory, which allows for a more strict connection with the fundamental interactions at the quark level. We investigate PVTC and PVTV effects induced by these potential on several few-nucleon observables, such as the longitudinal asymmetry in proton-proton scattering and radiative neutron-proton capture, and the electric dipole momentsof the deuteron and the trinucleon system.


INTRODUCTION
The interaction between nucleons is at the heart of nuclear physics and has been a subject of great scientific interest since many decades. The strong nuclear forces have their origin in the residual interaction between quarks and gluons inside colorless nucleons and are described by quantum chromodynamics (QCD). The resulting parity-conserving, time-reversal-conserving (PCTC) nuclear interactions are known to exhibit a complicated pattern with a delicate interplay of strongly state-dependent repulsive and attractive pieces. While the nucleon-nucleon (NN) scattering data below the pion production threshold can nowadays be accurately described by modern NN potentials, the (weaker) three-nucleon (3N) forces and the electromagnetic interactions (EM) between the nucleons, known to play an important role in the nuclear structure and dynamics, are not so well understood and subject of active research. The current status of PCTC nuclear forces is reviewed in other contributions to this topical issue.
In addition to the bulk PCTC interactions mentioned above, nuclear forces also feature much tinier components, which originate from the weak forces between quarks and/or physics beyond the standard model (BSM) and whose strength is smaller than that of the strong and EM interactions by many orders of magnitude. These tiny components are, nevertheless, extremely interesting since investigation of their effects may shed new light on fundamental symmetries and BSM physics. While effects of such exotic PCTC components are, of course, completely overwhelmed by the strong and EM nuclear forces, parity-(P) violating and/or time-reversal-(T) violating nuclear interactions can be determined by measuring specific observables which would vanish if these symmetries were conserved. In this contribution, we review the theory of parity-violating, time-reversal-conserving (PVTC) and parity-violating, time-reversal-violating (PVTV) nuclear forces and discuss selected applications.
Starting from the fifties of the last century, a wide variety of phenomenological models have been developed to describe nuclear forces, the most prominent utilizing the one-boson exchange picture, see Ref. [1] and references therein. More recently, the development of chiral effective field theory (χEFT) [2] has given a new impetus to the derivation of nuclear interactions [3,4,5]. The χEFT approach utilizes the spontaneously broken approximate SU(2) L ×SU(2) R chiral symmetry of QCD 1 in order to describe low-energy dynamics of pions, the (pseudo-) Goldstone bosons of the spontaneously broken axial generators, in a systematic and model-independent fashion within the framework of the effective chiral Lagrangian [6,7,8,9,10,11], see Refs. [12,13,14] for review articles. Owing to the derivative nature of the Goldstone boson interactions, the scattering amplitude in the pion-and single-baryon sectors can be calculated via a perturbative expansion in powers of Q/Λ χ , where Q refers to momenta of the order of the pion mass m π and Λ χ ∼ m ρ ∼ 1 GeV denotes the chiral symmetry breaking scale, with m ρ the ρ-meson mass. The effective Lagrangian involves (an infinite number of) all possible hadronic interactions compatible with the symmetries of QCD, which are naturally organized according to the number of derivatives and/or quark or pion mass insertions. 2 Every term in the effective Lagrangian is multiplied with a coefficient, whose strength is not fixed by the symmetry. These so-called low-energy constants (LECs) can be determined by fits to experimental data and/or obtained from lattice QCD simulations, see [13,14] and references therein. At every order in the Q/Λ χ -expansion, only a finite number of terms from the effective Lagrangian contributes to the scattering amplitude. The resulting framework, commonly referred to as chiral perturbation theory (χPT), is nowadays widely applied to analyze low-energy processes in the Goldstone boson and single-nucleon sectors. It has also been 1 Here and in what follows, we restrict ourselves to the two-flavor case of the light up and down quarks unless specified otherwise. 2 In the isospin limit, the quark and pion masses are related to each other via m 2 π = 2Bmq + O(m 2 q ), where B is a constant proportional to the quark condensate 0|ūu|0 = 0|dd|0 . generalized to study few-and many-nucleon systems, where certain resummations beyond perturbation theory are necessary in order to dynamically generate the ultrasoft scale associated with nuclear binding. According to [2], the breakdown of the perturbative expansion for the NN scattering amplitude is traced back to enhanced contributions of ladder diagrams, i.e. Feynman diagrams that become infrared divergent in the static limit of infinitely heavy nucleons. The simplest and natural way to resum enhanced ladder diagrams is provided by solving the nuclear Schrödinger equation. The framework therefore essentially reduces to the conventional quantum mechanical A-body problem. The corresponding nuclear forces and current operators are defined in terms of non-iterative parts of the scattering amplitude, which are free from the above mentioned enhancement. They can be derived from the effective chiral Lagrangian in a systematically improvable way via a perturbative expansion in powers of Q/Λ χ [4,5]. Assuming the scaling of few-nucleon contact operators according to naive dimensional analysis 3 , the PCTC interactions are dominated by the pairwise NN force, which receives its dominant contribution at order (Q/Λ χ ) ν with ν = 0, defined to be the leading order (LO). Parity conservation forbids the appearance of nuclear forces at order ν = 1, so that the next-to-leading order (NLO) contribution to the PCTC NN potential appears at order ν = 2. Next-to-next-to-leading order (N 2 LO) has ν = 3 and so on. PCTC three-and four-nucleon forces are suppressed and start contributing at orders ν = 3 (N 2 LO) and ν = 4 (N 3 LO), respectively. Presently, the chiral expansion of the PCTC NN force has been pushed to order ν = 5 (N 4 LO) [16,17,18,19], while many-nucleon interactions have been worked out up through N 3 LO, see [4,5] and references therein. We further emphasize that a number of alternative formulations of χEFT for nuclear systems have been proposed [20,21,22,23,24,25], see also Refs. [26,27,28,29,30,31] for a related discussion.
Another framework to analyze nuclear systems at very low energies is based on the so-called pionless formulation of EFT, see Refs. [32,33,31] for review articles. It is valid at momenta well below the pion mass, at which the pionic degrees of freedom can be integrated out. In the resulting picture, nucleons interact with each other solely through short-range contact two-and many-body forces. This formulation is considerably simpler than χEFT both at the conceptual and practical levels, and has been successfully applied to study e.g. Efimov physics and universality in few-body systems near the unitary limit, lowenergy properties of halo-nuclei and reactions of astrophysical relevance, see Refs. [32,33,31] and references therein.
In this paper we focus on the PVTC and PVTV interactions in the frameworks of χEFT and pionless EFT. We also outline various meson-exchange models frequently adopted to analyze the results for some PVTC and PVTV observables. In the subsections below, we briefly discuss the origin of the PVTC and PVTV interactions and summarize the current experimental and theoretical status of research along these lines.

The PVTC interaction
The PVTC component of the nuclear force is governed by the weak interaction between the quarks inside the nucleons (and pions). Studying such effects, therefore, opens a window on the so-called "pure" hadronic weak interaction (HWI) [34,35,36,37,38]. This part of the weak interaction is far less known experimentally.
A number of experiments aimed at studying PVTC in low-energy processes involving few-nucleon systems have been completed/are being planned at cold-neutron facilities, such as the Los Alamos Neutron Science Center (LANSCE), the National Institute of Standards and Technology (NIST) Center for Neutron Research, the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory, and the European Spallation Source (ESS) in Lund. The primary objective of this experimental program is to determine the LECs which appear in the PVTC nuclear potentials. For a recent review of the current status of experiments along this line and the impact of anticipated results see Ref. [39].
PVTC nuclear forces have already been analyzed in the framework of χEFT [40,41,42]. The LO PVTC NN force is driven by the one-pion-exchange term with ν = −1, while the NLO terms with ν = 1 emerge from two-pion-exchange diagrams and NN contact interactions 4 . In Ref. [44], it was shown that the PVTC NN potential involves only five independent contact operators at this order corresponding to five S-P transition amplitudes at low energies [45]. Including the PVTC pion-nucleon coupling constant h 1 π , the NN potential at NLO thus contains six LECs which need to be determined from experimental data. At N 2 LO one has to take into account five additional LECs, which determine the strength of the subleading PVTC pion-nucleon interactions [46].
In pionless EFT, the LO PVTC NN potential is completely described in terms of the already mentioned five contact terms [36,47]. The large-N c scaling of PVTC NN contact interactions was analyzed in Refs. [48,47]. These studies suggest that three out of five PVTC contact interactions are suppressed by a factor of (1/N c ) 2 or by the factor sin 2 θ W ≈ 0.23, see also a related discussion in [49]. If the large-N c scaling persists to the physically relevant case of N c = 3, the pionless potential at LO should be dominated by only 2 LECs [39]. Unfortunately, the currently available experimental data do not allow one to draw definitive conclusions on whether the suggested large-N c hierarchy of PVTC contact interactions is indeed realized in Nature.
Regarding the various meson-exchange models developed to describe the PVTC interaction, we will mainly discuss the model proposed by Desplanques, Donoghue, and Holstein (DDH) [50] which includes pion and vector-meson exchanges with seven unknown meson-nucleon PVTC coupling constants.

The PVTV interaction
PVTV nuclear forces originate from more exotic sources at the fundamental level, which include the so-called θ-term in the Standard Model (SM) Lagrangian [51], or even BSM interactions [52]. Due to the CPT theorem, any PVTV interaction also violates the CP symmetry, where C refers to charge conjugation. CP violation is a key ingredient for the dynamical generation of a matter-antimatter asymmetry in the Universe [53]. The SM with three generations of quarks has a natural source of CP-violation in the phase of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. This mechanism is however not sufficient for explaining the observed asymmetry [54].
The phase of the CKM matrix also does not contribute sizably to the nuclear PVTV interaction. For example, let us consider the electric dipole moment (EDM) of a system of particles. A nonzero permanent EDM of a particle or a system of particles necessarily involves the breaking of the parity and time-reflection symmetries. EDMs of the electron, nucleons and nuclei are mostly sensitive to P-and T-violating flavor-diagonal interactions. To induce a non-zero EDM, on the other hand, the phase of the CKM requires contributions from all three generations of quarks, including heavy quarks, leading to a large suppression [52,55,56,57]. For example, the expected size of the nucleon EDM based on the CKM mechanism in the SM is |d CKM N | ∼ 10 −18 e fm [58,59]. Therefore, any observed permanent EDM of an atomic or nuclear system larger in magnitude than the expected size within the SM would highlight PVTV effects beyond the CKM mixing matrix. The present experimental upper bounds on the EDMs of neutron and proton are |d n | < 3.0 · 10 −13 e fm [60,61] and |d p | < 2.0 · 10 −12 e fm, where the proton EDM has been inferred from a measurement of the diamagnetic 199 Hg atom [62] using a calculation of the nuclear Schiff moment [63]. For the electron, the most recent upper bound is |d e | < 1.1 · 10 −16 e fm [64], derived from the EDM of the ThO molecule. In all cases, the current experimental sensitivities are orders of magnitude away from CKM predictions.
χEFT allows to derive PVTV nuclear forces in a systematic and model independent way. To this aim, the PCTC effective chiral Lagrangian has to be extended to include all possible PVTV terms classified according to their chiral dimension. Some of these terms are induced, at the microscopic level, by the SM mechanisms discussed above. The effective chiral Lagrangian induced by the θ-term is discussed in Refs. [65,66]. BSM theories such as supersymmetry, multi-Higgs scenarios, left-right symmetric model etc. would give rise to additional PVTV sources of dimension six (and higher) in the quark-gluon Lagrangian [67]. The χEFT Lagrangians originating from these sources were derived in Refs. [68,69]. Various terms in the resulting effective chiral Lagrangian possess different scaling with respect to the underlying microscopic PVTV sources. χEFT can thus be used to establish relations between the fundamental PVTV mechanisms and specific terms in the nuclear potentials and, accordingly, specific pattern in the corresponding nuclear observables [65,68,69]. In principle, this offers the possibility to identify the fundamental sources of time-reversal violation and to shed light on some of the BSM scenarios, provided the corresponding LECs in the effective Lagrangian can be determined from Lattice QCD calculations or experimental data [70,71].
In the framework of χEFT, the PVTV NN potential was derived up to N 2 LO including one-and two-pion exchange contributions and the corresponding contact interactions [72,73]. Subsequent works showed the presence in the PVTV Lagrangian of a three-pion term [68], which was for the first time included in the calculations in Ref. [66]. This term also generates a PVTV 3N force at NLO, which contributes to the 3 H and 3 He EDM. The calculation reported in Ref. [66] was also the first one carried out using solely the interactions derived in χEFT. More precisely, the PVTV potential at NLO was used in combination with the N 2 LO PCTC potentials from Ref. [74]. Finally, in Ref. [75], the EDM of deuteron and trinucleons was studied using the χEFT PVTV potential up to N 2 LO along with the N 4 LO PCTC potential of Ref. [18]. In this paper, it was also shown that the N 2 LO contribution to the PVTV 3N force generated by the three-pion interaction vanishes. The LO χEFT PVTV potential has also been applied in combination with many-body methods to calculate Schiff moments of heavy nuclei [76].
Currently, no direct limits on EDMs of light nuclei have been established. However, experiments are planned to measure the EDM of protons and light nuclei in dedicated storage rings [77,78,79,80,81,82]. This new approach could reach an accuracy of ∼ 10 −16 e fm, although this has to be established in practice. If successful, these experiments would lead to a great improvement in the hadronic sector of EDM searches. A measurement of a non-vanishing EDM of this magnitude would provide evidence of a PVTV source beyond the CKM mechanism. However, a single measurement would be insufficient to identify the source of PVTV. For this reason, experiments with various light nuclei such as 2 H, 3 H and 3 He are planned. Such measurements would provide a complementary information needed to impose constraints on PVTV sources at the fundamental level.
A brief discussion of the PVTV potentials derived in the framework of the one-meson exchange model and in the pionless EFT approach will also be reported in this review.

Outline of the article
Our paper is organized as follow. In Section 2, we discuss the origins of PVTC and PVTV interactions at the fundamental level and list the relevant terms in the quark-gluon Lagrangian. In Section 3, we give the corresponding terms in the effective chiral Lagrangian and discuss the derivation of the PVTC and PVTV potentials in χEFT. In Section 4, we specifically focus on the contact few-nucleon interactions which enter the potentials in both chiral and pionless EFT formulations. We also discuss the expected hierarchy of the corresponding LECs as suggested by the large-N c analysis. Next, in Section 5, the various mesonexchange models developed to describe the PVTC and PVTV interactions will be summarized. Then, in Section 6, we report on a selected set of results for PVTC and PVTV observables in light nuclei up to A = 3. Finally, the main conclusions of this paper and future perspectives are summarized in section 7.

PARITY VIOLATION AND TIME-REVERSAL VIOLATION AT THE MICROSCOPIC LEVEL
Parity is violated in the SM of particle physics because of the different gauge interactions of left-and righthanded fermion fields. Only left-handed particles interact via SU(2) L gauge interactions such that this part of the SM violates parity maximally. The remaining color and electromagnetic interactions conserve parity modulo the QCD vacuum angle which is discussed below. Parity violation was first observed in semileptonic charged current interactions in 1957 [83]. Twenty years later, in the late '70s, PVTC was observed in neutral current electron-nucleus scattering [84], providing a strong confirmation of the SM. Subsequent PVTC electron scattering experiments have quantitatively confirmed the SM picture [85]. In addition to PVTC in β decays and semileptonic neutral current processes, the SM predicts PVTC in weak interactions between quarks. At energies smaller than the masses of the W and Z bosons, these interactions can be represented by four-fermion operators. Right below the electroweak (EW) scale, and limiting ourselves to the lightest u and d quarks, the four-fermion Lagrangian is where G F is the Fermi coupling constant and s 2 w ≡ sin 2 θ W = 0.231, with θ W the Weinberg mixing angle. q L and q R denote the left-handed and right-handed doublets q T L = (u L , d L ) and q T R = (u R , d R ), and the dots denote terms that conserve parity 5 . Eq. (1) was obtained assuming the CKM matrix to be the identity, that is V ud = 1. The three operators in Eq. (1) all break parity, but have different transformation properties under chiral symmetry and isospin. We note that the isovector and isotensor terms (the second and third operators) given in Eq. (1) are suppressed by a factor s 2 w with respect to the isoscalar one.
The operators in Eq. (1) need to be evolved using the renormalization group equations (RGE) from the EW scale down to the QCD scale, and in this process they mix with additional PVTC operators [86]. After 5 Here u and d denote the uand d-quark Dirac fields, respectively. Moreover u R,L = 1±γ 5 2 u, etc.
the RGE evolution, the PVTC Lagrangian assumes the form where in the SM, the coefficients C SM i are known functions of SM parameters as s w , the strong coupling constant g s , etc. Greek indices α and β appearing as superscript in some of the quark fields in Eq. (2) specify color indices. They are only shown for cases where the color contractions are not obvious. Notice that the QCD evolution does not remedy the s 2 w suppression of the isospin-one and -two operators [86]. BSM physics that arises at scales well above the EW can be represented at the EW scale via gaugeinvariant higher-dimensional operators [87,67]. This framework is usually called the SM Effective Field Theory (SM-EFT). SM-EFT operators can induce new PVTC couplings of the W and Z bosons to left-and right-handed quarks, and new PVTC four-fermion operators. After evolving the effective operators from the EW to the QCD scale, the net effect of BSM PVTC SM-EFT operators is to modify the coefficients C SM i in Eq. (2) with respect to their SM values, namely in Eq.
We have focused so far on operators involving only the u and d quarks. Flavor-conserving (∆F = 0) operators involving the s quark can also generate interesting contributions to hadronic P violation [88,86], such as contributions to isospin-one operators that are not suppressed by s 2 w .
While P and C are maximally broken by the V −A structure of the SM, the breaking of CP is much more delicate. In the SM with three generations of quarks, CP is broken by the phase of the CKM matrix, which explains all the observed CP violation in the kaon [89,90,91], and B meson systems [92,93]. Theoretical uncertainties are at the moment too large to definitively conclude whether the recently discovered CP violation in D decays [94] is compatible with the SM. The phase of CKM gives, on the other hand, unobservable contributions to flavor-diagonal CP violation, in particular to the neutron [95,55,59] and electron EDMs [96,97,98].
The second source of CP violation in the SM is the QCD θ term [99,51,100] where g s is the strong coupling constants and G a µν the gluon field tensors (a is a color index). The θ term is a total derivative, but it contributes to physical processes through extended, spacetime-dependent field configurations known as instantons. CP violation from the QCD θ term is intimately related to the quark masses. All phases of the quark mass matrix can be eliminated through non-anomalous SU(2) vector and axial rotations, except for a common phase ρ. The mass plus QCD θ terms which are left are where M = diag(m u , m d ). The parameters ρ and θ are not independent. In χEFT, it is convenient to rotate L mass+θ P V T V into a complex mass term with an anomalous U(1) A rotation, obtaining, after vacuum alignment [101], n f = 2 is the number of light flavors, and the combinations of light quarks massesm and ǫ are 2m = Eqs. (5) and (6) can be easily generalized to include strangeness. θ is a free parameter in the QCD Lagrangian, and one would expectθ = O(1). This would however lead to a large neutron EDM |d n | ∼ 10 −3θ e fm [102,103], ten orders of magnitude larger than the current limits, d n < 3.0 · 10 −13 e fm [61]. Thereforeθ 10 −10 , the so-called strong CP problem.
The phase of the CKM and the QCDθ term are the only CP-violating parameters in the SM Lagrangian. They are however not sufficient to explain the observed matter-antimatter asymmetry of the Universe [104,105,106,107], and it is therefore natural to think about CP-violating sources induced by BSM physics. The low-energy CP-violating operators relevant for EDMs have been cataloged in several works, e.g. Refs. [52,108,109,110]. Ref. [68] considered all the low-energy operators that are induced by SM-EFT operators at tree level, retaining the two lightest quarks. Generalization to three flavors are given, for example, in Refs. [111,112]. The most relevant SU(3) c × U(1) em -invariant purely hadronic operators induced by dimension-six SM-EFT operators are L 6,hadr where f abc are the structure constants of the Lie algebra of the color 1,2 and Ξ (ud) 1,2 are dimensionless and scale as (v/Λ X ) 2 , where v = 246 GeV is the Higgs vacuum expectation value, and Λ X is the scale of new physics. The Weinberg three-gluon, the quark EDM (qEDM), and the chromo-EDM (qCEDM) operators (the first, second, and third term given in Eq. (7), respectively) have received the most attention in the literature [52,113]. They can be written directly in terms of SU(3) c × SU(2) L × U(1) Y -invariant operators at the EW scale, and receive corrections by a variety of CP-violating operators in the SM-EFT, involving heavy SM fields. The fourquark operators, given in the second line of Eq. (7), can also be expressed in terms of gauge-invariant operators at the EW scale, and they arise, for example, in leptoquark models, see [114,115]. The fourquark operators, given in the third line of Eq. (7), are on the other hand induced by right-handed couplings of quarks to the W boson [116,68], and are generated, for example, in left-right symmetric models. While all operators in Eqs. (5) and (7) violate P and CP symmetry, they transform differently under isospin and chiral rotations. As such, the operators induce different χEFT Lagrangians at lower energies, and different hierarchies of CP-violating hadronic and nuclear observables such as EDMs or scattering observables.

PVTC AND PVTV CHIRAL POTENTIALS
In this section, we discuss the derivation of the PVTC and PVTV NN and 3N potentials within the framework of χEFT. In the first and second subsections we briefly review the properties of the PVTC and PVTV chiral Lagrangians. In subsection 3.3, we present briefly two methods used to derive the potentials starting from a Lagrangian. Finally, in the last two subsections, we present the PVTC and PVTV chiral potentials, respectively.
In order to discuss hadronic observables such as nuclear EDMs or PVTC asymmetries in pp scattering, the quark-level PVTC and PVTV Lagrangians of Eqs. (2) and (7) need to be matched onto nuclear EFTs, such as chiral EFT and pionless EFT. Due to the nonperturbative nature of QCD at low energy, this matching cannot be done in perturbation theory. Nevertheless, the approximate chiral and isospin symmetries of the QCD Lagrangian provide an organizing principle for low-energy interactions, see [12,13,14] for review articles.
In some cases, we will perform a non-relativistic reduction of the nucleon field N(x) and use N r (x) where p r (x) (n r (x)) is the two component Pauli spinor representing the static proton (neutron) field. Effects of the anti-nucleon degrees of freedom are taken into account in the form of 1/M relativistic corrections to the vertices. The coefficient of the annihilation operator reduces to χ s exp(ip · x), where χ s is a spinor describing a spin state with z-projection s = ± 1 2 .
The main "building block" to construct the chiral Lagrangian is the SU(2) pionic matrix field U(x), often written as (but its definition is not unique) [12] U(x) = e i fπ π(x)· τ , where f π ≈ 92.4 MeV is the pion decay constant. Another low energy constant frequently entering the chiral Lagrangian is the axial coupling constant g A ≈ 1.29. Following the standard convention, we give here the effective value that takes into account the Goldberger-Treiman discrepancy and is extracted from the empirical value of the pion-nucleon coupling constant. The effective chiral Lagrangian is constructed in terms of N(x) and U(x) and therefore contains vertices with arbitrary number of pion fields. In the following, we will retain explicitly only relevant terms with the minimum number of pion fields, obtained by expanding U(x) in powers of the pion field. Additional terms with a larger number of pion fields will only contribute to the PVTC and PVTV potential at higher orders in the chiral expansion. For an introduction to the chiral Lagrangians and their building blocks, the reader is referred to Refs. [12,13,14] and references therein.
Each term of the chiral Lagrangian will be classified by the so-called "chiral order". Each four-gradient of the pion matrix field or a multiplication by a pion mass increases the order of the term by one. Fourgradients acting on nucleon fields are more difficult to classify, since the time derivative brings down a factor proportional to the nucleon mass. An easier counting is obtained using the non-relativistic heavy baryon perturbation theory [117,12], which was used in the derivation of the PVTC potential in [46] and of the PVTV potential in [72]. In the following, we will use both the relativistic and nonrelativistic nucleon fields.
For the sake of completeness, we report first of all the terms of the PCTC Lagrangian that contribute to the PVTC and PVTV potentials up to the order we are interested in. In SU(2) χPT, the PCTC Lagrangian can be conveniently organized in sectors with different numbers of pions and nucleons (below we give the explict expression for the relevant terms in the πN Lagrangian only).
L P CT C = L P CT C,πN + L P CT C,N N + L P CT C,ππ + · · · , where "· · · " in the previous expression denotes terms of higher order and/or more pions fields of no interest here. Above The parameters c i=1−4 are LECs appearing in the Lagrangian of order Q 2 . They have dimension of the inverse of an energy. For a complete discussion of the terms appearing in the Lagrangians L P CT C,πN , L P CT C,N N , and L P CT C,ππ , etc., see Refs. [12,118].

The PVTC chiral Lagrangian
The effective chiral Lagrangian that involves contributions from the weak sector of the SM was first discussed in the seminal paper by Kaplan and Savage [88] and subsequently revisited in Refs. [119,40,46,42]. Also the PVTC Lagrangian can be conveniently organized in sectors with different numbers of pions and nucleons, explictly where the superscript (n) denotes the chiral order of each piece. The pion-nucleon interaction terms are collected in L P V T C,πN and those entering the PVTC potential up to Q 1 are the following [88,46] L (0) where The parameters h 1 π and h ∆I V,A are unknown LECs. The superscript ∆I labels the rank of the corresponding isospin tensor. The LECs can be estimated by naive dimensional analysis (NDA) [88,40,46,42] where Λ χ = 4πf π ∼ 1.2 GeV is the typical scale of the strong interaction. Eq. (21) shows the order-ofmagnitude estimates of the PVTC interactions. These estimates do not take into account factors of s 2 w and N c that could modify the expected scaling of the LECs.
The contact terms entering the Lagrangian L P V T C,N N are products of a pair of bilinears of nucleon fields that are odd under P and even under CP. The most general bilinear product reads where Γ A and Γ B are elements of the Clifford algebra with the possible addition of 4-gradients and F ab are unknown parameters. To violate P but conserve CP, at least one 4-gradient is required. We must build isoscalar, isovector and isotensor terms as discussed in Section 2. The operators moreover have to conserve the electric charge and thus commute with the third component of the isospin operator. The terms with only one gradient operator are collected in L (1) P V T C,N N (i.e. of chiral order 1). Only five independent terms can be written [44], corresponding to the five possible S ↔ P transitions in NN scattering [45]. It is more convenient to give the Lagrangian using the non-relativistic reduction of the nucleon fields N r : The factor 1 Λ 2 χ fπ has been chosen to ensure C i are dimensionless and for convenience in the power counting. The construction of Eq. (23) and the elimination of redundancies will be discussed in more details in Section 4. The operators multiplying the LECs C 1,2 are isoscalar, those multiplying C 3,4 change isospin by one unit, while that multiplying C 5 is an isotensor. The scaling of the LECs from naive dimensional analysis [120] is given by which once again does not take into account the suppression by s 2 w affecting, for example, the isovector operators. The operators in Eq. (23) contribute to the PVTC potential at NLO (suppressed by (Q/Λ χ ) 2 with respect to LO), and we will give the potential derived from them in Eq. (60). The terms appearing in L P V T C,N N contain two additional gradients and contribute to the PVTC potential at higher order. They have not been considered so far.
Finally, there are some terms with 3π vertices appearing in L (2) P V T C,πππ as discussed in Ref. [42]. These terms would contribute to the Q 2 PVTC potential, but their contributions at the end vanishes as discussed in Subsect. 3.4.

Connection to the underlying PVTC sources
Attempts to estimate the values of the coupling constants were performed mainly in the framework of the meson exchange models (which will be discussed in Section 5). However, since in both χEFT and meson exchange frameworks the lowest order pion-nucleon Lagrangian term is the same as given in Eq. (18), we can report here the values for h 1 π estimated from the underlying fundamental theory also before the advent of χEFT [121,122,123,124,125,126,127]. One of the most comprehensive calculation including all previous results was performed in 1980 by Desplanques, Donoghue, and Holstein (DDH) [50] using the valence quark model. Additional calculations have been performed later [128,129,130], using similar or other methods and finding qualitatively similar results. These estimates, however, are based on a series of rather uncertain assumptions (see, for example, Ref. [131]). For example, DDH presented not a single value for h 1 π but rather a range inside of which it was extremely likely that this parameter would be found [50]. In addition they presented also a single number called the "best value" but this is described simply as an educated guess in view of all the uncertainties. The values of h 1 π were [50] DDH: h 1 π = 4.56 × 10 −7 ("best value") , h 1 π = 0 − 11.4 × 10 −7 ("reasonable range") . (25) Some years ago, a lattice QCD calculation of h 1 π has also been made [132], resulting in the following estimate Lattice: where the theoretical uncertainty is related to the statistical Monte Carlo error. While the systematic errors are expected to be within the quoted statistical error [132], we stress that the calculation was performed at a heavy pion mass and not extrapolated to the physical point, disconnected diagrams were not included, and operator renormalization was neglected.
Regarding the other LECs entering the contact Lagrangian given in Eq. (23), no direct estimates were reported in literature. These LECs were estimated by comparing the expression of contact potential with the potential developed using the exchanges of heavy mesons, as for example, in the DDH potential [46,42] (this issue will be considered in more detail in Sect. 5). However, since also the DDH estimates are rather uncertain, here we will not discuss this issue.

The PVTV Lagrangian
The PVTV chiral Lagrangian taking into account the QCDθ term was first considered in the seminal paper by Crewther, di Vecchia, Veneziano and Witten [102], and consequently revisited in Refs. [133,134,135,136,137]. Subleading terms in the chiral expansion were systematically constructed in Refs. [65,69]. The chiral Lagrangian induced by the dimension-six operators in Eq. (7) were derived in Refs. [68,69].
As before, in SU(2) χPT, the PVTV Lagrangian can be organized in sectors with different numbers of pions and nucleons As in the previous subsection, we report here only the most important interactions for each sector, focusing on the terms with the minimum number of pion fields entering in the final expression of the potential. Terms with additional pions are not universal for the different PVTV sources at the quark level, but instead depend on their chiral-symmetry breaking pattern. These differences only enter at higher order in the potentials than we consider here.
In the PVTV case, the simultaneous violation of P, T, and isospin symmetry allows for a pion tadpole linear in the pion field ∼ π 3 with a corresponding LEC proportional to the symmetry-violating source terms at the quark level. Such tadpoles can always be removed by appropriate field redefinitions of the pion and nucleon fields [65,68,69]. At LO in the chiral expansion, the tadpole removal is the same as the vacuum alignment procedure at the quark level [101]. While tadpoles can be removed, the corresponding field redefinitions affect other couplings in the chiral Lagrangian. In particular, for chiral-symmetrybreaking CP sources that do not transform as a quark mass term, a PVTV three-pion vertex of chiral order Q 0 is left behind [68,69].
where∆ is a LEC. Other three-pion vertices will appear at N 2 LO, but they will contribute to high orders of the PVTV potential.
Arguably the most important interactions appear in the pion-nucleon sector. Simultaneous violation of P, T, and chiral symmetry allows for non-derivative single-pion-nucleon interactions, something which is not possible in the PCTC Lagrangian. In principle, three different interactions can be written corresponding, respectively, to an isospin singlet, vector, and tensor interaction. As discussed below, the relative size of the LECsḡ 0,1,2 strongly depends on the quark-level PVTV source under consideration. In the case of CP-violation from chiral invariant operators, such as the gCEDM,ḡ i are suppressed by powers of the pion masses, and the pion-nucleon Lagrangian contains chiral-invariant, derivative couplings as important as those in Eq. (32) [68]. These can however always be absorbed into a shift ofḡ 0 and of the ∆I = 0 NN operators discussed below.
The NLO Lagrangian contains several two-pion two-nucleon PVTV interactions [65,68,69,75], but, for all CP-violating sources, they contribute to the two-and three-body PVTV potentials at N 3 LO and N 2 LO, respectively. We therefore ignore these couplings. Isospin-breaking sources also generate a singlepion-nucleon NLO coupling. The coupling involves a time derivative of the pion field, thus inducing a relativistic correction in the O(Q) PVTV potential. At N 2 LO the number of interactions proliferates significantly and there are also new pure pionic interactions. These contributions can either be absorbed into LO LECs or appear at high orders in the PVTV potential considered here.
Apart from pionic and pion-nucleon interactions, there appear PVTV NN contact interactions. As in the PVTC case, at least one gradient is required such that these operators start at order Q. Terms with three or more gradients have not been considered so far. At order Q, only five independent interactions of this kind can be written, corresponding to the five possible S ↔ P transitions (see Section 4 for a general discussion of this kind of interaction terms). Neglecting terms with multiple pions, the Lagrangian reads (again, it is convenient to write it in terms of the non-relativistic nucleon field N r ) As suggested by the factor of Λ 2 χ which we pulled out of the definition of the LECs, in χEFT these operators contribute in general at N 2 LO and are suppressed with respect to the PVTV one-pion exchange (OPE) potential. The only exception, as discussed in Sec. 3.2.1, are quark-level operators that do not break chiral symmetry, for whichC 1,2 are as important as the contributions fromḡ 0,1 .
Finally, the calculation of EDMs or other PVTV electromagnetic moments requires the inclusion of electromagnetic currents. Nucleon EDMs are induced by pion loops involving the interactions in L (0,1) The renormalization of these loops requires the inclusion of short-distance counter terms contributing to Table 1. Scaling of the LECs in the chiral Lagrangian in dependence of the microscopic CP violation sources. We introduced the counting parameters With ǫ mπ ∼ ǫ χ , we introduced two different parameters to explicitly track insertions of the light quark masses from the QCD Lagrangian. ε is the isospin breaking parameter ε = (m d − m u )/(m d + m u ) ≃ 1/3. The scaling of the LECs induced by dimension-six sources assume a Peccei-Quinn mechanism. A "−" implies the interaction is only induced at higher order than considered here. The parametersC 1,2 ,C 3,4 , andC 5 are the LECs entering the contact PVTV potential, respectively of isoscalar, isovector, and isotensor type. the nucleon EDMs. Such counter terms indeed appear in the chiral Lagrangian where F αβ is the electromagnetic field strength andd 0 andd 1 are LECs related to the proton and neutron EDMs, respectively. The above interactions are sufficient for calculations of hadronic and nuclear PVTV scattering observables and EDMs up to NLO in the chiral expansion. Calculations of higher PVTV moments, such as magnetic quadrupole moments, can depend on additional LECs [138].

Connection to the underlying PVTV sources
In the previous section we listed the PVTV hadronic interactions relevant for observables of experimental interest. However, for a given PVTV source at the quark-gluon level, a specific hierarchy among the various interactions appear. The relative importance of the LECs in Eqs. (31), (32), (33) and (34) for the different microscopic sources of CP violation is summarized in Table 1. These estimates are based on NDA [120]. NDA is valid in the regime in which the strong coupling g s is non-perturbative, and, as done for NDA estimates of the chiral-invariant PCTC interactions, we will take g s ≃ 4π. In addition, for dimension-six sources, we assumed that a Peccei-Quinn mechanism [139] relaxesθ to an inducedθ ind , which depends on the coefficients and vacuum matrix elements of the operators in Eq. (7) [52,140,141]. The scaling of the couplings without this assumption can be found in Ref. [68]. To make the power counting explicit, we introduced three ratios of scales Numerically, ǫ χ ∼ ǫ mπ , but we define two different parameters to track the dependence of the LECs on the quark masses. To assess the size of the contribution of different CP violating sources to the nucleon and nuclear EDMs, the scaling of the LECs in Table 1 can be combined with a naive estimate of these observables. As we will discuss in detail in Sections 3.2.2 and 6.4, the nucleon EDM receives tree level contributions fromd 0,1 and loop contributions byḡ 0 andḡ 1 , leading to where e is the electric charge and the coefficients of the loops α 0,1 will be given explicitly in Section 3.2.2. The additional suppression ofḡ 1 is due to the fact that this coupling only involves neutral pions, which do not interact with a single photon at LO. Nuclear EDMs, on the other hand, receive tree level contributions from the single nucleon EDM, and from pion-nucleon and nucleon-nucleon couplings, The coefficients a n,p , a ∆,0,1,2 and A 1,...,5 depend on the nucleus under consideration, and in Section 6.4 we will present results for their calculation in chiral EFT for the deuteron, 3 H and 3 He. By power counting, they are expected to be O(1) (measured in units of fm in the case of the dimensionful a ∆,0,1,2 and A 1,...,5 ), barring isospin selection rules, which for example suppress the contributions of the isoscalar operatorsḡ 0 andC 1,2 in nuclei with N = Z, such as the deuteron [142,143] 6 .
The reader should be aware that the dimensionless Wilson coefficients of the dimension-six operators, c 1,2 also come with intrinsic suppression factors. These arise from the typical loop and chiral factors that appear in BSM models. For example, quark and gluon dipole operators are typically induced at the one-loop level, and the quark EDM and chromo-EDM coefficients come with explicit factors of the quark mass (already included in Eq. (7)). This implies that one can expect {c Of course this is just an estimate and certainly models exist where these operators appear only at the two-or higher-loop level. On the other hand, the four-quark operators Ξ and Σ can be induced at tree level, so that {Ξ, Σ} = O(ǫ Λ ). Once the matching coefficients are calculated in a given model, Table 1 and Eqs. (36)-(37) allow to identify the dominant low-energy operator and to get a rough idea of the EDM constraints. Table 1 highlights that the chiral and isospin properties of the quark-level CP-violating sources induce very specific hierarchies between different low-energy couplings. These hierarchies in turn imply different relations between the EDMs of the nucleon, deuteron, and three-nucleon systems, which, if observed, would allow to disentangle the various CP-violating sources. From Table 1 we see that chiral-symmetrybreaking sources, such asθ,c (u,d) g and Ξ (u,d) 1,2 , induce relatively large PVTV pion-nucleon couplings. These couplings appear in the table with entry 1, indicating no further suppression. In particular, the isoscalar θ term and isovector Ξ (u,d) predominantly induce, respectively,ḡ 0 andḡ 1 , while a qCEDM would yield both couplings with similar strengths. The consequence is that for these sources light nuclear EDMs are enhanced with respect to the nucleon EDM. For these chiral-symmetry-breaking sources, the contact nucleon interactions proportional toC i are suppressed in the chiral expansion because these operators involve an explicit derivative. The suppression can be explicitly seen combining the scaling in Table 1 with the explicit factor of ǫ χ in Eqs. (33) and (37).
Chiral invariant sources such as the Weinberg operatorC G and the four-quark operators Σ (u,d) 1,2 , on the other hand, require additional chiral-symmetry breaking to generateḡ 0,1 , as indicated by extra powers 6ḡ 0 andC 1,2 contribute to the deuteron EDM in conjunction with isospin breaking in the strong interaction, or via the spin-orbit coupling of the photon to the nucleons [143]. Both contributions are beyond the accuracy we work at in this paper. of ǫ mπ . In this case, EDMs of light-nuclei are expected to be of similar size as the nucleon EDM. Furthermore, the contact nucleon operators proportional toC 1,2 now contribute to the PVTV potential at the same order asḡ 0,1 . Finally, the qEDM mostly inducesd 0,1 , all other couplings being suppressed by O(α em ), where α em is the fine structure constant ∼ 1/137. In this case one expects nuclear EDMs to be dominated by the constituent nucleon EDMs.
While most statements are source-dependent, there is an important general message hidden in Table 1. There is no PVTV source for which the couplingsḡ 2 andC 3,4,5 appear at LO. For all sources they appear with a relative suppression of εǫ mπ or ǫ χ compared to other PVTV interactions. For most calculations one can simply neglect the associated interactions, reducing the number of LECs entering the expression of hadronic and nuclear observables. The suppression of the LECsḡ 2 andC 3,4,5 ultimately is a consequence of imposing gauge invariance on the dimension-six PVTV sources. Table 1 relies on NDA estimates for hadronic matrix elements [120]. A more quantitative assessment of the discriminating power of EDM experiments necessitates to replace the NDA estimates in Table 1 with solid nonperturbative calculations of the LECs. At the moment, there exist controlled estimates only of a few LECs. The pion-nucleon couplingsḡ 0 induced by the QCDθ term is related by chiral symmetry to modifications in the baryon spectrum [102]. In particular, in SU(2) χPTḡ 0 is related to the quark mass contribution to the nucleon mass splitting [65,144], up to N 2 LO corrections. Using Lattice QCD evaluations of the nucleon mass splitting [145,146], one finds where the 15% error includes both the Lattice QCD error on m n − m p , and an estimate of the error from N 2 LO chiral corrections. Unfortunately, chiral-symmetry-based relations do not allow to extractḡ 1 andd 0,1 .ḡ 1 has been estimated with resonance saturation leading toḡ 1 (θ)/ḡ 0 (θ) ≃ −0.2, somewhat larger than expected from NDA [73]. The LECsd 0,1 are usually estimated by naturalness arguments and considered to be of similar size to non-analytic contributions to the isoscalar and isovector nucleon EDM, see the Subsect. 3.2.2.
The relation between PVTV pion-nucleon couplings and corrections to the nucleon and pion masses is not specific to the QCDθ term, but can be generalized to all chiral-symmetry-breaking sources, such as for example the qCEDM [68,147] and Ξ (ud) 1,2 [148,141]. Since corrections to spectroscopic quantities should be easier to compute on the lattice, these chiral relations allow a calculation ofḡ 0,1 in Lattice QCD. While promising, this strategy has yet to lead to controlled results. The best estimate ofḡ 0,1 induced by the qCEDM comes from QCD sum rules [149,52] These estimates agree with NDA, especially forḡ 1 . However,ḡ 0 seems to be slightly suppressed, in agreement with large-N c expectations [150].
Only for the four quark operators proportional to Ξ (ud) 1,2 of Eq. (7) does the three-pion vertex with LEC ∆ appear at LO in the chiral Lagrangian. For this case, the LEC∆ is related by SU(3) symmetry to K → ππ matrix elements and K −K matrix elements that have been calculated on the lattice. We obtain with The matrix elements in Eq. (41) are in good agreement with NDA. The value of∆ also determines the tadpole component ofḡ 1 , which again is in line with NDA.
Most of the remaining LECs are undetermined at present. The focus of the Lattice QCD community has been on the matrix elements connecting the nucleon EDMs to theθ term [151,152,103], the qEDMs [153,154], the qCEDMs [152,155], and the Weinberg operator [156]. Some results are given in next subsection.

The nucleon EDM in chiral perturbation theory
The PVTV LECs defined in the previous section can be used to calculate the nucleon PVTV electric dipole form factor (EDFF). At zero momentum transfer, the EDFFs are identified with the nucleon EDMs. In dimensional regularization with modified minimal subtraction up to NLO in the chiral expansion, the EDMs are given by [137,157] where µ is the dimensional regularization scale. The leading loops proportional toḡ 0 are divergent and renormalized by the µ-dependent LECsd 0,1 . The NLO corrections proportional to m π /M are finite. The LEC∆ does not contribute at this order for any of the PVTV sources. As standard in χPT, the loops are associated to inverse powers of (4πf π ) 2 = Λ 2 χ . Combined with the scaling of the LECs in Table 1, we conclude that for theθ term and the qCEDMs the leading loop proportional toḡ 0 and the counter terms d 0,1 appear at the same order. For all other PVTV sources, the short-range counter termsd 0,1 are expected to dominate the nucleon EDMs. In no scenario can the EDMs be calculated solely from the pion-nucleon LECsḡ 0,1 as is often assumed in the literature. Estimates for the nucleon EDMs are often obtained by setting µ = M andd 0,1 (µ = M) = 0 such that EDMs depend on the value ofḡ 0,1 , which for some PVTV sources is better known.
The separation between the short-range and loop contributions is scheme dependent and therefore not physical. Lattice QCD calculations can therefore only calculate the total nucleon EDMs d n and d p . In recent years, significant efforts have been made towards calculating the nucleon EDMs in terms of the underlying PVTV sources. Most efforts have focused on the QCDθ term and the qEDM. The most recent results for theθ term [103] give in good agreement, but with sizeable uncertainties, with expectations from the chiral logarithm in Eq. (42) using Eq. (38). In the case of the qEDM, the nucleon EDM is related to the tensor charges, which have been computed with good accuracy [153,154]. Using the FLAG average [154], we get where Q u,d are the u and d-quark charges in units of the electric charge, and g u,d T the u and d-quark tensor charges of the proton, and the error on the r.h.s. of Eq. (45) is dominated by the uncertainty on the light quark masses.
On a longer time-scale, calculations of the qCEDMs and the Weinberg operator are also targeted. For now, the best results come from calculations using QCD sum rules [52,158].

From the Lagrangian to the potential
In this subsection, we briefly present two methods that have been used to derive nucleon-nucleon potentials starting from a Lagrangian. We first introduce the notation used here and in the next subsections.
The process under consideration is the scattering of two nucleons from an initial state |p 1 p 2 to the final state |p ′ 1 p ′ 2 (hereafter the dependence on the spin-isospin quantum numbers is understood). It is convenient to define the momenta where p j and p ′ j are the initial and the final momenta of the nucleon j. Furthermore it is useful to define which are the spin (isospin) matrix element between the final state s ′ j (t ′ j ) and the initial state s j (t j ) of the nucleon j.
Because k 1 = −k 2 ≡ k from the overall momentum conservation p 1 + p 2 = p ′ 1 + p ′ 2 , the momentumspace potential V is a function of the momentum variables k, K 1 and K 2 , namely Moreover, we can write in general where K = (K 1 − K 2 )/2, P = p 1 + p 2 = K 1 + K 2 , and the term V (P ) (k, K) represents a boost correction to V (CM) (k, K), the potential in the center-of-mass frame (CM). Below we will ignore the boost correction and provide expressions for V (CM) (k, K) only. Note that in the CM we define also p 1 = −p 2 ≡ p and p ′ 1 = −p ′ 2 ≡ p ′ . So we have k = p ′ − p and K = (p ′ + p)/2, so in the following we also write V (CM) as V (CM) (p, p ′ ). From now on, we will suppress the superscript "(CM)" for simplicity.
In order to derive the potential, two methods have been frequently used, the method of unitarity transformation (UT), and the method of the time-ordered perturbation theory (TOPT). They are briefly introduced below.
The time-ordered perturbation theory method. Let us consider the matrix element of the T -matrix, the "amplitude" of a process of scattering of two nucleons. Its square modulus |T f i | 2 is directly related to the cross section of the process. The conventional perturbative expansion for this matrix element is given as where E i is the energy of the initial state, H 0 is the Hamiltonian describing free pions and nucleons, and H I is the Hamiltonian describing interactions among these particles. These operators are defined to be in the Schrödinger picture and they can be derived from the Lagrangian constructed in terms of pions and nucleons as described, for example, in Refs. [159,118]. The evaluation of T f i is carried out in practice by inserting complete sets of H 0 eigenstates between successive H I factors. Power counting is then used to organize the expansion in powers of Q/Λ χ ≪ 1, where Q stands for either an external momenta or the pion mass. We will use the "naive" Weinberg counting rules [2], namely, we will count simply the powers of both the external momenta and pion mass insertions (we will consider low energy processes only). Each term will be of some order (Q/Λ χ ) ν . The terms with the lowest power of ν will be the LO, and so on.
In the perturbative series given in Eq. (50), a generic contribution will be characterized by a certain number of vertices coming from the interaction Hamiltonian H I and energy denominators, and it can be visualized also as a diagram (hereafter referred to as a TOPT diagram). Each vertex will give a "vertex function" and a δ conservation of the momenta of the particles involved in the vertex. The vertex functions are the results of the matrix elements of terms appearing in H I and are given as products of Dirac fourspinors, momenta, etc. A sum over the momenta of the particles entering the intermediate states is also present. When a diagram includes one or more loops, the δ's are not sufficient to eliminate all the sums over the momenta of the intermediate states. The energy denominators come from the factors 1/(E i − E α + iǫ), where E α is the (kinetic) energy of a specific intermediate state entering the calculation. The chiral order of each diagram can be calculated as follows. One needs to consider: 1. The chiral order of the vertex functions, which can be calculated from the non-relativistic (NR) expansion of the nucleon Dirac four-spinors (1/M expansion), and from various other factors. Typically, the powers of p/M coming from the NR expansion of the nucleon Dirac four-spinors are counted as ∼ Q 2 [2,4,160]. In other approaches however they are considered to be of order Q [20,161,162]. In this paper, we will follow the first prescription. 2. The energy denominators. We note that typical momenta p of the nucleons are much smaller than the mass of the nucleons, so we can treat them non relativistically. Namely . Usually in the energy denominator all the nucleon masses M cancel out and therefore we have two cases: • if there are no pions in the intermediate state, the energy denominator has only nucleon energy terms so it results of order 1/Q 2 .
• if there are pions in the intermediate states, the energy denominator reads where the term ∆E = E 1 + E 2 + · · · − E i where E 1 , . . . are the energies of the nucleons in the intermediate state and E i is the initial scattering energy. In the Taylor expansion the first term is of order Q −1 , while the other terms are usually called "recoil corrections". For the sake of consistency with the choice discussed above regarding the NR expansion of the Dirac 4-spinors, here we will count the p/M terms coming from recoil corrections as Q 2 as well. 3. The number of loops, or better the number of the sums over the intermediate state momenta that remain after using the conservation δ's. Each loop at the end will give a contribution of order Q 3 . 4. The number of disconnected parts of the diagram. For each of these parts, a δ factor expressing the momentum conservation of each part is present. Then, if there are N D disconnected parts, one of the δ simply gives the total momentum conservation, a factor common to all diagrams and therefore not relevant. Each of the remaining N D − 1 δ's at the end will "block" a sum over an external three-momentum, each one therefore reducing the chiral order by 3 units.
Once the T -matrix has been calculated, one would obtain in general where T (n) In all cases the sum starts from a minimum value n min , n min = 0 for the PCTC and n min = −1 for the PVTC and PVTV amplitudes. The idea now is to "define" the potential acting between the two nucleons so that it can reproduce the same amplitude T f i , namely, so that (for more details, see Ref. [118]) where H is the non-interacting Hamiltonian of two nucleons. Clearly, this procedure is not unique, since usually one imposes the relation T V = T f i to hold "on shell", namely by requiring the conservation of the energy between initial and final states. This induces an ambiguity, as discussed for example in Ref. [162]. However, the obtained potentials are expected to be equivalent by means of a unitary or at least a similarity transformation [163].
Finally, to invert Eq. (53), one assumes that V has the same Q expansion as the T matrix, and Eq. (53) can be solved for V (n) order-by-order (see, for example, Ref. [118] for more details). This procedure can be generalized to the A = 3 case to define a three-nucleon potential and so on.
The method of unitarity transformation. The method of unitary transformation (MUT) has been pioneered in the fifties of the last century to derive nuclear potentials in the framework of pion field theory [164,165]. In the context of chiral EFT, this approach was formulated in Refs. [166,167]. Similarly to TOPT, the MUT is applied to the pion-nucleon Hamiltonian which can be obtained from the effective Lagrangian in a straightforward way using the standard canonical formalism. Let η and λ denote the projection operators on the purely nucleonic subspace and the rest of the Fock space involving pion states with the usual properties η 2 = η, λ 2 = λ, ηλ = λη = 0 and η + λ = 1. To derive nuclear forces and/or current operators, the Hamiltonian needs to be brought into a block-diagonal form with no coupling between the ηand λ-subspaces, which can be achieved via a suitably chosen unitary transformation U. Following Okubo, a unitary operator can be conveniently parametrized in terms of the operator A = λAη that mixes the two subspaces via One then obtains the nonlinear decoupling equation for the operator A: The  50), and the corresponding matrix elements can also be interpreted in terms of TOPT-like diagrams. Notice that contrary to Eq. (50), the expressions in the MUT do, per construction, not involve energy denominators that vanish in the static limit of infinitely heavy nucleons and correspond to iterative contributions to the scattering amplitude. As explained in [167], in order to implement the chiral power counting in the algebraic approach outlined above it is convenient to rewrite it in terms of different variables. Using the rules given in the description of the TOPT approach and counting the powers of the soft scale Q for a given irreducible (i.e. of non-iterative type) connected N-nucleon TOPT-like diagram without external sources, one obtains for the chiral order n [2, 167] where L is the number of loops, V i is the number of vertices of type i. Further, the vertex dimension ∆ i is given by ∆ i = d i + 1/2n i − 2 with d i and n i being the number of derivatives and/or m π -insertions and the number of nucleon fields, respectively. The above expression is convenient to use for estimating the chiral dimension of TOPT-like diagrams. For the MUT, it is, however, advantageous to rewrite it in the equivalent form where p i is the number of pionic fields. The parameter κ i obviously corresponds to the inverse overall mass dimension of the coupling constant(s) accompanying a vertex of type i. In this form, chiral expansion becomes formally equivalent to the expansion in powers of the coupling constants, and it is straightforward to employ perturbation theory for solving the decoupling equation (56) and deriving the nuclear potentials ηHη. One nontrivial issue that emerges when applying chiral EFT to nuclear potentials concerns their renormalization. While on-shell scattering amplitudes, calculated in chiral EFT, can always be made finite by including the counterterms from the effective Lagrangian (provided one uses a chiral-symmetry preserving regularization scheme such as dimensional regularization), nuclear potentials represent scheme-dependent quantities, which correspond to non-iterative parts of the scattering amplitude. There is no a priori reason to expect all ultraviolet divergences emerging from TOPT-like diagrams, which give rise to nuclear forces, to be absorbable into a redefinition of the LECs. Indeed, it was found that the static PCTC three-nucleon force at order Q 4 of the two-pion-one-pion exchange type cannot be renormalized if one uses the unitary transformation given in Eq. (55) [168]. On the other hand, the employed parametrization of the operator U is clearly not the most general one and represents just one possible choice. The freedom to change the off-shell behavior of the nuclear potentials, already mentioned in the context of TOPT, has been exploited in a systematic way in the PCTC sector in order to enforce renormalizability of nuclear forces (using dimensional regularization) [167,169,170,171,172]. The MUT has also been successfully applied to the effective Lagrangian in the presence of external classical sources in order to derive the corresponding nuclear current operators, see [160] and references therein.

The PVTC potential up to order Q 2
In this subsection we will discuss in detail the derivation of the PVTC potential up to N 2 LO using the TOPT approach. We consider diagrams contributing to the T -matrix with one vertex coming from the PVTC Lagrangian, while all other vertices coming from the PCTC interaction. Diagrams with two or more PVTC vertices can be safely neglected.
The TOPT diagrams contributing to the PVTC T-matrix up to N 2 LO are shown in Fig. 1 in panels (a)-(l). The one pion exchange diagram (a) gives a contribution to the T -matrix of order Q −1 (that will be our LO). The diagram (b) represents a PVTC contact interaction of order Q; also the diagrams (c) and (d) with the PCTC contact vertex and one pion exchange are of order Q. The triangle diagram (e) with a PCTC ππNN vertex is of order Q, while if we consider the PVTC ππNN vertex as in panel (l) the diagram is of order Q 2 . The box diagrams (f) and (g) includes contribution of order Q 0 and Q; the contribution of order Q 0 is exactly canceled when inverting Eq. (53). Finally, the "bubble" diagram (h), the three-pion vertex (2) (4) (1) (3) Figure 2. Other diagrams that would contribute at NLO. These diagrams contribute to the renormalization of the LECs (panels (1), (2) and (3)) or give a vanishing contribution to the potential (panel (4)) due to the inversion of Eq. (53). Notation as in Fig. 1. diagram (i), the box diagram (j) with the πNN vertex coming from the subleading PVTC Lagrangian terms proportionals to the LECs h i V , and also the diagram (k) with the ππNN vertex coming from the subleading PCTC Lagrangian terms proportionals to the LECs c i , are of order Q 2 . These latter diagrams were considered for the first time in [173] using the MUT, and using TOPT in [174].
Contributions proportional to 1/M coming from the NR expansion of the vertex functions or from recoil corrections in this work are considered to be at least of order N 3 LO.
Other types of diagrams like those shown in panels (1), (2), (3) of Fig. 2 simply contribute to a renormalization of the coupling constants and masses, see Ref. [42] for more details. In the following, we will disregard these diagrams, but it should be taken into account that the formulas below are given in terms of the renormalized (physical) LECs and masses. The contribution of diagram (4) is cancelled when inverting Eq. (53).
Let us now consider each kind of diagram separately: • One pion exchange (OPE) diagram. Diagram (a) of Fig. 1 gives the LO contribution (Q −1 ) to the potential where ω k = k 2 + m 2 π , and it comes directly from the LO expansion of the vertices and energy denominators. Derived from the same diagram, there are terms coming from the NR expansion of the vertices, the first correction being of order (p/M) 2 . However, as discussed previously, they are counted to be of order Q 4 , and thus the corresponding terms are considered to be suppressed by four orders with respect to V P V T C,N N . They give a contribution to the potential of order Q 1 . As discussed in Chapter 4, this contribution can be written in various equivalent forms due to the Fierz identities [44]. We have chosen to write this part as follows [42] V (1) where Λ χ = 4πf π ≈ 1.2 GeV. The parameters C i , i = 1, . . . , 5 are LECs. Different (but equivalent) forms of this part were used in Refs. [41,173]. • Contact plus OPE diagrams. The diagrams (c) and (d) in Fig. 1 are representative of diagrams with a contact term and an OPE. However all these diagrams vanish after the integration over the loop variable. • NLO two pions exchange: triangle diagrams. There are 6 different time-orderings of diagrams given in panel (e) in Fig. 1. After summing them, the total contribution from the triangle diagrams results to be [40,175] V (1) where ω ± = (q ± k) 2 + 4m 2 π . The integral is singular and has to be regularized using some method. We will discuss this issue later. • NLO two pions exchange: box diagrams. There are 48 diagrams represented by the diagrams of type (f) and (g) of Fig. 1 when we consider all possible time orderings. The final contribution is [40,175] V (1) and it is of order Q 1 . Again the integral is singular. In this case, in the amplitude T f i there appears a term of order Q 0 coming from diagram (g), but it cancels out when inverting Eq. (53). • Bubble diagrams. We now turn to the diagrams contributing at order Q 2 , that is at N 2 LO. The sum of "bubble" diagrams depicted in panel (h) of Fig. 1 mutually cancel and these diagrams do not give any contribution to the PVTC potential. • Diagrams with three pion vertices. The expansion of the PVTC Lagrangian in terms of pions gives rise to two terms proportional to ( π) 3 which would contribute to T f i via the diagram depicted in panel (i) of Fig. 1 Figure 3. TOPT diagrams that would contribute to the PVTC 3N force For the notation see Fig. 1.
• N 2 LO two pion exchanges: triangle diagrams. The diagram depicted in panels (k) derives from a subleading ππNN vertices in the PCTC Lagrangians [173,75], see Eq. (13), Note in Eq. (64) the presence of the LEC c 4 , which belong to the PCTC sector [12]. The expression for the diagrams (l) comes from the LO PCTC and PVTC vertex functions. The final result is [173,75] V (2) where h 1 A and h 2 A are two of the LECs that appear in the Lagrangian terms given in Eq. (19).
Finally, we conclude this section by mentioning that at N 2 LO, one should also include PVTC 3N forces. Examples of diagrams contributing to this 3N force are reported in Fig. 3. The chiral order of diagrams with more than two nucleons is discussed in detail in Ref. [167]. The diagram depicted in panel (a) with a LO PCTC ππNN vertex would contribute at NLO but summing over all the time-orderings it vanishes. The other three diagrams (the one in panel (b) has a subleading PCTC ππNN vertex proportional to c i , i = 1, . . . , 4 [12]) are N 2 LO and therefore they have to be considered in order to perform fully consistent calculations in A ≥ 3 systems. These kind of diagrams have not yet been considered in literature. Note that diagrams with a 3N PVTC contact vertex are highly suppressed, so no new LEC has to be introduced in this case.

Regularization of the PVTC potential
In this section we deal with the divergences in the loop diagrams. We will briefly present three methods frequently used in literature, namely the dimensional regularization (DR) method used e.g. in [161], the spectral function regularization (SFR) [176], and the novel (semi-)local momentum-space regularization approach of Ref. [19].
• Dimensional regularization method. This technique is well known for dealing with divergences of loop integrals present in Feynman diagrams, where the integration is performed over four-momenta. In case of time-ordering diagrams, the loops involve integrations over 3-dimensional momenta. To deal with the singularities, the integrals are re-defined in d dimensions and successively one takes the limit d → 3. The singular part is singled out by terms ∼ 1/(3 − d), which then can be reabsorbed in some of the LECs. As usual, we define ǫ = 3 − d, and we assume that ǫ → 0. When we use the DR, it is better to "rescale" all the dimensional quantities with an energy scale µ. Therefore we define q =qµ, m =mµ, etc., where the "tilde" quantities are adimensional. We can now go to d dimension and manipulate the integrals as discussed in detail in Ref. [161], see also Ref. [177]. Here we limit ourselves to list the results needed to regularize the loop integrals we have encountered. Regarding the loop integrals appearing at NLO in Eqs. (61) and (62), we have where and which contains the divergent part (above γ is the Euler-Mascheroni constant). The loop integrals appearing in the N 2 LO diagrams as in Eqs. (64) and (65) are of the form The first integral is finite, but the second integrand diverges linearly for q → ∞. The finite contribution to potential can be obtained using the DR method. Alternatively, one can impose an ultraviolet cut-off Λ C on the integrals. The integrals then yield diverging parts as Λ C → ∞, which can be again reabsorbed by some LECs, finite parts independent on Λ C that are exactly the same as obtained using the DR method, and a number of other terms which can be expressed in power series of Q/Λ C , where Q is either k or m π . Taking the limit Λ C to infinity these latter parts would disappear. Since, in general we must fix Λ C at a value greater than the typical energies of the χEFT, then these additional terms carry at least an additional power of Q which means they give contributions at N 3 LO (or beyond) to the potential. Therefore, for the integral in Eq. (71), we have followed the prescription to absorb the divergent parts in some LEC's, disregard the parts depending on Q/Λ C , and retaining the finite parts as those given by the DR method. Explicitly, the two integrals are given by where • Spectral function regularization method. Pion loop integrals appearing in the two-pion exchange contributions discussed in the previous subsection can be generally expressed using a dispersive representation. Writing the momentum-space potentials in the general form V = i O i W i (k) with O i being spin-isospin-momentum operators and W i the corresponding structure functions that depend only on the momentum transfer k ≡ |k|, the unsubtracted dispersion relations for the functions W i (k) have the form [178] W where the spectral functions ρ i (µ) are given by ρ i = ℑ W i (0 + − iµ . Notice that the spectral integrals in Eq. (75) do not converge for potentials derived in chiral EFT since ρ i (µ) generally grow with µ, and must be subtracted the appropriate number of times. The subtractions introduce terms which are polynomial in k 2 and can be absorbed into the corresponding contact interactions. It was shown in Ref. [176] that even at fairly large internucleon distances, the potentials receive significant contributions from the spectral function in the region of µ Λ χ , where the chiral expansion cannot be trusted. It was, therefore, proposed in that paper to employ an ultraviolet cutoff Λ in the spectral integrals. This can be shown to be equivalent to introducing a particular ultraviolet cutoff in the loop integrals over the momentum q. Using a sharp cutoff Λ in the spectral integrals over µ leads to the following modification of the loop functions L(k) and A(k): where we have introduced l = Λ 2 − 4m 2 π . The resulting approach is referred to as the spectral function regularization. The limit of an infinitely large cutoff Λ corresponds to the previously considered case of dimensional regularization with L ∞ (k) = L(k) and A ∞ (k) = A(k). The spectral function regularization approach with a finite value of Λ was employed in the PCTC potentials of Ref. [74] and the more recent work [18], and in the derivation of the N 2 LO PVTC potential in Ref. [46]. • Local regularization in momentum space. The previously introduced spectral function regularization approach has an unpleasant feature of inducing long-range finite-Λ artifacts as can be seen by expanding the functions L Λ (k) and A Λ (k) in inverse powers of Λ. This feature may affect the applicability of chiral EFT for softer cutoff choices. Recently, local regulators in coordinate [179,17] and momentum space [19] were introduced, which do not affect the analytic structure of the pionexchange interactions and thus maintain the long-range part of the nuclear force. The approach of Ref. [19] amounts to replacing the static propagators of pions exchanged between different nucleons via 1 with q ≡ |q|. Such a regulator obviously does not induce any long-range artifacts at any order in the 1/Λ-expansion. This regularization approach can be easily implemented for two-pion exchange NN potentials with no need to recalculate the various loop integrals. Using the feature that the regulator does not affect long-range interactions, it is easy to show that the regularization of a generic two-pion exchange contribution simply amounts to introducing a specific cutoff in the dispersive representation (modulo short-range interactions), namely [19] 2 π ∞ 2mπ dµ µ ρ i (µ) In Ref. [19], the regularized two-pion exchange contributions were defined using the requirement (i.e. a convention) that the corresponding potentials in coordinate space and derivatives thereof vanish at the origin. This is achieved by adding to the right-hand side of Eq. (78) a specific combination of (locally regularized) contact interactions allowed by the power counting. For more details and explicit expressions see Ref. [19]. This local regularization scheme has not been used for PVTC or PVTV nuclear potentials.

The regularized PVTC potential
Once the loop integrals have been manipulated as discussed previously, we can now write the PVTC potential up to N 2 LO derived from χEFT. In the following, some of the LEC's have been further redefined to absorb the singular parts coming from the loop integrals. If one has chosen to regularize the loop integral using the SFR method, then the functions L(k) and A(k) below have to be substituted with L Λ (k) and A Λ (k), the spectral regularized functions, see Eq. (76). In summary, where V (−1) The NLO term V P V T C (l). Let us note that we have in total 11 LECs that must be determined from the experimental data: one in the LO term, five in the subleading order and five at N 2 LO. This potential is the same as the one derived using the MUT in Ref. [173].
Finally, the potential to be used in calculation of PVTC observables has to be regularized for large values of p, p ′ . The frequently used procedure is to multiply by a cutoff function containing a parameter Λ C Usual choices for f Λ C are [74] f where usually n = 6, adopted for example in Ref. [173], or adopted in Ref. [42]. The value of the cutoff Λ C is chosen to be around 400-600 MeV, and consistent with the analogous parameter used to regularize the PCTC potential.
The currently most accurate and precise PCTC NN potentials of Ref. [19] employ the local momentum-space regularization approach for pion-exchange contributions as described in section 3.4.1 in combination with a nonlocal Gaussian regulator given in Eq. (85) with n = 2 and Λ C = Λ for contact interactions (Λ is the cutoff used in the local regulator in Eqs. (77), (78)). The superior performance of the momentum-space regulator in Eq. (78) as compared with both the spectral-function regularization and a local multiplicative regularization as defined in Eq. (86) manifests itself in exponentially small distortions at large distances as visualized in Fig. 5 of [19].
Last but not least, we emphasize that using different regulators when calculating loop integrals in the nuclear potentials/currents and solving the Schrödinger equation to compute observables is generally incorrect. This issue becomes relevant at the chiral order, at which one encounters the first loop contributions to the 3N potentials and to the NN exchange current operators (i.e. at order Q 4 or N 3 LO in the PCTC sector) [180,181], which is beyond the accuracy of the calculations described in this review article. For more details and a discussion of a possible solution to this problem see Ref. [182].

Relevant PCTC and PVTC electromagnetic currents
Electromagnetic currents can be calculated in the χEFT expansion. For our purposes we require currents for the longitudinal asymmetry in radiative neutron capture on a proton target at thermal energies. As we deal with a real outgoing photon, the LO PCTC current is induced by the nucleon magnetic moment. At NLO there are contributions from the convection currents and one-pion-exchange currents proportional to g 2 A . At NLO the relevant currents become where κ 0 = −0.12 and µ v = 3.71 are the isoscalar and isovector anomalous nucleon magnetic moments. p j and p ′ j denote the incoming and outgoing momenta of nucleon j interacting with a photon of outgoing momentum q. The intermediate pions Ref. [173] used these currents in combination with N 3 LO χEFT potentials from Ref. [17] to calculate the total np → dγ capture cross section. Using just the LO currents gives a cross section of 305 ± 4 mb, which grows to 319 ± 5 at NLO. The remaining 4% discrepancy to the experimental cross section 334.2 ± 0.5, indicates that N 2 LO currents should probably be included.
A consistent calculation of PVTC observables as the photon asymmetry in the np → dγ radiative capture also requires the inclusion of PVTC currents. There is no one-body current in this case, as the anapole moment vanishes for on-shell photons [183]. As such, the leading PVTC currents arises from . one-pion-exchange currents where we stress the dependence on the PVTC pion-nucleon LEC h 1 π . Higher-order PVTC currents have not been developed.

The PVTV potential up to order Q
In this section, we discuss the derivation of the PVTV NN and 3N potentials at N 2 LO. The final expressions are given in terms of a sum of diagrams, which can be obtained either using the MUT [166,184,4], standard dimensional regularization [72] or the TOPT method [75]. In the following, we briefly report the derivation of the PVTV potential in the framework of TOPT approach.
The TOPT diagrams that give contribution to the NN PVTV potential up to N 2 LO (order Q 1 ) are shown in Fig. 4. We do not consider diagrams which give contributions only to the renormalization of the LECs. In this section we write the final expression of the NN PVTV potential V P V T V by having already taken into account the singular parts coming from loops. Note that for the PVTV potential the LO term is of order Q −1 as for the PVTC case. However, now there will be terms of order Q 0 , which will be denoted as NLO terms, etc. We have namely coming from OPE diagrams at LO, TPE at N 2 LO, three-pion vertices (3π) at NLO and at N 2 LO, and contact contributions (CT). From now on we defineḡ * 0 =ḡ 0 +ḡ 2 /3. In this case, we report here the final form of the potential, namely, the LECs appearing in the expressions below are the physical ones, having reabsorbed the various infinities generated by loops and diagrams like those shown in panels (1)-(3) of Fig.2.
• One pion exchange diagram. The OPE term, depicted in diagram (a) of Fig. 4, gives a contribution at LO, namely of order Q −1 , coming from the NR expansion of the vertices where there are an isoscalar, an isovector and an isotensor components. Contributions coming from the 1/M expansion are considered to be suppressed at least by four orders with respect to V (−1) Notice that the above LECsC 1 ,C 2 ,C 3 ,C 4 andC 5 have been redefined to absorb various singular terms coming from the TPE and 3π diagrams. It is possible to write ten operators which can enter V (1) P V T V (CT) at order Q but only five of them are independent as discussed in Chapter 4. In this work we have chosen to write the operators in terms of k, so that the r-space version of V  Fig. 4. This term has no isovector component, as shown for the first time in [73]. It reads where the loop functions L(k) and H(k) are defined in Eq. (68).
• Diagrams with three pion vertices The 3π-exchange term gives a NLO contribution through the diagram (i) of Fig. 4, where A(k) is given in Eq. (74). Additional contributions coming from diagram (i) deriving from the 1/M expansion of the energy denominators and vertex functions are here neglected since we count them as N 3 LO. The diagram in panel (j) of Fig. 4 contributes to V (3π) Note in Eq. (94) the presence of the c 1 , c 2 and c 3 LECs, which belong to the PCTC Lagrangian given in Eq. (13). In Eqs. (93) and (94),∆ is a renormalized LEC. The 3π PVTV vertex gives rise to a three body interaction through the diagram (k) in Fig. 4. The lowest contribution appears at NLO while at N 2 LO the various time orderings cancel out [75]. The final expression for the NLO of the 3N PVTV potential is, where k i = p ′ i − p. This expression is in agreement with that reported in Ref. [68,69].

The PVTV current
The PVTV current up to now has been considered coming from the LO one-body contribution where d p (d n ) is the proton (neutron) EDM. In the nuclear physics applications, it is customary to consider d p and d n as unknown parameters, although they in principle can be estimated in terms of the LECs entering the χEFT, as we have seen in Sect. 3.2.2. The complete derivation of PVTV two-body currents has not been completed. Partial results have been given in Refs. [143,73].

PVTC AND PVTV POTENTIALS IN PIONLESS EFT
In this section, we specifically focus on the contact few-nucleon interactions which enter the potentials in both chiral and pionless EFT formulations. We also discuss the expected hierarchy of the corresponding LECs as suggested by the large-N c analysis.

Effective Lagrangians
At distances much larger than the range of the interactions mediated by pions, the latter degrees of freedom can be integrated out of the effective theory, and the relevant effective Lagrangian can be written in terms of nucleon fields only, that interact through contact vertices.
At leading order these vertices involve one spatial derivative of fields, responsible for parity violation. Time derivatives can be eliminated by recursively using the equations of motion order by order in the low-energy expansion. This reflects our freedom in choosing the nucleon interpolating field, and amounts to a definite choice of the off-shell behavior of amplitudes. The theory can be formulated in terms of non-relativistic nucleon fields represented by two-component Pauli spinors N r (x). The relativistic 1/M corrections, which can in principle be worked out (see e.g. [185]) will be of no interest here. Relativistic covariance requires that the interactions depend on the relative momenta only (momentum-dependent "drift" corrections, which vanish in the center-of-mass frame of two nucleon systems, are part of the above mentioned relativistic corrections). Thus, gradients of nucleon fields in two-nucleon contact operators may only enter in the combinations where (ai ← → ∇ b) ≡ a(i∇b) − (i∇a)b and the factor i, meant to ensure the hermiticity, makes it odd under time-reversal.
Since the underlying mechanism of parity violation in the SM may induce ∆I = 0, 1, 2 transitions (at least to order G 2 F ), the effective Lagrangian will contain contact operators which transform as isoscalars or the neutral components of isovector and isotensors. In the two-nucleon case all these flavor structures are real, and therefore unaffected by the time-reversal operation, except for ( τ 1 × τ 2 ) z , which changes sign.

PVTC Lagrangian
Following the general considerations outlined above, there are ten possible structures entering the twonucleon contact Lagrangian in the PVTC case, The Fermi statistics of nucleon fields, together with Fierz's reshuffling of spin-isospin indices allow to establish linear relations between primed and unprimed operators, thus reducing the number of independent operators to five, so that the effective Lagrangian can be written as where C i are LECs. This Lagrangian is identical to that reported in Eq. (23). From this Lagrangian, one can derive the potential given in Eq. (60).
The five LECs are in a one-to-one correspondence with the possible S-P transitions in two-nucleon systems [45], namely 1 S 0 -3 P 0 (∆I = 0, 1, 2), 3 S 1 -1 P 1 (∆I = 0) and 3 S 1 -3 P 1 (∆I = 1). This may be shown explicitly by using the spin-isospin projection operators [45,186,187,188] normalized according to such that the operator (N T r P µ,α N r ) † creates a correctly normalized two-nucleon state with the appropriate spin-isospin quantum numbers. The relevant operators [188] O are related to the original basis via Fierz's transformations as follows, whence one can read the relation between the partial-waves projected LECs and the C i . The potential derived from the operators given in Eq. (103) has been often used in studies of PVTC observables. It is given explicitly as [37,39] V (1) where the five LECs Λ (...) ∆T are in one-to-one correspondence with C 1−5 . Explicitly where κ = 2Mm 2 ρ /f π Λ 2 χ .

PVTV Lagrangian
The T -odd sector is very similar (see also Ref. [189]): one starts with a list of ten redundant operators, and uses Fierz's identities to establish the linear relations, so that the Lagrangian only depends on five LECs, from which one can derive the potential given in Eq. (91).
The five S-P transition operators only differ from the T -even case by a factor of i, related to the original basis as follows,

Constraints from the large-N c limit
In 't Hooft combined large-N c and small coupling limit, with g 2 s N c fixed [190], QCD considerably simplifies, while maintaining many of the features of the actual theory, becoming a theory of stable hadrons. The baryons emerge as dense systems of many quarks, subjected to a mean field potential [191]. Nucleon-nucleon interactions exhibit in this limit a spin-flavor symmetry [192,193,194]. Indeed, due to the fact that nucleons carry definite spin and isospin of O(1), interactions inducing a change in either spin or isospin are suppressed relative to the dominant O(N c ) ones, that are either spin-isospin independent (∼ 1) or dependent on both (∼ στ ). The large-N c counting of momenta follows from the observation that the nucleon-nucleon scattering amplitude is in this limit a sum of meson exchange poles, each one depending only on the relative momentum transfer. The average relative momenta can only appear as relativistic corrections, which are suppressed by inverse powers of M ∼ O(N c ).
Apparently the resulting scaling laws do not conform with the operator identities (99) and (108) and seem to imply a dependence on the choice of operator basis. However, one can start with the redundant set of operators, pertinent to a theory of distinguishable nucleons, since the large-N c arguments outlined above are completely general and do not rely on the statistics of the interacting baryons (the only assumption is that they both carry spin and isospin of O(1)). As a result one obtains the large-N c scaling of the LECs in the PVTV contact potential, and in the PVTV one, Therefore we have only two leading LECs in the PVTC potential (C 2 and C 5 corresponding to ∆I = 0, 2 respectively) and only one in the PVTV potential (C 3 with ∆I = 1) [47,150]. This largely increases the predictive power for low-energy hadronic parity violation, and allows to put more severe constraints on the forthcoming experimental results. Notice however that the above results are obtained by simply projecting the Hartree Hamiltonian in the nucleon-nucleon sector. A consistent treatment would require to consider the induced effect on NN contact vertices of ∆ exchanges, since the latter are enhanced, in the large-N c limit, due to the degeneracy between nucleon and delta masses implied by the spin-flavor symmetry.
Moreover, for the PVTV case, this picture is obscured by the fact that the magnitude of the five contact LECs depends strongly on the particular type of the CP-violating source at the quark level. For example, the QCDθ term conserves isospin symmetry such thatC 3,4,5 are suppressed by powers of εǫ mπ compared toC 1,2 (see Table 1). Despite the possible 1/N c suppression ofC 1,2 compared toC 3 the former are still expected to dominate.

ONE-MESON EXCHANGE MODELS
In the past, a simple and rather efficient description of the strong PCTC NN interaction was obtained in terms of a sum of single meson exchanges [195,196]. These models started to be popular since the discovery of various meson resonances during the sixties. The potentials were generally constructed taking into account the exchanges of pions (J P = 0 − , m π = 138 MeV), η-mesons (J P = 0 − , m η = 550 MeV), and ρand ω-mesons (J P = 1 − , m ρ,ω = 770, 780 MeV), but clearly, the number of mesons to be included is somewhat arbitrary. This picture has been extended also to describe PVTC and PVTV interactions, simply considering single meson exchanges where one of the vertex is strong and PCTC, while the other violates P and conserves T or violates both P and T. Then, all the dynamics of such interactions is contained in a number of PVTC and PVTC nucleon-nucleon-meson (NNM) coupling constants. Table 2. Transformation properties of fermion bilinears with different elements of the Clifford algebra and various meson fields under hermitian conjugation (H), parity (P), and charge conjugation (C). Note that the pion and rho-meson fields are isospin triplets, a = 1, 2, 3.
One starts by writing the Lagrangian consistent of Yukawa-like NNM vertices, invariant under the proper Lorentz transformations, and either conserving or violating the discrete P, C, T symmetries. The building blocks of the Lagrangian are therefore nucleon bilinears multiplied by a meson field arranged so that Lorentz symmetry is satisfied. For the construction of the PCTC Lagrangian, one usually includes only isospin-conserving terms, however, for the PVTC and PVTV Lagrangians also isospin-violating terms have to be included as the underlying operators at the quark level are not necessarily isospin symmetric. A summary of the transformation properties of nucleon bilinears with different elements of the Clifford algebra and the various meson fields under hermitian conjugation (H), parity P, and charge conjugation C are reported in Table 2.
Using these properties it is not difficult to write the Lagrangians. For example, the strong L PCTC Lagrangian constructed with these mesons is given by (here we list only isospin -conserving terms) DDH [50] DDH [50] Table 3. Weak NNM couplings as estimated in Ref. [50]. All numbers are quoted in units of the value 3.8 × 10 −8 .
where q µ is the meson momentum 7 , π a , ρ µ a , η and ω µ are meson fields and g π , . . . PCTC coupling constants. Above, χ V and χ S are the ratios of the tensor to vector coupling constant for ρ and ω, respectively. Assuming vector-meson dominance [197], they can be related to the iso-vector and isoscalar magnetic moments of a nucleon (χ V = 3.70 and χ S = −0.12). Note that the pion and rho-meson are isospin triplets, therefore the fields have the isospin index a = 1, . . . , 3. Moreover, the rho-and omega-mesons have spin 1, and their fields correspondingly are vector fields with index µ = 0, . . . , 3.
Let us now consider the PVTC Lagrangian constructed in terms of the same mesons. In this case one has to take into account Barton's theorem [198], which asserts that exchange of neutral and spinless mesons between on-shell nucleons is forbidden by CP invariance, and therefore they cannot enter in a PVTC Lagrangian. Therefore only π ± , ρ and ω vertices need to be considered and the form of the PVTC effective Lagrangian is [131] where h 1 π , . . . are PVTC coupling constants to be determined. As discussed also in Section 3, where we focused in particular on the pion-nucleon PVTC constant h 1 π , attempts to estimate the magnitude of these couplings from the fundamental theory were reported in several papers [121,122,123,124,125,126,127]. In particular, in the DDH paper [50], the authors presented reasonable ranges inside of which these parameters were extremely likely to be found, together with a set of "best values" (see Table 3). Clearly, these values have to be considered as educated guesses in view of all the uncertainties of their evaluation. Of the seven unknown weak couplings h 1 π , h 0 ρ , . . ., there are estimates that indicate that h ′ 1 ρ is quite small [199] and this term was generally omitted, leaving PVTC observables to be described in terms of six constants. 7 More appropriately, these Lagrangian terms should be written in terms of four-gradients. For examplē where [, ] denotes the commutator.
The PVTV potential was derived in Ref. [202,203,204,142]. The momentum space version reads Also in this case, cut off functions can be applied in order to regularize the large k behavior of V PVTV . It is worthwhile to stress that the PVTV meson-exchange potential involves significantly more parameters than the LO PVTV chiral potential which depends in principle only on 4 LECsḡ 0,1 andC 1,2 , withḡ 2 ,∆, andC 3,4,5 appearing at subleading orders. While the meson-exchange potential can be mapped onto the short-distanceC i operators, the dynamics from the 3-pion∆ interaction is not captured in this way.

SELECTED RESULTS FOR VARIOUS PVTC AND PVTV OBSERVABLES
In this section we present a selection of results obtained with the chiral EFT potentials and currents described in Sect. 3 for various PVTC and PVTV observables. We will discuss first in the next three subsections the parity violation in the radiative neutron capture on the proton, then the longitudinal asymmetry in pp scattering, and subsequently the n-p and n-d spin rotations. Finally, in the last subsection, we present some results for the EDM of light nuclei. Our motivation to include these results in the review is mainly to establish benchmarks to help future applications. We include also a "minimum" analysis how the current experimental data constrain some of the values of the LECs entering the χEFT interactions.
Results obtained using the pionless EFT can be found, for example, in Refs. [36,37,39]. The mesonexchange potentials (in particular the DDH model) were used to analyze the results of several experiments of PVTC observables also in medium and heavy nuclei. For a summary of the obtained results, see, for example, Refs. [131,34,39]. Calculations of the EDM of light nuclei using the meson exchange potential were performed in Refs. [142,205,206].

Parity violation in radiative neutron capture on the proton
The radiative neutron capture on the proton np → dγ, where d denotes the deuteron and n a longitudinally polarized neutron, represents a very interesting process to study PVTC effects in nuclear physics. The longitudinal analyzing power for this process is defined as where dσ ± (θ) is the differential cross section for positive/negative helicity neutrons, and θ is defined as the angle between the neutron spin and the outgoing photon momentum. a γ has been targeted by several experiments in the last decades. A first nonzero signal was reported last year for incoming neutrons of thermal energies [207], a γ = (−3.0 ± 1.4 ± 0.2) · 10 −8 .
albeit only two standard deviations away from a null measurement.
The theoretical asymmetry is given by where X ℓ ( 2S+1 S J ) are reduced matrix elements (RMEs) either of electric (X = E) or magnetic (X = M) type, of multipolarity ℓ, and describing the EM transition from the n − p system in the scattering state 2S+1 S J [208].
Compared to the PVTC longitudinal analyzing power in proton-proton scattering discussed later, a γ comes with a big advantage. The initial neutron-proton state can be in a 3 S 1 state, such that the process is sensitive to 3 S 1 ↔ 3 P 1 transitions and thus depends on the LO PVTC NN potential. In chiral EFT the LO potential only depends on the LEC h 1 π , such that measurements of a γ provide a unique chance to pin down the value of this LEC, something which is much more difficult in proton-proton scattering, where the contribution of the LO potential vanishes. The disadvantage is that np → dγ is an electromagnetic process and therefore depends on P-conserving and P-violating electromagnetic currents.
As can be seen from Eq. (129), nonzero a γ requires an interference between electric and magnetic dipole currents. As such, just including the leading magnetic moment current in the presence of the LO PVTC NN potential leads to a vanishing result and NLO currents are necessary. There are then three relevant contributions that consist of interference between the isovector nucleon magnetic moment and 1. the one-body convection current in combination with the PVTC NN potential, 2. the two-body PCTC currents in combination with the PVTC NN potential, 3. the two-body PVTC currents.
Each of these contributions is sizeable: a 1 γ = (−0.27±0.03)h 1 π , a 2 γ = (−0.53±0.02)h 1 π , and a 3 γ = (0.72± 0.03)h 1 π where the theoretical error bands are obtained from cut-off variations in the strong NN potential and do not reflect uncertainties from higher-order contributions [173]. While these uncertainties are small on the individual contributions, they lead to a sizeable uncertainty on the total analyzing power [173] The cancellations between the different contributions are related to gauge invariance [209,208,173] and this explains the relatively large total theoretical uncertainty. While the electromagnetic currents given above are explicitly gauge invariant as they result from the gauge-invariant χEFT Lagrangian, explicit gauge invariance is lost due to applied regulator when solving the NN scattering and bound-state equations. Future calculations can probably reduce the uncertainty by using regulators that do not violate explicit gauge invariance, but such schemes have not been applied to PVTC processes. Alternatively, it is possible to apply the Siegert theorem to relate part of the electric dipole currents to the one-body charge density. Ref. [210] applied the Siegert theorem in combination with phenomenological strong potentials to calculate a γ finding a result in good agreement with the central value in Eq. (130). Such calculations however do not include an uncertainty estimate, for instance from missing transverse currents that are not included when applying the Siegert theorem. In this light, Eq. (130) can be interpreted as a conservative result. It would be interesting to redo the calculation of a γ in an updated framework to reduce the theoretical uncertainty.
The contribution to a γ of the short range components of the potential is considered to be negligible. For example, using the meson-exchange model, the calculations have shown that a γ is essentially unaffected by short-range contributions [211,212,213,208], represented in this case by ρ and ω exchanges. Within χEFT, a resonance saturation estimate of the short-distance LECs contributing to the asymmetry led to short-distance contributions to a γ of roughly 5 · 10 −9 and thus very small [173]. Therefore, considering the theoretical expression given in Eq. (129) and the experimental value given in Eq. (128), we obtain an estimate for the LEC h 1 Note that the large experimental error and the large theoretical uncertainty only allow to establish the positive sign and that the magnitude of this LEC is consistent with the Lattice QCD preliminary evaluation reported in Eq. (26) [132].

Parity violation in pp scattering
PVTC effects in proton-proton scattering can be studied by looking at the longitudinal analyzing power A z (E, θ) defined as, where θ is the scattering angle and E the energy of the protons in the laboratory frame, and σ + (θ, E)(σ − (θ, E)) the cross section when the polarization of the incoming proton is parallel (antiparallel) to the beam direction. Actually the experiments detect the particles scattered in angular range [θ 1 , θ 2 ] and the measured quantity is an "average" of the asymmetry over the total cross-section in this range, explicitly where is the unpolarized differential cross-section for the process. There exist three measurements of the angleaveraged pp longitudinal asymmetry A z (E), see Eq. (133), obtained at different laboratory energies E [214,215,216]. The measurements and the angle ranges are reported in Table 4.
The isospin state of two proton system is |pp ≡ |T = 1, T z = 1 , therefore the LO contribution that comes from the OPE vanishes and the LEC h 1 π will contribute to the observable only via the TPE box diagrams that appear at NLO and N 2 LO. Taking into account the isospin selection rules, the longitudinal asymmetry can be written as  Table 4. Values of A z and angle ranges for the three measurements of the pp longitudinal analyzing power [214,215,216].  [18] at the three energies corresponding to the experimental data points. The PVTC potential has been regularized as in Eq. (86) adopting the value Λ C = 500 MeV for the cutoff parameter. The PCTC potential has been regularized with the same value of the cutoff parameter. For the coefficient a where the first two terms are NLO contributions and the third term enters at N 2 LO. We defined and a  [18] are reported in Table 5. The only coefficient which receives contributions from both the NLO and N 2 LO potentials is a (pp) 0 . In the table, we report separately the two contributions and also the total contribution, given simply as The value of a (pp) 0 (N 2 LO) has been obtained assuming a value c 4 = 3.56 GeV −1 [217]. This correction to a (pp) 0 is of the order of ∼ 50% with respect to the NLO value, somewhat larger than expected. This is related by the unnaturally large value of the πNN LEC c 4 appearing in the PCTC Lagrangian (13). This value has been obtained from the Roy-Steiner analysis of πN scattering data at N 2 LO performed in Ref. [217].
Unfortunately, of the performed measurements, the two at the lowest energy do not give independent information. In fact, the observable A z at low energy scales as √ E, since its energy dependence in this energy range is driven solely by that of the S-wave (strong interaction) phase shift [218]. Because of this scaling, it is not possible to fit from these data all three LECs h 1 π , C, andh at the same time. If we fix the value h 1 π = 2.7 × 10 −7 from the central value as extracted from the np → dγ observable, see Eq. (131), then we can perform a χ 2 analysis of the three data points listed in Table 4 in order to fix the values of C andh. Note that this value of h 1 π was obtained from the np → dγ calculation performed in Ref. [173]  reported in Table 5 assuming the value h 1 π = 2.7 × 10 −7 .
using a different PCTC potential than that one used compute the a (pp) i coefficients . However, since the np → dγ experiment depends mainly on the peripheral regions of the process, the value of a γ is not very sensitive to the PCTC interaction (see also the calculations reported in Ref. [219]).
First of all, if we restrict ourselves to an NLO analysis, using h 1 π = 2.7 × 10 −7 we would obtain C = (49 ± 2) · 10 −7 . If we take into account also the N 2 LO LEC, we report in Fig. 5 the C andh values for which χ 2 ≤ 2, which form an elliptic region. As can be seen, there appears to be a strong correlation between C andh and the range of allowed values of the LECs is rather large 5 × 10 −7 < C < 67 × 10 −7 and −1.5 × 10 −7 <h < 2.5 × 10 −7 . Note that the ellipse is rather narrow and almost coincides with a straight line. See also Refs. [46,42] for a similar analysis performed at NLO for the LECs h 1 π and C only.
The previous discussion did not take into account the large uncertainty of the h 1 π coupling constant after the fit of the n -p radiative capture asymmetry. In Table 6, we report representative values of C andh giving the minimum value of χ 2 corresponding to range of values for h 1 π as given in Eq. (131). In the fourth column we report values for C if we neglect the N 2 LO contributions (settingh = 0). We conclude that the combination of the pp and np → dγ asymmetries allows for a rough extraction of the LO and NLO LECs h 1 π and C, but is insufficient to also pinpoint the N 2 LO LECh. The uncertainty of the extractions of h 1 π and C is dominated by theoretical and experimental uncertainties related to the PVTC asymmetry in the radiative neutron capture process.  Table 6. Values for C andh corresponding to different values of h 1 π (all LECs are given in units of 10 −7 ) giving the minimum value of the χ 2 in the fit of the three experimental p -p data points. For example, for h 1 π = 2.7×10 −7 , the C,h values are those lying in the center of the elliptical contour shown in Fig. 5. The fourth column corresponds to an analysis where we ignore the N 2 LO contributions and thus seth = 0.

The n-p and n-d spin rotation
The spin rotation of neutron traversing a slab of matter in a plane transverse to the beam direction induced by the PVTC potential is given by where ρ is the density of hydrogen or deuterium nuclei for X = p or d, | pẑ; m n , m X (±) are the n-X scattering states with outgoing-wave (+) and incoming-wave (−) boundary conditions and relative momentum p = pẑ taken along the spin-quantization axis (theẑ-axis), S X is the X spin, and v rel = p/µ is the magnitude of the relative velocity, µ being the n-X reduced mass. The expression above is averaged over the spin projections m X ; however, the phase factor ǫ mn = (−) 1/2−mn is ±1 depending on whether the neutron has m n = ±1/2. We consider the n-p and n-d spin rotations for vanishing incident neutron energy (measurements of this observable are performed using ultracold neutron beams). In the following, we assume ρ = 0.4 × 10 23 cm −3 . The rotation angle depends linearly on the PVTC LECs, as higher-order weak corrections are negligible. We write where the a The values of these coefficients for the n-p case and the cut-off value Λ = 500 MeV are listed in Table 7. From that Table, it is possible to appreciate the chiral convergence for the coefficients a P V T C (TPE) multiplies the operator (τ 1z + τ 2z ). The n-p spin rotation is sensitive to all the LECs except for the LECs C 4 and h 1 A multiplying again the isospin term (τ 1z + τ 2z ); in particular, there is a large sensitivity to C 5 and h 2 A , which multiply the isotensor terms of the PVTC potential.  Table 7. Values of the coefficients entering the expression of the n -p spin rotation in units of Rad m −1 calculated for the χEFT N 2 LO PVTC potential described in Sect. 3.4 and the N 4 LO PCTC potential derived in Ref. [18] at vanishing neutron beam energy. The PVTC potential has been regularized as in Eq. (86) adopting the value Λ C = 500 MeV for the cutoff parameter. The PCTC potential has been regularized with the same value of the cutoff parameter. For a  Table 8. The same as in Table 7 but for the n -d spin rotation and using the χEFT NLO PVTC potential and the N 3 LO PCTC potential derived in Ref. [5].
Regarding the n-d spin rotation, the coefficients, as reported in Table 8, are calculated by using only the NLO PVTC potential. We note the large sensitivity to h 1 π (this fact is well known [201,220]), and to the LEC's C 2 and C 3 .
At present there are no measurements of these quantities, however their experimental knowledge could be very useful in isolating certain combinations of LECs.

EDM of light nuclei
The EDM operatorD is composed by two parts, D PCTC is the electric dipole operator derived from the current J P CT C given in Eq. (87), after using the long wavelength approximation and the continuity equation [221], explicitlŷ where e > 0 is the electric unit charge, τ z (i) and r i are the z component of the isospin and the position of the i-th particle. This operator implicitly takes into account also the main part of the two-body PCTC currents. TheD PVTV contribution comes from the PVTV current at LO given in Eq. (96) and it readŝ where d p and d n are the EDM of proton and neutron, respectively and σ i is the spin operator which act on the i-th particle. As discussed in Sect. 3.5.1 and in Refs. [143,73] theD PVTV should also include contributions from transition currents at N 2 LO. These are not considered in this review.
The EDM of an A nucleus can be expressed as where |ψ A + (|ψ A − ) is defined to be the even-parity (odd-parity) component of the wave function. In general, due to the smallness of the LECs, the EDM depends linearly on the PVTV LECs where the a i for i = 0, 1, 2, A i for i = 1, . . . , 5, a ∆ , and a p , a n are coefficients independent on the LEC values (all coefficients except a p and a n have the unit of a length). For the deuteron, d 2 PVTV is dominated by one-body components, proportional to the neutron and proton EDM. The coefficients a p and a n multiplying the intrinsic neutron and proton EDM, as already pointed out first in Ref. [222] and then in Ref. [66], are given by, where P D is the percentage of D-wave present in the deuteron wave function. d 2 PCTC , in the case of the deuteron, receives contribution only from the LECsḡ 1 ,∆,C 3 andC 4 . The coefficients calculated with the χEFT N 2 LO PVTV potential described in Sect. 3.5 and the N 4 LO PCTC potential derived in Ref. [18] are reported in Table 9. The cutoff for both the PCTC and PVTV potentials has been chosen to be Λ C = 500 MeV. The coefficients a 1 , A 3 and A 4 agree well with the power counting expectation in Eq. (37). The slight suppression of a 1 compared with the naive estimate a 1 ∼ 1 is in very good agreement with the perturbative pion power counting [70]. The LO perturbative pion calculation of a 1 agrees with the value in Table 9 at the 20% level [70]. Results obtained in chiral EFT with N 2 LO PCTC potentials [66], and with "hybrid" approaches [143,222] based on chiral PVTV and phenomenological PCTC potentials, also agree well with the results reported in Table 9. The contribution of the three-pion coupling a ∆ is a bit more problematic. We find in this case that the contribution of the N 2 LO term is of the order of ∼ 60% of the NLO term. We will discuss the issue of these large N 2 LO corrections more in detail below.
Depending on the source of CP violation at the quark level, the deuteron EDM can be dominated by different LECs. For sources such as quark chromo-EDMs and four-quark operators Ξ, for whichḡ 1 is induced without any chiral suppression, the pion-exchange contribution proportional toḡ 1 is expected to dominate the deuteron EDM. For sources such as quark EDMs or the Weinberg operator, however, the deuteron EDM is well approximated by the sum of the nucleon EDMs. For the θ-term, the pion-exchange a n (a p ) 0  Table 9. Values of the coefficients entering the expression of the deuteron EDM calculated for the χEFT N 2 LO PVTV potential described in Sect. 3.5 and the N 4 LO PCTC potential derived in Ref. [18]. The PVTC potential has been regularized as in Eq. (86) adopting the value Λ C = 500 MeV for the cutoff parameter. The PCTC potential has been regularized with the same value of the cutoff parameter. For a ∆ we give explicitly the contribution of the different orders, the sum of the two contributions is given in the last row.  Table 10. The same as in Table 9 but for the 3 H and 3 He EDM.
contributions are expected to be minor as well. Given measurements of the deuteron and nucleon EDMs, one can, therefore, identify the underlying source of CP violation [223,70].
As regarding the 3 H and 3 He EDMs, the results are summarized in Table 10. The coefficients a 0 and a 1 are again a bit smaller than the O(1) expectation. Note that the value for a 0 reported in Table 9 is approximately 50% smaller than that reported in Ref. [66]. This difference can be traced back to the contribution of the TPE, which was not included in that work. Performing the calculations at LO, namely including only the OPE term, the a 0 coefficient results to agree with that reported in Ref. [66]. The values of the numerical coefficients are mostly equal in modulus between 3 H and 3 He except a p and a n . The coefficients associated to isovector terms have the same sign while all the others are opposite. Again the contribution of the N 2 LO potential term to a ∆ is significant, about 60%. This issue is discussed below.
Let us now consider in more detail the issue of the NLO and N 2 LO contributions to a ∆ . We have seen that in all cases the N 2 LO correction to a ∆ is of the order of 60%, a bit larger than expected. Explicitly, the coefficient a ∆ can be written as [75] a ∆ = a ∆ (NLO) + a ∆ (N 2 LO) ,  Table 11. Values of the various components coefficients a ∆ as given in Eq. (149) in units of e fm for the different nuclei. The coefficients have been evaluated using the N 4 LO PC potential derived in Ref. [18] and using a cutoff parameter of value Λ C = 500 MeV.
where a ∆ (0) comes from the NLO potential V  Tables 9 and 10, the following values were adopted: c 1 = −1.10 GeV −1 , c 2 = +3.57 GeV −1 , and c 3 = −5.54 GeV −1 as reported in Refs. [217,224]. The large N 2 LO corrections are caused by the large values of these LECs. 8 For more detail, see Ref. [75]. For the trinucleon systems, the values of a ∆ (3N) give a correction to a ∆ (NLO) of the order of ∼ 25%, which is in line with the chiral perturbation theory prediction because these contributions appear at the same order.
Similarly to the deuteron EDM, the trinucleon EDMs can be dominated by different terms. As the isoscalar interaction proportional toḡ 0 andC 1,2 now gives a sizable contribution, the trinucleon EDMs are noticeably different from the nucleon EDMs for the QCD θ-term, the quark chromo-EDMs, the four-quark operators Ξ, and potentially the Weinberg operator and the four-quark operators Σ. These EDMs therefore provide complementary information to the deuteron and nucleon EDMs. Combined measurements of all these EDMs would allow one to unravel various BSM models of new CP violation [71].

CONCLUSIONS AND PERSPECTIVES
In this paper we have discussed the current status of the PVTC and PVTV nuclear interactions using the traditional approach based on phenomenological boson exchange models and utilizing the modern frameworks of pionless and chiral EFT. The study of PVTC signals in nuclei is interesting since it derives from the nonleptonic weak interactions between quarks. Furthermore, a solid understanding of the manifestation of PVTC interactions at the nuclear level would give us confidence in the analysis of the more exotic PVTV case and other BSM nuclear observables. In fact, PVTV observables provide very valuable information since they are sensitive to interactions originating from the θ-term in the SM and even to more exotic mechanisms appearing in BSM theories.
As discussed in this review, the theoretical understanding of the PVTC and PVTV interactions is already rather advanced. Interactions in χEFT have been developed up to N 2 LO. The convergence of the χEFT appears to be problematic only for the contributions proportional to the ππNN LECs c i , due to the large values of those coefficients as it results from πN scattering data [217]. Given that the LECs c 2,3,4 are largely driven by the ∆(1232) [225], one may expect a better convergence in the formulation of chiral EFT that includes the ∆ as an explicit degree of freedom. Furthermore, large N c analysis may help in reducing the number of contact LECs. Also Lattice QCD calculations start to give valuable information [132,226]. 8 Notice that the values of the LECs c i obtained from the pion-nucleon amplitude at NLO, which would be appropriate for V (1) P V T V , are considerably smaller in magnitude. We, however, decided to adopt the larger values to be consistent with the employed PCTC potential.
We have also reported the results of the theoretical calculations of several observables performed using the potentials derived within the χEFT framework. The PVTC observables considered include i) the longitudinal asymmetry in n-p radiative capture, ii) the longitudinal asymmetry in proton-proton elastic scattering, and iii) the spin rotation of a neutron beam passing through a hydrogen and deuterium gas. As an example of PVTV observable, we have studied the EDMs of some light nuclei. The main motivation to study these observables is that for such light systems, the theoretical analysis can be carried out without invoking any uncontrolled approximations. Thus, the comparison with the experimental data can be performed unambiguously. The analyses of PVTC and PVTV observables using meson exchange models can be found in other review articles [131,34,39] and are not reported here.
As discussed previously, there exists a first measurement of the parameter a γ of the radiative neutron capture on the proton np → dγ. The large error derives from the smallness of this parameter which makes this measurement very awkward [207]. This observable is directly connected to the LO pionnucleon PVTC coupling constant a γ ∼ h 1 π . However, as we have seen, the theoretical estimate of the proportionality coefficient has been obtained with a relatively large theoretical uncertainty due to sizeable cancellations between different contributions. Therefore, to infer information from this observable, it will be necessary to make progress in both the experimental and theoretical analyses.
Other important information is brought forth by the three measurements at different energies of the p-p longitudinal asymmetry. This observable is sensitive to h 1 π via the TPE component of the PVTC potential and also to other LECs. In fact, owing to the isospin quantum numbers T = 1, T z = 1 of the p-p system, the LO contribution vanishes. Moreover, at NLO (N 2 LO), this observable depends on two (three) combinations of the LECs. Unfortunately, only two of the performed measurements give independent information. These two data have not been obtained with enough accuracy, so the constraints to the (combinations of) LECs which can be obtained are not so stringent [46,42], as discussed in Subsect. 6.2. For this observable the wave functions are easily obtained, however the vanishing of the LO contribution makes the χPT convergence more uncertain. On the other hand, it would be very useful to have more accurate experimental measurements.
Regarding the spin rotation observables, no experiments to measure the n-p and n-d spin rotation angles, which could provide useful information on some of the contact term LECs, are planned at present. The experimental detection of a non-vanishing n-p spin rotation would be rather important for two reasons: i) the theoretical treatment of the two-nucleon system does not present any difficulty numerically, while ii) this observable is sensitive to the LO term and therefore the chiral expansion of the potential is well under control, as discussed in Subect. 6.3. Regarding the n-d spin rotation, the same is not completely true since, as discussed in Sect. 3.4 one has to include also the PVTC 3N interaction terms which start to appear at N 2 LO. This is an interesting extension of χEFT which will be considered in future. From the experimental point of view, we note that there is an experiment in progress trying to measure the n-4 He spin rotation at NIST [227]. Some years ago there was a measure of the longitudinal asymmetry in p-4 He scattering, but this experiment has been performed at a rather high energy of the proton beam (46 MeV) [228] and this makes the theoretical treatment very difficult and impossible without some approximations.
Also a measurement of the n-3 He longitudinal asymmetry at the SNS facility is in progress, and the experimental value should be published soon. For this A = 4 system it is possible to perform accurate calculations of the wave functions [42], and therefore this observable may give valuable information in the near future.
To have the possibility to pin down all the LECs (6 LECs at NLO and 5 more at N 2 LO) more experimental information will be necessary in any case. In particular, an interesting possibility would be to measure PVTC observables in the A = 3 system, such as the longitudinal asymmetry of p-d elastic scattering and the photon asymmetry in n-d radiative capture. For both reactions, the theoretical treatment would be straightforward, once the PVTC 3N force has been taken into account.
Regarding the PVTV observables, the measurement of EDMs of particles is the most promising observable for studying CP violation beyond CKM mixing matrix effects. Currently, there are proposals for the direct measurement of EDMs of electrons, single nucleons and light nuclei in dedicated storage rings [77,78,81,229,82]. This new approach plans to reach an accuracy of ∼ 10 −16 e fm, improving the sensitivity in particular in the hadronic sector. Any measurement of a non-vanishing EDM of this magnitude would provide evidence of PVTV beyond CKM effects [52,55,56,57]. However, a single measurement will be insufficient to identify the source of PVTV, only the availability of the measurement of EDM of various light nuclei such as 2 H, 3 H, and 3 He can impose constrains on all the LECs. Other light nuclear EDMs have been discussed in Refs. [206,230]. EDMs of heavy diamagnetic systems provide very important information as well, but the systems are too large for chiral EFT calculations.
Other observables sensitive to PVTV effects are the transmission of polarized neutrons through a polarized target [231,232], in particular for heavy nuclei the PVTV effects can be enhanced by factors as large as 10 6 [233,234], see also Ref. [235]. In order to exploit this enhancement, some experiments are being planned, such as the NOPTREX experiment at RIKEN [236,237]. Also polarized nucleonpolarized deuteron scattering has been proposed as a way to detect PVTV signals [220,238]. Finally, searching for a large PVTV violations in polarized β-decay of 8 Li via measurements of the triple vector correlation is under consideration [239]. Clearly, it would be important to be able to detect a nonzero PVTV signal in all these experiments to be able to pin down the values of all the LECs.
From the theoretical point of view, calculations of the EDM of 2 H, 3 H and 3 He can be performed very accurately, also taking into account the contributions of the PVTV 3N force. The robustness of the calculation has been checked by evaluating the EDMs of the nuclei using different chiral orders in the PCTC potential. The discrepancy between the use of the N 2 LO and the N 4 LO PCTC potential has been found to be approximately 5% [75].
Currently, the only missing ingredient is the two-body PVTV N 2 LO currents [143,73]. Once this problem will be solved, one will achieve a fully consistent calculation of the EDM of light nuclei up to N 2 LO. There are also plans to perform theoretical studies of PVTV observables in n-p and n-d scattering in order to have independent and complementary information of PVTV effects.
The PVTV χEFT interaction developed in the previous sections depends on 11 coupling constants that should be determined by comparing with experimental data. As already pointed out by many authors [65,68,69] and discussed in Subsect. 3.2.1, the LECsḡ 2 ,C 3 ,C 4 andC 5 are suppressed for all CP-violation sources. However, for certain sources, the suppression is not too severe. For example, in Ref. [66], an analysis of the nuclear EDM in the minimal left-right scenario is presented in which the Lagrangian terms with LECsC 3 andC 4 appear at N 2 LO. In any case, since the CP-violation sources are not known, the only way to determine them is to fit all possible LECs and then compare the results with predictions for various scenarios.
Most of the observables discussed so far were obtained (or they are planned to be studied) at low energies, where also the pionless EFT framework is valid. The advantage of this framework is related to the fact that the resulting potentials depends on five LECs. Then, assuming the validity of the large N c analysis [47,49], the number of dominant LECs could be further reduced. This new paradigm is advocated for the PVTC case in Ref. [39]. For this case, only two LECs are expected to be dominant, the other three demoted to be subleading. Unfortunately, the photon asymmetry of np → dγ depends on the subleading LECs (this could explain its relative smallness) and therefore cannot be used to give information on the two leading LECs. Moreover, only one of the measured p-p longitudinal asymmetry at low energy may be used to test if this hierarchy is realized in Nature (the other measurement is taken at too high energy to be used in the pionless EFT framework). The other observable which could give valuable information could be n-3 He longitudinal asymmetry, for which experimental data should appear soon. However, no theoretical calculations of this observable performed in the framework of pionless EFT are available at present. Additional information could be obtained by calculations of these LECs using Lattice QCD, presently in progress. Regarding the PVTV observables in pionless EFT, here the large N c analysis predicts that only one of the LECs should be dominant, while the other four being suppressed. However, this picture is partially obscured by the fact that the magnitude of the five contact LECs would depend very much on the particular type of the CP-violating source.
In conclusions, the study of PVTC and PVTV observables is an active area of research that provide important tests of the SM and hopefully future evidence for BSM physics.

CONFLICT OF INTEREST STATEMENT
The Authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.