Parity Doubling in QCD Thermodynamics †

: Motivated by the recent lattice study by FASTSUM collaboration, effective masses of the baryon parity-doublers are shown for various pion masses. A general trend of the nucleon and delta parity-doublers is consistent with the lattice Quantum Chromodynamics (QCD) observation, whereas the hyperon masses exhibit a qualitatively different behavior, traced back to the lattice set-up with the heavy pion comparable to the kaon. As an application to hot QCD, we demonstrate the ﬂuctuations and correlations involving baryon number in hot hadronic matter with modiﬁed masses of negative-parity baryons, in the context of the hadron resonance gas. Confronting the baryon number susceptibility, baryon–charge and baryon–strangeness correlations as well as their ratios with the lattice QCD data for the physical pion mass, we ﬁnd that the strong downward mass shift in the hyperons can accidentally reproduce some correlation ratios, however it also tends to overshoot the individual ﬂuctuations and correlations of lattice simulations. Another application of nucleon parity doubling is the physics of neutron stars. Under beta equilibrium and charge neutrality, hadronic matter with unbroken chiral symmetry can be favored in the core of the neutron stars.


Introduction
Spectral properties of hadrons are expected to change in a medium. A search for those modifications because of partial restoration of the chiral symmetry has been one of the central subjects in heavy-ion collisions where hot and/or dense matter is supposed to be created. Substantial medium modifications of the lowest-lying hadrons have been predicted in a large class of chiral models, whereas their reliable estimates in lattice QCD are thus far limited to the screening masses. Recently, the thermal masses of non-strange and strange baryons with positive and negative parity have been extracted from the temporal correlation functions by the FASTSUM collaboration [1,2]. The negative-parity states clearly show downward mass shifts, whereas the positive-parity baryons are rather insensitive to temperature. The obtained spectra follow qualitatively an expectation from parity doubling of the chiral symmetry, indicating that the states with the same spin but opposite parity tend to become degenerate when the symmetry gets partially restored.
In this contribution, we briefly summarize our recent study on the thermodynamics of parity doubling: the in-medium masses of baryon octet and decuplet and their potential influence over fluctuations and correlations in the context of a chiral effective theory [3] and in-medium Hadron Resonance Gas (HRG) model [4,5], as well as the physics of neutron stars in a hybrid approach for hadrons and quarks [6].

Parity Doubling
The baryon octet and decuplet with opposite parity can be modeled in an effective chiral approach with N f = 2 + 1 [3], where the masses are described as certain functions of the light-quark σ q and strange-quark σ s condensates as well as the degenerate mass of parity doublers m 0 , which is not anchored to dynamical chiral symmetry breaking. Their modifications are therefore driven entirely by the in-medium quark condensates. We use the thermal profiles of the condensates for the physical pion and kaon masses measured at 0.6 < T/T c < 1.2 in lattice QCD by the HotQCD collaboration [7].
In Figure 1, we present the results from this effective chiral model as well as the masses based on a parameterization motivated by the lattice results in [1,4].  Table 1.
for a given baryonic state i with negative parity, and with M i − (T c ) and b i the fitting parameters to be fixed for each channel. The data are given as a function of T/T c ; the value of T c does not affect the fit but later it is set to the physical value of 154 MeV. The parameter b i corresponds to the width of the chiral crossover. A complete set of those values is given in Table 1. Table 1. Assignment of the negative-parity states in in-medium HRG [5].
LGT [1,4]  We note that the simulations [1,2] were performed for heavier light-quark mass than the physical one, leading to m π = 384 MeV, whereas the strange quark is set to the physical one. To correct the unphysical effect from the heavy up and down quarks, M i − (T = 0) is set to its PDG value in the following calculation. We need to correct the value of M − (T c ) as well and re-scale it by multiplying the factor M PDG . It is uncertain how to assign the negative-parity states to the observed ones because of unknown quantum numbers of some of the candidates. In fact, one finds, e.g. three low-lying Lambda baryons with negative parity on the PDG table: Λ − (1405), Λ − (1670), Λ − (1800). Hence, we follow the suggestion in [4] (Set C) and further adopt two different assignments, Sets A and B. Set A includes the lighter state for Λ − , Ξ − , and Σ * − while Set B the heavier ones.
When the SU(3) limit with a common quark-mass is taken, the chiral model yields a quite similar trend to the lattice observation up to T c . In the realistic setup with the physical m π , one observes a sizable difference between the model and Equation (1), which sets in at a rather low temperature, T/T c ∼ 0.6. The distinct features are understood as remnants of the underlying flavor symmetries and their universality [3]; the SU(3) leads to a first-order chiral phase transition, whereas the SU(2) to a second-order. The physics with N f = 2 + 1 near T c is to a large extent governed by the O(4) universality class, so that such a drastic decrease based on Equation (1) is rather an SU(3) artifact and would not be expected in the realistic QCD thermodynamics.
This becomes more striking in strange baryons. The hyperon masses from the chiral model and from Equation (1) are shown in Figure 2. It is clearly seen that the discrepancy increases with larger strangeness. It would be intriguing to confirm whether such a strong mass-reduction of the hyperons still persists in simulations with a lighter pion.  Figure 2. Temperature dependence of the strange baryons with negative parity for ∆ − , for the Assignment C in Table 1.

