Reflection positivity and complex analysis of the Yang-Mills theory from a viewpoint of gluon confinement

TL;DR

This paper investigates the confining properties of the Yang-Mills theory by analyzing the massive Yang-Mills model, demonstrating reflection positivity violation, and exploring the complex structure of the gluon propagator through both numerical and analytical methods.

Contribution

It establishes the equivalence of the massive Yang-Mills model with a gauge-invariant extension and confirms its ability to reproduce lattice results and explain gluon confinement features.

Findings

Gluon and ghost propagators match lattice data with one-loop corrections.

Reflection positivity is violated in the model due to complex conjugate poles.

The gluon propagator's complex structure explains the Gribov-Stingl form fit.

Abstract

In order to understand the confining decoupling solution of the Yang-Mills theory in the Landau gauge, we consider the massive Yang-Mills model which is defined by just adding a gluon mass term to the Yang-Mills theory with the Lorentz-covariant gauge fixing term and the associated Faddeev-Popov ghost term. First of all, we show that massive Yang-Mills model is obtained as a gauge-fixed version of the gauge-invariantly extended theory which is identified with the gauge-scalar model with a single fixed-modulus scalar field in the fundamental representation of the gauge group. This equivalence is obtained through the gauge-independent description of the Brout-Englert-Higgs mechanism proposed recently by one of the authors. Then, we reconfirm that the Euclidean gluon and ghost propagators in the Landau gauge obtained by numerical simulations on the lattice are reproduced with good accuracy from the massive Yang-Mills model by taking into account one-loop quantum corrections. Moreover, we demonstrate in a numerical way that the Schwinger function calculated from the gluon propagator in the Euclidean region exhibits violation of the reflection positivity at the physical point of the parameters. In addition, we perform the analytic continuation of the gluon propagator from the Euclidean region to the complex momentum plane towards the Minkowski region. We give an analytical proof that the reflection positivity is violated for any choice of the parameters in the massive Yang-Mills model, due to the existence of a pair of complex conjugate poles and the negativity of the spectral function for the gluon propagator to one-loop order. The complex structure of the propagator enables us to explain why the gluon propagator in the Euclidean region is well described by the Gribov-Stingl form.

Figures

Equations

$${\mathcal{L}}_{\mathrm{mYM}}^{\mathrm{tot}}=$$

$${\mathcal{L}}_{\mathrm{YM}}=$$

$${\mathcal{L}}_{\mathrm{GF}}=$$

$${\mathcal{L}}_{\mathrm{FP}}=$$

$$=$$

$${\mathcal{L}}_{\mathrm{m}}=$$

$$⟨\mathcal{O}[\mathcal{A}]{⟩}_{\mathrm{mYM}}:=$$

$$⟨\mathcal{O}[\mathcal{A}]{⟩}_{\mathrm{mYM}}=$$

$$=$$

$${\Delta }_{\mathrm{FP}}[\mathcal{A}]:=\det ({\partial }^{\mu }{\mathcal{D}}_{\mu }[\mathcal{A}{]}^{AB}).$$

$${\mathcal{L}}_{m\mathrm{YM}}^{\mathrm{tot}}=$$

$${\mathcal{L}}_{\mathrm{YM}}=$$

$${\mathcal{L}}_{\mathrm{GF}+\mathrm{FP}}=$$

$$+\mathcal{N}\cdot {\partial }^{\mu }{\mathcal{A}}_{\mu }-\frac{\beta }{2}g\mathcal{N}\cdot (i\bar{\mathcal{C}}\times \mathcal{C})$$

$$+i\bar{\mathcal{C}}\cdot {\partial }^{\mu }{\mathcal{D}}_{\mu }[\mathcal{A}]\mathcal{C}+\frac{\beta }{4}{g}^2(i\bar{\mathcal{C}}\times \mathcal{C})\cdot (i\bar{\mathcal{C}}\times \mathcal{C})$$

$$=$$

$$+\frac{\beta }{4}(\bar{\mathcal{N}}\cdot \bar{\mathcal{N}}+\mathcal{N}\cdot \mathcal{N}),$$

$${\mathcal{L}}_m=$$

$${\mathcal{L}}_{\mathrm{RF}}=$$

$${\mathcal{L}}_{\mathrm{YM}}=$$

$${\mathcal{L}}_{\mathrm{kin}}:=$$

$$f(\Phi (x)):=\Phi (x{)}^{†}\cdot \Phi (x)-\frac{1}{2}{v}^2=0,$$

$$\Phi (x)=$$

$${\mathcal{A}}_{\mu }(x)$$

$$\Phi (x)$$

$$\Theta (x):=$$

$$=$$

$${\mathcal{L}}_{\mathrm{RF}}=$$

$${\mathcal{L}}_{\mathrm{YM}}[\mathcal{A}]=$$

$${\mathcal{L}}_{\mathrm{kin}}[\mathcal{A},\Theta ]=$$

$${\mathcal{L}}_{\mathrm{c}}[u,\Theta ]=$$

$$f(\Theta ):=\mathrm{tr}\left({\Theta }^{†}\Theta -\frac{1}{2}{v}^{2}1\right)/\mathrm{tr}(1),$$

$${\mathcal{A}}_{\mu }(x)$$

$$\hat{\Theta }(x)$$

$$\hat{\Theta }(x)=\Theta (x)/\left(\frac{v}{\sqrt{2}}\right),v>0.$$

$$\hat{\Theta }(x)\in G=SU(2).$$

$${\mathcal{W}}_{\mu }(x):=$$

$$=$$

$${\mathcal{L}}_{\mathrm{kin}}[\mathcal{A},\hat{\Theta }]=$$

$$=$$

$${\mathcal{W}}_{\mu }(x)={\mathcal{A}}_{\mu }(x)-i{g}^{-1}\hat{\Theta }(x){\partial }_{\mu }\hat{\Theta }(x{)}^{†}.$$

$${\mathcal{W}}_{\mu }(x)=\hat{\Theta }(x){\mathcal{A}}_{\mu }^{{\hat{\Theta }}^{†}}(x)\hat{\Theta }(x{)}^{†},$$

Full text

††preprint: CHIBA-EP-238, 2020.02.06

I Introduction

It is still a challenging problem in particle physics to explain quark and gluon confinement in the framework of quantum gauge field theories YM54 . The very first question to this problem is to clarify what criterion should be adopted to understand confinement. For quark confinement, there is a well-established gauge-invariant criterion given by Wilson Wilson74 , namely, the area law falloff of the Wilson loop average leading to the linear static quark potential with a non-vanishing string tension. For gluon confinement, on the other hand, there is no known gauge-invariant criterion to the best of the authors’ knowledge. This is also the case for more general hypothesis of color confinement including quark and gluon confinement as special cases. Once the gauge is fixed, however, there are some proposals. For instance, the Kugo-Ojima criterion for color confinement is given for the Lorentz covariant Landau gauge KO79 . Indeed, it is rather difficult to prove the color confinement criterion even in a specific gauge, although there appeared an announcement for a proof of the Kugo-Ojima criterion for color confinement in the covariant Landau gauge CF18 . Even if color confinement is successfully proved in a specific gauge, this does not automatically guarantee color confinement in the other gauges. Therefore the physical picture for confinement could change gauge by gauge.

The information on confinement is expected to be encoded in the gluon and ghost propagators which are obtained by fixing the gauge. Recent investigations have confirmed that in the Lorentz covariant Landau gauge the decoupling solution decoupling-lattice ; decoupling-analytical is the confining solution of the Yang-Mills theory in the three- and four-dimensional spacetime, while the scaling solution is realized in the two-dimensional spacetime. Therefore, it is quite important to understand the decoupling solution in the Lorentz covariant Landau gauge. Of course, there are so many approaches towards this goal. In this paper, we focus on the approach TW10 ; RSTW17 ; Machado16 ; Kondo15 which has been developed in recent several years and has succeeded to reproduce some features of the decoupling solution with good accuracy. We call this approach the mass-deformed Yang-Mills theory with the gauge fixing term or the massive Yang-Mills model in the covariant gauge for short.

However, the reason why this approach is so successful is not fully understood yet in our opinion. In the original works TW10 the massive Yang-Mills model in the Landau gauge was identified with a special parameter limit of the Curci-Ferrari model CF76b . However, the Curci-Ferrari model is not invariant under the usual Becchi-Rouet-Stora-Tyutin (BRST) transformation, but invariant just under the modified BRST transformation which does not respect the usual nilpotency.

In this paper we show based on the previous works Kondo16 ; Kondo18 that the mass-deformed Yang-Mills theory with the covariant gauge fixing term has the gauge-invariant extension which is given by a gauge-scalar model with a single fixed-modulus scalar field in the fundamental representation of the gauge group, provided that a constraint called the reduction condition is satisfied. We call such a model the complementary gauge-scalar model. This equivalence is achieved based on the gauge-independent description Kondo16 ; Kondo18 of the Brout-Englert-Higgs (BEH) mechanism Higgs1 ; Higgs2 ; Higgs3 which does not rely on the spontaneous breaking of gauge symmetry NJL61 ; Goldstone61 . This description enables one to give a gauge-invariant mass term of the gluon field in the Yang-Mills theory which can be identified with the gauge-invariant kinetic term of the scalar field in the complementary gauge-scalar model.

In this paper, we first confirm that the massive Yang-Mills model with one-loop quantum corrections being included in the Euclidean region reproduces with good accuracy the gluon and ghost propagators of the decoupling solution of the Yang-Mills theory in the Landau gauge obtained by numerical simulations on the lattice. In fact, the resulting gluon and ghost propagators in the massive Yang-Mills model can be well fitted to those on the lattice by adjusting the parameters, namely, the coupling constant $g$ and the gluon mass parameter $M$.

For gluon confinement, the violation of reflection positivity is regarded as a necessary condition for confinement. In fact, it is known that the gluon propagator in the Yang-Mills theory exhibits the violation of reflection positivity. This fact was directly shown by the numerical simulations on the lattice, e.g., in the covariant Landau gauge CMT05 ; Bowmanetal07 . In this paper, by using the relevant gluon propagator in the massive Yang-Mills model, we calculate the Schwinger function in a numerical way to demonstrate that the reflection positivity is violated at the physical point of parameters reproducing the Yang-Mills theory.

In order to understand these facts and consider the meaning of gluon confinement, we perform the analytic continuation of the gluon and ghost propagators in the Euclidean region to those in the Minkowski region on the complex momentum squared plane. The consideration of the complex structure of the propagator enables us to give an analytical proof that the reflection positivity is violated for any choice of the parameters without restricting to the physical point of the Yang-Mills theory in the massive Yang-Mills model with one-loop quantum corrections being included. For this proof, it is enough to show that the Schwinger function necessarily becomes negative in some region, which is achieved by calculating separately the contributions to the gluon Schwinger function from the pole part and the continuous (branch cut) part of the gluon propagator based on the generalized spectral representation in the massive Yang-Mills model to one-loop order. It turns out that the violation of reflection positivity is an immediate consequence of the facts that the gluon propagator has a pair of complex conjugate poles and that the spectral function of the gluon propagator has negative value on the whole range, see HK18 . See e.g., DORS19 ; BT19 for the construction of the spectral function from the Euclidean data of numerical simulations on the lattice.

The complex structure of the propagator enables us to explain why the gluon propagator in the Euclidean region is well described by the Gribov-Stingl form Stingl86 , as demonstrated in the numerical simulations on the lattice DOS18 . Indeed, the pole part of the gluon propagator due to a pair of complex conjugate poles exactly reproduces the Gribov-Stingl form which is fitted to the numerical simulations to very good accuracy, after subtracting the small contribution coming from the continuous part represented by the spectral function obtained from the discontinuity across the branch cut on the positive real axis on the complex momentum plane. See also Kondo11 for another explanation for the occurrence of the gluon propagator of the Gribov-Stingl form.

The above result suggests that gluon confinement is not restricted to the confinement phase of the ordinary Yang-Mills theory, and can be extended into more general situations, namely, anywhere represented by the massive Yang-Mills model, which includes the Higgs phase in the complementary gauge-scalar model. In the lattice gauge theory, it is known that the confinement phase in the pure Yang-Mills theory is analytically continued to the Higgs phase in the relevant gauge-scalar model, which is called the Fradkin-Shenker continuity FS79 as a special realization of the Osterwalder-Seiler theorem OS78 . There are no local order parameters which can distinguish the confinement and Higgs phases. There is no thermodynamic phase transition between confinement and Higgs phases lattice-gauge-scalar-fund , in sharp contrast to the adjoint scalar case lattice-gauge-scalar-adj where there is a clear phase transition between the two phases. Therefore, confinement and Higgs phases are just subregions of a single confinement-Higgs phase FMS80 ; tHooft80 ; Maas17 . Therefore, permanent violation of positivity can be understood in light of the Fradkin-Shenker continuity between confinement-like and Higgs-like regions in a single confinement phase in the gauge-scalar model.

