# Radial Flow in a Multiphase Transport Model at FAIR Energies

Authors: Soumya Sarkar; Provash Mali; Somnath Ghosh; Amitabha Mukhopadhyay
Journal: Advances in High Energy Physics
 (2018)

---

## Abstract

Azimuthal distributions of radial velocities of charged hadrons produced in nucleus-nucleus(AB) collisions are compared with the corresponding azimuthal distribution of charged hadron multiplicity in the framework of a multiphase transport (AMPT) model at two different collision energies. The mean radial velocity seems to be a good probe for studying radial expansion. While the anisotropic parts of the distributions indicate a kind of collective nature in the radial expansion of the intermediate “fireball,” their isotropic parts characterize a thermal motion. The present investigation is carried out keeping the upcoming Compressed Baryonic Matter (CBM) experiment to be held at the Facility for Antiproton and Ion Research (FAIR) in mind. As far as high-energy heavy-ion interactions are concerned, CBM will supplement the Relativistic Heavy-Ion Collider (RHIC) and Large Hadron Collider (LHC) experiments. In this context our simulation results at high baryochemical potential would be interesting, when scrutinized from the perspective of an almost baryon-free environment achieved at RHIC and LHC.


---

## Body

## 1. Introduction

When two heavy nuclei collide with each other at high-energy it is expected that a color deconfined state composed of strongly coupled quarks and gluons is formed. The properties of such a state, formally known as Quark-Gluon Plasma (QGP) [1], are governed by the rules of quantum chromodynamics (QCD). In order to understand the bulk properties of this extended QCD state and to understand the dynamical processes that might be involved in its formation and subsequent decay, over last three decades or so QGP has been widely searched in many high-energy experiments [2]. Of all the efforts in this regard, the study on the emission of final state particles with respect to the reaction plane of an AB collision, a plane spanned by the beam direction and the impact parameter vector, has been a popular technique that can characterize the thermodynamic and hydrodynamic properties of the QGP matter [3, 4]. In this technique the Fourier decomposition of the anisotropic azimuthal distribution has been widely employed to explore the collective behavior of the final state particles. More specifically, the second harmonic coefficient, traditionally known as the elliptic flow parameter (v2), is of special interest [5]. The v2 results obtained from the RHIC and LHC experiments show considerable hydrodynamical behavior of the matter present in the overlapping zone of the colliding nuclei, an intermediate “fireball” that gets thermalized within a very short time interval (<1 fm/c) and subsequently expands almost like a “perfect fluid” having a very small ratio of shear viscosity to entropy density [6–8]. In RHIC [9–13] and LHC [14–17] experiments the v2 parameter has been widely studied as a function of the centrality of the collision, transverse momentum (pT), and rapidity (y) or pseudorapidity (η) of the produced particles, for different colliding systems and at varying collision energies. Using the AMPT model, presence of anisotropy in the azimuthal distribution of transverse rapidity (yT) has been investigated in Au + Au collision at sNN=200 GeV [18]. Recently, we have reported some simulation results on the anisotropy present in the azimuthal distribution of pT of the emitted charged hadrons that has relevance to the radial flow of charged hadrons produced at FAIR energies [19]. The CBM experiment is dedicated to explore the deconfined QCD matter at high baryon density and low to moderate temperature. It is a fixed target experiment being designed with incident beam energy range Elab=10–40 GeV per nucleon, which is expected to produce partonic matter of density 6 to 12 times the normal nuclear matter density at the central rapidity region [20]. But at the same time it should be kept in mind that our present understanding of the collective flow of hadronic/partonic matter at this energy region is constrained by the availability of only a few experimental results [21]. Therefore, to get an idea about the expected behavior of any variable or parameter that is relevant in this regard, we have to rely mostly upon the event generators and models. In this article we present some basic simulated results on the (an)isotropy of the radial velocity of charged hadrons produced in Au + Au collisions at Elab=10A and 40A GeV using the AMPT model [22, 23]. The paper is organized as follows: a brief description of the methodology used in this analysis is given in Section 2, the AMPT model is summarily described in Section 3, our simulation based results are discussed in Section 4, and finally in Section 5 our observations are listed.


