**1. Introduction**

The phase diagram of matter is one of the most fascinating subjects in physics, which also has important implications on chemistry and biology. Its phase boundaries and (possibly) critical points have been the focus of physics research for centuries. Apart from the traditional phase diagram in the plane of temperature *T* and pressure *P*, its transport properties like the shear and bulk viscosities, the electric conductivity, etc. are also of fundamental interest. These transport coefficients emerge from the stationary limit of correlators and provide additional information on the systems in thermal and chemical equilibrium apart from the equation of state. In this context, the phase diagram of strongly interacting matter has been the topic of most interest for decades and substantial experimental and theoretical efforts have been invested to shed light on this issue. It contains the information about the properties of our universe from early beginning—directly after the Big Bang—when the matter was in the QGP phase at very large temperature *T* and about zero baryon chemical potential *μ<sup>B</sup>*, to the later stages of the universe, where in the expansion phase stars and Galaxy have been formed. Here, the matter is at low temperature and large baryon chemical potential. Relativistic and

#### *Particles* **2020**, *3*

ultra-relativistic heavy-ion collisions (HICs) nowadays offer the unique possibility to study some of these phases, in particular a QGP phase and its phase boundary to the hadronic one. Furthermore, the phase diagram of strongly interacting matter in the (*<sup>T</sup>*, *μB*) plane can also be explored in the astrophysical context at moderate temperatures and high *μB* [1], i.e., in the dynamics of supernovae or—more recently—in the dynamics of neutron star merges.

In order to reproduce the mini Big Bangs in laboratories, heavy-ion accelerators are built which allow for investigating the creation of the QGP under controlled conditions. Hadronic spectra and relative hadron abundances from these experiments reflect important aspects of the dynamics in the hot and dense zone formed in the early phase of the reaction and collective flows provide information on the transport properties of the medium generated on short time scales. Whereas heavy-ion reactions at Relativistic Heavy-Ion Collider (RHIC) and Large Hadron Collider (LHC) energies probe a partonic medium at small baryon chemical potential *μ<sup>B</sup>*, the current interest is in collisions at lower bombarding energies where the net baryon density is higher and *μB* accordingly. Such conditions will be realized in future accelerators at the Facility for Antiproton and Ion Research (FAIR) in Darmstadt and the Nuclotron-based Ion Collider fAcility (NICA) in Dubna.

Current methods to explore QCD in Minkowski space for non-vanishing quark (or baryon) densities (or chemical potential) are effective approaches. Using effective models, one can study the dominant properties of QCD in equilibrium, i.e., thermodynamic quantities as well as transport coefficients. To this aim, the dynamical quasiparticle model (DQPM) has been introduced [2–6], which is based on partonic propagators with sizeable imaginary parts of the self-energies incorporated. Whereas the real part of the self-energies can be attributed to a dynamically generated mass (squared), the imaginary parts contain the information about the interaction rates in the system. Furthermore, the imaginary parts of the propagators define the spectral functions of the degrees of freedom which might show narrow (or broad) quasiparticle peaks. A further advantage of a propagator based approach is that one can formulate a consistent thermodynamics [7] as well as a causal theory for non-equilibrium configurations on the basis of Kadanoff–Baym equations [8].

In order to explore the properties of the QGP close to equilibrium, transport coefficients are calculated such as shear *η* and bulk *ζ* viscosities, electric conductivity *σ*0, etc. While basically all of the effective models have similar equations of state (EoS), which match well with available lattice QCD data at *μB* = 0, the transport coefficients can vary significantly for different models (cf. [9]). Exploration of transport coefficients of the hot and dense QGP can provide useful information for simulations of heavy-ion collisions (HIC) based on hydrodynamical models for which they are used as input parameters. The experimental data for elliptic flow can be well reproduced by hydrodynamical simulations with a small value for the shear viscosity over entropy density [10,11].