This paper is organized as follows. In sec. II, we introduce the massive Yang-Mills model in the covariant gauge. In sec. III, we show that the massive Yang-Mills model with quantum corrections to one-loop order well reproduces the gluon and ghost propagators of the decoupling solution. In sec. IV, we show that the gluon propagator exhibits violation of reflection positivity through the calculation of the Schwinger function. In sect. V, we perform the analytic continuation of the propagator to the complex momentum to examine the complex structure. In the final section we draw the conclusion and discuss the future problems to be tackled. In Appendix A, we give a recursive construction of the transverse and gauge-invariant gluon field to show the gauge-invariant extension of the massive Yang-Mills model. In Appendix B, we give another way for solving the reduction condition.

II Gauge-invariant extension of the mass-deformed Yang-Mills theory in the covariant Landau gauge

II.1

Mass deformation of the Yang-Mills theory in the covariant Landau gauge

We introduce the mass-deformed Yang-Mills theory in the covariant gauge which is defined just by adding the naive mass term $\mathscr{L}_{\rm m}$ to the ordinary massless Yang-Mills theory in the (manifestly Lorentz) covariant gauge fixing. The total Lagrangian density $\mathscr{L}_{\rm mYM}^{\rm tot}$ of the massive Yang-Mills model consists of the Yang-Mills Lagrangian $\mathscr{L}_{\rm YM}$, the gauge-fixing (GF) term $\mathscr{L}_{\rm GF}$, the associated Faddeev-Popov (FP) ghost term $\mathscr{L}_{\rm FP}$, and the mass term $\mathscr{L}_{\rm m}$,

[TABLE]

where $\mathscr{A}_{\mu}^{A}$ denotes the Yang-Mills field, $\mathscr{N}^{A}$ the Nakanishi-Lautrup field, and $\mathscr{C}^{A},\mathscr{\bar{C}}^{A}$ the Faddeev-Popov ghost and antighost fields, which take their values in the Lie algebra $\mathscr{G}$ of a gauge group $G$ with the structure constants $f_{ABC}$ ($A,B,C=1,...,{\rm dim}G$). We call this theory the massive Yang-Mills model in the covariant gauge for short.

The expectation value of an operator $\mathscr{O}[\mathscr{A}]$ of $\mathscr{A}_{\mu}^{A}$ is given according to the path integral quantization using the total action $S_{\rm mYM}^{\rm tot}[\mathscr{A},\mathscr{C},\mathscr{\bar{C}},\mathscr{N}]$ and the integration measure $\mathcal{D}\mathscr{A}\mathcal{D}\mathscr{C}\mathcal{D}\mathscr{\bar{C}}\mathcal{D}\mathscr{N}$

[TABLE]

In the Landau gauge $\alpha=0$, especially, the average is cast into a simpler form by integrating the Nananishi-Lautrup field $\mathscr{N}^{A}$ and subsequently the ghost and antighost field $\mathscr{C}^{A},\mathscr{\bar{C}}^{A}$ as

[TABLE]

with the Faddeev-Popov determinant,

[TABLE]

In this paper we do not intend to take into account the Gribov problem. The reasons are as follows. In this paper we deal with the massive Yang-Mills model as a low-energy effective model of the Yang-Mills theory and perform the perturbative analysis based on this model.

In the ultraviolet region the perturbative analysis of the Yang-Mills theory is valid due to the ultraviolet asymptotic freedom and is free from the Gribov problem, since the perturbative analysis can be done in the neighborhood of the origin of the configuration space of the gauge field within the first Gribov region and therefore does not reach the Gribov horizon where the Gribov problem becomes serious. This is also the case for the massive Yang-Mills model, since the effect of mass term can be ignored in the ultraviolet region.

Of course, in the usual perturbative treatment of the Yang-Mills theory, we encounter the Landau pole at which the gauge coupling constant diverges and the perturbative analysis breaks down at an intermediate momentum scale before reaching the deep infrared region. For the massive Yang-Mills model, however, we can adopt the infrared safe renormalization scheme in which the perturbation theory does not break down and remains valid from the large momentum all the way down to the zero momentum, as can be seen from the fact that the gauge coupling constant remains finite without divergence in the whole momentum region, and even vanishes in the zero momentum limit, as reviewed in section III. Therefore, we think that the massive Yang-Mills model can be treated in the whole region without seriously worrying about the Gribov problem, although there is no rigorous proof on this claim.

We regard the massive Yang-Mills model adopted in this paper as a low-energy effective model of the Yang-Mills theory where the mass term is generated in the dynamical way due to quantum corrections, for instance, according to the Wilsonian renormalization group. The mass term plays also the role of an infrared regulator and the massive Yang-Mills model is thereby free from the infrared divergence even in the vanishing momentum limit. Of course, the generation of the gluon mass term originates from non-perturbative effects and should be investigated from the first principles, which is however beyond the scope of this paper. Incidentally, we tried to show the existence of such mass term in WMNSK18 .

The massive Yang-Mills model just defined is a special case of a massive extension of the massless Yang-Mills theory in the most general renormalizable gauge having both BRST and anti-BRST symmetries given by Baulieu85

[TABLE]

where $\beta$ is a parameter which correspond to the gauge-fixing parameters in the $M\rightarrow 0$ limit, $\mathscr{D}_{\mu}[\mathscr{A}]\mathscr{C}(x):=\partial_{\mu}\mathscr{C}(x)+g\mathscr{A}(x)\times\mathscr{C}(x)$, and $\bar{\mathscr{N}}:=-\mathscr{N}+gi\bar{\mathscr{C}}\times\mathscr{C}$. This model is called the Curci-Ferrari model CF76b with the coupling constant $g$, the mass parameter $M$, and the parameter $\beta$. [In the Abelian limit with vanishing structure constants $f^{ABC}=0$, the FP ghosts decouple and the Curci-Ferrari model reduces to the Nakanishi model Nakanishi72 .] For $M\not=0$, the physics depends on the parameter $\beta$. This result should be compared with the $M=0$ case, in which $\beta$ is a gauge fixing parameter and hence the physics should not depend on $\beta$. In the $M=0$ case, indeed, any choice of $\beta$ gives the same physics. However, this is not the case for $M\not=0$. See e.g., Kondo13 for more details. The massive Yang-Mills model is regarded as a $\beta=0$ case of the Curci-Ferrari model. This point of view taken in the preceding works TW10 is good from the viewpoint of renormalizability, since the Curci-Ferrari model is known to be renormalizable. However, the Curci-Ferrari model lacks the physical unitarity at least in the perturbation theory CF76b ; Kondo13 . Indeed, the massive Yang-Mills model does not have the nilpotent BRST symmetry, although it has the modified BRST symmetry which does not respect the usual nilpotent property and reduces to the ordinary BRST symmetry only in the massless limit $M\to 0$. In this paper we try to find an extended theory with the ordinary nilpotent BRST symmetry, which reproduces the massive Yang-Mills model under an appropriate prescription. As a candidate for such a theory we investigate a specific gauge-scalar model.

In what follows we show that the massive Yang-Mills model in a covariant gauge has the gauge-invariant extension which is given by the gauge-scalar model with a single radially fixed (or fixed modulus) scalar field in the fundamental representation of a gauge group if the theory is subject to an appropriate constraint which we call the reduction condition. We call such a gauge-scalar model the complementary gauge-scalar model. In other words, the complementary gauge-scalar model with a single radially fixed scalar field in the fundamental representation reduces to the mass-deformed Yang-Mills theory in a fixed gauge if an appropriate reduction condition is imposed.

For $G=SU(2)$, the complementary gauge-scalar model is given by

[TABLE]

with a single fundamental scalar field $\Phi$ subject to the radially fixed condition,

[TABLE]

where $v$ is a positive constant $(v>0)$ and $\Phi(x)$ is the $SU(2)$* doublet* formed from two complex scalar fields $\bm{\phi}_{1}(x),\bm{\phi}_{2}(x)$,

[TABLE]

where $D_{\mu}[\mathscr{A}]$ is the covariant derivative in the fundamental representation $D_{\mu}[\mathscr{A}]:=\partial_{\mu}-ig\mathscr{A}_{\mu}$.

This gauge-scalar model is invariant under the gauge transformation,

[TABLE]

It is more convenient to convert the scalar field into the gauge group element. For this purpose, we introduce the matrix-valued scalar field $\Theta$ by adding another $SU(2)$ doublet $\tilde{\Phi}:=\epsilon\Phi^{*}$ as

[TABLE]

Then the complementary $SU(2)$ gauge-scalar model with a single radially fixed scalar field in the fundamental representation is defined by

[TABLE]

where $u$ is the Lagrange multiplier field to incorporate the holonomic constraint (7) written in the matrix form $f(\Theta)=0$. The radially fixed gauge-scalar model with the Lagrangian density (11) is invariant under the gauge transformation,

[TABLE]

Then we introduce the normalized matrix-valued scalar field $\hat{\Theta}$ by

[TABLE]

The above constraint (7) implies that the normalized scalar field $\hat{\Theta}$ obeys the conditions: $\hat{\Theta}(x)^{\dagger}\hat{\Theta}(x)=\hat{\Theta}(x)\hat{\Theta}(x)^{\dagger}=\bm{1},$ and $\det\hat{\Theta}(x)=1$. Therefore, $\hat{\Theta}$ is an element of $SU(2)$:

[TABLE]

This is an important property to provide a gauge-independent BEH mechanism.

The massive vector boson field $\mathscr{W}_{\mu}\in\mathscr{G}=su(2)$ is defined in terms of the original gauge field $\mathscr{A}_{\mu}\in\mathscr{G}=su(2)$ and the normalized scalar field $\hat{\Theta}\in G=SU(2)$ as shown in a previous paper Kondo18 ,

[TABLE]

According to the gauge-independent BEH mechanism Kondo16 ; Kondo18 , the kinetic term of the scalar field $\Theta$ is identical to the mass term of $\mathscr{W}_{\mu}$,

[TABLE]

The massive vector field $\mathscr{W}_{\mu}$ is rewritten using $\hat{\Theta}(x)\hat{\Theta}(x)^{\dagger}=\bm{1}$ into

[TABLE]

Then it is shown that the massive vector boson field $\mathscr{W}_{\mu}$ has the expression, 111 In Kondo18 $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}(x)$ was written as $\tilde{\mathscr{W}}_{\mu}(x)$ .

[TABLE]

where $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}$ denotes the gauge transform of $\mathscr{A}_{\mu}$ by $\hat{\Theta}\in G$. Notice that $\mathscr{W}_{\mu}$ transforms according to the adjoint representation under the gauge transformation,

[TABLE]

whereas $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}$ is gauge invariant,

[TABLE]

Therefore, the mass term can be written in terms of the gauge-invariant field $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}$ as

[TABLE]

This theory is supposed to obey the reduction condition for the massive vector field mode $\mathscr{W}_{\mu}(x)$. The stationary form of the reduction condition is given by

[TABLE]

where $\mathscr{D}_{\mu}[\mathscr{A}]$ is the covariant derivative in the adjoint representation $\mathscr{D}_{\mu}[\mathscr{A}]:=\partial_{\mu}-ig[\mathscr{A}_{\mu},\cdot]$. The stationary reduction condition is cast into

[TABLE]

This implies that imposing the reduction condition $\chi(x):=\mathscr{D}^{\mu}[\mathscr{A}]\mathscr{W}_{\mu}(x)=0$ is equivalent to imposing the “Landau gauge condition” $\partial^{\mu}\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}(x)=0$ or transverse condition for the gauge-invariant field $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}(x)$. Therefore, we can use the (gauge-transformed) reduction condition $\chi^{\hat{\Theta}^{\dagger}}$ written as

[TABLE]

and the associated Faddeev-Popov determinant $\Delta_{\rm FP}^{\rm red}$ reads

[TABLE]

Notice that the reduction condition $\chi^{\hat{\Theta}^{\dagger}}$ and the associated FP determinant $\Delta^{\rm red}$ are written in terms of $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}$ alone, $\chi^{\hat{\Theta}^{\dagger}}=\chi[\mathscr{A}^{\hat{\Theta}^{\dagger}}]$ and $\Delta^{\rm red}=\Delta^{\rm red}[\mathscr{A}^{\hat{\Theta}^{\dagger}}]$, and hence they are gauge invariant.

