Extracting topological information from momentum space propagators

TL;DR

This paper introduces a new gauge-invariant topological invariant based on momentum space propagators, capable of classifying complex poles and potentially detecting phenomena like chiral symmetry breaking and deconfinement.

Contribution

A novel topological invariant sensitive to propagator structures is proposed, extending to complex poles and applicable in non-perturbative quantum field theory contexts.

Findings

Invariant can classify propagators with complex poles

Potential to detect chiral symmetry breaking

May identify deconfinement transitions

Abstract

A new topological invariant quantity, sensitive to the analytic structure of both fermionic and bosonic propagators, is proposed. The gauge invariance of our construct is guaranteed for at least small gauge transformations. A generalization compatible with the presence of complex poles is introduced and applied to the classification of propagators typically emerging from non-perturbative considerations. We present partial evidence that the topological number can be used to detect chiral symmetry breaking or deconfinement.

Equations

$${\mathcal{N}}_3=N\mathrm{Tr}{\int }_{\Sigma }KG\mathrm{d}{G}^{-1}\wedge G\mathrm{d}{G}^{-1}\wedge G\mathrm{d}{G}^{-1},$$

$$\delta {\mathcal{N}}_3=3N{\int }_{\Sigma }\text{mbox}Tr\left[K\delta [G\mathrm{d}{G}^{-1}]\wedge G\mathrm{d}{G}^{-1}\wedge G\mathrm{d}{G}^{-1}\right],$$

$$3N{\int }_{\partial \Sigma }\text{mbox}Tr\left[K(G\delta G^{-1})G\mathrm{d}{G}^{-1}\wedge G\mathrm{d}{G}^{-1}\right]=0,$$

$$G^{AB}(p)=\mathcal{G}(p){\delta }^{AB}=\frac{Z(p^2)}{i/p+\mathcal{M}(p^2)}{\delta }^{AB},$$

$$\underset{|\vec{p}|\to +\infty }{\lim }\frac{\delta \mathcal{M}}{|\vec{p}|{\left(1+\frac{{\mathcal{M}}^2}{p^2}\right)}^2}=0.$$

$${\mathcal{N}}_3$$

$${\epsilon }_{ijk}\mathrm{Tr}\left[{\gamma }^5{\gamma }^0{\gamma }_i{\gamma }_j{\gamma }_k\right]$$

$$3p^l{p}^k{\epsilon }_{ijk}\mathrm{Tr}\left[{\gamma }^5{\gamma }^0{\gamma }_l{\gamma }^i{\gamma }^j\right]$$

$${\mathcal{N}}_3$$

$$\theta (p)=4\left[\frac{\mathcal{M}(p)/p}{1+{\left(\mathcal{M}(p)/p\right)}^2}+\mathrm{Arctan}\left(\mathcal{M}(p)/p\right)\right].$$

$$\mathcal{M}(p^2)=\frac{M^3}{p^2+m^2}+\mu .$$

$${\mathcal{N}}_{\Gamma }=\frac{1}{2\pi \mathrm{i}}{\oint }_{\Gamma }\mathrm{d}z4\mathrm{i}\sqrt{z}\frac{\left[\mathcal{M}(z)-2z{\partial }_z\mathcal{M}(z)\right]}{{\left(z+{\mathcal{M}}^2(z)\right)}^2},$$

$$D(p^2)=\frac{Z(p^2)}{p^2+{\mathcal{M}}^2(p^2)}$$

$$C(\alpha )=R{\mathrm{e}}^{i\alpha }\text{with}\pi -\frac{\epsilon }{R}\le \alpha \le \frac{\epsilon }{R}$$

$$L_{-}(t)=t+i\epsilon \text{with}-R\le t<0$$

$$L_{+}(t)=t+i\epsilon \text{with}0\le t\le R.$$

$${\mathcal{N}}_{\Gamma }=\frac{1}{2\pi \mathrm{i}}{\oint }_{\Gamma }dz\frac{f^{\prime }(z)}{f(z)}=-\frac{1}{2\pi }{\int }_{\theta \left(f(\Gamma )\right)}d\theta (z)$$

$$f(z)={\mathrm{e}}^{-\mathrm{i}\theta (z)},\theta (z)=2\left[\frac{\sqrt{z}\mathcal{M}(z)}{z+{\mathcal{M}}^2(z)}+\mathrm{Arctan}\left(\frac{\mathcal{M}(z)}{\sqrt{z}}\right)\right].$$

$$\delta {\mathcal{N}}_{\Gamma }$$

$$-\frac{2\sqrt{z}z}{(z+{\mathcal{M}}^2{)}^2}{\partial }_z\delta \mathcal{M}=\frac{3\sqrt{z}\delta \mathcal{M}}{(z+{\mathcal{M}}^2{)}^2}-\frac{4\sqrt{z}z(1+2\mathcal{M}{\partial }_z\mathcal{M})\delta \mathcal{M}}{(z+{\mathcal{M}}^2{)}^3}-\frac{\partial }{\partial z}\left(\frac{2\sqrt{z}z\delta \mathcal{M}}{(z+{\mathcal{M}}^2{)}^2}\right).$$

