Next Article in Journal
Tripartite Entanglement: Foundations and Applications
Next Article in Special Issue
Study of Very Forward Neutrons with the CMS Zero Degree Calorimeter
Previous Article in Journal
Driving Interactions Efficiently in a Composite Few-Body System
Previous Article in Special Issue
Open-Charm Hadron Measurements in Au+Au Collisions at √sNN = 200 GeV by the STAR Experiment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temperature Dependence of the Axion Mass in a Scenario Where the Restoration of Chiral Symmetry Drives the Restoration of the UA(1) Symmetry

1
Physics Department, Faculty of Science-PMF, University of Zagreb, Bijenička cesta 32, 10000 Zagreb, Croatia
2
Ruđer Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Universe 2019, 5(10), 208; https://doi.org/10.3390/universe5100208
Submission received: 2 August 2019 / Revised: 30 September 2019 / Accepted: 1 October 2019 / Published: 8 October 2019
(This article belongs to the Special Issue The Zimányi School and Analytic Hydrodynamics in High Energy Physics)

Abstract

:
The temperature (T) dependence of the axion mass is predicted for T s up to 2.3 × the chiral restoration temperature of QCD. The axion is related to the U A ( 1 ) anomaly. The squared axion mass m a ( T ) 2 is, modulo the presently undetermined scale of spontaneous breaking of Peccei–Quinn symmetry f a (squared), equal to QCD topological susceptibility χ ( T ) for all T. We obtain χ ( T ) by using quark condensates calculated in two effective Dyson–Schwinger models of nonperturbative QCD. They exhibit the correct chiral behavior, including the dynamical breaking of chiral symmetry and its restoration at high T. This is reflected in the U A ( 1 ) symmetry breaking and restoration through χ ( T ) . In our previous studies, such χ ( T ) yields the T-dependence of the U A ( 1 ) -anomaly-influenced masses of η and η mesons consistent with experiment. This in turn supports our prediction for the T-dependence of the axion mass. Another support is a rather good agreement with the pertinent lattice results. This agreement is not spoiled by our varying u and d quark mass parameters out of the isospin limit.

1. Introduction

The axion, one of the oldest hypothetical particles beyond the Standard Model, intensely sought for by many experimentalists already for 40 years now, still escapes detection [1]. It was introduced theoretically [2,3,4,5] to solve the so-called Strong CP problem of QCD. The problem is that no experimental evidence of CP-symmetry violation has been found in strong interactions, although the QCD Lagrangian L QCD ( x ) can include the so-called θ -term L θ ( x ) = θ Q ( x ) where gluon field strengths F μ ν b ( x ) form the CP-violating combination Q ( x ) named the topological charge density:
Q ( x ) = g 2 32 π 2 F μ ν b ( x ) F ˜ b μ ν , where F ˜ b μ ν 1 2 ϵ μ ν ρ σ F ρ σ b ( x ) .
Whereas Q ( x ) can be re-cast in the form of a total divergence μ K μ , discarding L θ is not justified even if F μ ν b ( x ) vanish sufficiently fast as | x | . Specifically, F F ˜ = μ K μ can anyway contribute to the action integral, since in QCD there are topologically nontrivial field configurations such as instantons. They are important for, e.g., obtaining the anomalously large mass of the η meson. Also, precisely the form (1) from the θ -term appears in the axial anomaly, breaking the U A ( 1 ) symmetry of QCD - see Equation (2).
For these reasons, one needs L θ ( x ) = θ Q ( x ) in the QCD Lagrangian, as reviewed briefly in Section 1 of Ref. [6]. Moreover, the Strong CP problem cannot be removed by requiring that the coefficient θ = 0 , since QCD is an integral part of the Standard Model, where weak interactions break the CP symmetry. This CP violation comes from the complex Yukawa couplings, yielding the complex CKM matrix [7,8] and the quark-mass matrix M which is complex in general. To go to the mass–eigenstate basis, one diagonalizes the mass matrix, and the corresponding chiral transformation changes θ by arg det M . Hence, in the Standard Model the coefficient of the Q F F ˜ term is in fact θ ¯ = θ + arg det M [6]. Therefore, to be precise, we change our notation to θ ¯ -term, L θ L θ ¯ .
Since CP is not a symmetry of the Standard Model, there is no a priori reason θ ¯ , which results from the contributions from both the strong and weak interactions, should vanish. And yet, the experimental bound on it is extremely low, | θ ¯ | < 10 10 [9], and in fact consistent with zero. Therefore, the mystery of the vanishing strong CP violation is: why is θ ¯ so small?
The most satisfactory answer till this very day has been provided by axions, even though the original variant has been ruled out [1]. In the meantime, they turned out to be very important also for cosmology, as promising candidates for dark matter—see from relatively recent references such as [10,11] to the earliest papers [12,13,14]. (For an example of a broader review of axion physics, see [15].) It is thus no wonder that ever since the original proposal of the axion mechanism [2,3,4,5] in 1977–1978, many theorists kept developing various ideas on this theoretically much needed object, trying to pinpoint the properties of this elusive particle and increase chances of finding it.
However, to no avail. There have even been some speculations that the axion is hidden in plain sight, by being experimentally found, paradoxically, already years before it was conjectured theoretically: namely, that the axion should in fact be identified with the well-known η meson with a minuscule admixture of a pseudoscalar composite of neutrinos [16]. Nevertheless, while an intimate relation between the axion and η doubtlessly exists, reformulations of the axion theory, let alone so drastic ones, are in fact not needed to exploit this axion- η relationship: thanks to the fact that both of their masses stem from the axial anomaly and are determined by the topological susceptibility of QCD, in the present paper we show how our previous study [17] of the temperature (T) dependence of the η and η mesons give us a spin-off in the form of the T-dependence of the axion mass, m a ( T ) . It is given essentially by the QCD topological susceptibility χ ( T ) , which is rather sensitive to changes of the lightest quark masses m q : Equation (9) vanishes linearly when m q 0 even for just one flavor q. We thus examine the effect of their values on χ ( T ) also out of the isosymmetric limit, and find that such a variation can be accommodated well. The agreement with lattice results on χ ( T ) is reasonably good.

2. Connection with the Complex of the η and η Mesons

2.1. Some Generalities on the Influence of the Anomaly on η and η

In this paper, we neglect contributions of quark flavors heavier than q = s and take N f = 3 as the number of active flavors.
At vanishing and small temperatures, T 0 , the physical η meson is predominantly1 η 0 , the singlet state of the flavor S U ( 3 ) group, just like its physical partner, the lighter isospin-zero mass eigenstate η is predominantly the octet state η 8 . Unlike the S U ( 3 ) octet states π , K and η 8 , the singlet η 0 is precluded from being a light (almost-)Goldstone boson of the dynamical breaking of the (only approximate) chiral symmetry of QCD (abbreviated as DChSB). Namely η 0 receives a relatively large anomalous mass contribution from the non-Abelian axial ABJ2 anomaly, or gluon anomaly for short. An even better name for it is the U A ( 1 ) anomaly, since it breaks explicitly the U A ( 1 ) symmetry of QCD on the quantum level.
The breaking of U A ( 1 ) by the anomaly makes the flavor singlet ( a = 0 ) axial current of quarks, A 0 μ ( x ) = q = u , d , s q ¯ ( x ) γ μ γ 5 q ( x ) , not conserved even in the chiral limit:
μ A 0 μ ( x ) = i q = u , d , s 2 m q q ¯ ( x ) γ 5 q ( x ) + 2 N f Q ( x ) , ( N f = 3 ) ,
unlike the corresponding octet currents A a μ ( x ) , a = 1 , 2 , , N f 2 1 . In the chiral limit, the current masses of non-heavy quarks all vanish, m q 0 ( q = u , d , s ) , but the divergence of the singlet current A 0 μ is not vanishing due to the U A ( 1 ) anomaly contributing no other but the topological charge density operator Q ( x ) (1)—that is, precisely the quantity responsible for the strong CP problem.
The quantity related to the U A ( 1 ) -anomalous mass in the η - η complex is the QCD topological susceptibility χ ,
χ = d 4 x 0 | T Q ( x ) Q ( 0 ) | 0 ,
where T denotes the time-ordered product.
Figure 1 shows how anomaly contributes to the mass matrix (in the basis of quark–antiquark ( q q ¯ ) pseudoscalar bound states P) by depicting how hidden-flavor q q ¯ pseudoscalars mix, transiting through anomaly-dominated gluonic intermediate states.

2.2. On Some Possibilities of Modeling the U A ( 1 ) Anomaly Influence

