A unified quark-nuclear matter equation of state from the cluster virial expansion within the generalized Beth-Uhlenbeck approach

We consider a cluster expansion for strongly correlated quark matter where the clusters are baryons with spectral properties that are described within the generalized Beth-Uhlenbeck approach by a medium dependent phase shift. We employ a simple ansatz for the phase shift which fulfils the Levinson theorem by describing an on-shell bound state with an effective mass and models the continuum by an anti-bound state located at the mass of the three-quark threshold. The quark and baryon interactions are accounted for by the coupling to scalar and vector meson mean fields modelled by density functionals. At increasing density and temperature, due to the different medium-dependence of quark and baryon masses, the Mott dissociation of baryons occurs and the nuclear cluster contributions to the thermodynamics vanish. It is demonstrated on this simple example that this unified approach to quark-nuclear matter is capable of describing crossover as well as first order phase transition behaviour in the phase diagram with a critical endpoint. Changing the meson mean field, the case of a"crossover all over"in the phase diagram is also obtained.


Introduction
The investigation of the phase diagram of quantum chromodynamics (QCD) is one of the major goals of nuclear and particle physics. Monte-Carlo simulations of the QCD partition function on space-time lattices have reached a precision stage where discretization errors and the continuum limit are under control and calculations are performed for physical quark masses providing the pseudocritical temperature of the hadron-to-quark matter crossover transition at vanishing baryo-chemical potential µ as T (µ = 0) = 156.5 ± 1.5 MeV [1]. Unfortunately, these calculations face the sign problem for µ > 0 and cannot elucidate the QCD phase structure in the whole T − µ plane where a central aim is to identify the position of one (or several) critical endpoint(s) beyond that of the nuclear gas-liquid transition. Eventually, there is no such critical endpoint and the hadron-to-quark matter transition is a crossover all over the QCD phase diagram.
In order to make reliable predictions for the QCD phase structure, effective non-perturbative approaches to lowenergy QCD at finite T and µ have been developed, such as the Dyson-Schwinger equation approach [2]. Within this approach remarkable progress has been achieved towards a theory of the QCD phase diagram and a unified equation of state (EoS) of quark-hadron matter [3]. However, the inclusion of baryons into this scheme is still an open task. Another systematic non-perturbative approach to low-energy QCD is based on applying the functional renor-malization group methods. For an overview, see [4]. On the basis of the Polyakov-Quark-Meson model, the phase diagram with a critical endpoint could be obtained, albeit without baryons [5,6]. Recently, this approach has been developed towards a formulation of the hadronization problem, i.e. to describe hadrons as bound states of quarks [7].
While these fundamental approaches are being further developed it is worthwhile to construct effective dynamical models for QCD thermodynamics which are more practical to use and at the same time provide a scheme for the unified description of hadronic and quark-gluon degrees of freedom where hadrons appear as composites (clusters) of quarks and gluons.
In the description of the hadron-quark matter transition it became customary to employ the so-called two-phase approaches where quark and hadronic matter EoS are modelled separately and subsequently joined by a phase transition construction as, e.g., the Maxwell construction. For an early review on such constructions, see [8]. A more recent overview on the thermodynamics of the deconfinement transition with relevant references is given in [9]. This approach has the disadvantage that it ignores the fact that hadrons are bound states of quarks and their thermodynamics should be directly related to the underlying quark dynamics [10]. The goal would be to describe the hadrons as solutions of the equations of motion of multi-quark states in medium so that their dissociation within the Mott effect is obtained naturally. The quark interactions shall capture the most important aspects of confinement, namely, that they allow a sufficient number of hadronic eigenstates up to masses of, e.g., m max ≈ 1.5 GeV, and that no free quark states are excited before the chiral transition is reached.
The simultaneous description of the thermodynamics of bound states (clusters) and their constituents is the subject of the physics of non-ideal plasmas, in particular of particle clustering and Mott transitions in nuclear matter [11,12]. This description has recently been advanced to a cluster virial expansion [13] on the level of the generalized Beth-Uhlenbeck approach for nuclear matter [14]. The extension of such a description to quark matter systems where the elementary degrees of freedom are the quarks and the clusters are the hadrons has been advanced in Ref. [15] (see also references therein). It meets the problem that quark confinement at low densities has to be taken into account. This can effectively be achieved by coupling the chiral quark dynamics to the Polyakov loop [16,17]. However, at low temperatures and finite baryon densities such a description lacks the knowledge of the Polyakovloop potential, because no lattice QCD calculations exist in that domain for fitting it. An alternative has been suggested in the form of a relativistic density functional approach that suppresses quarks at low densities by a diverging scalar self-energy (mass term), see Ref. [18].
In the present work we will employ the concept, which is described in [19], for the implementation of quark confinement/deconfinement and describe the appearance/disappearance of baryons within a simplified realization of the generalized Beth-Uhlenbeck approach that works with generic ansatz for the medium dependence of the baryon phase shifts that encodes their Mott dissociation. As a result we shall obtain the thermodynamics of the non-ideal quark-baryon plasma in the form of a cluster expansion and will be able to relate the existence of a critical endpoint of a first order deconfinement transition to the microphysics of the model.