We show that the massive Yang-Mills (mYM) model in the Landau gauge can be converted to the complementary gauge-scalar (CGS) model, namely, radially fixed gauge-scalar model subject to the reduction condition. In fact, the vacuum expectation value of a gauge-invariant operator $\mathscr{O}[\mathscr{A}]$ of $\mathscr{A}_{\mu}^{A}$ reads

[TABLE]

where the normalized matrix-scalar field $\hat{\Theta}$ is introduced and the integration over the gauge volume $\int\mathcal{D}\hat{\Theta}^{\dagger}$ is inserted in the second equality, the integration variable $\mathscr{A}$ is renamed to $\mathscr{A}^{\hat{\Theta}^{\dagger}}$ in the third equality, the gauge invariance of the Yang-Mills action $S_{\rm YM}[\mathscr{A}^{\hat{\Theta}^{\dagger}}]=S_{\rm YM}[\mathscr{A}]$, the integration measure $\mathcal{D}\mathscr{A}^{\hat{\Theta}^{\dagger}}=\mathcal{D}\mathscr{A}$ and the operator $\mathscr{O}[\mathscr{A}^{\hat{\Theta}^{\dagger}}]=\mathscr{O}[\mathscr{A}]$ is used in the fourth equality, and the FP determinant $\Delta_{\rm FP}[\mathscr{A}]$ for the Landau gauge $\partial^{\mu}\mathscr{A}_{\mu}=0$ in the massive Yang-Mills model is identified with the FP determinant (25) for the reduction condition (24) in the fifth equality. In the last step, the delta function $\delta(\hat{\Theta}^{\dagger},h[\mathscr{A}])$ on the group $G$ satisfying $\int\mathcal{D}\hat{\Theta}^{\dagger}\delta(\hat{\Theta}^{\dagger},h[\mathscr{A}])=1$ is used to rewrite

[TABLE]

which is valid when the following equation for a given $\mathscr{A}_{\mu}$ has a unique solution of $h=h[\mathscr{A}]\in G$,

[TABLE]

This uniqueness of the solution corresponds to assuming that there are no Gribov copies if $\partial_{\mu}\mathscr{A}_{\mu}^{h[\mathscr{A}]}(x)=0$ is regarded as the gauge-fixing condition. Notice that we have taken into account the radially fixed constraint (7) in replacing the scalar field ${\Theta}^{\dagger}$ by the normalized matrix-valued (or group-valued) scalar field $\hat{\Theta}^{\dagger}$ in the last step.

We have assumed that the solution is unique in showing the equivalence in the above. Therefore, the equivalence is valid up to the Gribov copies. As mentioned already, however, we do not intend to seriously consider the Gribov problem in this paper, since we take the same standpoint as before explained in the above.

Incidentally, by adopting the absolute Landau gauge for $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}$ as the reduction condition, we can extract the gauge field configuration as the unique solution without Gribov copies. Then we can show the exact equivalence between the massive Yang-Mills model and a specific gauge-scalar model. Consequently, the resulting theory inevitably becomes nonloal as expected from the effective theory, which however does not affect the perturbative analysis done in this paper.

II.2

Solving the reduction condition

In the complementary gauge-scalar model, the scalar field $\Phi$ and the gauge field $\mathscr{A}$ are not independent field variables, because we intend to obtain the massive pure Yang-Mills theory which does not contain the scalar field $\Phi$. Therefore, the scalar field $\Phi$ is to be eliminated in favor of the gauge field $\mathscr{A}$. This is in principle achieved by solving the reduction condition as an off-shell equation, which is different from solving the field equation for the scalar field $\Phi$ as adopted in the preceding studies Stueckelberg38 ; KG67 ; SF70 ; Cornwall74 ; Cornwall82 ; DT86 . 222 See e.g., DTT88 ; RRA04 for reviews of the Stückelberg field. Notice that the reduction condition is an off-shell condition. Therefore, solving the reduction condition is different from solving the field equation for the Stückelberg field as done in the preceding works Capri-etal05 . This means that the solution of the reduction condition does not necessarily satisfy the field equation, while the solution of the field equation of the complementary gauge-scalar model automatically satisfies the reduction condition Kondo18 .

Consequently, the resulting massive Yang-Mills model with the covariant gauge-fixing term and the associated Faddeev-Popov ghost term becomes power-counting renormalizable in the perturbative framework, as demonstrated to one-loop order in the next section.

Moreover, the entire theory is invariant under the usual Becchi-Rouet-Stora-Tyutin (BRST) transformation $\delta_{BRST}$. The nilpotency $\delta_{BRST}\delta_{BRST}=0$ of the usual BRST transformations ensures the unitarity of the theory in the physical subspace of the total state vector space as the BRST-invariant sector according to Kugo and Ojima KO79 . This situation should be compared with the Curci-Ferrari model CF76b which is not invariant under the ordinary BRST transformation, but instead can be made invariant under the modified BRST transformation $\delta_{BRST}^{\prime}$. Nevertheless, this fact does not guarantee the unitarity of the Curci-Ferrari model due to the lack of usual nilpotency of the modified BRST transformation satisfying $\delta_{BRST}^{\prime}\delta_{BRST}^{\prime}\delta_{BRST}^{\prime}=0$, see e.g., Kondo13 .

We proceed to eliminate the scalar field $\Phi$ or $\Theta$ by solving the reduction condition to obtain the massive Yang-Mills model from the complementary gauge-scalar model

[TABLE]

Notice that introducing the reduction condition does not break the original gauge symmetry. The general form of the transverse and gauge-invariant Yang-Mills gauge field $\mathscr{A}_{\mu}^{h[\mathscr{A}]}$ satisfying (24) can be obtained explicitly by order by order expansion in powers of the gauge field $\mathscr{A}$ up to the Gribov copies. Indeed, $\mathscr{A}_{\mu}^{h[\mathscr{A}]}$ satisfying the transverse condition,

[TABLE]

is obtained as a power series in $\mathscr{A}$,

[TABLE]

where we have defined the transverse field $\mathscr{A}_{\mu}^{T}$ in the lowest order term linear in $\mathscr{A}$ as

[TABLE]

Then we find that the transverse field $\mathscr{A}_{\mu}^{h[\mathscr{A}]}$ is rewritten into

[TABLE]

Under an infinitesimal gauge transformation $\delta_{\Lambda}$ defined by $\delta_{\Lambda}\mathscr{A}_{\mu}=\mathscr{D}_{\mu}[\mathscr{A}]\Lambda:=\partial_{\mu}\Lambda-ig[\mathscr{A}_{\mu},\Lambda],$ $\Psi_{\nu}$ transforms as

[TABLE]

Therefore, $\mathscr{A}_{\mu}^{h}$ given by (33) is left invariant by infinitesimal gauge transformations order by order of the expansion,

[TABLE]

In Appendix A, we give a recursive construction of the transverse field $\mathscr{A}_{\mu}^{h[\mathscr{A}]}$ and the proof of gauge invariance of the resulting $\mathscr{A}_{\mu}^{h[\mathscr{A}]}$.

The mass term of $\mathscr{W}_{\mu}$ is equal to that of $\mathscr{A}_{\mu}^{\hat{\Theta}^{\dagger}}$,

[TABLE]

Therefore, the “mass term” of gauge-invariant field $\mathscr{A}_{\mu}^{h}$ is used to rewrite the kinetic term of the scalar field:

[TABLE]

In this way, we have eliminated the scalar field by solving the reduction condition.

Only when we adopt the covariant Landau gauge $\partial\cdot\mathscr{A}=0$ as the gauge-fixing condition, the infinite number of nonlocal terms disappear so that $S_{\rm kin}^{*}$ reduces to the naive mass term of $\mathscr{A}$,

[TABLE]

In the Landau gauge, thus, the complementary gauge-scalar model with the reduction condition reduces to the massive Yang-Mills model with the naive mass term.

The explicit expression of the massive vector field $\mathscr{W}_{\mu}$ in terms of $\mathscr{A}_{\mu}$ is given in Appendix B. Notice that $\mathscr{W}_{\mu}$ agrees with $\mathscr{A}_{\mu}^{T}=\mathscr{A}_{\mu}$ in the Landau gauge $\partial\cdot\mathscr{A}=0$.

III

Massive Yang-Mills model and decoupling solutions

In this section we give a review of the pertubative results TW10 ; RSTW17 obtained for the massive Yang-Mills model and reconfirm them from our viewpoint for later convenience.

In order to reproduce the decoupling solution of the Yang-Mills theory in the covariant Landau gauge, we calculate one-loop quantum corrections to the gluon and ghost propagators in the massive Yang-Mills model. The Nakanishi-Lautrup field $\mathscr{N}^{A}$ can be eliminated so that the gauge-fixing term reduces to

[TABLE]

The results in the Landau gauge is obtained by taking the limit $\alpha\to 0$ in the final step of the calculations. Only in the Landau gauge $\alpha=0$ the massive Yang-Mills model with a mass term $\mathscr{L}_{m}$ has the gauge-invariant extension. In order to obtain the gauge-independent results in the other gauges with $\alpha\not=0$, we need to include an infinite number of non-local terms in addition to the naive mass term $\mathscr{L}_{m}$ for gluons, as shown in the previous section.

III.1 Feynman rules for the massive Yang-Mills model

The Feynman rules for the massive Yang-Mills model are given as follows. The diagrammatic representations of the Feynman rules are given in Fig. 1.

(P1) gluon propagator $\langle\mathscr{AA}\rangle$

[TABLE]

(P2) ghost propagator $\langle\mathscr{C\bar{C}}\rangle$

[TABLE]

(V1) three-gluon vertex function $\langle\mathscr{AAA}\rangle$

[TABLE]

(V2) gluon-ghost-antighost vertex function $\langle\mathscr{AC\bar{C}}\rangle$

[TABLE]

(V3) four-gluon vertex function $\langle\mathscr{AAAA}\rangle$

[TABLE]

Here the momentum conservation is omitted and the momentum flow at each vertex is regarded as incoming, while the momentum of antighost as outgoing. Notice that the Feynman rules are the same as those of the ordinary Yang-Mills theory in the Lorenz gauge except for the gluon propagator which was replaced by the massive propagator (40).

The gluon propagator (40) has the same form as that in the renormalizable $R_{\xi}$ gauge where unitarity is not manifest. For any finite values of $\alpha$, the gluon propagator has good high-energy behavior, namely, the asymptotic behavior $O(1/k^{2})$ as $k\to\infty$, and hence the theory is renormalizable by power counting. For example, the choice $\alpha=1$ leads to the propagator $\frac{-1}{k^{2}-M^{2}}g_{\mu\nu}$. In the limit $\alpha\to\infty$, the gluon propagator reduces to the standard form for a massive spin-one particle, as can be easily seen in the second form. In the unitary gauge particle content is manifest, since there are no unphysical fields, and hence unitarity is transparent, while renormalizability is not transparent.

For any finite values of $\alpha$, the gluon propagator has an extra unphysical pole at $k^{2}=\alpha M^{2}$ besides the physical pole (massive gauge bosons) at $k^{2}=M^{2}$, as can be seen in the second form of (40). In order to preserve unitarity, the unphysical poles must be eliminated or mutually cancel in the $S$-matrix element involving only physical particles. In the spontaneously broken gauge theory, the would-be Nambu-Goldstone boson field has the propagator with the unphysical pole at $k^{2}=\alpha M^{2}$, and this unphysical pole of the would-be Nambu-Goldstone particle cancels one of the gauge boson in order to preserve unitarity. This is not the case in our model, since there are no Nambu-Goldstone particles without spontaneous symmetry breaking. The above type of cancellation of unphysical poles can be proven to all orders in perturbation theory by using the generalized Ward-Takahashi identities which are a consequence of the gauge invariance of the theory.

In the limit $\alpha\to 0$, however, the gluon propagator reduces to the simple form for a massive spin-one particle with the transverse projector $\frac{-1}{k^{2}-M^{2}}\left(g_{\mu\nu}-\frac{k_{\mu}k_{\nu}}{k^{2}}\right)$, as can be seen in the third form of (40), and the contribution from the unphysical pole at $k^{2}=0$ disappears in this limit. Therefore, the Landau gauge is the very special gauge which guarantees renormalizability and allows the existence of the gauge-invariant extension as demonstrated for the massive Yang-Mills model in the previous section.

III.2 One-loop quantum corrections and renormalization

