Critical exponents and transport properties near the QCD critical endpoint from the statistical bootstrap model

We present an estimate of the behaviour of the shear and bulk viscosity coefficients when the QCD critical point is approached from the hadronic side, describing hadronic matter within the statistical bootstrap model (SBM) of strong interactions. The boostrap model shows critical behavior near the quark-hadron transition temperature if the parameter characterizing the degeneracy of Hagedorn states is properly chosen. We calculate the critical exponents and amplitudes of relevant thermodynamic quantities near the QCD critical point and combine them with an Ansatz for the shear and bulk viscosity coefficients to derive the behavior of these coefficients near the critical point. The shear viscosity to entropy density ratio is found to decrease when the temperature is increased, and to approach the Kovtun-Son-Starinets bound 1/(4{\pi}) faster near the critical point, while the bulk viscosity coefficient is found to rise very rapidly.


I. INTRODUCTION
The phase diagram of strongly interacting matter has attracted a lot of attention in the past couple of decades [1][2][3]. At low temperature T and low baryon chemical potential µ B , quantum chromodynamics (QCD) matter consists of colorless hadrons, while at high temperature and baryon density the degrees of freedom are colored quarks and gluons. Lattice QCD simulations indicate that, at µ B ∼ 0, the transition from the hadronic phase to the quark-gluon plasma phase (QGP) is actually an analytic cross over [4][5][6][7][8][9]. At finite baryon densities µ B /T 1, however, lattice QCD calculations are affected by the sign problem [10][11][12][13], and, even though very recent works, based on the analytical continuation from imaginary to real µ B , are now probing values of the chemical potential up to real µ B ∼ 300 MeV [14], results at large values of µ B are still scarce. Hence, in this regime of the phase diagram one has to resort to effective models of QCD, such as the Nambu-Jona-Lasinio model [15,16], the quark-meson-coupling model [17] etc. Finally, at low temperature and sufficiently large baryon density many phenomenological models predict the transition between the hadronic phase and the deconfined phase to be of first order [18][19][20][21][22][23][24][25][26]: see also the discussion in ref. [27].
Apart from the effective QCD models mentioned above, a simple model to describe the hadronic phase of QCD is the hadron resonance gas model (HRG). This model is based on the S-matrix formulation of statistical mechanics [28]. At low density, as it turns out, the thermodynamics can be approximately modeled in terms of a non-interacting gas of hadrons and resonances [29][30][31]. The predictions of this model have been compared with lattice QCD simulations, finding good agreement for temperatures up to T ∼ 150 MeV except for some discrepancies in the trace anomaly [32][33][34]. Later studies found that the agreement can be improved, if the contribution of a continuous density of states is included in the mass spectrum of the HRG [35][36][37][38]. Remarkably, analogous results have been obtained also in lattice simulations of SU(N ) Yang-Mills theories without dynamical quarks [39,40], and even in three spacetime dimensions [41].
A description of the density of hadron states in terms of a continuous distribution is the basis of the statistical bootstrap model (SBM) [42,43], which attracted a lot of attention in the particle physics community in the pre-QCD era. The mass spectrum of abundant formation of heavy resonances and higher angular momentum states can be consistently described by a self-similar structure of hadrons through the bootstrap condition. These high-mass resonances have an interesting effect on the strong interaction thermodynamics: in the thermodynamic system dominated by exponentially rising resonance states there is a finite limiting temperature T H , called Hagedorn temperature. The existence of this limiting temperature indicates that the hadron resonance gas cannot exist at physical temperatures T > T H , and suggests that strongly interacting matter should then enter a different phase. The bootstrap condition of the SBM requires the density of states to be of the form ρ(m) ∼ m a exp(m/T H ) [44][45][46], where a is a constant. Interestingly, the string model (or dual resonance model) of strong interactions [47] also predicts this type of density of states. The a constant plays an important role in determining the thermodynamics of the SBM near the Hagedorn temperature. In fact, for the choice a = −4 both the energy density and the entropy density remain finite near T H and one expects a phase transition to take place [24,44,48], so that T H can be interpreted as a critical temperature, T c .
A particularly interesting point in the QCD phase diagram is the conjectured critical end point (CEP). It should be remarked that, so far, the existence of the CEP has neither been proven theoretically, nor has it been observed experimentally. However, its existence is strongly suggested by the aforementioned model calculations investigating the phase diagram region at low temperatures and baryon densities larger than that of nuclear matter, which predict a first-order transition line separating the hadronic phase from a deconfined phase: since that line is known not to extend all the way to the µ B = 0 axis (where the transition is actually a crossover), it should end at a CEP, where the transition should be a continuous one [49]. A lot of theoretical investigation has been carried out, and is still going on, to locate the CEP and predict possible experimental signatures, see refs. [50][51][52][53] for reviews. On the experimental side, an entire experimental program, namely the Beam Energy Scan (BES) program, has been devoted at the Super Proton Synchrotron (SPS) and at the Relativistic Heavy Ion Collider (RHIC) to search for the CEP [54,55]. In particular, as suggested in ref. [49], the existence and the location of the CEP could be revealed by the observation of a suppression in temperature and chemical potential fluctuations on an event-by-event basis, and by large fluctuations in the multiplicity of low-energy pions.
A very important feature of the critical point is the emergence of a universal critical mode. As the system approaches the critical point, this mode rises very rapidly with some power of the correlation length ξ, which eventually diverges at the critical point. For instance, the variance, skewness and kurtosis of the non-Gaussian fluctuation of the critical mode grow as ξ 2 , ξ 9/2 and ξ 7 , respectively [56][57][58]. In the experimental search for the critical endpoint, these critical fluctuations can be accessed by measuring event-by-event fluctuations of particle multiplicities [55].
While "static" critical phenomena have been extensively studied theoretically, an avenue that has been explored less is the one of "dynamical" critical phenomena. As it turns out, critical singularities can also occur in quantities encoding the dynamical properties of the medium, like transport coefficients. Away from the critical point, the dynamic properties of a system can be characterized by hydrodynamics, which provides an effective description of the fluid in the lowfrequency, long-wavelength limit. Hydrodynamics describes fluctuations of conserved quantities at equilibrium and any additional slow variable that occurs due to the existence of a spontaneously broken symmetry. In the hydrodynamic effective theory the dynamical critical fluctuations are described by coupling the order-parameter field with the conserved momentum density. In this model, which is called the H model in the classification of dynamical critical phenomena [59] by Hohenberg and Halperin, the transport coefficients depend on the correlation length as This behavior suggests that the transport coefficients would affect the bulk hydrodynamic evolution of the matter created in heavy-ion collisions near the QCD critical point [60][61][62][63][64][65][66][67].
It is worth emphasizing that, while lattice calculations remain the tool of choice for theoretical first-principle studies of various quantities relevant for strong-interaction matter, their applicability in studies of transport coefficients in the proximity of the QCD critical endpoint is severely limited, for a two-fold reason. On one hand, as we remarked above, the existence of the sign problem poses a formidable barrier to lattice simulations at finite baryon-number density: a barrier that might even be impossible to overcome with classical computers, if it is related to fundamental computational-complexity issues [68]. On the other hand, even at zero baryon-number density, the lattice determination of transport coefficients of QCD matter involves its own difficulties, due to the fact that lattice QCD calculations are done in a Euclidean spacetime, and typically the extraction of quantities involved in real-time dynamics requires a Wick rotation back to Minkowski signature, with the reconstruction of a continuous spectral function from a finite set of Euclidean data, which is an ill-posed numerical problem [69]. Despite some recent progress (see, e.g., refs. [70][71][72][73][74][75]), a general solution to this type of problems is still unknown.
For these reasons, phenomenological models remain a useful theoretical tool to get some insight into the physics near the QCD critical endpoint. In this work, we extract critical exponents [44] and amplitudes of thermodynamic quantities relevant near the critical point within the statistical bootstrap model. We then derive the singularities characterizing shear and bulk viscosity coefficients, starting from an Ansatz for viscosity coefficients [76] that is suitable for a hydrodynamic system with conserved baryon charge. We then estimate viscosity coefficients near the critical point from the hadronic side using the critical exponents of this model.
We organize the paper as follows. In section II we review the derivation of the critical exponents (and amplitudes) close to the critical point in the critical bootstrap model, that was first worked out in ref. [44], with a few additional remarks and comments. In section III we derive the singularities of shear and bulk viscosity near the critical point. In section IV we present our results. Finally, in section V we summarize our findings and conclude. Throughout the paper, we work in natural units ( = c = k B = 1).