Since relativistic heavy-ion collisions start with impinging nuclei in their groundstates, a proper non-equilibrium description of the entire dynamics through possibly different phases up to the final asymptotic hadronic states—eventually showing some degree of equilibration—is mandatory. To this aim, the Parton–Hadron–String Dynamics (PHSD) transport approach [5,12–15] has been formulated more then a decade ago (on the basis of the Hadron-String-Dynamics (HSD) approach [16]), and it was found to well describe observables from p+A and A+A collisions from SPS to LHC energies including electromagnetic probes such as photons and dileptons [5]. In order to explore the partonic systems at higher *μ<sup>B</sup>*, the PHSD approach has been recently extended to incorporate partonic quasiparticles and their differential cross sections that depend not only on temperature *T* as in the previous PHSD studies, but also on chemical potential *μB* explicitly—cf. [17]. Within this extended approach, we have studied the 'bulk' observables in HIC for different energies—from AGS to RHIC, and systems—strongly asymmetric C+Au and symmetric Au+Au/Pb+Pb collisions. We have found only a small influence of *μB* dependences of parton properties (masses and widths) and their interaction cross sections in bulk observables [17].

In this work, we extend our study from Ref. [17] to the collective flow (*<sup>v</sup>*1, *v*2) coefficients and their sensitivity to the *μB* dependences of partonic cross sections. In addition, we explore the relations between the in and out-of equilibrium QGP by means of transport coefficients and collective flows. Additionally, we show explicitly the 'bulk' results for asymmetric heavy-ion collisions such as Cu+Au and discuss which hadronic species and observables are more sensitive to such effects.

#### **2. The PHSD Approach**

We start with reminding the basic ideas of the PHSD transport approach and the DQPM. The Parton–Hadron–String Dynamics transport approach [5,12–15] is a microscopic off-shell transport approach for the description of strongly interacting hadronic and partonic matter in and out-of equilibrium. It is based on the solution of Kadanoff–Baym equations in first-order gradient expansion [13] employing 'resummed' propagators from the dynamical quasiparticle model (DQPM) [2–4] for the partonic phase.

The Dynamical Quasiparticle Model (DQPM) has been introduced in Refs. [2–4] for the effective description of the properties of the QGP in terms of strongly interacting quarks and gluons with properties and interactions which are adjusted to reproduce lQCD results on the thermodynamics of the equilibrated QGP at finite temperature *T* and baryon (or quark) chemical potential *μq*. In the DQPM, the quasiparticles are characterized by single-particle Green's functions (in propagator representation) with complex self-energies. The real part of the self-energies is related to the mean-field properties, whereas the imaginary part provides information about the lifetime and/or reaction rates of the particles.

In PHSD, the partons (quarks and gluons) are strongly interacting quasiparticles characterized by broad spectral functions *ρj* (*j* = *q*, *q*¯, *g*), i.e., they are off-shell contrary to the conventional cascade or transport models dealing with on-shell particles, i.e., the *δ*-functions in the invariant mass squared. The quasiparticle spectral functions have a Lorentzian form [5] and depend on the parton mass and width parameters:

$$\rho\_{\dot{j}}(\omega, \mathbf{p}) = \frac{\gamma\_{\dot{j}}}{E\_{\dot{j}}} \left( \frac{1}{(\omega - E\_{\dot{j}})^2 + \gamma\_{\dot{j}}^2} - \frac{1}{(\omega + E\_{\dot{j}})^2 + \gamma\_{\dot{j}}^2} \right) \tag{1}$$

separately for quarks/antiquarks and gluons (*j* = *q*, *q*¯, *g*). With the convention *<sup>E</sup>*<sup>2</sup>(**p**<sup>2</sup>) = **p**2 + *M*<sup>2</sup> *j* − *γ*2 *j* , the parameters *M*<sup>2</sup> *j* and *γj* are directly related to the real and imaginary parts of the retarded self-energy, e.g., Π*j* = *M*<sup>2</sup> *j*− <sup>2</sup>*iγj<sup>ω</sup>*.

The actual parameters in Equation (1), i.e., the gluon mass *Mg* and width *<sup>γ</sup>g*—employed as input in the present PHSD calculations—as well as the quark mass *Mq* and width *γq*, are depicted in Figure 1 as a function of the temperature *T* and baryon chemical potential *μ<sup>B</sup>*. These values for the masses and widths have been fixed by fitting the lattice QCD results from Refs. [18,19] in thermodynamic equilibrium. One can see that the masses of quarks and gluons decrease with increasing *μ<sup>B</sup>*, and a similar trend holds for the widths of partons.