## 2. Methodology

Before the collision, the nucleons belonging to individual nucleus possess only longitudinal degrees of freedom. Transverse degrees of freedom are excited into them only after an interaction takes place. In mid-central collisions the overlapping area of the colliding nuclei is almond shaped in the transverse plane. This initial asymmetry in the geometrical shape gives rise to different pressure gradients along the long and the short axis of the overlapping zone and correspondingly to a momentum space asymmetry in the final state. As a result, if the matter present in the intermediate “fireball” exhibits a fluid-like behavior, then a collective flow of final state particles is observed, which is reflected in the azimuthal distribution of the particle number as well as in the azimuthal distribution of an appropriate kinematic variable likepT, yT and the transverse or radial velocity vT [24]. The radial velocity has two components, the radial flow velocity and the velocity due to the random thermal motion of the particles constituting the intermediate fireball. For an ideal fluid the radial flow velocity should be isotropic. However, for a nonideal viscous fluid, the shear tension is proportional to the gradient of the radial velocity along the azimuthal direction, which again is related to the anisotropy of radial velocity [18]. An analysis of the v2 results over the energy range Elab≅1–160 GeV has shown that the observed v2 values are lower than what is expected from a phenomenology based on the three-fluid dynamics [25]. The difference has been attributed to dissipative effects like viscosity. A single parameter (the Knudsen number) fit of the v2 results over a wide range of collision energy suggests that the upper limit of the shear viscosity to specific entropy ratio η/s~1-2, a value much higher than what is estimated for an almost ideal fluid created at RHIC or LHC energies. However, to understand the exact nature of the flow characteristics or the nature of the fluid created at FAIR energies, we shall have to wait till the CBM results in this regard become available.We introduce the transverse (radial) velocity as(1)vT=pTE=pTmTcosh⁡y,where E=mTcosh⁡y is the energy of the particle, mT=m02+pT2 is its transverse mass, m0 is the particle rest mass, and y is its rapidity. For a large sample of events the total radial velocity 〈VT(ϕm)〉 of all particles falling within the mth azimuthal bin is defined as (2)VTϕm=1Nev∑j=1Nev∑i=1nmvT,iϕm,where vT,i(ϕm) is the radial velocity of the ith particle, nm is the total number of particles present in the mth bin, Nev is the number of events under consideration, and  denotes an averaging over events. In this paper we have chosen the transverse velocity as the basic variable in terms of which the azimuthal asymmetry has been studied and compared the results obtained thereof with those of the azimuthal asymmetry associated with the charged particle multiplicity distribution. An azimuthal distribution of 〈VT(ϕm)〉 contains information of the asymmetry in the multiplicity distribution as well as that of the radial expansion. By taking an average over the particle number the mean transverse velocity 〈〈vT(ϕm)〉〉 is introduced as(3)vTϕm=1Nev∑j=1Nev1Nm∑i=1NmvT,iϕm,where  represents first an average over all particles present in the mth azimuthal bin and then over all events present in the sample. This double averaging reduces the multiplicity influences significantly, and the corresponding distribution measures only the radial expansion. In this context we must mention that the mean radial velocity actually consists of contributions coming from three different sources, the average isotropic radial velocity, the average anisotropic radial velocity, and the average velocity associated with thermal motion. It should be noted that both radial and thermal motion contribute to the isotropic velocity of the distribution. Like the azimuthal distribution of charged particle multiplicity dNch/dϕ, it is also possible to expand the azimuthal distributions of the total and mean transverse velocities in Fourier series as(4)dVTdϕ≈v0VT1+2v2VTcos⁡2ϕ,dvTdϕ≈v0vT1+2v2vTcos⁡2ϕ.In these expansions only the leading order terms (n= 0 and 2) are retained. The anisotropy present in any of the distributions [see (4)] is quantified by the second Fourier coefficient v2, whereas v0 is a measurement of the isotropic flow.


## 3. AMPT Model