$$\delta {\mathcal{N}}_{\Gamma }=-\frac{4}{\pi }{\oint }_{\Gamma }\frac{\partial }{\partial z}\left(\frac{\sqrt{z}z\delta \mathcal{M}}{(z+{\mathcal{M}}^2{)}^2}\right)\mathrm{d}z.$$

$$F(z)=\sum _{N=-m}^{+\infty }{C}_n(z-a{)}^n=\sum _{n=1}^m\frac{C_{-n}}{(z-a{)}^n}+\sum _{n=0}^{+\infty }{C}_n(z-a{)}^n.$$

$$F^{\prime }(z)$$

$$F^{\prime }(z)=\sum _{n=-m-1}^{+\infty }{E}_n(z-a{)}^n,$$

$$\underset{|z|\to \infty }{\lim }\mathcal{M}(z)=\mu \text{and}\underset{|z|\to \infty }{\lim }\delta \mathcal{M}(z)=\delta \mu .$$

$$D(p^2)=\frac{c(1+d{p}^2)}{p^4+2r^2{p}^2+r^4+b^2}=\frac{cd}{p^2+{\mathcal{M}}^2(p^2)}.$$

$${\mathcal{M}}^2(z)=\frac{z(2r^{2}d-1)+(r^4+b^2)d}{1+dz}.$$

$$\underset{|z|\to \infty }{\lim }\mathcal{M}(z)=\sqrt{\frac{2r^{2}d-1}{d}}\text{and}\underset{|z|\to \infty }{\lim }\delta \mathcal{M}(z)=\frac{4r{d}^{\frac{5}{2}}\delta r+{\mathrm{d}}^{\frac{1}{2}}\delta d}{d^2(2r^{2}d-1{)}^{\frac{1}{2}}},$$

$$C(t)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }dp{e}^{ipt}D(p^2),$$

$${\mathcal{N}}_{\Gamma }=-\frac{1}{2\pi }\underset{R\to \infty }{\lim }\left[{\int }_{C(\alpha )}\mathrm{d}\theta (z)+{\int }_{L_{-}(t)}\mathrm{d}\theta (z)+{\int }_{L_{+}(t)}\mathrm{d}\theta (z)\right].$$

$${\cal N}_{\Gamma}~{}=~{}-\frac{1}{2\pi}\Bigg{[}\int_{\pi-\frac{\varepsilon}{R}}^{\frac{\varepsilon}{R}}\frac{\mathrm{d}}{\mathrm{d}\alpha}\theta(R\mathrm{e}^{i\alpha})\,\mathrm{d}\alpha+\int_{-R}^{0}\frac{\mathrm{d}}{\mathrm{d}t}\theta(t+i\varepsilon)\,\mathrm{d}t+\int^{R}_{0}\frac{\mathrm{d}}{\mathrm{d}t}\theta(t+i\varepsilon)\,\mathrm{d}t\Bigg{]}\;,$$

Full text

I Introduction

Topology has been used since a long time in the study of condensed matter physics and, nowadays, it has become the main theoretical tool for the description and classification of topological insulators (see Qi and Zhang (2011) for a comprehensive review). This is because it is possible to construct topological invariant quantities, i.e. , quantities that are invariant under “smooth deformations” (homeomorphisms). In this manner the space of parameter-dependent results is naturally divided into disjoint sectors, invariant with respect to homeomorphisms, that are characterized by different values for the topological invariants. Also in a particle physics, or more general, (quantum) field theory context, the importance of topology cannot be underestimated, see Rajaraman (1982) for the typical illustrative examples.

From the physical point of view, different sectors correspond to different physical phases. Similarly, momentum space topology (MST) has been applied to the classification of the ground state of relativistic quantum field theories (see So (1985), where MST is first applied to lattice fermions) into universality classes, however without the same resonance as in condensed matter Horava (2005); Volovik (2013, 2010, 2007); Zubkov and Volovik (2012); Zubkov (2018). The MST framework might be a very powerful tool and has also been applied to the investigation of emergent gravity, Klinkhamer and Volovik (2008); Volovik (2009). Indeed relativistic quantum fields share some topological properties with topological materials Horava (2005). The Standard Model (SM) below the Electroweak scale, for example, belongs to the same universality class as of the superfluid 3He and of the three-dimensional topological insulators obeying time-reversal symmetry Volovik (2013). However, up to now, topological invariants constructed in MST have not been applied to bosonic systems. In a sense, this could be interpreted as one of the major drawbacks of MST, which only applies to “half of the world”, for reasons that will become clear in the next sections. On the contrary, it would be desirable to apply MST to bosonic fields, especially considering that the analytic structure in momentum space of bosonic propagators might encode extremely valuable physical information.

