Thermodynamics and susceptibilities of isospin imbalanced QCD matter

TL;DR

This paper investigates the thermodynamic properties and susceptibilities of isospin imbalanced quark matter using the NJL model, revealing phase transition behaviors and agreement with lattice and chiral perturbation theory data.

Contribution

It provides new insights into the phase transition dynamics and susceptibility behaviors of isospin imbalanced QCD matter within the NJL model framework.

Findings

Susceptibilities peak at the critical isospin chemical potential.

Temperature smooths the phase transition.

Pion susceptibility is zero in the normal phase and finite in the superfluid phase.

Abstract

We study the thermodynamics and the susceptibilities of quark matter in the framework of two-flavor Nambu-Jona-Lasinio model at finite isospin chemical potential and temperature. Isospin number density, normalized energy density and trace anomaly are shown to be in good agreement with the available lattice data as well as with the results from chiral perturbation theory at zero temperature. We also study how susceptibilities depend on the isospin chemical potential and on temperature. We find a peak for the chiral, pion, and isospin susceptibilities at the critical isospin chemical potential, $\mu_I^c(T)$, at the boundary of the phase transition between the normal and pion superfluid phase. Moreover, temperature makes the transition from normal to pion condended phase smoother. We also note that the pion susceptibility always remains zero in the normal phase while it is finite in the superfluid phase.

Equations

$$\mathcal{L}=\bar{q}(i{\gamma }^{\mu }{\partial }_{\mu }-m+\hat{\mu }{\gamma }_0)q+G[(\bar{q}q{)}^2+(\bar{q}i{\gamma }_5\tau q{)}^2],$$

$$(\bar{q}q{)}^2\approx 2(\bar{q}q)⟨\bar{q}q⟩-⟨\bar{q}q{⟩}^2,$$

$$(\bar{q}i{\gamma }_5\tau q{)}^2\approx 2(\bar{q}i{\gamma }_5\tau q)⟨\bar{q}i{\gamma }_5\tau q⟩-⟨\bar{q}i{\gamma }_5\tau q{⟩}^2.$$

$$⟨\bar{q}q⟩=\sigma ={\sigma }_u+{\sigma }_d,$$

$${\pi }^{\pm }=⟨\bar{q}i{\gamma }_5{\tau }_{\pm }q⟩\equiv \Pi e^{\pm i\theta },$$

$$\Pi =⟨\bar{q}i{\gamma }_5{\tau }_{1}q⟩.$$

$$\mathcal{L}$$

$$\Omega$$

$$E_p$$

$$E_p^{\pm }$$

$$\frac{\partial \Omega }{\partial \sigma }=\frac{\partial \Omega }{\partial \Pi }=0.$$

$$n_I$$

$$P$$

$$E$$

$$P^{\text{LO}}$$

$$E^{\text{LO}}$$

$$n_I=\frac{\partial P^{\text{LO}}}{\partial {\mu }_I}=\frac{f_{\pi }^2}{{\mu }_I^3}({\mu }_I^4-m_{\pi }^4).$$

