Updates on the QCD phase diagram from lattice

Different aspects of the phase diagram of strongly interacting matter described by quantum chromodynamics (QCD), which have emerged from the recent studies using lattice gauge theory techniques, are discussed. A special emphasis is given on understanding the role of the anomalous axial U(1) symmetry in determining the order of the finite temperature chiral phase transition in QCD with two massless quark flavors and tracing its origin to the topological properties of the QCD vacuum.


Introduction
Understanding the phase diagram of strongly interacting hadronic matter described by quantum chromodynamics (QCD) is of fundamental importance, as it enables us to explain the origin of mass of 99.9% of the visible matter in the present universe. It has motivated large scale experiments from the Large Hadron Collider at CERN to the Relativistic Heavy Ion Collider at Brookhaven. Unraveling the nature and phases of QCD is one of the most challenging problems in modern theoretical physics, since much of it is driven by non-perturbative interactions. Lattice QCD is one of the most successful non-perturbative methods available to us. Using lattice techniques, it has been conclusively demonstrated that the phase transition in 2 + 1-flavor QCD, at vanishingly small baryon densities that may have occurred in the early epoch of cosmological evolution of the universe, is a smooth crossover [1][2][3]. Bulk thermodynamic quantities like entropy density, pressure and the equation of state (EoS) at vanishingly small baryon densities are now known to a very high precision [4,5] and more recently the continuum estimates for the EoS available at baryon densities as large as μ B /T ∼ 2.5 [6]. Lattice techniques are now becoming sophisticated enough to provide fundamental insights on the more microscopic details of different phases of QCD, in particular, on the origins of chiral symmetry The Institute of Mathematical Sciences, HBNI, Chennai 600113 India breaking and confinement that drive the phase transition in QCD. Understanding the mechanism behind the confinement of color degrees of freedom in gauge theories is one the most intriguing problems of contemporary research and any progress made towards understanding this phenomenon will definitely be an important milestone.
In this review article, updates and latest results on the QCD phase diagram at vanishingly small baryon densities, from very recent lattice studies will be discussed. Different phases of QCD as a function of the quark masses have been summarized in the famous Columbia plot. I will discuss about the current status of the Columbia plot and how a detailed understanding of it can in turn enrich our understanding of the phase diagram of QCD with physical quark masses. Moreover, I will further explain how the different phases in the Columbia plot provide us with insights about the intimate connection between chiral symmetry breaking and confinement. Both these phenomena are believed driven by topological properties of QCD, thus providing us with further insights about the vacuum structure of QCD. In recent years, the Columbia plot is studied on the lattice including a new axis to it, defined by an imaginary quark chemical potential. I will conclude this article with the new results coming from these studies. Throughout this article, the connection of the phase transition in QCD to the topological properties of the vacuum will be stressed.

