Introduction

The exploration of quantum chromodynamics (QCD) phase diagrams and nuclear matter properties has been one of the most important motivations for relativistic heavy-ion collisions[1-6]. The region of high temperature and vanishing baryon chemical potential of the QCD phase diagram, where the phase transition from hadronic matter to quark-gluon plasma (QGP) [7] is a smooth crossover, has been well studied in experiments at the RHIC [8, 9, 10] and LHC [11, 12]. Theoretical calculations suggest that there may exist a critical end point in the finite baryon chemical potential region [13]. It is well known that several observables, such as directed flow [14-17] and moments of distributions of conserved charges [19], have been proposed to be sensitive to critical behavior. In the RHIC Beam Energy Scan phase I program (BES-I), the non-monotonic behavior of the net proton directed flow slope with respect to rapidity as a function of collision energy is observed, and the minimum of the distribution suggests the softest point of the equation of state (EoS) or prediction of the critical end point [18].

The collective flow [20, 21] has been extensively used to study the transverse properties of hot and dense matter created in heavy-ion collisions owing to their sensitivity to expansion dynamics. They are defined by the coefficients of the Fourier expansion of the azimuthal distribution of the emitted particles with respect to the reaction plane: [22] $E\frac{{\text{d}}^{3}N}{{\text{d}}^{3}p}=\frac{1}{2\pi }\frac{{\text{d}}^{2}N}{p_{\text{T}}\text{d}{p}_{\text{T}}\text{d}y}\left[1+2{\displaystyle \sum _{n=1}^{\infty }}{v}_n\cos \left[n\left(\varphi -{\Psi }_{\text{RP}}\right)\right]\right]\mathrm{.}$ (1) Here, ϕ is the azimuthal angle of the emitted particles, and Ψ_RP is the reaction plane angle. In this simulation, the reaction plane angle, Ψ_RP, is zero. The coefficients, vn, can then be determined by $v_n=\langle \cos [n(\varphi )]\rangle$, where the average runs over all particles in all events. The first two coefficients, v₁ and v₂, are of particular interest as they are established early during the system evolution. Recently, the STAR experiment has reported new results for collective flow at $\sqrt{s_{\text{NN}}}$ = 3 GeV [23, 24], which corresponds to a region with a high baryon density of 750 MeV [25]. The data suggest that hadronic interactions dominate the EoS in collisions at such energies. Exploring the phase structure in the high baryon density region is the main task of the second phase of the beam energy scan at the RHIC (BES-II) program [26], as well as the motivation for future facilities: such as the facility for proton and ion research (FAIR) [27], nuclotron-based ion collider facility (NICA) [28], and high-intensity heavy-ion accelerator facility (HIAF) [29, 30]. Regarding collectivity concerns, the responses of baryons and mesons are quite different. For example, we will focus on the collective measurements of protons, representative of baryons, and π⁺s, representative of produced mesons. In the case of pions, if the effect of isospin is neglected, there should be no differences in the flow of π⁺ and π^-. Calculations of Au+Au collisions from a hadronic transport model UrQMD [31, 32] will be used in this study.

This paper is organized as follows: The basic features of the UrQMD model are briefly discussed in Sect. 2. In Sect. 3, we present a v₁ and v₂ comparison between the UrQMD calculations and recent STAR preliminary results at $\sqrt{s_{\text{NN}}}$ = 3 GeV. The results of the energy dependence of v₁ and v₂ from UrQMD are also discussed. Finally, a summary is provided in Sect. 4.

The UrQMD Model

The ultra-relativistic quantum molecular dynamics model (UrQMD) [31, 32] is a microscopic transport model for simulating the process of relativistic heavy-ion collisions. It has been widely and successfully used for studying pp, pA, and AA collisions in heavy-ion collision physics over a broad energy range from a few GeV to CERN LHC energies. The particle production in UrQMD is similar to other transport models, which include resonance excitation and decay, as well as string dynamics and strangeness exchange reactions. The propagation of these hadrons and their exacted states is based on binary elastic and inelastic scattering, in which the scattering cross-sections are obtained from experimental data and models. The cascade version of the UrQMD model successfully describes the particle production, as well as directed flow measurements in heavy-ion collisions at collision energies $\sqrt{s_{\text{NN}}}$ > 7 GeV [33, 34]. However, at lower collision energies, it is necessary to include the mean-field potential to describe the collective flow in the high baryon density region [36, 35, 37, 38, 39]. A baryonic skyrme-potential [36, 35], based on the relative baryon density, is popularly used for heavy-ion collisions, which defines the nuclear EoS: $U=\alpha \left(\frac{\rho }{{\rho }_0}\right)+\beta {\left(\frac{\rho }{{\rho }_0}\right)}^{\gamma }\mathrm{.}$ (2) Here, ρ and ρ₀ are the baryon density is the ground-state and baryon density, respectively. By changing the parameters α, β, and γ, the stiffness of the nuclear EoS can be controlled. In this study, the parameters α, β, and γ are -124 MeV, 71 MeV, and 2, respectively, corresponding to a nuclear incompressibility of κ = 380 MeV. This value of κ was also used in this study [38], which is consistent with the HADES collective flow measurements at $\sqrt{s_{\text{NN}}}$ = 2.4 GeV [44]. Note that the value of κ depends on the collision energy and can be model dependent. As discussed in Ref. [39-43], the extracted values of κ are below 300 MeV for Au+Au collisions at lower energies. The momentum-dependent potential is not included in the mean field in this study.