A scalar mean-field *Us*(*ρs*) for quarks and antiquarks can be defined by the derivative of the potential energy density with respect to the scalar density *ρ<sup>s</sup>*,

$$\mathcal{U}\_s(\rho\_s) = \frac{dV\_p(\rho\_s)}{d\rho\_s},\tag{2}$$

which is evaluated numerically within the DQPM. Here, the potential energy density is defined by

$$V\_p(T, \mu\_q) = T\_{\mathcal{K}^-}^{00}(T, \mu\_q) + T\_{q-}^{00}(T, \mu\_q) + T\_{q-}^{00}(T, \mu\_q), \tag{3}$$

where the different contributions *T*<sup>00</sup> *j*− correspond to the space-like part of the energy-momentum tensor component *T*<sup>00</sup> *j* of parton *j* = *g*, *q*, *q*¯ (cf. Section 3 in Ref. [3]). The scalar mean-field *Us*(*ρs*) for quarks and antiquarks is repulsive as a function of the parton scalar density *ρs* and shows that the scalar mean-field potential is in the order of a few GeV for *ρs* > 10 fm−3. The mean-field potential (2) is

employed in the PHSD transport calculations and determines the force on a partonic quasiparticle *j*, i.e., ∼ *Mj*/*Ej*<sup>∇</sup>*Us*(*x*) = *Mj*/*Ej dUs*/*dρs* ∇*ρs*(*x*), where the scalar density *ρs*(*x*) is determined numerically on a space-time grid.

**Figure 1.** The effective quark (**left**) and gluon (**right**) masses *M* (**upper row**) and widths *γ* (**lower row**) as a function of the temperature *T* and baryon chemical potential *μB*.

Furthermore, a two-body interaction strength can be extracted from the DQPM as well from the quasiparticle width in line with Ref. [2]. On the partonic side, the following elastic and inelastic interactions are included in the latest version of PHSD (v. 5.0) *qq* ↔ *qq*, *q*¯*q*¯ ↔ *q*¯*q*¯, *gg* ↔ *gg*, *gg* ↔ *g*, *qq* ¯ ↔ *g*, *qg* ↔ *qg*, *gq* ¯ ↔ *gq* ¯ exploiting 'detailed-balance' with cross sections calculated from the leading order Feynman diagrams employing the effective propagators and couplings *g*<sup>2</sup>(*T*/*Tc*) from the DQPM [17]. In Ref. [17], the differential and total off-shell cross sections have been evaluated as a function of the invariant energy of colliding off-shell partons √*s* for each *T*, *μ<sup>B</sup>*. We remind that in the previous PHSD studies (using v. 4.0 and below) the cross sections depend only on *T* as evaluated in Ref. [20].

When implementing the differential cross sections and parton masses into the PHSD5.0 approach, one has to specify the 'Lagrange parameters' *T* and *μB* in each computational cell in space-time. This has been done by employing the lattice equation of state and a diagonalization of the energy-momentum tensor from PHSD as described in Ref. [17].

The transition from partonic to hadronic d.o.f. (and vice versa) is described by covariant transition rates for the fusion of quark–antiquark pairs or three quarks (antiquarks), respectively, obeying flavor current–conservation, color neutrality as well as energy–momentum conservation [14]. Since the dynamical quarks and antiquarks become very massive close to the phase transition, the formed resonant 'prehadronic' color-dipole states (*qq*¯ or *qqq*) are of high invariant mass, too, and sequentially decay to the groundstate meson and baryon octets, thus increasing the total entropy.

On the hadronic side, PHSD includes explicitly the baryon octet and decouplet, the 0<sup>−</sup>- and 1<sup>−</sup>-meson nonets as well as selected higher resonances as in the Hadron–String–Dynamics (HSD) approach [16]. Note that PHSD and HSD (without explicit partonic degrees-of-freedom) merge at low energy density, in particular below the local critical energy density *εc* ≈ 0.5 GeV/fm3.

## **3. Transport Coefficients**