The status of the Columbia plot of QCD at μ B = 0
The up and down quark masses in QCD are much lighter than its intrinsic scale, i.e m l = m u,d << QCD ; hence, the U L (2) × U R (2) symmetry of the QCD action is very mildly broken. (1) and QCD with two light and one heavier strange quark flavors, has to a very good approximation, a SU (2) This symmetry is broken to SU(2) V × U B (1) as the temperatures are lowered below ∼ QCD , leading to chiral symmetry breaking. The anomalous U A (1) part is always broken due to quantum effects. Though the chiral symmetry is exact in the limit m u , m d → 0, however, remnants of it exist in chiral observables. Thus even for QCD with physical quark masses, the temperature at which an inflection point exist for the subtracted chiral condensate is consistent with the one at which the chiral susceptibility or its disconnected part shows a maxima. An unweighted average of these temperature estimates in the continuum limit has allowed for a very precise determination of the pseudo-critical temperature T c = 156.5 ± 1.5 MeV [7], to within a 1% precision.
The U A (1) part of the chiral/flavor symmetry is an anomalous symmetry; thus, there is no corresponding order parameter. From renormalization group studies of model quantum field theories with the same flavor symmetries as QCD, it has been observed that the order of phase transition for two flavor QCD depends on the effective magnitude of the U A (1) breaking at T c [8]. Further studies using the epsilon expansion and renormalization group [9] as well as the conformal bootstrap techniques [10] have revealed a possibility of a first-order or even a second-order phase transition of is effectively restored near T c . This is in contrast to an O(4) second-order transition, if it remains broken. Magnitude of the effective breaking of U A (1), can only be determined non-perturbatively. Lattice techniques can and have contributed towards a more systematic understanding on this issue.
Before elaborating on it further, I would like to explain its consequences for the QCD phase diagram in the context of the Columbia plot. The current status of the Columbia plot is shown in Fig. 1, taken from Ref. [11]. The plot summarizes the fate of chiral phase transitions when Fig. 1 The current status of Columbia plot from various lattice studies, most recent of which are shown as arrows on the diagram, taken from reference [11] the masses of the light and strange quark flavors are varied. QCD with physical quark masses lie in the crossover region which is extended over a range of m u , m s . The upper right corner of the plot is now well understood since for quark masses infinitely large, it corresponds to a SU (3) gauge theory which has a first order transition [12]. This first-order region is separated from the crossover region by a Z(2) second-order line. The lower left corner of the plot is yet to be understood comprehensively. From model quantum field theories with same symmetries as N f = 3 QCD, it is expected that a first-order region exist which should again be separated from the crossover region by a second-order Z(2) line. Scaling studies of chiral susceptibilities along the very left part of the diagonal (which represents m s = m u,d ) on N τ = 6 lattices with a particular fermion discretization called the highly improved staggered quarks (HISQ), constrain the Z(2) line to exist for pion masses m π < 50 MeV [13]. With yet another fermion discretization on the lattice known as clover improved Wilson fermions, the corresponding critical pion mass has been measured to be m π < 170 MeV. This study at present is performed on rather coarse lattices [14]. However, it is important to note that irrespective of the choice of the fermion discretization on the lattice, the firstorder region tends to shrink when the lattice spacings are made finer [15], i.e., when one approaches the continuum limit. This picture seems to be realized also for N f = 4 QCD [16].
The second important issue that is still remains unresolved is whether this first-order region in the lower left corner of the plot ends in a tri-critical point in the chiral limit, i.e., m u,d = 0 and a finite m s or continues all the way to the m s → ∞ axis. Which of these two scenarios are realized in the continuum is ultimately related to the effective magnitude of the anomalous U A (1) symmetry. If U A (1) is effectively restored, the first-order region in the lower left hand corner of the plot may extend as a tiny strip parallel to the m s axis all the way to m s → ∞, separated from the crossover region by the Z(2) line. On this Z(2) line, the values of m u,d is much smaller than physical quark mass. To introduce the readers to the latest developments, arrows are marked on the Columbia plot in Fig. 1, to indicate the directions which are pursued by the most recent lattice studies in this regard. In summary, • The green arrow shows the approach for N f = 2 QCD, followed by one of the recent lattice studies [17,18]. With the current lattice volume (2.4 fm) 3 and inverse lattice spacing a −1 = 2.6 GeV, the results seem to suggest that U A (1) is restored at ∼ 1.1 T c for m u,d 5 MeV. On the other hand, expectations from N f = 3 QCD suggest that in the continuum, the first-order region, if it survives and continues from the lower left corner to the N f = 2 axis, will very narrow, characterized by m u,d << 5 MeV. It is thus a challenge to reconcile both these results; perhaps, it will be resolved in the continuum limit. • The blue line on the Columbia plot shows the other approach that was pursued by an independent lattice study [19], where the strange quark mass m s is fixed to its physical value and m u,d is successively reduced to check whether one approaches the Z(2) line or an O(4) second-order line. New results on chiral susceptibility for N τ = 8, 12 lattices with HISQ fermions suggest that its peak height decreases with increasing lattice volumes, ruling out first-order phase transition for M π > 80 MeV. Scaling studies of the chiral condensate seems to rule out Z(2) scaling for M π > 55 MeV, implying that the Z(2) line if it exists would be further left towards the m u,d = 0 axis of the Columbia plot. • Yet, another approach pursued was to measure the eigenvalue distribution of the QCD Dirac operator and understand the fate of U A (1) from the characteristics of the eigenvalue spectrum. A recent lattice study has also followed along the path marked by the blue arrow on the Columbia plot, to verify if U A (1) remains broken as light quark mass is successively reduced from its physical value [20]. If indeed U A (1) remains strongly broken, one will directly encounter the O(4) line instead of a Z(2) second-order line when moving towards the chiral limit m u,d → 0. The eigenvalue densities as observed different pion masses M π ∼ 160, 140, 110 MeV seem to support the existence of the O(4) line [20]. • A recent study was performed by considering N f as a continuous parameter and study the fate of the chiral phase transition as a function of N f [21]. Starting with the first-order region corresponding to N f = 3 QCD with finite quark masses, one can zoom in to the tri-critical scaling regime to constrain the minimum number of light quark flavors for which the tri-critical line exist, N tric f . Using the scaling relation , this present study on N τ = 4 lattices concludes that the N tric f < 2, which seems to suggest a first-order transition for two flavor QCD. These results are being further verified in the continuum limit.
Emphasizing again that the U A (1) is not the symmetry of the QCD partition function, thus, there is no corresponding order parameter. However, one can indirectly quantify its effective magnitude from different observables. One of the earliest suggestions was to look for the degeneracy of the integrated two-point correlation functions of iso-triplet pseudo-scalar and scalar mesons [22] in the chiral symmetry restored phase. In two flavor QCD, the possible meson correlation functions are: Here, the meson operators are defined as , τ being the isospin operator and ψ l denotes the field operator corresponding to the two degenerate light quark flavors. When the chiral symmetry is restored then χ π = χ σ and similarly χ η = χ δ . If one further demands anomalous U A (1) restoration, then it would result in χ π = χ η which due to chiral symmetry restoration would result in χ π = χ δ . Thus the difference between the integrated correlation functions χ π − χ δ is sensitive to the U A (1).
The integrated correlators can in turn be written in terms of the eigenvalues λ and density ρ(λ) of the QCD Dirac operator since χ π − χ δ = dλ Therefore, typical characteristics of the eigenvalue spectrum as a function of temperature can provide us with information about the fate of the U A (1) [23]. One way to trivially realize U A (1) restoration along with the chiral symmetry is to have a gap in the infrared region of the eigenvalue spectrum, i.e., ρ(λ → 0) = 0 [24]. A recent theoretical study has suggested that is important to look at even higher order correlation functions in all these mesonic quantum number channels [25]. In the chiral limit of two-flavor QCD, U A (1) breaking effects can be invisible in up to 6-point correlation functions in the scalarpseudo-scalar channel if the eigenvalue density goes as ρ(λ) ∼ λ 3 [25].
Recent measurements of χ π −χ δ in two flavor QCD with fermions with nearly exact chiral symmetry on the lattice (along the direction of the green arrow in the Columbia plot) report that U A (1) may be restored at about 1.1 T c [18]. Current systematics are still large to make any predictions closer to T c . For physical quark masses, the U A (1) breaking is still finite on lattices of size 48 3 × 12 which tend to decrease to zero in the limit m l → 0 [17,18]. On the other hand, pursuing along the blue arrow on the Columbia plot, the eigenvalue density of the Dirac operator in 2 + 1 flavor QCD (by fixing the strange quark mass m s to its physical value and successively reducing the light quark masses towards the chiral limit) does not show any trends towards opening up of a gap in the infrared region [20]. In fact, the infrared part of the spectrum has a similar behavior ρ(λ) ∼ λ at around 1.1 T c , when the light quark masses are reduced from m l = m s /20 [26] to m s /80. Even though for the given lattice sizes studied, the QCD Dirac operator has millions of eigenvalues; the first ∼ 100 of them contribute maximally to the U A (1) breaking observable, described earlier. An appropriately renormalized observable m 2 s (χ π − χ δ )/T 4 , reported in the same work, remains finite even for light quark mass m l = m s /40 at temperatures up to 1.1 T c . This suggests U A (1) breaking may as well survive towards the chiral limit [20]. An observable that measures the localized topological fluctuations of QCD vacuum at any temperature T is the topological susceptibility defined as χ t = T V Q 2 . Here, Q is the topological charge of a typical gauge configuration and the angular bracket denotes quantum averaging and also over the Euclidean space-time volume. The topological susceptibility can be further related to the disconnected part of the meson correlators χ η introduced earlier in Eq. 1, We recall here that the susceptibilities in the different quantum number channels are related as: where the π and δ-meson correlators only have connected parts. In the chiral symmetry restored phase of two-flavor QCD, Ward identities thus ensure that χ t = 1 2 m 2 l (χ π − χ δ ). This identity points to the fact that the axial U(1) symmetry breaking or its restoration as a function of temperature is of topological origin. This identity has been verified to hold even for 2+1 flavor QCD with nearly physical quark masses [27], showing that χ t is an observable that can be used to determine the extent of effective breaking of the anomalous U(1) symmetry in the chiral symmetry restored phase. The temperature dependence of the continuum extrapolated values for χ t in QCD matches exactly with the expectations from a dilute instanton gas approximation (DIGA) for temperatures T 3 T c [27][28][29][30], whereas a non-trivial temperature dependence is observed for T c < T < 3 T c [27,31]. New results for χ t with a different fermion discretization in 2 + 1 + 1 QCD have also confirmed this overall picture [32].
Just above T c , the temperature dependence for χ t is unambiguously different from the expectations of DIGA [27,31], the origin of which is not yet understood. This reemphasizes the fact that the axial U(1) symmetry is not effectively restored above T c , and thus, the QCD vacuum does not immediately go over to a trivial dilute instanton gas regime. Can this give us some hints about the topological origin of confinement? We recall that at finite temperatures, the eigenvalues of Polyakov loop at spatial infinity (the holonomy values) characterize the instanton solution. For trivial holonomy, the finite action solution for the non-Abelian gauge fields at non-zero temperatures is known and termed as caloron [33]. Calorons with non-trivial holonomy [34,35] in SU(N) gauge theory consists of N constituent instanton-dyons or simply dyons, which carry a fraction 1/N of the net topological charge of the caloron. Moreover, dyons carry both color electric and magnetic charges. Whereas calorons with trivial holonomy cannot explain confinement in gauge theories, mean-field studies of dyon gas hints to the fact that they interact with the holonomy potential and thus driving towards confinement [36].
Earlier lattice studies [37,38] have reported on the observation of dyons in pure gauge theory as well as in QCD respectively. On the lattice, the gauge ensembles are generated during a Monte-Carlo evolution with (anti)periodic boundary conditions along the temporal direction for (fermion) gauge fields respectively; thus, an isolated dyon which carries a non-trivial chromo-magnetic charge cannot exist. However, zero modes of a probe or valence Dirac operator, with generalized temporal periodicity phase ψ(τ +β) = e iφ ψ(τ ), will detect the dyon whose action is characterized by the difference between those eigenvalues of the Polyakov loop within which the phase φ lies. Recently, this technique has been utilized to not only detect dyons but also identify the different species of dyons in the temperature range T c < T < 1.2 T c [39,40]. The density profiles of zero-mode wave functions of the valence QCD Dirac operator and their characteristic fall-off at large distances have been suggested as one of the observables to identify different species of dyons on the lattice. Moreover, the eigenvalue spectrum of the nearzero modes of the valence Dirac operator with different periodicity phases [40] show distinct features corresponding to the different species of dyons. To illustrate this point, we refer to Fig. 2 from [40]. The density of near-zero eigenvalues corresponding to the L-dyons (anti-periodic boundary condition or φ = π) depletes at T ∼ 1.1 T c , signaling the restoration of chiral symmetry. However, the infrared part of the eigenvalue density corresponding to M-dyons (φ = ±π/3) is significantly larger than the Ldyons, which would imply presence of a large number of closely situated M-dyon pairs at high temperatures as compared to the L-dyons.
At higher temperatures T > 2 T c , the holonomy is trivial but there may be localized fluctuations of the Polyakov loop, which has been conjectured to provide the disordered landscape required to localize bulk eigenfunctions of the QCD Dirac operator [41]. Such localization of the Dirac eigenvalues have been reported in a recent lattice study in 2+1+1 flavor QCD [42].