A. Critical exponents
The analysis of critical phenomena is based on the assumption that, in the T → T c limit, any relevant thermodynamic quantity can be separated into a regular part and a singular part. The singular part may be divergent or it can have a divergent derivative. It is further assumed that the singular part of all the relevant thermodynamic quantities is proportional to some power of t, where t = (T − T c )/T c . These powers, called critical exponents, characterize the nature of singularity at the critical point. The critical exponents,α,β,γ andν are defined through the following power laws [77] (in the limit t → 0 − ): where C V , n B,c , k T and ξ respectively denote the specific heat, the critical baryon density, the isothermal compressibility and the correlation length, while C − , N − , K − and Ξ − are the corresponding amplitudes from the hadronic side (T < T c ). Note that eq. (3) is an equation of state, relating baryon density n B and pressure p near the critical point.

B. Formulation of the model
We follow ref. [44] to extract the amplitudes and critical exponents within the statistical bootstrap model. We first discuss the case of vanishing baryonic chemical potential, µ B = 0. Consider an ideal gas of hadrons and all possible resonance states as non-interacting constituents: the partition function of this system can be written as [78] where ρ(m) is the hadron spectrum included in the HRG model. In the simplest formulation of the model, that was discussed in ref. [44], only pions were considered as the "basic" hadrons. More recently, however, it has become customary to include all the hadrons and resonances that have been detected experimentally up to some energy scale M and take ρ(m) = j δ(m − m j ). Such discrete mass spectrum leads to the physical hadron resonance gas model. In the physical HRG, if g i is the degeneracy of the i-th hadronic species, then for spin degrees of freedom the degeneracy factors turn out to be g i ∼ m 2 i [48]. Thus, one sees that the spin multiplicity already can result in an unbounded increase in resonances. The upshot of the m 2 dependence of resonance degeneracy is that the partition function of the physical resonance gas and all of its higher-order derivatives remain finite at T c . Thus, the required degeneracy structure is absent in the physical resonance gas and hence it does not show critical behavior.
It turns out that the degeneracy structure required to show critical behavior is present in the Hagedorn density of states, which can be used to model the spectral density above M in terms of a continuous distribution. Consider a density of states of the form where the sum ranges over all hadrons species with mass m i ≤ M , g i denotes the degeneracy of each species, while ρ H is the continuous contribution to the density of states. For our analysis, we included all hadronic states reported in ref. [79] with masses not larger than M = 2.25 GeV.
It should be noted that choosing a different M value in the same ballpark would not lead to significant differences e.g. in the equilibrium thermodynamic quantities in the low-temperature phase. The reason for this is that at low temperatures the thermodynamics is dominated by the lightest hadrons, and including or not including the contribution from some discrete heavy states does not have significant impact on the equation of state at low T . Note that, in the simplest possible formulation of the model, the discrete part of the spectrum could include only pions, and one could model all the states of the spectrum above the two-pion threshold (setting M = 2m π ) in terms of a Hagedorn density of states: where g π denotes the pion degeneracy, which is equal to 3. While such picture is clearly a very crude model of the hadron spectrum, it still captures some interesting finite-temperature features, at least at a qualitative level, and is useful to highlight some general consequences for the thermodynamics and transport properties. The logarithm of the partition function (6) with the density of states in eq. (7) is written as in which the first summand on the right-hand side, which does not depend on the Ansatz for the continuous part of the density of states, encodes the contribution from a gas of non-interacting hadrons in the discrete part of the spectrum (i.e. hadrons whose masses are not larger than M ). In particular, the contribution to ln Z discrete due to pions can be written in the form where K n (z) is a modified Bessel function of the second kind of order n. For large real values of its argument, one has The contributions to ln Z discrete from the other hadron species in the discrete part of the spectrum can be derived in a similar way, and one obtains where the sum over i ranges over all hadrons with mass m i ≤ M , as in eq. (7), and η i = −1 for bosons, while η i = 1 for fermions. The second summand appearing on the right-hand side of eq. (9) represents the contribution due to the continuous part of the spectrum: where, as above, M is the threshold separating the discrete (for m ≤ M ) and the continuous (for m > M ) parts of the spectrum. Performing the momentum integration, one obtains Using eq. (11), for m/T 1 one gets All the thermodynamic functions can be readily obtained from the partition function in eq. (13) once the continuous part of the mass spectrum ρ(m) is specified.
In the statistical bootstrap model (see refs. [80,81] for reviews), hadronic matter at high temperature is dominated by formation of resonances whose number grows exponentially. The bootstrap condition leads to a solution for the mass spectrum of the form [42,43] ρ H (m) = A m a e bm (17) where A, a, and b are constant parameters. In particular, the parameter A provides the normalization of the resonance contributions relative to that of the pions. The parameter a specifies the nature of the degeneracy of high-mass resonances, and also determines the critical behavior of hadronic matter. One possible solution of the bootstrap condition was derived in ref. [82], yielding a 3. Finally, the parameter b turns out to be the inverse of the Hagedorn temperature at which thermodynamic functions show singular behavior.
Restoring the dependence on the fugacity z B = exp(µ B /T ), the contribution to the partition function associated with the continuous spectrum (17) can be written as At this point, we should stress an important observation: in order to obtain eq. (18), in which z B is factorized on the right-hand side, it has been implicitly assumed that b is independent from µ B , i.e. that the critical temperature T c does not depend on the fugacity z B . Strictly speaking, however, this is not fully justified: as has been discussed in detail in the literature [83][84][85], in the presence of arbitrary fugacity z B , the bootstrap equation takes the form where φ is an input function, receiving contributions from the physical hadrons of the theory, while G, which encodes their interactions in terms of the bootstrap picture (whereby strongly interacting systems of particles form larger clusters of particles, which in turn form larger clusters, etc.) is the Laplace transform of the mass spectrum. Eq. (19) has a square-root branch point singularity for φ = 2 ln 2 − 1 (or, equivalently, for G = ln 2), which defines the boundary of the hadronic phase in this model through a non-trivial relation between T and z B . In other words, strictly speaking, the critical temperature T c is a non-trivial function of z B . In eq. (18) and in the rest of this work, however, we assume that the dependence of T c on the fugacity is mild, i.e. we work in the approximation in which b = 1/T c is constant. While this simplification may appear to be crude, it is worth noting that during the past few years lattice QCD calculations have conclusively proven that the change of state between the hadronic, broken-chiral-symmetry phase and the deconfined, chirally symmetric phase at zero chemical potential is a crossover [86][87][88], and that at small but finite values of µ B the curvature of the line describing the crossover in the QCD phase diagram is very small [89][90][91][92][93]. As a consequence, it is not unreasonable to expect that, even within the approximation of a critical temperature independent from the chemical potential, the statistical bootstrap model may still capture the physics close to a possible critical endpoint of QCD at finite chemical potential. Assuming T c to be approximately independent from µ B simplifies the expression of the partition function, and allows one to get more analytical insight into the physical quantities of interest. In a nutshell, the fact that z B factors out in the expression of the logarithm of Z H implies that the dependence on the chemical potential in this model is somewhat "trivial". While the validity of this approximation at large values of µ B is not obvious, lattice results lead us to think that its use at least for small and intermediate values of µ B should be a reasonable approximation. With these caveats in mind, in the next section we shall calculate the critical exponents by taking appropriate derivatives of the partition function (18) and then taking the T → 1/b limit.