Light pseudoscalar mesons can be studied by various methods. We have preferred using [17,21,22,23,24,25,26,27,28,29,30,31,32] the relativistic bound-state approach to modeling nonperturbative QCD through Dyson–Schwinger equations (DSE), where, if approximations are consistently formulated, model DSE calculations also reproduce the correct chiral behavior of QCD. This is of paramount importance for descriptions of the light pseudoscalar mesons, which are quark–antiquark bound states but simultaneously also the (almost-)Goldstone bosons of DChSB of QCD. (For general reviews of the DSE approach, see, e.g., Refs. [33,34,35,36]. About our model choice at T = 0 and T > 0 , in the further text see especially Appendix A).
Figure 1 illustrates how hard computing would be “in full glory” the U A ( 1 ) -anomalous mass and related quantities, such as the presently all-important topological susceptibility (3) in the DSE approach with realistically modeled QCD interactions—especially if the calculation should be performed in a consistent approximation with the calculation of the light pseudoscalar bound states, to preserve their correct chiral behavior. (For this reason, they have most often been studied in the rainbow-ladder approximation of DSE, which is inadequate for the anomalous contributions [33,34,35,36], as also Figure 1 shows.)
However, our DSE studies of pseudoscalar mesons have been able to address not only pions and kaons, but also η and η mesons, for which it is essential to include the anomalous U A ( 1 ) symmetry breaking at least at the level of the masses. This was done as described in Refs. [21,24,29,30,37], namely exploiting the fact that the U A ( 1 ) anomaly is suppressed in the limit of large number of QCD colors N c [38,39]. This allows treating the anomaly contribution formally as a perturbation with respect to the non-anomalous contributions to the η and η masses [21,24,29]. This way we avoid the need to compute the anomalous mass contribution together, and consistently, with the non-anomalous, chiral-limit-vanishing parts of the masses. The latter must be evaluated by some appropriate, chirally correct method, and our preferred tool—the relativistic bound-state DSE approach [33,34,35,36]—is just one such possibility. The point is that they comprise the non-anomalous part of the η - η mass matrix, to which one can add, as a first-order perturbation, the U A ( 1 ) -anomalous mass contribution M U A ( 1 ) —and it does not have to be modeled, but taken from lattice QCD [29]. Specifically, at T = 0 , M U A ( 1 ) can be obtained from χ YM , the topological susceptibility of the (pure-gauge) Yang-Mills theory, for which reliable lattice results have already existed for a long time3 [40,41,42].
This can be seen from the remarkable Witten–Veneziano relation (WVR) [38,39] which in a very good approximation relates the full-QCD quantities ( η , η and K-meson masses M η , M η and M K respectively, and the pion decay constant f π ), to the pure-gauge quantity χ YM :
M η 2 + M η 2 2 M K 2 = 2 N f χ YM f π 2 M U A ( 1 ) 2 .
The right-hand-side must be the total U A ( 1 ) -anomalous mass contribution in the η - η complex, since in the combination on the left-hand-side everything else cancels at least to the second order, O ( m q 2 ) , in the current quark masses of the three light flavors q = u , d , s . This is because the non-anomalous, chiral-limit-vanishing parts M q q ¯ of the masses of pseudoscalar mesons4 P q q ¯ composed of sufficiently light quarks, satisfy the Gell–Mann–Oakes–Renner (GMOR) relation with their decay constants f q q ¯ and the quark–antiquark ( q q ¯ ) condensate signaling DChSB:
M q q ¯ 2 = q ¯ q 0 ( f q q ¯ ch . lim ) 2 ( m q + m q ) + O ( m q 2 , m q 2 ) ( q , q = u , d , s ) .
Here f q q ¯ ch . lim = f q q ¯ ( m q , m q 0 ) , and q ¯ q 0 denotes the massless-quark condensate, i.e., the q q ¯ chiral-limit condensate, or “massless” condensate for short. (In the absence of electroweak interactions, the “massless” condensates have equal values for all flavors: q ¯ q 0 = q ¯ q 0 .) It turns out that even s-flavor is sufficiently light for Equation (5) to provide reasonable approximations.
Using WVR and χ YM to get the anomalous part of the η and η masses is successful [29] only for T 0 , or at any rate, T’s well below T c , the pseudocritical temperature of the chiral transition. In the absence of a systematic re-derivation of WVR (4) at T > 0 , its straightforward extension (simply replacing all quantities by their T-dependent versions) is tempting, but was found [31] unreliable and with predictions in a drastic conflict with experiment when T starts approaching T c . This is because the full-QCD quantities M η ( T ) , M η ( T ) , M K ( T ) and f π ( T ) have very different T-dependences from the remaining quantity χ YM ( T ) , which is pure-gauge and thus much more resilient to increasing temperature: the critical temperature of the pure-gauge, Yang-Mills theory, T Y M , is more than 100 MeV higher than QCD’s T c = ( 154 ± 9 ) MeV [44,45]. The early lattice result T Y M 260 MeV [40,41] is still accepted today [46], and lattice groups finding a different T Y M claim only it is even higher, for example T Y M = ( 300 ± 3 ) MeV of Gattringer et al. [42]. (There are even some claims about experimentally established T Y M = 270 MeV [47].)
We thus proposed in 2011. [48] that the above mismatch of the T-dependences in WVR (4) can be removed if one invokes another relation between χ YM and full-QCD quantities, to eliminate χ YM , i.e., substitute pertinent full-QCD quantities instead of χ YM at T > 0 . This is the Leutwyler–Smilga (LS) relation, Equation (11.16) of Ref. [49], which we used [17,37,48] in the inverted form (and in our notation):
χ YM = χ 1 + χ ( 1 m u + 1 m d + 1 m s ) 1 q ¯ q 0 ( χ ˜ ) ,
to express (at T = 0 ) pure-gauge χ YM in terms of the full-QCD topological susceptibility χ χ QCD , the current quark masses m q , and q ¯ q 0 , the condensate of massless, chiral-limit quarks. The combination which these full-QCD quantities comprise, i.e., the right-hand-side of the LS relation (6), we denote (for all T) by the new symbol χ ˜ for later convenience - that is, for usage at high T, where the equality (6) with χ YM does not hold.
The remarkable LS relation (6) holds for all values of the current quark masses. In the limit of very heavy quarks, it correctly yields χ χ quenched = χ YM for m q , it but it also holds for the light m q . In the light-quark sector, the QCD topological susceptibility χ can be expressed as [49,50,51]:
χ = q ¯ q 0 1 m u + 1 m d + 1 m s + C m ,
where C m represents corrections of higher orders in light-quark masses m q . Thus, it is small and often neglected, leaving just the leading term as the widely used [52] expression for χ in the light-quark, N f = 3 sector. However, setting C m = 0 in the light-quark χ (7) returns us χ YM = through Equation (6) [48]. Or conversely, setting χ YM = in the LS relation (6), gives the leading term of χ (7) (see also Ref. [43]). This can be a reasonable, useful limit considering that in reality χ YM / χ 40 . Nevertheless, in our previous works on the η - η complex at T > 0 [17,48] we had to fit C m (and parameterize it with Ansätze at T > 0 ), since we needed the realistic value of χ YM ( T = 0 ) = χ ˜ ( T = 0 ) from lattice to reproduce the well-known masses of η and η at T = 0 . (However, just for χ ( T ) this is not necessary.)
Replacing [48] χ YM ( T ) by the full-QCD quantity χ ˜ ( T ) obviously keeps WVR at T = 0 , but avoids the ’YM vs. QCD’ T-dependence mismatch with f π ( T ) and the LHS of Equation (4), so it is much more plausible to assume the straightforward extension of T-dependences. The T-dependences of χ ( T ) (7) and χ ˜ ( T ) , and thus also of the anomalous parts of the η and η masses, are then obviously dictated by q ¯ q 0 ( T ) , the “massless” condensate.
General renormalization group arguments suggest [53] that QCD with three degenerate light-quark flavors has a first-order phase transition in the chiral limit, whereas in QCD with (2+1) flavors (where s-quark is kept significantly more massive) a second-order chiral-limit transition5 is also possible and even more likely [62,63,64]. What is important here, is that in any case the chiral-limit condensate q ¯ q 0 ( T ) drops sharply to zero at T = T c . (The dotted curve in Figure 2 is just a special example thereof, namely q ¯ q 0 ( T ) calculated in Ref. [48] using the same model as in [17] and here.) This causes a similarly sharp drop of χ ˜ ( T ) and χ ( T ) . (We may be permitted to preview similar dotted curves in Figures 5 and 7 in the next section and anticipate that in this case one would get a massless axion at T = T c .) This was also the reason, besides the expected [65,66] drop of the η mass, Ref. [48] also predicted so drastic drop of the η mass at T = T c , that it would become degenerate with the pion. However, no experimental indication whatsoever for a decreasing behavior of the η mass, and much less for such a conspicuous sharp mass drop, has been noticed to this day, which seems to favor theoretical descriptions with a smooth crossover. Also, recent lattice QCD results (see [44,67,68] and their references) show that the chiral symmetry restoration is a crossover around the pseudocritical transition temperature T c .
To describe such a crossover behavior of the chiral transition, we incorporated [17,37] into our approach Shore’s generalization [69] of the Witten–Veneziano relation. To be precise, we studied it at T = 0 already in 2008 [30] and adapted it to our DSE bound-state context by applying some very plausible simplifications [20]. The more recent reference [37] presented the analytic, closed-form solutions to Shore’s equations for the pseudoscalar meson masses. These solutions showed that Shore’s approach is then actually quite similar to the original WVR, leading to a similar η - η mass matrix [37].
Presently, the most important advantage is that Shore’s generalization leads to the crossover T-dependence. As obtained in Section 3 for two specific interactions modeling the nonperturbative QCD interaction, these condensates exhibit a smooth, crossover chiral symmetry transition around T c . Here Figure 2 illustrates this generic behavior by displaying the results obtained in Section 3 for a specific DSE model: the higher current quark mass, the smoother the crossover behavior, which then results in the crossover behavior also of other quantities, like the presently all-important quantity, the QCD topological susceptibility χ ( T ) .
This comes about as follows: the quantity which in Shore’s mass relations [69] has the role of χ YM in the Witten–Veneziano relation, is called the full-QCD topological charge parameter A. (Shore basically took over this quantity from Di Vecchia and Veneziano [50].) At T = 0 , it is approximately equal to χ YM in the sense of 1 / N c expansion. Shore uses A to express the QCD susceptibility χ through a relation similar to the Leutwyler–Smilga relation (see Equations (2.11) and (2.12) in Ref. [69]), but using the condensates u ¯ u , d ¯ d , s ¯ s of realistically massive u , d , s quarks. The inverse relation, yielding A (with the opposite sign convention), is the most illustrative for us:
A = χ 1 + χ ( 1 m u u ¯ u + 1 m d d ¯ d + 1 m s s ¯ s ) ( A = χ YM + O ( 1 N c ) at T = 0 ) .
Obviously, it is analogous to the inverted LS relation (6) defining χ ˜ , except that A is expressed through “massive” condensates. (If they are all replaced by q ¯ q 0 , then A χ ˜ .) They are in principle different for each flavor, but in the limit of usually excellent isospin symmetry, u ¯ u = d ¯ d .
One can examine the limiting assumption A = in analogy with taking the limit χ YM = compared to χ = χ Q C D . Then, for A = (be it in our Equation (8) or Shore’s Equations (2.11) and (2.12) for χ ), one recovers the leading term of the QCD topological susceptibility expressed by the “massive” condensates. However, if one needs a finite A, as in η - η calculations [17] where one needs to reproduce A χ Y M , one also needs the appropriate correction term C m , just as C m in Equation (7), so that:
χ ( T ) = 1 1 m u u ¯ u ( T ) + 1 m d d ¯ d ( T ) + 1 m s s ¯ s ( T ) + C m .
Again, C m is a very small correction term of higher orders in the small current quark masses m q ( q = u , d , s ), and we can neglect it in the present context, where we actually have a simpler task than finding T-dependence of the η and η masses in Ref. [17]. Since it turns out that for determining the T-dependence of the mass of the QCD axion we do not need to find A, we set C m = 0 in this paper throughout. One needs just the topological susceptibility χ ( T ) for that, and just the leading term of (9) will suffice to yield the crossover behavior found on lattice ( e.g., in Refs. [70,71,72]).

3. The Axion Mass from the Non-Abelian Axial Anomaly of QCD

Peccei and Quinn (PQ) introduced [2,3] a new global symmetry U ( 1 ) PQ which is broken spontaneously at some very large, but otherwise still unknown scale f a > 10 8 GeV [1,73], which determines the absolute value of the axion mass m a . Nevertheless, this constant cancels from ratios such as m a ( T ) / m a ( 0 ) , where T is temperature. Thus, useful insights and applications, such as those involving the nontrivial part of axion T-dependence, are possible in spite of f a being presently unknown.
The factor in the axion mass which carries the nontrivial T-dependence, is the QCD topological susceptibility χ ( T ) . This quantity is also essential for our description of the η - η complex at T > 0 , since it relates the T-dependence of the anomalous breaking of U A ( 1 ) symmetry.

3.1. The Axion as the Almost-Goldstone Boson of the Peccei–Quinn Symmetry