Transport models are best suited to studyAB collisions at the energy range under our consideration. Since transport models treat chemical and thermal freeze-out dynamically, they have the ability to describe the space-time evolution of the hot and dense “fireballs” created in collision between two heavy nuclei at relativistic energy. As mentioned before, in this simulation study we use the AMPT model with partonic degrees of freedom, the so-called string melting version, with the expectation that under the FAIR conditions transitions from the initial nuclear matter to the QGP state (if any) and then from the QGP state to the final hadronic state will take place at high baryon density and low to moderate temperature. Previous calculations have shown that flow parameters consistent with experiment can be developed through AMPT and the model can successfully describe different aspects of the collective behavior of hadronic/partonic matter produced in AB interactions [26–32]. The string melting version of AMPT should be even more appropriate to model particle emission data, where a transition from nuclear matter to deconfined QCD state is expected. AMPT is a hybrid model where the primary particle distribution and other initial conditions are taken from the heavy-ion jet interaction generator (HIJING) [33], and Zhang’s parton cascade (ZPC) formalism [34] is used in subsequent stages. Note that the ZPC model includes only parton-parton elastic scattering with an in-medium cross section derived from pQCD, the effective gluon screening mass being taken as an adjustable parameter. In the string melting version of the AMPT model, all hadrons are produced from string fragmentation like that in the HIJING model. The strings are converted into valence quarks and antiquarks. They are subsequently allowed to interact through the ZPC formalism and propagate according to a relativistic transport model [23]. Finally, the quarks and antiquarks are converted to hadrons via a quark coalescence formalism.


## 4. Results and Discussion

In this section we describe our results obtained from the Au + Au minimum bias events simulated by the AMPT model (string melting version) atElab=10A and 40A GeV. A representative value of the parton scattering cross section (σ=3 mb) is used in this analysis. The σ value is chosen so as to match with a previously studied collective behavior at FAIR energies [31]. We have indeed compared the NA49 results [21] on the pT dependence of elliptic flow parameter v2 by varying σ over a range of 0.1 to 6 mb. We have seen that even though the σ values differ almost by two orders of magnitude, the corresponding differences in the simulated v2 values are not that significant [35]. The size of each sample of Au + Au events used in this analysis is one million. We begin with the multiplicity distribution of charged hadrons, represented schematically in Figure 1. The nature of the distributions is more or less similar at both energies considered. However, the average and the highest multiplicities are naturally quite larger at Elab=40A GeV. In Figure 2 we plot the pT distributions of charged hadrons. As expected, with increasing pT we observe an approximately exponential fall in the particle number density. It is interesting to note that, at low pT values, up to pT≈1.5 GeV/c, the slopes of the distributions at both energies hardly differ, but at high pT, beyond pT=2.0 GeV/c, the slope values are considerably different, this being stiffer at Elab=10A GeV. The inverse slope can be related to the temperature of the intermediate “fireball” as and when it achieves a thermal equilibrium. With varying collision energies, the particles produced in soft processes, therefore, correspond to almost same source temperature, irrespective of the collision energy. A Monte Carlo Glauber (MCG) model [36] is employed to characterize the geometry of an AB collision. Using the MCG model the average transverse momentum 〈pT〉 of particles produced in AB collisions belonging to a particular centrality can be determined. In Figure 3 such a plot of 〈pT〉 against the number of participating nucleons (Npart), a measure of the centrality of the collision, is graphically shown. At both the collision energies considered in this analysis, at low Npart the 〈pT〉 values are significantly different, 〈pT〉 increases almost linearly with increasing centrality, and each distribution saturates at similar value, 〈pT〉≈0.365 GeV/c. The fact that at the highest centrality the saturation value of pT of produced particles is almost independent of the incident beam energy is perhaps due to kinematic reasons. The transverse degree of freedom that was absent before the collision took place is excited into the interacting AB system due to multiple nucleon-nucleon (NN) scattering and rescattering. Our results indicate that the degree of such excitations, which predominantly depends on the number of binary collisions Ncoll, appears to remain almost same for the most central collisions of the Au + Au system in the FAIR energy region. In Figure 4 we present the azimuthal distributions of (a) the total radial velocity 〈VT〉, (b) the multiplicity (Nch), and (c) the mean radial velocity 〈〈vT〉〉 of charged hadrons produced at Elab=40A GeV in the mid-rapidity region, within Δy=±1.0 symmetric about the central value y0 and within the 0–80% centrality range. Presence of anisotropy in all three distributions is clearly visible. It is also observed that while all three distributions exhibit same periodicity, their amplitudes are quite different. In order to show that all three distributions can analytically be described by a single function like N1+αcos⁡(2ϕ) without significant contributions coming from other harmonics, we fit the distributions with exactly the same relative vertical axis range with respect to the value of the parameter α centred around the same value of the other parameter N (here N=1.0) and plot them together in Figure 4(d) along with the respective fitted lines. When appropriately scaled, we find that the elliptic anisotropy present in the distribution of total radial velocity is almost equal in magnitude to that coming from the anisotropy in multiplicity distribution. In comparison, corresponding anisotropy in the mean radial velocity is quite small. The results at Elab=10A and 40A GeV are qualitatively similar.Figure 1
Charged hadron multiplicity distribution in Au + Au collisions atElab=10A and 40A GeV.Figure 2
Charged hadronpT distribution in Au + Au collisions at Elab=10A and 40A GeV.Figure 3
Average transverse momentum of charged hadrons as a function ofNpart in Au + Au collisions at Elab=10A and 40A GeV.Figure 4
Azimuthal distribution of (a) total radial velocity, (b) multiplicity, (c) mean radial velocity, and (d) all the aforesaid quantities properly normalized for charged hadrons produced in Au + Au collisions atElab=40A GeV.
(a) (b) (c) (d)