In this work, we use version 3.4 of the UrQMD model to generate Monte Carlo event samples of Au+Au collisions at center of mass energies of $\sqrt{s_{\text{NN}}}$ = 2.11, 2.22, 2.32, 2.4, 2.51, 2.7, 2.86, 3.0, 3.5, 4.5 and 4.9 GeV, corresponding to beam kinetic energies of 0.5, 0.7, 1.0, 1.23, 1.5, 2.0, 2.5, 2.91, 4.65, 8.9 and 11 GeV per nucleon, respectively. Both cascade and mean-field modes of UrQMD are used to systematically study the collective flow in the high baryon-density region.

Results and Discussions

Following the STAR experimental analysis procedures, the collision centrality in this study is determined by the reference multiplicity, which is the number of charged pions, kaons, and protons within the pseudorapidity range $\left|\eta \right|$ < 0.5 [45]. In Fig. 1, the reference multiplicity distribution in Au+Au collisions with the UrQMD cascade version from the center-of-mass energy $\sqrt{s_{\text{NN}}}$ = 2.11 – 4.9 GeV is shown. These hadrons are produced more copiously at higher collision energies. The collision centrality bins for each energy can be determined by fitting the Glauber model [46]. In the following discussion, we report results from the 10–40% centrality bin where the collectivity, is expected to have the maximum strength, especially the second order of the Fourier coefficient v₂.

Figure 1Full-screen View

(Color online) Reference multiplicity distribution in UrQMD with cascade mode in Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 2.11–4.9 GeV. Similar to the experimental setup, the charged particles were measured within $\left|\eta \right|<$ 0.5.

According to the recent STAR collective flow measurements at $\sqrt{s_{\text{NN}}}$ = 3 GeV [23, 24], the UrQMD model results from both mean-field and cascade modes are compared with experimental data for the same centrality interval and kinematic selection criteria. Figure 2 presents the rapidity dependence of v₁ and v₂ for proton and π⁺ in 10–40% centrality Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 3 GeV. The symbols denote the STAR data. The red and blue bands represent the results from the UrQMD mean-field and cascade mode, respectively. Due to the experimental acceptance, only the negative rapidity range results are shown [23].

Figure 2Full-screen View

(Color online) Rapidity dependence of v₁ (top panel) and v₂ (bottom panel) of the proton and π⁺ in the 10–40% centrality Au+Au collision at $\sqrt{s_{\text{NN}}}$ = 3 GeV. The symbols denote STAR data. The red and blue bands represent the results from the UrQMD mean-field and cascade mode, respectively.

Negative v₂ values were observed at mid-rapidity for both protons and π⁺ at $\sqrt{s_{\text{NN}}}$ = 3 GeV, indicating the shadowing effect of the passing spectators. Meanwhile, the proton v₂ is positive in the forward rapidity region. For the v₁ results, a strong v₁ signal was observed compared to the high-energy results [47]. The proton v₁ increases from the forward to mid-rapidity region, whereas π⁺ v₁ exhibited the opposite trend. The standard UrQMD cascade mode failed to describe the experimental data. However, by including the mean-field potential, the UrQMD results reproduce the rapidity dependence and negative v₂ values at mid-rapidity for both protons and π⁺. These agreements between the experimental data and UrQMD calculations with mean-field potential imply that hadronic interactions play an important role in the collective flow development at such collision energies. Similar to Fig. 2, we present the p_T dependence of v₁ and v₂ for protons and π⁺ in 10–40% centrality Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 3 GeV in Fig. 3. The mid-rapidity range of -0.5 < y < 0 for protons and π⁺ are used in this analysis, which is consistent with the STAR experimental selection. Both v₁ and v₂ are negative and decrease with increasing p_T. Again, the red bands shown in the figure indicate the UrQMD results with mean-field potential, which can reproduce the p_T dependence and qualitatively describe the experimental data. However, quantitative differences are visible in some p_T bins for both protons and π⁺.