The transport properties of the QGP close to equilibrium can be characterized by various transport coefficients. The shear viscosity *η* and bulk viscosity *ζ* describe the fluid's dissipative corrections at leading order. Both coefficients are generally expected to depend on the temperature *T* and baryon chemical potential *μ<sup>B</sup>*. In the hydrodynamic equations, the viscosities appear as dimensionless ratios, *η*/*s* and *ζ*/*<sup>s</sup>*, where *s* is the fluid entropy density. Such specific viscosities are more meaningful than the unscaled *η* and *ζ* values because they describe the magnitude of stresses inside the medium relative to its natural scale.

In our recent studies [17,21,22], we have investigated the transport properties of the QGP in the (*<sup>T</sup>*, *μB*) plane based on the DQPM. One way to evaluate the viscosity coefficients of partonic matter is the Kubo formalism [23–26], which was used to calculate the viscosities for a previous version of the DQPM within the PHSD transport approach in a box with periodic boundary conditions in Ref. [27] as well as in the latest study with the DQPM model in Refs. [17,21]. Another way to calculate transport coefficients (explored also in [17,21]) is to use the relaxation–time approximation (RTA) [28–31].

The shear viscosity based on the RTA (cf. [32]) reads as:

$$\eta^{\rm RTA}(T,\mu\_{\rm q}) = \frac{1}{15T} \sum\_{i=q, p\_{\rm q\xi}} \int \frac{d^3 p}{(2\pi)^3} \frac{\mathbf{p}^4}{E\_i^2} \mathbf{r}\_i(\mathbf{p}, T\_\nu \mu) \, d\_i (1 \pm f\_i) f\_i. \tag{4}$$

where *dq* = 2*Nc* = 6 and *dg* = <sup>2</sup>(*N*2*c* − 1) = 16 are degeneracy factors for spin and color in case of quarks and gluons, *τi* are relaxation times. Equation (4) includes the Bose enhancement and Pauli-blocking factors, respectively. The pole energy is *E*2*i* = *p*2 + *<sup>M</sup>*2*i* , where *Mi* is the pole mass given in the DQPM. The notation ∑*j*=*q*,*q*¯,*<sup>g</sup>* includes the contribution from all possible partons which in our case are the gluons and the (anti-)quarks of three different flavors (*<sup>u</sup>*, *d*,*<sup>s</sup>*).

We consider two cases for the relaxation time for quarks and gluons: (1) *<sup>τ</sup>i*(**<sup>p</sup>**, *T*, *μ*) = 1/Γ*i*(**<sup>p</sup>**, *T*, *μ*) and (2) *<sup>τ</sup>i*(*<sup>T</sup>*, *μ*) = 1/2*γi*(*<sup>T</sup>*, *μ*), where <sup>Γ</sup>*i*(**<sup>p</sup>**, *T*, *μ*) is the parton interaction rate, calculated microscopically from the collision integral using the differential cross sections for parton scattering, while *<sup>γ</sup>i*(*<sup>T</sup>*, *μ*) is the width parameter in the parton propagator (1).