### 4.1. Centrality Dependence ofv2 and v0

Elliptic flow results from interactions between particles comprising the intermediate “fireball,” and hence it is a useful probe for the identification of local thermodynamic equilibrium. Thev2 values are smaller for the extreme central and peripheral collisions, which can be explained in terms of the initial geometric effects and the pressure gradient produced thereof [37]. In the hydrodynamical limit v2 is proportional to the elliptic eccentricity (ε2) of the overlapping region of the colliding nuclei, whereas in the low density limit v2/ε2 is proportional to the product of the rapidity density of charged particles dNch/dy and inverse of the overlapping area of the colliding nuclei. It is believed that the centrality dependence of elliptic flow provides valuable information regarding the degree of equilibration achieved by the intermediate “fireball” and also regarding the characteristics of (re)scattering effects present therein [38]. Some model based results at FAIR energies can be found in [19, 31, 35].In Figure5 we compare the centrality dependence of the v2 parameter obtained from distributions of all three variables under consideration. The overall centrality dependence is found to be similar for v2(VT) and v2(Nch). However, v2vT, which in magnitude is quite small in comparison with the v2 values obtained from the other two variables, behaves quite differently. All three variations are consistent with our observations of Figure 4(d). It is to be noted that the anisotropy in mean radial velocity, which describes the radial expansion, is significantly smaller than that of the corresponding multiplicity distribution in the mid-central region. In this regard we also intend to scrutinize the effects of the collision energy involved. It is observed that v2VT and v2(〈Nch〉) at Elab=40A GeV are marginally higher than those obtained at Elab=10A GeV, a general feature of any v2 result, which has been confirmed over a wide energy range. The v2(〈〈vT〉〉) values are not significantly different at the two collision energies involved either. We expect that the isotropy parameter (v0) of all aforesaid distributions is also of certain importance and we graphically plot the results in Figure 6. The v0 values associated with 〈VT〉 and 〈Nch〉 distributions show a linear dependence with increasing Npart, being highest in the most central events. This feature of v0 can be ascribed to the fact that the azimuthally integrated magnitude of transverse flow increases with increasing centrality of the collisions. On the contrary, an increasing trend in the v0vT values with increasing Npart is restricted only to the peripheral collisions, and beyond Npart≈80 the v0(〈〈vT〉〉) values achieve a saturation, being nearly independent of the centrality of the collisions. A significant energy dependence of v0 is also observed for all three variables considered in this analysis. We do not see any significant energy dependence in the variation of v0(〈〈VT〉〉) with Npart. The v0(〈Nch〉) values are however consistently higher at Elab=40A GeV than those at Elab=10A GeV, the difference becoming larger with increasing Npart. Once again v0(〈〈vT〉〉) behaves quite differently in this regard. The values at Elab=10A GeV are consistently higher than those at Elab=40A GeV. We may recall that the mean radial velocity has been defined in a way such that the multiplicity effects are removed. Therefore, we conclude that the particle multiplicity plays a dominant role to determine the total transverse flow, and a higher energy input results in a lower amount of azimuthally integrated transverse flow.Figure 5
Centrality dependence of anisotropy parameterv2 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.Figure 6
Centrality dependence of isotropic flow coefficientv0 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.
(a) (b) (c)


