CP conservation in the strong interactions
Wen-Yuan Ai Bjorn Garbrecht Carlos Tamarit









empty




KCL-PH-TH/2023-51



MITP-24-043

TUM-HEP-1508/24








CP conservation in the strong interactions








Wen-Yuan Ai,^a Björn Garbrecht^b and Carlos Tamarit^c

0.3cm
^aTheoretical Particle Physics and Cosmology, King’s College London,
 Strand, London WC2R 2LS, UK

 ^bPhysik-Department T70, Technische Universität München,

James-Franck-Straße, 85748 Garching, Germany

^cPRISMA+ Cluster of Excellence and MITP,
Johannes Gutenberg University, D-55099 Mainz, Germany


1.4cm







We discuss matters related to the point that topological quantization in the strong interaction is a consequence of an infinite spacetime volume. Because of the ensuing order of limits, i.e. infinite volume prior to summing over topological sectors, CP is conserved. Here, we show that this reasoning is consistent with the construction of the path integral from steepest-descent contours. We reply to some objections that aim to support the case for CP violation in the strong interactions that are based on the role of the CP-odd theta-parameter in three-form effective theories, the correct sampling of all configurations in the dilute instanton gas approximation and the volume dependence of the partition function. We also show that the chiral effective field theory derived from taking the volume to infinity first is in no contradiction with analyses based on partially conserved axial currents.






.85cm




§ INTRODUCTION


The strong interactions conserve charge–parity (CP). This has been established through many observations, to the greatest precision in searches for the permanent electric dipole moment (EDM) of the neutron <cit.>, and calls for an explanation of certain theoretical aspects.

Since only a few years after the discovery of quantum chromodynamics (QCD) being the theory of the strong interactions <cit.> it has been suggested that due to the generic presence of a topological term in the action, the charge–parity symmetry CP should be violated <cit.>. To the present day, there has been no such observation, so arguably the coefficient θ  of the topological term (more precisely θ̅ as we introduce below) would need to be zero, corresponding to unnatural tuning and therefore a problem.

The present paper adds to the discussion following Ref. <cit.>. There, it has been brought up that the effective decomposition (quantization) of gauge field configurations of finite action into topological sectors of integer winding number without imposing ad hoc boundary conditions can only be derived when taking the volume of Euclidean spacetime to infinity. In contrast, there is no physical motivation for fixed boundary conditions on a finite surface in Euclidean space. As a consequence, the limit of infinite spacetime volume must be taken before summing over topological sectors, and it turns out that CP remains conserved this way.

One of the main points of the present paper is to extend this line of reasoning by demonstrating that this order of limits corresponds to a well-defined integration contour in the path integral constructed from the steepest-descent flows. In contrast,  the opposite conventional order of limits results in an integration that is inequivalent in the sense of the Cauchy theorem. Another purpose of this article is to address some criticisms regarding  Ref. <cit.>. This article also partly serves as a review of the related topics, though without including a complete list of references. For some other reviews, see e.g. Refs. <cit.>.

Without going into extensive technical detail, the basic reasons for CP conservation in the strong interaction are given in the summary of Section <ref>. There, we refer to the sections in the remainder of the paper where the statements of Section <ref> are supported at a more technical level.





§ SUMMARY AND OUTLINE




Strong interactions are described by a Yang–Mills theory which generally involves CP-odd parameters through the masses of the quarks as well as through the topological term. The Lagrangian in the Euclidean spacetime is 

    L=1/2g^2 trF_μνF_μν+ψ̅(γ̂_μ D_μ+ m  e^ iαγ_5)ψ
    -
    /16π^2θ  tr F_μνF_μν ,

where we use the convention Tr (T^a T^b)=δ ^ab/2, [T^a,T^b]= f^abcT^c for the Lie algebra generators T^a and the structure constants f^abc. Above, F_μν=F_μν^a T^a with F_μν^a≡∂_μ A^a_ν-∂_ν A_μ^a+f^abcA^b_μ A^c_ν being the field strength tensors. F^a_μν≡1/2ε_μναβ F^a_αβ (ε_1234=1) is the Hodge dual of F^a_μν. The covariant derivative takes the form

    D_μψ_i=(∂_μ- A_μ^a T^a)ψ_i

when ψ_i lives in the fundamental representation of the gauge group and 

    D_μψ_i=∂_μψ_i- A_μ^a[T^a,ψ_i]

when ψ_i lives in the adjoint representation. 

The Euclidean gamma matrices γ̂_μ are obtained from the Minkowskian counterparts

    γ^μ=[    0  σ^μ; σ̅^μ    0 ] ,

(where σ^μ=(1_2,σ⃗) and σ̅^μ=(1_2,-σ⃗) with σ⃗^i the Pauli matrices)
via 

    γ̂_4=γ_0=γ^0 ,  γ̂_i=γ_i=-γ^i .



These matrices satisfy the Clifford algebras

    {γ^μ,γ^ν}=2g^μν1_4 ,    {γ̂_μ,γ̂_ν}=2δ_μν1_4 .

Following Ref. <cit.>,
we use the same γ^5 in Euclidean spacetime as for Minkowski spacetime[In Ref. <cit.>, it is mistakenly stated that γ^5=γ̂_1γ̂_2γ̂_3γ̂_4.]

    γ^5≡γ^0γ^1γ^2γ^3=-γ̂_1γ̂_2γ̂_3γ̂_4 .

Since in this paper we mostly work in Euclidean space, from now on we remove the hat on the Euclidean gamma matrices. Note that the γ^5 in the Euclidean spacetime used here may differ from that used in some other papers, e.g. Ref. <cit.>, by a minus sign. The chirality in Euclidean spacetime in these two notations is thus defined oppositely. The only effect of this appears in applying the Atiyah–Singer index theorem <cit.> when e.g. deriving the anomalous axial current by counting the zero modes of the massless Dirac operator.