In the left part (a) of Figure 2, we show the ratio of the shear viscosity to entropy density as a function of the scaled temperature *T*/*Tc* for *μB* = 0 calculated using the Kubo formalism (green solid line) and RTA approach with the interaction rate Γon (red solid line) and the DQPM width 2*γ* (dashed green line). The RTA approximation (4) of the shear viscosity with the DQPM width 2*γ* and with the interaction rate Γon are quite close to each other at *μB* = 0 and also very close to the one from the Kubo formalism [17] indicating that the quasiparticle limit (*γ M*) holds in the DQPM.

**Figure 2.** Plot (**a**): the ratio of shear viscosity to entropy density as a function of the scaled temperature *T*/*Tc* for *μB* = 0 calculated using the Kubo formalism (green solid line) and the RTA approach with the interaction rate Γon (red solid line) and the DQPM width 2*γ* (dashed green line). The dashed gray line demonstrates the Kovtun–Son–Starinets bound [33] (*η*/*s*)KSS = 1/(<sup>4</sup>*π*), and the symbols show lQCD data for pure SU(3) gauge theory taken from Ref. [34] (pentagons). The solid blue line shows the results from a Bayesian analysis of experimental heavy-ion data from Ref. [35]. Plot (**b**): the ratio of the bulk viscosity to entropy density *ζ*/*s* as a function of the scaled temperature *T*/*Tc* for *μB* = 0 calculated using the RTA approach with the on-shell interaction rate Γon (red solid line) and the DQPM width 2*γ* (dashed green line). The symbols correspond to the lQCD data for pure SU(3) gauge theory taken from Refs. [36] (pentagons) and [37] (circles). The solid blue line shows the results from a Bayesian analysis of experimental heavy-ion data from Ref. [35]. The dot-dashed and dashed lines correspond to the results from the non-conformal holographic model for *φM* = 2 and 3, correspondingly, from Ref. [38].

The ratio *η*/*s* increases with an increase of the scaled temperature. The actual values for the ratio *η*/*s* are in a good agreemen<sup>t</sup> with the gluodynamic lattice QCD calculations at *μB* = 0 from Ref. [34]. Moreover, our DQPM results are in qualitative agreemen<sup>t</sup> with the results from a Bayesian analysis of experimental heavy-ion data from Ref. [35]. We mention that the DQPM result differs from the recent calculations for the shear viscosity at *μB* = 0 in the quasiparticle model in Ref. [39] where the width of quasiparticles is not considered which leads to a high value for the *η*/*s* ratio. This shows the sensitivity of this ratio to the modelling of partonic interactions and the properties of partons in the hot QGP medium. We remind also that in Refs. [17,21,22] we find that the ratio *η*/*s* shows a very weak dependence on *μB*and has a similar behavior as a function of temperature for all *μB*≤ 400 MeV.

The expression for the bulk viscosity of the partonic phase derived within the RTA reads (following Ref. [30])

$$\zeta^{\rm REA}(T,\mu) = \frac{1}{9T} \sum\_{i=q,\theta,\emptyset} \int \frac{d^3 p}{(2\pi)^3} \text{tr}(\mathbf{p}\_i T, \mu) \left. \frac{d\_i (1 \pm f\_i) f\_i}{E\_i^2} \left( \mathbf{p}^2 - 3c\_s^2 \left( E\_i^2 - T^2 \frac{dm\_i^2}{dT^2} \right) \right)^2 \right|\_{\mathbf{p}=\mathbf{p}\_i, \mathbf{p}\_i = \mathbf{p}\_i} \tag{5}$$

where *c*2*s* is the speed of sound squared, and *dm*2*i dT*<sup>2</sup> is the DQPM parton mass derivative which becomes large close to the critical temperature *Tc*.

On the right side (b) of Figure 2, we show the ratio of the bulk viscosity to entropy density *ζ*/*s* as a function of the scaled temperature *T*/*Tc* for *μB* = 0 calculated using the RTA approach with the interaction rate Γon (red solid line) and the DQPM width 2*γ* (dashed green line). The symbols correspond to the lQCD data for pure SU(3) gauge theory taken from Refs. [36] (pentagons) and [37] (circles). The solid blue line shows the results from a Bayesian analysis of experimental heavy-ion data from Ref. [35]. The dot-dashed and dashed lines correspond to the results from the non-conformal holographic model [38] for *φM* = 2 and 3, correspondingly, where *φM* is the model parameter which characterizes the non-conformal features of the model. We find that the DQPM result for *ζ*/*s* is in very good agreemen<sup>t</sup> with the lattice QCD results and shows a rise closer to *TC* contrary to the holographic results, which show practically a constant behavior independent of model parameters. This rise is attributed to the increase of the partonic mass closer to *TC* as shown in Figure 1, thus the mass derivative term in Equation (5) also grows. The Bayesian result also shows a peak near *TC*; however, the ratio drops to zero while lQCD data indicate the positive *ζ*/*s* as found also in the DQPM. The *μB* dependence of *ζ*/*s* has been investigated within the DQPM in Refs. [17,21,22], where it has been shown that it is rather weak for *μB* ≤ 400 MeV, similar to *η*/*<sup>s</sup>*. As follows from hydrodynamical calculations, the results for the flow harmonic *vn* is sensitive to the transport coefficients [10,11,35]. Thus, there are hopes to observe a *μB* sensitivity of *v*1, *v*2.

## **4. Heavy-Ion Collisions**

In our recent study [17], we have investigated the sensitivity of 'bulk' observables such as rapidity and transverse momentum distributions of different hadrons produced in heavy-ion collisions from AGS to top RHIC energies on the details of the QGP interactions and the properties of partonic degrees-of-freedom. For that, we have considered the following three cases:

(1) **'PHSD4.0'**: the masses and widths of quarks and gluons depend only on *T*. The cross sections for partonic interactions depend only on *T* as evaluated in the 'box' calculations in Ref. [20] in order to merge the QGP interaction rates from all possible partonic channels to the total temperature dependent widths of the DQPM propagator. This has been used in the PHSD code (v. 4.0 or below) for extended studies of many hadronic observables in p+A and A+A collisions at different energies [5,12–15,40].

(2) **'PHSD5.0 -** *μB* = 0**'**: the masses and widths of quarks and gluons depend only on *T*; however, the differential and total partonic cross sections are obtained by calculations of the leading order Feynman diagrams employing the effective propagators and couplings *g*<sup>2</sup>(*T*/*Tc*) from the DQPM at *μB* = 0 [17]. Thus, the cross sections depend explicitly on the invariant energy of the colliding partons √*s* and on *T*. This is realized in the PHSD5.0 by setting *μB* = 0, cf. [17].