We now take into account quantum corrections to the gluon and ghost propagators to one-loop order. In Fig. 2, we enumerate the one-loop diagrams which contribute to the gluon and ghost propagators to one-loop order.

In the massive Yang-Mills model we introduce the renormalization factors $Z_{\mathscr{A}},Z_{\mathscr{C}}=Z_{\mathscr{\bar{C}}},Z_{g},Z_{M^{2}},\tilde{Z}_{\alpha}$ to connect the bare unrenormalized fields (gluon $\mathscr{A}_{B}$, ghost $\mathscr{C}_{B}$ and antighost $\mathscr{\bar{C}}_{B}$) and bare parameters (the coupling constant $g_{B}$, the mass parameter $M_{B}$ and the gauge-fixing parameter $\alpha_{B}$) to the renormalized fields ($\mathscr{A}_{R}$, $\mathscr{C}_{R}$ and $\mathscr{\bar{C}}_{R}$) and renormalized parameters ($g_{R}$, $M_{R}$ and $\alpha_{R}$) respectively Taylor71 ; BSNW96 ; Wschebor08 :

[TABLE]

For comparison with the lattice data, we move to the Euclidean region and use $k_{E}$ to denote the Euclidean momentum so that $k^{2}=-k_{E}^{2}$.

For gluons, we introduce the two-point vertex function $\Gamma_{\mathscr{A}}^{(2)}$ as the inverse of the transverse part $\mathscr{D}_{\rm T}$ of the propagator 333 In this paper we focus on the Landau gauge. For the gluon propagator, therefore, we discuss the transverse part alone.

and the vacuum polarization function $\Pi_{T}$ as

[TABLE]

where $\delta_{Z}$ and $\delta_{M^{2}}$ are counterterms to cancel the divergence coming from the vacuum polarization function $\Pi_{T}$ to obtain the finite renormalized one ${\Pi}_{T}^{\rm fin}$

[TABLE]

under the suitable renormalization conditions to be discussed shortly, and they are related to the renormalization factors as

[TABLE]

We define the dimensionless versions $\hat{\mathscr{D}}_{\rm T}(s)$ and $\hat{\Pi}(s)$ of $\mathscr{D}_{\rm T}(k_{E})$ and $\Pi_{T}(k_{E}^{2})$ with the hat respectively

[TABLE]

with the dimensionless squared momentum

[TABLE]

and

[TABLE]

The gluon vacuum polarization function in the covariant Landau gauge $\alpha=0$ calculated using the dimensional regularization in Euclidean space is given to one-loop order as the power-series Laurent expansion in $\epsilon:=2-\frac{D}{2}$ 444 These expressions are obtained by taking the limit $\alpha\to 0$ of those with an arbitrary $\alpha$ given in Suda18 .

[TABLE]

where $C_{2}(G)$ is the quadratic Casimir operator of a gauge group $G$, $\gamma$ is the Euler constant, and $\eta$ is the value of $s$ at the scale $\tilde{\mu}$ introduced through the dimensional regularization for dimensional reasons

[TABLE]

Here we have defined the functions of $s$,

[TABLE]

Notice that there are no singular term in the finite part $\hat{\Pi}_{T}^{\rm fin}(s)$ even at $s=0$, because there does not exist $O(s^{-2})$ term in the bracket […] of (52), since the expansion of $h(s)$ around $s=0$ reads

[TABLE]

which follows from

[TABLE]

Thus, we have the finite part of the gluon vacuum polarization to one-loop

[TABLE]

which has the $s=0$ limit,

[TABLE]

III.3 Naive (zero-momentum) renormalization conditions

For gluons, we can take a naive vanishing-momentum renormalization condition such that

[TABLE]

The first renormalization condition adopted by Tissier and Wschebor TW10 is the vanishing-momentum renormalization condition which is written in terms of $\Gamma_{\mathscr{A}}^{(2)}$ or equivalently $\hat{\Pi}_{T}^{\rm fin}$ as

[TABLE]

where we have introduced the dimensionless ratio of the renormalization scale $\mu$ to the mass defined by

[TABLE]

Adopting the renormalization condition [TW1], we obtain the renormalized gluon vacuum polarization function,

[TABLE]

Note that constant terms in $[...]$ are canceled by the subtraction: $-(s\rightarrow\nu)$.

However, it has been shown TW10 that the vanishing-momentum renormalization condition (59): $\Gamma_{\mathscr{A}}^{(2)}(k_{E}=0)=M^{2}$ or $\hat{\Pi}^{\rm fin}(s=0)=0$ yields the infrared Landau pole, namely, the coupling constant diverging at a certain momentum in the infrared region. Therefore, we use another renormalization condition given in the next subsection.

III.4 Infrared safe renormalization condition

For ghost, we introduce the two-point vertex function $\Gamma_{gh}^{(2)}$, the propagator $\Delta_{gh}$ and the self-energy function $\Pi_{gh}$,

[TABLE]

where $\delta_{C}$ is a counterterm to cancel the divergence coming from the ghost self-energy function $\Pi_{gh}$ to obtain the finite one $\Pi_{gh}^{\rm fin}$.

[TABLE]

and is related to the renormalization factor as

[TABLE]

We also define the dimensionless versions $\hat{\Delta}_{gh}(s)$ and $\hat{\Pi}_{gh}(s)$ of $\Delta_{gh}(k_{E}^{2})$ and $\Pi_{gh}(k_{E}^{2})$ as

[TABLE]

with

[TABLE]

The ghost self-energy function $\Pi_{gh}(k)$ in the covariant Landau gauge $\alpha=0$ is also calculated using the dimensional regularization and the dimensionless version $\hat{\Pi}_{gh}(s)$ is given to one-loop order by

[TABLE]

For ghosts, we impose the renormalization condition

[TABLE]

The renormalization condition (69) determines the counterterm $\delta_{C}$ as

[TABLE]

Then we obtain the renormalized ghost self-energy function under the renormalization condition (69)

[TABLE]

We now return to the gluon renormalization. To avoid the infrared Landau pole for the coupling, we replace the vanishing-momentum renormalization condition (59) by the second one:

[TABLE]

There is a well-known non-renormalization for the coupling in the Taylor scheme Taylor71 which also holds in the massive Yang-Mills model in the Landau gauge: The identity

[TABLE]

implies in the Landau gauge

[TABLE]

since in the Landau gauge,

[TABLE]

The implication of the first renormalization condition of (72) is explained as follows. For the massive Yang-Mills model in the Landau gauge $\alpha=0$ as a special limit of the Curci-Ferrari model, the non-renormalization theorem holds in the sense that a combination of renormalization factors is finite to all orders in the loop expansions BSNW96 ; Wschebor08 : The identity

[TABLE]

implies in the Landau gauge

[TABLE]

As $Z_{M^{2}}Z_{\mathscr{A}}=1+\delta_{M^{2}}$ from (48) and $Z_{\mathscr{C}}=1+\delta_{C}$ from (65), the non-renormalization theorem (76) in the Landau gauge reduces to the relation between the counterterms

[TABLE]

which means in the one-loop level

[TABLE]

This is the result of the first renormalization condition of (72).

Then the remaining $\delta_{Z}$ is determined from the second renormalization condition of (72): $\hat{\Pi}_{T}^{\rm fin}(s=\nu)=\hat{\Pi}_{T}(s=\nu)+\nu\delta_{Z}^{(1)}-\delta_{C}^{(1)}=0$ by using $\delta_{C}^{(1)}=(\ref{delta_C})$ as

[TABLE]

Then, by substituting (80) and (79) into (51): $\hat{\Pi}_{T}^{\rm fin}(s)=\hat{\Pi}_{T}(s)+s\delta_{Z}^{(1)}+\delta_{M^{2}}^{(1)}$, the renormalized gluon vacuum polarization function is modified into HK18

[TABLE]

The gluon vacuum polarization at $s=0$ has a positive value

[TABLE]

where we have used the fact that $f(s)$ is a monotonically increasing function of $s$ with $f(0)=\frac{5}{2}$.

We enumerate the obtained renormalization factors as functions of $g^{2}$ and $\nu$

[TABLE]

and

[TABLE]

We can obtain the renormalization group functions using these renormalization factors. For instance, the anomalous dimension of the field $\Phi$ is obtained from the renormalization factor $Z_{\Phi}=1+Z_{\Phi}^{(1)}+\cdots$ according to

[TABLE]

where the replacement of the derivative with respect to $\mu^{2}$ by $\nu=\mu^{2}/M^{2}$ is valid to one-loop order, since $M$ is the renormalized mass which depends on the renormalization scale $\mu$. Therefore, the ghost field has the anomalous dimension to one-loop order

[TABLE]

Similarly, the anomalous dimension of the gluon field is calculated to one-loop order as

[TABLE]

Notice that $\gamma_{\mathscr{C}}$ is always negative ($\gamma_{\mathscr{C}}=0$ at $\nu=0$). We find that $\gamma_{\mathscr{A}}$ is negative for $\nu>0.28$, becomes zero at $\nu\sim 0.28\sim 0.53^{2}$ and positive for $\nu<0.28$ ($\gamma_{\mathscr{A}}=1/3$ at $\nu=0$).

The $\beta$ function for the gauge coupling constant is obtained from

[TABLE]

which is indeed calculated to one-loop as

[TABLE]

This equation is rewritten into a differential equation with respect to $\nu$

[TABLE]

Thus, by introducing the indefinite integral $W$ of $w$ which has the closed form

[TABLE]

the running gauge coupling constant is given by

[TABLE]

Notice that $W$ has the asymptotic expansions for $\nu\gg 1$ and $\nu\ll 1$ respectively

[TABLE]

Hence, the running gauge coupling constant behaves in the ultraviolet region $\nu\gg 1$ and infrared one $\nu\ll 1$ respectively as

[TABLE]

In the ultraviolet region $\nu\gg 1$ or $\mu\gg 1$, the beta function $\beta_{\tilde{g}^{2}}$ in the massive Yang-Mills model is negative for $\nu\gg 1$, since (91) has the expansion for $\nu\gg 1$

[TABLE]

This result is in agreement with the standard, universal beta function of the usual Yang-Mills theory reflecting the ultraviolet asymptotic freedom

[TABLE]

In the infrared region $\nu\ll 1$ or $\mu\ll 1$, on the other hand, the beta function $\beta_{\tilde{g}^{2}}$ of the massive Yang-Mills model becomes positive in the deep infrared regime, since (91) has the expansion for $\nu\ll 1$

[TABLE]

This implies that the running coupling constant $g^{2}(\mu)$ decreases towards the infrared region and vanishes as $\mu\to 0$

[TABLE]

Therefore the RG flow drives the system towards a weak coupling region as $\mu$ goes to zero. This fact justifies the use of the one-loop approximation to study the Yang-Mills theory even in infrared region. See Fig. 3.

We find that the beta function $\beta_{g^{2}}$ is negative for $\nu>0.07$, becomes zero at $\nu\sim 0.07\sim 0.26^{2}$ and positive for $\nu<0.07$. This implies that the running coupling constant $g^{2}(\mu)$ of the decoupling solution increases monotonically in decreasing the scale $\mu$ until $\mu$ reaches the value $\mu/M\sim 0.26$, and it turns over at $\mu/M\sim 0.26$ and decreases towards the infrared limit $g^{2}(\mu)\to 0$ as $\mu\to 0$.

Finally, we study the RG flow in the two-dimensional parameter space ($\nu:=\mu^{2}/M^{2}$, ${\tilde{g}^{2}}:=\frac{g^{2}C_{2}(G)}{16\pi^{2}}$) of the four-dimensional massive Yang-Mills model. See Fig. 4. First, we fix the value of $\nu$ to a relatively large value $\nu_{0}$ (which is equivalent to set $\tilde{m}^{2}:=M^{2}/\mu^{2}$ to a relatively small value), e.g., $\nu_{0}=100^{2}$ and varies the value of $\tilde{g}^{2}(\nu)$ to see the differences of the resulting trajectories. Then we find that the running coupling constant $\tilde{g}^{2}(\nu)$ remains finite for all $\nu$ if the initial value $\tilde{g}^{2}(\nu_{0})$ at $\nu_{0}$ is smaller than and equal to a certain value $\tilde{g}_{*}^{2}(\nu_{0})$, while it diverges at a finite $\mu$ if $\tilde{g}^{2}(\nu_{0})$ is greater than the value $\tilde{g}_{*}^{2}(\nu_{0})$. Therefore, the decoupling solution exists for $\tilde{g}^{2}(\nu_{0})<\tilde{g}_{*}^{2}(\nu_{0})$, while the scaling solution is realized at the critical value $\tilde{g}^{2}(\nu_{0})=\tilde{g}_{*}^{2}(\nu_{0})$ RSTW17 . For $\tilde{g}^{2}(\nu_{0})>\tilde{g}_{*}^{2}(\nu_{0})$, we have an infrared Landau pole. Therefore, the coupling constant behaves in decreasing $\nu$ from $\nu_{0}$, $\nu_{0}>\nu=\frac{\mu^{2}}{M^{2}}\searrow 0$ as