C. Critical exponents in the statistical bootstrap model
With the change of variable having expressed the integral in terms of the upper incomplete Γ function. The energy density can then be written as and for T → 1/b one finds that Hence for a < −7/2 the energy density remains finite (and approaches some critical value ε c ) as T → T H , implying that the system cannot exist in this state for ε > ε c and suggesting that a phase transition must take place.
The specific heat at constant volume can then be written as and for T → 1/b one gets Comparing eq. (24) with eq. (2) we deduce the amplitude C − as while the critical exponentα readŝ The baryon number density n B can be evaluated as hence for T close to 1/b we get the critical density as The inverse of the isothermal compressibility is defined as and for a non-interacting resonance gas it takes the following, very simple form: For temperatures close to 1/b, one obtains, from which it is straightforward to deduce the amplitude K − for a < −5/2 (32) and the critical exponentγ asγ We note that a continuous density of states with a < −7/2 makes the energy and entropy densities finite, while all higher-order derivatives diverge near T H . In our analysis of transport coefficients we shall consider the a = −4 case [24,48,94] which leads to normal behavior of the hadronic system near the boundaries of the quark-hadron phase transition line, since it does not allow the energy density to become infinite even for pointlike particles.
At this point, an important observation is in order. Hadrons are not elementary, pointlike particles: rather, they arise as color-singlet bound states of the strong interaction, and, for this reason, they can be associated with a characteristic finite size, of the order of the fm. As a consequence of the very nature of hadrons as complex bound states of relativistic, strongly interacting constituents (which defies a description in terms of sufficiently simple phenomenological models), the measurement and even the definition of hadron sizes are, in general, non-trivial (see, for example, ref. [95] for an experimental determination of the radius of a well-known hadron: the proton). It is worth noting that, if corrections related to the finiteness of the particles' physical size are taken into account in our model, the restriction on the admissible values of a become milder, in the sense that finite-particle-size corrections make some of the divergent quantities obtained in the pointlike approximation finite. The fact that finite-particle-size effects can have even a qualitative impact on the details of the description of the thermodynamics of the confining phase of QCD is hardly surprising, as it is well known that they have a significant role in fits of particle multiplicities produced in heavy-ion collisions [96][97][98], and even in the interpretation of non-perturbative theoretical predictions from lattice simulations [99]. For this reason, in a more complete discussion, a priori one should not discard the a values that lead to unphysical infinities for a system of pointlike particles. However, a fully systematic discussion of finite-particle-size effects would involve a non-trivial amount of additional technicalities (and a certain degree of arbitrariness in the way to define these effects), and lies beyond the scope of our present work. For this reason, in the following we restrict our attention to the simpler, idealized case of pointlike particles, which is nevertheless expected to provide a reasonable approximation of the physics that is studied in currents experiments, especially in view of the fact that the typical sizes of the systems produced in nuclear collisions are significantly larger than hadron sizes [100], and which does not introduce additional parameters in the description.