The pseudoscalar axion field a ( x ) arises as the (would-be massless) Goldstone boson of the spontaneous breaking of the PQ symmetry U ( 1 ) PQ [4,5]. The axion contributes to the total Lagrangian its kinetic term and its interaction with fermions of the Standard model. Nevertheless, what is important for the resolution of the strong CP problem, is that the axion also couples to the topological charge density operator Q ( x ) defined in Equation (1) and generating the U A ( 1 ) -anomalous term in Equation (2). The θ ¯ -term in L Q C D thus changes into
L θ ¯ L axion θ ¯ + = θ ¯ + a f a g 2 64 π 2 ϵ μ ν ρ σ F μ ν b F ρ σ b .
Because of this axion-gluon coupling, the U ( 1 ) PQ symmetry is also broken explicitly by the U A ( 1 ) anomaly (gluon axial anomaly). This gives the axion a nonvanishing mass, m a 0 [4,5].
Gluons generate an effective axion potential, and its minimization leads to the axion expectation value a which makes the modified coefficient of Q ( x ) in Equation (10) vanish: θ ¯ + a / f a θ ¯ = 0 .
Obviously, the experiments excluding the strong CP violation, such as [9], have in fact been finding that consistent with zero is θ ¯ , the coefficient of Q ( x ) in the QCD Lagrangian when the axion exists. The strong CP problem is thereby solved, irrespective of the initial value of θ ¯ . (Relaxation from any θ ¯ -value in the early Universe towards the minimum at θ ¯ = a / f a is called misalignment production. The resulting axion oscillation energy is a good candidate for cold dark matter [10,11,12,13,14,15].)

3.2. Axion Mass from the Topological Susceptibility from Condensates of Massive Quarks

Modulo the (squared) Peccei–Quinn scale f a 2 , the axion mass squared is at all temperatures T given by the QCD topological susceptibility [11,43,70,71,72,74] very accurately [75,76] (up to negligible corrections of the order ( p i o n m a s s ) 2 / f a 2 ):
m a 2 ( T ) = 1 f a 2 χ ( T ) ,
as revealed by the quadratic term of the expansion of the effective axion potential.
We explained in Section 2 how the U A ( 1 ) symmetry-breaking quantity χ ( T ) can be obtained through Equation (9) as a prediction of any method which can provide the quark condensates q ¯ q ( T ) ( q = u , d , s ). Thus, one can get the T-dependence of the axion mass (11) through the mechanism where DChSB drives the U A ( 1 ) symmetry breaking. (And conversely, of course: the chiral restoration then drives the restoration of U A ( 1 ) symmetry.)
An excellent tool to study DChSB, and in fact “produce” it in the theoretical sense, is one of the basic equations of the DSE approach - the gap equation. The most interesting thing it does for nonperturbative QCD is explaining the notion of the constituent quark mass around 1 3 of the nucleon mass M N by generating them via DChSB, in the same process which produces the q q ¯ condensates. Thanks to this, q q ¯ condensates can be evaluated from dressed quark propagators. Specifically, hadronic-scale large ( M N / 3 at small momenta p) dressed quark-mass functions M q ( p 2 ) B q ( p 2 ) / A q ( p 2 ) are generated despite two orders of magnitude lighter current quark masses m q , and in fact even in the chiral limit, when m q = 0 ! This happens in low-energy QCD thanks to nonperturbative dressing via strong dynamics, making strongly dressed quark propagators S q ( p ) out of the free quark propagators S q free :
S q free ( p ) = 1 i γ · p + m q S q ( p ) = 1 i γ · p A q ( p 2 ) + B q ( p 2 ) ( Euclidean space expressions ) .
The solution for the dressed quark propagator S q ( p ) of the flavor q, i.e., the dressing functions A q ( p 2 ) and B q ( p 2 ) , are found by solving the gap equation
S q 1 ( p ) = S q free ( p ) 1 Σ q ( p ) , ( q = u , d , s ) ,
where Σ q ( p ) is the corresponding DChSB-generated self-energy, for example, Equation (14) if the rainbow-ladder truncation is adopted.
In the present work, all we want to model of nonperturbative QCD are the condensates q ¯ q at all temperatures T, and for that the solutions of the quark-propagator gap Equation (13) are sufficient, i.e., we do not need the Bethe–Salpeter equation (BSE) for the q q ¯ pseudoscalar bound states. However, we want the same condensates, and basically the same (leading term of) the topological susceptibility as we had in our related η - η paper [17], and in a number of earlier papers, such as Refs. [31,48]. Thus, we now use the same model interaction we have been using then in the consistent rainbow-ladder (RL) truncation of DSE’s to produce chirally correctly behaving pseudoscalar mesons - that is, with the non-anomalous parts of their masses given by GMOR (5).
Thus, the quark self-energy in the gap Equation (13) in the RL truncation is
Σ q ( p ) = d 4 ( 2 π ) 4 g 2 D μ ν a b ( p ) eff λ a 2 γ μ S q ( ) λ b 2 γ ν ,
where D μ ν a b ( k ) eff is an effective gluon propagator, which should be chosen to model the nonperturbative, low-energy domain of QCD. This can be done in varying degrees of DSE modeling, depending on the variety of problems one wants to treat [33,34,35,36,77]. For example, in the context of low-energy meson phenomenology, if one does not aim to address problems of perturbative QCD, it is better not to include the perturbative part of the QCD interaction. Otherwise, in the words of very authoritative DSE practitioners, “the logarithmic tail and its associated renormalization represent an unnecessary obfuscation" [78].
In medium, the original O ( 4 ) symmetry is broken to O ( 3 ) symmetry. The most general form of the dressed quark propagator then has four independent tensor structures and four corresponding dressing functions. At nonvanishing temperature, T > 0 , we use the Matsubara formalism, where four-momenta decompose into three-momenta and Matsubara frequencies: p = ( p 0 , p ) p n = ( ω n , p ) . Therefore, the (inverted) dressed quark propagator S q ( p ) (13) becomes
S q 1 ( p n ) = S q 1 ( p , ω n ) = i γ · p A q ( p 2 , ω n ) + i ω n γ 4 C q ( p 2 , ω n ) + B q ( p 2 , ω n ) + i ω n γ 4 γ · p D q ( p 2 , ω n ) .
(The T-dependence of the propagator dressing functions is understood and, to save space, is not indicated explicitly, except in Appendix A).
Nevertheless, the last dressing function D q ( p 2 , ω n ) is so very small that it is quite safe and customary to neglect it—e.g., see Refs. [34,79]. Thus, also we set D q 0 , leaving only A q , C q and B q .
For applications in involved contexts, such as calculations at T > 0 , appropriate simplifications are very welcome for tractability. This is why in Refs. [17,31,48] and presently, we adopted relatively simple, but phenomenologically successful [30,31,32,77,80] separable approximation [77]. The details on the functional form and parameters of the presently used model interaction can be found in Appendix A.
As already pointed out in the original Ref. [77], the model Ansätze for the nonperturbative low-energy interaction (“interaction form factors”) are such that they provide sufficient ultraviolet suppression. Therefore, as noted already in Ref. [77], no renormalization is needed and the multiplicative renormalization constants, which would otherwise be needed in the gap Equations (13) with (14), are 1. The usual expression for the condensate of the flavor q then becomes
q ¯ q = N c p Tr S q ( p ) N c T n Z d 3 p ( 2 π ) 3 Tr S q ( p , ω n ) ,
where Tr is the trace in Dirac space, and the combined integral-sum symbol indicates that when the calculation is at T > 0 , the four-momentum integration decomposes into the three-momentum integration and summation over fermionic Matsubara frequencies ω n = ( 2 n + 1 ) π T , n Z .
As is well known, the condensates (16) are finite only for massless quarks, m q = 0 , i.e., only q ¯ q 0 is finite, while the “massive” condensates are badly divergent, and must be regularized, i.e., divergences must be subtracted. Since the subtraction procedure is not uniquely defined, the chiral condensate at nonvanishing quark mass is also not uniquely defined. However, the arbitrariness is in practice slight and should rather be classified as fuzziness. It should not be given too large importance in the light of small differences between the results of various sensible procedures.
Our regularization procedure is subtracting the divergence-causing m q (∼ several MeV) from the scalar quark-dressing function B q ( p 2 ) (∼ several hundred MeV) whenever it is found in the numerator of the condensate integrand. To justify our particular regularization of massive condensates as physically meaningful and sensible, we have examined its consistency with two different subtractions used on lattice [81,82,83] and in a recent DSE-approach paper [84].
We shall now test our massive condensates obtained from the separable rank-2 DSE model (see Appendix A), whose regularized versions have already been shown in Figure 2.
Let us first consider the subtraction on lattice (normalized to 1 for T = 0 ) first proposed in Ref. [81] in their Equation (17), rewritten in our notation and applied to our condensate of u-quarks (and of course d-quarks in the isospin limit):
R ψ ¯ ψ ( T ) = R u ¯ u ( T ) = u ¯ u ( T ) u ¯ u ( 0 ) + q ¯ q 0 ( 0 ) q ¯ q 0 ( 0 ) .
In Figure 3, the upper, red curve shows (normalized) u-quark condensate u ¯ u ( T ) / u ¯ u ( 0 ) when regularized in the usual way, by subtracting m q from B q ( p 2 ) in the numerator of the condensate integrand. It agrees very well with the lattice regularization R u ¯ u (17) of our condensate u ¯ u ( T ) , represented by the green curve. The agreement with the lattice data points taken (if pertinent) from Table 6 of Ref. [82] is also rather good.
Next, we examine the consistency of our subtraction with the most usual condensate subtraction on the lattice, which combines the light and strange quark condensates and their masses like this:
Δ ¯ l , s ( T ) = l ¯ l l ( T ) m l m s s ¯ s s ( T ) .
Following Isserstedt et al. [84], in Figure 4 we make comparison of the normalized version thereof
Δ l , s ( T ) = l ¯ l ( T ) m l m s s ¯ s ( T ) l ¯ l ( 0 ) m l m s s ¯ s ( 0 ) ( l = u or d in the isospin symmetric limit )
with the lattice data of Ref. [83]. The agreement is very good, which implies also the agreement with the subtracted condensates in the recent DSE paper [84], which made this successful comparison first (in its Figure 3).
To conclude: results shown in Figure 3 and Figure 4 demonstrate that certain arbitrariness in the choice of regularization does not disqualify our massive condensates from useful applications, such as using them in Equation (9) to make predictions on the topological susceptibility.

3.3. Axion Mass and Topological Susceptibility—Results from the Rank-2 Separable Model in the Isosymmetric Limit