Cluster Expansion
It is shown in [19,20], that the Φ-derivable approach [21,22] can be generalized to resemble a cluster decomposition of the thermodynamic potential where i labels the clusters and stands for quark, meson, diquark and baryon contributions, c i = +1(−1) for bosons (fermions). For the description of the system within this approach, one would need to include the diagrams, shown in the upper row of Fig. 1, as contributions to the Φfunctional. Note here, that we restrict ourselves to dia-grams of the sunset type. Before solving the entire system, we want to do a first step, by replacing the meson and diquark propagators by their mean fields, as it is also shown in the lower row of diagrams in Fig. 1. The remaining particles, baryons and quarks, are affected by a mean field, coming from the diquark and meson contributions, which we describe effectively by density-dependent self-energies. Diagrammatically, the quark (baryon) selfenergy is obtained by a functional derivative of the Φfunctional w.r.t. the quark (baryon) propagator, i.e. by cutting a quark (baryon) propagator line in the lower row of diagrams in Fig. 1. The third term describes the coupling of quark and baryon sectors and the last term is a remnant of the meson-diquark interaction which in our approximation becomes a constant of the thermodynamic potential and therefore its derivatives, such as particle density, vanish. As it is shown in [23,24], for both types of diagrams in Fig. 1, a cancellation holds which allows to write the cluster expansion in the form of the generalized Beth-Uhlenbeck approach [19,20] with the phase shift δ i (E) as a medium-dependent quantity, which already includes all properties of the cluster i, corresponding to its spectrum of bound and scattering states. The matrix A ij is formed by the number of constituents j in the cluster i, while j ⊂ i. The species i has degeneracy factor g i and obeys the Fermi distribution The chosen diagrams in the lower row of Fig. 1 are of Hartree type. The resulting self energies thus are real and represent shifts of the poles of the quark and nucleon quasiparticle propagators. Hence the phase shifts defined as δ i = arctan Im Gi Re Gi , degenerate to functions which can only attain values of nπ, where n is an integer standing for the number of bound states.
In chemical equilibrium the chemical potentials µ i of all particles i can be expressed by their conserved charges as where B i is the baryon number of the particle i, C i is its electrical charge number and µ B , µ C are the associated chemical potentials. It should be noted, in this work we do not consider the possible occurrence of other charges like strangeness or lepton number. Using the quasiparticle energy-momentum dispersion relation with the scalar (vector) self energy S i (V i ), we can substitute the energy of the quasiparticle by its effective mass + + + + Figure 1: Contributions of the baryon-meson-quark-diquark system to the Φ-functional. Upper row shows the full set of sunset type diagrams. Lower row shows the resulting diagrams to be taken into account, after collapsing the diquark and meson-propagators. The first and second one represent the mean field contributions to the quark and the baryon propagators. Here the contributions result from effective interactions with the meson and diquark fields. The third diagram represents the coupling of baryons to quarks, which is elementary done by diquarks, but here described as effective coupling term. The last term, which is the remaining of the meson-diquark contributions is a constant of the thermodynamic potential, whose derivatives (like particle density) vanish.
Here we replaced the energy-dependent phase shift δ i (E) by a (effective) mass dependent representation δ i (M ), which absorbs the contributions of scalar self energy.
In the system we assume now up-and down-quarks j = {u, d} as constituent particles and additionally protons and neutrons as composites so that i = {u, d, p, n}. Hence the degeneracy factors are g i = (6, 6, 2, 2) and quark constituent of the clusters is We define the phase shifts as for nucleons and quarks, respectively. The simple ansatz for the phase shift as depicted in Fig. 2 has already been utilized before in the description of composite pions in quark matter in [25]. It features a bound state (jump up from 0 to π) at the respective nucleon mass as well as an anti-bound state (jump from π to zero) at a threshold mass M thr i as a degenerate three-quark continuum phase shift satisfying in this manner Levinson's theorem. As elementary particles, quarks are here described as bound states without substructure, i.e. with infinite binding energy.
The threshold masses of nucleons are defined as the sum In this way, the bound state is density-dependent via the self-energy shifts of the quarks as well as the nucleons and, once the constituent masses drop below the nucleon mass, the Mott transition occurs. Due to the simple shape of the phase shifts, the mass integration in Eq. (6) can be performed analytically. Note here, that the derivative of the Heavyside distributions Θ(x) results in Dirac distributions and sin(π/2) = 1 at the steps from zero to π. As result we get the usual quasiparticle density for (free) quarks and an altered expression for (bound) nucleons which contains one usual "free" contribution, subtracted by another contribution, based on the mass of the threequark continuum threshold. One can see that once the threshold mass drops below the nucleon mass at the Mott transition, the entire nucleon particle density vanishes.