Inspired by the power of MST in the classification of the vacuum of quantum field theories, we propose a new topologically invariant object that is not only sensitive to the full analytic structure of fermionic propagators, but that can also be applied to bosonic field propagators. Particularly, we will consider here the case of certain classes of fits to lattice fermion propagators Furui and Nakajima (2006); Dudal et al. (2016); Burgio et al. (2012); Rojas et al. (2013) and gluon propagators Dudal et al. (2010); Cucchieri et al. (2012a); Dudal et al. (2012, 2018), although the proposed tool can be applied to any relativistic fermionic or bosonic propagator. Hence, the present results give rise to a considerable enhancement of the power, as well as applicability, of MST. The paper is organized as follows: in Section II we recall the definition and properties of the momentum space topological invariant. In Section III we extract, departing from the standard definition, a specific representation for the MST invariant that is very useful for a generalization to the bosonic sector. In Section IV, we generalize this definition to make it compatible with the potential presence of complex poles and we apply this latest definition to several concrete propagators in Section V. Finally, in Section VI, our conclusions can be found.

II The momentum space topological invariant

Let us start with the following topological quantity, to our knowledge first introduced by So So (1985), see also Zubkov (2018); Katsura and Koma (2018); Volovik and Zubkov (2017),

[TABLE]

where $K$ stands for the matrix representation of a vacuum symmetry; $G$ is the two-point Green function of the fermion field; the integral is on a three/dimensional hypersurface $\Sigma$, in our case defined by imposing $p_{0}=0$. From now on, we assume the following notation: $p^{2}=p_{i}p^{i}$ and $\not{p}=\gamma_{i}p^{i}$, with $i=1,2$ and $3$ (the spatial indices); finally, $N$ is a normalization factor. Notice that the foregoing implies that $[K,G]=0$. A rather similar, albeit not exactly the same, topological invariant was also studied in the lattice context of Coste and Luscher (1989).

At first glance, the topological invariant defined above can be constructed only for fermion propagators, since the $\gamma$-matrices are indispensable to get a non-trivial result. A direct computation reveals that if one tries to make sense of the above expression for bosonic propagators with internal indices, the result vanishes identically while for fermions the $\gamma$-matrices save the day by generating the completely antisymmetric tensor.

Our main goal now for the following sections is to propose a new, and more general, topological invariant that can also be applied to bosonic fields.

If there would be no interactions between the fermions, the quantity (1) becomes nothing else than the number of massive flavors of Dirac fermions Zubkov and Volovik (2012). Adding interactions in such a way that these do not spoil the condition $[G,K]=0$, the change of Eq. (1) under a continuous deformation $G\to G+\delta G$, where $\delta G$ encodes the variations of the parameters in $G$, can be written as

[TABLE]

at least, if the three-form $\delta[G\mathrm{d}G^{-1}]\wedge G\mathrm{d}G^{-1}\wedge G\mathrm{d}G^{-1}$ is continuous inside $\Sigma$, according to the Leibniz integration rule (Courant and John, 1989, p. 466). Thus, to ensure $\delta\mathcal{N}_{3}=0$, it is necessary to have

[TABLE]

where, as $\Sigma\equiv S^{3}$ in momentum space, $\partial\Sigma$ is the spherical two-dimensional surface with radius $|\vec{p}|$, whilst $|\vec{p}|\to+\infty$. Eq. (2) is the most general condition that any Green function $G$ must fulfill so that $\mathcal{N}_{3}$ is invariant under small (continuous) deformations of $G$. In Section III, we shall derive a more specific condition for a particular, but general enough, class of fermionic propagators. These have been used to fit lattice data for the quark propagator in QCD; see Furui and Nakajima (2006); Dudal et al. (2016); Burgio et al. (2012); Rojas et al. (2013).

A priori, one might question the physical relevance of the topological invariants discussed here as they are explicitly based on gauge variant input, viz. propagators that will depend on a chosen gauge. In practice, this will boil down to a possible dependence on a gauge parameter. Indeed, changing the gauge parameter can always be achieved via consecutive application of small (infinitesimal) gauge variations (or BRST variation, if you wish), which is just a special class of continuous deformations leaving the topological number untouched. We will have nothing to say about large gauge transformations and invariance of ${\cal N}_{3}$ w.r.t. those.

III Topological invariant for fermionic case: old and new results

In this Section we will apply the previous ideas to the case of a fermion propagator in a $SU(N)$ gauge theory, with most general parameterization (assuming Lorentz invariance of course):

[TABLE]

where $AB$ are internal Dirac fermion indices111Usually, fermions belong to the fundamental representation of the internal $SU(N)$ group.. It can be shown that the boundary condition (2) is recovered if Zubkov and Volovik (2012)

[TABLE]

For this kind of fermionic propagators and up to powers of logarithms as dictated by the renormalization group equation, the dynamical mass $\mathcal{M}(p^{2})$ tends to a constant mass $\mu$ (possibly zero) in the UV limit ($p^{2}\rightarrow\infty$), while the renormalization function $Z(p^{2})$ goes to $1$ in such a limit. Therefore, when $p^{2}\rightarrow\infty$ the standard UV fermionic propagator $\frac{\delta^{AB}}{\mathrm{i}\not{p}+\mu}$ is restored, up to logarithmic terms. On the other hand, in the IR limit ($p^{2}\rightarrow 0$) one usually assumes that $\mathcal{M}(p^{2})$ continuously tends to a constant value $\mathcal{M}_{0}$. This can be appreciated from e.g. lattice, functional or fitting approaches, Braun et al. (2016); Aguilar and Papavassiliou (2011); Alkofer et al. (2009); Bashir et al. (2012); Oliveira et al. (2019); Bowman et al. (2005); Burgio et al. (2012); Dudal et al. (2016); Furui and Nakajima (2006).