Our result for χ ( T ) 1 / 4 = m a ( T ) f a is presented in Figure 5 as a solid curve and compared, up to T 2.3 T c , with the corresponding results of two lattice groups [70,71], rescaled to the relative-temperature T / T c . (Table 1 gives numerical values of our results at T = 0 .)
In our case, the results for χ ( T ) and condensates u ¯ u ( T ) , d ¯ d ( T ) and s ¯ s ( T ) needed to obtain it, are predictions of the dynamical DSE model used in the T > 0 study of η - η [17]. This is the same modeling of the low-energy, nonperturbative QCD interactions as we have already employed in our earlier studies of light pseudoscalar mesons at T 0 [30,31,32,48]: the separable model interaction—see, e.g., [77,85], and references therein. We have adopted the so-called rank-2 variant from Ref. [77]. The adopted model with our choice of parameters is defined in detail in the Appendix A of the present work, after the subsection II.A of Refs. [31,86]. It employs the model current-quark-mass parameters m u = m d m l = 5.49 MeV and m s = 115 MeV. The model prediction for condensates at T = 0 are s ¯ s = ( 238.81 MeV ) 3 for the heaviest quark, while isosymmetric condensates of the lightest flavors, u ¯ u = d ¯ d l ¯ l = ( 218.69 MeV ) 3 are quite close to the “massless” one, q ¯ q 0 = ( 216.25 MeV ) 3 .
Contrary to, e.g., Ref. [48], where the condensate of massless quarks q ¯ q 0 ( T ) was used, in Ref. [17] and here we follow Shore [69] in using condensates of light quarks with nonvanishing current masses. The smooth, crossover behavior around the pseudocritical temperature T c for the chiral transition (now confirmed at vanishing baryon density by lattice studies such as [44,67,68,70,71]), is obtained thanks to the DChSB condensates of realistically massive light quarks—i.e., the quarks with realistic explicit chiral symmetry breaking [17].
In contrast, using in Equation (9) the massless-quark condensate q ¯ q 0 (which drops sharply to zero at T c ) instead of the “massive” ones, would dictate a sharp transition of the second order at T c [17,48] also for χ ( T ) , illustrated in Figure 5 by the dotted curve. This would of course imply that axions are massless [87] for T > T c . It is of academic interest to know what consequences would be thereof for cosmology, but now it is clear that only crossover is realistic [70,71].
The rather good agreement with lattice in Figure 5 resulted without any refitting of this model, either in Ref. [17] for η and η , or in this subsection. The model is in the isosymmetric limit, m u = m d m l , which is perfectly adequate for most purposes in hadronic physics. Nevertheless, the QCD topological susceptibility χ in its version (9) contains the current quark masses in the form of harmonic averages of m q q ¯ q ( q = u , d , s ). A harmonic average is dominated by its smallest argument, and presently this is the lightest current-quark-mass parameter, motivating us to investigate the changes occurring beyond the isospin symmetric point.

3.4. Axion Mass and Topological Susceptibility from Rank-1 and Rank-2 Models out of the Isosymmetric Limit

The previous isosymmetric case, pertinent also for the η - η study [17], has the current-quark-mass model parameters m u = m d = 5.49 MeV. This is above the most recent PDG quark-mass values [1], but anyway yields χ ( T = 0 ) = ( 72.73 MeV ) 4 already a little below the lattice results [70,71], and below the most recent chiral perturbation theory result χ ( T = 0 ) = ( 75.44 MeV ) 4 [76].
This seems not to bode well for the attempts out of the isosymmetric limit, because lowering the values of the current masses seems to threaten yielding unacceptably low values of the topological susceptibility. Indeed, taking the central values from the current quark masses m u = 2.2 0.4 + 0.5 MeV and m d = 4.70 0.3 + 0.5 MeV and m s = 95 3 + 9 MeV recently quoted by PDG [1], yields just ( 62.50 MeV ) 4 for the leading term of Equation (9) at T = 0 .
However, our model m u , m d and m s are phenomenological current-quark-mass parameters, and cannot be quite unambiguously and precisely related to the somewhat lower PDG values of the current quark masses. The better relation is through the ratios of quark masses, for which PDG gives m u / m d = 0.48 0.08 + 0.07 [1].
We thus require that m u fit / m d fit = 0.48 be satisfied by the new non-isosymmetric mass parameters m u fit and m d fit when they are varied to reproduce the recent most precise value χ ( T = 0 ) = ( 75.44 MeV ) 4 [76]. We get m u fit = 4.66 MeV, resulting in the condensate u ¯ u ( T = 0 ) = ( 218.35 MeV ) 3 and m d fit = 9.71 MeV, resulting in d ¯ d ( T = 0 ) = ( 220.33 MeV ) 3 . (The s-mass parameter is not varied, i.e., m s m s fit . The rest of model parameters, namely those in the Ansatz functions F 0 ( p 2 ) and F 1 ( p 2 ) modeling the strength of the rank-2 nonperturbative interaction (see Appendix), are also not varied.)
The T-dependence of the resulting χ ( T ) 1 / 4 is given by the short-dashed black curve in Figure 6. Except its better agreement with the lattice results [70,71] at low T, the new (dashed) χ ( T ) 1 / 4 curve is very close to the isosymmetric (solid) curve.
Now we will check the model dependence by comparing our results presented so far (obtained in the rank-2 model) with those we get in the separable rank-1 model of Ref. [77]. It is similar to the previously considered rank-2 one by modeling the low-energy, nonperturbative QCD interaction with an Ansatz separating the momenta p a , p b of interacting constituents, but is of a simpler form, proportional to just F 0 ( p a 2 ) F 0 ( p b 2 ) . Its presently interesting feature is that for similar quark-mass parameters, it yields significantly larger condensates than those in the separable rank-2 model, and thus also larger χ . (This also holds at low and vanishing T even for χ ( T ) calculated using only the “massless” condensate q ¯ q 0 ( T ) . This case is depicted in Figure 7 as the dotted curve.)
The original rank-1 model employs the light-quark current mass parameters in the isosymmetric limit: m u = m d m l = 6.6 MeV [77]. However, in Equation (9) we also need the s-flavor. The fit to the kaon mass yields m s = 142 MeV. The model prediction for condensates at T = 0 are then s ¯ s = ( 251.49 MeV ) 3 for the heaviest quark, while isosymmetric condensates of the lightest flavors, u ¯ u = d ¯ d l ¯ l = ( 249.27 MeV ) 3 are quite close to the “massless” one, q ¯ q 0 = ( 248.47 MeV ) 3 . This gives too large topological susceptibility at T = 0 , namely χ ( 0 ) = ( 84.08 MeV ) 4 . Nevertheless, for large T, it also falls with T somewhat faster than the rank-2 χ ( T ) , since rank-1 condensates fall with T somewhat faster than the rank-2 ones.
The isosymmetric rank-1 χ ( T ) is depicted by the solid black curve in Figure 7, showing that it actually falls with T faster even than χ ( T ) ’s from lattice [70,71] for practically all T’s high enough to induce changes. Then, comparing Figure 6 and Figure 7 shows that the lattice high-T results are in between high-T results of the two separable models.
To go out of the isospin limit with the rank-1 model, we again require that the changed parameters m u fit and m d fit (together with the condensates resulting from them) fit the currently most precise T = 0 value of the topological susceptibility, χ = ( 75.44 MeV ) 4 [76], while also obeying m u fit / m d fit = 0.48 , i.e., the central value of the PDG [1] m u / m d ratio. (Again, other model parameters including m s are not varied: m s m s fit .)
In our rank-1 model, these requirements yield m d fit = 6.56 MeV, i.e., it practically remained the same as in the originally fitted model [77]. Of course, its condensate d ¯ d also remains the same. The lightest flavor has the mass parameter lowered to m u fit = 3.15 MeV. Now it has only slightly lower condensate u ¯ u ( T = 0 ) = ( 248.87 MeV ) 3 , which is even closer to the “massless” q ¯ q 0 ( T = 0 ) . Nevertheless, u ¯ u ( T ) retains the crossover behavior for T > 0 , although it falls with T steeper than more “massive” condensates.
The resulting non-isosymmetric χ ( T ) 1 / 4 = m a ( T ) f a is in Figure 7 shown as the short-dashed black curve, which is everywhere consistently the lowest (among the “massive”, crossover curves).

4. Summary and Discussion

In the DSE framework, we have obtained predictions for the nontrivial part of the T-dependence of the axion mass m a ( T ) = χ ( T ) / f a , Equation (11), by calculating the QCD topological susceptibility χ ( T ) , since the unknown Peccei–Quinn scale f a is just an overall constant. We have used two empirically successful dynamical models of the separable type [77] to model nonperturbative QCD at T > 0 . We also studied the effects of varying the mass parameters of the lightest flavors out of the isospin limit, and found that our χ ( T ) , and consequently m a ( T ) , are robust with respect to the non-isosymmetric refitting of m u and m d = m u / 0.48 .
All these results of ours on χ ( T ) , and consequently the related axion mass, are in satisfactory agreement with the pertinent lattice results [70,71], and in qualitative agreement with those obtained in the NJL model [88]. Everyone obtains qualitatively similar crossover of χ ( T ) around T c , but it would be interesting to speculate what consequences for cosmology could be if χ ( T ) , and thus also m a ( T ) , would abruptly fall to zero at T = T c due to a sharp phase transition of the “massless” condensate q ¯ q 0 ( T ) . Of course, dynamical models of QCD can access only much smaller range of temperatures than lattice, where T 20 T c has already been reached [71]. (On the other hand, the thermal behavior of the U A ( 1 ) anomaly could not be accessed in chiral perturbation theory [89].)
Since it is now established that (at vanishing and low density) the chiral transition is a crossover, it is important that one can use massive-quark condensates, which exhibit crossover behavior around T T c . In the present work, they give us directly, through Equation (9), the crossover behavior of χ ( T ) . However, these are regularized condensates, because a nonvanishing current quark mass m q makes the condensate q ¯ q plagued by divergences, which must be subtracted. In Section 3, we have shown that our regularization procedure is reasonable and in good agreement with at least two widely used subtractions on the lattice.
To discuss our approach from a broader perspective, it is useful to recall that JLQCD collaboration [90] has recently pointed out how the chiral symmetry breaking and U A ( 1 ) anomaly are tied, and stressed the importance of the q q ¯ chiral condensate in that. The axion mass presently provides a simple example thereof: through Equation (11), m a ( T ) is at all temperatures directly expressed by the QCD topological susceptibility χ ( T ) , which is a measure of U A ( 1 ) breaking by the axial anomaly. We calculate χ ( T ) through Equation (9) from the quark condensates, which in turn arise from DChSB. In addition, conversely: melting of condensates around T T c signals the restoration of the chiral symmetry. Therefore, the U A ( 1 ) symmetry breaking and restoration being driven by the chiral ones is straightforward.
The relation of χ ( T ) to the η mass is, however, a little less straightforward [17] because it involves several other elements, but the topological susceptibility remains the main one. Since our present results on the axion are, in a way, a by-product of the framework which was initially formulated to understand better the T-dependence of η and η masses, we have explained it in detail in Section 2.2. Therefore, here in the Summary, we should just stress that the topological susceptibility χ ( T ) is the strong link between the QCD axion and the η - η complex. It so also in the case of the present paper regarding our η - η reference [17]: specifically, we should note that the actual T-dependence of η and η is rather sensitive to the behavior of χ ( T ) , and rather accurate χ ( T ) is needed to get acceptable M η ( T ) and M η ( T ) . Thanks to its crossover behavior, our χ ( T ) gives in Ref. [17] empirically allowed T-dependence of the masses in the η - η complex. However, even a crossover, if it were too steep, would lead to the unwanted (experimentally never seen) drop of the η mass, just as a too slow one would not yield the drop of the η mass required according to some experimental analyses [65,66].
In that sense, our present predictions on m a ( T ) are thus supported by the fact that our calculated topological susceptibility χ ( T ) gives the T-dependence of the U A ( 1 ) anomaly-influenced masses of η and η mesons [17] which is consistent with experimental evidence [65,66].