Fluctuations and Correlations
In Ref. [4], the observed medium modification was used to possibly explain a missing contribution in correlations among the conserved charges, in the context of the in-medium Hadron Resonance Gas (HRG) model. In the following, we examine the fluctuations and correlations of the net-baryon with net-charge and net-strangeness on the basis of the standard HRG, the in-medium HRG and the chiral approach [5]. In particular, we focus on the quantities χ B 2 , χ BQ 11 , and χ BS 11 at µ B = µ Q = µ S = 0. Since mesons do not contribute to those susceptibilities, they are good measures of in-medium effects in the baryonic sector. We neglect an intrinsic chemical-potential dependence in the baryon masses in the present calculations for simplicity, which is well justified at small chemical potential. By the standard thermodynamic calculations, we obtain the thermal properties of the χ B 2 , χ BQ 11 and −χ BS 11 , as shown in Figure 3.   [8,9] and Budapest-Wuppertal collaborations [10].
The in-medium HRG apparently overshoots individual susceptibilities of lattice QCD because of the strong decrease in the masses of the negative-parity states. We see a moderate enhancement in all the three assignments, Lattices with assignments A, B and C in Table 1, and it becomes stronger when M(T c ) is corrected, since the correction further reduces the masses near T c . The results from the chiral effective model follows the same trend, but the amount of the enhancement is much smaller than the tanh-parameterization (Equation (1)), because of the much weaker mass shifts particularly in the hyperon sectors.
In Figure 4, we present the ratios χ BQ 11 /χ B 2 and −χ BQ 11 /χ BS 11 together with corresponding lattice results.  The results with Lattices (A) and (B) follow the trend seen in the HotQCD data. As becomes clear in Figure 3, however, the individual susceptibilities cannot be reproduced by the in-medium HRG. The coincidence comes from the strong mass reduction in the charge asymmetric states. The chiral approach leads to the opposite trend that the ratio is above the standard HRG, due to the much milder mass-reduction in the hyperon sector. Therefore, we conclude that the lattice data cannot be explained solely by the mass reduction of the negative-parity baryons.

The Structure of Neutron Stars
Parity doubling is expected to play a central role in cold but dense matter as well. We recently applied the chiral effective model with nucleon parity doubling to the physics of neutron stars under β equilibrium and charge neutrality [6]. We utilize a hybrid model where both quarks and hadrons are treated as dynamical degrees of freedom together with a mechanism to suppress the quarks (nucleons) at low (high) density [11,12]. The model yields the chiral symmetry restoration as a first-order, a second-order phase transition or a smooth crossover, depending on two major parameters; one is the degenerate mass of nucleon parity doublers with restored chiral symmetry, m 0 , and the other is the infra-red cutoff of the quark thermal distribution function in matter-free space, αb 0 [11]. The most characteristic feature is that at increasing baryon density, the chiral symmetry is restored within the hadronic phase by lifting the mass splitting between chiral partner states, before quark deconfinement takes place, as depicted in Figure 5 where the energy-density profiles of the stars are shown as a function of the radius. The two parameters m 0 and αb 0 can be restricted by imposing the following astrophysical constraints: 1. The maximum mass of the neutron stars [13]; 2. the onset of the direct Urca process relevant for the cooling of neutron stars [14,15]; and 3. the dimensionless tidal deformability in the binary merger [16].
In Figure 6, we compile those three constraints in the (αb 0 , m 0 )-plane. The region where they overlap gives the most preferable sets of the parameters. The observational neutron-star data provide useful constraints on the structure of strongly interacting matter. They may constrain the phase diagram of isospin-symmetric QCD matter, which is of major relevance for the heavy-ion physics. In Figure 7 In view of the constraints given in Figure 6, the scenarios with αb 0 = 350 MeV and αb 0 = 370 MeV are rather excluded. Hence, either rather low temperature for the critical end point or even its absence in the phase diagram for isospin-symmetric matter is favored.

Conclusions
Motivated by the recent lattice QCD result on the baryon-octet and -decuplet masses at finite temperature, we formulated thermal masses of baryon parity doublers as functions of the chiral order parameter. We examined the consequences of several scenarios with in-medium masses for the negative-parity octet and decuplet baryons, adopting a lattice-motivated parameterization and that from a chiral effective theory. We showed that the reproduction of the lattice QCD results of the ratios χ BQ 11 /χ B 2 and −χ BQ 11 /χ BS 11 is rather accidental when the in-medium masses are naively introduced in the conventional HRG approach. Since the strong mass shifts of the hyperons observed by FASTSUM is a consequence of the approximate flavor SU(3) due to a heavy pion mass, it is critical to confirm if the large hyperon mass-shift persists for the physical pion mass.
We emphasize that the treatment of the resonances in the conventional HRG model is insufficient and a proper inclusion of the widths improves the thermodynamics (see, e.g., [17]). Hence, our results do not rule out the manifestation of parity doubling. To correctly account for the influence of the chiral symmetry restoration over the thermodynamic quantities, a consistent framework beyond the mass shifts is required.
We also presented the consequences of a recently developed model for the equation of state of dense QCD matter under neutron-star conditions and the phenomenology of compact stars. It was shown that the chiral symmetry restoration by parity doubling within the hadronic phase can be archived in the core of massive neutron stars. The high-mass stars, such as the PSR J0348 + 0432 pulsar with the mass 2.01(4)M , can be realized in three different ways: (1) single-phase nuclear matter with broken chiral symmetry; (2) its core is made of a mixed phase surrounded by chirally broken and confined nuclear matter; or (3) two-phase nuclear matter where the core with restored chiral symmetry but still confined is surrounded by chirally broken and confined matter.