Self-energies
Recently, we developed a relativistic density functional (RDF) formalism [18], which is capable of handling more complex types of interaction, such as confinement.

Quarks
The model for confinement is motivated by the stringflip model (SFM) [26], and already achieved great success in the studies of astrophysical phenomena [18,[27][28][29]. In this work, we restrict ourself to only the important interactions in order to demonstrate the effect of deconfinement and Mott dissociation It captures the aspects of (quark) confinement through its resulting density-dependent scalar self-energy contribution with the effective string tension which results in an effective suppression of quark degree of freedom (dof) at low densities and temperatures. Furthermore, we include a minimal description of vector repulsion, which is for the discussion of this paper not relevant, but necessary to move the transition to reasonable densities. Note that the formulation of SFM is based on quark quantities. To have a consistent description in this paper, we explicitly added an index q for the quark densities and introduced the scalar quark density n q,s = 3n s as well as the vector quark density n q,v = 3n v , what is important in order to compare the parameters with the original work.

Hadrons
For description of hadronic degrees of freedom we use the DD2 model [30,31]. It is a nucleonic model with mean field meson interactions of the Walecka type. The model has proven to be very accurate in describing properties of nuclei and nuclear matter. In the present work it is only important, that it describes hadronic degrees of freedom and itself already features the nucleonic liquid-gas transition. The DD2 model can be expressed in the RDF formalism with the potential which results in the scalar self-energy and the species dependent vector self-energy All DD2 coupling terms and parameters can be found in [31].

Results
For the scenario with a critical endpoint, the parameters are chosen to be √ D 0 = 290 MeV, α = 0.15 fm 6 and a = 12 MeVfm 3 . This choice is not due to any constraints or physical application. It just states the most simple example to demonstrate the feature of the cluster expansion in this scenario.
The resulting EoS can be seen in Fig. 3, where we compare the common two-phase approach with the clustervirial expansion, utilizing the Maxwell construction in both cases. Once there is an occupation of quarks in the system, the EoS gets altered from the pure hadronic one. In our case this leads directly into thermodynamic instability ∂µ/∂n < 0, most likely due to choice of parameters. This instability can be seen as a second Van-der-Waals wiggle and hence indicates a first-order phase transition as it can be seen in the right panel of Fig. 3. Now we have two first-order phase transitions, one for the pure hadronic liquid-gas transition and one for the quark-hadron transition, connected to the deconfinement of quarks. Due to the extended persistence of hadrons as bound states in the quark-gluon plasma, this phase gets softer and the phase transition is at much lower densities than in the two-phase approach. Fig. 4 illustrates the phase transition from the hadronic to the quark EoS in the pressure versus chemical potential plane, obtained by Maxwell constructions for both, the cluster expansion and the two-phase approach. The inset magnifies the region of phase transitions of the cluster expansion. The pressure was obtained via the free energy density by numerical integration.
In Fig. 5 we show a phase diagram, based on Maxwell constructions. One can see, that besides the liquid-gas transition, the quark-hadron phase transition has a critical endpoint, which separates the region of first-order transitions from that of smooth crossover transitions. The crossover was obtained from the inflection point in the µ − n plot. This calculation features a very tiny region of hadronic liquid and the critical temperature of the quark-hadron transition is comparably small. Once again we want to   mention, that the parameter choice is purely academic and a thorough parameter scan needs to be performed. An interesting outcome of our model is the fact that the thermodynamic phase transition does not coincide with the microscopic Mott transition (nucleon dissociation). Besides the liquid-gas transition (n B < 0.15 fm −3 ), we obtain 4 different regions along the density axis (in brackets we indicate the borders at T = 0): The appearance of quarks at n B = 0.19 fm −3 is triggering the thermodynamic instability (spinodal). The existence of the third phase is a crucial difference to the two-phase approach. This phase is on the high-density end of the Maxwell construction. Due to the additional (hadronic) degrees of freedom it is much softer (compared to the pure quark phase), resulting in the exceptionally early onset. Another parameter set, featuring the possibility of a crossover-allover scenario, one gets by changing the vol-ume fraction parameter to α = 0.10 fm 6 , which controls how fast deconfinement is reached along density [18]. The corresponding line of the crossover in the phase diagram can be seen in Fig. 5 as well. With these parameters, the two-phase approach would still give a first-order phase transition, at slightly higher densities.
Even though one can see in the right panel of Fig. 5 a slight increase of the spinodal region at moderate temperatures, there was no scenario found with two critical endpoints. Such scenario is suggested as a possibility of the quark-hadron continuity [33][34][35] due to the appearance of nonvanishing diquark condensate signalling a colorsuperconducting phase of quark matter. The inclusion of the diquark correlations is straightforward in the present approach but goes beyond the scope of this exploratory work.

Conclusions
We have demonstrated that a unified quark-nuclear EoS can be obtained within a cluster expansion for strongly correlated quark matter where the clusters are baryons with spectral properties that are described within the generalized Beth-Uhlenbeck approach by a medium-dependent phase shift. As an instructive example we have discussed here the simple and generic model for the phase shift of baryons with a step-up at the effective mass of the baryon describing an on-shell bound state and a step-down located at the mass of the three-quark continuum threshold which models the continuum as an anti-bound state. This simple ansatz fulfils the Levinson theorem by construction. The quark and baryon interactions are accounted for by the coupling to scalar and vector meson mean fields modelled by density functionals. At increasing density and temperature, due to the different medium dependence of quark and baryon masses, the Mott dissociation of baryons occurs and the contributions of nuclear clusters to the thermodynamics vanish. It is demonstrated on this simple example that this unified approach to quark-nuclear matter is capable of describing a first order phase transition with a critical endpoint as well as the case of crossover all over.
The next step would be to include a larger set of particles like diquarks, mesons and baryons and to include a more sophisticated functional for the quark interaction. With a sufficient amount of particles it should be possible to adjust the parameter set in such a way that it follows Lattice QCD calculations at zero density and fulfils known constraints for low temperatures.