### 4.2. Transverse Momentum Dependence ofv2 and v0

It is well known that the anisotropy coefficientv2 depends on pT of charged hadrons. Hydrodynamics as well as resonance decay are expected to dominate at low pT, whereas high pT particles are expected to stem out from the fragmentation of jets modified in the hot and dense medium of the intermediate “fireball” [15]. At FAIR energies the production of high pT hadrons would be rare, and owing to statistical reasons we restrict our analysis up to pT=2.0 GeV/c. v2 arising from multiplicity distributions of the produced hadrons has been studied widely as a function of pT using the data available from the experiments held at RHIC [39] and LHC [40]. Simulation results under FAIR-CBM conditions utilizing the UrQMD, AMPT (default), and AMPT (string melting) models can be found in [19, 31]. Figure 7 depicts that the anisotropies present in 〈Nch〉, 〈VT〉, and 〈〈vT〉〉 distributions rise monotonically with increasing pT. At Elab=40A GeV, beyond pT=1.5GeV/c, there is a trend of saturation in the v2 values extracted from all three variables. Once again we conclude that the multiplicity dominates over the radial velocity at a particular pT bin, and v2(〈Nch〉) and v2(〈VT〉) are found to be almost equal in the 0⩽pT⩽2.0 GeV/c range. Once we get rid of the multiplicity effects, the actual anisotropy present in the radial velocity comes out, which we can see in the plot of v2(〈〈vT〉〉) against pT shown in the same diagram. As a result, within 0.25≤pT≤1.25 GeV/c the v2(〈〈vT〉〉) values are slightly lower at higher Elab. At FAIR energies, however, we do not find any noticeable deviation in the trend of this pT dependence of v2 from its nature observed at RHIC energies [24]. Comparing Figure 7(a) with Figure 7(b), we see a very weak (almost insignificant) energy dependence of v2 in terms of all three variables concerned. We may reckon that FAIR-CBM energy range may not provide us with a proper platform to study the energy dependence of anisotropy, but it may be suitable for studying the issues related to the isotropy measure v0. The pT dependence of v0 has been shown in Figure 8. It is observed that the v0 coefficients associated with 〈Nch〉, 〈VT〉, and 〈〈vT〉〉 while being plotted against pT exhibit similar nature. In the low pT region, the v0 values extracted from each variable rise with increasing pT, attain a maximum, and then fall off to a very small saturation value (almost zero) at both incident energies beyond pT=1.25 GeV/c. Once again, while v0(〈VT〉) values at Elab=10A and 40A GeV are almost identical, the v0(〈Nch〉) values at Elab=40A GeV are higher in the low pT region (pT≲0.7 GeV/c) than those at 10A GeV. On the contrary, the v0(〈〈vT〉〉) values obtained at Elab=40A GeV are lower in the low pT region (pT≲0.5 GeV/c) than those at 10A GeV. At FAIR energy the random thermal motion of particles perhaps dominates over their collective behavior, which at high pT leads to a very small amount of azimuthally integrated magnitude of net flow.Figure 7
Transverse momentum dependence of anisotropy parameterv2 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.
(a)
10 A GeV (b)
40 A GeVFigure 8
Transverse momentum dependence of isotropic coefficientv0 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.
(a) (b) (c)


## 4.1. Centrality Dependence ofv2 and v0