Assuming the above mentioned asymptotic behavior of $Z(p^{2})$ and $\mathcal{M}(p^{2})$, one can verify directly that propagators of the kind (3) fulfil condition (4).

Then, following Volovik (2010), the matrix element $K$ will be considered as being the one that accounts for the CT-symmetry, representable by $\gamma^{5}\gamma^{0}$ (the four-dimensional Euclidean Dirac matrices). A nice overview of CPT symmetry can be found in Oblak . Notice that $\gamma_{5}\gamma_{0}$ does commute with the general Dirac fermion propagator given at Eq. (3). Then, since we already showed that ${\cal N}_{3}$ is a topological invariant, we can set $Z(p^{2})=1$, as this $Z(p^{2})$ for $p^{2}>0$ is a smooth deformation of $1$. One can again appreciate this from non-perturbative functional or lattice computations; see e.g. Braun et al. (2016); Aguilar and Papavassiliou (2011); Alkofer et al. (2009); Bashir et al. (2012); Oliveira et al. (2019); Bowman et al. (2005); Burgio et al. (2012); Furui and Nakajima (2006). Therefore, after some algebraic manipulations, we can reduce

[TABLE]

to

[TABLE]

by making use of the Euclidean identities

[TABLE]

Let us now observe that Eq. (6) still carries the “fingerprint” of the topological invariance of the original expression (5). Indeed, even after performing the integral over the momentum in order to arrive at Eq. (6), this expression can be interpreted as a one-dimensional topological invariant, associated to the analytic structure of $\mathcal{M}(p^{2})$ itself. A fundamental observation is that one can also forget about how we did arrive at Eq. (6), since such an expression is a perfectly well-defined winding number associated to the dynamical mass $\mathcal{M}(p^{2})$ through the curve $\theta(p)$ defined in Eq. (8) below, providing there is no pathological behavior of the integrand. In other words, if one is presented the expression (6) for the first time, without any previous knowledge of Eq. (5), one is still able to show that this quantity (6) is a well-defined one-dimensional topological invariant. More precisely,

[TABLE]

where we have defined the function $\theta(p)$ such that

[TABLE]

The significant advantage of this observation is that Eq. (6) can be considered as a legitimate invariant in itself for any dynamical (fermionic or bosonic) mass functions (that satisfy boundary conditions).

In the following we present a detailed analysis of the topological classification of a Dirac fermionic system whose dynamical mass reads

[TABLE]

For the obtained values of parameters coming from a lattice fit Furui and Nakajima (2006); Dudal et al. (2016); Rojas et al. (2013), $(M^{3},m^{2},\mu)$, (always keeping that order) all of them are positive and, therefore, the result for the topological invariant is always $\mathcal{N}_{3}=1$, as one can readily verify. However, depending on possible other configurations of these parameters, other values of ${\cal N}_{3}$ can emerge. This “pathological” behavior could imply phase transitions, provided in such a case $\delta\mathcal{N}_{3}=0$ cannot be guaranteed for a small variation of the parameters. This is translated into the fact that there are configurations, characterized by an equation on the space of parameters $f(M^{3},m^{2},\mu)=0$, which are unstable under small variation of the parameters. These surfaces can be realized as phase boundaries Goldenfeld (1992). One concrete example of such phase boundary is shown in Fig. 1, where the two-dimensional surface is characterized by the equation $M^{3}+m^{2}\mu=0$ on the three-dimensional space of parameters222One can check that this condition translates into the Green function having a zero mass pole.. On such surface the topological invariant is $\mathcal{N}_{3}=0$. Notice that this is also the topological invariant’s value in the chirally invariant case of a simple massless Dirac quark propagator. As such, the topological number $\mathcal{N}_{3}$ can discriminate between chirally symmetric and chirally broken phases.

Away from this critical surface, the $\mathcal{N}_{3}$ jumps. To set the mind, working in unspecified mass units, for $(-1,1,1)$ we have $\mathcal{N}_{3}=0$; for $(-1.01,1,1)$, $\mathcal{N}_{3}=-1$ and for $(-0.99,1,1)$, $\mathcal{N}_{3}=+1$. For all of these parameter sets, the propagator poles are located on the negative real axis, whilst the mass function remains positive for $p^{2}>0$. Stated otherwise, all of them could in principle be used to try to fit the non-perturbative quark mass function as found on lattice; see e.g. Bowman et al. (2005); Oliveira and Silva (2019). In future work, it would be interesting to compute $\mathcal{N}_{3}$ directly based on lattice data, after which all kinds of approximations to it, for instance born out of Dyson-Schwinger/functional renormalization group equations mentioned already before, or approaches likes Alkofer et al. (2004); Peláez et al. (2014), can be tested if they belong to the same homotopy class or not.