[TABLE]

III.5 Fitting to the numerical simulations

We utilize the data obtained by the numerical simulations on the lattice for the Yang-Mills theory in the covariant Landau gauge to determine the parameters, the coupling constant $g$ and the gluon mass parameter $M$, in the massive Yang-Mills model.

In fitting the data of numerical simulations for the gluon propagator on the lattice DOS16 to the analytical expression $\mathscr{D}$ for the gluon propagator with one-loop quantum corrections, we need to take into account the fact that the renormalization conditions adopted in the lattice simulations DOS16 are different from those adopted in this paper, leading to the different scale for the gluon propagator. Otherwise, the fitting does not work so well and the appropriate parameters cannot be obtained. For this purpose, we introduce an overall scale factor $Z$ which can scale the gluon propagator as a whole to absorb the difference of the renormalization conditions. In DOS16 , indeed, such a scaling of data obtained by numerical simulations for the gluon propagator was adopted to satisfy the renormalization condition $\mathscr{D}_{T}(k_{E}^{2}=\mu^{2})=1/\mu^{2}$ at $\mu=4$ GeV. This kind of rescaling was also adopted in TW10 . Consequently, the fitting works surprisingly well to give the precise values for the parameters $g$, $M$ and $Z$ as shown in Fig. 5 in the fitting range $0<k_{E}\leq 4$GeV at $\mu=1$GeV for $G=SU(3)$ where the fitting parameter with errors are given by

[TABLE]

We use these parameters to plot the ghost propagator using the analytical expression by including quantum corrections to one-loop order in the massive Yang-Mills model, as shown in Fig. 6.

Both gluon propagator and ghost propagator in the decoupling solution of the Yang-Mills theory are well reproduced by the values (102) of parameters $g$ and $M$. In what follows we call these values of the parameters the physical point for the Yang-Mills theory.

As a side remark, let us add some comments on the validity of the massive Yang-Mills model in reproducing the infrared behaviors of the Yang-Mills theory. There is no guarantee in advance that such a specific model with a “phenomenological” mass term for gluons being just included captures the intricacies of the real Yang-Mills dynamics. We acknowledge that the surprising agreement between the numerical lattice data of the Yang-Mills theory and the simple one-loop propagator of the massive Yang-Mills model could be accidental, and that the gluon mass term will, at best, only capture some aspects, not all aspects, of the intricate dynamics of the original Yang-Mills theory or QCD. In fact, this type of the massive model for the real QCD is shown to give a poor agreement for the quark sector of QCD with numerical lattice results PTW14 . Nevertheless, we can still claim that this model gives a gluon propagator showing excellent agreement with the lattice data. Indeed, it is shown GPRT19 that the two-loop calculations for the gluon and ghost propagators considerably improve the one-loop result to show more excellent agreement with the lattice data. In these investigations, it is also confirmed that the pure Yang-Mills sector indicates the infrared-safety, namely, the finiteness of the running gauge coupling constant in all scales, which makes the perturbative method more feasible. Incidentally, the one-loop calculation for the three-point gluon vertex functions gives a “satisfying” agreement with the available lattice data PTW13 . In view of these works, the massive Yang-Mills model will be valid to capture some aspects of the gluon sector of QCD relevant to our investigation, even though the other important aspects may be missing. At least for the gluon, therefore, it will be worthwhile to study the analytic structure of the propagator of this model, which is one of our purposes in this paper.

IV

Reflection positivity violation in the massive Yang-Mills model

In this section, we observe that the Euclidean gluon propagator in the massive Yang-Mills model exhibits violation of reflection positivity. This result suggests gluon confinement in the Yang-Mills theory.

Usually the quantum field theory (QFT) obeying the Wightman axioms Wightman56 ; AQFT is first defined in the Minkowski region and then analytically continued to the Euclidean region to obtain the Euclidean QFT which consequently obeys the Osterwalder-Schrader (OS) axioms OS73 . However, we want to start from the Euclidean QFT obeying the OS axioms (or better axioms if any) and check which kinds of QFT can be defined in the Minkowski spacetime which is to be obtained by analytic continuation from the Euclidean region.

In our opinion, only the Euclidean QFT can be rigorously defined as the QFT. Probably, QFT describing only non-confining particles will be defined both in the Euclidean and the Minkowski space in the equivalent way. However, we have no evidences that the QFT describing confining particles can be formulated in the Minkowski spacetime in the same way as QFT for non-confining particles. In contrast, we know some examples of Euclidean QFT which exhibit confinement, e.g., the linear potential for the static quark potential is observed in the Euclidean Yang-Mills theory on the lattice. Therefore, the validity of the Euclidean QFT for confining particles is tested everyday on the lattice in the non-perturbative manner. In view of these, we examine the validity of the reflection positivity as an axiom or one of the general properties to be satisfied by the Euclidean QFT.

IV.1 Reflection positivity and the Schwinger function

The OS axioms OS73 are general properties to be satisfied for the QFT formulated in the Euclidean space, which are the Euclidean version of the Wightman axioms for the relativistic QFT formulated in the Minkowski spacetime. A relativistic QFT described by a set of the Wightman functions satisfying the Wightman axioms can be constructed from a set of Schwinger functions (Euclidean Green’s functions) if they obey the OS axioms. In particular, the axiom of reflection positivity is the Euclidean counterpart to the positive definiteness of the norm in the Hilbert space of the corresponding Wightman QFT. If the reflection positivity is violated, a particular Euclidean correlation function cannot have the interpretation in terms of stable particle states, which is regarded as a manifestation of confinement. To demonstrate the violation of reflection positivity in the OS axioms, one counterexample suffices.

For the special case of a single propagator, the reflection positivity reads

[TABLE]

where $\mathscr{S}_{+}(\mathbb{R}^{D})$ denotes a complex-valued test (Schwartz) function with support in $\{(\bm{x},x_{D});x_{D}>0\}$. The reflection positivity is rewritten as

[TABLE]

where we defined $\Delta(\bm{p},x_{D}-y_{D})$ by

[TABLE]

In what follows we call $\Delta(\bm{p},x_{D}-y_{D})$ the Schwinger function. For this inequality to hold for any test function $f\in\mathscr{S}_{+}(\mathbb{R}^{D})$, the Schwinger function $\Delta$ must satisfy the positivity

[TABLE]

We consider a particular Schwinger function in the $D$-dimensional spacetime defined by the Fourier transform of the Euclidean propagator $\tilde{\mathscr{D}}(\bm{p},p_{E}^{D})$,

[TABLE]

If $\tilde{\mathscr{D}}(\bm{0},p_{E}^{D})$ is even in $p_{E}^{D}$, namely, $\tilde{\mathscr{D}}(\bm{0},-p_{E}^{D})=\tilde{\mathscr{D}}(\bm{0},p_{E}^{D})$, the Schwinger function reduces to

[TABLE]

To demonstrate the violation of reflection positivity, one counterexample suffices. Therefore, non-positivity of the Schwinger function $\Delta(t)$ at some value of $t$ leads to the violation of reflection positivity. Consequently, the reflection positivity is violated for the gluon propagator. The corresponding states cannot appear in the physical particle spectrum. This is consistent with gluon confinement.

For the free massive propagator,

[TABLE]

we find $\Delta(t)$ is positive for any $t$:

[TABLE]

Therefore, there is no reflection-positivity violation for the free massive propagator, as expected. For unconfined particles, the reflection positivity should hold.

IV.2 Positivity violation for the decoupling solution of the Yang-Mills theory

In order to examine the violation of the reflection positivity through the behavior of the gluon Schwinger function, we first construct a set of gluon and ghost propagators in such a way that they are renormalized to satisfy the renormalization conditions TW2 and (69) in the massive Yang-Mills model to reproduce the decoupling solution in the Yang-Mills theory to one-loop order. The integral in obtaining the Schwinger function as the Fourier transform of the gluon propagator is not so easy to be performed analytically, hence we resort to the numerical calculations for this definite integral.

In Fig. 7, we give the plot for the gluon propagator and the associated Schwinger function in the Landau gauge $\alpha=0$ for the $SU(3)$ massive Yang-Mills model at the physical point of parameters $g=4.1$ and $M/\mu=0.454$. We observe that the Schwinger function takes negative values for $\mu t>6$ and hence the reflection positivity is violated. Therefore, this result suggests that the reflection positivity is violated for the decoupling solution in the Yang-Mills theory. The more detailed analysis of the reflection positivity will be given in the next section from the viewpoint of the complex structure of the gluon propagator.

IV.3 Positivity violation in the complementary gauge-scalar model

In what follows, we examine how the gluon propagator and the Schwinger function are modified if the parameters $g$ and $M$ deviate from the physical point. In this case the massive Yang-Mills model is no longer regarded as a low-energy effective theory of the original Yang-Mills theory. However, the resulting model can be regarded as the gauge-scalar model with the complementarity between Higgs and confinement in the sense that the confinement phase in the Yang-Mills theory is analytically connected with no phase transition to the Higgs phase in the gauge-scalar model through the BEH mechanism, which is called the Fradkin-Shenker continuity.

IV.3.1 Smaller coupling constant

First, we take smaller values for the coupling constant $g$ than the physical value $g=4.1$ and keep the mass parameter $M$ fixed to the physical value $M/\mu=0.454$. In Fig. 8, the gluon propagator and the associated Schwinger functions are given for a smaller value $g=2.3$. For a further smaller value $g=1$, they are given in Fig. 9.

For smaller coupling constant $g$, the gluon propagator ${\mathscr{D}}$ seems to be monotonically decreasing in $k_{E}$. The Schwinger function falls off very slowly from $t=0$ value and keeps its positivity until very large value of $t$, although it is difficult to see the difference from the graphs. Consequently, the smallest value of $t$ giving the negative value of the Schwinger function shifts to larger values of $t$, and eventually goes to infinity as $g\to 0$. This result is reasonable, since, in the vanishing coupling limit $g\to 0$, the gluon propagator must reduce to the free massive propagator in the tree level. Therefore, the reflection positivity must be recovered and the Schwinger function keeps positivity everywhere in the limit $g\to 0$. As far as the results of the numerical calculations are concerned, the positivity seems to be not violated and restored for relatively smaller coupling constants.

However, this observation turns out to be wrong. In fact, we can prove analytically that the reflection positivity of the gluon Schwinger function is violated for any value of the parameters $g$ and $M$ in the massive Yang-Mills model with one-loop quantum corrections being included. The proof will be given in the next section. The Schwinger function $\Delta$ is an oscillating and exponentially fall-off function of $t$ approaching zero finally as $t\to\infty$. Therefore, it is difficult to examine the violation of positivity in the large $t$ region in the numerical way due to the restriction on the precision of numerical calculations. For smaller coupling constant $g$, therefore, the Schwinger function takes a smaller but negative value for larger $t$, until the negativity disappears only in the limit $g\to 0$.

IV.3.2 Smaller mass parameter

Next, we keep the coupling constant fixed to the physical value $g=4.1$, and take smaller gluon mass parameter $M/\mu$ than the physical value $M/\mu=0.454$. For a smaller value $M/\mu=0.2$, the gluon propagator and the associated Schwinger functions are given in Fig. 10. For a further smaller value $M/\mu=0.141$, they are given in Fig. 11.

As the value of mass parameter $M/\mu$ is chosen to be smaller and smaller than the physical value for the Yang-Mills theory, the gluon propagator $\tilde{\mathscr{D}}(p)$ exhibits sizable non-monotonic behavior and the Schwinger function exhibits more enhanced negativity, leading to the clearer violation of reflection positivity.

For smaller mass $M$ or larger coupling constant $g$ than the physical value for the Yang-Mills theory, the gluon propagator $\tilde{\mathscr{D}}_{\rm T}(p)$ exhibits stronger non-monotonic behavior.

IV.3.3 Presence of Euclidean poles