(3) **'PHSD5.0 -** *μB***'**: the masses and widths of quarks and gluons depend on *T* and *μB* explicitly; the differential and total partonic cross sections are obtained by calculations of the leading order Feynman diagrams from the DQPM and explicitly depend on invariant energy √*<sup>s</sup>*, temperature *T* and baryon chemical potential *μ<sup>B</sup>*. This is realized in the full version of PHSD5.0, cf. [17].

The comparison of the 'bulk' observables for A+A collisions within the three cases of PHSD in Ref. [17] has illuminated that they show a very low sensitivity to the *μB* dependences of parton properties (masses and widths) and their interaction cross sections such that the results from PHSD5.0 with and without *μB* were very close to each other. Only in the case of kaons, antiprotons *p*¯ and antihyperons Λ¯ + Σ¯ 0, a small difference between PHSD4.0 and PHSD5.0 could be seen at top SPS and top RHIC energies. A similar trend has been found for very asymmetric collisions of C+Au: a small sensitivity to the partonic scatterings was found in the kaon and antibaryon rapidity distributions too. This could be understood as following: at high energies such as top RHIC where the QGP volume is very large in central collisions, the *μB* is very low, while, when decreasing the bombarding energy—in order to increase *μ<sup>B</sup>*, the fraction of the QGP is decreasing such that the final observables are dominated by the hadronic phase, i.e., the probability for the hadrons created at the QGP hadronization to re-scatter, decay, or be absorbed in hadronic matter increases strongly; accordingly, the sensitivity to the properties of the QGP is washed out.

## *4.1. Asymmetric Systems*

In Ref. [17], we have investigated the sensitivity to *μB* of the 'bulk' observables in asymmetric heavy-ion collisions for C+Au. The spectra for C+Au indicated that they show a slightly larger sensitivity to *μB* for antiprotons and strange hadrons - kaons and antihyperons than for pions and protons. Here, we present the results for the asymmetric Cu+Au collisions.

In Figures 3 and 4, we show the rapidity distributions (left plot) and *pT*-spectra at midrapidity (|*y*| < 0.5) (right plot) for *<sup>π</sup>*±, *<sup>K</sup>*±, *p*, *p*¯, Λ + Σ0, Λ¯ + Σ¯ 0 for 10% central Cu+Au collisions at 30 AGeV and √*sNN* = 200 GeV for three cases: (1) PHSD4.0 (green dot-dashed lines), (2) PHSD5.0 with partonic cross sections and parton masses/widths calculated for *μB* = 0 (blue dashed lines) and (3) with cross sections and parton masses/widths evaluated at the actual chemical potential *μB* in each individual space-time cell (red lines). Similar to C+Au collisions, we find for Cu+Au collisions a small difference in the rapidity distributions of antiprotons and in the strangeness sector - in kaon and especially in Λ ¯ + Σ ¯ 0 *y*-distributions. Similar statements hold for the *pT* spectra which show a slightly different slope at low and high momenta of anti-strange baryons. This suggests that the strange degree-of-freedom might be experimentally explored in asymmetric systems to obtain additional information on the partonic interactions.