More general, if the parameters take (critical) values such that the denominator of the integrand in Eq. (6) has singular behavior for positive values of $p^{2}$, there will be a transition $\Delta\mathcal{N}_{3}\neq 0$ of the topological invariant for a small variation of the critical parameters. As $\mathcal{N}_{3}$ is an integer number, this transition must come with another integer value for the topological number.

IV A generalized topological invariant

In this Section, we propose to extend the momentum space topological invariant (6) in order to make it sensitive to complex poles of the mass function $\mathcal{M}$, of which we always assume to know the analytic structure333Evidently, this is a big assumption. For example, a Monte Carlo based lattice computation will never give direct access to the full analytic structure in the complex momentum squared plane. But, using dedicated inversion schemes can shed light on this anyhow, see e.g. Fischer et al. (2018), in some cases also (numerical) solutions over the complex momentum squared plane of the Dyson-Schwinger equations Maris (1995)Fischer et al. (2018) can be obtained. . Such an extension is a topological invariant, applicable either to fermionic or to bosonic systems.

This quantity is defined by the following integral in the complex plane:

[TABLE]

where $\Gamma$ is the contour defined to enclose all the possible complex poles of the integrand (10) lying in the upper half of the complex plane (see Fig. 2), while avoiding every possible real pole, branch point and branch cut. This new topological object is a relative of the topological object $N_{W}$ defined in Hayashi and Kondo (2019), with the difference that their object is defined in terms of the propagator itself, while our ${\cal N}_{\Gamma}$ depends on the mass function ${\cal M}(z)$, which in general can be defined by writing in full generality a fermion propagator as in (3), or for a bosonic propagator (stripping off all possible color/Lorentz tensorial structures)

[TABLE]

where the same comments as before apply to the wave-function normalization function $Z(p^{2})$ and mass function ${\cal M}(p^{2})$.

Despite the similarity, it is not clear to us if both objects encode the same physical information. Furthermore, notice that the integrand of Eq. (10) has at least two branch points, [math] and $\infty$, but depending on the analytical structure of ${\cal M}(z)$, the number of non-removable singularities can increase. Likewise, it seems that also the quantity defined in Hayashi and Kondo (2019) is susceptible to the appearance of new branch points, regarding the analytical expression of the propagator.

Therefore, we choose the closed contour $\Gamma$ as being composed of three paths with the following parametrization (see Fig. 2):

[TABLE]

It is important to notice that Eq. (10) can be rewritten as the winding number of some function $f(z)$, indeed,

[TABLE]

with

[TABLE]

Notice that the $\mathrm{Arctan\!}$-term gives the relevant contribution to the integral (15), in a sense that an equally well defined topological object (15) could be given by $\theta(z)=2\mathrm{Arctan\!}\left(\frac{{\cal M}(z)}{\sqrt{z}}\right)$. From Eq. (15) one can clearly see that our ${\cal N}_{\Gamma}$ is the winding number of $f(z)$ given at Eq. (16), and as such it represents the difference between the number of zeros and poles of $f(z)$ as the latter function is per assumption meromorphic within the contour $\Gamma$ and has no poles or zeros on the contour $\Gamma$.

IV.1 Invariance of $\mathcal{N}_{\Gamma}$ under small deformations of the parameters

Now, in order to derive the boundary condition for ${\cal N}_{\Gamma}$ to be invariant under smooth variations of ${\cal M}$ (with respect to the parameters of ${\cal M}$), we will rely on two procedures. In this Subsection, we use general complex analysis results, while the second argument, $cf.$ Appendix A, makes direct use of expressing Eq. (10) through the function $\theta(z)$ defined in Eq. (16).

If we consider the variation of $\mathcal{N}_{\Gamma}$, we get444The contour $\Gamma$ itself depends on the integrand, in the sense that $\Gamma$ is chosen in such way that it does not encircle or pierce any cuts or has poles on it. As before, this means that the renormalization group related logarithms affecting the integrand are analytic within $\Gamma$ and as such can be considered as smooth deformations. Stated otherwise, these logarithmic terms can again be neglected being irrelevant for determining the topological number ${\cal N}_{\Gamma}$.

[TABLE]

where $\mathcal{M}$ is a function of $z$; we did not write it explicitly to avoid notational clutter. In order to have something proportional to $\delta\mathcal{M}$, we need to integrate by parts the term with $\partial_{z}\delta\mathcal{M}$. This is achieved via

[TABLE]

Plugging this directly into $\delta\mathcal{N}_{\Gamma}$, we obtain

[TABLE]