For quite small mass parameter $M^{2}/\mu^{2}$ or large coupling constant $g$, the gluon propagator $\tilde{\mathscr{D}}(k_{E}^{2})$ becomes singular at two values of $k_{E}^{2}$ and takes negative values in between. In Fig. 12, the gluon propagator is given for the parameters $g=4.1$ and $M/\mu=0.08$. This result is consistent with the statement HK18 that the gluon propagator has poles in the Euclidean region (namely, tachyonic poles) with multiplicity two or a pair of complex conjugate poles under some assumptions on the propagator and the spectral function. The related issue will be discussed in the next section.

Therefore, this singular behavior affects the associated Schwinger function $\Delta(t)$. This feature will be an artifact due to the limitation of one-loop calculations. Therefore, we exclude the relevant region of parameters from the following considerations.

IV.4 Magnitude of positivity violation and the complementary gauge-scalar model

Finally, we investigate to what extent the reflection positivity is violated depending on the choice of the parameters $g$ and $M$, although the reflection positivity is everywhere broken. We examine the magnitude of positivity violation in the massive Yang-Mills model which could be regarded as the complementary gauge-scalar model. To estimate the violation of positivity of the Schwinger function $\Delta(t)$, we adopt the ratio $\min_{0<t<\infty}\Delta(t)/\Delta(t=0)$ between the smallest value $\min_{0<t<\infty}\Delta(t)$ of $\Delta(t)$ and the value at the origin $\Delta(t=0)$. Fig. 13 gives the 3D plot and the contour plot of $\min_{0<t<\infty}\Delta(t)/\Delta(t=0)$ on the two-dimensional parameter plane ($\frac{M^{2}}{\mu^{2}}$,$\lambda$)=($\frac{M^{2}}{\mu^{2}}$,$\frac{g^{2}C_{2}(G)}{16\pi^{2}}$). Fig. 14 gives the same plot with larger range of parameters. Note that the left-upper $(\frac{M^{2}}{\mu^{2}}\ll 1,\lambda\gg 1)$ and right-upper $(\frac{M^{2}}{\mu^{2}}\gg 1,\lambda\gg 1)$ regions in the contour plot correspond to the region to be excluded where the Euclidean poles occur. Note that the ratio $\min_{0<t<\infty}\Delta(t)/\Delta(t=0)$ must be negative. However, there are spikes showing positive values in Fig. 14, which are artifacts of our numerical calculations due to the algorithm used for looking for the very small negative value in the very large $t$ as the minimum. These spikes are to be ignored.

If $g^{2}\rightarrow 0$, the theory has no interaction and the propagator approaches the free massive propagator $\mathscr{D}(k)=\frac{1}{k^{2}+M^{2}}$. In this limit, the Schwinger function is positive for any value of $M$ and there is no violation of reflection positivity. For small $g^{2}$ and large $M^{2}/\mu^{2}$, namely, for large $1/g^{2}$ and large $v^{2}\simeq(M^{2}/\mu^{2})/g^{2}$, the Schwinger function exhibits small violation of positivity. This region corresponds to the Higgs-like region in the complementary gauge-scalar model. For large $g^{2}$ and small $M^{2}/\mu^{2}$, namely, for small $1/g^{2}$ and small $v^{2}\simeq(M^{2}/\mu^{2})/g^{2}$, the Schwinger function exhibits large violation of positivity. This region corresponds to the confinement-like region in the complementary gauge-scalar model.

However, there is no phase transition between the positivity violation and restoration. There is just a smooth crossover separating large and small violation of positivity. The massive Yang-Mills model has only one confinement phase. This result is interpreted as the Fradkin-Shenker continuity in the complementary gauge-scalar model from the viewpoint of the gauge-invariant extension from the massive Yang-Mills model to the gauge-invariant complementary gauge-scalar model explained in section II.

V

Complex analysis of the gluon propagator

In the previous section we have investigated the propagator in the Euclidean region. We have shown the violation of reflection positivity in the massive Yang-Mills model. However, this result is obtained only in the numerical way. In this section, we study the propagator on the complex plane of the squared momentum $k^{2}$, which follows from the analytic continuation of the propagator from the Euclidean region to the Minkowski region. We find that the violation of the reflection positivity in the Euclidean region is understood from the existence of a pair of complex conjugate poles and the discontinuity across the branch cut yielding the negative spectral function represented by the generalized spectral representation of the gluon propagator. As a consequence of the complex structure, we give an analytical proof that the reflection positivity is always violated for any choice of the parameters $M$ and $g$ in the massive Yang-Mills model to one-loop order.

V.1 Spectral representation of a propagator

It is well-known that a propagator $\mathscr{D}(k^{2})$ in the Minkowski region $k^{2}>0$ (for the time-like momentum $k$) has the spectral representation of the Källén–Lehmann form under assumptions of the general principles of the QFT such as the spectral condition, the Poincaré invariance and the completeness of the state space spectral_repr_UKKL : The full propagator $\mathscr{D}(k^{2})$ of the field $\phi$ is written as the weighted sum of the free propagator,

[TABLE]

with the weight function $\rho(\sigma^{2})$ called the spectral function being obtained from the state sum

[TABLE]

where $d$ is the space dimension, $D$ is the spacetime dimension, the sum is over all the intermediate states with the total momentum $P_{n}$, and $\theta(k_{0})$ is a step function ensuring the positivity $k_{0}\geq 0$. The spectral function $\rho$ has contributions from a stable single-particle state with physical mass $m_{P}$ (pole mass) and intermediate many-particle states $|p_{1},...,p_{n}\rangle$ with a continuous spectrum, such as two-particle states, three-particle states, and so on,

[TABLE]

Then the spectral representation is written as the sum of the contributions from the real pole $k^{2}=m_{P}^{2}$ and the branch cut

[TABLE]

This spectral representation can be extended to the complex momentum $k^{2}\in\mathbb{C}$. See the left panel of Fig. 15. A propagator $\mathscr{D}(k^{2})$ as a complex function of the complex variable $z=k^{2}\in\mathbb{C}$ has the spectral representation with the spectral function $\rho$,

[TABLE]

This representation (115) is applied to an arbitrary $k^{2}$ in the complex plane except for the singularities located on the positive real axis $[s_{\rm min},\infty)$. The spectral function $\rho$ (116) known as the dispersion relation is obtained from the discontinuity across the branch cut, $\mathscr{D}(z+i\epsilon)-\mathscr{D}(z-i\epsilon)=\mathscr{D}(z+i\epsilon)-\mathscr{D}(z+i\epsilon)^{*}=2i\operatorname{Im}\mathscr{D}(z+i\epsilon).$ It is explicitly checked that the two definitions of the spectral functions (113) and (116) agree with each other once the theory is specified. The representation (115) is obtained under the following assumptions Kondo03 :

$\mathscr{D}(z)$ is holomorphic except singularities on the positive real axis. 2. 2.

$\mathscr{D}(z)\rightarrow 0$ as $|z|\rightarrow\infty$. 3. 3.

$\mathscr{D}(z)$ is real on the negative real axis.

This is indeed the case of the quantum Yang-Mills theory, see e.g., HK18 .

The spectral representation has a straightforward generalization in the presence of complex simple poles, see e.g., Siringo17a ; HK18 . Suppose that the propagator has simple complex poles at $z=z_{\ell}$ $(\ell=1,\cdots,n)$. See the right panel of Fig. 15. Then the propagator $\mathscr{D}(k^{2})$ has the generalized spectral representation,

[TABLE]

where $\gamma_{\ell}$ is a small contour circulating clockwise around the pole at $z_{\ell}$. Here we have separated the propagator $\mathscr{D}$ into the contribution from the complex poles $\mathscr{D}_{p}$ and that from the branch cut $\mathscr{D}_{c}$. This is derived from the following assumptions HK18 :

$\mathscr{D}(z)$ is holomorphic except singularities on the positive real axis and a finite number of simple poles. 2. 2.

$\mathscr{D}(z)\rightarrow 0$ as $|z|\rightarrow\infty$. 3. 3.

$\mathscr{D}(z)$ is real on the negative real axis.

Note that the poles must appear as real poles or pairs of complex conjugate poles as a consequence of the Schwarz reflection principle $\mathscr{D}(z^{*})=[\mathscr{D}(z)]^{*}$.

From now on, we focus on a propagator with a pair of complex conjugate simple poles. This is indeed the case for the gluon propagator of the massive Yang-Mills model as will be shown in the next subsection. For a propagator with one pair of complex conjugate simple poles at $k^{2}=v\pm iw$, the generalized spectral representation (117) reduces to

[TABLE]

V.2 Gluon propagator on the complex momentum plane

We first perform the analytic continuation of the propagator $\mathscr{D}$ in the Euclidean region $k^{2}=-k_{E}^{2}<0$ to the entire complex plane $k^{2}\in\mathbb{C}$. Fig. 16 is the plot of the real and imaginary parts of the complex-valued gluon propagator $\mathscr{D}(k^{2})$ on the complex momentum plane $k^{2}\in\mathbb{C}$, at the physical point of the parameters (102) in the massive Yang-Mills model. Note that the gluon propagator $\mathscr{D}(k^{2})$ is real-valued on the negative real axis (space-like momentum) $k^{2}=-k_{E}^{2}<0$, since the imaginary part $\operatorname{Im}\mathscr{D}(k^{2})$ is zero on the negative real axis (space-like momentum). The real part $\operatorname{Re}\mathscr{D}(k^{2})$ on the negative real axis $k^{2}=-k_{E}^{2}<0$ is identical to the Euclidean propagator. We observe that the gluon propagator has a pair of complex conjugate poles and the imaginary part has discontinuities across the branch cut on the positive real axis $\mathscr{D}(k^{2}+i\epsilon)\not=\mathscr{D}(k^{2}-i\epsilon)$ ($k^{2}>0$, $\epsilon\downarrow 0$), while there are no discontinuities on the negative real axis $\mathscr{D}(k^{2}+i\epsilon)=\mathscr{D}(k^{2}-i\epsilon)$ ($k^{2}<0$, $\epsilon\downarrow 0$). Therefore, in discussing the behavior of the propagator on the positive real axis, we must specify which side is used. In what follows we use the limit $\mathscr{D}(k^{2}+i\epsilon)$ ($\epsilon\downarrow 0$).

Next, we focus on the real axis $k^{2}\in\mathbb{R}$ to see the behavior of the complex-valued gluon propagator $\mathscr{D}(k^{2})$ as a function of a real-valued momentum $k^{2}\in\mathbb{R}$. Fig. 17 is the plot of the real and imaginary parts of the complex-valued gluon propagator on the real axis $k^{2}\in\mathbb{R}$ at the physical point of the parameters (102) in the massive Yang-Mills model. On the negative real axis $k^{2}=-k_{E}^{2}<0$ (the Euclidean region), we find that the real part $\operatorname{Re}\mathscr{D}(k^{2})$ is always positive, and the imaginary part $\operatorname{Im}\mathscr{D}(k^{2})$ is identically zero. On the positive real axis $k^{2}>0$ (the Minkowski region), $\operatorname{Re}\mathscr{D}(k^{2})$ changes the sign such that it is positive for small $k^{2}$, and negative for large $k^{2}$, which implies the existence of (at least one) zeros of $\operatorname{Re}\mathscr{D}(k^{2})$ in the Minkowski region $k^{2}>0$. The scaled imaginary part $\operatorname{Im}\mathscr{D}(k^{2}+i\epsilon)/\pi$ is identical to the spectral function $\rho(k^{2})$. Therefore, the spectral function is identically zero in the Euclidean region,

[TABLE]

However, it is non-trivial in the Minkowski region. It is remarkable that the spectral function is always negative,

[TABLE]

in the massive Yang-Mills model to one-loop order.

For a given propagator $\mathscr{D}(k^{2})$, we can decompose it into the contribution from the branch cut $\mathscr{D}_{c}(k^{2})$ and that from the poles $\mathscr{D}_{p}(k^{2})$. Fig. 18 gives this decomposition of the gluon propagator for the Euclidean momentum $\mathscr{D}(k^{2})=\mathscr{D}_{p}(k^{2})+\mathscr{D}_{c}(k^{2})$ for $k^{2}<0$.

According to the separation of the propagator, the Schwinger function is also separated into the two parts: the continuous cut part $\Delta_{c}(t)$ coming from the spectral function and the pole part $\Delta_{p}(t)$ coming from the pole part $\mathscr{D}_{p}$ of the propagator $\mathscr{D}$,

[TABLE]

Especially, the cut part $\Delta_{c}(t)$ is directly written as an integral of the spectral function as follows

[TABLE]