**Figure 3.** The rapidity distributions (left plot) and *pT*-spectra at midrapidity (|*y*| < 0.5) (right plot) for *<sup>π</sup>*±, *<sup>K</sup>*±, *p*, *p*¯, Λ + Σ0, Λ¯ + Σ¯ 0 for 10% central Cu+Au collisions at 30 A GeV for PHSD4.0 (green dot-dashed lines), PHSD5.0 with partonic cross sections and parton masses calculated for *μB* = 0 (blue dashed lines) and with cross sections and parton masses evaluated at the actual chemical potential *μB* in each individual space-time cell (red lines).

**Figure 4.** The rapidity distributions (left plot) and *pT*-spectra at midrapidity (|*y*| < 0.5) (right plot) for *<sup>π</sup>*±, *<sup>K</sup>*±, *p*, *p*¯, Λ + Σ0, Λ¯ + Σ¯ 0 for 10% central Cu+Au collisions at √*sNN* = 200 GeV for PHSD4.0 (green dot-dashed lines), PHSD5.0 with partonic cross sections and parton masses calculated for *μB* = 0 (blue dashed lines) and with cross sections and parton masses evaluated at the actual chemical potential *μB* in each individual space-time cell (red lines).

## *4.2. Directed Flow*

Now, we test the traces of *μB* dependences of the QGP interaction cross sections in collective observables such as directed flow *v*1 considering again three cases of the PHSD as discussed above.

Figure 5 depicts the directed flow *v*1 of identified hadrons (*K*±, *p*, *p*¯, Λ + Σ0, Λ¯ + Σ¯ 0) versus rapidity for √*sNN* = 27 GeV. One can see a good agreemen<sup>t</sup> between PHSD results and experimental data from STAR collaboration [41]. However, the different versions of PHSD for the *v*1 coefficient show a quite similar behavior; only antihyperons show a slightly different flow. This supports again the finding that strangeness, and in particular anti-strange hyperons, are the most sensitive probes for the QGP properties.

**Figure 5.** Directed flow of identified hadrons as a function of rapidity at √*sNN* = 27 GeV for PHSD4.0 (green lines), PHSD5.0 with partonic cross sections and parton masses calculated for *μB* = 0 (blue dashed lines) and with cross sections and parton masses evaluated at the actual chemical potential *μB* in each individual space-time cell (red lines) in comparison to the experimental data of the STAR Collaboration [41].

## *4.3. Elliptic Flow*

As follows from the hydrodynamic simulations [10,11] and from the Bayesian analysis [35], an elliptic flow *v*2 is sensitive to the transport properties of the QGP characterized by transport coefficients such as shear *η* and bulk *ζ* viscosities. In this section, we present the results for the elliptic flow of charged hadrons from HIC within the PHSD5.0 with and without *μB* dependences and compare the results with PHSD4.0, again.

The left plots '(a)' in Figures 6 and 7 display the actual results for charged hadron elliptic flow as a function of pseudo-rapidity *η* (Figure 6) and of transverse momentum *pT* (Figure 7) for minimum bias Au+Au collisions at √*sNN* = 200 GeV for PHSD4.0 (green lines), PHSD5.0 with partonic cross sections and parton masses calculated for *μB* = 0 (blue dashed lines), and with cross sections and parton masses evaluated at the actual chemical potential *μB* in each individual space-time cell (red lines) in comparison to the experimental data from the STAR collaboration [42] (solid stars) and PHOBOS [43] (solid dots). One can see the difference for *<sup>v</sup>*2(*pT*) in case of charged hadrons for high *pT* > 0.5 GeV between PHSD4.0 and PHSD5.0.

*Particles* **2020**, *3*