The issue to prove that $\delta\mathcal{N}_{\Gamma}=0$ is the following: we must ascertain that the function $\frac{\partial}{\partial z}\left(\frac{\sqrt{z}z\delta\mathcal{M}}{(z+\mathcal{M}^{2})^{2}}\right)$ has a zero residue inside $\Gamma$, otherwise the zero value of the integral cannot be guaranteed. Therefore we consider that $\mathcal{M}$ has the form of a rational function (fermionic case), or the square root of some rational polynomial (bosonic case); we thus ignore, without any loss of generality, the logarithmic tails, as explained before. Then $\frac{\sqrt{z}z\delta\mathcal{M}}{(z+\mathcal{M}^{2})^{2}}$ is a meromorphic function times a possible square root of some rational polynomial. As the contour $\Gamma$ explicitly avoids the branch points and branch cuts coming from such a square root, we can say that $\frac{\sqrt{z}z\delta\mathcal{M}}{(z+\mathcal{M}^{2})^{2}}$ is a meromorphic function inside the contour $\Gamma$. Therefore, we are under the hypothesis of Lemma IV.1, at least inside $\Gamma$. A direct result of the lemma, with the proof sketched below, is that our $\frac{\partial}{\partial z}\left(\frac{\sqrt{z}z\delta\mathcal{M}}{(z+\mathcal{M}^{2})^{2}}\right)$ indeed has a zero residue inside $\Gamma$.

Lemma IV.1**.**

*Let $F(z)$ be a complex function with a pole at $z=a$ with finite multiplicity $m$ and analytic in the punctured disc $D^{\prime}(a,r)=\{z\in\mathbb{C},\;0<|z-a|<r\}$.

Then $F^{\prime}(z)$ has a pole with multiplicity $(m+1)$ and zero residue at $z=a$.*

Proof.

As the multiplicity of the pole is $m$, we may consider the Laurent series of $F(z)$ in $D^{\prime}(a,r)$ (see for instance Ablowitz and Fokas (2003))

[TABLE]

The coefficient $C_{-1}$ is the residue of $F(z)$ at $z=a$. Therefore,

[TABLE]

where $B_{n}=-(n-1)C_{-n+1}$ and $D_{n}=(n+1)C_{n+1}$. Thus

[TABLE]

with $E_{n}=B_{-n}$ for $n\leq-2$, $E_{n}=D_{n}$ for $n\geq 0$ and $E_{-1}=0$. Eq. (18) is, by construction, the Laurent series of $F^{\prime}(z)$ at $z=a$. Then, $z=a$ is a pole of multiplicity $(m+1)$, as the series starts with a term proportional to $(z-a)^{-m-1}$. It also has a residue equal to zero, as $E_{-1}=0$. ∎

V Applications

In this Section we will apply our complex topological object (10) on both fermionic and gluonic sectors as a matter of exercise, but, being the form of the propagator for Dirac fermions at finite temperature more complicated because of the absence of Lorentz invariance, the Dirac fermions will be treated at zero temperature while the gluonic case, only, will be treated at finite temperature.

V.1 A Dirac quark propagator that fits lattice data

Here we assume the analytic continuation $p^{2}\to z$ of the Dirac quark propagator that fits lattice data, whose $\mathcal{M}(p^{2})$ mass function is given by Eq. (9). The boundary conditions (27) for $\delta{\cal N}_{\Gamma}=0$ are satisfied, namely,

[TABLE]

We must emphasize that, because of the shape of our contour $\Gamma$, every possible complex pole with positive real part lies within the surrounded region (see Fig. 2). For the propagator on the right hand side of Eq. (3) and with dynamical mass given by Eq. (9), there will always be three poles, one of these being negatively real, and two complex conjugate ones. Notice that for the present quark propagator there will not appear new branch points due to the specific structure of ${\cal M}(z)$.

For general values of the parameters $(M^{3},m^{2},\mu)$, we find ${\cal N}_{\Gamma}=-2$ and ${\cal N}_{\Gamma}=+2$. More precisely, ${\cal N}_{\Gamma}=-2$ was associated to a negative mass at zero momentum, and also with a massless quark at zero momentum evolving to a strictly negative mass. It is possible to find combinations of the parameters so that the mass function is divergent at some point. In these cases, if the mass starts from negative values (considering $p^{2}$ from zero to positive values), ${\cal N}_{\Gamma}$ also acquires the value $-2$. On the other hand, we verified that ${\cal N}_{\Gamma}=+2$ is associated to configurations of the parameters such that the dynamical mass is positive at zero momentum, and also if it is a massless quark at zero momentum with an increasing positive mass. For example, for the specific values of $(M^{3},m^{2},\mu)$ that makes the quark propagator fit to lattice data, Furui and Nakajima (2006); Dudal et al. (2016), our topological invariant acquires the positive value $+2$.

V.2 A gluon propagator that fits lattice data

In order to explicitly show that our generalized topological object (10) can in fact be applied to the topological analysis of bosonic relativistic quantum fields within the MST framework, we perform here the topological classification of the space defined by the two-point Green function of gluons, whose analytic expression is one that has been used in literature to fit lattice data at zero and finite temperature, according to Aouane et al. (2012); Cucchieri et al. (2012b); Dudal et al. (2016, 2018). Such analytic expression can be derived by means of the refined Gribov-Zwanziger framework, see e.g. Dudal et al. (2008, 2019a).

