{'fraction_non_alphanumeric': 0.05609033691225472, 'fraction_numerical': 0.016882636060718253, 'mean_word_length': 4.090652091971354, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 15, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
\eta' Multiplicity and the Witten-Veneziano relation at finite temperature

We discuss and propose the minimal generalization of the Witten-Veneziano relation to finite temperatures, prompted by STAR and PHENIX experimental results on the multiplicity of eta' mesons. After explaining why these results show that the zero-temperature Witten-Veneziano relation cannot be straightforwardly extended to temperatures T too close to the chiral restoration temperature T_Ch and beyond, we find the quantity which should replace, at T>0, the Yang-Mills topological susceptibility appearing in the T=0 Witten-Veneziano relation, in order to avoid the conflict with experiment at T>0. This is illustrated through concrete T-dependences of pseudoscalar meson masses in a chirally well-behaved, Dyson-Schwinger approach, but our results and conclusions are of a more general nature and, essentially, model-independent.

We discuss and propose the minimal generalization of the Witten-Veneziano relation to finite temperatures, prompted by STAR and PHENIX experimental results on the multiplicity of η ′ mesons. After explaining why these results show that the zero-temperature Witten-Veneziano relation cannot be straightforwardly extended to temperatures T too close to the chiral restoration temperature T Ch and beyond, we find the quantity which should replace, at T > 0, the Yang-Mills topological susceptibility appearing in the T = 0 Witten-Veneziano relation, in order to avoid the conflict with experiment at T > 0. This is illustrated through concrete T -dependences of pseudoscalar meson masses in a chirally well-behaved, Dyson-Schwinger approach, but our results and conclusions are of a more general nature and, essentially, model-independent.

I. INTRODUCTION
The ultrarelativistic heavy-ion collider facilities like RHIC at BNL and LHC at CERN strive to produce a new form of hot QCD matter. The experiments show [1,2] that it has very intricate properties and presents a big challenge especially for theoretical understanding. While above the (pseudo)critical temperature T c ∼ 170 MeV this matter is often called the quark-gluon plasma (QGP), it cannot be a perturbatively interacting quarkgluon gas (as widely expected before RHIC results [1,2]) until significantly higher temperatures T ≫ T c . Instead, the interactions and correlations in the hot QCD matter are still strong (e.g., see Refs. [3,4]) so that its more recent and more precise name is strongly coupled QGP (sQGP) [4]. One of its peculiarities seems to be that strong correlations in the form of quark-antiquark (qq) bound states and resonances still exist [3,5] in the sQGP well above T c . In the old QGP paradigm, even deeply bound charmonium (cc) states such as J/Ψ and η c were expected to unbind at T ≈ T c , but lattice QCD simulations of mesonic correlators now indicate they persist till around 2T c [6,7] or even above [8]. Similar indications for light-quark mesonic bound states are also accumulating from lattice QCD [9] and from other methods [3,10,11]. This agrees well with the findings on the lattice (e.g., see Ref. [12] for a review) that for realistic explicit chiral symmetry breaking (ChSB), i.e., for the physical values of the current quark masses, the transition between the hadron phase and the phase dominated by quarks and gluons, is not an abrupt, singular phase transition but a smooth, analytic crossover around the pseudocritical temperature T c . It is thus not too surprising that a clear experimental signal of, e.g., deconfinement, is still hard to find and identify unambiguously.
The most compelling signal for production of a new * Electronic address: klabucar@oberon.phy.hr; Corresponding author; Senior Associate of Abdus Salam ICTP, Trieste, Italy form of QCD matter, i.e., sQGP, would be a restoration -in hot and/or dense matter -of the symmetries of the QCD Lagrangian which are broken in the vacuum. One of them is the [SU A (N f ) flavor] chiral symmetry, whose dynamical breaking results in light, (almost-)Goldstone pseudoscalar (P ) mesons -namely the octet P = π 0 , π ± , K 0 ,K 0 , K ± , η, as we consider all three lightquark flavors, N f = 3. The second one is the U A (1) symmetry. Its breaking by the non-Abelian axial Adler-Bell-Jackiw anomaly ('gluon anomaly' for short) makes the remaining pseudoscalar meson of the light-quark sector, the η ′ , much heavier, preventing its appearance as the ninth (almost-)Goldstone boson of dynamical chiral symmetry breaking (DChSB) in QCD.
The first experimental signature of a partial restoration of the U A (1) symmetry seems to have been found in the √ s NN = 200 GeV central Au+Au reactions at RHIC. Namely, Csörgő et al. [13,14] analyzed combined data of PHENIX [15] and STAR [16] collaborations very robustly, through six popular models for hadron multiplicities, and found that at 99.9% confidence level, the η ′ mass, which in the vacuum is M η ′ = 957.8 MeV, is reduced by at least 200 MeV inside the fireball. It is the sign of the disappearing contribution of the gluon axial anomaly to the η ′ mass, which would drop to a value readily understood together with the (flavor-symmetrybroken) octet of qq ′ (q, q ′ = u, d, s) pseudoscalar mesons. This is the issue of the "return of the prodigal Goldstone boson" predicted [17] as a signal of the U A (1) symmetry restoration.
Nevertheless, the T -dependence of χ YM is such [24] that the straightforward extension of Eq. (1) to T > 0 [24], i.e., replacement of all quantities 1 therein by their respective T -dependent versions M η ′ (T ), M η (T ), M K (T ), f π (T ) and χ YM (T ), leads to a conflict with experiment [13,14]. Since this extension of Eq. (1) to T > 0 was studied in Ref. [24] before the pertinent experimental analysis [13,14], one of the purposes of this paper is to revisit the implications of the results of Ref. [24] for WVR at T > 0, and demonstrate explicitly that they are practically model-independent. The other, more important purpose is to propose a mechanism which can enable WVR to agree with experiment at T > 0.

II. THE RELATIONS CONNECTING TWO THEORIES, QCD AND YM
Both issues pointed out before Eq. (1) and around it, are best understood in a model-independent way if one starts from the chiral limit of vanishing current quark masses (m q = 0) for all three light flavors, q = u, d, s. Then not only pions and kaons are massless, but is also η, which is then (since the situation is also SU(3)-flavor-symmetric) a purely SU(3)-octet state, η = η 8 . In contrast, η ′ is then purely singlet, η ′ = η 0 ; since the divergence of the singlet axial quark currentqγ µ γ 5 1 2 λ 0 q is nonvanishing even for m q = 0 due to the gluon anomaly, the η ′ mass squared receives the anomalous contribution ∆M 2 η ′ (= λ 4 /f 2 η ′ in the notation of Ref. [17]) which is nonvanishing even in the chiral limit: However, λ 4 and f η ′ are known accurately 2 only in the large N c limit. There, in the leading order in 1/N c , λ 4 is given by the YM (i.e., "pure glue") topological susceptibility χ YM times 2N f = 6 [18,19], and the "η ′ decay constant" f η ′ is the same as f π [27]. Thus, keeping only the leading order in 1/N c , the last equality is WVR in the chiral limit. The consequences of Eq. (2) remain qualitatively the same realistically away from the chiral limit. This will soon become clear on the basis of, e.g., Eq. (3) below. Namely, due to DChSB in QCD, for relatively light current quark masses m q (q = u, d, s), the qq ′ bound-state pseudoscalar meson masses (including the nonanomalous 1 Throughout this paper, all quantities are for definiteness assumed at T = 0 unless their T -dependence is specifically indicated in formulas or in the text. 2 Also note that a unique "η ′ decay constant" f η ′ is, strictly speaking, not a well-defined quantity, as two η ′ decay constants are actually needed: the singlet one, f 0 η ′ , and the octet one, f 8 η ′ ; e.g., see an extensive review [25] or the short Appendix of Ref. [26]. parts of the η ′ and η masses) behave as The pseudoscalar mesons (including η ′ ) thus obtain relatively light nonanomalous contributions M qq ′ to their masses M P , allowing them to reach the empirical values. That is, instead of the eight strictly massless Goldstone bosons, π 0 , π ± , K 0 ,K 0 , K ± and η are relatively light almost-Goldstones. Among them, in the limit of isospin symmetry (m u = m d ), only η now receives also the gluon-anomaly contribution since the explicit SU (3) flavor breaking between the nonstrange (N S) u, d-quarks and s-quarks causes the mixing between the isoscalars η and η ′ . For m q = 0, Eq. (2) is replaced by the usual WVR (1) containing also the nonanomalous contributions to meson masses. Nevertheless, these contributions largely cancel due to the approximate SU(3) flavor symmetry and to DChSB [i.e., Eq. (3)]. This can be seen assuming the usual SU(3) qq content of the pseudoscalar meson nonet with well-defined isospin 3 quantum numbers, in particular the isoscalar (I = 0) octet and singlet etas, 3, whose mixing yields the physical particles η and η ′ . Since the nonanomalous parts of the η 0 and η 8 masses squared, , the nonanomalous parts of the η and η ′ masses are canceled by 2M 2 K in WVR (1). Another way of seeing this is expressing the nonanomalous parts of , showing again that already WVR's chiral-limitnonvanishing part (2) reveals the essence of the influence of the gluon anomaly on the masses in the η ′ -η complex. This is important also for the presently pertinent finite-T context because thanks to this, below it will be shown model-independently that WVR (1) containing the YM topological susceptibility χ YM implies T -dependence of η ′ mass in conflict with the recent experimental results [13,14].
Namely, the gluon-anomaly contribution (2) is established at T = 0 but it is not expected to persist at high temperatures. Ultimately, η ′ should also become a massless Goldstone boson at sufficiently high T , where χ YM (T ) → 0. However, according to WVR, ∆M η ′ (T ) falls only for T where f π (T ) 2 does not fall faster than 6χ YM (T ), as stressed in Ref. [24].
The WVR's chiral-limit version (2) manifestly points out the ratio χ YM (T )/f π (T ) 2 as crucial for the anomalous η ′ mass, but the above discussion shows that this remains essentially the same away from the chiral limit.
In the present context, it is important for practical calculations to go realistically away from the chiral limit, in which the chiral restoration is a sharp phase transition at its critical temperature T Ch where the chiral-limit pion decay constant vanishes very steeply, i.e., as steeply as the chiral quark condensate. In contrast, for realistic explicit ChSB, i.e., m u and m d of several MeV, this transition is a smooth crossover (e.g., see Ref. [12]). For the pion decay constant, this implies that f π (T ) still falls relatively steeply around pseudocritical temperature T Ch , but less so than in the chiral case, and even remains finite, enabling the usage of WVR (1) for the temperatures across the chiral and U A (1) symmetry restorations.
WVR is very remarkable because it connects two different theories: QCD with quarks and its pure-gauge, YM counterpart. The latter, however, has much higher characteristic temperatures than QCD with quarks: the "melting temperature" T YM where χ YM (T ) starts to decrease appreciably was found on lattice to be, for example, T YM ≈ 260 MeV [31,32] or even higher, T YM ≈ 300 MeV [33]. In contrast, the pseudocritical temperatures for the chiral and deconfinement transitions in the full QCD are lower than T YM by some 100 MeV or more (e.g., see Ref. [12]) due to the presence of the quark degrees of freedom.
This difference in characteristic temperatures, in conjunction with χ YM (T ) in WVRs (1) and (2) would imply that the (partial) restoration of the U A (1) symmetry (understood as the disappearance of the anomalous η 0 /η ′ mass) should happen well after the restoration of the chiral symmetry. But, this contradicts the RHIC experimental observations of the reduced η ′ mass [13,14] if WVRs (1), (2) hold unchanged also close to the QCD chiral restoration temperature T Ch , around which f π (T ) decreases still relatively steeply 4 [24] for realistic explicit ChSB, thus leading to the increase of 6χ YM (T )/f π (T ) 2 and consequently also of M η ′ .
There is still more to the relatively high resistance of χ YM (T ) to temperature: not only does it start falling at rather high T YM , but χ YM (T ) found on the lattice is falling with T relatively slowly. In some of the applications in the past (e.g., see Refs. [34,35]), it was customary to simply rescale a temperature characterizing the pure-gauge, YM sector to a value characterizing QCD with quarks. (For example, Refs. [34,35] rescaled T YM = 260 MeV found by Ref. [31] to 150 MeV). However, even if we rescale the critical temperature for melting of the topological susceptibility χ YM (T ) from T YM down to T Ch , the value of 6χ YM (T )/f π (T ) 2 still increases a lot [24] for the pertinent temperature interval starting 4 Relative to decay constants of mesons containing a strange quark; e.g., compare fss(T ) of the unphysical ss pseudoscalar with fπ(T ) in Fig. 1.
already below T Ch . This happens because χ YM (T ) falls with T more slowly than f π (T ) 2 . (It was found [24] that the rescaling of T YM would have to be totally unrealistic, to less than 70% of T Ch , in order to achieve sufficiently fast drop of the anomalous contribution that would allow the observed enhancement in the η ′ multiplicity.) These WVR-induced enhancements of the η ′ mass for T ∼ T Ch were first noticed in Ref. [24]. This reference used a concrete dynamical model (with an effective, rank-2 separable interaction, convenient for computations at T ≥ 0) [36] of low energy, nonperturbative QCD to obtain mesons as qq ′ bound states in Dyson-Schwinger (DS) approach [37][38][39], which is a bound-state approach with the correct chiral behavior (3) of QCD. Nevertheless, this concrete dynamical DS model was used in Ref. [24] to get concrete values for only the nonanomalous parts of the meson masses, but was essentially not used to get model predictions for the mass contributions from the gluon anomaly, in particular χ YM (T ). On the contrary, the anomalous mass contribution was included, in the spirit of 1/N c expansion, through WVR (1). Thus, the T -evolution of the η ′ -η complex in Ref. [24] was not dominated by dynamical model details, but by WVR, i.e., the ratio 6χ YM (T )/f π (T ) 2 . Admittedly, f π (T ) was also calculated within this model, causing some quantitative model dependence of the anomalous mass in WVR, but this cannot change the qualitative observations of Ref.
[24] on the η ′ mass enhancement. Namely, our model f π (T ), depicted as the dash-dotted curve in Fig. 1, obviously has the right crossover features [12]. It also agrees qualitatively with f π (T )'s calculated in other realistic dynamical models [10,38]. Various modifications were tried in Ref. [24] but could not reduce much the η ′ mass enhancement caused by this ratio, let alone bring about the significant η ′ mass reduction found in the RHIC experiments [13,14].
One must therefore conclude that either WVR breaks down as soon as T approaches T Ch , or that the Tdependence of its anomalous contribution is different from the pure-gauge χ YM (T ). We will show that the latter alternative is possible, since WVR can be reconciled with experiment thanks to the existence of another relation which, similarly to WVR, connects the YM theory with full QCD. Namely, using large-N c arguments, Leutwyler and Smilga derived [27], at T = 0, the relation (in our notation) between the YM topological susceptibility χ YM , and the full-QCD topological susceptibility χ, the chiral-limit quark condensate qq 0 , and m, the harmonic average of N f current quark masses m q . That is, m is N f times the reduced mass. In the present case of N f = 3, q = u, d, s, so that Equation (4) is a remarkable relation between the two pertinent theories. For example, in the limit of all very heavy quarks (m q → ∞, q = u, d, s), it correctly leads to the result that χ YM is equal to the value of the topological susceptibility in quenched QCD, χ YM = χ(m q = ∞). This holds because χ is by definition the vacuum expectation value of a gluonic operator, so that the absence of quark loops would leave only the pure-gauge, YM contribution. However, the Leutwyler-Smilga relation (4) also holds in the opposite (and presently pertinent) limit of light quarks. This limit still presents a problem for getting the full-QCD topological susceptibility χ on the lattice [40], but we can use the light-quark-sector result [27,41] where C m stands for corrections of higher orders in small m q , and thus of small magnitude. The leading term is positive (as qq 0 < 0), but C m is negative, since Eq. (4) shows that χ ≤ min(−m qq 0 /N f , χ YM ).
Although small, C m should not be neglected, since C m = 0 would imply, through Eq. (4), that χ YM = ∞. Instead, its value (at T = 0) is fixed by Eq. (4): All this starting from Eq. (4) has so far been at T = 0. If the left-and right-hand side of Eq. (4) are extended to T > 0, it is obvious that the equality cannot hold at arbitrary temperature T > 0. The relation (4) must break down somewhere close to the (pseudo)critical temperatures of full QCD (∼ T Ch ) since the pure-gauge quantity χ YM is much more temperature-resistant than the righthand side, abbreviated as χ. The quantity χ, which may be called the effective susceptibility, consists of the full-QCD quantities χ and qq 0 , the quantities of full QCD with quarks, characterized by T Ch , just as f π (T ). As T → T Ch , the chiral quark condensate qq 0 (T ) drops faster than the other DChSB parameter in the present problem, namely f π (T ) for realistically small explicit ChSB. (See Fig. 1 for the results of the dynamical model adopted here from Ref. [24], and, e.g., Refs. [10,38] for analogous results of different DS models). Thus, the troublesome mismatch in T -dependences of f π (T ) and the pure-gauge quantity χ YM (T ), which causes the conflict of the temperature-extended WVR with experiment at T > ∼ T Ch , is expected to disappear if χ YM (T ) is replaced by χ(T ), the temperature-extended effective susceptibility. The successful zero-temperature WVR (1) is, however, retained, since χ YM = χ at T = 0.
Extending Eq. (6) to T > 0 is something of a guesswork as there is no guidance from the lattice for χ(T ) [unlike χ YM (T )]. Admittedly, the leading term is straightforward as it is plausible that its T -dependence will simply be that of qq 0 (T ). Nevertheless, for the correction term C m such a plausible assumption about the form of T -dependence cannot be made and Eq. (7), which relates YM and QCD quantities, only gives its value at T = 0. We will therefore explore the T -dependence of the anomalous masses using the following Ansatz for the T ≥ 0 generalization of Eq. (6): where the correction-term T -dependence is parametrized through the power δ of the presently fastest-vanishing (as T → T Ch ) chiral order parameter qq 0 (T ).
The T ≥ 0 extension (8) of the light-quark χ, Eq. (6), leads to the T ≥ 0 extension of χ: (9) We now use χ(T ) in WVR instead of χ YM (T ) used by Ref. [24]. This gives us the temperature dependences of the masses in the η-η ′ complex, such as those in Fig. 2 illustrating the cases δ = 0 and δ = 1.
It is clear that χ(T ) (9) blows up as T → T Ch if the correction term there vanishes faster than qq 0 (T ) squared. Thus, varying δ between 0 and 2 covers the cases from the T -independent correction term, to (already experimentally excluded) enhanced anomalous masses for δ noticeably above 1, to even sharper mass blow-ups for δ → 2 when T → T Ch . On the other hand, it does not seem natural that the correction term vanishes faster than the fastest-vanishing order parameter qq 0 (T ). Indeed, already for the same rate of vanishing of the both terms (δ = 1), one can notice in Fig. 2 the start of the precursors of the blow-up of various masses in the η ′ -η complex as T → T Ch although these small mass bumps are still experimentally acceptable. Thus, in Fig. 2 we depict the δ = 1 case, and δ = 0 (T -independent correction term) as the other acceptable extreme. Since they turn out to be not only qualitatively, but also quantitatively so similar that the present era experiments cannot discriminate between them, there is no need to present any 'in-between results', for 0 < δ < 1.
Next we turn to completing the explanation how the above-mentioned results in Fig. 2 were obtained. Using χ(T ) in WVR instead of χ YM (T ) used by Ref. [24], does not change anything at T = 0, where χ(T ) = χ YM (0), which remains an excellent approximation even well beyond T = 0. Nevertheless, this changes drastically as T approaches T Ch . For T ∼ T Ch , the behavior of χ(T ) is dominated by the T -dependence of the chiral condensate, tying the restoration of the U A (1) symmetry to the chiral symmetry restoration.
As for the nonanomalous contributions to the meson masses, we use the same DS model (and parameter values) as in Ref. [24], since it includes both DChSB and correct QCD chiral behavior as well as realistic explicit ChSB. That is, all nonanomalous results (M π , f π , M K , f K , as well as M ss and f ss , the mass and decay constant of the unphysical ss pseudoscalar meson) in the present paper are, for all T , taken over from Ref. [24]. We used this same model also for computing the chiral quark condensate qq 0 , including its T -dependence displayed in Fig. 1.
This defines completely how the results displayed in Fig. 2 were generated. For details, see Ref. [24] (and Ref. [30] for M η0 and M η8 ). Here we list only the formulas which, in conjunction with Fig. 1, enable the reader to understand easily the T -dependences of the masses in Fig. 2: The theoretical η ′ and η mass eigenvalues are In all expressions after Eq. (11), the T -dependence is understood.
In both cases considered for the topological susceptibility (8) [δ = 0, i.e., the constant correction term, and δ = 1, i.e., the strong T -dependence ∝ qq 0 (T ) of both the leading and correction terms in χ(T )], the results are consistent with the experimental findings on the decrease of the η ′ mass of Csörgő et al. [13,14].

III. SUMMARY, DISCUSSION AND CONCLUSIONS
In the light of the recent experimental results on the η ′ multiplicity in heavy-ion collisions [13,14], we revis-ited the earlier theoretical work [24] concerning the thermal behavior of the η ′ -η complex following from WVR straightforwardly extended to T > 0. We have confirmed the results of Ref. [24] on WVR where the ratio χ YM (T )/f π (T ) 2 dominates the T -dependence, and clarified that these results are practically model-independent. It is important to note the difference between our approach and those that attempt to give model predictions for topological susceptibility, such as Refs. [42,43]. By contrast, in Refs. [24] and here, as well as earlier works [26,30,44,45] at T = 0, a DS dynamical model is used (as far as masses are concerned) to obtain only the nonanomalous part of the light pseudoscalar meson masses (where the model dependence is however dominated by their almost-Goldstone character), while the anomalous part of the masses in the η ′ -η complex is, through WVR, dictated by 6χ YM /f 2 π . In this ratio, f π (T ) is admittedly model-dependent in quantitative sense, but other realistic models yield qualitatively similar crossover behaviors [46] of f π (T ) for m q = 0, as exemplified by our Fig. 1, and Fig. 2 in Ref. [24], and by Fig. 6 in Ref. [10]. Such f π (T ) behaviors are also in agreement with the T -dependence expected of the DChSB order parameter on general grounds: a pronounced fall-off around T Ch -but exhibiting, in agreement with lattice [12], a smooth crossover pattern for nonvanishing explicit ChSB, a crossover which gets slower with growing m q [e.g., compare f π (T ) with f ss (T ) in Fig. 1]. In contrast to the QCD topological susceptibility χ, the YM topological susceptibility and even its T -dependence χ YM (T ), including its "melting" temperature T YM , can be extracted [24] reasonably reliably from the lattice [22,31]. Thus, it was not modeled in Ref. [24]. Hence our assertion that the results of Ref. [24] unavoidably imply that the straightforward extension of WVR to T > 0 is falsified by experiment [13,14], especially if one recalls that even the sizeable T -rescaling [34,35] T YM → T Ch was among the attempts to control the η ′ mass enhancement [24]. Nevertheless, we have also shown that there is a plausible way to avoid these problems of the straightforward, naive extension of WVR to T > 0, and this is the main result of the present paper. Thanks to the existence of another relation, Eq. (4), connecting the YM quantity χ YM with QCD quantities χ and qq 0 , it is possible to define a quantity, χ, which can meaningfully replace χ YM (T ) in finite-T WVR. It remains practically equal to χ YM up to some 70% of T Ch , but beyond this, it changes following the T -dependence of qq 0 (T ). In this way, the successful zero-temperature WVR is retained, but the (partial) restoration of U A (1) symmetry [understood as the disappearing contribution of the gluon anomaly to the η ′ (η 0 ) mass] is naturally tied to the restoration of the SU A (3) flavor chiral symmetry and to its characteristic temperature T Ch , instead of T YM .
It is very pleasing that this fits in nicely with the recent ab initio theoretical analysis using functional methods [47], which finds that the anomalous breaking of U A (1) symmetry is related to DChSB (and confinement) in a self-consistent manner, so that one cannot have one of these phenomena without the other.
Of course, the most important thing is that this version of the finite-T WVR, obtained by χ YM (T ) → χ(T ), is consistent with experiment [13,14] for all reasonable strengths of T -dependence [0 ≤ δ < ∼ 1 in Eq. (8)]. Namely, the both tablets in Fig. 2 show, first, that η ′ mass close to T Ch suffers the drop of more than 200 MeV with respect to its vacuum value. This satisfies the minimal experimental requirement abundantly. Second, Fig.  2 shows an even larger drop of the η 0 mass, to some 400 MeV, close to the "best" value of the in-medium η ′ mass (340 MeV, albeit with large errors) obtained by Csörgő et al. [13,14]. This should be noted because the η 0 mass inside the fireball is possibly even more relevant. Namely, although it is, strictly speaking, not a physical meson, η 0 is the state with the qq content closest to the qq content of the physical η ′ in the vacuum. Thus, among the isoscalar qq states inside the fireball, η 0 has the largest projection on, and thus the largest amplitude to evolve, by fireball dissipation, into an η ′ in the vacuum.