Figure 3Full-screen View

(Color online) Transverse momentum dependence of v₁ (top panel) and v₂ (bottom panel) of the proton and π⁺ in the 10–40% centrality Au+Au collision at $\sqrt{s_{\text{NN}}}$ = 3 GeV. The symbols denote STAR data. The red and blue bands represent the results from the UrQMD mean-field and cascade mode, respectively.

We systematically studied the mean-field potential effect on the particle production, directed flow, and elliptic flow development in the high baryon density region. The particle production yield as a function of kinetic freeze-out time for protons and π⁺ in 3 GeV Au+Au central collisions (0–10%) at mid-rapidity ($\left|y\right|<$ 0.5) is presented in Fig. 4. The red and blue lines represent the results from the UrQMD mean-field and cascade modes, respectively. Note that these two distributions were normalized by the number of events. Overall, hadrons are produced mostly during kinetic freeze-out in the time interval 10–25 fm/c. The mean value of the distribution is shown in each panel using the corresponding colors. It can be seen that the mean value from the mean-field version is smaller than that of the cascade version. This can be explained by the larger pressure generated by the density dependent fields, which drives the strong expansion. Thus, these hadrons are pushed outside the system and freeze out earlier than those from the cascade mode.

Figure 4Full-screen View

(Color online) Particle production yield as a function of freeze-out time for proton (left panel) and π⁺ (right panel) from UrQMD 3 GeV 0–10% Au+Au collisions. The red and blue lines represent the results from the mean-field and cascade modes, respectively.

To quantify the strength of directed flow at mid-rapidity, a linear fit is employed to extract the slope of the v₁(y) distribution [47, 18], as shown in Fig. 2. The directed flow is exactly zero at y = 0 owing to the momentum conservation effect. Figure 5 presents the v₁ slopes $\text{d}{v}_1/\text{d}y{|}_{y=0}$ and elliptic flow v₂ at mid-rapidity as a function of kinetic freeze-out time (t) for proton and π⁺ in 3 GeV Au+Au mid-central collisions (10–40%), the red and blue bands represent the results from the UrQMD mean-field and cascade mode, respectively. Both v₁ and v₂ are established in the very early stage of the system evolution and show a rapid increase for t ≲ 15 fm/c. They then increased slightly with increasing freeze-out time. As can be seen in the figure, the v₁ values for the mid-rapidity proton and pion start from below zero before turning into positive values at a later time, as experimentally observed. Owing to the baryonic mean-field, at t $\tilde{1}$ 5 fm/c, the v₁ slope for the proton is much larger than that of π⁺, especially in the results for the mean-field mode. Conversely, while the cascade predicts an in-plane expansion, for example, positive v₂ for both protons and pions, the mean-field tends to keep the incoming nucleons together. As a result, shadowing seems to be sustained throughout the collision period, leading to negative values of v₂ and out-of-plane expansion at this energy. In short, compared to the cascade mode, the mean-field potential provides a stronger expansion in the x-z plane (v₁) and an effective blocking that leads to the expansion in the y-z plane (v₂) in Au+Au collisions at 3 GeV.

Figure 5Full-screen View

(Color online) v₁ slopes ($\text{d}{v}_1/\text{d}y{|}_{y=0}$) and v₂ at mid-rapidity as a function of freeze-out time for proton (left panel) and π⁺ (right panel) from UrQMD 3 GeV 10–40% Au+Au collisions. The red and blue bands represent the results from the mean-field and cascade modes, respectively.

The mean transverse momentum $\langle p_{\text{T}}\rangle$ as a function of collision energy $\sqrt{s_{\text{NN}}}$ at a mid-rapidity of $\left|y\right|<$ 0.5 for protons and π⁺ in UrQMD 10–40% Au+Au collisions is shown in Fig. 6. The red and blue symbols represent the results from the cascade and mean-field modes, respectively. It can be observed that $\langle p_{\text{T}}\rangle$ increases with increasing collision energy for both particle species owing to the higher energies distributed in the transverse direction. The mean-field potential slightly increases $\langle p_{\text{T}}\rangle$ for all particles as a result of the increased radial flow. The increase in $\langle p_{\text{T}}\rangle$ is greater for heavier particle protons. The enhancement in the mean-field mode is due to the more repulsive interactions in the high-baryon-density region.