Assuming the particular notation of Aouane et al. (2012), the gluon propagator reads,

[TABLE]

Performing the analytic continuation $p^{2}\to z$ to the complex plane of the mass function, ${\cal M}$, one has

[TABLE]

Notice that the boundary conditions (27) are satisfied,

[TABLE]

so that ${\cal N}_{\Gamma}$ is indeed a topological invariant under smooth variation of ${\cal M}$ with respect to the parameters $(r,\ b,\ d,\ c)$.

As mentioned before, since our proposed ${\cal N}_{\Gamma}$ is constructed in the complex plane and explicitly depends on the expression of the (analytic continuation of the) mass function of the propagator ${\cal M}(p^{2})$, it is quite possible that non-removable singularities appear, others than [math] and $\infty$, as in the case of the present example. Given the expression of the complex mass function ${\cal M}(z)$, at Eq. (21), one can check that the integrand of Eq. (10) has also the new branch points $-\frac{1}{d}$ and $\frac{-b^{2}d-dr2^{2}}{-1+2dr^{2}}$, which are real and, therefore, lie outside the region bounded by $\Gamma$.

In this example we restrict ourselves to the finite temperature results of Aouane et al. (2012). The values of $(r,~{}b,~{}d,~{}c)$ are taken from their Table II. Namely, the authors of Aouane et al. (2012) assessed the following values of temperature (in unities of the critical temperature): $T/T_{c}=0.65$; $T/T_{c}=0.74$; $T/T_{c}=0.86$; $T/T_{c}=0.99$; $T/T_{c}=1.20$; $T/T_{c}=1.48$; $T/T_{c}=1.98$; and $T/T_{c}=2.97$. For all the temperature values with the exception of the highest three, they set $b=0.0$. Notice that for $T/T_{c}=1.20$ (clearly above the critical temperature) the parameter $b$ is also set to zero.

For the values reported in Aouane et al. (2012), our topological invariant acquires the following values: ${\cal N}_{\Gamma}=0$ whenever $b=0$ , since the propagator then displays only a (negative) real (double) pole. Such configuration appears for temperatures $T/T_{c}=0.65$; $T/T_{c}=0.74$; $T/T_{c}=0.86$; $T/T_{c}=0.99$; $T/T_{c}=1.20$. On the other hand, we find ${\cal N}_{\Gamma}=+2$ for the three highest temperatures $T/T_{c}=1.48$; $T/T_{c}=1.98$; and $T/T_{c}=2.97$ (where $b\neq 0$). Therefore, these results suggest a phase transition at a critical temperature somewhere between $T/T_{c}=1.20$ and $T/T_{c}=1.48$, different from the deconfinement phase transition found by the authors of Aouane et al. (2012).

Interestingly, the authors of Hayashi and Kondo (2019) found $N_{W}=0$ for gauge propagators displaying only real poles; and $N_{W}=-2$ for gauge propagators with a set of complex conjugate poles. Here, again, we point out that our topological invariant assumes the same absolute values as the one of Hayashi and Kondo (2019), despite that they do not necessarily represent the same topological quantity. Therefore, following a completely different approach from the one in Hayashi and Kondo (2019) (although both are aimed to the momentum space topological analysis of the gauge field two-point Green’s function), we could find a phase transition between different regimes of the gauge propagator, in agreement with Hayashi and Kondo (2019).

VI Conclusions

We developed a topological classification of relativistic quantum systems within the momentum space topology framework, starting with the usual topological invariant ${\cal N}_{3}$ as applied to Dirac fermions. This topological invariant can experience integer jumps when one crosses different open regions, which are limited by phase boundaries. This change can be related to the non-analyticities appearing in the integrand of (6).

Afterwards, we proposed a new topological invariant, ${\cal N}_{\Gamma}$, that is sensitive to the existence of complex poles of the propagator of the relativistic quantum field, be it fermion or boson. This topological invariant was constructed by performing the analytic continuation $p^{2}\to z$ to the complex plane of the usual ${\cal N}_{3}$. As a result, ${\cal N}_{\Gamma}$ depends on the closed contour $\Gamma$ and on the analytic expression of ${\cal M}(z)$. In order to ensure the invariance of ${\cal N}_{\Gamma}$ with respect to smooth variations of ${\cal M}(z)$, we found that ${\cal M}(z)$ and $\delta{\cal M}(z)$ must be constant in the limit $|z|\to\infty$. We applied this topological number ${\cal N}_{\Gamma}$ to a model Dirac quark propagator whose analytic expression can fit quite well corresponding lattice data. In addition, we also applied it to a gluon propagator form that can fit zero and finite temperature lattice data. We made use of the rational function fits according to Aouane et al. (2012), and we could find a phase transition at a temperature within the range $T/T_{c}=1.20$ and $T/T_{c}=1.48$, which is slightly above the thermodynamic critical (deconfinement) temperature found by the authors.