Adding a new dimension to the Columbia plot
Till now, we have been discussing the conventional twodimensional Columbia plot. Imposing an additional third dimension in a form of an imaginary chemical potential [43] iμ q can impose further constraints on the conventional Columbia plot. The QCD partition function in presence of iμ q is free from the infamous sign-problem and has the property that Z( where n ∈ Z. The center of the gauge group is a symmetry of the pure gauge theory partition function but remains broken when quark fields are introduced. However, in the presence of an imaginary quark chemical potential, the center symmetry is a good symmetry of QCD with finite quark masses again. The phase of the Polyakov loop is an observable in this case, which will identify the different Z(3) sectors as one varies Fig. 2 The eigenvalue distribution of a valence Dirac operator (overlap operator) with temporal periodicity phases π, −π/3, π/3 respectively on 2+1 flavor QCD gauge ensembles. It is evident from the plot that whereas the chiral condensate in QCD due to the L-dyons (corresponding to the usual anti-periodic boundary condition given by φ = π) vanishes at around 1.1 T c , the contributions from M-dyons survive. This figure has been taken from Ref. [40] iμ q . The Roberge-Weiss (RW) points characterized by μ q = (2n + 1) iπT 3 [44] represent the transition between adjacent center sectors, which is first order for high temperatures and a smooth crossover for lower temperatures. From the continuity of the free energy, the first-order lines should end in a second-order RW end-point.
How do the deconfinement and chiral transitions at μ q = 0 connect to the RW points? From reflection symmetry of the partition function, the first-order lines for QCD with very heavy quark masses m q are expected to meet at the RW point, which will be a triple point. For intermediate values of m q , the crossover curve at μ q = 0 may meet at the RW end-point, expected to be in the Z(2) universality class. In the chiral limit, the scenario is still far from being understood. Numerical simulations on relatively small N τ = 4 lattices with staggered fermions have shown a first order RW transition for both N f = 2, 3 [43,45], which seems to survive in the chiral limit [46]. This scenario was later confirmed in studies with Wilson fermions [47]. The discussion in this section is succinctly summarized in the modified Columbia plot at μ q = μ B /3 = iπT/3, shown in Fig. 3, taken from [11]. Now, I would like to elaborate on the consequences of these results obtained in the imaginary chemical potential plane to the nature of the N f = 2 chiral transition at μ q = 0. If the N f = 2 chiral transition at iμ q = 0 is second (first) order, then the first-order RW transition will end in a tri-critical point for μ 2 q < 0 (μ 2 q > 0) respectively. This opens up a new possibility to understand the nature of chiral phase transition and the role of U A (1) for N f = 2 QCD. Indeed, the first lattice study along this line [46], performed with staggered fermions on N τ = 4 for different lattice volumes with N s = 8, 12, 16, reported a tri-critical point at μ 2 q = 0.85(5)T 2 which seemed to suggest that N f = 2 chiral transition at μ q = 0 is first order at least on coarser lattices. Subsequently improved versions of staggered fermions have been used to reduce lattice cutoff effects. By fixing the m s to its physical value and reducing the m l = m u,d at μ q = iπT/3, i.e., along the blue line shown on the lower RW plane of Fig. 3, the recent results [48] report that a first-order RW transition is not observed for M π ≥ 50 MeV. In the continuum limit, it would imply that this first-order region, when it continues to the μ q = 0 plane would be a very narrow strip parallel to the m s axis. Remarkably, the chiral and RW transition seem to follow each other as one reduces the m u,d [48].
Furthermore, if the end-point of the line of first-order RW transitions is of second order, then it would belong to the three-dimensional Ising universality class. The imaginary part of the Polyakov loop L would behave like the Fig. 3 Current status of the Columbia plot with an additional direction added to it which is characterized by the imaginary baryon chemical potential (which is a third of the quark chemical potential μ q ) normalized by the temperature, μ B /T, taken from [11]. The lower plane is the current understanding we have of the generalized Columbia plot at μ B /T = 3μ q /T = iπ Fig. 4 The Z(2) scaling of the expectation value of the imaginary part of Polyakov loop near the chiral crossover transition temperature T c for different lattice volumes, taken from [49] magnetization order parameter of the Ising model. This is due to the fact that under Z(2) transformation, the real part of Polyakov loop does not change sign whereas its imaginary part does. Scaling studies performed [49] with the expectation values of |ImL| , for different lattice volumes, N τ = 4 and N σ = 8−24 around the chiral crossover transition temperature T c ∼ 200 MeV and M π = 135−90 MeV and on the 3μ q /T = iπ plane, confirm a secondorder Z(2) scaling, essentially supporting the previous findings ( Fig. 4) [48].

Summary and outlook
In this review article, the most recent findings about the QCD phase diagram arising from lattice studies is outlined in the context of the Columbia plot. Very new insights about the order of the chiral phase transition for N f = 2 QCD has been obtained both from the study of its eigenvalue spectrum (and integrated meson correlators) and indirectly from the scaling studies of imaginary part of the Polyakov loop expectation values at an imaginary baryon chemical potential μ B = iπT. Though most studies argue for a second-order O(4) phase transition, however, no final conclusion has yet been achieved. The order of chiral phase transition at N f = 2 is intimately connected to the effective magnitude of the anomalous axial U(1) symmetry breaking. It has been stressed in this article that the origin of this effective breaking of axial U(1) is related to the topological fluctuations of the QCD vacuum hence to quantify it is one of the most challenging problems in lattice field theory. Ultimately, a comprehensive understanding of this phenomenon will in turn advance our knowledge about the mechanism of confinement and chiral transition in QCD with physical quark masses.