The same procedure is also applied to the Schwinger function. Fig. 19 shows the respective ratio $\Delta_{c,p}(t)/\Delta(t)$ of the pole or cut part $\Delta_{c,p}(t)$ to the total Schwinger function $\Delta(t)$. Using the already known spectral function $\rho(k^{2})$ calculated according to $\rho(k^{2})=\operatorname{Im}\mathscr{D}(k^{2})/\pi$, the cut part $\Delta_{c}(t)$ of the Schwinger function is obtained by integrating $\rho(k^{2})$ according to (124). Then the pole part $\Delta_{p}(t)$ of the Schwinger function is obtained as the difference $\Delta_{p}(t)=\Delta(t)-\Delta_{c}(t)$ from the total Schwinger function $\Delta(t)$. Note that the ratio can become divergent at a zero $t_{0}$ of the Schwinger function $\Delta(t_{0})=0$, which should be ignored as an artifact of this procedure, see also the caption of Fig. 19.

V.3 A pair of complex conjugate poles and Gribov-Stingl form

If the propagator has no complex poles besides the singularities on the real positive axis, the complex pole part vanishes $\mathscr{D}_{p}(k^{2})=0$ and the Euclidean gluon propagator obeys the usual spectral representation

[TABLE]

Then the Schwinger function is calculated from the cut part alone

[TABLE]

In this case, we find that the positivity of the spectral function $\rho$ implies the positivity of the Schwinger function $\Delta$

[TABLE]

which implies that non-positivity of the Schwinger function $\Delta$ yields non-positivity of the spectral function $\rho$, namely, violation of reflection positivity,

[TABLE]

Thus, when the propagator has no singularities other than the positive real axis, the positivity of the spectral function is directly related to the positivity of the Schwinger function, or the reflection positivity. The violation of reflection positivity can be seen as the non-positivity of the spectral function. However, this is not the case for the Yang-Mills theory, as demonstrated in the massive Yang-Mills model shortly.

Suppose that the propagator has a pair of complex conjugate poles at $k^{2}=v\pm iw$ $(v,w\in\mathbb{R},w>0)$ with the respective residues $Z,Z^{*}\in\mathbb{C}$. Then the pole part of the propagator in the Euclidean region is represented as

[TABLE]

This pole part of the propagator agrees with the *Gribov-Stingl form555 If we apply the definition of the spectral function $\rho$ given in eq.(118) to the Gribov-Stingl propagator (130), we obtain the identically vanishing spectral function, since the Gribov-Stingl propagator does not have the branch cut on the real $k^{2}$ axis across which there is a discontinuity: $\rho(k^{2}):=\frac{1}{\pi}\operatorname{Im}\mathscr{D}(k^{2}+i\epsilon)=\frac{1}{2i\pi}[\mathscr{D}(k^{2}+i\epsilon)-\mathscr{D}(k^{2}-i\epsilon)]$ for $k^{2}\in\mathbb{R}$. Therefore, the Gribov-Stingl propagator has only the pole part and does not have the continuous cut part. * Stingl86 with real parameters $c_{0},c_{1},c_{2},d_{0},d_{1}\in\mathbb{R}$,

[TABLE]

Note that all the coefficients $c_{0},c_{1},c_{2},d_{0},d_{1}$ are not independent. The Gribov-Stingl form actually has four independent parameters, since one of them is eliminated by the rescaling. This number of independent parameters agrees with that of the pole part of the propagator with a pair of complex conjugate poles characterized by the four parameters $v,w,\operatorname{Re}(Z),\operatorname{Im}(Z)$.

For instance, the correspondence between two sets of parameters is given as

[TABLE]

which has the inverse relation

[TABLE]

For $v+iw$ to be a complex number (namely, $w$ to be a real number), the parameters of the Gribov-Stingl form must satisfy the restriction

[TABLE]

Assuming this condition, we can obtain the closed form for the pole part of the Schwinger function

[TABLE]

Indeed, the contribution from one of the poles is exactly evaluated as

[TABLE]

where we have defined

[TABLE]

where $\alpha$ must be located on the upper half plane of the complex $k_{E}$ plane. Therefore, the pole part of the Schwinger function coming from a pair of complex conjugate poles is exactly obtained as Kondo11

[TABLE]

At $t=0$, $\Delta_{p}(0)$ has the value,

[TABLE]

We find that $\Delta_{p}(t)$ is oscillating between positive and negative values, although the absolute value $|\Delta_{p}(t)|$ becomes smaller for larger $t>0$. 666 The pole part of the Schwinger function becomes positive only when the poles become real ones $w=0$ (or $\varphi=0$)

$\displaystyle\Delta_{p}^{\rm GS}(t)=\frac{|\operatorname{Re}(Z)|}{\sqrt{|v|}}\exp[-t\sqrt{|v|}].$

(139)

On the other hand, the cut part $\Delta_{c}(t)$ of the Schwinger function is estimated using the integral representation (124):

[TABLE]

This representation is an exact relation between the spectral function and the cut part of the Schwinger function, which holds irrespective of the existence or non-existence of complex poles.

At least to one-loop order in the massive Yang-Mills model, the spectral function $\rho(\sigma^{2})$ takes the negative value $\rho(\sigma^{2})<0$ for all $\sigma^{2}>0$:

[TABLE]

as demonstrated numerically in Fig. 17 (at the physical point) in this paper and shown analytically for any value of the parameters $g$ and $M$ in HK18 . According to (140), therefore, the cut part $\Delta_{c}(t)$ of the Schwinger function takes negative value for any value of $t$,

[TABLE]

although $\Delta_{c}(t)$ takes smaller and smaller negative value for larger and larger $t>0$.

Moreover, the negative spectral function (141) yields the existence of one pair of complex conjugate poles or two real poles in the Euclidean region as shown in HK18 . According to (137), therefore, the pole part $\Delta_{p}(t)$ of the Schwinger function due to a pair of complex conjugate poles takes negative value for a certain value of $t$:

[TABLE]

Thus, the Schwinger function $\Delta(t)$ obtained as a sum of two parts, $\Delta(t)=\Delta_{p}(t)+\Delta_{c}(t)$ has necessarily negative value at a certain value of $t$,

[TABLE]

Thus we complete the analytical proof that the reflection positivity is always violated irrespective of the choice of the parameters $g$ and $M$ in the massive Yang-Mills model to one-loop order.

In particular, the propagator of the Gribov type is a special case corresponding to $c_{1}=0$ and $d_{0}=0$

[TABLE]

which has a pair of pure imaginary poles

[TABLE]

with the real-valued residue

[TABLE]

The pole part of the Schwinger function for the propagator of the Gribov type (145) is given by

[TABLE]

where we have used $\varphi=\frac{\pi}{4}$ and $r^{1/2}=\left(c_{0}/c_{2}\right)^{1/4}=\sqrt{|w|}$.

V.4 Fitting of the pole part at the physical point

The general analysis given in the above can be substantiated by choosing the physical point for the parameters of the massive Yang-Mills model. We have obtained the pole part of the gluon propagator, as demonstrated in Fig. 18. This pole part $\mathscr{D}_{p}(k^{2})$ is obtained by the difference $\mathscr{D}_{p}(k^{2})=\mathscr{D}(k^{2})-\mathscr{D}_{c}(k^{2})$ from the propagator $\mathscr{D}(k^{2})$ once the cut part $\mathscr{D}_{c}(k^{2})$ of the propagator is specified according to (120) from the imaginary part of the propagator on the positive real axis through the spectral function $\rho(k^{2})$. We have fitted the resulting pole part $\mathscr{D}_{p}(k^{2})$ of the gluon propagator to the Gribov-Stingl form. Fig. 20 is the result of fitting of the pole part of the gluon propagator to the Gribov-Stingl form at the physical point of the parameters. Here we have introduced dimensionless versions of the parameters $\hat{c_{0}},\hat{c_{1}},\hat{c_{2}},\hat{d_{0}},\hat{d_{1}}$ for $c_{0},c_{1},c_{2},d_{0},d_{1}$ and squared momentum $s$ for $k^{2}$, which are scaled by appropriate powers of the gluon mass $M$ to make the dimensionless Gribov-Stingl form

[TABLE]

The fitting parameters are determined as Watanabe19

[TABLE]

which is subject to the restriction (133)

[TABLE]

This result is translated into the complex pole and the residue of the gluon propagator

[TABLE]

It is ensured that this data reproduces the location of the poles given in Fig. 16. We find that the fitting errors are very small and to good accuracy the pole part of the gluon propagator is identical to the Gribov-Stingl form. This result strongly suggests that the pole part of the gluon propagator indeed stems from a pair of complex conjugate poles. For the other argument for understanding the Gribov-Stingl form, see Kondo11 .

The pole part of the Schwinger function can be calculated according to (137) once all the parameters of the Gribov-Stingl form are determined. If our analysis of the complex structure of the propagator is correct, the result should agree with the pole part of the Schwinger function given in Fig. 19. In fact, Fig. 21 shows excellent agreement between pole part of the Schwinger function obtained from the spectral function and the expression (137) with the parameters obtained through the fit of the gluon propagator to the Gribov-Stingl form. This result supports the validity of our arguments.

Thus we have shown that the gluon propagator consists of the pole part due to a pair of complex conjugate poles and the cut part due to the branch cut on the positive real axis, in agreement with the generalized spectral representation (120). This is also the case for the associated Schwinger function. In this way we can conclude that the reflection positivity is violated in the massive Yang-Mills model at the physical point. 777 Notice that the existence of complex poles in the momentum representation of the two-point function does not necessarily violate spacelike commutativity. For instance, it was shown by Nakanishi Nakanishi71 that the existence of complex pole is compatible with spacelike commutativity in a complex scalar field theory with indefinite metric. Notice that this theory is manifestly Lorentz covariant in a finite duration of time, against the title of the paper.

V.5 Parameter dependence other than the physical point

We investigate the gluon propagator and the associated Schwinger function at choices of the parameters other than the physical point.

V.5.1 Smaller gauge coupling

For smaller value of the coupling constant $\lambda:=Ng^{2}/(4\pi)^{2}=0.1$ with the physical value for $M$, we obtain the gluon propagator in Fig. 22 and the Schwinger function in Fig. 23. The cut part is relatively large and has the opposite sign to the pole part to cause cancellation. The fall-off of both parts of the Schwinger function is slow for smaller value of the coupling constant. Therefore, the large $t$ behavior of the Schwinger function must be investigated to see the violation of positivity due to the cancellation between two parts for this choice of parameters other than the physical point which can be identified with an effective model of the pure Yang-Mills theory.

V.5.2 Smaller gluon mass

For smaller value of the gluon mass at $M^{2}/\mu^{2}=0.02$, we obtain the gluon propagator in Fig. 24 and the Schwinger function in Fig. 25. The cut part is relatively rather small. Therefore, violation of reflection positivity largely comes from the pole part.

We find that the cut parts $\mathscr{D}_{c}$ and $\Delta_{c}$ of the gluon propagator and the Schwinger function are always negative. This result reconfirms in a numerical way that the spectral function is negative irrespective of the choice of the parameters. Thus, the reflection positivity is always violated for any choice of the parameters.

VI

Conclusion and discussion

We have examined the mass-deformed Yang-Mills theory or the massive Yang-Mills model in the covariant Landau gauge with two parameters, the coupling constant $g$ and the mass parameter $M$, in order to reproduce the confining decoupling solution of the pure Yang-Mills theory. By choosing appropriate values for a set of two parameters $g$ and $M$, we have shown that the massive Yang-Mills model well reproduces simultaneously the gluon and ghost propagators of the decoupling solution obtained by the numerical simulations on the lattice (at least) in the low-momentum region. Such a choice of the parameters is called the physical point for the Yang-Mills theory.

Then we have shown that the reflection positivity is violated in the massive Yang-Mills model at the physical point of the parameters by observing the negativity of the Schwinger function which is obtained in a numerical way as the Fourier transform of the gluon propagator.

The violation of reflection positivity was also confirmed by examining the complex structure of the complex-valued gluon propagator obtained by performing the analytic continuation of the Euclidean propagator to the entire complex squared momentum plane. We have verified that the violation of reflection positivity in the Euclidean region detected by the Schwinger function associated with the Euclidean gluon propagator is a consequence of the complex structure of the complex-valued gluon propagator: (i) the negativity of the spectral function obtained from the discontinuity of the gluon propagator across the branch cut on the positive real axis on the complex squared momentum plane, (ii) the existence of a pair of complex conjugate poles in the gluon propagator. At the physical point, the contribution from the cut part to the gluon propagator in the Euclidean region is relatively small compared with that from the pole part. Therefore, the propagator in the Euclidean region is well approximated by the contribution from a pair of complex conjugate poles in the complex region, which implies that the propagator in the Euclidean region is well described by the Gribov-Stingl form, in agreement with the lattice result DOS18 . The violation of reflection positivity is regarded as a necessary condition for gluon confinement. Therefore, our results of reflection positivity violation at the physical point of the massive Yang-Mills model support strongly gluon confinement in the Yang-Mills theory.