While Eq. (1) has the clear physical meaning of being related to a non-dissipative conductivity (Hall, spin Hall, etc, $cf.$ Zubkov and Wu (2019); Zhang and Zubkov (2019)) at this stage the physical interpretation of our invariant is less obvious. What is clear is that it is sensitive to the presence of complex conjugate poles thus, as shown in refs. Hayashi and Kondo (2019) and Cyrol et al. (2018), it “measures” the possible occurrence of a violation of the “reflection positivity condition” (the Euclidean counterpart to the positive definiteness of the norm in the Hilbert space of the corresponding Wightman quantum field theory, see the Osterwalder-Schrader axioms Osterwalder and Schrader (1973, 1975)) in terms of the Schwinger function. Violation of the reflection positivity condition is regarded as a signal for gluon confinement (see Section VI of Hayashi and Kondo (2019) and Section V and VI of Cyrol et al. (2018), as well as Krein et al. (1992).).

For the particular propagator (20), depending on the value of the parameters, we will encounter either two (possibly degenerate) real poles or two complex conjugate ones. Only in the latter case, we will have a contribution to the topological invariant as defined in expression (10). Moreover, in the presence of a set of such poles, it has been shown explicitly in e.g. Dudal et al. (2008); Roberts and Williams (1994); Gao et al. (2018) that the Schwinger function,

[TABLE]

is no longer positive definite, in accordance with lattice measurement of the same quantity, see for example Cucchieri et al. (2005); Bowman et al. (2007).

For these reasons we believe that our topological invariant proposed in equation (10) may be interpreted as an indication of whether the quantum system is in a confined or deconfined phase. However, for the moment, this is only an intuition without a rigorous analytic proof, which deserves further investigation.

For future work, it would be interesting to use improved fits to more recent finite temperature gluon data, to see if the phase transition coincides with the deconfinement transition. That such is possible can be inferred from the preliminary fitting report of Cucchieri et al. (2012b), where in contrast to the here used numbers of Aouane et al. (2012), there are complex poles below (up to $T=0$) and around $T_{c}$. More takes on the analytic structure of gluon and/or quark propagators can be found in $e.g.$ Dudal et al. (2014); Siringo (2016); Lowdon (2018); Silva et al. (2017); Ilgenfritz et al. (2018); Cyrol et al. (2018); Kondo et al. (2019); Dudal et al. (2019b); Binosi and Tripolt (2019). Although these rational function fits never capture the renormalization group-controlled logarithmic tails, we explained how our construct is insensitive to these structures anyhow.

In addition, it would be rather interesting to follow the temperature-evolution, if any, of the quark topological number based on the rather recent lattice output of Oliveira and Silva (2019), to see if one can find a topological signal of the deconfinement and/or chiral transition in the quark sector. We hope to come back to these issues in the foreseeable future.

Acknowledgments

P.P. , D.D. and L.R. acknowledge the warm hospitality of the Centro de Estudios Científicos (CECs) in Valdivia, Chile, during different stages of this work. P. P. is supported by the project High Field Initiative (CZ.02.1.01/0.0/0.0/15_003/0000449) from the European Regional Development Fund. I. F. J. is supported by Fondecyt Grant No. 3170278. This work was partially supported by the research Grant Number 2017W4HA7S “NAT-NET: Neutrino and Astroparticle Theory Network” under the program PRIN 2017 funded by the Italian Ministero dell’ Istruzione, dell’ Università e della Ricerca (MIUR) and also by the Initiative “TAsP: Theoretical Astroparticle Physics” funded by the Italian Istituto Nazionale di Fisica Nucleare (INFN). This work has been funded by the Fondecyt Grant No. 1160137. The Centro de Estudios Científicos (CECs) is funded by the Chilean Government through the Centers of Excellence Base Financing Program of Conicyt.

Appendix A Alternative proof that $N_{\Gamma}$ is a homotopic invariant

For a second, more explicit, demonstration of the invariance of Eq. (10), we split our contour as $\Gamma=C(\alpha)+L_{-}(t)+L_{+}(t)$ at the limit $R\to\infty$, as depicted in Fig. 2, writing ${\cal N}_{\Gamma}$ as

[TABLE]

Given the parameterizations (12), (13) and (14), one has,

[TABLE]

with $R\to\infty$. Let us focus on the integral over the contour piece $C(\alpha)$. As we avoid poles and branch points/cuts on $\Gamma$, such an integral reduces to

[TABLE]

if ${\cal M}(z)$ is at most constant at $|z|\to\infty$. This means that our topological invariant reduces to

[TABLE]

The variation $\delta_{{}_{\cal M}}{\cal N}_{\Gamma}$, due to $\delta{\cal M}(z)$, on Eq. (25), lead us to

[TABLE]

The variation $\delta_{{}_{\cal M}}\theta(z)$ reads

[TABLE]

just as Eq. (17). Thus, as there are no singularities along the contour $\Gamma$, in order to ensure $\delta_{{}_{\cal M}}{\cal N}_{\Gamma}=0$, the analytical continuation of the mass function ${\cal M}(z)$ and its small variation $\delta{\cal M}(z)$ must be constant at the boundary $|z|\to\infty$, i.e.,

[TABLE]