Elliptic flow results from interactions between particles comprising the intermediate “fireball,” and hence it is a useful probe for the identification of local thermodynamic equilibrium. Thev2 values are smaller for the extreme central and peripheral collisions, which can be explained in terms of the initial geometric effects and the pressure gradient produced thereof [37]. In the hydrodynamical limit v2 is proportional to the elliptic eccentricity (ε2) of the overlapping region of the colliding nuclei, whereas in the low density limit v2/ε2 is proportional to the product of the rapidity density of charged particles dNch/dy and inverse of the overlapping area of the colliding nuclei. It is believed that the centrality dependence of elliptic flow provides valuable information regarding the degree of equilibration achieved by the intermediate “fireball” and also regarding the characteristics of (re)scattering effects present therein [38]. Some model based results at FAIR energies can be found in [19, 31, 35].In Figure5 we compare the centrality dependence of the v2 parameter obtained from distributions of all three variables under consideration. The overall centrality dependence is found to be similar for v2(VT) and v2(Nch). However, v2vT, which in magnitude is quite small in comparison with the v2 values obtained from the other two variables, behaves quite differently. All three variations are consistent with our observations of Figure 4(d). It is to be noted that the anisotropy in mean radial velocity, which describes the radial expansion, is significantly smaller than that of the corresponding multiplicity distribution in the mid-central region. In this regard we also intend to scrutinize the effects of the collision energy involved. It is observed that v2VT and v2(〈Nch〉) at Elab=40A GeV are marginally higher than those obtained at Elab=10A GeV, a general feature of any v2 result, which has been confirmed over a wide energy range. The v2(〈〈vT〉〉) values are not significantly different at the two collision energies involved either. We expect that the isotropy parameter (v0) of all aforesaid distributions is also of certain importance and we graphically plot the results in Figure 6. The v0 values associated with 〈VT〉 and 〈Nch〉 distributions show a linear dependence with increasing Npart, being highest in the most central events. This feature of v0 can be ascribed to the fact that the azimuthally integrated magnitude of transverse flow increases with increasing centrality of the collisions. On the contrary, an increasing trend in the v0vT values with increasing Npart is restricted only to the peripheral collisions, and beyond Npart≈80 the v0(〈〈vT〉〉) values achieve a saturation, being nearly independent of the centrality of the collisions. A significant energy dependence of v0 is also observed for all three variables considered in this analysis. We do not see any significant energy dependence in the variation of v0(〈〈VT〉〉) with Npart. The v0(〈Nch〉) values are however consistently higher at Elab=40A GeV than those at Elab=10A GeV, the difference becoming larger with increasing Npart. Once again v0(〈〈vT〉〉) behaves quite differently in this regard. The values at Elab=10A GeV are consistently higher than those at Elab=40A GeV. We may recall that the mean radial velocity has been defined in a way such that the multiplicity effects are removed. Therefore, we conclude that the particle multiplicity plays a dominant role to determine the total transverse flow, and a higher energy input results in a lower amount of azimuthally integrated transverse flow.Figure 5
Centrality dependence of anisotropy parameterv2 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.Figure 6
Centrality dependence of isotropic flow coefficientv0 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.
(a) (b) (c)


## 4.2. Transverse Momentum Dependence ofv2 and v0