$${\mu }_I^{\text{peak}}\simeq \{\begin{array}{ll}1.274m_{\pi },& \text{NJL},\\ 1.276m_{\pi },& \text{CHPT}\text{cite[cite]@@bibrefAuthors Phrase1YearPhrase216Carignano.Mammarella.ea51503-51503PRD@@citephrase(@@citephrase)},\\ \{1.20,1.25,1.275\}{m}_{\pi },& \text{Lattice data}\text{cite[cite]@@bibrefAuthors Phrase1YearPhrase212Detmold.Orginos.ea54507-54507PRD@@citephrase(@@citephrase)}.\end{array}$$

$${\bar{\mu }}_I\simeq 1.754m_{\pi },$$

$${\chi }_{\sigma }=-\frac{\partial \sigma }{\partial m}.$$

$${\chi }_{\pi }=\frac{\partial \Pi }{\partial m}.$$

$${\chi }_I=\frac{d{n}_I}{d{\mu }_I}.$$

Full text

I Introduction

There is compelling evidence that Quantum Chromodynamics (QCD) is the correct theory describing the strong interaction between quarks and gluons. QCD is a non-Abelian gauge theory, which has a nontrivial vacuum state and provides a rich phase structure and exciting physical phenomena in various environments Andersen, Naylor, and Tranberg (2016); Kharzeev (2015); Miransky and Shovkovy (2015); Philipsen (2019). At low energy densities quarks are confined inside hadrons, while they form a weakly coupled plasma in the very high energy limit. Therefore, the strongly interacting QCD matter is expected to undergo a hadron-quark phase transition at high temperature and/or chemical potential. Moreover, the chiral symmetry is spontaneously broken at low temperature/baryon density due to the nonvanishing chiral condensate, while it is expected to restore at high temperature and high density Stephanov (2006); Fukushima and Hatsuda (2011). Lattice calculations show that at small baryon chemical potential the transition is actually a smooth crossover Borsanyi et al. (2010a); Bazavov et al. (2012); Borsanyi et al. (2010b, 2014); Cheng et al. (2010). On the other hand, many studies Alford, Rajagopal, and Wilczek (1998); Halasz et al. (1998); Costa, Ruivo, and de Sousa (2008) suggest that at large chemical potential the smooth crossover meets a first order transition at a particular point in the $(T,\mu)$ plane named the critical endpoint (CEP), where the crossover becomes a second order phase transition. The quest for the possible existence of CEP as well as its precise location are still on-going, and have attracted a lot of attention, see e.g. Ferreira, Costa, and Providência (2018); Qin et al. (2011); Costa et al. (2014); Ruggieri et al. (2014); Costa et al. (2015); Lu et al. (2015); Ruggieri (2011).

In QCD, fluctuations of the conserved charges like the baryon and isospin numbers, exhibit critical behavior around a phase transition; these fluctuations are encoded into the corresponding susceptibilities as well as into higher order moments Jeon and Koch (2000); Asakawa, Heinz, and Muller (2000); Hatta and Ikeda (2003); Hatta and Stephanov (2003). For example, the chiral susceptibility, which represents the modification of the chiral condensate to a small perturbation of the current quark mass, is often used to describe the chiral phase transition Wang, Sun, and Zong (2013); Min et al. (2008); Chakraborty, Mustafa, and Thoma (2003); Zhao et al. (2008); Bernard et al. (2005); Aoki et al. (2006); Cheng et al. (2007). Besides, it has been suggested that the moments of the baryon number can be used to identify the CEP in the QCD phase diagram Harada, Sasaki, and Takemoto (2010); Ding et al. (2015); Fujii (2003); Fujii and Ohtani (2004); Kunihiro (1991); Stephanov, Rajagopal, and Shuryak (1998); Gavai and Gupta (2001); Gavai, Gupta, and Majumdar (2002); Cui, Takeuchi, and Wu (2011); Jiang, Luo, and Zong (2011); He et al. (2009); Ratti, Roessner, and Weise (2007).

In the two-flavor case the quark chemical potentials $\mu_{f}$ can be expressed in term of the baryon chemical potential $\mu_{B}=N_{c}(\mu_{u}+\mu_{d})/2$, and the isospin chemical potential $\mu_{I}=\mu_{u}-\mu_{d}$111Note that another definition of the isospin chemical potential $\mu_{I}=(\mu_{u}-\mu_{d})/2$ is also frequently used in the literature, resulting in the critical point as $\mu_{I}^{c}(T=0)=m_{\pi}/2$ (see, e.g., Refs. Brandt, Endrődi, and Schmalzbauer (2018); Adhikari and Andersen (2017)). with $N_{c}=3$ the number of colors. Along the $\mu_{B}$ axis, the knowledge from first principle lattice simulation on the thermodynamics of strongly interacting matter is limited, due to the notorious sign problem Muroya et al. (2003). Thus, effective field theories of QCD as well as phenomenological models are necessary, and they might provide powerful tools for us to have a better understanding of QCD physics in the nonperturbative regime. On the other hand, lattice simulations are feasible for $\mu_{B}=0$ and $\mu_{I}\neq 0$, therefore it is possibe to study QCD at finite $\mu_{I}$ by means of first principle calculations and compare these with the predictions of effective models. This is an useful exercise to test the effectiveness of QCD models. For example, an analysis based on the leading order (LO) Son and Stephanov (2001); Splittorff, Son, and Stephanov (2001) and the next-to-leading order (NLO) Adhikari, Andersen, and Kneschke (2019); Adhikari and Andersen (2019a) chiral perturbation theory (CHPT) shows that a second order phase transition from the normal to pion superfluid phase is expected to happen at $\mu_{I}^{c}(T=0)=m_{\pi}$, where $m_{\pi}$ is the pion mass. This has been confirmed by lattice simulations Kogut and Sinclair (2002a, b, 2004); Detmold, Orginos, and Shi (2012) as well as by the Nambu$-$Jona-Lasinio (NJL) model analytically He, Jin, and Zhuang (2005); He and Zhuang (2005); Warringa, Boer, and Andersen (2005).

Many aspects of pion condensate have been extensively studied in the effective field theories Loewe and Villavicencio (2003, 2005); Adhikari, Cohen, and Sakowitz (2015); Mammarella and Mannarelli (2015); Carignano et al. (2017); Carignano, Mammarella, and Mannarelli (2016); Adhikari (2019); Adhikari and Andersen (2019b), lattice simulations Brandt and Endrodi (2016); Shi (2013); Scior, von Smekal, and Smith (2018); Brandt, Endrodi, and Schmalzbauer (2018) and effective models Mao, Petropoulos, and Zhao (2006); Matsuzaki (2010); Ueda et al. (2013); Stiele, Fraga, and Schaffner-Bielich (2014); Nishihara and Harada (2014); Andersen and Kyllingstad (2009); Xia, He, and Zhuang (2013); Wu, Ping, and Zong (2017); Chao, Huang, and Radzhabov (2018); Xiong, Jin, and Li (2009); Abuki et al. (2009); Avancini et al. (2019) (for a review, see Mannarelli (2019)). We also mention that a new type of compact star made of a charged pion condensate has been proposed recently Mao (2014); Brandt et al. (2018); Andersen and Kneschke (2018). While we refer to the original articles for the detailed picture for both two-flavor and three-flavor quark matter, we remind here the well established picture for the two-flavor case with finite bare quark masses. At zero temperature there is a second order phase transition from the normal phase, in which the pion condensate is vanishing, to the phase with nonzero pion condensate: this phase transition happens when $|\mu_{I}|=m_{\pi}$; increasing $|\mu_{I}|$ there is eventually another second order phase transition back to the normal phase. At finite temperature the gap between the two critical $|\mu_{I}|$ becomes smaller and eventually the window for the pion condensate phase closes: at very high temperature there is room for the normal phase only, see for example Fig. 4 of Wu, Ping, and Zong (2017).

In this paper, we study several susceptibilities of isospin inbalanced QCD matter at vanishing $\mu_{B}$ and zero and/or finite temperature, which to the best of our knowledge has not yet been explicitly discussed in the literature. As mentioned above, doing this study is useful to test how effective models predictions compare with first principle calculations. The thermodynamic properties of pion superfluid phase at vanishing temperature has been studied at the lowest order of CHPT at zero temperature Carignano, Mammarella, and Mannarelli (2016) and it has been found a good agreement with the lattice QCD simulations for $\mu_{I}\lesssim 2m_{\pi}$, but also that the LO approximation breaks down at higher $\mu_{I}$. The NJL model can be employed in a larger range of chemical potential and temperature with respect to LO CHPT and it is the purpose of this work to show how this model describes the thermodynamics and the fluctuations for the transition to the pion superfluid phase. We will use the two-flavor NJL model in this article for simplicity, leaving the three-flavor case to a future study.

The article is organized as follows. In Sec. II, we give a brief introduction for the theoretical framework of NJL model at finite $\mu_{I}$ and temperature. In Sec. III, after comparing the NJL model results for several thermodynamic quantities with LO CHPT calculation and/or lattice data at vanishing temperature, we give the numerical results for several susceptibilities of our interest at finite $\mu_{I}$ and temperature. Finally, we draw our conclusions in Sec. IV

II The NJL model

The standard NJL model Lagrangian density for two flavors of quarks is given by Nambu and Jona-Lasinio (1961a, b); Vogl and Weise (1991); Klevansky (1992); Volkov (1993); Hatsuda and Kunihiro (1994)

[TABLE]

where $q$ represent the quark fields, $\tau$ the Pauli matrices in flavor space, $G$ the four-fermion interaction coupling, $m$ the degenerate quark mass of up and down quarks and $\hat{\mu}=\mathrm{diag}\{\mu_{u},\mu_{d}\}$ denotes the diagonal matrix for quark chemical potentials.

In the grand canonical ensemble the thermodynamics of strongly interacting matter can be considered as a function of the isospin and baryon chemical potentials and temperature. In this article we limit ourselves to the case $\mu_{B}=0$ which leads to $\mu_{u}=\mu_{I}/2$, $\mu_{d}=-\mu_{I}/2$. For simplicity, we perform the calculations in the mean field approximation, i.e.,

[TABLE]

We introduce the chiral condensates as

[TABLE]

with $\langle\bar{u}u\rangle=\sigma_{u}$ and $\langle\bar{d}d\rangle=\sigma_{d}$. At $\mu_{I}\neq 0$ a charged pion condensate is also expected: in order to account for this we introduce

[TABLE]

with $\theta$ a real number and $\tau_{\pm}=\tau_{1}\pm i\tau_{2}$; the effective potential does not depend on $\theta$ so we choose $\theta=0$ in the above equation and we are left with the pion condensates, $\Pi$, in the $\tau_{1}$-direction of isospin, namely

[TABLE]

The Lagrangian density can thus be written as

[TABLE]

with the effective quark mass defined as $M_{q}=m-2G\sigma$. The thermodynamic potential is

[TABLE]

where $\beta=1/T$ and

[TABLE]

The ground state is given by the values of $\sigma$, $\Pi$ that minimize $\Omega$: this is equivalent to require that $\sigma$ and $\Pi$ are solution of the gap equations,

[TABLE]

Given $\Omega$ at the ground state it is straightforward to obtain the isospin number density, the pressure and the energy density as

[TABLE]

where $\Omega_{0}$ is the thermodynamic potential at $T=\mu_{I}=0$.

The two-flavor NJL model has three parameters: the current quark mass $m=0.005$ GeV, the four-fermion coupling strength $G=5.01$ GeV*-2*, and the hard three-momentum cutoff $\Lambda=0.653$ GeV, which are fixed by reproducing the empirical values of the pion mass $m_{\pi}=0.134$ GeV, the physical pion decay constant $f_{\pi}=0.093$ GeV, and the chiral condensate $\sigma_{0}=2(-0.25~{}\text{GeV})^{3}$ in vacuum Zhuang, Hufner, and Klevansky (1994).

Before going on it is useful to remind the results of LO CHPT for the pressure and the energy density in the pion superfluid phase, namely

[TABLE]

where $f_{\pi}$ is the pion decay constant. The result for isospin density will be also useful in the following section:

[TABLE]

III Results and discussion

III.1 Isospin number density, energy density, and trace anomaly

In Fig. 1, we plot the isospin number density normalized to $m_{\pi}^{3}$ at zero temperature as a function of $\mu_{I}/m_{\pi}$. For comparison we also include the results from LO CHPT222Recently, the contribution from Lee-Huang-Yang term is also included Lepori and Mannarelli (2019). (purple dotted line) as well as recent lattice data (magenta circles) Brandt et al. (2018). The results obtained from these methods agree quantitatively with each other in the considered range of $\mu_{I}$. When $\mu_{I}<m_{\pi}$ the isospin number density is zero, which corresponds to the normal phase of the system. In the pion superfluid phase with $\mu_{I}>m_{\pi}$, the isospin number density becomes nonzero, which increases monotonically with $\mu_{I}$.

In Fig. 2, we plot the energy density normalized to $\mu_{I}^{4}$ as a function of $\mu_{I}/m_{\pi}$ at $T=0$. The results from the LO CHPT (purple dotted line) in Eq. (16) are also shown for comparison. The result from the NJL model agrees with the one from CHPT at zero temperature though the energy density from the former is slightly larger than that from the latter at high $\mu_{I}$. The quantity $E/\mu_{I}^{4}$ becomes nonzero and positive when $\mu_{I}$ exceeds the critical value $\mu_{I}^{c}(T=0)=m_{\pi}$. In particular, it develops a peak at $\mu_{I}^{\text{peak}}\simeq 1.274m_{\pi}$. Recently, a calculation of this peak based on the LO CHPT has been done in Ref. Carignano, Mammarella, and Mannarelli (2016); moreover, this quantity has also been computed on the lattice:

[TABLE]

In the above equation the three lattice results quoted have been obtained with three different lattice volumes, namely $L^{3}$=$\{16^{3}$, $20^{3}$, $24^{3}\}$ respectively. The extrapolation of the lattice results to the continuum limit gives $\mu_{I}^{\text{peak}}=1.30(7)m_{\pi}$. We notice the nice agreement of the NJL model with LO CHPT and lattice data.

In Fig. 2 we also show the normalized trace anomaly, $(E-3P)/\mu_{I}^{4}$, as a function of $\mu_{I}/m_{\pi}$ at $T=0$. Again we notice the agreement between the NJL model and the LO CHPT results in the entire isospin chemical potential range considered. In more detail, we observe that the blue dashed line (lower one), which represents the result calculated in the NJL model at zero temperature, keeps zero at small $\mu_{I}$, and becomes nonzero and positive when $\mu_{I}$ slightly larger than $m_{\pi}$. Similar to the behavior of $E/\mu_{I}^{4}$, the normalized trace anomaly also develops a peak at an intermediate $\mu_{I}$ although the peak of the former appears at a slightly smaller value of $\mu_{I}$ with respect to the latter. As $\mu_{I}$ increases inside the domain of the pion condensed phase, the normalized trace anomaly decreases and becomes negative. In fact, the point $\bar{\mu}_{I}$ fulfilling the conformal relation $E-3P=0$ can be evaluated in the NJL model at zero temperature as

[TABLE]

which is in good agreement with the LO CHPT result given by $\bar{\mu}_{I}=\sqrt{3}m_{\pi}$ Carignano, Mammarella, and Mannarelli (2016). This point separates the $E>3P$ and $E<3P$ regions. In fact, this value $\bar{\mu}_{I}\simeq 1.754m_{\pi}$ is very close to the value of $\mu_{I}$ corresponding to the BEC-BCS crossover, which is estimated to be $1.702m_{\pi}$ at zero temperature according to the analyses in Refs. He, Jin, and Zhuang (2006); Sun, He, and Zhuang (2007); He and Zhuang (2007); Abuki et al. (2010). However, although these two values almost coincide with each other, whether $\bar{\mu}_{I}$ can be identified with the isospin chemical potential corresponding to the BEC-BCS crossover or not need to be investigated in more detail in the future.

III.2 Chiral, pion, and isospin susceptibilities

The transition to the pion condensate phase can be well identified by means of the susceptibilities. In principle, they can also be measured on the lattice since there is no sign problem at finite $\mu_{I}$ due to the real and positive fermionic determinant. Chiral condensate is the order parameter of spontaneous chiral symmetry breaking, while pion condensate indicates the spontaneous isospin symmetry breaking. The chiral susceptibility, $\chi_{\sigma}$, which corresponds to the zero-momentum projection of the scalar propagator and encodes all fluctuations of the order parameter, is defined as Karsch and Laermann (1994); Bernard et al. (2005); Chakraborty, Mustafa, and Thoma (2003); Aoki et al. (2006)

[TABLE]

We can also define the pion susceptibility, which is obtained as the first derivative of pion condensate with respect to the current quark mass, i.e.,

[TABLE]

We also define the isospin number susceptibility as

[TABLE]

Note that we write the derivative operator above as a total derivative instead of a partial one because one might take into account the fact that the condensates may have a dependence on $\mu_{I}$.

The chiral, pion, and isospin susceptibilities as functions of $\mu_{I}/m_{\pi}$ for four values of the temperature are illustrated in Fig. 3. The temperatures are chosen as $T=0$ (black short-dashed line), $T=0.1$ GeV (blue dashed line), $T=0.15$ GeV (red dotted line), and $T=0.22$ GeV (magenta dot-dashed line). We notice that the zero temperature susceptibilities exhibit a discontinuity at $\mu_{I}=\mu_{I}^{c}(T)$, signaling the boundary of a second order phase transition. At nonzero temperatures, e.g., $T=0.1$ GeV and 0.15 GeV, the situation is similar except that the discontinuity shifts to larger $\mu_{I}$. For $T=0.22$ GeV instead the susceptibilities are continuous in the whole range of $\mu_{I}$, which can be easily understood with the absence of a transition to the pion condensed phase in agreement with He, Jin, and Zhuang (2005). We notice that the chiral susceptibility stays finite in the pion condensed phase, while the pion susceptibility vanishes in the normal phase: this is obviously related to the fact that the bare quark mass is finite, so chiral symmetry is always broken explicitly and this leads to a finite chiral condensate. In the chiral limit this would not happen and the chiral susceptibility would vanish in the pion condensed phase.

In Fig. 4, we show the thermal behavior of the chiral susceptibility at several values of $\mu_{I}$. The chiral susceptibility at $\mu_{I}=0$, which is denoted by the black short-dashed line, firstly grows up smoothly at low temperature, then it develops a smooth peak at $T=0.2$ GeV in correspondence of the chiral crossover. At higher temperatures, the chiral susceptibility decreases with temperature. We also include several cases with $\mu_{I}>\mu_{I}^{c}(T)$ for comparison. For $\mu_{I}=0.3$ GeV, which is represented by the blue dashed line, the chiral susceptibility has a much larger value at $T=0$ compared to the other cases. This line first slowly increases with temperature and reaches a smooth peak at $\mu_{I}/m_{\pi}=0.175$, signaling a crossover. After the peak the line drops and show a kink at $\mu_{I}/m_{\pi}=0.183$, in correspondence of the transition to the normal phase. Similarly, the lines for $\mu_{I}=0.5$ GeV and 0.7 GeV, which are represented by red dotted and orange dot-dashed lines respectively, also show a modest kink peak at the critical temperature.

In Fig. 5, we plot $\chi_{\pi}$ versus temperature for several values of $\mu_{I}$, i.e., $\mu_{I}=0.2,~{}0.3,~{}0.5,~{}0.7$ GeV. We do not illustrate the result for $\chi_{\pi}$ with $\mu_{I}$ smaller than the critical $\mu_{I}^{c}(T)$ because in the normal phase $\chi_{\pi}=0$. For $\mu_{I}=0.2$ GeV, which is labeled by the black short-dashed line, $\chi_{\pi}$ increases very smoothly at low temperature, and then grows up rapidly when $T$ is close to the critical temperature; clearly it vanishes above this temperature because the pion condensate is zero there.

In Fig. 6 We plot $\chi_{I}$ versus the temperature for several values of $\mu_{I}$ . For the black short-dashed line, which represents the result obtained at $\mu_{I}=0$, the isospin number susceptibility increases smoothly and monotonously with temperature. At higher $\mu_{I}$, e.g., $\mu_{I}=0.3$ GeV (blue dashed line), $\mu_{I}=0.5$ GeV (red dotted line), and $\mu_{I}=0.7$ GeV (magenta dot-dashed line), $\chi_{I}$ keep almost unchanged in the entire considered temperature range expect for the cusp at the phase transition point.

IV Conclusions

In this work, we have studied the isospin chemical potential and temperature dependence of several thermodynamic quantities, in particular of several susceptibilities, of isospin-imbalanced QCD matter at vanishing baryon chemical potential within the two-flavor NJL model. We have compared the isospin number density at zero temperature with recent lattice data as well as with LO CHPT, finding a good agreement with the latter results, see Fig. 1. In addition, the result for the normalized energy density $E/\mu_{I}^{4}$ as a function of $\mu_{I}/m_{\pi}$ is also in fair agreement with that from LO CHPT at zero temperature. In particular, the location of the peak of $E/\mu_{I}^{4}$ as well as the trace anomaly computed within the NJL model agree with LO CHPT and/or lattice simulations: we have found $\mu_{I}^{\text{peak}}\simeq 1.274m_{\pi}$ for the peak of the energy density and $\bar{\mu}_{I}\simeq 1.754m_{\pi}$ for the condition $E-3P=0$.

We have then considered the chiral, pion and isospin susceptibilities at finite temperature. These are interesting quantities because they allow to define the critical lines in the $(T,\mu_{I})$ plane and can also be computed on the lattice, so they can be used as a test of the effective models of QCD. The qualitative behavior of these susceptibilities is in agreement with naive expectations: in particular, the chiral susceptibility is always finite because of the nonvanishing bare quark mass that induces a finite chiral condensate even at large $T/\mu_{I}$. On the other hand, the pion susceptibility vanishes in the normal phase and is nonzero in the pion condensed phase. The transition from the normal phase to the pion condensed phase is of the second order with a divergent pion susceptibility. The isospin number susceptibility also keep finite and almost unchanged in the normal phase except that it is zero for the zero temperature case, and shows a discontinuity at the phase transition.

In this study we have ignored the role of baryon and strangeness chemical potentials Carignano, Mammarella, and Mannarelli (2016). These external sources are known to play an important role and we leave the study of the more complete problem to future works.

Acknowledgement

The authors acknowledge Jingyi Chao and John Petrucci for inspiration and comments on the first version of this article. ZYL acknowledges the warm hospitality of the Department of Nuclear Science and Technology of Lanzhou University where part of this work was done. This work is supported in part by the China Postdoctoral Science Foundation under Grant No. 2017M620920, by DFG and NSFC through funds provided to the Sino-German CRC 110 “Symmetries and the Emergence of Structure in QCD” (NSFC Grant No. 11621131001 and DFG Grant No. TRR110), by the NSFC under Grant Nos. 11747601, 11835015, and 11705163, by the Chinese Academy of Sciences (CAS) under Grant No. QYZDB-SSW-SYS013 and No. XDPB09, and by the CAS Center for Excellence in Particle Physics (CCEPP). The work of M. R. is supported by the National Science Foundation of China (Grants No. 11805087 and No. 11875153) and by the Fundamental Research Funds for the Central Universities (Grants No. 862946).