**Figure 6.** Plot (**a**): elliptic flow of charged hadrons as function of pseudo-rapidity *η* for minimum bias Au+Au collisions at √*sNN*=200 GeV for PHSD4.0 (green lines), PHSD5.0 with partonic cross sections and parton masses calculated for *μB* = 0 (blue dashed lines), and with the actual *μB* (red lines) in comparison to the experimental data from STAR [42] (solid starts) and PHOBOS [43] (solid dots). Plot (**b**): individual contributions to *v*2 without their relative weights to the total *v*2, which are indicated by a green solid line for PHSD5.0 with *μB*: the red dotted line corresponds to the final hadrons coming from the QGP without rescattering in the hadronic phase, the blue dashed line indicates the *v*2 of hadrons coming from strings while the brown dot-dashed line shows the *v*2 of hadrons coming from mesonic and baryonic resonance decays. Plot (**c**): individual contributions to *v*2 including their relative weights to the total *v*2.

**Figure 7.** Elliptic flow of charged hadrons as a function of *pT* for 0–50% central Au+Au collisions at √*sNN* = 200 GeV. The line description is similar to Figure 6.

The channel composition of *v*2 for PHSD5.0—with cross sections and parton masses evaluated at the actual chemical potential *μB* in each individual space-time cell—is shown in the middle plots '(b)' of Figures 6 and 7. We sorted the particles according to their production channels into three parts: the red dotted line corresponds to the final hadrons coming from the QGP without rescattering in the hadronic phase, the blue dashed line indicates the *v*2 of hadrons coming from strings (without further rescattering) while the brown dot-dashed line shows the *v*2 of hadrons coming from mesonic and baryonic resonance decays. One can see a large difference between the averaged elliptic flow for the different channels: the *v*2 of hadrons from string decay is the lowest since string production occurs dominantly at the initial phase of the heavy-ion collision; the *v*2 of hadrons from the QGP is the largest versus *η* as follows from the middle plot '(b)' of Figure 6. However, this is mainly due to the low *pT* hadrons which give a larger contribution to *<sup>v</sup>*2(*η*)—cf. the middle part '(b)' of Figure 7. Here, the high *pT* hadrons from the QGP show a lower *v*2 than those coming from strings or resonance decays.

The right parts '(c)' of Figures 6 and 7 present the individual contributions to *v*2 including their relative weights to the total *v*2. It shows that the properly weighted channel decomposition of *v*2 looks rather different—the contribution of the hadrons from the QGP is now small since most of them rescatter in the hadronic phase, i.e., the relative fraction of hadrons directly coming from QGP

hadronization is very small. The total *v*2 is dominated by the hadrons coming from the decay of resonances. The fraction of hadrons from string decays is very small due to the fact that strings are formed mainly in the beginning of collisions, and a very small fraction of hadrons can survive directly. Thus, the information in *v*2 about the QGP properties is washed out to a large extent by final hadronic interactions.

In Figure 8, we present the elliptic flow of identified hadrons (*K*±, *p*, *p*¯, Λ + Σ0, Λ¯ + Σ¯ 0) as a function of *pT* at √*sNN* = 27 GeV for PHSD4.0 (green lines), PHSD5.0 with partonic cross sections and parton masses calculated for *μB* = 0 (blue dashed lines) and with cross sections and parton masses evaluated at the actual chemical potential *μB* in each individual space-time cell (red lines) in comparison to the experimental data of the STAR Collaboration [44]. Similar to the directed flow shown in Figure 5, the elliptic flow from all three cases for PHSD shows a rather similar behavior, the differences are very small (within the statistics achieved here). Only antiprotons and antihyperons show a small decrease of *v*2 at larger *pT* for PHSD5.0 compared to PHSD4.0, which can be attributed to the explicit √*s*-dependence and different angular distribution of partonic cross sections in the PHSD5.0. We note that the underestimation of *v*2 for protons and Λ's we attribute to the details of the hadronic vector potential involved in this calculations which seems to underestimate repulsion.

**Figure 8.** Elliptic flow of identified hadrons (*K*±, *p*, *p*¯, Λ + Σ0, Λ¯ + Σ¯ 0) as a function of *pT* at √*sNN* = 27 GeV for PHSD4.0 (green lines), PHSD5.0 with partonic cross sections and parton masses calculated for *μB* = 0 (blue dashed lines) and with cross sections and parton masses evaluated at the actual chemical potential *μB* in each individual space-time cell (red lines) in comparison to the experimental data of the STAR Collaboration [44].