It is well known that the anisotropy coefficientv2 depends on pT of charged hadrons. Hydrodynamics as well as resonance decay are expected to dominate at low pT, whereas high pT particles are expected to stem out from the fragmentation of jets modified in the hot and dense medium of the intermediate “fireball” [15]. At FAIR energies the production of high pT hadrons would be rare, and owing to statistical reasons we restrict our analysis up to pT=2.0 GeV/c. v2 arising from multiplicity distributions of the produced hadrons has been studied widely as a function of pT using the data available from the experiments held at RHIC [39] and LHC [40]. Simulation results under FAIR-CBM conditions utilizing the UrQMD, AMPT (default), and AMPT (string melting) models can be found in [19, 31]. Figure 7 depicts that the anisotropies present in 〈Nch〉, 〈VT〉, and 〈〈vT〉〉 distributions rise monotonically with increasing pT. At Elab=40A GeV, beyond pT=1.5GeV/c, there is a trend of saturation in the v2 values extracted from all three variables. Once again we conclude that the multiplicity dominates over the radial velocity at a particular pT bin, and v2(〈Nch〉) and v2(〈VT〉) are found to be almost equal in the 0⩽pT⩽2.0 GeV/c range. Once we get rid of the multiplicity effects, the actual anisotropy present in the radial velocity comes out, which we can see in the plot of v2(〈〈vT〉〉) against pT shown in the same diagram. As a result, within 0.25≤pT≤1.25 GeV/c the v2(〈〈vT〉〉) values are slightly lower at higher Elab. At FAIR energies, however, we do not find any noticeable deviation in the trend of this pT dependence of v2 from its nature observed at RHIC energies [24]. Comparing Figure 7(a) with Figure 7(b), we see a very weak (almost insignificant) energy dependence of v2 in terms of all three variables concerned. We may reckon that FAIR-CBM energy range may not provide us with a proper platform to study the energy dependence of anisotropy, but it may be suitable for studying the issues related to the isotropy measure v0. The pT dependence of v0 has been shown in Figure 8. It is observed that the v0 coefficients associated with 〈Nch〉, 〈VT〉, and 〈〈vT〉〉 while being plotted against pT exhibit similar nature. In the low pT region, the v0 values extracted from each variable rise with increasing pT, attain a maximum, and then fall off to a very small saturation value (almost zero) at both incident energies beyond pT=1.25 GeV/c. Once again, while v0(〈VT〉) values at Elab=10A and 40A GeV are almost identical, the v0(〈Nch〉) values at Elab=40A GeV are higher in the low pT region (pT≲0.7 GeV/c) than those at 10A GeV. On the contrary, the v0(〈〈vT〉〉) values obtained at Elab=40A GeV are lower in the low pT region (pT≲0.5 GeV/c) than those at 10A GeV. At FAIR energy the random thermal motion of particles perhaps dominates over their collective behavior, which at high pT leads to a very small amount of azimuthally integrated magnitude of net flow.Figure 7
Transverse momentum dependence of anisotropy parameterv2 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.
(a)
10 A GeV (b)
40 A GeVFigure 8
Transverse momentum dependence of isotropic coefficientv0 obtained from the azimuthal distributions of total radial velocity, multiplicity, and mean radial velocity in Au + Au collision at Elab=10A and 40A GeV.
(a) (b) (c)


## 5. Conclusion

In this paper we present some basic results on the elliptic and radial flow of charged hadrons. The study is based on the azimuthal distributions of total transverse velocity, mean transverse velocity, and multiplicity of charged hadrons produced in Au + Au collisions atElab=10A GeV and 40A GeV. We have used the AMPT model (string melting version) to generate the events. We observe that azimuthal asymmetries are indeed present in all three distributions. However, we also note that in our simulation results the azimuthal anisotropy of the final state particles is predominantly due to the asymmetry of particle multiplicity distribution, and only a small fraction of this asymmetry is due to kinematic reasons. The overall nature of the dependence of the elliptic anisotropy parameter on the centrality of the collision and transverse momentum of produced particles are similar for the three variables considered in the present analysis. The elliptic flow parameter is highest in the mid-central collisions, and within the interval 0⩽pT⩽2.0 GeV/c it is highest at the highest pT. From our simulated results in the FAIR energy range we find a very small energy dependence of the elliptic flow parameter. On the other hand, the azimuthally integrated magnitude of the radial flow is maximum for most central collisions and its values are high in the low pT region. From this analysis we see that the contribution to v0 from the asymmetry in multiplicity distribution and that coming from the asymmetry in kinematic variable vT exhibit an opposite incident beam energy dependence. While the former is higher at higher Elab, the latter is higher at lower Elab. Our simulated results are consistent with those obtained from RHIC and LHC energies and do not require any new dynamics to interpret. However, in future there is enough scope to appropriately model these results in terms of relevant thermodynamic and hydrodynamic parameters associated with the intermediate “fireball” produced in AB collisions.



---

*Source: 7453752-2018-04-22.xml*