We have regarded the massive Yang-Mills model at the physical point as the low-energy effective model of the pure Yang-Mills theory. However, the massive Yang-Mills model with the parameters $g$ and $M$ other than the physical point has another meaning. We have discussed that the massive Yang-Mills model in the covariant Landau gauge has the gauge-invariant extension, which is identified with the complementary gauge-scalar model with a radially fixed fundamental scalar field which is subject to an appropriate reduction condition. In other words, the gauge-scalar model with a radially fixed fundamental scalar field subject to the reduction condition can be gauge-fixed to becomes the massive Yang-Mills model in the covariant Landau gauge. The gauge-invariant extension of a non-gauge theory is performed through the gauge-independent description of the BEH mechanism Kondo18 without relying on the spontaneous symmetry breaking which was first proposed for the adjoint scalar field Kondo16 .

Therefore, the Yang-Mills theory in the confinement phase characterized by the decoupling solution is identified with the massive Yang-Mills model with the physical point of the parameters $g$ and $M$. This physical point of the massive Yang-Mills model represents a point in the parameter space of the complementary gauge-scalar model obtained as a gauge-invariant extension of the massive Yang-Mills model. Thus, the violation of reflection positivity in the massive Yang-Mills model for any value of the parameters $g$ and $M$ is consistent with the Fradkin-Shenker continuity in the sense that the massive Yang-Mills model describes both confinement-like and Higgs-like regions in the single confinement phase of the complementary gauge-scalar model. Our result seems to be consistent with the other approaches CPRW18 ; RSTW17 .

Let us make comments on the sum rule for the spectral function called the superconvergence relation OZ80 ,

[TABLE]

It is obvious that this sum rule cannot be satisfied for the negative spectral function. In HK18 , remarkably, the generalized sum rule for the spectral function has been derived in the presence of a pair of complex conjugate poles

[TABLE]

provided that the propagator has the asymptotic behavior

[TABLE]

in the region far from the origin of the complex $k^{2}$ plane. In fact, it is shown that the gluon propagator fulfills this condition in Yang-Mills theories in the Landau gauge due to the asymptotic freedom and the negativity of the anomalous dimension OZ80 . The assumption (155) is enough to obtain

[TABLE]

The real part of (156) leads to a generalized sum rule for the spectral function (154) by taking into account the relation (118):

[TABLE]

while the imaginary part leads to another relation,

[TABLE]

Therefore, the usual superconvergence relation (153) does not hold unless the residue of the complex pole is pure imaginary $\operatorname{Re}Z=0$. The preliminary results for the massive Yang-Mills model to one-loop order at the physical point are Watanabe19

[TABLE]

which leads to

[TABLE]

It will be interesting to examine whether the generalized sum rule holds or not, and to what extent it is satisfied beyond one-loop level when the gluon propagator has a pair of complex conjugate poles.

Moreover, it is desirable to extend the results obtained in this paper to a finite temperature to see whether or not the reflection positivity violated in the low-temperature confinement phase is recovered in the high-temperature deconfinement phase even in the pure Yang-Mills theory. Then we can ask whether or not the transition detected by the positivity violation/restoration agrees with the confinement/deconfinement transition detected by the Polyakov loop average. It is also interesting to examine how the relevant complex structure changes depending on the temperature. These issues will be discussed in subsequent papers.

Finally, we give some comments on the obstructions stemming from the presence of complex poles in the gluon propagator to the formal field theoretic issues such as loss of locality. In local QFT, it is recognized that any 2-point correlation function is an analytic function in the cut complex $p^{2}$(squared momentum)-plane with singularities along the time-like (positive) real axis only. The assumptions to establish this analytic property are AQFT : Lorentz covariance (covariance under space-time translations), the spectrum condition, local (or space-like) commutativity, and uniqueness and cyclicity of the vacuum. For any other singularity structure of 2-point correlations, at least one of these assumptions must be violated.

In order to consider this issue in the Yang-Mills theory, we can take into account the observation OZ80 ; Oehme95 that the correlation functions of the Yang-Mills field vanish in the limit $p^{2}\to\infty$ in all directions of the complex $p^{2}$-plane due to the ultraviolet asymptotic freedom of the Yang-Mills theory. However, non-trivial entire functions with that property do not exist. Therefore, they must have singularities somewhere. Usually, the singularities are supposed to exist on the positive real $p^{2}$ axis. However, this does not deny the existence of complex conjugate poles discussed in this paper. See also Oehme95 and section 2.5 and 5.4 of AS01 for the review.

Indeed, the existence of complex poles does not only play a crucial role in the violation of the reflection positivity but also breaks the spectral representation, which is a fundamental implication of local field theories. In the standard point of view, e.g., from the Jost-Lehmann-Dyson (JLD) representation JL57 ; Dyson58 , complex poles must violate the local spacelike commutativity.

One might claim that the non-locality of the Yang-Mills theory in a gauge-fixed picture is rather “natural” due to the Gribov-Singer obstruction, see Gribov78 ; Singer78 ; Maas13 and Zwanziger89 ; Baulieu-etal . The problem of locality is discussed in Stingl96 ; HKRSW90 , in which they assert that complex poles describe short-lived excitations, and the locality is broken in the level of propagators, but the corresponding $S$-matrix remains causal. However, their way of reconstructing the Minkowskian propagator from the Euclidean propagator is questionable. By a straightforward reconstruction of the analytic continuation in the complex time plane, the resulted Minkowskian propagator differs from the previous one Stingl96 and predicts the Lee-Wick type propagator LW69 . In fact, without the positive definiteness of the state space and the spectral condition, which are not guaranteed for confined degrees of freedom, complex spectra can appear. The Lee-Wick type theories can yield complex poles without the loss of the spacelike commutativity, see, e.g. Nakanishi71b as the simplest example for the propagator with complex poles. Notice that the complex spectra enable a theory to evade the restriction of the axiomatic or analytic theorems like the JLD representation because it deviates from the framework of tempered distributions. In this scenario, complex poles are not an indication of the non-locality but just a reflection of unphysical degrees of freedom, such as timelike photons. This issue will be further discussed elsewhere.

Another important issue to be addressed in the presence of the complex poles in the gluon propagator is to answer the question how the correlation function of color singlet composite operators can have the real poles, since such composite states must be observed. This issue was investigated in Zwanziger89 and Baulieu-etal . It is further argued in Stingl96 that complex singularities with time-like real part might be acceptable for the propagators of unphysical colored fields, by the reason that such singularities might conspire to cancel with singularities or zeros in other unphysical correlation functions so as to be absent from physical amplitudes. This will give rise to an infinite hierarchy of constraints on such unphysical singularities in arbitrary high $n$-point functions. An example of such compensating singularities are those in the non-perturbative expansion scheme, see Stingl96 ; HKRSW90 and sections 2.5 and 5.4 of AS01 . More serious discussion on the related issues will be given in future works.

Acknowledgements

The authors would like to thank Professor Oliveira for providing us with the data of numerical simulation for gluon and ghost propagators presented in this paper. K.-I. K. thanks Matthieu Tissier, Nicolas Wschebor, Julien Serreau, and Urko Reinosa for the discussion on positivity violation. This work was supported by Grant-in-Aid for Scientific Research, JSPS KAKENHI Grant Number (C) No.19K03840 and No.15K05042. R. M. was supported by Grant-in-Aid for JSPS Research Fellow Grant Number 17J04780.

Appendix A Recursive construction of the gauge-invariant transverse field

In this section, we derive the recurrence relation for obtaining the power-series solution of the reduction condition which is equivalent to the transverse condition,

[TABLE]

First, we expand $h[\mathscr{A}]$ into the power series in the gauge field $\mathscr{A}$,

[TABLE]

where the superscript $(k)$ denotes $k$-th term in the power series. Note that $(h^{-1})^{(k)}=(h^{(k)})^{\dagger}$ because $h^{-1}=h^{\dagger}$. Then the $k$-th term of $\mathscr{A}_{\mu}^{h}=h\mathscr{A}_{\mu}h^{\dagger}-ig^{-1}\partial_{\mu}hh^{\dagger}$ is given by

[TABLE]

By substituting (164) into 161 we obtain the recurrence relation as

[TABLE]

By using this recurrence relation, we can derive several features of $\mathscr{A}^{h}_{\mu}$. First, we show that $\mathscr{A}_{\mu}^{h}$ can be written as

[TABLE]

This is shown as follows. Indeed, by substituting the recurrence relation 165 into 164, we obtain

[TABLE]

where we have defined

[TABLE]

By taking into account the fact that $F^{(k)}_{\mu}$ is not Hermitian, the $k$-th term of $\Psi_{\mu}$ in 166 can be written as

[TABLE]

Next, we show that $A_{\mu}^{h}$ is gauge invariant by using 165. Now let the gauge transformation of $\mathscr{A}$ by $V$ be

[TABLE]

We observe that $\mathscr{A}_{\mu}^{h}$ is indeed gauge invariant if the gauge transformation of $h$ obeys

[TABLE]

whose infinitesimal form for $V=e^{-ig\Lambda}$ is given by

[TABLE]

In the following, we show that 172 holds order by order.

Note that because $\delta\mathscr{A}_{\mu}=\partial_{\mu}\Lambda-ig[\mathscr{A}_{\mu},\Lambda]$, $\delta h^{(k)}$ contains $k$-th term and $(k-1)$-th term of power series of $\delta h$. For this reason, we separate $\delta h^{(k)}$ into two parts as

[TABLE]

where $\delta h^{(k)}_{-}$ has the order of $k-1$ and $\delta h^{(k)}_{=}$ has the order of $k$. Therefore, the $k$-th order term of 172 is written as

[TABLE]

This relation is shown to hold with mathematical induction as follows. For $k=0$, it follows from 165 that

[TABLE]

By taking the variation under the gauge transformation, we obtain

[TABLE]

Thus we have shown 174 holds for $k=0$. Next, suppose that 174 holds for $k-1$. Then we proceed to show that it holds for $k$. By taking variation of 165 for $k+1$ under the gauge transformation we obtain

[TABLE]

By taking the variation of 165 for $k$ we obtain

[TABLE]

By summing up the underlined part, double-underlined part, wavy-lined part and broken-lined part respectively, we obtain

[TABLE]

where we have used the assumption of induction. The sum of 179 and 180 cancels the third term in the parentheses of 178. The sum of 181 and 182 cancels the first term in the parentheses of 177. Therefore 174 is satisfied for $k$. Thus we have shown that $\mathscr{A}^{h}$ is invariant under a gauge transformation.

Appendix B

Solving the reduction condition in another way

By using the massive vector field mode $\mathscr{W}_{\mu}$ (17), the reduction condition reads

[TABLE]

For the scalar field $\hat{\Theta}(x)$, we introduce the Lie algebra $\mathscr{G}$-valued field $\theta(x)$ as

[TABLE]

In the following, we solve the reduction condition by expressing the scalar field $\hat{\Theta}(x)$ as a power series in the gauge field $\mathscr{A}_{\mu}$. The Lie algebra form of the pure gauge reads

[TABLE]

The more general expression is given as

[TABLE]

By substituting this result into the reduction condition $\chi=0$, we have

[TABLE]

This is recast into

[TABLE]

which yields

[TABLE]

Substituting recursively for $\theta$, we obtain a power series,

[TABLE]

The massive vector field mode $\mathscr{W}_{\mu}$ (17) is written as

[TABLE]

where we have defined the transverse field $\mathscr{A}_{\mu}^{T}$ in the lowest order term linear in $\mathscr{A}$ as

[TABLE]

Notice that $\mathscr{W}_{\mu}$ agrees with $\mathscr{A}_{\mu}$ in the Landau gauge $\partial\cdot\mathscr{A}=0$.

Thus, by substituting $\mathscr{W}$ of (191) into (16), the term $S_{\rm kin}$ reads

[TABLE]

By performing integration by parts and taking into account the transversality $\partial_{\mu}\mathscr{A}_{\mu}^{T}=0$, the action $S_{\rm kin}^{*}$ takes the form,

[TABLE]

This indeed agrees with the expression (37).