Author Contributions

Conceptualization, D.K. (Dubravko Klabučar); methodology, D.K. (Dubravko Klabučar), D.H. and D.K. (Dalibor Kekez); software, D.H.; validation, D.K. (Dubravko Klabučar), D.H. and D.K. (Dalibor Kekez); formal analysis, D.H., D.K. (Dubravko Klabučar) and D.K. (Dalibor Kekez); investigation, D.H., D.K. (Dubravko Klabučar) and D.K. (Dalibor Kekez); data curation, D.H.; writing–original draft preparation, D.K. (Dubravko Klabučar); writing–review and editing, D.K. (Dubravko Klabučar), D.K. (Dalibor Kekez) and D.H.; visualization, D.H., D.K. (Dalibor Kekez); supervision, D.K. (Dubravko Klabučar).

Funding

This work was supported in part by STSM grants from COST Actions CA15213 THOR and CA16214 PHAROS.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DChSBDynamical chiral symmetry breaking
DSEDyson–Schwinger equations
QCDQuantum chromodynamics
ABJAdler-Bell-Jackiw
CKMCabibbo-Kobayashi-Maskawa
CPcharge conjugation parity

Appendix A. Separable Interaction Models for Usage at T≥0

At T = 0 , the Dyson–Schwinger equation (DSE) approach in the rainbow-ladder approximation (RLA) tackles efficiently solving Dyson–Schwinger gap equation and Bethe–Salpeter equations, but extending this to T > 0 is technically quite difficult. We thus adopt a simple model for the strong dynamics from Ref. [77], namely the model we already used in Refs. [17,30,31,48]. For the effective gluon propagator in a Feynman-like gauge, we use the separable Ansatz:
g 2 D μ ν a b ( p ) eff = δ a b g 2 D μ ν eff ( p ) δ μ ν D ( p 2 , 2 , p · ) δ a b ,
whereby the dressed quark-propagator gap Equations (13) with (14) yields
B q ( p 2 ) = m q + 16 3 d 4 ( 2 π ) 4 D ( p 2 , 2 , p · ) B q ( 2 ) 2 A q 2 ( 2 ) + B q 2 ( 2 )
A q ( p 2 ) 1 p 2 = 8 3 d 4 ( 2 π ) 4 D ( p 2 , 2 , p · ) ( p · ) A q ( 2 ) 2 A q 2 ( 2 ) + B q 2 ( 2 ) .
More specifically, the so-called rank-2 separable interaction entails:
D ( p 2 , 2 , p · ) = D 0 F 0 ( p 2 ) F 0 ( 2 ) + D 1 F 1 ( p 2 ) ( p · ) F 1 ( 2 ) .
Then, the solutions of Equations (A2) and (A3) for the dressing functions are of the form
B q ( p 2 ) = m q + b q F 0 ( p 2 ) and A q ( p 2 ) = 1 + a q F 1 ( p 2 ) ,
reducing Equations (A2) and (A3) to the nonlinear system of equations for the constants b q and a q :
b q = 16 D 0 3 d 4 ( 2 π ) 4 F 0 ( p 2 ) B q ( 2 ) 2 A q 2 ( 2 ) + B q 2 ( 2 )
a q = 2 D 1 3 d 4 ( 2 π ) 4 2 F 1 ( 2 ) A q ( 2 ) 2 A q 2 ( 2 ) + B q 2 ( 2 ) .
If one chooses that the second term in the interaction (A4) is vanishing, by simply setting to zero the second strength constant, D 1 = 0 , one has a still simpler rank-1 separable Ansatz, where A q ( p 2 ) = 1 .
The analytic properties of these model interactions are defined by the choice of the interaction “form factors” F 0 ( p 2 ) and F 1 ( p 2 ) . In the present work we will use the functions [32,86]
F 0 ( p 2 ) = exp ( p 2 / Λ 0 2 ) and F 1 ( p 2 ) = 1 + exp ( p 0 2 / Λ 1 2 ) 1 + exp ( ( p 2 p 0 2 ) / Λ 1 2 ) ,
which satisfy the constraints F 0 ( 0 ) = F 1 ( 0 ) = 1 and F 0 ( ) = F 1 ( ) = 0 .
For the numerical calculations we fix the free parameters of the model at T = 0 as in Refs. [32,86], to reproduce in particular the vacuum masses of the pseudoscalar and vector mesons, M π = 140 MeV, M K = 495 MeV, M ρ = 770 MeV, the pion decay constant f π = 92 MeV, and decay widths, Γ ρ 0 e + e = 6.77 keV, Γ ρ π π = 151 MeV as basic requirements from low-energy QCD phenomenology.
We thus use the same parameter set as in Refs. [32,86], namely m u = m d = m l = 5.49 MeV, m s = 115 MeV, D 0 Λ 0 2 = 219 , D 1 Λ 0 4 = 40 , Λ 0 = 0.758 GeV, Λ 1 = 0.961 GeV and p 0 = 0.6 GeV for the rank-2 model.
For fixing the parameters in the rank-1 model, we use only the masses of pion and kaon, the pion decay constant, and GMOR as one additional constraint. This gives m u = m d = m l = 6.6 MeV, m s = 142 MeV, D 0 Λ 0 2 = 113.67 , and Λ 0 = 0.647 GeV for our values of the rank-1 parameters.
At T > 0 , p p n = ( ω n , p ) . Presently, pertinent are the fermion Matsubara frequencies ω n = ( 2 n + 1 ) π T . Due to loss of O ( 4 ) symmetry in medium, the dressed quark propagator (12) is at T > 0 replaced by
S q 1 ( p , ω n ; T ) = i γ · p A q ( p 2 , ω n ; T ) + i γ 4 ω n C q ( p 2 , ω n ; T ) + B q ( p 2 , ω n ; T ) .
For separable interactions, the dressing functions A q , C q and B q depend only on the sum p n 2 = ω n 2 + p 2 . In the separable models (A4), with their characteristic form (A5) of the propagator solutions at T = 0 , the dressing functions obtained as solutions of the gap equation at T > 0 are:
A q ( p n 2 ; T ) = 1 + a q ( T ) F 1 ( p n 2 ) , C q ( p n 2 ; T ) = 1 + c q ( T ) F 1 ( p n 2 ) , B q ( p n 2 ; T ) = m q + b q ( T ) F 0 ( p n 2 ) .
That is, the former gap constants a q and b q become temperature-dependent gap functions a f ( T ) , b f ( T ) and c f ( T ) obtained from the nonlinear system of equations:
a q ( T ) = 8 D 1 9 T n d 3 p ( 2 π ) 3 F 1 ( p n 2 ) p 2 A q ( p n 2 , T ) d q 1 ( p n 2 , T ) ,
c q ( T ) = 8 D 1 3 T n d 3 p ( 2 π ) 3 F 1 ( p n 2 ) ω n 2 C q ( p n 2 , T ) d q 1 ( p n 2 , T ) ,
b q ( T ) = 16 D 0 3 T n d 3 p ( 2 π ) 3 F 0 ( p n 2 ) B q ( p n 2 , T ) d q 1 ( p n 2 , T ) ,
where the denominator function is d q ( p n 2 , T ) = p 2 A q 2 ( p n 2 , T ) + ω n 2 C q 2 ( p n 2 , T ) + B q 2 ( p n 2 , T ) .