Figure 6Full-screen View

(Color online) Mean transverse momentum $\langle p_{\text{T}}\rangle$ as a function of collision energy for proton (squares) and π⁺ (circles) from UrQMD 10–40% Au+Au collisions. The red and blue symbols represent the results from the mean-field and cascade modes, respectively.

We will now discuss the energy dependence of directed and elliptic flows from UrQMD and compare this with world experimental data. Figure 7 shows the energy dependence of the mid-rapidity slope of directed flow ($\text{d}{v}_1/\text{d}y{|}_{y=0}$) in 10–40% centrality Au+Au collisions. The transverse momentum range in the UrQMD calculations for proton and π⁺ is 0.4 < p_T < 2.0 GeV/c and 0.2 < p_T < 1.6 GeV/c, respectively, which is consistent with the STAR 3 GeV results [23]. The black symbols denote the world experimental data [48, 50, 23]. The red and blue symbols represent the calculations from the UrQMD mean-field and cascade modes, respectively. For the π⁺ v₁ and v₂ measurements, only the STAR results from $\sqrt{s_{\text{NN}}}$ = 3 and 4.5 GeV are available. The proton $\text{d}{v}_1/\text{d}y{|}_{y=0}$ values increase with decreasing collision energy for both the UrQMD cascade and mean-field calculations. While the mean-field potential significantly enhances the $\text{d}{v}_1/\text{d}y{|}_{y=0}$ values, allowing comparison with the experimental data for both protons and π⁺. However, there are still some differences relating to the choice of the p_T range. As we can see in Fig. 2, v₁ has a strong p_T dependence at such a low energy.

Figure 7Full-screen View

(Color online) v₁ slopes ($d{v}_1/dy{|}_{y=0}$) at mid-rapidity as a function of collision energy for proton (left panel) and π⁺ (right panel) from UrQMD 10–40% Au+Au collisions. The red and blue symbols represent the results from the mean-field and cascade modes, respectively. World experimental data are represented by black symbols.

The energy dependence of the elliptic flow at mid-rapidity in 10–40% centrality Au+Au collisions is shown in Fig. 8. The black symbols denote the world experimental data [49, 50, 23]. Both the sign and absolute value of v₂ reflect the medium properties. v₂ increases as a function of the collision energy. The results from the cascade mode are always above the results from the mean-field mode, but both converge at a high energy of approximately 5 GeV. At approximately 2.4 GeV, the proton v₂ has already changed from a negative to positive value in the cascade mode, and the change occurs at a much higher energy in the case of the mean-field calculation. According to world data, the proton elliptic flow transits from in-plane expansion to out-of-plane expansion at $\sqrt{s_{\text{NN}}}$ ∼ 3.6 GeV as a result of the shadowing effect, where particles are blocked in the reaction plane by the spectators and are emitted mainly in the out-of-plane direction. In Fig. 8, the UrQMD calculations with the mean-field potential qualitatively describe the experimental data for both protons and π⁺, and are quantitatively consistent with the data for $\sqrt{s_{\text{NN}}}<$ 3.6 GeV.

Figure 8Full-screen View

(Color online) Similar to Fig. 7: v₂ at mid-rapidity as a function of the collision energy for protons (left panel) and π⁺ (right panel) from UrQMD to 10–40% Au+Au collisions.

Summary

A hadronic transport model, UrQMD, was employed to study the directed and elliptic flows in the high baryon density region. The model provides a good description of the recent STAR results of 10–40% Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 3 GeV by including the mean-field potential. This indicates that hadronic interactions play an important role in the collective flow development in the high-baryon-density region. Other measurements from the 3 GeV Au+Au collisions, including strangeness production [51] and high moments of protons [52], also suggest that hadronic interactions dominate the dynamics in the high baryon density region.

The results of the time evolution of directed and elliptic flows show that the collective flow is fully developed at $\tilde{1}$ 0-15 (fm/c).

The mean-field potential slightly enhances $\langle p_{\text{T}}\rangle$ because of the stronger radial expansion, and shifts the freeze-out to an earlier time. In addition, compared to the cascade mode, the mean-field option delays the proton v₂ crossing zero as a function of collision energy, from $\sqrt{s_{\text{NN}}}$ = 2.5 to 3.5 GeV. A negative or positive v₂ corresponds to an out-of-plane or in-plane expansion, respectively.

These results will be useful for understanding experimental data from upcoming experimental facilities focused on the high baryon density region.