III. TRANSPORT COEFFICIENTS NEAR THE CRITICAL POINT
Approaching the critical point, the thermodynamic quantities relevant for the computation of transport coefficients are: energy density (ε), baryon number density (n B ), specific heats (C V and C p ), isothermal compressibility (k T ), speed of sound (C s ) and correlation length (ξ). A set of Ansätze for the transport coefficients near the critical point can be written in terms of thermodynamic quantities as [76] Near the critical point, the correlation length ξ is the only relevant length scale. Further, longitudinal perturbations can be assumed to be those of the non-equilibrium modes near T c . A particularly simple form of the functions F η,ζ , namely Here, h c and s c respectively denote the enthalpy and entropy densities at T c , both of which are finite when one sets a = −4 in the Hagedorn density of states, while λ c = (∂p/∂T ) V at T = T c . The amplitudes f η and f ζ are free parameters, which can be fixed by imposing some constraint on the viscosity coefficients near T c . For instance, as we already mentioned, the gauge-string duality [101][102][103] suggests a universal lower bound 1/(4π) for the η/s ratio [104]. Similar constraints can be imposed on the ζ/s ratio, too.  [36,48]. According to the discussion in section II, the parameter b is set to the inverse of T c , whose value is T c = 0.160 GeV.

density of states a b [GeV
Note that the value of A 1 chosen for the ρ 1 model corresponds to A 1 = 15T 3 c , as discussed in the text.
Before discussing the behavior of viscosity coefficients near T c , it is instructive to point out a few remarks about the thermodynamics of the model. In table I, we report the parameters of the continuous part of the density of states, taken from refs. [36,48]. Note that, as discussed in The quark chemical potential is assumed to be µ B = 0 for the blue curves, while the red curves are obtained at µ B = 220 MeV. The parameter b is set to the inverse of the critical Hagedorn temperature, as discussed in the paragraph after eq. (17). In addition, we also plot the curves representing the contribution due to an ideal pion gas (dotted green curves), i.e. to the lightest states in the discrete part of the spectrum in eq. (7), which does not depend on the functional form that is assumed to model the continuous part of the spectrum. those references, the continuous part of the density of states is assumed to start at mass values corresponding to the pion-pair threshold, and that, in addition to the continuous part, the density of states also includes a δ-like contribution at the pion mass. From the density of states constructed using the parameters in table I, one obtains the equilibrium thermodynamic quantities shown in fig. 1, namely the pressure (p), the energy density (ε) and the entropy density (s), in units of T 4 (for p and ε) and T 3 (for s). We have also plotted the trace of the energy-momentum tensor ∆ in units of the fourth power of the temperature, ∆/T 4 = (ε − 3p)/T 4 . The solid blue curves correspond to a continuous density of states of the form ρ 1 (m) = A 1 m −4 exp(bm) at vanishing chemical potential, whereas the solid red ones are obtained at µ B = 220 MeV. To give an idea of the dependence of equilibrium thermodynamic quantities on a, we also show the results that one would obtain for a different value of a, i.e. for a spectral density with continuous part ρ 1 (m) = A 1 m −17/4 exp(bm), which are shown, with the same color code, by the dashed curves. One immediately realizes that, as compared with the solid curves, the dashed ones exhibit only small quantitative differences. The reason for the choice a = −17/4 stems from the fact that, as was discussed in subsection II C, the specific heat exhibits power-law behavior only if a is larger than −9/2. On the other hand, we also remarked that a is constrained to be less than −7/2, because in this range the energy density remains finite when T tends to the critical temperature. This leaves us with (−9/2, −7/2) as the most interesting interval of values for a. Thus, a = −17/4 is a value which is exactly equidistant from our choice a = −4 and the lower end of the interval of interesting values, and as such is expected to reveal some information on the dependence of our results on the choice of a. As the plots in fig. 1 clearly show, this dependence is very mild, indicating that our predictions for these quantities are robust (at least within the interval of a values, i.e. −9/2 ≤ a ≤ −7/2).
Finally, the dotted green curves show the contributions from the ideal pion gas, i.e. the lightest hadrons included in the discrete and model-independent part of the density of states in eq. (7), which can be directly derived from eq. (10): for example, the pion-gas contribution to the pressure (that one can denote as p π ) can be written as It is known from comparison with lattice QCD results (as reviewed, for instance, in the recent ref. [105]) that the hadron resonance gas model provides a very accurate description of the equation of state for all temperatures below T c . The contribution to thermodynamics from the part of the hadronic spectrum that is modelled in terms of a continuous density of states becomes significant when the temperature is sufficiently large. Nevertheless, in the case of ρ 1 with a = −4 both the energy and the entropy densities remain finite for T → T − c . This reflects the fact that, for a = −4, the second derivative of the partition function is divergent, but the first is not. In fact, setting A 1 = 15T 3 c = 0.06144 GeV 3 corresponds to ε c /T 4 c 4 [48]. It is worth noting that the bootstrap model predicts the existence of a phase transition at the finite critical temperature T c . This can be interpreted by saying that this phenomenological model, which provides a description for the thermodynamics of hadronic matter in rather simple terms (e.g. neglecting hadron-hadron interactions) and without reference to the microscopic QCD Lagrangian, is able to capture the existence of a finite temperature, above which hadrons cannot exist anymore. To draw an analogy with the description of physics at the electro-weak scale within and beyond the Standard Model, the statistical bootstrap model can be interpreted as an "effective field theory" describing the thermal properties of nuclear matter in terms of its "low-energy degrees of freedom" (i.e. those that manifest themselves at energy scales below the characteristic hadronic scale, O(10 2 ) MeV), and its breakdown at a finite temperature T c hints at the existence of "new physics" above that scale. In this case, the "new physics" above that temperature is the quark-gluon plasma, whose existence could be argued (and reconciled with the bootstrap model [42,43]) after the introduction of QCD [106]. In this analogy, QCD plays the role of the "more fundamental theory", which holds up to higher energies (being, in fact, a renormalizable, asymptotically free and ultraviolet-complete theory) and at the same time reduces to the "effective model" at low energies, by predicting the existence of massive hadrons through the mechanisms of color confinement and dynamical chiral symmetry breaking [107]. One should remark that, despite the remarkable qualitative prediction of a finite maximal temperature at which hadrons exist, the bootstrap model does not capture all quantitative details of the change of phase between hadronic matter and the quark-gluon plasma: in particular, non-perturbative lattice calculations based on the QCD Lagrangian show that, for zero or nearly zero values of the baryonic chemical potential, this change of phase is actually an analytical crossover, rather than an actual phase transition (see refs. [8,108] and references therein). As a consequence, the statistical bootstrap model prediction (for a = −4 and at zero net baryon density) of a phase transition with critical exponentsα = 1/2,β = 1, andγ = 0 is disproven by lattice QCD. Still, the statistical bootstrap model remains a useful phenomenological model, in particular when studying regions of the phase diagram at large baryonic densities, where a critical endpoint might exist, and in which, as we already pointed out in section I, lattice QCD calculations are hampered by particularly severe computational challenges. to entropy density ratio with those from other models [76,104,109,110]. Our predictions for a = −4, denoted by the solid blue line, are also compared with those for a = −17/4 (shown by the dotted blue curve, which is nearly indistinguishable from the former), and the difference δ η/s is plotted in the inset figure.
The right-hand-side panel shows the prediction of the statistical bootstrap model for the bulk viscosity to entropy density ratio (for a = −4 and for a = −17/4, and the difference between the two, denoted by δ ζ/s and displayed in the inset figure) and its comparison with other works [60,61,76]. Our results correspond to f η = 0.5 and f ζ = 0.85, which are fixed by requiring consistency with other models in the vicinity of the critical point, as discussed in the text. Fig. 2 shows the predictions for the shear and bulk viscosities near the QCD critical point based on eqs. (36) and (37). We take the correlation length amplitude to be Ξ − = 1 fm and the estimate for the critical point location to be (T c , µ B,c ) = (160 MeV, 220 MeV) [111]. For a = −4 one can easily derive the critical exponents and amplitudes needed for the estimate of the viscosity coefficients near T c . The critical exponents are not independent but are constrained by scaling laws. In particular, the exponentsα andν are related by the Josephson scaling lawνd = 2 −α, where d is the number of space dimensions [77].
In the left-hand-side panel of fig. 2, the solid blue curve shows the shear viscosity to entropy density ratio within the statistical bootstrap model, with density of states specified by ρ 1 , and with the f η parameter appearing on the right-hand side of eq. (36) set to 0.5. Note that the choice of this amplitude value, which we have done with the procedure discussed below, introduces some systematic uncertainties. On the other hand, to give an idea of the dependence of this prediction on the parameter a, we also present the prediction that one would obtain for a = −17/4 (again with f η = 0.5), which is displayed by the dotted blue curve, and which is nearly indistinguishable from the latter. Hence, in the inset plot we show the quantity δ η/s , defined as the prediction of the bootstrap model for η/s for a = −17/4, minus the one for a = −4: the relative difference between the predictions corresponding to the two a values is at the per mille level. We conclude that the dependence of our prediction on a (within the range of values of a of our interest) has a negligible impact on the uncertainties affecting the prediction for the η/s ratio. In the larger plot in the figure, we also compare the critical solution for η/s with those obtained from various other models: in particular, the dashed black curve corresponds to the conjectured universal lower bound 1/(4π) for this ratio, that was derived in ref. [104], while the brown curve describes the result that one would obtain for a pion gas [109], and the magenta curve shows the result that can be derived assuming the medium to be described in terms of a hadronic mixture [110] at low temperature and density. Finally, the red curve corresponds to the same solution for the viscosity coefficients as in this work, but with the critical exponents of the three-dimensional Ising model and the amplitudes constrained by universality arguments [76]. In our case, we chose to fix the f η amplitude to optimize the consistency with the other predictions shown in the figure at temperatures −0.2 t −0.1: in particular, in that temperature interval our choice yields an almost perfect consistency with the curve predicted in ref. [110], which is the one that is intermediate among those predicted in those works. We should remark, however, that in general the choice of the f η amplitude remains a source of systematics that are difficult to quantify (and, hence, the value that we quote should be taken cum grano salis). Note, however, that, as shown by the two curves derived in ref. [76] with two different choices for f η , i.e. the solid and dash-dot-dotted red lines, the choice of the numerical value of the amplitude has a strong impact at temperatures far from the critical point, but this discrepancy is already reduced to small values for reduced temperatures between approximately −0.2 and −0.1. We note that the critical behavior of the statistical bootstrap model leads to a linear decrease in η/s as a function of the temperature, and that at low temperatures the estimated magnitude of η/s is in agreement with that of a pion gas, or of the hadron gas mixture. Near T c there is a mild violation of the bound conjectured in ref. [104] (which could make it problematic to fix f η through some constraint in a region of temperatures very close to T c ). Such violation has also been noted for one of the solutions discussed in ref. [76], shown by the dashed red curve in fig. 2.
The right-hand-side panel of fig. 2 shows the bulk viscosity to entropy density ratio, in which one notes that the statistical bootstrap model predicts a rapid increase in the bulk viscosity as a function of the temperature. Also in this case, we present our results both for a = −4 (solid blue line) and for a = −17/4 (dotted blue curve), and the difference between the latter and the former, which is denoted by δ ζ/s and shown by the dashed blue line in the inset plot. In this case, the relative difference between the predictions corresponding to the two a values is below 10 −2 , meaning that also for ζ/s the dependence on a induces a very mild systematic uncertainty. Near T c , our results, with the amplitude coefficient appearing in eq. (37) fixed to f ζ = 0.85 by requiring an approximate match with those of refs. [60,61] at t −0.1, are in remarkable agreement with those from that work (shown by the dashed green curve), where the bulk viscosity has been estimated under the assumption of a QCD critical point belonging to the dynamical universality class of the so-called H model [59]. Remarkably, this agreement between the two curves is observed for essentially all negative values t −0.1, which is non-trivial, as that is the region in which the ζ/s ratio grows rapidly to very large values. Nevertheless, also in this case the readers should be warned that there is no obvious method to fix the value of f ζ in a unique, completely rigorous way from first principles, and the systematic uncertainties associated with any choice remain difficult to assess. For comparison, in the plot we also show the prediction for ζ/s from ref. [76]. Coming to the interpretation of fig. 2, we note that a large bulk viscosity should manifest itself in heavyion collisions through the decrease of the average transverse momentum of final-state hadrons. Moreover, due to the increase in entropy associated with the dissipation through large bulk viscosity, this effect should be accompanied by an increase in total multiplicity for final-state hadrons. The large bulk viscosity near the critical point would play a particularly important role in the elliptic flow measurement of the matter produced in the BES program.
Note that the features of the transport coefficients predicted by our model are only expected to hold close to T c , and there is no reason to expect the curves plotted in fig. 2 to be quantitatively accurate predictions at temperatures much smaller than the critical one. The reason for this was already discussed in ref. [76], in which it was remarked that the extrapolation of power-law behavior beyond the critical region can be, at best, a crude approximation. Indeed, by definition, the critical exponents only capture the "universal" critical features of the system, not its full dynamics. Nevertheless, it is interesting to plot these quantities in a range of temperatures similar to the one that was used for the equilibrium thermodynamics quantities (for which, as we pointed out above, the predictions of our model are instead expected to extend to all temperatures below T c ), which allows one to highlight, in particular, the monotonically decreasing dependence of η/s as a function of the temperature for T ≤ T c , and the dramatic increase of ζ/s close to the critical point.

A. Discussion
In this work we derived the predictions of the statistical bootstrap model for thermodynamic quantities and transport coefficients near the critical endpoint of QCD. While it is well known that equilibrium thermodynamic quantities at temperatures below the (pseudo-)critical one are described well in terms of a gas of non-interacting hadrons, when all experimentally observed hadronic states with masses up to approximately 2 GeV [38,112] are included, the introduction of a continuous, Hagedorn-like, density of states for heavier states in the spectrum leads to the manifestation of critical behavior, without substantially altering the predictions for the equation of state at low temperatures. Moreover, as we remarked, the phenomenological implications of the model do not depend on the precise value of M , which in our computation was set to 2.25 GeV.
Even though the derivation of the critical exponents for this model is based on the assumption of a spectral density valid at zero chemical potential, and the dependence on µ B of the logarithm of the partition function is simply encoded in a fugacity factor, we argued that it may still capture the correct physics close to a possible critical point at finite µ B . For a continuous density of states of the form ρ(m) ∼ m a exp (bm) with a < −7/2, the energy and entropy densities remain finite even for point-like hadrons, while all higher-order derivatives diverge near T c . For a = −4 the energy density remains finite as T → T c , signalling the existence of a phase transition.
In passing, it is worth mentioning that a continuous spectral density of the form required for self-consistency in the statistical bootstrap model, eq. (17), also arises if one models hadrons in terms of confining strings of glue (as was done, for example, in ref. [113]).
Next, we studied the critical behavior of shear and bulk viscosities within the statistical bootstrap model. Identifying the thermodynamic quantities whose singular parts would contribute to the viscosity coefficients it is possible to write down Ansätze for the viscosity coefficients valid near the critical point. Using the Ansätze in eqs. (34) and (35) together with the singular part of the relevant thermodynamic quantities in eqs. (2)-(5), one can obtain the dominating contributions for the viscosity coefficients in eqs. (34) and (35). We found that the statistical bootstrap model predicts the shear viscosity to entropy density ratio η/s to decrease quite rapidly near T c . We observe that the magnitude of η/s away from the critical point is in good agreement with the predictions of non-critical models, and that there is a mild violation of the η/s ≥ 1/(4π) bound [104] near T c . It is worth emphasizing that this (slight) violation of the η/s ≥ 1/(4π) bound might be unphysical, i.e. an artifact of the model. In fact, it is also worth remarking that, while the conjecture of the η/s ≥ 1/(4π) bound was first derived in a holographic context [104] (and is expected to be saturated in strongly coupled gauge theories with a known gravity dual [114], such as the N = 4 supersymmetric Yang-Mills theory [115]), its origin is, in fact, much more general, being related to the uncertainty principle of quantum systems. As for the ζ/s ratio, we found it to rise very rapidly near T c , in remarkable agreement with refs. [60,61].
The anomalous behavior of shear and bulk viscosity coefficients near the critical endpoint might be very important for heavy-ion collision experiments. In particular, an enhanced bulk viscosity should manifest itself in heavy ion collisions through a decrease of the average transverse momentum of final state hadrons, and a corresponding increase in entropy. This feature may be particularly important for the experimental search for the critical endpoint of QCD through the BES program.
It is interesting to discuss a comparison of our findings with other related works. In particular, a study similar to ours was recently reported in ref. [116], which also predicts the transport properties of hot QCD matter within a hadron resonance model, finding a very low shear viscosity to entropy density ratio near T H , in agreement with our results. In contrast to our work, however, the focus of that article is not on the behavior near criticality (where, as we have discussed in detail above, much information can be derived with purely analytical calculations and general universality arguments), but instead on their numerical determination in a wider range, using the Gießen Boltzmann-Uehling-Uhlenbeck transport model [117] and Monte Carlo calculations. The fact that the numerical approach used in ref. [116] reproduces our analytical results close to criticality is a non-trivial cross-check of the results.
We should emphasize again that there is no fundamental proof that the statistical bootstrap model described in this work should necessarily provide a complete description of the fundamental properties at the CEP. Based on very general arguments (including the continuous nature of the transition at the critical endpoint, spacetime dimensionality, and the underlying symmetriesor lack thereof-of the theory), one may instead argue that the universality class of the critical endpoint of QCD should instead be the one of the Ising model in three dimensions (for a review, see ref. [118]). The critical exponents in this model have recently been computed to very high precisions using conformal bootstrap techniques [119,120], findingα = 0.11008(1),γ = 1.237075(10), etc., which are clearly incompatible with those predicted by the model that we considered here. If the critical endpoint of QCD exists, it may well be that its actual critical exponents are those of the three-dimensional Ising model, and that deviations from the description in terms of the statistical bootstrap model start to occur when one approaches the CEP. In this respect, it would be interesting to study theoretically how these deviations start to manifest themselves when the system is off, but close to, the critical point-perhaps using the analytical tools of conformal perturbation theory [121] (see also ref. [122], for an explicit example of application), as recently proposed in ref. [123]. An interesting issue associated with the description of the critical endpoint of QCD in terms of the three-dimensional Ising model concerns the identification of the lines, in the QCD phase diagram, that describe relevant deformations of the model (see also ref. [124]): what are the directions that correspond to a "thermal" and to a "magnetic" perturbation of the critical point? How do they affect the reliability of the description of the thermodynamics of the hadronic phase in terms of a hadron resonance gas with a spectrum of the form in eq. (7)? The answers to these questions may have important phenomenological implications for the evolution of the medium in energy scans going through or close to the critical endpoint, since they could directly affect the dynamics of hadrons before freeze-out.
Finally, it should be noted that, by construction, the statistical bootstrap model does not allow one to describe the approach to the critical endpoint of QCD from "above", i.e. from the deconfined phase. Perturbative computations show that the η/s ratio is generally large for a weakly coupled quark-gluon plasma [125,126] (a seemingly counter-intuitive result, which, in fact, reflects the fact that suppression of interactions makes "transverse" propagation of momentum difficult), but it is well known that thermal weak-coupling expansions are affected by non-trivial divergences, which are not present at T = 0 (see, for example, ref. [127] and references therein), and require a sophisticated treatment [128][129][130][131].

B. Conclusions
To summarize, in this work we derived the theoretical predictions of the statistical bootstrap model in the vicinity of the critical endpoint of QCD. Working in the approximation in which the critical temperature does not depend on the value of the chemical potential (which, as we remarked in section II B, has support from lattice QCD calculations showing that the crossover line in the phase diagram of the theory has very small curvature [89][90][91][92]), we showed that, for a suitable choice of its parameters, this model "predicts" the existence of a phase transition, and allows one to derive the associated critical exponents in an analytical way. Moreover, the model also gives predictions for the transport properties near the CEP, which are encoded in the shear and bulk viscosities. Both for the equation of state and for these transport coefficients, the dependence of the predictions of the model on the parameter a (within the rather narrow interval of physical interest) is very mild. In spite of the relative simplicity of the model, these results are qualitatively and quantitatively very similar to those obtained from other approaches. These findings may hopefully guide the future experimental identification of the CEP in heavy-ion collision experiments and the determination of the physical properties of QCD matter in the proximity of the critical endpoint.