References

  1. Tanabashi, M.; Hagiwara, K.; Hikasa, K.; Nakamura, K.; Sumino, Y.; Takahashi, F.; Tanaka, J.; Agashe, K.; Aielli, G.; Amsler, C.; et al. [Particle Data Group]. Review of Particle Physics. Phys. Rev. 2018, D98, 030001. [Google Scholar] [CrossRef]
  2. Peccei, R.D.; Quinn, H.R. CP Conservation in the Presence of Instantons. Phys. Rev. Lett. 1977, 38, 1440–1443. [Google Scholar] [CrossRef]
  3. Peccei, R.D.; Quinn, H.R. Constraints Imposed by CP Conservation in the Presence of Instantons. Phys. Rev. 1977, D16, 1791–1797. [Google Scholar] [CrossRef]
  4. Weinberg, S. A New Light Boson? Phys. Rev. Lett. 1978, 40, 223–226. [Google Scholar] [CrossRef]
  5. Wilczek, F. Problem of Strong P and T Invariance in the Presence of Instantons. Phys. Rev. Lett. 1978, 40, 279–282. [Google Scholar] [CrossRef]
  6. Peccei, R.D. The Strong CP problem and axions. Lect. Notes Phys. 2008, 741, 3–17. [Google Scholar] [CrossRef]
  7. Cabibbo, N. Unitary Symmetry and Leptonic Decays. Phys. Rev. Lett. 1963, 10, 531–533. [Google Scholar] [CrossRef]
  8. Kobayashi, M.; Maskawa, T. CP Violation in the Renormalizable Theory of Weak Interaction. Prog. Theor. Phys. 1973, 49, 652–657. [Google Scholar] [CrossRef]
  9. Baker, C.A.; Doyle, D.D.; Geltenbort, P.; Green, K.; van der Grinten, M.G.D.; Harris, P.G.; Iaydjiev, P.; Ivanov, S.N.; May, D.J.R.; Pendlebury, J.M.; et al. An Improved experimental limit on the electric dipole moment of the neutron. Phys. Rev. Lett. 2006, 97, 131801. [Google Scholar] [CrossRef]
  10. Wantz, O.; Shellard, E.P.S. Axion Cosmology Revisited. Phys. Rev. 2010, D82, 123508. [Google Scholar] [CrossRef]
  11. Berkowitz, E.; Buchoff, M.I.; Rinaldi, E. Lattice QCD input for axion cosmology. Phys. Rev. 2015, D92, 034507. [Google Scholar] [CrossRef]
  12. Preskill, J.; Wise, M.B.; Wilczek, F. Cosmology of the Invisible Axion. Phys. Lett. 1983, 120B, 127–132. [Google Scholar] [CrossRef]
  13. Abbott, L.F.; Sikivie, P. A Cosmological Bound on the Invisible Axion. Phys. Lett. 1983, 120B, 133–136. [Google Scholar] [CrossRef]
  14. Dine, M.; Fischler, W. The Not So Harmless Axion. Phys. Lett. 1983, 120B, 137–141. [Google Scholar] [CrossRef]
  15. Kim, J.E.; Nam, S.; Semetzidis, Y.K. Fate of global symmetries in the Universe: QCD axion, quintessential axion and trans-Planckian inflaton decay-constant. Int. J. Mod. Phys. 2018, A33, 1830002. [Google Scholar] [CrossRef]
  16. Dvali, G.; Funcke, L. Domestic Axion. arXiv 2016, arXiv:1608.08969. [Google Scholar]
  17. Horvatić, D.; Kekez, D.; Klabučar, D. η and η mesons at high T when the UA(1) and chiral symmetry breaking are tied. Phys. Rev. 2019, D99, 014007. [Google Scholar] [CrossRef]
  18. Adler, S.L. Axial vector vertex in spinor electrodynamics. Phys. Rev. 1969, 177, 2426–2438. [Google Scholar] [CrossRef]
  19. Bell, J.S.; Jackiw, R. A PCAC puzzle: π0γγ in the σ model. Nuovo Cim. 1969, A60, 47–61. [Google Scholar] [CrossRef]
  20. Feldmann, T. Quark structure of pseudoscalar mesons. Int. J. Mod. Phys. 2000, A15, 159–207. [Google Scholar] [CrossRef]
  21. Kekez, D.; Klabučar, D.; Scadron, M.D. Revisiting the U(A)(1) problems. J. Phys. 2000, G26, 1335–1354. [Google Scholar] [CrossRef]
  22. Kekez, D.; Klabučar, D. Two photon processes of pseudoscalar mesons in a Bethe-Salpeter approach. Phys. Lett. 1996, B387, 14–20. [Google Scholar] [CrossRef]
  23. Kekez, D.; Bistrović, B.; Klabučar, D. Application of Jain and Munczek’s bound state approach to gamma gamma processes of Pi0, eta(c) and eta(b). Int. J. Mod. Phys. 1999, A14, 161–194. [Google Scholar] [CrossRef]
  24. Klabučar, D.; Kekez, D. eta and eta-prime in a coupled Schwinger-Dyson and Bethe-Salpeter approach. Phys. Rev. 1998, D58, 096003. [Google Scholar] [CrossRef]
  25. Kekez, D.; Klabučar, D. Gamma* gamma —> pi0 transition and asymptotics of gamma* gamma and gamma* gamma* transitions of other unflavored pseudoscalar mesons. Phys. Lett. 1999, B457, 359–367. [Google Scholar] [CrossRef]
  26. Bistrović, B.; Klabučar, D. Anomalous gamma —> 3 pi amplitude in a bound state approach. Phys. Lett. 2000, B478, 127–136. [Google Scholar] [CrossRef]
  27. Kekez, D.; Klabučar, D. Eta and eta-prime in a coupled Schwinger-Dyson and Bethe-Salpeter approach: The gamma* gamma transition form-factors. Phys. Rev. 2002, D65, 057901. [Google Scholar] [CrossRef]
  28. Kekez, D.; Klabučar, D. Pseudoscalar q anti-q mesons and effective QCD coupling enhanced (A**2) condensate. Phys. Rev. 2005, D71, 014004. [Google Scholar] [CrossRef]
  29. Kekez, D.; Klabučar, D. Eta and eta’ mesons and dimension 2 gluon condensate < A**2 >. Phys. Rev. 2006, D73, 036002. [Google Scholar] [CrossRef]
  30. Horvatić, D.; Blaschke, D.; Kalinovsky, Y.; Kekez, D.; Klabučar, D. eta and eta-prime mesons in the Dyson-Schwinger approach using a generalization of the Witten-Veneziano relation. Eur. Phys. J. 2008, A38, 257–264. [Google Scholar] [CrossRef]
  31. Horvatić, D.; Klabučar, D.; Radzhabov, A.E. eta and eta-prime mesons in the Dyson-Schwinger approach at finite temperature. Phys. Rev. 2007, D76, 096009. [Google Scholar] [CrossRef]
  32. Horvatić, D.; Blaschke, D.; Klabučar, D.; Radzhabov, A.E. Pseudoscalar Meson Nonet at Zero and Finite Temperature. Phys. Part. Nucl. 2008, 39, 1033–1039. [Google Scholar] [CrossRef]
  33. Alkofer, R.; von Smekal, L. The Infrared behavior of QCD Green’s functions: Confinement dynamical symmetry breaking, and hadrons as relativistic bound states. Phys. Rept. 2001, 353, 281. [Google Scholar] [CrossRef]
  34. Roberts, C.D.; Schmidt, S.M. Dyson-Schwinger equations: Density, temperature and continuum strong QCD. Prog. Part. Nucl. Phys. 2000, 45, S1–S103. [Google Scholar] [CrossRef]
  35. Holl, A.; Roberts, C.D.; Wright, S.V. Hadron physics and Dyson-Schwinger equations. In Proceedings of the 20th Annual Hampton University Graduate Studies Program (HUGS 2005), Newport News, VA, USA, 31 May–17 June 2005. [Google Scholar]
  36. Fischer, C.S. Infrared properties of QCD from Dyson-Schwinger equations. J. Phys. 2006, G32, R253–R291. [Google Scholar] [CrossRef]
  37. Benić, S.; Horvatić, D.; Kekez, D.; Klabučar, D. A UA(1) symmetry restoration scenario supported by the generalized Witten-Veneziano relation and its analytic solution. Phys. Lett. 2014, B738, 113–117. [Google Scholar] [CrossRef]
  38. Witten, E. Current Algebra Theorems for the U(1) Goldstone Boson. Nucl. Phys. 1979, B156, 269–283. [Google Scholar] [CrossRef]
  39. Veneziano, G. U(1) Without Instantons. Nucl. Phys. 1979, B159, 213–224. [Google Scholar] [CrossRef]
  40. Alles, B.; D’Elia, M.; Di Giacomo, A. Topological susceptibility at zero and finite T in SU(3) Yang-Mills theory. Nucl. Phys. 1997, B494, 281–292. [Google Scholar] [CrossRef]
  41. Boyd, G.; Engels, J.; Karsch, F.; Laermann, E.; Legeland, C.; Lutgemeier, M.; Petersson, B. Thermodynamics of SU(3) lattice gauge theory. Nucl. Phys. 1996, B469, 419–444. [Google Scholar] [CrossRef]
  42. Gattringer, C.; Hoffmann, R.; Schaefer, S. The Topological susceptibility of SU(3) gauge theory near T(c). Phys. Lett. 2002, B535, 358–362. [Google Scholar] [CrossRef]
  43. Di Vecchia, P.; Rossi, G.; Veneziano, G.; Yankielowicz, S. Spontaneous CP breaking in QCD and the axion potential: An effective Lagrangian approach. J. High Energy Phys. 2017, 12, 104. [Google Scholar] [CrossRef]
  44. Dick, V.; Karsch, F.; Laermann, E.; Mukherjee, S.; Sharma, S. Microscopic origin of UA(1) symmetry violation in the high temperature phase of QCD. Phys. Rev. 2015, D91, 094504. [Google Scholar] [CrossRef]
  45. Bazavov, A.; Ding, H.-T.; Hegde, P.; Kaczmarek, O.; Karsch, F.; Laermann, E.; Maezawa, Y.; Mukherjee, S.; Ohno, H.; Petreczky, P.; et al. The QCD Equation of State to O( μ B 6 ) from Lattice QCD. Phys. Rev. 2017, D95, 054504. [Google Scholar] [CrossRef]
  46. Bazavov, A. An overview of (selected) recent results in finite-temperature lattice QCD. J. Phys. Conf. Ser. 2013, 446, 012011. [Google Scholar] [CrossRef] [Green Version]
  47. Stoecker, H.; Zhou, K.; Schramm, S.; Senzel, F.; Greiner, C.; Beitel, M.; Gallmeister, K.; Gorenstein, M.; Mishustin, I.; Vasak, D.; et al. Glueballs amass at RHIC and LHC Colliders! - The early quarkless 1st order phase transition at T=270 MeV - from pure Yang-Mills glue plasma to GlueBall-Hagedorn states. J. Phys. 2016, G43, 015105. [Google Scholar] [CrossRef]
  48. Benić, S.; Horvatić, D.; Kekez, D.; Klabučar, D. eta-prime Multiplicity and the Witten-Veneziano relation at finite temperature. Phys. Rev. 2011, D84, 016006. [Google Scholar] [CrossRef]
  49. Leutwyler, H.; Smilga, A.V. Spectrum of Dirac operator and role of winding number in QCD. Phys. Rev. 1992, D46, 5607–5632. [Google Scholar] [CrossRef]
  50. Di Vecchia, P.; Veneziano, G. Chiral Dynamics in the Large n Limit. Nucl. Phys. 1980, B171, 253–272. [Google Scholar] [CrossRef]
  51. Dürr, S. Topological susceptibility in full QCD: Lattice results versus the prediction from the QCD partition function with granularity. Nucl. Phys. 2001, B611, 281–310. [Google Scholar] [CrossRef]
  52. Bernard, V.; Descotes-Genon, S.; Toucas, G. Topological susceptibility on the lattice and the three-flavour quark condensate. JHEP 2012, 06, 051. [Google Scholar] [CrossRef]
  53. Pisarski, R.D.; Wilczek, F. Remarks on the Chiral Phase Transition in Chromodynamics. Phys. Rev. 1984, D29, 338–341. [Google Scholar] [CrossRef]
  54. Höll, A.; Maris, P.; Roberts, C.D. Mean field exponents and small quark masses. Phys. Rev. 1999, C59, 1751–1755. [Google Scholar] [CrossRef]
  55. Kiriyama, O.; Maruyama, M.; Takagi, F. Current quark mass effects on chiral phase transition of QCD in the improved ladder approximation. Phys. Rev. 2001, D63, 116009. [Google Scholar] [CrossRef]
  56. Ikeda, T. Dressed quark propagator at finite temperature in the Schwinger-Dyson approach with the rainbow approximation: Exact numerical solutions and their physical implication. Prog. Theor. Phys. 2002, 107, 403–420. [Google Scholar] [CrossRef]
  57. Blank, M.; Krassnigg, A. The QCD chiral transition temperature in a Dyson-Schwinger-equation context. Phys. Rev. 2010, D82, 034006. [Google Scholar] [CrossRef]
  58. Fischer, C.S.; Mueller, J.A. Quark condensates and the deconfinement transition. PoS 2009, CPOD2009, 023. [Google Scholar] [CrossRef]
  59. Qin, S.x.; Chang, L.; Chen, H.; Liu, Y.x.; Roberts, C.D. Phase diagram and critical endpoint for strongly-interacting quarks. Phys. Rev. Lett. 2011, 106, 172301. [Google Scholar] [CrossRef]
  60. Gao, F.; Chen, J.; Liu, Y.X.; Qin, S.X.; Roberts, C.D.; Schmidt, S.M. Phase diagram and thermal properties of strong-interaction matter. Phys. Rev. 2016, D93, 094019. [Google Scholar] [CrossRef]
  61. Fischer, C.S. QCD at finite temperature and chemical potential from Dyson–Schwinger equations. Prog. Part. Nucl. Phys. 2019, 105, 1–60. [Google Scholar] [CrossRef]
  62. Ejiri, S.; Karsch, F.; Laermann, E.; Miao, C.; Mukherjee, S.; Petreczky, P.; Schmidt, C.; Soeldner, W.; Unger, W. On the magnetic equation of state in (2+1)-flavor QCD. Phys. Rev. 2009, D80, 094505. [Google Scholar] [CrossRef]
  63. Ding, H.T.; Hegde, P.; Karsch, F.; Lahiri, A.; Li, S.T.; Mukherjee, S.; Petreczky, P. Chiral phase transition of (2+1)-flavor QCD. Nucl. Phys. 2019, A982, 211–214. [Google Scholar] [CrossRef]
  64. Ding, H.T.; Hegde, P.; Kaczmarek, O.; Karsch, F.; Lahiri, A.; Li, S.T.; Mukherjee, S.; Petreczky, P. Chiral phase transition in (2 + 1)-flavor QCD. PoS 2019, LATTICE2018, 171. [Google Scholar] [CrossRef]
  65. Csorgo, T.; Vertesi, R.; Sziklai, J. Indirect observation of an in-medium η mass reduction in s N N = 200 GeV Au+Au collisions. Phys. Rev. Lett. 2010, 105, 182301. [Google Scholar] [CrossRef]
  66. Vertesi, R.; Csorgo, T.; Sziklai, J. Significant in-medium η ’ mass reduction in s N N = 200 GeV Au+Au collisions at the BNL Relativistic Heavy Ion Collider. Phys. Rev. 2011, C83, 054903. [Google Scholar] [CrossRef]
  67. Aoki, S.; Fukaya, H.; Taniguchi, Y. Chiral symmetry restoration, eigenvalue density of Dirac operator and axial U(1) anomaly at finite temperature. Phys. Rev. 2012, D86, 114512. [Google Scholar] [CrossRef]
  68. Buchoff, M.I.; Cheng, M.; Christ, N.; Ding, H.-T.; Jung, C.; Karsch, F.; Lin, Z.; Mawhinney, R.D.; Mukherjee, S.; Petreczky, P.; et al. QCD chiral transition, U(1)A symmetry and the dirac spectrum using domain wall fermions. Phys. Rev. 2014, D89, 054514. [Google Scholar] [CrossRef]
  69. Shore, G.M. Pseudoscalar meson decay constants and couplings, the Witten-Veneziano formula beyond large N(c), and the topological susceptibility. Nucl. Phys. 2006, B744, 34–58. [Google Scholar] [CrossRef]
  70. Petreczky, P.; Schadler, H.P.; Sharma, S. The topological susceptibility in finite temperature QCD and axion cosmology. Phys. Lett. 2016, B762, 498–505. [Google Scholar] [CrossRef]
  71. Borsanyi, S.; Fodor, Z.; Kampert, K.H.; Katz, S.D.; Kawanai, T.; Kovacs, T.G.; Mages, S.W.; Pasztor, A.; Pittler, F.; Redondo, J.; et al. Calculation of the axion mass based on high-temperature lattice quantum chromodynamics. Nature 2016, 539, 69–71. [Google Scholar] [CrossRef] [Green Version]
  72. Bonati, C.; D’Elia, M.; Mariti, M.; Martinelli, G.; Mesiti, M.; Negro, F.; Sanfilippo, F.; Villadoro, G. Axion phenomenology and θ-dependence from Nf = 2 + 1 lattice QCD. J. High Energy Phys. 2016, 03, 155. [Google Scholar] [CrossRef]
  73. Takahashi, F.; Yin, W.; Guth, A.H. QCD axion window and low-scale inflation. Phys. Rev. 2018, D98, 015042. [Google Scholar] [CrossRef]
  74. Bonati, C.; D’Elia, M.; Martinelli, G.; Negro, F.; Sanfilippo, F.; Todaro, A. Topology in full QCD at high temperature: A multicanonical approach. J. High Energy Phys. 2018, 11, 170. [Google Scholar] [CrossRef]
  75. Grilli di Cortona, G.; Hardy, E.; Pardo Vega, J.; Villadoro, G. The QCD axion, precisely. J. High Energy Phys. 2016, 01, 034. [Google Scholar] [CrossRef]
  76. Gorghetto, M.; Villadoro, G. Topological Susceptibility and QCD Axion Mass: QED and NNLO corrections. J. High Energy Phys. 2019, 03, 033. [Google Scholar] [CrossRef]
  77. Blaschke, D.; Burau, G.; Kalinovsky, Y.L.; Maris, P.; Tandy, P.C. Finite T meson correlations and quark deconfinement. Int. J. Mod. Phys. 2001, A16, 2267–2291. [Google Scholar] [CrossRef]
  78. Alkofer, R.; Watson, P.; Weigel, H. Mesons in a Poincare covariant Bethe-Salpeter approach. Phys. Rev. 2002, D65, 094026. [Google Scholar] [CrossRef]
  79. Contant, R.; Huber, M.Q. Phase structure and propagators at nonvanishing temperature for QCD and QCD-like theories. Phys. Rev. 2017, D96, 074002. [Google Scholar] [CrossRef]
  80. Blaschke, D.; Kalinovsky, Y.L.; Radzhabov, A.E.; Volkov, M.K. Scalar sigma meson at a finite temperature in a nonlocal quark model. Phys. Part. Nucl. Lett. 2006, 3, 327–330. [Google Scholar] [CrossRef]
  81. Burger, F.; Ilgenfritz, E.M.; Kirchner, M.; Lombardo, M.P.; Muller-Preussker, M.; Philipsen, O.; Urbach, C.; Zeidlewicz, L. Thermal QCD transition with two flavors of twisted mass fermions. Phys. Rev. 2013, D87, 074508. [Google Scholar] [CrossRef]
  82. Kotov, A.Y.; Lombardo, M.P.; Trunin, A.M. Fate of the η in the quark gluon plasma. Phys. Lett. 2019, B794, 83–88. [Google Scholar] [CrossRef]
  83. Borsanyi, S.; Fodor, Z.; Hoelbling, C.; Katz, S.D.; Krieg, S.; Ratti, C.; Szabo, K.K. Is there still any Tc mystery in lattice QCD? Results with physical masses in the continuum limit III. J. High Energy Phys. 2010, 09, 073. [Google Scholar] [CrossRef]
  84. Isserstedt, P.; Buballa, M.; Fischer, C.S.; Gunkel, P.J. Baryon number fluctuations in the QCD phase diagram from Dyson-Schwinger equations. arXiv 2019, arXiv:1906.11644. [Google Scholar] [CrossRef]
  85. Cahill, R.T.; Gunner, S.M. The Global color model of QCD for hadronic processes: A Review. Fizika 1998, B7, 171. [Google Scholar]
  86. Blaschke, D.; Horvatić, D.; Klabučar, D.; Radzhabov, A.E. Separable Dyson-Schwinger model at zero and finite T. In Proceedings of the Mini-Workshop Bled 2006 on Progress in Quark Models, Bled, Slovenia, 10–17 July 2006; pp. 20–26. [Google Scholar]
  87. Klabučar, D.; Horvatić, D.; Kekez, D. T-dependence of the axion mass when the UA(1) and chiral symmetry breaking are tied. In Proceedings of the Workshop on Excited QCD 2019, Schladming, Austria, 30 January–3 February 2019. [Google Scholar]
  88. Lu, Z.Y.; Ruggieri, M. Effect of the chiral phase transition on axion mass and self-coupling. Phys. Rev. 2019, D100, 014013. [Google Scholar] [CrossRef]
  89. Gu, X.W.; Duan, C.G.; Guo, Z.H. Updated study of the η-η mixing and the thermal properties of light pseudoscalar mesons at low temperatures. Phys. Rev. 2018, D98, 034007. [Google Scholar] [CrossRef]
  90. Fukaya, H. Can axial U(1) anomaly disappear at high temperature? EPJ Web Conf. 2018, 175, 01012. [Google Scholar] [CrossRef]