We note that the terms tr F_μνF_μν and ψ̅ iγ^5 ψ are both parity-odd and charge-conjugation even <cit.>. The CP-odd parameters are therefore α, the phase pertaining to the mass m of the fermion ψ, and θ, the coefficient of the topological term. The fermion ψ in the fundamental representation is referred to as quark. To focus on the principal aspects, for the most part of the discussion we take here just one quark flavour and the gauge group to be SU(2). This is the minimal setup that allows to study the interplay of the CP-odd parameters α and θ. In particular, we are interested in how  α and θ appear in the effective interaction (known as `t Hooft operator) that captures nonperturbative effects associated with the chiral anomaly.  In the single-flavour model, this operator has the same form as a quark mass term, and therefore both would be hard to discern. Nonetheless, the main question of whether and how α and θ appear in the effective `t Hooft operator can be answered within this setup. The generalization to the phenomenologically relevant case of several flavours is presented in Ref. <cit.>. In the strong interactions, the group SU(2) can be viewed as embedded within the SU(3) colour group. The technical details of this construction are reviewed in Ref. <cit.>.

It is well known that the presence of a CP-odd Lagrangian term does not readily imply CP-violating physical effects. A necessary condition is the existence of a CP-odd combination of the Lagrangian terms that is invariant under field redefinitions. In the present case, this condition is met as the parameter

    θ̅=θ+α

is invariant under redefinitions of the quark fields, in particular through anomalous chiral transformations. The chiral anomaly <cit.> also implies that θ is an angular variable, i.e. all observables must be 2π-periodic in θ. Therefore, the integral  1/(16 π^2)∫ d^4 x  tr FF that multiplies θ must be an integer in order to contribute to the action and thereby to the partition function.


Here, we ask the question of whether a nonvanishing θ̅ is also a sufficient condition for CP violation in strong interactions. Does θ̅ have physical effects, in particular, is there a neutron EDM depending on its value? Evidently, θ̅ has no impact on the classical equations of motion since the topological term is a total derivative. Nonetheless, under certain assumptions, based on the fact that the third homotopy group of the gauge group or one of its SU(2) subgroups is π_3( SU(2))= Z, the energy functional is periodic under so-called large gauge transformations. The situation is therefore reminiscent of a periodic quantum mechanical potential in a crystal, and θ̅ would then correspond to the crystal momentum <cit.>.

One way to see how far the analogy goes is to study canonical quantization of the gauge theory <cit.>. Since non-Abelian gauge theories are typically handled through functional quantization this possibility has not yet been investigated in all aspects pertinent to the present questions. We briefly comment on canonical quantization in Section <ref> and in more detail in a separate paper <cit.>.

For the time being, we focus on functional quantization since it has been the principal method to carry out calculations on CP violation in the strong interactions ever since the matter was brought up <cit.>. In the functional approach, we take the partition function

    Z[η,η̅]=
    ∫ Dψ̅  Dψ  DA  e^-lim_Ω→∞∫_Ω d^4x ( L-η̅ψ-ψ̅η)

as the defining point of the quantum field theory, where A is the gauge potential. We write this as a functional of external fermionic sources η(x), η̅(x) as a provision in order to derive quark correlation functions that may or may not exhibit CP invariance. Furthermore, we make explicit that the integral over Euclidean spacetime is understood as a limiting procedure, taking the spacetime volume Ω to infinity. It is this limit that allows us to state Eq. (<ref>) without specifying boundary conditions on the path integral and that, moreover, this way we obtain the vacuum correlation functions of the theory <cit.>. We shall review this point in Section <ref>.

Now, as argued in Ref. <cit.>, the partition function (<ref>) in infinite spacetime volume Ω→∞ receives its nonvanishing contributions from saddle points of finite action and fluctuations about these. For these field configurations, the winding number Δ n is an integer that labels the topological sector:

    Δ n=lim_Ω→∞∫_Ω d^4x 1/16π^2 tr FF∈ Z  for nonvanishing contributions to Z .

This is the desired outcome because it is consistent with θ being an angular variable, as required by the chiral anomaly. As topological quantization, i.e. integer Δ n, is a consequence of Ω→∞, we must carry out this limit before summing over topological sectors.

Suppose now that it is valid to organize the calculation of the path integral by adding contributions from the individual topological sectors. Then, Eq. (<ref>) implies that the partition function should be evaluated as

    Z[η,η̅]=
    ∑_Δ n∫ Dψ̅  Dψ  DA_Δ n e^-lim_Ω→∞∫_Ω^4x ( L-η̅ψ-ψ̅η) .

The subscript Δ n on DA_Δ n indicates that the path integral is supposed to cover the configurations with given Δ n, i.e. to sweep over the given topological sector. The contributions of the topological sectors are evaluated with the limit Ω→∞ before they are added together. A rearrangement of the limits will in general lead to different results and is therefore not justified.

In Section <ref>, we expand on this argument and put it on a more formal footing. The path integrals within individual topological sectors correspond to steepest-descent contours for the exponent of the Euclidean path integral. For different Δ n, these contours can only be connected by configurations of infinite action. It thus follows that the arrangement of limits in Eq. (<ref>) indeed corresponds to a good contour for the path integral in Eq. (<ref>). Further, this formally establishes that the decomposition of the path integral into topological sectors is valid in the first place.


This brings us to the salient point: In Ref. <cit.>, it has been shown that the limit Ω→∞ does not commute with the sum over Δ n in Eq. (<ref>), with the consequence that the quark correlations do not exhibit CP violation. While we review this technical argument in Section <ref>, the basic reason is that  lim_Δ n→±∞lim_Ω→∞Δ n/Ω=0, see also Section <ref>. As explained in Section <ref>, one can thus conclude that there is no EDM for the neutron, no matter what the value of θ̅ is. Since in the sum over Δ n all integer values are taken, this is where the analogy with the quantum-mechanical crystal breaks down as for the latter the number of potential minima may be large but remains finite. Therefore, the order of the path integral and the limit of an infinite spacetime volume is not an issue in that case.




After all, the strong interactions are complete without the necessity of tuning the parameter θ̅ to be small or extending the theory by additional scalar fields and nonrenormalizable operators. While this is a gratifying conclusion, a scrutiny of the argument is warranted, not least because the prevalent line of reasoning arrives at the contrary verdict: In order to deduce CP-violation in the strong interactions, one would have to impose that the limit Ω→∞ is taken last, i.e.

    Z[η,η̅^†]?=lim_Ω→∞∫ Dψ̅  Dψ  DA  e^-∫_Ω^4x ( L-η̅ψ-ψ̅η) ,

and at the same time specify boundary conditions on the finite surfaces ∂Ω  such as

    A_μ= iω∂_μω^-1 on some finite ∂Ω ,

where ω(x)∈ SU(2), which corresponds to a pure gauge and implies topological quantization, i.e. Δ n ∈ Z. Equations (<ref>) and (<ref>) together are either directly or indirectly implied in the bulk of the existing literature, including the initial papers on the topic <cit.>. Although Eq. (<ref>) may be motivated by considering fields in the classical ground state, i.e. of vanishing classical energy, on the initial and final spatial hypersurfaces, the quantum ground state also receives contributions from other field configurations. Therefore, the boundary condition (<ref>) does not follow from the partition function (<ref>) (note that Eq. (<ref>) moreover assumes that three-dimensional space is finite), unlike the topological quantization (<ref>) that is implied by the partition function (<ref>), see Section <ref> for more detail.



To our knowledge, the published papers do not provide a conclusive reason for how the procedure given by Eqs. (<ref>) and (<ref>) might be deduced from the functional (<ref>) that defines the theory. Given that the limits do not commute, i.e., that Eq. (<ref>) is not equivalent with Eqs. (<ref>) and (<ref>), as shown explicitly in Section <ref>, there also cannot be such a derivation. (Note that while in Ref. <cit.> it is shown that the correlators in a large but finite box depend only on the boundary conditions through the Chern–Simons flux, the issue that the limits do not commute is not addressed in that work.) Neither are we aware of an argument why Eq. (<ref>), which is the standard textbook expression  (up to the fact that we write the infinite-spacetime limit explicitly, which is a purely notational matter), might be incorrect to start with. Unless taking Ω→∞ there is also no apparent reason why the boundary condition (<ref>) should be physical.



0.9
0.95
In the simplest terms, the reason for CP-conservation in the strong interactions can thus be stated as follows:


  * For θ to be physical, we must have Δ n∈ Z since θ is an angular variable. Since fixed boundary conditions on a finite surface are not physical, this topological quantization can only follow from Ω→∞.

  * Given the order of limits that is thus implied, there is no CP-violation for Ω→∞ since lim_Δ n→±∞lim_Ω→∞Δ n/Ω=0.






While the main conclusions and technicalities on the absence of CP violation in the strong interactions have been presented in Ref. <cit.>, one objective of the present paper is to add a more formal interpretation of the difference between Eqs. (<ref>) and (<ref>) versus Eqs. (<ref>) and (<ref>). In Section <ref>, we recall to that end the reason for taking Euclidean time to infinity in the first place. Then, we show that Eq. (<ref>) corresponds to a contour integration that can be derived and assembled from steepest-descent flows, while the prescription of Eqs. (<ref>) and (<ref>)  does not correspond to a connected integration contour. Since the reasoning in the present work is based on infinite Euclidean spacetime as the analytic continuation of Minkowski spacetime, we briefly comment in Section <ref> on how calculations in finite Euclidean spacetimes with and without boundaries can be made meaningful, but we leave a detailed discussion to a separate paper. Next, as it allows for an explicit demonstration of CP conservation as a consequence of Eq. (<ref>) and for an intuitive interpretation of the matter, we review in Section <ref> the dilute instanton gas calculation of Ref. <cit.>. This lays the ground to address objections concerning the volume dependence of the partition function in Section <ref>. Since there is no physical interpretation of fixed boundary conditions on finite Euclidean surfaces, we show that the partition function in fact shows the expected behaviour when evaluated in finite volumes with open boundary conditions. Further objections are based on effective field theory (EFT) descriptions, which is why we review the role of the θ parameter in the effective `t Hooft vertex as well as in chiral perturbation theory in Section <ref>. With this preparation, we can reply in Section <ref> to an objection using the topological term in hadronic matrix elements and the role of θ in the effective description of the dynamics of the topological current. After some additional comments on the sampling of topological configurations in the different orders of limits and the recent literature, we wrap up and conclude in Section <ref>.
Except for Section <ref> where we work in Minkowski spacetime, we work in Euclidean spacetime throughout all other sections.




§ PATH INTEGRAL AND TOPOLOGICAL QUANTIZATION





  
Euclidean partition function for infinite volume
As a defining point of the quantum field theory, one may take the partition function (<ref>). Such a partition function based on an infinite Euclidean volume is a common starting point for a wide range of calculations. Here, we go for variations with respect to the sources η(x) and η̅(x) that yield Euclidean correlation functions for the quarks. 
Eventually, these can be analytically continued to Minkowski spacetime and interpreted as vacuum correlations.

Nonetheless, as some of these matters are contested in the context of CP conservation in the strong interactions, we shall briefly revisit here the reasons for taking an infinite Euclidean spacetime volume, or, more precisely, why we take imaginary time to infinity. For more discussion, see Ref. <cit.>.  In short, it allows one to obtain vacuum correlation functions without specifying the vacuum in terms of a wave functional (which appears to be practically impossible in quantum field theory at the nonlinear level).

Imaginary time arises from the analytic continuation of real time. The corresponding Wick rotation is straightforward for any spacetime that is stationary in real time. Note that taking the spacetime volume Ω to infinity comes as a consequence of the limit of infinite Euclidean time. The spatial geometry is not the decisive reason for this even though we take here unbounded Cartesian space R^3, for definiteness. The reason for taking time to infinity is as follows: By Wick rotation, the correlations derived from Z for infinite imaginary time correspond to the analytically continued expectation values for the state of the lowest energy that is accessible given the conservation laws. This remains true also for clockwise rotations of the real time axis by an angle 0<ϑ≤π/2 in the complex plane provided the infinite time limit is applied prior to taking ϑ to zero. (Do not confuse ϑ here with the angle θ as the coefficient of the topological term.) That is, for given initial and final states |i,f⟩ (which here may be taken as the usual linear combinations of field eigenstates of integer Chern–Simons number to comply with gauge invariance under large gauge transformations, i.e. the so-called θ-vacua) the path integral corresponds to

    Z_ϑ=lim_t→∞⟨ f| e^- iH t exp(- iϑ)| i⟩ ,

where H is the Hamiltonian. It therefore projects on the lowest energy eigenstate that is accessible. For ϑ=π/2, one gets the Euclidean path integral. In this sense, the Minkowskian vacuum correlation functions are analytic continuations of the Euclidean ones with a discontinuity on the real time axis.

If instead we were keeping time purely real or still complex in general but finite, then we would have to weigh each path contributing to the partition function by the ground state wave functional Ψ evaluated at the endpoints of these paths, as we discuss in Section <ref>. The main reason for taking complex time to infinity is therefore to avoid this complication because then we can evaluate the path integral without explicit use of the vacuum state. In particular, since the ground state wave functional in Yang–Mills theory is not known to an approximation that addresses the present purpose, taking the imaginary part of time to infinity is the main method of making the analytic approximation of vacuum correlators feasible.



  
Evaluation of the partition function and topological quantization
Having specified the problem through the partition function (<ref>), we now turn to its evaluation. In the following, we show that Eq. (<ref>) implies an integration contour that joins together the steepest-descent paths for each topological sector Δ n. The field configurations on these steepest-descent contours are not bound by any finite spacetime volume, i.e. there is no value R so that F_μν(x) F_μν(x)=0 (at least to some approximation) for |x|>R for all field configurations on a given steepest descent. Therefore, we must evaluate the path integrals for the different topological sectors Δ n in the limit Ω→∞ before interfering these <cit.>. This also implies that organizing the evaluation of the path integral in terms of a sum over contributions from the individual topological sectors is valid in the first place. Further, we argue that the order of limits in Eqs. (<ref>) and (<ref>) is opposite to what is implied by the form of the correct integration contour.

To start, we state why it is necessary to specify an appropriate integration contour. Since the integrand in Eq. (<ref>) for θ̅≠0 2π is not positive definite, explicitly because of the topological term and implicitly because of fermion determinants, we must determine the integration contours to leave the integral well defined. These specify the order in which the integral ∫ DA must be carried out so that we can derive how to sum over the topological sectors. In particular, it will allow us to discern whether it is Eq. (<ref>) or Eqs. (<ref>) and (<ref>) that correspond to the correct procedure.

To determine the contours, we first note that in general, the real part of the Euclidean Yang–Mills action 

    S⊃ S_ YM=1/2g^2∫_Ω d^4 x  trF_μνF_μν

should be bounded from below on its domain and thus have global or local minima that may or may not exhibit degeneracies. If Ω is infinite, this implies that the vector potential at such minima reduces to pure gauge configurations at infinity, i.e.

    A_μ(x)→ω(x) ∂_μω^-1(x) for |x|→∞  for local minima of  Re[S] ,

where ω(x)∈ SU(2). Because the surface at infinity is homeomorphic to S^3, and the third homotopy group of the gauge symmetry is π_3( SU(2))= Z, this immediately implies topological quantization as in Eq. (<ref>) for these minima, i.e. integer winding number Δ n. In addition, for each of these minima, there are degeneracies, i.e. flat directions of the action, parameterized by moduli <cit.>. Note that this reasoning applies without imposing boundary conditions ad hoc (as in Eq. (<ref>)) because the limit Ω→∞ appears inside the path integral (<ref>). This is in line with standard introductions of path integrals in infinite spacetimes that do not impose particular boundary conditions, see e.g. Ref. <cit.> where this point is mentioned explicitly.




Given the minima of the real part of the classical Euclidean action with integer winding number, the contributions to the path integral for the individual topological sectors Δ n can then be evaluated on steepest-descent contours 
(of the negative action -S)
passing through these minima. The contours are Lefschetz thimbles and are determined through flow equations <cit.>. By this reasoning of steepest-descent contours, upon dealing with the usual ultraviolet divergences and the vacuum contributions in the infinite spacetime volume, the path integrals in the individual topological sectors are convergent. The steepest-descent contours for different Δ n do not intersect for finite S because, in the infinite spacetime volume, solutions of different winding number are separated by infinite action barriers. The integration contours over the different sectors Δ n can therefore only be connected via configurations of infinite action S that give no contribution to the partition function (<ref>). Note that these infinite action configurations connecting the sectors are allowed precisely because we do not impose boundary conditions when evaluating the partition function. For comparison, in finite volumes with fixed boundary conditions (<ref>) the field configurations in different topological sectors are not continuously connected, not even via a path through configurations of infinite action since paths with a noninteger winding number are forbidden by the boundary conditions (<ref>), no matter whether S is finite or infinite. In Figure <ref>, we schematically illustrate these integration paths and their crucial differences.

Therefore, we first carry out the integration over the entire infinite spacetime in a given topological sector Δ n individually. Then, we can connect this integration with the steepest-descent contour for a different Δ n via configurations of infinite action that do not contribute to the path integral. This determines the integration contour that should be used in order to evaluate Eq. (<ref>) and that is given by the order of limits specified in Eq. (<ref>). Any alternative contour must be a continuous deformation in compliance with the prerequisites of the Cauchy theorem, a criterion that is not met with Eqs. (<ref>) and (<ref>), which consequently lead to a different result. In practice, this means that, without specifying ad hoc boundary conditions, we must not interfere the different sectors before taking Ω→∞. Otherwise, we would partition and rearrange the full integration contour in a non-continuous way that in general leads to an inequivalent result because the integrand, or more specifically, the sum over the topological sectors, is not positive definite and not absolutely convergent.

We can therefore write the path integral as in Eq. (<ref>) and as indicated in Figure <ref>. We emphasize that the decomposition into topological sectors follows from Ω→∞ and is not a consequence of the saddle point approximation (which only comes into play in the dilute instanton gas approximation in Section <ref>). In turn, finite surfaces ∂Ω imply that topological charge is not quantized, i.e. it can flow in and out of the volume Ω, unless imposing unphysical constraints.

Since the winding number density Δ n/Ω apparently is a measure of CP violation, it is already clear that for each term of the series (<ref>), there is no CP-violating contribution, as lim_Ω→∞Δ n /Ω=0. For the quark system, the corresponding calculation is explicitly presented in Section <ref>. 





  
Commuting the order of limits
The decisive point in the present discussion regarding the CP symmetry of the strong interactions is whether Eqs. (<ref>) and (<ref>) are consistent with Eq. (<ref>) and the integration contour that it implies. As we shall review in Section <ref>, the conclusion that the CP symmetry is violated when α+θ≠0π relies on a partition function as in Eq. (<ref>) in conjunction with the ad hoc boundary conditions (<ref>). In Eq. (<ref>) the order of limits is therefore opposite to Eq. (<ref>) that we have derived from the starting point given by Eq. (<ref>).

Now to further (beyond the apparent contradiction with Eq. (<ref>) regarding CP) assess the validity of Eqs. (<ref>) and (<ref>), we take  Ω to be finite and first assume that the boundary conditions from Eq. (<ref>) are not imposed. In particular, one can then move topological charge (i.e. instantons in the weak coupling limit) across the boundary ∂Ω so that there is no topological quantization into sectors with integer winding number Δ n and also no conservation of topological charge within Ω. Without these, nontrivial minima of the action do not exist also. To bring topology back into the picture, one might therefore impose the boundary conditions given in Eq. (<ref>). These boundary conditions can be derived as a consequence of the 

infinite spacetime volume Ω→∞ in the partition function (<ref>) <cit.>. However, imposing these boundary conditions and computing the partition function according to  Eqs. (<ref>) and (<ref>) requires commuting the limit of infinite spacetime volume with the sum over infinitely many topological sectors which is however not justified.

Equations (<ref>) and (<ref>) therefore do not follow from Eq. (<ref>) so that one may keep looking for alternative arguments. However, while topological quantization also emerges for fixed boundary conditions on a given compact ∂Ω as in Eq. (<ref>), there does not appear to be a valid reason for imposing such configurations. In fact, the vacuum wave functional has nonvanishing support on configurations that do not observe Eq. (<ref>) and thus have a nonvanishing classical energy. Since the field operators do not commute with the Hamiltonian, the configurations obeying Eq. (<ref>) are as good or as bad as any other field configuration subject to a different boundary condition on ∂Ω. (For the boundary conditions given in Eq. (<ref>), Δ n are integers whereas, for more general boundary conditions specified up to gauge transformations, Δ n are given by a fixed real number plus any integer.) So there is no preference for choosing pure gauges as a boundary condition on a sequence of finite surfaces, even as these surfaces are taken to infinity. For example, if ∂Ω is spherical, instantons could be placed at certain angles and close to the radius of ∂Ω, what defines boundary conditions with noninteger Δ n. Except that there is not a valid reason to impose Eqs. (<ref>) and (<ref>), whose consequence for the outcome of the calculation is material, there also is no a priori justification why Ω should be taken to infinity at the same rate for all topological sectors in Eq. (<ref>).

Finally, note that for examples that do not involve noncommuting limits, correlations from fixed boundary conditions on finite ∂Ω converge in imaginary time to the vacuum correlators as Ω→∞. However, this does not imply by analogy that in the present case, where we must sum over infinitely many topological sectors, Eqs. (<ref>) and (<ref>) yield the correct vacuum correlation functions.



§ FINITE EUCLIDEAN SPACETIMES



Above, we have reviewed the reasoning for computing the path integral without specifying boundary conditions in favour of taking complex time to infinity. While yielding the physical correlation functions, taking time to infinity clearly is a mathematical trick. However, as we have discussed in Section <ref>, simply using Eqs. (<ref>) and (<ref>) is not a valid procedure. Nonetheless, it should still be possible, at least in principle, to carry out the calculation in a finite spacetime volume.

To this end, we see three ways of doing this. All of these turn out to require the replacement of the fixed boundary conditions (<ref>) with different configurations that again lead to the same conclusion of CP conservation. These particular possibilities are:


  * We can take a finite time interval at the price of having to project on a vacuum wave functional (see the present section).

  * In order to avoid the projection on the wave functional, we can stay within functional quantization and consider a finite subvolume of infinite Euclidean space. Then, we have to integrate over all possible boundary configurations on the subvolume (see Section <ref>).

  * We can take compact spacetimes without boundaries. For definiteness, consider here a four-torus with a finite Euclidean time interval of length β=1/T, where T is the temperature. The relation with Minkowski-spacetime is given by its correspondence with the canonical thermodynamic partition function. This once again requires canonical quantization that restricts the form of the wave functionals (see the present section).





  
Projecting on the wave functional
Regarding the first option, unless taking complex time to infinity or assuming finite temperature, we need to specify the ground state in order to get the correct boundary conditions for the path integral. That is, when restricting to a real time interval from t^' to t, we must weigh each path contributing to the partition function

    Z_ M(t_2,t_1)
    =   ∫ Dψ̅  Dψ  DA Ψ^*[A_μ(t_2,𝐱),ψ(t_2,𝐱)] Ψ[A_μ(t_1,𝐱),ψ̅(t_1,𝐱)]
    
    ×    e^ i∫_Ω_12^4x ( L_ M-η̅ψ-ψ̅η) ,

by the unknown ground state wave functional

    |Ψ, t⟩=Ψ[A_μ(𝐱,t),ψ̅(𝐱,t))
    ] |A_μ(𝐱,t)⟩⊗|ψ(𝐱,t)⟩

evaluated at the endpoints of these paths. Here, L_ M is the Minkowskian Lagrangian, |A_μ(𝐱,t)⟩, and |ψ(𝐱,t)⟩ are field eigenstates and Ω_12 is the four-volume bounded by the three-volumes at the times t_1 and t_2. 


In turn, for physical boundary conditions that do account for fluctuations about the classical minimal-energy states, topological quantization cannot be assumed. Note while this means that the action generally receives contributions from field configurations of noninteger Δ n, this does not preclude the wave functional from transforming with a phase factor under large gauge transformations. However one may wonder whether one can still assume topological quantization in an approximate sense, so that one might still obtain a good result from the sum over path integrals with boundary conditions corresponding to pure gauges imposed on some finite volume. To settle this question, we would have to find the ground state wave functional in canonical quantization of Yang–Mills theory <cit.>, which does not appear to be practically possible. Finite temperature field theory also relies on canonical quantization, but tangible conclusions may be drawn in that context, as we discuss next.



  
Compact spacetimes without boundary
As for introducing finite spacetime volume through temperature, a notable example is de Sitter space, where the Euclidean counterpart is a sphere, and hence Ω is finite for this spacetime. The latter can be interpreted as the representation of a canonical ensemble of states over a static patch of de Sitter space with a Lorentzian signature. Another important situation where Ω is finite is when space is a three torus and we consider a canonical ensemble as well. Then, the Euclidean representation is a four-torus and corresponds to the continuum limit of computations in lattice QCD. Both of these very relevant examples of finite Euclidean volume therefore require first the canonical quantization in the respective background geometries in order to connect these with physical observables in Lorentzian spacetimes. This is mandatory in order to evaluate the trace of the canonical density matrix from which the path integral representation can then be derived.

In canonical quantization of gauge theory, large gauge transformations on the spatial section can play a special role. These are gauge transformations that are not continuously connected with the identity but give rise to equivalence classes that are a representation of the homotopy group. This also happens for infinite volumes (provided that the gauge field configurations approach a unique value at spatial infinity) and leads to the well-known θ-vacua <cit.>. However, the infinite volume limit taken inside the path integral readily implies CP conservation so that it is interesting to further focus on finite spatial volumes.

On finite spatial volumes, large gauge transformations are only singled out when fixing the gauge corresponding to periodic (i.e. single-valued) gauge potentials on the torus or single-valued gauge potentials on the sphere. Imposing single-valuedness, we find that the canonical quantization on these finite spatial volumes only admits states that are invariant under large gauge transformations, i.e. without any phase incurred, and that there hence can be no CP violation. This also resolves the matter of renormalizability of the states raised in Refs. <cit.>. On the other hand, without imposing single-valuedness, all gauge transformations on the spatial sections are continuously connected with the identity transformation, and again no CP-violating effects can be deduced. The details of this argument shall be published elsewhere <cit.>.




§ DILUTE INSTANTON GAS APPROXIMATION



The most sensitive probe of possible CP violation associated with the strong interactions is the EDM of the neutron. At the relevant energy scale, QCD is deeply in the nonperturbative regime. This is a well-known and obvious drawback for any analytical approximation. Yet, one can observe from semiclassical calculations that instantons play a central role in the spontaneous breaking of chiral symmetry as well as in mediating the effects from the anomalous axial U(1)_ A symmetry that notably explains the large mass of the η^' meson <cit.>. It is therefore also strongly indicated that the role of the topological term can be understood from a semiclassical evaluation of the effective fermion interaction mediated by instantons, i.e. the `t Hooft operator <cit.>. This corresponds to the expectation that the presence or absence of CP violation should prevail when crossing between the strongly and weakly coupled regimes at low and high energies, respectively. Moreover, the generic arguments in Section <ref> as well as in Refs. <cit.>, where cluster decomposition and the index theorem are used, do not refer to the semiclassical approximation.

The semiclassical approximation therefore remains of substantial interest, being the only analytic procedure to make quantitative statements about CP violation in the strong interactions. Also, it offers a very useful perspective on the central issues with this topic.

In the present context, the semiclassical approach is given by the dilute instanton gas approximation. Stationary and quasi-stationary points of the action are described in terms of instantons and their individual collective coordinates <cit.>. Stationary points are the classical solutions.  These are the minima of the action for each topological sector characterized by winding number Δ n.  For Δ n=± 1, they are given by Belavin–Polyakov–Schwarz–Tyupkin (BPST) (anti-)instanton solutions <cit.> whose classical Yang–Mills action is 

    S_ BPST=8π^2/g^2 .

Explicitly, the BPST instanton reads in the regular gauge

    A^a_μ=2η_aμν(x-x_0)_ν/(x-x_0)^2+ρ^2,

where x_0 and ρ are free parameters corresponding to the center location of the instanton and its size, respectively. Here η_aμν are the `t Hooft symbols <cit.>

    η_aμν=ε_aμν,     μ,ν=1,2,3 
     -δ_aν,    μ=4 
    δ_aμ,    ν=4 
     0,    μ=ν=4  .

Similarly one can define η̅_aμν by a change in the sign of δ in the above equation. For the anti-instanton, we should replace η_aμν by η̅_aμν in Eq. (<ref>).

To visualize the BPST instanton (in analogy to Ref. <cit.> for one-dimensional instantons), we consider some explicit expressions with x_0=0, i.e., with the center set at the origin. For example,

    A_1^3(x)=2 x_2/x^2+ρ^2 ,
       A_2^3(x)=-2 x_1/x^2+ρ^2 .

A_1^3(x) is symmetric in the hyperplane {x_1,x_3,x_4} and A_2^3 symmetric in the hyperplane {x_2,x_3,x_4}. Without loss of generality, in Figure <ref>, we show these as a function of x_1,x_2 by taking x_3=x_4=0 in arbitrary units. The field strength components read 

    F_μν^a=-4η_aμνρ^2/x^2+ρ^2 .

As an example, we plot F_21^3 as a function of x_1 and x_2 in Figure <ref>. These quantities are gauge dependent. The gauge-independent quantities are

    tr F_μν F^μν= tr F_μνF^μν =-96ρ^4/[x^2+ρ^2]^4 .

For the anti-instanton, one would have tr F_μν F^μν=- tr F_μνF^μν.
We plot these in Figure <ref>. From the graph, one can see that the instanton indeed has a radius characterized by the value of ρ (ρ=1 in the plot).












For |Δ n|>1, they should be obtained from the Atiyah–Drinfeld–Hitchin–Manin (ADHM) construction <cit.>. In the dilute gas picture, they correspond to Δ n instantons, no anti-instantons for Δ n>0 and -Δ n anti-instantons, no instantons for Δ n<0. The collective coordinates for the individual instantons describe their size, their gauge orientation as well as their position in Euclidean spacetime. As we follow the steepest-descent contours, the action S evolves towards larger values, and we can encounter quasi-stationary points. These can be described in terms of the number n of instantons and n̅ of anti-instantons, where both of these objects can coexist within such configurations. In the sector (i.e. on the thimble) characterized by Δ n, it must hold that n-n̅ =Δ n. Each of these individual instantons and anti-instantons is again parametrized in terms of the aforementioned collective coordinates. In Figure <ref>, we illustrate the typical structure of the thimbles in the semiclassical approximation.




Now, we aim to integrate out the gluon fields in order to see explicitly what quark correlations breaking the anomalous chiral symmetry U(1)_A they leave behind. We follow Ref. <cit.>, but here we work with Euclidean time for simplicity.

 The relevant quark correlation function is given by

    ⟨ψ(x)ψ̅(x^')⟩
    =-.δ^2 log Z[η,η̅]/δη̅(x) δη(x^')|_η=η̅=0=
    ∫ Dψ̅  Dψ  DA ψ(x)ψ̅(x^')
     e^-lim_Ω→∞∫_Ω^4x  L/∫ Dψ̅  Dψ  DA  e^-lim_Ω→∞∫_Ω^4x  L .

While this is a standard expression one should note that 
the numerator and denominator in this equation are not well defined in the thermodynamic limit Ω→∞, even when ultraviolet divergences have been renormalized. However, this does not force us to keep Ω finite. Rather, divergent extensive contributions in the numerator and denominator from spacetime regions far away from x and x^' cancel. In standard perturbation theory, these contributions are represented by vacuum diagrams where the divergence results from their overall invariance under spacetime translations. We go into more detail regarding this point in Section <ref>. In the present semiclassical evaluation, we shall see how to deal with these extensive contributions a bit further down the line of argument. 

To proceed with the evaluation of Eq. (<ref>), we approximate the Green's function of the quarks in the background of one anti-instanton (Δ n=-1 as)

    S(x,x^')
    =ψ̂_λ̅(x)ψ̂^†_λ̅(x^')/λ̅≈ψ̂_0 L(x)ψ̂_0 L^†(x^')/m  e^- iα+_λ≠0ψ̂_λ(x)ψ̂^†_λ(x^')/λ .

The middle expression is the exact spectral sum representation in terms of the eigenvalues λ̅ of the Dirac operator of massive quarks in the anti-instanton background, and ψ̂_λ̅ are the corresponding eigenfunctions. As for the approximation on the right, by ψ̂_0 L,R we denote the `t Hooft zero modes of the massless Dirac operator in the corresponding one anti-instanton or instanton background <cit.>, which are purely chiral and where their handedness is indicated by L,R. The nonzero eigenvalues of the massless Dirac operator are given by λ and the pertaining eigenmodes by ψ_λ. Note that the contribution breaking chiral symmetry, i.e. the first term in the approximate expression, aligns with the CP phase α pertaining to the quark mass and not with the angle θ. For real masses, this approximation has been used in Refs. <cit.>.

In the semiclassical approximation, we carry out the path integral by taking the quasistationary configurations of the action, i.e. with n instantons and n̅ anti-instantons and evaluate the leading fluctuations, i.e. the functional determinants corresponding to one-loop order. For such a quasistationary background, the Green's function for the quarks should be well approximated by <cit.>

    S_n,n̅(x,x^')≈ S_0(x,x^')+   ∑_ν̅=1^n̅ψ̂_0 L(x-x_0,ν̅)ψ̂^†_0 L(x^'-x_0,ν̅)/m  e^- iα
    +   ∑_ν=1^nψ̂_0 R(x-x_0,ν)ψ̂^†_0 R(x^'-x_0,ν)/m  e^ iα .

Here x_0,ν and x_0,ν̅ are the locations of instantons and anti-instantons, respectively, S_0 is the Green's function of a Dirac fermion with mass m exp( iαγ^5) in a translation-invariant (i.e., void of instantons) background. This approximation neglects contributions from overlapping instantons which are more suppressed as the instanton gas becomes more dilute. While the Green's function close to the individual instantons and anti-instantons is dominated and therefore approximated by the `t Hooft zero modes, sufficiently far away from the points x_0,ν and x_0,ν̅ the Green's function is given by the form in the background without instantons, i.e. 

    S_0(x,x^')=(-γ_μ∂_μ+m  e^- iαγ^5)∫ d^4p/(2π)^4  e^- ip(x-x^')1/p^2+m^2 .

In Eq. (<ref>), we note the alignment of the instanton-induced breaking of chiral symmetry with the quark masses so that there is no indication of CP violation at this level but also note that θ has not yet entered into the calculation.


Given the Green's functions (<ref>), we can proceed with evaluating the fermion correlation on the thimble (or equivalently in a fixed topological sector) characterized by the winding number Δ n:  

    ⟨ψ(x) ψ̅(x^')⟩_Δ n
    
    =   ∑_n̅,n≥ 0  n-n̅=Δ n∫ DA_n̅,n Dψ̅ Dψ ψ(x)ψ̅(x^')  e^-S[A,ψ̅,ψ]
    
    =   ∑_n̅,n≥ 0  n-n̅=Δ n1/n̅! n!(∏_ν̅=1^n̅ ∫_Ω d^4x_0,ν̅ dΣ_ν̅ J_ν̅)
    (∏_ν=1^n ∫_Ω d^4x_0,ν dΣ_ν J_ν)
    S_n̅,n(x,x^')
       5cm× 
    
     e^-S_ BPST (n̅ +n) e^- i(n̅ -n)(α+θ)
    (-Θϖ)^(n̅ +n) .

The symbol DA_n̅, n implies that the path integral is evaluated in terms of fluctuations and moduli about the classical background, i.e. (quasi-)stationary point, made up from n instantons and n̅ anti-instantons. The integration over collective coordinates other than the locations of the instantons and anti-instantons are denoted by dΣ_ν,ν̅, and the Jacobians from the transformation of the zero modes in the path integral in favour of the collective coordinates are denoted by J_ν,ν̅. The one-loop determinant of the gauge field about a single instanton or anti-instanton (denoted by A̅ below) with the zero modes omitted and divided by the gauge field determinant in the background A=0 is given by

    ϖ≡1/√(^'_A̅/_A=0) ,

where the prime on the determinant indicates the omission of zero eigenvalues.
In an analogous manner, Θ represents the modulus of the ratio of the fermionic determinants in the one-(anti)instanton and A=0 backgrounds, 

    Θ≡| (- D-m e^ iαγ^5)/ (-∂-m e^ iαγ^5)| .


As usual, the partition function diverges in the thermodynamic limit so that we keep the spacetime volume Ω finite for now. Nonetheless, we need to take Ω→∞ before eventually summing over the topological sectors Δ n as the latter are only a consequence of infinite spacetime volume and to remain true to the integration contour implied by Eq. (<ref>), cf. the discussion in Section <ref>.

In order to normalize, i.e., to divide out vacuum contributions, we also need the partition function in a fixed topological sector. Proceeding as for the fermion correlation, we obtain

    Z_Δ n=   ∑_n̅,n≥ 0  n-n̅=Δ n∫ DA_n̅,n Dψ̅ Dψ  e^-S[A,ψ̅,ψ]
    
    =   ∑_n̅,n≥ 0  n-n̅=Δ n1/n̅! n!(-Ω J  Θ ϖ  e^-S_ BPST∫ dΣ)^(n̅ +n) e^- i(n̅ -n)(α + θ) .


Next, we turn to the collective coordinates and integrate out the location of a single anti-instanton as

    ∫_Ω d^4 x_0,ν̅ S(x,x^')
    ≈    ∫_Ω d^4 x_0,ν̅[S_0(x,x^')
    +ψ̂_0 L(x-x_0,ν̅)ψ̂_0 L^†(x^'-x_0,ν̅)/m  e^- iα+⋯]
    
    
    =   Ω (S_0(x,x^')+⋯) +m^-1 e^ iα h(x,x^')P_ L .

The dots above represent contributions from the zero modes of the (anti)-instantons whose centers were not integrated over. This expression defines the overlap function h(x,x^')—a rank-two tensor in spinor space:

    h(x,x^')P_ L=    ∫_Ω d^4 x_0,ν̅ ψ̂_0 L(x-x_0,ν̅) ψ̂_0  L^†(x^'-x_0,ν̅) ,
    
    h(x,x^')P_ R=    ∫_Ω d^4 x_0,ν ψ̂_0 R(x-x_0,ν) ψ̂_0 R^†(x^'-x_0,ν) .

Further, we integrate over the remaining collective coordinates as

    h̅(x,x^')≡   ∫Σ h(x,x^')/∫Σ .

Notice that we ignore here the fact that for the classical instanton, the integral over the dilatational mode is divergent. The running coupling will however render the correlations finite in a more complete calculation.


Strictly speaking, the dilute instanton gas approximation is only applicable when the integral over the dilatational mode converges, i.e. when contributions from both, large and small instantons are cut off. This is naturally the case in the ultraviolet, where small instantons are suppressed for asymptotically free theories as g→0. In the infrared, one may attempt to keep g perturbatively small by a bespoke particle content that controls the running coupling. As a matter of principle, an infrared cutoff can also be enforced when the gauge symmetry is spontaneously broken so that the size of the instantons is limited by the inverse gauge boson mass <cit.>. While none of this applies to the strong interactions, such considerations show that the dilute instanton gas is a meaningful concept. Note that an infrared cutoff for the instanton size has no implications for the integration over the locations x_0,ν, x_0,ν̅ of the instanton centers. The preservation of Poincaré symmetry, and in particular Lorentz invariance, demands that the values of  x_0,ν, x_0,ν̅ should remain unconstrained. Hence, the dependence of the results on the spacetime volume is unchanged in the presence of a cutoff for the instanton size. As a consequence, the former has no consequence for the order of limits of infinite spacetime volume and infinite maximal absolute value of the topological charge.  Note that even with the aforementioned size cutoff it is clear that either expression (<ref>) or (<ref>) can be technically evaluated. Further, the presence or absence of divergences from infrared instantons does not decide which order of limits must be taken because the presence of sectors of integer Δ n is a topological argument that does not depend on the validity of the semiclassical expansion. We eventually note here that the validity of the dilute instanton gas and its generalization toward the inclusion of interactions between instantons has been addressed in Refs. <cit.>.

The present point of view is that the saddle point approximation in the dilute instanton gas approach, while not quantitatively applicable to the strong interactions, yields information about the symmetries that are respected by the theory. This does not only apply to the present work that argues in favour of the evaluation of the partition function according to Eq. (<ref>) but also to Refs. <cit.> that assume Eqs. (<ref>) and (<ref>). To our knowledge, Ref. <cit.> is the only paper that explicitly evaluates the 't Hooft operator for nonzero θ̅. While the saddle point approximation is an important cross-check, we note that the conclusion about the absence of CP violation does not rely on it, cf. the boxed argument in Section <ref> and in the part of Section <ref> on the evaluation of the partition function and topological quantization as well as Refs. <cit.>. Furthermore, as will be reviewed at the end of Section <ref>, the results obtained with the dilute instanton gas can be recovered from general arguments based on cluster decomposition and the index theorem, without making use of the dilute gas approximation.


Integrating now over all locations of instantons and anti-instantons, we obtain the correlation function for fixed Δ n:

    ⟨ψ(x) ψ̅(x^')⟩_Δ n
    
    =   ∑_n̅,n≥ 0  n-n̅=Δ n1/n̅! n![
     h̅(x,x^')(n̅/m  e^- iα P_ L+n/ m  e^ iαP_ R) Ω^n̅+n -1
    + S_0(x,x^') Ω^n̅+n]
       8cm×κ^n̅ +n
    (-1)^n+n̅ e^ iΔ n(α + θ)
    
    =   [( e^ iα I_Δ n+1(2 κΩ)P_ L+ e^- iα I_Δ n-1(2 κΩ) P_ R)κ/m h̅(x,x^')
    +I_Δ n(2 κΩ) S_0(x,x^')]
       8cm×
    (-1)^Δ n e^ iΔ n(α + θ) ,

where
κ=∫ dΣ J  Θ ϖ  e^-S_ BPST is the instanton density per spacetime volume[Note that in Minkowskian spacetime, we define κ=∫ dΣ  JΘϖ e^-S_ BPST in Ref. <cit.>. The κ in both cases is the same and is real due to the fact that the Jacobian J in Minkowski spacetime contains an additional factor of  compared to its Euclidean counterpart.] and I_ν(x) is the modified Bessel function.

The terms involving the overlap function h̅ are due to the instanton effects on the quarks and break chiral symmetry. While we should expect that these scale in the same way with the spacetime volume Ω as the term with S_0, i.e. the contribution from regions between instantons, the explicit dependence on Ω in Eq. (<ref>) is different. However, we see that the scaling after all is the same. Relax for the moment the constraint of fixed Δ n and use that κ may be interpreted as the likelihood for finding an instanton in a unit four-volume.  Then, for large Ω the sum is dominated by particular value of n≈n̅:

    ⟨ n ⟩=∑_n=0^∞ n(κΩ)^n/n!/∑_n=0^∞(κΩ)^n/n!=κΩ .

Moreover, the relative fluctuation vanishes in the infinite-volume limit <cit.>:

    √(⟨ (n -⟨ n⟩)^2⟩)/⟨ n⟩=1/√(κΩ) .

This means that in the coefficients in front of the chiral projection operators within the middle expression in Eq. (<ref>), we can replace n,n̅→κΩ. This basic behaviour, i.e. that the central value for the number of instantons is given by κΩ, is also reflected by the fact that for large arguments, the modified Bessel functions become independent of their index, i.e. lim_x→∞I_Δ n(x)/I_Δ n^'(x)=1. Since all the modified Bessel functions in Eq. (<ref>) tend to the same value, we see directly from this expression that there is no relative CP phase between the terms from the quark masses and instanton-induced breaking of chiral symmetry in the infinite-volume limit. Correspondingly, the partition function for fixed Δ n turns out as <cit.>

    Z_Δ n=I_Δ n(2 κΩ) (-1)^Δ n e^ iΔ n(α+θ) .


Now, when calculating the correlation function as the sum over the topological sectors, we have to take the limit Ω→∞ first for the reasons explained in Section <ref>. Because of the divergence in the thermodynamic limit, the numerator and denominator have to be treated together, and we obtain

    ⟨ψ(x)ψ̅(x^')⟩
     
     =   lim_N→∞ N∈ Nlim_Ω→∞∑_Δ n=-N^N⟨ψ(x) ψ̅(x^')⟩_Δ n/∑_Δ n=-N^N Z_Δ n
    =   S_0inst(x,x^')
    +κh̅(x,x^') m^-1 e^- iαγ^5 .

In Section <ref>, we show that this procedure amounts to dropping the divergent extensive contributions that correspond to the vacuum diagrams in standard perturbation theory. In this final result, the phase from the quark mass in S_0 inst, cf. Eq. (<ref>), is aligned with the phase from the instanton-induced effects in the term with the overlap function h̅, so that there are no CP-violating effects.

One may wonder about Eq. (<ref>) why we take the limit N→∞ in front of the fraction whereas by Eq. (<ref>), it appears that it should hold for numerator and denominator separately in the first place. As we have noted though, without normalization by vacuum contributions, the partition function is not well defined in the thermodynamic limit. The present procedure is necessary to divide out the extensive contributions causing the divergence. It is unique in the sense that we carry out the integrals over each steepest-descent contour before interfering them. Doing otherwise would correspond to a partitioning and reordering of the full integration contour that consists of the steepest-descent contours connected via configurations of infinite action, see Figure <ref> for illustration. This amounts to an incorrect manipulation of a path integral that is not absolutely convergent.


Now we consider what happens when the limits are ordered the other way around, i.e. sum over the topological sectors before taking Ω→∞, according to Eqs. (<ref>) and (<ref>). We reiterate though that this procedure is not valid because topological quantization can only be deduced in infinite spacetime volume.
As for the fermion correlation, one obtains

    ∑_n̅,n≥ 01/n̅! n![
     h̅(x,x^')(n̅  m^-1 e^ iα P_ L+n  m^-1 e^- iα P_ R) (Ω)^n̅+n -1+ S_0(x,x^') (Ω)^n̅+n]
       8cm×
    (- κ)^n̅ +n e^ iΔ n(α + θ)
    
    =   [-( e^-iθ P_ L+ e^ iθ  P_ R)κ/mh̅(x,x^')
    + S_0(x,x^')] e^-2  κΩcos(α+θ) ,

and for the partition function

    ∑_n,n̅1/n!n̅!(-κΩ)^n̅+n e^- i  (n̅-n)(α+θ)= e^-2κΩcos(α+θ) .

Taking the ratio, the overall exponential factors cancel but now there is a misalignment between the phases in S_0 inst and in the instanton-induced term. This means that as Ω→∞, there is an infinite amount of destructive interference that suppresses the statistically more likely contributions with approximately equal numbers of n and n̅ (see Eq. (<ref>)) in favour of outliers for which Δ n/Ω does not go to zero. Equations (<ref>) and (<ref>), if they were correct, would signal CP-violating effects. Note that in either result, terms that break the U(1)_A symmetry from both, instanton-mediated effects and the quark mass m are present. When we turn to the phenomenology of the strong interactions and generalization to several flavours in Section <ref>, we shall recall that both, the breaking of chiral symmetry through the quark masses and from instantons are necessary in order to explain the spectrum of mesons, in particular why the η^'-meson is much heavier than the pions <cit.>. In either order of limits, this phenomenology is explained. So the meson spectrum alone cannot be used in order to conclude the the correct order of limits.



§ THERMODYNAMIC LIMIT AND CLUSTER DECOMPOSITION



We establish here that with the limiting procedure in Eq. (<ref>), contributions that are divergent due to the infinite spacetime volume cancel between numerator and denominator. This corresponds to the usual cancellation of vacuum diagrams when evaluating connected correlation functions in standard perturbation theory (i.e. without expanding around nontrivial classical solutions).

The present argument is also interesting for what concerns Eq. (<ref>). The partition function is defined in the limit Ω→∞ in the first place. This appears as an obstacle to using log Z as an extensive, volume-dependent quantity in line with what is familiar from thermodynamics. We shall see here that an expression with such a property can nonetheless be defined when restricting Z to some subvolume of Ω. Note that as such a restriction is arbitrary, no boundary conditions on the subvolume can be placed.

While we are working here at zero temperature, we note that one can use the Polyakov line at finite temperature in order to control and study the deconfinement phase transition, including contributions from the gradient expansion of the quark determinant <cit.>.

We use a well known line of reasoning <cit.> and consider the expectation value of an operator O in an infinite spacetime volume Ω, and interfere different topological sectors Δ n as

    ⟨ O⟩= lim_N→∞ N∈ Nlim_Ω→∞∑_Δ n=-N^N f(Δ n)∫_Δ n Dϕ  O  e^-S_Ω[ϕ]/∑_Δ n=-N^N f(Δ n)∫_Δ n Dϕ  e^-S_Ω[ϕ] ,

where Dϕ is the path integral measure over all fields involved. Now, let O be an operator corresponding to a correlation function evaluated for some spacetime points. For example, in Eq. (<ref>), O=ψ(x)ψ̅(x^'). For the action, we write S_Ω to indicate that it is obtained from integrating the Lagrangian over the spacetime volume Ω. As for the Lagrangian, we take it not to include the topological term -θ tr F F/(16π^2). Rather, we have the function f(Δ n) taking care of the dependence on the topological sector.

Now, consider partitioning the spacetime volume as Ω=Ω_1∪Ω_2 so that Δ n(Ω)=Δ n_1(Ω_1)+Δ n_2(Ω_2). We further assume that the spacetime arguments of the operator fall within Ω_1, and we write O_1 in favour of O to indicate this. We can thus write

    ⟨ O_1⟩= lim_N_1, N_2→∞ N_1, N_2∈ N
     lim_Ω→∞∑_Δ n_1=-N_1^N_1∑_Δ n_2=-N_2^N_2 f(Δ n_1+Δ n_2)∫_Δ n_1 Dϕ  O_1  e^-S_Ω_1[ϕ]∫_Δ n_2 Dϕ  e^-S_Ω_2[ϕ]/∑_Δ n_1=-N_1^N_1∑_Δ n_2=-N_2^N_2 f(Δ n_1+Δ n_2)∫_Δ n_1 Dϕ  e^-S_Ω_1[ϕ]∫_Δ n_2 Dϕ  e^-S_Ω_2[ϕ] .

Since there may be instantons sitting right at the boundaries of the two subvolumes, Δ n_1,2 will not be strictly integer. However, if the instanton gas is sufficiently dilute, integer winding numbers may still correspond to an adequate approximation.

Now, as required by the cluster decomposition principle, provided Ω_1 is chosen large enough, ⟨ O_1⟩ must not depend on contributions from Ω_2 to the path integral. This is generally the case when the numerator and the denominator decompose into factors that only depend on Δ n_1, Ω_1 or  Δ n_2, Ω_2, respectively. Then the contributions from the volume Ω_2 can be reduced from the fractions. This generally happens when

    f(Δ n_1+Δ n_2)=f(Δ n_1)f(Δ n_2)⇒ f(Δ n)= e^ iΔ n θ .

So the contributions from the topological term, that we have left aside thus far, can indeed be accounted for through f(Δ n). Note that the argument holds for either order of limits, i.e. the one from Eqs. (<ref>), that implies Eq. (<ref>) which is imposed here, as well as for the commuted version from Eq. (<ref>).

To carry out the limits, we write Eq. (<ref>) as

    ⟨ O_1⟩=   lim_N_1, N→∞ N_1,N∈ N
     lim_Ω→∞ ∑_Δ n=-N^N ∑_Δ n_1=-N_1^N_1 f(Δ n)∫_Δ n_1 Dϕ  O_1  e^-S_Ω_1[ϕ]∫_Δ n_2=Δ n-Δ n_1 Dϕ  e^-S_Ω_2[ϕ]/∑_Δ n=-N^N ∑_Δ n_1=-N_1^N_1 f(Δ n)∫_Δ n_1 Dϕ  e^-S_Ω_1[ϕ]∫_Δ n_2=Δ n -Δ n_1 Dϕ  e^-S_Ω_2[ϕ] .

Corresponding to Eq. (<ref>), the integrations over the volume Ω_2 lead to

    ⟨ O_1⟩=   lim_N_1, N→∞ N_1,N∈ N
     lim_Ω_2→∞
       -.4cm∑_Δ n=-N^N∑_Δ n_1=-N_1^N_1 f(Δ n) I_Δ n-Δ n_1(2κΩ_2)(-1)^N_f(Δ n-Δ n_1) e^ i α(Δ n-Δ n_1)∫_Δ n_1 Dϕ  O_1  e^-S_Ω_1[ϕ]/∑_Δ n=-N^N∑_Δ n_1=-N_1^N_1 f(Δ n) I_Δ n-Δ n_1(2κΩ_2)(-1)^N_f(Δ n-Δ n_1) e^ i α(Δ n-Δ n_1)∫_Δ n_1 Dϕ  e^-S_Ω_1[ϕ] .

The explicit exponential factors here are phases from the fermion determinants.
Note that these have not been absorbed in κ, which we have defined to be real.

Now we are aiming for an expression for ⟨ O_1 ⟩ with finite Ω_1, without making reference to Ω_2. This proves to be possible because the contributions from Ω_2 can be interpreted as vacuum factors that reduce out from the normalized expectation value.

Since we must take Ω→∞ to have well-defined integer Δ n, we also have to take here Ω_2→∞. The Bessel functions with a factor Ω_2 in their argument then go to a common limit so that we can factorize out the sum over Δ n. We are left with

    ⟨ O_1⟩=   ∑_Δ n_1=-∞^∞∫_Δ n_1 Dϕ (-1)^-N_fΔ n_1 e^- i αΔ n_1 O_1  e^-S_Ω_1[ϕ]/∑_Δ n_1=-∞^∞∫_Δ n_1 Dϕ (-1)^-N_fΔ n_1 e^- i αΔ n_1 e^-S_Ω_1[ϕ] .

We therefore see that taking the limits as in Eq. (<ref>) leads to the correct cancellation of “disconnected” terms, in particular those that originate from regions that are far separated from the spacetime arguments of the observable O_1.

Moreover, in Eq. (<ref>) the θ-angle from the function f(Δ n) does not occur anymore. We can see this as a consequence of phases incurred in Ω_1 being canceled against complementary phases from Ω∖Ω_1. The remaining explicit dependence on the unphysical phase α cancels when the fermionic part of the path integral Dϕ is carried out. Since the path integral here is restricted to Ω_1, which is finite, we can compute the expectation values in finite volumes after all from a partition function in the form of Eq. (<ref>) but with the parameter θ̅ set to zero. This way, the logarithm of the partition function can be taken as an extensive quantity.

In the previous derivation, when going from Eq. (<ref>)  to (<ref>) we made use of the result for the partition function in the dilute instanton gas approximation, Eq. (<ref>). However, it is worth pointing out that the latter result can be derived from the cluster decomposition principle alone, without making use of the dilute instanton gas approximation. One can start by noting  that the factorization of the path integrals in  the denominator in Eq. (<ref>) can be written in terms of the following relations between the partition functions Z_Δ n(Ω) in the full volume and their counterparts Z_Δ n_1(Ω_1), Z_Δ n_2(Ω_2) for the subvolumes Ω_1, Ω_2,

    Z_Δ n(Ω) = ∑_Δ n_1 Z_Δ n_1(Ω_1)Z_Δ n-Δ n_1(Ω_2) .

Equation (<ref>) is an infinite set of identities that can be used to solve for Z_Δ n(Ω) from a set of minimal assumptions. First, we note that Z_Δ n(Ω) are complex. For starters, they receive a phase e^ iθΔ n due to the-θ term.
Further complex phases in Z_Δ n(Ω) can only come from the phases α_i of the fermion masses. At least at the leading order, the fermionic path integration yields determinants of the massive Dirac operator in a background of topological charge Δ n, which can be fully general and is not assumed to be precisely captured by the dilute instanton gas approximation. The phase of the total fermionic determinant is then fixed by the Atiyah-Singer index theorem <cit.>, and for a single fermion is given by e^αΔ n. As a consequence, one can write

    Z_Δ n(Ω)=  e^ (θ+α)Δ ng̃_Δ n(Ω) ,   g̃_Δ n(Ω)∈ℝ .

Parity considerations and appropriate limits of the cluster-decomposition relation (<ref>) can be used to motivate the simple ansatz <cit.>

    g̃_Δ n(Ω) = Ω^|Δ n|f_Δ n(Ω^2),   f_Δ n(0)≠ 0 .

Notably, the previous ansatz together with the assumption of analyticity in Ω give rise to a unique solution for the infinite tower of identities in Eq. (<ref>), which can be written as

    Z_Δ n(Ω) =  e^(θ+α)I_Δ n(2βΩ) ,  β≡ f_Δ n=1(0)∈ℝ .

As advertised, this recovers the result of Eq. (<ref>) without making use of the dilute instanton gas approximation.

Equation (<ref>) can be taken even further, as it allows one to rederive the phases of fermionic correlators and confirm the conclusions of this section without using the dilute instanton gas approximation. Defining a complex mass parameter 𝔪 as

    𝔪≡ m  e^α ,

the mass terms in the Lagrangian of Eq. (<ref>) can be written as

    L⊃𝔪ψ̅P_R ψ+𝔪^*ψ̅P_Lψ .
 
Then, one can view the complex mass parameters as sources for integrated correlators,

    ∂/∂𝔪Z_Δ n= -∫ d^4 x ⟨ψ̅P_ Rψ⟩_Δ n ,   ∂/∂𝔪^*Z_Δ n= - ∫ d^4 x ⟨ψ̅P_ Lψ⟩_Δ n .

As the partition functions of Eq. (<ref>) have been derived on general grounds, the previous correlators are meant to include nonperturbative effects.
Noting that the reality condition in the β parameter of Eq. (<ref>) implies β=β(𝔪𝔪^*) and writing α =-(/2)log(𝔪/𝔪^*) yields the following spacetime averaged correlators <cit.>

    1/Ω∫^4 x ⟨ψ̅P_ Rψ⟩ =     lim_N→∞lim_Ω→∞∑_|Δ n|<N∫^4x ⟨ψ̅P_ Rψ⟩_Δ n/Ω∑_|Δ m|<N Z_Δ m = -2𝔪^* ∂_𝔪𝔪^*β(𝔪𝔪^*) ,
    
      1/Ω∫^4 x ⟨ψ̅P_Lψ⟩ =    lim_N→∞lim_Ω→∞∑_|Δ n|<N∫^4x ⟨ψ̅P_ Lψ⟩_Δ n/Ω∑_|Δ m|<N Z_Δ m =-2𝔪 ∂_𝔪𝔪^*β(𝔪𝔪^*) .

It is readily seen that the total phase of the fermionic correlators, including nonperturbative effects, is aligned with the phases of the tree-level masses in the Lagrangian. This generalizes the result of Eq. (<ref>) and leads again to the conclusion of no CP violation. Again, the order of limits plays a crucial role in Eq. (<ref>).



 


§ EFFECTIVE THEORIES AND EFFECTIVE OPERATORS





We shall now draw the connection from the results of the semiclassical approximation corresponding to integrating out gluons from Section <ref>  with observables probing CP conservation or violation in the strong interactions. The main object of interest in that context is the `t Hooft vertex, which can be inferred from the correlation function (<ref>) as the Lagrangian term 

    -ψ̅(x) Γ e^ iαγ^5ψ(x) .

This vertex generates the same correlation functions as in Eq. (<ref>) for the EFT where gluons have been integrated out. Figure <ref> illustrates how such a model fits into the picture of the different EFTs discussed in the present context. In addition, there will also be in general nonlocal operators from the long-range interactions of the gluons, because with quark degrees of freedom still in the theory, there is no cutoff parameter that allows for a local expansion. The new operators appear in favour of the gluon kinetic term F_μνF_μν as well as the topological term F_μνF_μν which disappear together with the gluons.



In the dilute instanton gas approximation in Section <ref>, one has integrated out gluons in the semiclassical approximation. For this to be valid, the theory should be perturbative throughout, which can be achieved in principle by adding a bespoke matter content that controls the renormalization group evolution in this peculiar way. Certainly, however, this is neither the case for the theory specified in Eq. (<ref>) with the gauge group SU(2) and one quark flavour nor for QCD with SU(3) and three flavours of light quarks. As a consequence, one should expect substantial deviations from the correlation function given in Eq. (<ref>), in particular for large distances between x and x^'. Nonetheless, there should still be a small distance, high energy contribution of this form. When drawing conclusions about CP conservation or violation, one therefore must make the assumption that the CP-odd coefficient of the `t Hooft vertex appears in the same way within the extra operators that have to be added in principle to account for the low-energy behaviour. Note however that this shortcoming applies to the conclusions based on either order of limits when the calculation is carried out semiclassically.

To some extent, the above matter is addressed by the concluding argument from Section <ref>, where the leading fermion correlations are constrained without the dilute instanton gas approximation but using instead cluster decomposition and the index theorem. There, no assumption about the fermion correlation is made but for its U(1)_A-violating form. While the resulting fermion correlation then can only be stated in the coincident limit, the conclusions about CP conservation based on the order of the infinite volume limit and the sum over topological sectors should therefore extend to the nonperturbative low-energy regime as well.

The underlying theory that we are concerned with after all is QCD, which is specified (now with the gauge group SU(3),  N_f flavours of light quarks and in Minkowski spacetime) as (we choose ε_1230=+1 in Minkowski spacetime)

    L_ M⊃-1/2g^2 trF_μνF^μν+∑_i,j=1^N_fψ̅_i(γ^μ D_μδ_ij-M_ij P_ R-M^†_ij P_ L)ψ_j+
    1/16π^2θ  tr F_μνF^μν ,

where in the mass-diagonal basis

    M_ij=δ_ij m_j   (no sum over j) .

For the model as in Eq. (<ref>), the vertex corresponding to Eq. (<ref>) is

    -Γ_N_f e^- iα̅∏_j=1^N_f(ψ̅_j P_ Lψ_j)-Γ_N_f e^ iα̅∏_j=1^N_f(ψ̅_j P_ Rψ_j) ,

where

    α̅=∑_i^N_fα_i= M .



We need to sort in what way (cf. Figure <ref>) this is connected to the EFT of hadrons that is valid at low energies and should describe those possible CP-violating effects that are accessible by current precision experiments. A principal obstacle to systematically deriving quantitative predictions lies of course within the circumstance that perturbation theory is not valid anymore at low energies. 

Yet, the symmetries, even when realized approximately only, offer a standard method of constraining the EFT. In the fundamental theory as well as on the EFT side one can introduce operators of the physical fields coupled to external sources (sometimes called spurions) so that these operators are invariant under local symmetry transformations. On the side of the EFT, the coefficient of these operators has to be obtained through computational or experimental matching. Variation with respect to these sources then allows one to express matrix elements of the fundamental theory in terms of parameters of the EFT.

In the present case, we can apply this method by perceiving the quark masses that break the chiral flavour symmetries SU(N_f)_A as well as the operators breaking U (1)_A as external sources that transform according to these explicitly broken symmetries. In the following discussion, we occasionally let N_f=2, meaning that only up and down quarks are considered, for simplicity. But for expressions explicitly depending on N_f, we keep N_f general. First, we parametrize a chiral transformation as

    ψ→ L P_ Lψ + R P_ Rψ ,

where P_ L,R=(1∓γ^5)/2 and L,R are independent unitary matrices. For an axial transformation, R=L^-1 so that the SU(2)_A transformations are given by

    ψ_i→ [ e^ iγ⃗·σ⃗γ^5]_ijψ_j .

The Lagrangian (<ref>) would remain invariant if the mass matrix transformed as

    M→ L M R^†=  e^- iγ⃗·σ⃗ M e^- iγ⃗·σ⃗ .

In this transformation, M corresponds to a spurion field.

The corresponding EFT Lagrangian (cf. Figure <ref>) with the lowest-order terms is (see, e.g., Refs. <cit.> and Ref. <cit.> where the effective theory is derived from integrating out quark fields)

    L^ EFT_ M = f_π^2/4 Tr ∂_μ U∂^μ U^†+f_π^2 B_0/2  Tr (MU+U^† M^†)+|λ|  e^- iξ f_π^4  det  U+|λ|  e^ iξ f_π^4  det  U^† ,

where 



     

     




    U= U_0  e^ i/f_πΦ=U_0 Ũ ,   U_0=⟨ U ⟩=([ e^ iφ_u       0;       0 e^ iφ_d ]) ,  Φ=[[  π^0+η^' √(2) π^+; √(2) π^- -π^0+η^' ]] .

In the equations above, f_π is the pion decay constant and λ, B_0 are EFT coefficients to be determined experimentally or computationally and B_0 is directly related to the magnitude of the chiral quark condensate. The phases of the latter correspond to  φ_u,φ_d, and we have assumed a diagonal mass matrix M. The squared pion and η' masses are then given by

    m^2_π=    B_0 (m_ucos(α_u+φ_u)+m_dcos(α_d+φ_d)) ,
    
        m^2_η'≈    16|λ| f^2_πcos(ξ-φ_u-φ_d) ,


where we have made the phenomenologically valid approximation that m_π^2≪ m^2_η^'.

We leave the terms with the parameter λ aside just yet as U is invariant under SU(2)_A but transforms with U(1)_A in a way that we shall get to shortly. In correspondence with Eq. (<ref>), the meson fields behave under axial transformations as


    U→ R U L^† =  e^ iγ⃗·σ⃗ U  e^ iγ⃗·σ⃗

The term with the parameter B_0 should be matched so that the correct correlation functions are produced. Corresponding to the invariance of the underlying theory (<ref>), the Lagrangian (<ref>) is invariant under the simultaneous SU(2)_A transformations (<ref>) and (<ref>).

Now, after all, the (up and down) quark masses do not transform under SU(2)_A, they rather break this symmetry explicitly. We can still perceive these as local sources though, that perturb the correlators of the theory about the case with full SU(2)_A symmetry. For the local source, we can then take the fixed physical values of M so that Eq. (<ref>) accounts for the perturbation through the quark masses to linear order. In the EFT, one can continue this to higher orders pending on the precision that is aimed for.

Now, consider U(1)_A transformations 


    ψ_i→ e^ iβγ^5ψ_i

and recall the expression for α̅ from Eq. (<ref>).
The fundamental theory (<ref>) would remain invariant if the quark mass transformed as

    M→ e^-2 iβ M  ,  so that  α̅→α̅-2β N_f .

The chiral anomaly requires that the coefficient θ of the topological term goes as

    θ→θ+2 N_f β

in order to keep the Lagrangian invariant. Note that this implies that the combination

    θ̅=θ+α̅

is invariant under chiral rephasings and in general is nonzero. The presence of such an invariant does however not yet guarantee that it leads to physical effects.

We thus see that there are two local sources that transform under the symmetry U(1)_A: θ and -α. Noting that under this symmetry

    U→ e^2 iN_fβ U ,

the EFT Lagrangian (<ref>) remains invariant if either

    ξ=
       -α̅
      θ .

In principle, one may also allow linear combinations of the parameters -α̅ and θ. As this does not follow from either order of limits for the sum over topological sectors and spacetime volume that we discuss here, we do not consider this combination option further.

We also note that the operator with the coefficient B_0 breaks U(1)_A. So instead of the quark mass phase in M, one could also use θ to write this as an invariant operator with the help of chiral-variant source fields. However, the symmetric theory should respond to  U(1)_A-breaking perturbations through a quark mass term in the same way as it does for SU(2)_A-breaking. In this sense, the term with B_0 is unique to linear order in M. The explicit breaking of U(1)_A through instantons is independent of the quark masses, cf. Eq. (<ref>) together with the fact that κ= O(|M|), and therefore M does not appear in the terms with U.


Now recall that Eq. (<ref>) leads to the effective vertex (<ref>) in the theory where gluons have been integrated out. At this level, θ has disappeared so that the only option for the EFT Lagrangian (<ref>) is

    ξ=-α̅ .

The CP-odd coefficients can then be removed by an overall field redefinition. On the other hand, if it were ξ=θ, there would be a residual CP-odd term.



Further, note that Eq. (<ref>) shows that the mass of the η^' in general does not vanish in the limit of m_u,d→ 0, no matter which of the values ξ takes in Eq. (<ref>). In turn, the fact that the η^' is heavy compared to the pions as such does not lead to a conclusion about which is the correct order of limits.





Finally, the parameter ξ in the coefficient of the `t Hooft operator enters the calculation of the nucleon EDM as follows: Given the EFT Lagrangian (<ref>) and choosing a basis in which M is diagonal, the minimum of the field U is given by U_0 as in Eq. (<ref>) 




where in the limit of |λ|≫ B_0 m_d/f_π^4, and for ξ+α_u+α_d in the first quadrant, one has <cit.>

    m_u sin(φ_u+α_u)=m_d sin(φ_d+α_d)=sin(ξ+α_u+α_d)/√(1/m_u^2+1/m_d^2+2cos(ξ+α_u+α_d)/m_u m_d) .

Going beyond the assumption |λ|≫ B_0 m_d/f_π^4 leads to a mixing of the flavour eigenstates π^0 and η^' within the mass eigenstates.

In order to expand in terms of the meson fields, following  Eqs. (<ref>) and (<ref>) leads to

    Tr[M U +U^† M^†]=   1/2[m_u cos(φ_u+α_u)+m_d cos(φ_d+α_d)]
     Tr[Ũ+Ũ^†]
    +    i/2[m_u sin(φ_u+α_u)+m_d sin(φ_d+α_d)] Tr[Ũ-Ũ^†]
    +   terms mixing π^0 with η^' .

The operators in the second line are CP odd, and we note that

    Tr[Ũ-Ũ^†]=2 i/f_π  Tr[Φ]- i/3 f_π^3 Tr[Φ^3]=4 i/f_π η'- i/f_π^3(2/3η^'^3+2η^'[(π^0)^2+2π^+π^-]) .




Substituting this into the term with B_0 in Eq. (<ref>) and expanding in the meson fields, one generally would obtain CP-violating effects if ξ≠-α̅, the most immediate consequence of which would be η^'→ 2π (recall that the meson fields are CP-odd) through the interaction term

    L^ EFT_ M⊃   
    B_0sin(ξ+α_u+α_d)/f_π√(1/m_u^2+1/m_d^2+2cos(ξ+α_u+α_d)/m_u m_d)[(π^0)^2+2π^+π^-]η^'
    ≈   m_u m_d (ξ+α_u+α_d)/(m_u+m_d)^2m_π^2/ 6 f_π Tr[Φ^3]
     ,

The latter expression (which follows when assuming ξ+α_u+α_d≪ 1) is shown here for comparison with Eq. (8) of Ref. <cit.>. In the latter, the quark masses are taken as real, (see Eq. (5) in   <cit.>) which means that  the parameter θ of that reference should correspond to  ±θ̅ in Eq. (<ref>). Matching the resulting signs for the phases of the quark masses in Eq. (8) of Ref. <cit.> leads to an identification with -θ̅. It then follows that up to the central issue that which value in Eq. (<ref>) is taken by ξ, the results from the present EFT description and the partially conserved axial currents in Ref. <cit.> are therefore in agreement, as they should be. We further compare the two approaches given different values of ξ in Section <ref>. Finally, let us note also that the coefficient in Eq. (<ref>) is different in the approximation of three light flavours where an extra factor of √(2/3) occurs.

To see what the above CP-odd interactions of the pions and η' would imply for the nucleons, one can add their interactions to the EFT Lagrangian as

    L^π, p, n_ M⊃ iN̅∂ N
    -(m_NN̅ U P_ LN
    + icN̅ U^†∂ U P_ L N
    +d N̅ M^† P_ L N
    +e N̅ U M U P_L N
    + h.c.) ,

where the nucleon doublet transforms as

    N→ L  P_ L N+R  P_ R N .

Again, promoting M to a source that transforms under the axial symmetries rather than breaking these, this Lagrangian is invariant. Substituting the expectation value of the chiral condensate (<ref>), (<ref>) for small ξ+α̅, expanding in the meson field and applying field redefinitions N→ N so as to obtain the canonically normalized flavour eigenstates of the nucleons, one finds the interaction terms <cit.>

    L^neutron_ M⊃
    c_1∂_μπ^aN̅τ^a/2γ^μγ^5  N+c_2(ξ+α̅)N̅π^a τ^a/2 N .

The first of these is CP even, as it couples two axial currents, π^a being a pseudoscalar field. The second term is CP odd, as it couples a scalar density with a pseudoscalar field. At one loop level, if it were ξ≠-α̅, this would induce an EDM through the famous diagrams shown in Figure <ref>. Note that the weak interactions make an additional contribution to the neutron EDM <cit.>, which however is too small and usually neglected in the discussion of CP in the strong interactions.









§ SOME OBJECTIONS AND ANSWERS TO THESE


We review here and reply to some objections that we have been made aware of, mostly to the extent that these can be related to (partly earlier) articles or to conference talks that have been published online.



  
Instanton configurations in the different limits
As argued in Section <ref> the different orders of limits correspond to a partitioning and rearrangement of the integration contour that leads to inequivalent results for the path integration. Still, upon taking the limits, the path integral covers the same field configurations. Yet, arguments have been put forward that partitioning the contour should be harmless or that the order of limits as in Eqs. (<ref>) and (<ref>) does not include all relevant contributions as opposed to Eqs. (<ref>) and (<ref>) <cit.>.

Regarding the partitioning of the contour, one may state that for a given configuration of finite 
action, one can find a radius R so that <cit.>

    A_μ(x)= iω(x) ∂_μω^-1(x)+ O(1/R^2)

and, consequently,

    Δ n=1/16π^2∫_|x|<R d^4x  trFF=integer+ O(1/R) .

So in this sense, even for finite volumes one can at least approximately categorize certain configurations by an integer Δ n. However, R is not universal and not even a function of Δ n since the individual units of winding number can be separated arbitrarily far for given Δ n. So we cannot take this as an argument for the calculation in fixed volumes Ω to be equivalent to the full result.

Concerning Eq. (<ref>) perhaps not accounting for all relevant configurations, one may attempt to argue as follows: Consider the dilute instanton gas picture and let ρ be the radius of an instanton. Then, demanding that the instantons and anti-instantons do not overlap, the maximum number of instantons and anti-instantons n+n̅ satisfies |Δ n| ≤ N ≤ (n+n̅)_ max≈ Ω/ρ^4. Therefore, there is room to take N→∞ as Ω→∞ and it seems not to be appropriate to cap Δ n by some finite value in an infinite volume as Eq. (<ref>) suggests.

However, such a cap is not forced by Eq. (<ref>) in the following sense: Decompose the full spacetime into subvolumes, e.g., Ω=Ω_1∩Ω_2. For each topological sector Δ n, there is a constraint Δ n_1+Δ n_2=Δ n where Δ n_1,Δ n_2 are the winding numbers (which may not be integers precisely) in the subvolumes, but no constraint on Δ n_1 and Δ n_2 separately. Therefore, Δ n_1 can be arbitrarily large. For an observer in a subvolume, say Ω_1, the path integral (<ref>) therefore includes configurations of arbitrarily large winding number density Δ n_1/Ω_1 within Ω_1. Hence there is no cap on Δ n_1 in the finite subvolume Ω_1. Indeed, this is already clear from Eq. (<ref>) which is derived from the cluster decomposition principle.



  
Chiral limit

In a theory with at least one massless quark, the parameter α̅ can be chosen arbitrarily with no consequence to the Lagrangian. Without further ado, this implies that θ̅ in Eq. (<ref>) cannot be physical and that there is no CP violation in such a model. In the dilute instanton gas picture, this behaviour results from the suppression of single instantons through the zero-mode from the fermion determinant that makes the factor ϖ in Section <ref> and consequently κ vanishes proportionally to the absolute value of the quark mass determinant.

In this sense, as κ∼ m, we can take in Eq. (<ref>) m→ 0 and obtain a well-defined limit. Since the gluons and therefore the topological term have been integrated out, the parameter α in this expression is still arbitrary but unphysical, as it can be removed by a chiral rotation of the fermion field. The same reasoning applies to the effective operator (<ref>).



  
Effective theory for the topological current
We respond here to comments concerning Ref. <cit.> that have been made in Ref. <cit.>, see also Ref. <cit.>. It is argued there that finite spacetime volumes Ω lead to an unphysical breaking of the conservation of the winding number Δ n (or the density Δ n/Ω) so that the absence of CP violation would be an artifact of such regularization. But apparently, in Eq. (<ref>) and consequently Eq. (<ref>) that lead us here as well as in Ref. <cit.> to the conclusion of no CP violation, Ω→∞ is the first limit that is taken. This is in contrast to Eqs. (<ref>) and (<ref>), which lead to CP-violating observables but where the topological sectors are fixed prior to taking Ω→∞. So the criticism of producing finite volume artifacts would rather be an issue for the latter prescription, and in fact it is, as we have discussed in the previous sections.

While this comes as a rather immediate conclusion, it is of interest to see how the EFT for the topological current, which is introduced in Refs. <cit.>, fits into the present considerations of imposing boundary conditions and orders of limits. (See Figure <ref> for where it stands in relation to the other EFTs that are discussed here.) The topological current can be defined as 

    K_μ=1/4π^2ϵ_μναβ tr[1/2 A_ν∂_α A_β-  i/3A_ν A_α A_β] ,

where we recall tr(T^a T^b T^c)= f^abc/4 for SU(2). The topological charge density and hence the topological term can be explicitly written as a total divergence

    q=1/16π^2 tr F_μνF_μν=∂_μ K_μ .


One of the interesting points concerning the current K_μ is the form of its two-point correlation. Some information can be extracted from the chiral susceptibility

    χ_Ω=   1/Ω .⟨Δ n^2⟩|_θ̅=1/2(1-(-1)^N_f)π
    =   1/Ω .⟨(∫_Ω^4x  q(x))^2⟩|_θ̅=1/2(1-(-1)^N_f)π=
    ∫_Ω^4x .⟨ q(x) q(0)⟩|_θ̅=1/2(1-(-1)^N_f)π .

We evaluate this for the CP-even values of θ̅= 1/2(1-(-1)^N_f)π. When summing over topological sectors in a finite spacetime volume with fixed boundary conditions according to Eqs. (<ref>) and (<ref>), the vacuum energy is minimized at the chosen value of θ̅ and χ remains positive. Choosing instead the starting point Eq. (<ref>) and consequently Eq. (<ref>), the value of θ̅ is irrelevant by the arguments of Sections <ref> and <ref>. We attach a subscript on χ in order to indicate that it is important which volume is referred to. Significant differences can arise when Ω corresponds to a subvolume of the spacetime as opposed to the full spacetime, in particular when boundary conditions control the overall topological fluctuations. For expressions that apply to all volumes we omit the subscript.

One should also pay attention to the fact that χ will in general have connected and disconnected parts. If there is CP violation in a certain setup, e.g. imposed by unphysical boundary conditions, then ⟨Δ n⟩/Ω≠0. When one aims to characterize  the volume-scaled variance of the topological charge, one should then subtract the disconnected contributions, i.e. consider

    χ_Ω-⟨Δ n ⟩^2/Ω .

Of course, one may also define χ without the disconnected pieces to start with, which we do not do here for the sake of simpler expressions.

Applying the result (<ref>) to Eq. (<ref>) yields

    χ_Ω_1≡   1/Ω_1⟨(∫_Ω_1^4x  q(x))^2⟩
    
    =   1/Ω_1∑_Δ n_1=-∞^∞∫ DA_Δ n_1 Dψ̅ Dψ∫_Ω_1 d^4 z  q(z) ∫_Ω_1 d^4 z^'  q(z^') e^-S_Ω_1[A_μ](-1)^-N_f Δ n_1 e^- i α̅Δ n_1/∑_Δ n_1=-∞^∞∫ DA_Δ n_1 Dψ̅ Dψ e^-S_Ω_1[A_μ](-1)^-N_f Δ n_1 e^- i α̅Δ n_1
    
    =   1/Ω_1∑_n_1=0^∞∑_n̅_1=0^∞1/n_1!n̅_1!(n_1-n̅_1)^2 (κΩ_1)^n_1+n̅_1/∑_n_1=0^∞∑_n̅_1=0^∞1/n_1!n̅_1!(κΩ_1)^n_1+n̅_1=2κ .

That is, when calculating the susceptibility from the partition function (<ref>) but evaluating Eq. (<ref>) in a finite subvolume Ω_1 of infinite Euclidean spacetime, χ_Ω_1 is nonzero. On the other hand, as lim_Ω→∞Δ n/Ω=0, in infinite volume we have

    χ_Ω→∞=0 .


One next introduces the Fourier transforms of the correlation functions, which we indicate here by a tilde:

    ⟨ q q⟩ (p)=
    ∫ d^4 (x-y)  e^ ip(x-y)⟨ q(x) q(y)⟩
    =p_μ p_ν⟨ K_μ K_ν⟩(p) .

From Eq. (<ref>), one can conclude the following infrared behaviour <cit.>:

    lim_p→ 0p_μ p_ν⟨ K_μ K_ν⟩(p)=
    ∫ d^4 x ⟨∂_μ K_μ(x) ∂_ν K_ν(y)⟩=χ .

The former equation implies e.g. in transverse gauge <cit.>, where ∂_μϵ_μνρσ A_σ=0,

    ⟨K_μ K_ν⟩(p) = χ p_μ p_ν/p^4+ O(p^2)

so that one would expect for χ≠0 a simple massless pole in ⟨ K_μ K_ν⟩(p) for p_μ→ 0. Here we have further dropped the spacetime volume subscript on χ.

Now per Eq. (<ref>), K_μ can be understood as the Hodge dual of a three-form field. So it is of interest to compare it with the action for a massless three-form L_μ,

    S=∫_Ω d^4 x  (∂_μ L_μ)^2 .

The expression “massless” refers here to the absence of a mass term in the action and is not supposed to indicate that there is a massless propagating degree of freedom, which in fact there is not. The variation of this action is given by

    δ S=-2∫_Ω d^4x (∂_ν∂_μ L_μ)δ L_ν+2∫_∂Ω da_ν (∂_μ L_μ)δ L_ν ,

where  a_ν is a normal surface element on Ω. Discarding the boundary term for the moment, one would conclude that ∂_μ L_μ=const. But then again, dropping the boundary term is in general not justified because it in general contributes to the equations of motion. The situation here is different from the model specified by the Lagrangian (<ref>) in an infinite spacetime volume Ω: There, solutions of finite action exist while the topological term does not need to vanish and can be written per Eqs. (<ref>) and (<ref>) as a boundary term. However, within a given topological sector the boundary term gives a fixed contribution - iθΔ n to the action. Therefore it has no impact on the equations of motion.

Back to the action (<ref>), we must therefore make further assumptions about the boundary term. While there are many different options, we follow Ref. <cit.> that works with solutions

    ∂_μ L_μ=θ_L ,

where θ_L is an integration constant. While not stated explicitly in Ref. <cit.>, this appears to require boundary conditions δ L_μ=0 on ∂Ω and

    θ_L=1/Ω∫_∂Ω da_ν  L_ν .

Interpreting for the moment L_μ as a four-dimensional electric field in the background of a constant charge density θ_L, it should be clear that Eq. (<ref>) has solutions that satisfy the condition (<ref>). Note that for a given θ_L, L_μ∼θ_L Ω^1/4. To keep the boundary conditions well defined in terms of finite L_ν, we therefore have to impose these on finite Ω or we have to take θ_L=0.

Given the boundary condition δ L_μ=0, we can also derive the equation of motion

    ∂_μ^x∂_ν^x⟨ L_ν(x) L_ρ(y)⟩=-δ_μρδ^4(x-y)

for the correlation function. It further implies

    ⟨(∂_μ L_μ(x))(∂_ν L_ν(y)) ⟩= 4δ^4(x-y)+const.

To arrive at the above equation, one uses the
translation invariance to deduce that the correlator is a function of (x-y). Then one can take ∂^x_μ = -∂^y_μ and contract μ and ρ in Eq. (<ref>).

Having stated all this, it is interesting to follow Ref. <cit.> and consider the action (<ref>) as an EFT of the full model (<ref>), of which all degrees of freedom but K_μ have been integrated out. The effective action should generally feature higher order terms in ∂_μ L_μ beyond Eq. (<ref>) as well as nonlocal expressions that account for the rich phenomenology of QCD compared to the present simplistic model <cit.>. Yet, it appears that the following characteristic features of QCD are recovered.

Equation (<ref>) indicates that the topological fluctuations correspond to white noise. In that sense, the model (<ref>) without higher order and nonlocal terms captures the infrared behaviour on scales where the fluctuations are indeed expected to be uncorrelated. We may use the δ-function term in relation (<ref>) to identify

    L_μ=2 K_μ/√(χ) .

Indeed, substituting the above into Eq. (<ref>) (ignoring the constant term for the moment) and taking a Fourier transform one recovers Eq. (<ref>).


On the basis of these correspondences, the argument of Ref. <cit.> then goes that the integration constant θ_L selects a vacuum state with a nonzero value of the CP-odd observable ∂_μ L_μ  that it is therefore proportional, at least at the linearized level, to the parameter θ in the microscopic theory, which would thus turn out to be physical.

However, as we shall see next, we disagree with this interpretation. Because of the different boundary conditions, the  EFT for L_μ cannot correspond to either of the QCD partition functions (<ref>) or (<ref>), due to reasons that will also make clear why θ_L is not proportional to an angle.

Translating the boundary condition (<ref>) back to the fundamental theory implies that one fixes 

    Δ n=∫_∂Ω da_ν  K_ν=1/16π^2∫_Ω d^4 x  FF .

Moreover, these boundary conditions do not imply that the physical fields vanish on ∂Ω, and in fact per Eq. (<ref>) they do not for θ_L≠0, because θ_L∝∂_μ L_μ∝∂_μ K_μ, which is gauge invariant. The present setup is therefore not an EFT for the partition function specified through Eqs. (<ref>) and (<ref>). For the latter, one samples over all integer values of Δ n and the physical fields vanish on ∂Ω.

A consequence of this discrepancy is the fact that from Eqs. (<ref>) and (<ref>), one concludes that ⟨Δ n⟩/Ω=2 i(-1)^ N_fκsinθ̅ is purely imaginary <cit.> whereas in the EFT and for the boundary conditions that are assumed here, Δ n is given by the fixed real value (<ref>). The two setups therefore are not only different, but they also predict distinctly inequivalent results for the observables. There should also be no relation between θ from the model of Eqs. (<ref>) and (<ref>) and θ_L, which explains why the EFT in terms of the latter does not show any periodic behaviour as would be required for an angular variable.

One should notice that while Δ n is fixed over the full volume Ω, it will nonetheless fluctuate in any subvolume, cf. Eq. (<ref>). Still, the boundary conditions on ∂Ω are required, and these are very different from those in Eq. (<ref>). Further, if we interpret L_μ as an effective infrared field, nonvanishing ∂_μ L_μ on ∂Ω does not conflict with vanishing microscopic fields on that surface. But following the above remark concerning real versus imaginary ⟨Δ n⟩, the expectation values for Δ n in Ω still do not agree between the EFT and the setup specified through Eqs. (<ref>) and (<ref>).

It is also clear now that the present EFT (<ref>) with the boundary condition (<ref>) is also not equivalent to the partition function (<ref>) that we take as the object to theoretically define the strong interactions. The EFT corresponds to imposing a fixed flux (<ref>)—not necessarily an integer—on a finite surface. This in contrast to the case of Eq. (<ref>), for which we see from Eq. (<ref>) that all integer values of the topological flux at infinity are accounted for in the action. Imposing a fixed flux on a finite surface cannot yield the correct vacuum correlations, not least since a general finite Euclidean ∂Ω does not have a geometrically meaningful continuation to some domain in Minkowski space.

Though as just has been argued the EFT is not compatible with the usual QCD partition functions,
one may nevertheless point out that the topological susceptibility in the EFT can have a similar infrared behaviour as in the theory (<ref>). To see this, consider the case θ_L=0 for simplicity. The generalization can be carried out by accounting for disconnected contributions to the correlations in Eq. (<ref>). Since the topological flux Δ n vanishes, this means by Eq. (<ref>) that χ_Ω=0. In fact, this requirement fixes the constant term in Eq. (<ref>). This clearly resembles the behaviour exhibited in Eq. (<ref>). For either setup, there is no problem with χ_Ω=0 on the full volume Ω because there is no inhibition of topological fluctuations in subvolumes. Note that in the EFT for the topological current, for θ_L≠0 it follows that χ_Ω=⟨Δ n⟩^2/Ω is purely disconnected (because the winding number in Ω is fixed through the flux on ∂Ω) and still does not correspond to the value that is locally observed in a subvolume and includes connected contributions.

Now in Ref. <cit.>, it is suggested that the calculations in Ref. <cit.> (i.e. using the partition function (<ref>)) are equivalent to choosing boundary conditions on a finite surface that lift the mass of the three-form, supposedly preventing to capture the true massless dynamics. However, as discussed above, the vanishing of the chiral susceptibility from the partition function (<ref>) for the full volume in Eq. (<ref>) also occurs for the EFT for the topological current, up to disconnected contributions. In either case, this has the implication that the massless simple pole in ⟨ K_μ K_ν⟩(p) disappears when evaluating Eq. (<ref>) for the full spacetime. In this sense, both models behave very similarly and no conclusion can be drawn that the EFT captures dynamics that Eq. (<ref>) does not. In particular, one cannot support the statement from Ref. <cit.>  that with  Eq. (<ref>), one loses crucial information about CP-violating vacuum states based on some ill-behaved infrared behaviour of e.g. χ compared to the EFT.

Though also for the three-form EFT, Eq. (<ref>) does not yield a massless simple pole, it is of interest to give a more direct argument for why the presence of such a pole is not necessary for consistency with observed topological fluctuations to start with. First observe that Eq. (<ref>) is only true when the spacetime integral is taken over the full spacetime; otherwise, the Fourier transform is not simply projected into its zero momentum limit. To see this, consider a topological susceptibility defined by integrating over a finite subvolume Ω_1⊂Ω

    χ_Ω_1=∫_Ω_1 d^4 x ⟨∂_μ K_μ(x) ∂_ν K_ν(y)⟩=∫_Ω_1^4 x ∫^4 p/(2π)^4  e^- p(x-y)p_μ p_ν⟨K_μ K_ν⟩(p) .

In an infinite spacetime Ω, when Ω_1=Ω is also infinite, the spacetime integral yields a delta-function and one recovers the result of Eq. (<ref>). For a finite Ω_1⊂Ω, however the topological susceptibility will receive contributions from the correlator ⟨K_μ K_ν⟩(p) evaluated at nonzero momentum.



Lattice calculations in frozen topologies (that have χ_Ω=0 up to disconnected contributions) still observe nonzero topological susceptibility when sampling over a subvolume of the full lattice <cit.>. Hence, there is the logical possibility of a topological susceptibility that vanishes for the full spacetime, yet is nonzero when defined over subvolumes. This is exactly what happens with QCD defined as in Eq. (<ref>), as can be seen in Eqs. (<ref>), (<ref>). Now with χ vanishing in the full spacetime volume, Eq. (<ref>) implies that the correlator ⟨ K_μ K_ν⟩(p) cannot have a massless pole. In Ref. <cit.> it is stated without proof that this implies a massive pole instead. Let us adopt this assumption or more loosely introduce m_K as an infrared regulator that may or may not correspond exactly to the effect from the finite-volume cutoff on the correlation function.  We may also have m_K stipulate the absence of the simple pole in Eq. (<ref>). Either way,

    ⟨K_μ K_ν⟩(p) = χ_0  g_μν/p^2+m_K^2 .

Provided such a massive pole is consistent with nonzero χ_Ω_1, there is then no motivation for a massless pole for the three-form EFT of Ref. <cit.>. Hence any conclusion derived from insisting on a massless pole in the EFT does not apply.

One should therefore demonstrate how one can reconcile the topological susceptibility being zero at infinite volume while remaining nonzero at finite subvolumes. It turns out that such behaviour can actually be derived from the massive correlator of Eq. (<ref>). For estimating the susceptibility defined over a finite region Ω_R of radius R, one can consider an integration with a cutoff function e^-R^2/r^2 which suppresses the fluctuations for r>R,

    χ_Ω_Λ≡∫^4p/(2π)^4∫_ℝ^4 d^4 x  e^-R^2/r^2 e^- p(x-y) p_μ p_ν⟨K_μ K_ν⟩(p) .

Substituting Eq. (<ref>) and considering y=0, the result is

    χ_Ω_R= χ _0 (1-R ^4 m^4/16  e^R ^2 m_K^2/4Ei(-1/4 m_K^2 R ^2)-1/4  R^2 m_K^2)∼{[ χ_0,   R≪1/m_K,;   0,   R≫1/m_K. ].

Above, Ei(z)=-∫_-z^∞ e^-t/t   d t is the exponential integral function.

As advertised, one gets a zero topological susceptibility in the full volume, but also constant susceptibility for subvolumes with associated length scales L≪ 1/m_K. This behaviour is compatible with lattice results <cit.> and matches the result of Eqs. (<ref>), (<ref>), which are based on the dilute instanton gas approximation. 




  
Taking θ as a perturbation and PCAC
Sometimes it is argued that θ can be seen to be physical without making topological considerations by referring to calculations of possible CP-violating effects using identities for partially conserved axial currents (PCAC), e.g. in the discussion of Section 3 of Ref. <cit.> or in Ref. <cit.>. (Note though that these papers do not state explicitly that their calculations would allow them to come to this result without making topological arguments.) However, we show here that such conclusions ultimately rely on the assumption ξ=θ for Eq. (<ref>). The answer to the question of whether ξ=α̅ or ξ=θ relies however on topology so that PCAC methods or related arguments cannot be used in order to bypass this step.

To review such points, we express the anomalous divergence of the axial  U(1) current as 

    ∂_μ∑_j=1^N_f⟨ψ_j^†γ^5γ_μψ_j⟩
    =2 N_f/16π^2 trFF + 2 ∑_j=1^N_f⟨ψ_j^†γ^5 m_j  e^ iα_i γ^5ψ_j⟩ .

Using this relation, we can therefore see that 

    ⟨ A| θ/16π^2 trFF |B⟩
    = -θ/N_f⟨ A | ∑_j=1^N_fψ_j^†γ^5 m_j  e^ iα_i γ^5ψ_j | B ⟩ ,

where we have dropped the total divergence of the axial current. We have inserted a factor θ on each side in view of the subsequent discussion which is however arbitrary at this point.

Now, we first follow Ref. <cit.> and carry out chiral field redefinitions so that in Eq. (<ref>) the mass matrices are purely real (and diagonal) and in particular α̅=0 and θ=θ̅. Moreover, one can take the freedom of non-anomalous SU(2)_A chiral field redefinitions to set α_j=0 individually. Assuming θ to be small one may view the left-hand side of Eq. (<ref>) as a perturbative insertion of the CP-odd topological term into a hadronic transition matrix element. When |A⟩ and |B⟩ are eigenstates of CP with opposite eigenvalues, a nonzero result would then signal CP violation.  

Before proceeding, we note that we could just as well have set θ=0, α̅=+θ̅ and treated the axial U(1)_A phase of the quark masses as the perturbation. This is an equivalent point of view that is taken e.g. in Ref. <cit.> and perhaps more true to the circumstance that the gluons have been integrated out in the chiral EFT so that the operator θ FF has been removed as well. Either way, the right-hand side of Eq. (<ref>) turns out as 

    -θ̅/N_f⟨ A | ∑_j=1^N_fψ_j^†γ^5 m_j ψ_j | B ⟩ .


This way, all explicit phases have been pulled out in terms of a single factor in front of the matrix element (<ref>). If in addition, the chiral condensate has no phase, i.e. in Eq. (<ref>) U_0= 1_N_f or, in other words, the condensate points into the real direction, the matrix element (<ref>) can be evaluated using standard reduction formulae as presented e.g. in Ref. <cit.> leading to a nonzero result.

So taking θ=0 and α̅ as a perturbation and under the additional assumption that the chiral condensate points into the real direction, U_0= 1_N_f, one sees here CP violation. This extra assumption, however, is exactly what characterizes the choice ξ=θ in Eq. (<ref>)—as opposed to the correct choice in  Eq. (<ref>)—as the ground state implied by ξ=θ=0 can always be chosen as U_0= 1_N_f. This can be seen to follow from Eqs. (<ref>), (<ref>), provided that one uses the freedom to perform field redefinitions  without changing θ=0 (or equivalently, without changing α̅=α_u+α_d) to set α_u,α_d as follows,

    α_u=    m_d/m_u+m_d α̅,
    
        α_d=    m_u/m_u+m_d α̅.

If instead ξ=-α̅, U_0 can be calculated according to Eqs. (<ref>), (<ref>). Then, U_0 will be complex and aligned with -α̅, leading to additional perturbative corrections that will compensate for Eq. (<ref>). That this cancellation has to happen is just a consequence of the fact that we can dial α̅ to zero in the chiral effective Lagrangian (<ref>) by a redefinition of the field U.


We can illustrate this point explicitly for the example of η^'→2π discussed in Refs. <cit.> as well as in Section <ref> of the present paper. Suppose we start with the theory with θ=0 and α̅=0. Then, ξ=0 for either order of limits and there is no CP violation. Next, we reintroduce α̅≠0 with |α̅|≪ 1 through the Lagrangian term

    δ L_θ̅=- iθ̅m_u m_d/m_u+m_d∑_i=1^N_fψ̅_i γ^5 ψ_i ,

where here N_f=2 for simplicity and we have written α̅ as θ̅ since θ=0 per the present assumptions. Equation (<ref>) follows from expanding the mass terms in the Lagrangian (<ref>) to first order in α_i, restricting to two flavours and substituting Eq. (<ref>). We could have used the approximate SU(2)_A symmetry to attribute the phases differently among the flavours u and d but the present form allows most straightforwardly to connect with the relations from Refs. <cit.> and Section <ref>. In particular, inserting the Lagrangian of Eq. (<ref>) between the vacuum and |η'π^0π^0⟩, using PCAC relations, leads to

    ⟨ 0|δ L| η^'π^0π^0⟩
     =m_u m_d/(m_u+m_d)^2m_π^2/f_πθ̅ .

The previous result would give the total contribution to the CP-violating matrix element under the crucial assumption (which is not made explicit in references such as  <cit.>, neither do these articles contain a derivation of the expectation value ⟨ U⟩) that there are no additional CP-odd phases from the quark condensates when taking matrix elements on the physical states. Equation (<ref>) matches Eq. (5) of Ref. <cit.> up to a factor of  √(2/3) which can be understood from the fact that the former reference included a third light quark which was not accounted for here.
Of course, this result from PCAC relations also follows directly from the EFT Lagrangian term (<ref>) when assuming ξ=θ=0, so that ξ+α_u+α_d=α_u+α_d=θ̅.

Rather than relying on assumptions, the values of the phases in the quark condensate are calculable, as was shown in our effective theory analysis in  Section <ref>, in which the condensate phases were included in Eq. (<ref>) from the beginning and determined by solving the tadpole conditions, cf. Eq. (<ref>).
Now, when choosing the evaluation of the partition function according to Eqs. (<ref>) and (<ref>) the chiral condensate U_0 aligns with θ. This follows from the fact that in this order of limits, one should identify ξ=θ in Eq. (<ref>). Then Eqs. (<ref>) and (<ref>) relate φ_u,φ_d to θ.   As θ=0 in the present setup, one obtains φ_u=φ_d=0, so that the chiral condensate induces no additional phases and Eq. (<ref>) then is the only contribution to the total matrix element, as mentioned before. If instead the evaluation is done as in Eq. (<ref>) which corresponds to an evaluation of the partition function (<ref>) on a connected integration contour, ξ=-α̅ in Eq. (<ref>), so that Eqs. (<ref>), (<ref>) give U_0= e^- iα̅. So this deviation from an orientation of the chiral condensate in the real direction must be accounted for, which leads to additional complex phases in matrix elements, and thus extra contributions to Eq. (<ref>). The existence of additional phases to those leading to the PCAC result of Eq. (<ref>) is captured by our EFT result (<ref>) in the following manner: The angle α̅=α_u+α (which gives the part that can alternatively be evaluated through the PCAC relation (<ref>)) is compensated by ξ=-α̅, leading to a vanishing CP-violating amplitude.



In summary, the arguments from Refs. <cit.>, while technically correct if it were U_0= 1_N_f for θ=0, does not demonstrate that this assumption actually holds. To assess this, one, after all, has to consider topology in order to derive the effective `t Hooft vertex and to sort out the correct limiting procedure, i.e. either Eq. (<ref>) or Eqs. (<ref>) and (<ref>). As  Eq. (<ref>) turns out to be correct, U_0 aligns with -α̅ instead of θ so that there is no CP violation.



  
Relation of the present arguments to some recent literature There are other recent papers that argue for CP conservation in the strong interactions. We do not discuss these here comprehensively but state why they substantially differ from the line of reasoning in the present work, or, respectively, if there is some prospect to establish connections.

In Ref. <cit.>, the starting assumption is a finite toroidal four-dimensional geometry, where the different topological sectors are weighted as exp( iΔ n θ). Because of the finite volume, this decisively differs from Eq. (<ref>). It is further reported in that work that for θ≠0 no confinement occurs so that the experimentally observed confinement must result from θ=0. We can neither confirm this latter statement nor is it central to the discussion in the present work.

In Ref. <cit.>, it is stated that topological charge, i.e. winding number is not observable and that no `t Hooft operator can be derived. We neither reproduce this, and the `t Hooft operator (<ref>) here is of physical consequence and central to the present discussion.

The paper <cit.> makes the point that coherence between topological sectors must not lead to observable consequences in the presence of causal horizons, which arguably only happens when θ=0. In the present work, it turns out that θ has no material consequence in infinite Euclidean spacetime and hence neither in Minkowski spacetime. However, the present result (<ref>) is suggestive in that the partition function in a finite subvolume with free boundary conditions of Euclidean space, as derived from Eq. (<ref>), does not exhibit the parameter θ any more. It would therefore be interesting to understand the possible relation between the Euclidean subvolumes and the domains within a causal horizon of Refs. <cit.>.




§ CONCLUSIONS



Topological quantization appears central to the correct assessment of the CP-odd θ parameter in QCD. It is therefore important to deduce its origin and to account for the implications of the setup of the problem. When we take Euclidean spacetime as the analytic continuation of Minkowski spacetime and do not impose ad hoc boundary conditions, time must be taken to infinity. Then, integration contours can be deduced that imply that we are first to evaluate the path integral in the individual topological sectors in infinite volume and then to sum over these sectors. As a consequence, there is no CP violation present in QCD, particularly not in the effective `t Hooft vertex. While this corresponds to a stringent reasoning within zero-temperature QCD, beyond the present work we shall next address finite-temperature QCD as a setup with a clear physical meaning in a finite spacetime volume.



§ ACKNOWLEDGEMENTS


The work of W.Y.A. is supported by the UK Engineering and Physical Sciences Research Council (EPSRC), under Research Grant No. EP/V002821/1.  The research of C.T. was
supported by the Cluster of Excellence Precision Physics, Fundamental Interactions,
and Structure of Matter (PRISMA+, EXC 2118/1) within the German Excellence
Strategy (Project-ID 390831469)

utphys







