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 describes 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 continuum threshold, so that the Levinson theorem is fulfilled 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 its contributions to the thermodynamics vanish. It is demonstrated on this simple example that this unified approach to quark-hadron 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 c (μ = 0) = 156.5±1.5 MeV a e-mail: niels-uwe.bastian@uwr.edu.pl (corresponding author) [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. Perhaps 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 renormalization 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 quark-hadron 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]. By construction, you will get always a first-order phase transition, what is a contradiction to the Lattice results for low chemical potentials. A van-der-Waals model for hadrons allows the construction of a critical endpoint end hence a cross-over at high temperatures [11]. Similar results one obtains with density and temperature dependent excluded volume approaches [12]. These transitions are purely thermodynamic and therefore do not result in a transition to quark matter. An alternative approach is using a switching function to match a hadronic and quark EOS [13,14].
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. An approach based on a Hagedorn bag-like model is described in [15].
The simultaneous description of the thermodynamics of bound states (clusters) and their constituents is the subject of the physics of non-ideal plasmas as well as the particle clustering and Mott transitions in nuclear matter [16,17]. This description has recently been advanced to a cluster virial expansion [18] on the level of the generalized Beth-Uhlenbeck approach for nuclear matter [19]. 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 worked out in Ref. [20] (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 [21,22]. However, at low temperatures and finite baryon densities such a description lacks the knowledge of the Polyakov-loop 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. [23].
In the present work we will employ the concept, which is described in [24], 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. Within the scope of the current study we restrict ourselves to iso-spin symmetric matter.

Cluster expansion
It is shown in [24,25], that the -derivable approach [26,27] can be generalized to resemble a cluster decomposition of the thermodynamic potential where i labels the particle species and stands for quark, meson, diquark and baryon contributions, c i = +1(−1) for bosons (fermions). The full cluster Green's function G i fulfils the Dyson equation where the cluster self energy i is defined as the functional derivative of the -functional i = δ /δG i . In the cluster virial expansion scheme that we follow, we restrict ourselves to the choice of all two-loop diagrams of the "sunset" type that can be drawn with three cluster Green's functions for cluster sizes i, j and k that fulfil the relation i = j + k with i > j, k in order to preserve charge conservation at the vertices. 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 [24]. 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. Here we neglect the meson-meson couplings, which would lead to purely mesonic diagrams. These types of vertices would only play a role in the non-linear extension of the model, which is not considered at this point. 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. 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 [28,29], for both types of diagrams in Fig. 1, a cancellation holds which allows to write the cluster expansion in the form of the generalized + + + + Fig. 1 Contributions of the baryon-meson-quark-diquark system to the -functional. The upper row shows the full set of sunset type diagrams. The lower row shows the resulting diagrams to be taken into account, after collapsing the diquark and meson-propagators. The first and second diagram represent the mean field contributions to the quark and the baryon propagators. Here the contributions result from effective interac-tions 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 an effective coupling term. The last term, which is the remainder of the meson-diquark contributions is a constant of the thermodynamic potential, whose derivatives (like particle density) vanish Beth-Uhlenbeck approach [24,25] 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 i j 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 quasi-particle propagators. Hence the phase shifts defined as δ i = arctan(Im G i /Re G i ), degenerate to functions which can only attain values of nπ , where n is an integer standing for the number of bound states. Using the quasi-particle dispersion relation with the scalar (vector) self energy S i (V i ), we can substitute the energy of the quasi particle by its effective mass M i = m i + S i and obtain 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.
Within this letter, we consider isospin-symmetric matter and therefore can reduce the quark flavours to one degenerate one j = {q} as constituent particles and additionally nucleons (protons and neutrons) as degenerate composites so that i = {q, n}. Hence the degeneracy factors are g i = (12, 4) T and quark constituent of the clusters is A i j = (1, 3) T . The generalization to isospin-asymmetric or strange matter and the inclusion of other baryons (e.g. hyperons), can be done analogously. 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 [30]. It features a bound state (jump up 0 → π ) at the nucleon mass M i and an anti-bound state (jump down π → 0) at the continuum threshold mass M thr i , satisfying in this manner Levinson's theorem. Here the spectral density of the continuum states degenerates to that of a negative delta distribution at the continuum threshold and acts like the negative of the contribution from a bound state in the same channel. As elementary particles, quarks are here described as quasi particles without substructure. The analogue of the threshold in this case is at positive infinity.
The threshold masses of nucleons are defined as the sum of the masses of their constituents M thr n = 3M q . 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. At the quasi-particle mass M i of the nucleon it has a sharp jump from zero to π , referring to the bound state. At the continuum threshold mass M thr i it jumps back from π to zero, resembling a continuum phase shift that degenerates to an anti-bound state Due to the simple shape of the phase shifts, the mass integration in Eq. 5 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 quasi particle density for (unbound) quarks and an altered expression for (bound) nucleons which contains one usual "free" contribution, subtracted by another contribution, based on the mass of the three-quark 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. In chemical equilibrium the chemical potentials μ i of all particles i can be expressed by their conserved charges. For symmetric nuclear matter, where only the baryon chemical potential μ B is relevant, this leads to μ i = B i μ B , where B i is the baryon number of the particle i. Given the relation between densities and chemical potentials, one can obtain the free energy density Consequently, one can easily obtain the pressure as p = μn − f and construct the phase transitions by fulfilling the Gibbs conditions of phase equilibrium.

Self-energies
Recently, we developed a relativistic density functional (RDF) formalism [23], which is capable of handling more complex types of interaction (such as confinement) on the mean field level. Models within this formalism are defined by a density functional U (n s , n v , n vi ), representing the interaction contributions. Here n s , n v , n vi are the scalarisoscalar, vector-isoscalar and vector-isovector densities, which appear for quarks or nucleons as the expectation values of the corresponding bilinear expressions of quark and nucleon fields, respectively. Based on the density functional, one can derive the scalar S and vector V self energies, which are used in the dispersion relation Eq. 4. As the scalar self energy defines the effective masses of particles, it represents the imperative element for the Mott transition.

Quarks
The model for confinement is motivated by the string-flip model (SFM) [31], and already achieved great success in the studies of astrophysical phenomena [23,[32][33][34]. In this work, we restrict ourselves to only the important interactions in order to demonstrate the effect of deconfinement and Mott dissociation with the effective string tension D(n v ) = D 0 exp(−αn 2 q,v ). It captures the aspects of (quark) confinement through its resulting density-dependent scalar self-energy contribution which results in an effective suppression of quark degrees 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.
The values of the vacuum string tension D 0 , the available volume fraction α and the linear vector coupling parameter a are given in the results section.

Hadrons
For description of hadronic degrees of freedom we use the DD2 model [35,36]. 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 impor- All DD2 coupling terms and parameters can be found in [36].

Results
For the scenario with a critical endpoint, the parameters are chosen to be √ D 0 = 290 MeV, α = 0.15fm 6 and a = 12MeV fm 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 cluster-virial 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 upper 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.
In the lower panel we show the threshold mass M thr n and the effective mass of nucleons M n , which depends linearly on the scalar density n s of the system. Once quarks contribute to n s , M n gets lowered, compared to the model without quarks. The reduction of nucleons in the system has a stronger effects to n s than the increase of quarks, which causes an extremum in n s and M n . This effect is not seen in quarks, because of the vector-density-dependent effective string tension, see Eq. 12.
In Fig. 4 we show a phase diagram, based on Maxwell constructions. Besides the liquid-gas transition, the quarkhadron 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 as it is shown in Fig. 3  The appearance of quarks at n B = 0.19 f m −3 is triggering the thermodynamic instability (spinodal). The existence of the third phase is a crucial difference to the twophase 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-all-over scenario, one gets by changing the volume fraction parameter to α = 0.10 fm 6 , which controls how fast deconfinement is reached along density [23]. The corresponding line of the crossover in the phase diagram can be seen in Fig. 4 as well. With these parameters, the twophase approach would still give a first-order phase transition, at slightly higher densities.
Even though one can see in the inset of Fig. 4 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 [38][39][40] due to the appearance of non-vanishing diquark condensate signalling a coloursuperconducting 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 stepup at the effective mass of the baryon describing an on-shell bound state and a step-down located at the mass of the threequark 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 a 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. Furthermore, a generalization to arbitrary charge-fraction (isospin-asymmetry) is straightforward and would make the model applicable for astrophysical systems.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.