1.
The mass eigenstate η is approximated only roughly by the pure S U ( 3 ) singlet state η 0 , due to the relatively large explicit breaking of the flavor S U ( 3 ) symmetry by much heavier s-quark: 2 m s / ( m u + m d ) = 27.3 ± 0.7 [1].
2.
ABJ anomaly stands for names of Adler, Bell, and Jackiw, as a reminder of their pioneering work on anomalies [18,19] exactly half a century ago this year.
3.
This is in contrast with χ = χ QCD , the full-QCD topological susceptibility (3), which is much harder to find on the lattice because of the light quark flavors. χ = χ QCD approaches χ YM only if one takes the quenched limit of infinitely massive quarks, χ YM = χ quench , since quarks then disappear from the loops of Equation (3).
4.
The combinations P q q ¯ need not always pertain to physical mesons. The pseudoscalar hidden-flavor states u u ¯ , d d ¯ , s s ¯ are not physical as long as the U A ( 1 ) symmetry is not restored ( i.e., the anomaly effectively turned off, see around Equation (2.6) in Ref. [43] for example), but build the S U ( 3 ) states η 0 , η 8 and π 0 .
5.
This is a feature exhibited by DSE models, or at least by most of them, through the characteristic drop of their chiral-limit, massless q q ¯ condensates [54,55,56,57,58,59,60,61].
Figure 1. U A ( 1 ) anomaly-induced, hidden-flavor-mixing transitions from pseudoscalar quark–antiquark states P = q q ¯ to P = q q ¯ include both possibilities q = q and q q . Springs symbolize gluons. All lines and vertices are dressed in accord with the nonperturbative QCD. Nonperturbative configurations are essential for nonvanishing anomalous mass [20] contribution to η 0 η , since Q ( x ) is a total divergence. The gray blob symbolizes the infinity of all intermediate gluon states enabling such transitions, so that the three bold dots represent any even [21] number of additional gluons. Just one of infinitely many, but certainly the simplest realization thereof, is when such a transition is mediated by just two gluons (and no additional intermediate states), whereby the above figure reduces to the so-called “diamond graph”.
Figure 1. U A ( 1 ) anomaly-induced, hidden-flavor-mixing transitions from pseudoscalar quark–antiquark states P = q q ¯ to P = q q ¯ include both possibilities q = q and q q . Springs symbolize gluons. All lines and vertices are dressed in accord with the nonperturbative QCD. Nonperturbative configurations are essential for nonvanishing anomalous mass [20] contribution to η 0 η , since Q ( x ) is a total divergence. The gray blob symbolizes the infinity of all intermediate gluon states enabling such transitions, so that the three bold dots represent any even [21] number of additional gluons. Just one of infinitely many, but certainly the simplest realization thereof, is when such a transition is mediated by just two gluons (and no additional intermediate states), whereby the above figure reduces to the so-called “diamond graph”.
Universe 05 00208 g001
Figure 2. The relative-temperature T / T c dependences (where T Ch T c ) of (the 3rd root of the absolute value of) the q q ¯ condensates, and of (the 4th root of) the topological susceptibility χ ( T ) and the full QCD topological charge parameter A ( T ) . Everything was calculated in the isosymmetric limit using the separable rank-2 DSE model which we had already used in Refs. [17,31,48]. (See Appendix A for the model interaction form and parameters, and the first line of Table 1 for the numerical values of the condensates and χ at T = 0 .) Only the chiral-limit condensate q ¯ q 0 ( T ) falls steeply to zero at T = T Ch , indicative of the second-order phase transition. This sharp transition would through (7) and (6) be transmitted to, respectively, the chiral-limit χ ( T ) and χ ˜ ( T ) , and ultimately to the η and η masses in Ref. [48]. The highest curve (dash-dotted) and the second one from above (dashed) are (3rd roots of the absolute values of) the condensates s ¯ s ( T ) and u ¯ u ( T ) , respectively. Their smooth crossover behaviors carry over to χ ( T ) and A ( T ) through Equations (9) and (8), respectively, leading to the empirically acceptable predictions [17] for the T-dependence of the η and η masses. It turns out that χ ( T ) (9) also gives the smooth crossover behavior also to the T-dependence of the axion mass.
Figure 2. The relative-temperature T / T c dependences (where T Ch T c ) of (the 3rd root of the absolute value of) the q q ¯ condensates, and of (the 4th root of) the topological susceptibility χ ( T ) and the full QCD topological charge parameter A ( T ) . Everything was calculated in the isosymmetric limit using the separable rank-2 DSE model which we had already used in Refs. [17,31,48]. (See Appendix A for the model interaction form and parameters, and the first line of Table 1 for the numerical values of the condensates and χ at T = 0 .) Only the chiral-limit condensate q ¯ q 0 ( T ) falls steeply to zero at T = T Ch , indicative of the second-order phase transition. This sharp transition would through (7) and (6) be transmitted to, respectively, the chiral-limit χ ( T ) and χ ˜ ( T ) , and ultimately to the η and η masses in Ref. [48]. The highest curve (dash-dotted) and the second one from above (dashed) are (3rd roots of the absolute values of) the condensates s ¯ s ( T ) and u ¯ u ( T ) , respectively. Their smooth crossover behaviors carry over to χ ( T ) and A ( T ) through Equations (9) and (8), respectively, leading to the empirically acceptable predictions [17] for the T-dependence of the η and η masses. It turns out that χ ( T ) (9) also gives the smooth crossover behavior also to the T-dependence of the axion mass.
Universe 05 00208 g002
Figure 3. The relative-temperature T / T c dependence of the subtracted (and normalized) condensate R ψ ¯ ψ defined by Equation (17) and introduced by Ref. [81]. The lattice data points are from Figure 6 of Ref. [82], but scaled for the critical temperatures T χ from their Table 2, which is different for the “crosses” (data points [82] for m π 370 MeV) and “bars” (data points [82] for m π 210 MeV). The lower, green curve results from the R ψ ¯ ψ (17) subtraction of our u-quark condensate. The upper, red curve is the T-dependence of our u-quark condensate when regularized in the usual way (see text).
Figure 3. The relative-temperature T / T c dependence of the subtracted (and normalized) condensate R ψ ¯ ψ defined by Equation (17) and introduced by Ref. [81]. The lattice data points are from Figure 6 of Ref. [82], but scaled for the critical temperatures T χ from their Table 2, which is different for the “crosses” (data points [82] for m π 370 MeV) and “bars” (data points [82] for m π 210 MeV). The lower, green curve results from the R ψ ¯ ψ (17) subtraction of our u-quark condensate. The upper, red curve is the T-dependence of our u-quark condensate when regularized in the usual way (see text).
Universe 05 00208 g003
Figure 4. The relative-temperature T / T c dependence of the (normalized) subtracted quark condensate (19) from the lattice [83] (blue squares) and from our condensates. Slightly lower, green curve results from our unsubtracted condensates plugged in Equation (19), while the very slightly higher, red curve is from our already subtracted condensates.
Figure 4. The relative-temperature T / T c dependence of the (normalized) subtracted quark condensate (19) from the lattice [83] (blue squares) and from our condensates. Slightly lower, green curve results from our unsubtracted condensates plugged in Equation (19), while the very slightly higher, red curve is from our already subtracted condensates.
Universe 05 00208 g004
Figure 5. The relative-temperature T / T c dependence of (the leading term of) χ ( T ) 1 / 4 from our often adopted [17,30,31,32,48] isosymmetric DSE rank-2 separable model: solid curve for Equation (9) with massive-quark condensates, while the dotted curve results from using q ¯ q 0 instead. χ ( T ) 1 / 4 (with uncertainties) from lattice: dash-dotted curve extracted from Petreczky et al. [70] and long-dashed curve, from Borsany et al. [71]. (Colors online).
Figure 5. The relative-temperature T / T c dependence of (the leading term of) χ ( T ) 1 / 4 from our often adopted [17,30,31,32,48] isosymmetric DSE rank-2 separable model: solid curve for Equation (9) with massive-quark condensates, while the dotted curve results from using q ¯ q 0 instead. χ ( T ) 1 / 4 (with uncertainties) from lattice: dash-dotted curve extracted from Petreczky et al. [70] and long-dashed curve, from Borsany et al. [71]. (Colors online).
Universe 05 00208 g005
Figure 6. The short-dashed black curve shows the non-isosymmetric case of the leading term of χ ( T ) 1 / 4 , Equation (9), with m u fit = 4.66 MeV and m d fit = m u fit / 0.48 = 9.71 MeV, and appropriately recalculated condensates u ¯ u ( T ) and d ¯ d ( T ) . The vertical scale is zoomed with respect to Figure 5 to help resolve the short-dashed curve from the solid curve representing again the isosymmetric case of the same separable rank-2 model. Also, the lattice results [70,71] are again depicted as in Figure 5.
Figure 6. The short-dashed black curve shows the non-isosymmetric case of the leading term of χ ( T ) 1 / 4 , Equation (9), with m u fit = 4.66 MeV and m d fit = m u fit / 0.48 = 9.71 MeV, and appropriately recalculated condensates u ¯ u ( T ) and d ¯ d ( T ) . The vertical scale is zoomed with respect to Figure 5 to help resolve the short-dashed curve from the solid curve representing again the isosymmetric case of the same separable rank-2 model. Also, the lattice results [70,71] are again depicted as in Figure 5.
Universe 05 00208 g006
Figure 7. From the calculation in the separable rank-1 DSE model [77], the relative-temperature T / T c dependence of (the leading term of) χ ( T ) 1 / 4 is represented by: (i) the solid curve for the isosymmetric case with m u = m d = 6.6 MeV and m s = 142 MeV, (ii) the dotted curve is for the same mass parameters, but with all condensates approximated by the “massless” condensate q ¯ q 0 ( T ) , and (iii) the short-dashed curve for the non-isosymmetric case m u fit = 3.15 MeV, m d fit = 6.56 MeV, while m s fit m s = 142 MeV. The pertinent lattice results are presented in the same way as in the previous two Figures: the dash-dotted and long-dashed curves extracted, respectively, from Refs. [70,71]. Colors online.
Figure 7. From the calculation in the separable rank-1 DSE model [77], the relative-temperature T / T c dependence of (the leading term of) χ ( T ) 1 / 4 is represented by: (i) the solid curve for the isosymmetric case with m u = m d = 6.6 MeV and m s = 142 MeV, (ii) the dotted curve is for the same mass parameters, but with all condensates approximated by the “massless” condensate q ¯ q 0 ( T ) , and (iii) the short-dashed curve for the non-isosymmetric case m u fit = 3.15 MeV, m d fit = 6.56 MeV, while m s fit m s = 142 MeV. The pertinent lattice results are presented in the same way as in the previous two Figures: the dash-dotted and long-dashed curves extracted, respectively, from Refs. [70,71]. Colors online.
Universe 05 00208 g007
Table 1. For the both variants of the DSE separable model (with the rank-2 and rank-1 interaction Ansatz) used in the present paper, various sets of values of the model quark-mass parameters m q ( q = u , d , s ) are related to the model results for the topological susceptibility χ and the “massive" condensates q ¯ q at T = 0 . The topological susceptibility χ varies because of varying m q and (to a much lesser extent) because of the changes of “massive" condensates induced by changes of the quark-mass parameters m q . The massless-quark condensate, q ¯ q 0 , depends only on the dynamical DSE model: always q ¯ q 0 = 216.25 3 MeV 3 for the rank-2 model, and q ¯ q 0 = 248.47 3 MeV 3 for the rank-1 model. Thus, the topological susceptibility χ 0 calculated with the chiral-limit condensate, varies for a given model only because of varying values of m q . All values are in MeV (or the indicated 3rd or 4th powers of MeV).
Table 1. For the both variants of the DSE separable model (with the rank-2 and rank-1 interaction Ansatz) used in the present paper, various sets of values of the model quark-mass parameters m q ( q = u , d , s ) are related to the model results for the topological susceptibility χ and the “massive" condensates q ¯ q at T = 0 . The topological susceptibility χ varies because of varying m q and (to a much lesser extent) because of the changes of “massive" condensates induced by changes of the quark-mass parameters m q . The massless-quark condensate, q ¯ q 0 , depends only on the dynamical DSE model: always q ¯ q 0 = 216.25 3 MeV 3 for the rank-2 model, and q ¯ q 0 = 248.47 3 MeV 3 for the rank-1 model. Thus, the topological susceptibility χ 0 calculated with the chiral-limit condensate, varies for a given model only because of varying values of m q . All values are in MeV (or the indicated 3rd or 4th powers of MeV).
T = 0 m u m d m s χ 0 u ¯ u d ¯ d s ¯ s χ
(with q q ¯ 0 )
rank-2
m u = m d [17]5.495.49115 72.18 4 218.69 3 218.69 3 238.81 3 72.73 4
with constraint
m u = 0.48 m d [1],4.669.71115 74.64 4 218.35 3 220.33 3 238.81 3 75.44 4
fitted m u & m d
rank-1
m u = m d 6.66.6142 83.87 4 249.27 3 249.27 3 251.49 3 84.08 4
with constraint
m u = 0.48 m d [1],3.156.56142 75.31 4 248.87 3 249.21 3 251.49 3 75.43 4
fitted m u & m d

Share and Cite

MDPI and ACS Style

Horvatić, D.; Kekez, D.; Klabučar, D. Temperature Dependence of the Axion Mass in a Scenario Where the Restoration of Chiral Symmetry Drives the Restoration of the UA(1) Symmetry. Universe 2019, 5, 208. https://doi.org/10.3390/universe5100208

AMA Style

Horvatić D, Kekez D, Klabučar D. Temperature Dependence of the Axion Mass in a Scenario Where the Restoration of Chiral Symmetry Drives the Restoration of the UA(1) Symmetry. Universe. 2019; 5(10):208. https://doi.org/10.3390/universe5100208

Chicago/Turabian Style

Horvatić, Davor, Dalibor Kekez, and Dubravko Klabučar. 2019. "Temperature Dependence of the Axion Mass in a Scenario Where the Restoration of Chiral Symmetry Drives the Restoration of the UA(1) Symmetry" Universe 5, no. 10: 208. https://doi.org/10.3390/universe5100208

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop