Abstract

We have studied the equation of state and dissociation temperature of bottomonium state by correcting the full Cornell potential in isotropic medium by employing the effective fugacity quasi-particle Debye mass. We had also calculated the bottomonium suppression in an expanding, dissipative strongly interacting QGP medium produced in relativistic heavy-ion collisions. Finally we compared our results with experimental data from RHIC 200GeV/nucleon Au-Au collisions, LHC 2.76 TeV/nucleon Pb-Pb, and LHC 5.02 TeV/nucleon Pb-Pb collisions as a function of number of participants.

1. Introduction

At the Relativistic Heavy-Ion Collider (RHIC) situated at Brookhaven National Laboratory (BNL), heavy-ion collisions have been studied. After the pioneer work done in the direction of suppression by Matsui and Satz, and some other development of the potential models, suppression was observed by both SPS and RHIC [1]. Due to the Debye screening of the Quantum Chromo-Dynamic (QCD) potential between the two heavy quarks, quarkonia suppression was originally claimed to be an unambiguous signal of the formation of a quark-gluon plasma (QGP). Quarkonia suppression was suggested to be a signature of the QGP and we can measure the suppression ( as well as ), both at RHIC and at the LHC.

In heavy-ion collisions to determine the properties of the medium formed in A+A collisions and p + p collisions, the A+A collision deviates from simple superposition of independent p + p collisions. This deviation is quantified with the nuclear modification factor (). This factor is the ratio of the yield in heavy-ion collisions over the yield in p + p collisions, scaled by a model of the nuclear geometry of the collision. The value of =1 indicates no modification due to the medium. We can say that the probe of interest is suppressed in heavy-ion collisions if is less than 1. A quarkonia meson that forms on the outside surface will not dissociate regardless of the temperature of the medium because it does not have a chance to interact with it. This is why we never see a that is equal to zero. The suppression can also be affected by the QGP, the formation time of the quarkonia meson, and the QGP lifetime as well. For instance, a high quarkonia meson could have a formation time long enough that it actually does not see the QGP at all and thus is not suppressed.

In the early days most of the interests were focused on the suppression of charmonium states [13] of collider experiments at SPS and RHIC, but several observations are yet to be understood; namely, the suppression of (1S) does not increase from SPS to RHIC, even though the centre-of-mass energy is increased by fifteen times. The heavy-ion program at the LHC may resolve those puzzles because the beam energy and luminosity are increased by ten times that of the RHIC. Moreover the CMS detector has excellent capabilities for muon detection and provides measurements of (2S) and the family, which enables the quantitative analysis of quarkonia. That is why the interest may be shifted to the bottomonium states at the LHC energy.

A potential model for the phenomenological descriptions of heavy quarkonium suppression would be quite useful inspite of the progress of direct lattice QCD based determinations of the potential. The large mass of heavy quarks and their small relative velocity make the use of nonrelativistic quantum mechanics justifiable to describe the quarkonia in the potential models. This is one of the main goals of this present study that argues for the modification of the full Cornell potential as an appropriate potential for heavy quarkonium at finite temperature. QGP created at RHIC have a very low viscosity to entropy ratio, i.e., [49], and in the nonperturbative domain of QCD, with temperature close to , the quark matter in the QGP phase is strongly interacting.

In the present paper, we shall employ quasi-particle model for hot QCD equations of state [10, 11] to extract the Debye mass [12] which is obtained in terms of quasi-particle degrees of freedom. We first obtained the medium modified heavy quark potential in isotropic medium and estimate the dissociation temperature. Here, we have used the viscous hydrodynamics to define the dynamics of the system created in the heavy-ion collisions. We have included only the shear viscosity and not included the bulk viscosity. We will look the issue of bulk viscosity in near future.

Our work is organized as follows. In Section 2, we briefly discuss our recent work on medium modified potential in isotropic medium. In Sections 2.1 and 2.2 we study the real and imaginary part of the potential in the isotropic medium and effective fugacity quasi-particle model (EQPM) in Section 2.3. In Section 3 we studied binding energy and dissociation temperature of , , and state considering isotropic medium. Using this effective potential and by incorporating quasi-particle Debye mass, we have then developed the equation of state for strongly interacting matter and have shown our results on pressure, energy density, and speed of sound along with the lattice data. In Section 4, we have employed the aforesaid equation of state to study the suppression of bottomonium in the presence of viscous forces and estimate the survival probability in a longitudinally expanding QGP. Results and discussion will be presented in Section 5 and finally, we conclude in Section 6.

2. Medium Modified Effective Potential in Isotropic Medium

We can obtain the medium modification to the vacuum potential by correcting its both Coulombic and string part with a dielectric function encoding the effect of deconfinement [25]:

Here the functions, and , are the Fourier transform (FT) of the dielectric permittivity and Cornell potential, respectively. After assuming as distribution ( ) we evaluated the Fourier transform of the linear part asWhile putting , we can write the FT of the linear term asThus the FT of the full Cornell potential becomes

To obtain the real and imaginary parts of the potential, we put the temporal component of real and imaginary part in terms of retarded (or advanced) and symmetric parts in the Fourier space in isotropic medium which finally gives

Let us now discuss the real and imaginary part of the potential modified using the above define and along with effective fugacity quasi-particle model (EQPM) in the next subsections.

2.1. Real Part of the Potential in the Isotropic Medium

Now using the real part of retarded (advanced) propagator in isotropic medium, we getwhere the real part of the dielectric permittivity (also given in [2628]) becomesNow using (6) and real part of dielectric permittivity (7) in (1), we getSolving the above integral, we findwhere . In the limit , we have

2.2. Imaginary Part of the Potential in the Isotropic Medium

To obtain the imaginary part of the potential in the QGP medium, the temporal component of the symmetric propagator in the static limit has been considered, which reads [29, 30]Now the imaginary part of the dielectric function in the QGP medium isAfterwards, the imaginary part of the medium potential is easy to obtain owing to the definition of the potential (1) as done in [31]:After performing the integration, we findwhere .

2.3. Effective Fugacity Quasi-Particle Model (EQPM)

In our calculation, we use the Debye mass for full QCD:Here, is the QCD running coupling constant, () and is the number of flavors, the function has the form , and is the quasi-gluon effective fugacity and is quasi-quark effective fugacity. These distribution functions are isotropic in nature. These fugacities should not be confused with any conservations law (number conservation) and have merely been introduced to encode all the interaction effects at high temperature QCD. Both and have a very complicated temperature dependence and asymptotically reach to the ideal value unity [11]. The temperature dependence of and fits well to the form given below:Here and , , , and are fitting parameters, for both EOS1 and EOS2. Here, EoS1 is the hot QCD [1315] and EoS2 is the hot QCD EoS [16] in the quasi-particle description [10, 11], respectively. Now, the expressions for the Debye mass can be rewritten in terms of effective charges for the quasi-gluons and quarks aswhere and are the effective charges given by the equations:

In our present analysis we had used the temperature dependence of the quasi-particle Debye mass, , in full QCD with to determine charmonium suppression in an expanding, dissipative strongly interacting QGP medium. This quasi-particle Debye mass, , has the following form:

3. Binding Energy and Dissociation Temperature

To obtain the binding energies with heavy quark potential, we need to solve the Schrödinger equation numerically. In the limiting case discussed earlier, the medium modified potential resembles to the hydrogen atom problem [1]. The solution of the Schrödinger equation gives the eigenvalues for the ground states and the first excited states in charmonium (, , etc.) and bottomonium (, , etc.) spectra:where is the mass of the heavy quark.

In our analysis, we have fixed the critical temperature ( = ) and have taken the quark masses , as = 4.5 GeV, = 5.01 GeV, and GeV, as calculated in [32], and the string tension () is taken as . Let us now proceed to the computation of the dissociation temperatures for the above-mentioned quarkonia bound states.

As we know, dissociation of a quarkonia bound state in a thermal QGP medium will occur whenever the binding energy, , of the said state will fall below the mean thermal energy of a quasi-parton. In such situations, the thermal effect can dissociate the quakonia bound state. To obtain the lower bound of the dissociation temperatures of the various quarkonia states, the (relativistic) thermal energy of the partons will be 3 . The dissociation is supposed to occur whenever

’s for the sates , , and with the dissociation temperature are listed in Tables 1 and 2 for EoS1 and EoS2, respectively. We observe that (on the basis of temperature dependence of binding energy) dissociates at lower temperatures as compared to and for both the equations of state.

4. Formulation

In relativistic nucleus-nucleus collisions, the equation of state for the quark matter is an important observable and the properties of the matter are sensitive to it. The expansion of QGP is quite sensitive to EoS through the speed of sound and explores the sensitivity of the quarkonium suppression to the equation of state [33, 34].

For a strongly coupled QGP, Bannur [17] developed an equation of state by incorporating running coupling constant and did an appropriate modification to take account of color and flavor degrees of freedom and obtained a reasonably good fit to the lattice results. Now we will discuss briefly the equation of state which is expressed as a function of plasma parameter [35]:Plasma parameter is the ratio of average potential energy to average kinetic energy of particles, is assumed to be weak (), and is given byWe have studied the variation of plasma parameter with temperature and as well with the number of flavors that are present in the system and shown in Figure 1 for EoS1 and EoS2, respectively. As the temperature increases, potential becomes weaker and hence the plasma parameters have started waning; albeit at very large temperature it increases slightly due to the contribution coming from the (positive) finite-range terms in the potential, unlike the decreasing trend in Bannur model [17] always due to the presence of Coulomb interaction alone in the deconfined phase.

Let us consider that hadron exists for and goes to QGP for for strongly coupled plasma in QCD. As it was assumed that confinement interactions due to QCD vacuum have been melted [17] at and thus for , there are the strongly interacting plasma of quarks and gluons and no glue balls or hadrons. After inclusion of relativistic and quantum effects, the equation of state which has been obtained in the plasma parameter can be written asNow, the scaled-energy density is written as in terms of ideal contributionwhere is given byHere, is the number of flavors of quarks and gluons. Now, we will employ two-loop level QCD running coupling constant in scheme [36]:Here and . In scheme, and are the renormalization scale and the scale parameter, respectively. For the EoS to depend on the renormalization scale, the physical observables should be scale independent. We invade the problem by trading off the dependence on renormalization scale () to a dependence on the critical temperature . where =0.5772156 and , which is a constant depending on colors and flavors. There are several incertitude, associated with the scale parameter and renormalization scale , which occurs in the expression used for the running coupling constant . This issue has been considered well in literature and resolved by the BLM criterion due to Brodsky, Lepage, and Mackenzie [37]. is allowed to vary between and [38]. For our motive, we choose close to the central value [39] for =0 and for both =2 and =3 flavors the value is . If the factor is , then the above expression reduces to the expression used in [17, Eq. ], after neglecting the higher order terms of the above factor. However, this possibility does not hold good for the temperature ranges used in the calculation and causes an error in coupling which finally makes the difference in the results between our model and Bannur model [17]. First of all, we will calculate the energy density from (25) and using the thermodynamic relationwe calculated the pressure aswhere is the pressure at some reference temperature . Now, the speed of sound can be calculated once we know the pressure and energy density .

5. Survival of Bottomonium State

In order to derive the survival probability for an expanding QGP firstly, we explore the effects of dissipative terms up to first order in the stress-tensor. In the presence of viscous forces, the energy-momentum tensor is written aswhere the stress-energy tensor, , up to first order is given bywhere is the coefficient of the shear viscosity and is the symmetrized velocity gradient.

In Bjorken expansion, the equation of motion is given byThe solution of equation of motion (33) is given as

where the constant isand the symbols areandThe first term accounts for the contributions coming from the zeroth-order expansion (ideal fluid) and the second term is the first-order viscous corrections. We now have all the ingredients to write down the survival probability. Chu and Matsui [40] studied the transverse momentum dependence () of the survival probability by choosing the speed of sound (ideal EoS) and the extreme value . Instead of taking arbitrary values of , we tabulated the values of in Tables 1 and 2 corresponding to the dissociation temperatures for bottomonium states for EOS1 and EOS2. One can define initial energy density asHere, represents the proportionality of the deposited energy to the nuclear thickness where is the average initial energy density and will be given by the modified Bjorken formula [41, 42]:where is the transverse overlap area of the colliding nuclei and is the transverse energy deposited per unit rapidity. We use the experimental value of the transverse overlap area and the pseudo-rapidity distribution [43, 44] at various values of number of participants . These numbers are then multiplied by a Jacobian 1.25 to yield the rapidity distribution which will be further used to calculate the average initial energy density from Bjorken formula (39). After getting the value of average initial energy density we can obtain the initial energy density from formula (38). The scaling factor has been introduced in order to obtain the desired values of initial energy densities [45, 46] for most central collision which are consistent with the predictions of the self-screened parton cascade model [47] and also with the requirements of hydrodynamic simulation [45, 46] to fit the pseudo-rapidity distribution of charged particle multiplicity for various centralities observed in PHENIX experiments at RHIC energy. Let be the angle between the transverse momentum and position vector . Now assuming that is formed inside screening region at a point whose position vector is and moves with transverse momentum making an azimuthal angle, then the condition for escape of without forming bottomonium states is expressed aswhere is the position vector at which the bottom, anti-bottom quark pair is formed, is the proper formation time required for the formation of bound states of from correlated pair, and is the mass of bottomonia (, , for different resonance states of bottomonium). Assume the radial probability distribution for the production of pair in hard collisions at transverse distance as Here we take in our calculation as used in [40]. Then, in the color screening scenario, the survival probability for the bottomonium in QGP medium can be expressed as [40, 48, 49]where the maximum positive angle allowed by (26) becomes [50]since the experimentalists always measure the quantity, namely, integrated nuclear modification factor. We get the theoretical integrated survival probability as follows:

In nucleus-nucleus collisions, it is known that only about 60% of the observed originate directly in hard collisions while 30% of them come from the decay of and 10% from the decay of . Hence, the -integrated inclusive survival probability of in the QGP becomes [33, 51]

6. Results and Discussions

In our results, we had obtained the variation of plasma parameter with temperature and as well with the number of flavors that are present in the system and shown in Figure 1 for EoS1 and EoS2, respectively. After that, in Figure 2, we have plotted the variation of pressure () with temperature () using EoS1 and EoS2 for 3-flavor QGP along with Bannur EoS [17] and compared it with lattice results [1721]. For each flavor, and are adjusted to get a good fit to lattice results in Bannur model. Now, energy density , speed of sound , and so forth can be derived since we had obtained the pressure, . In Figure 3, we had plotted the energy density () with temperature () using EoS1 [1315] and EoS2 for 3-flavor QGP along with Bannur EoS [17] and compared it with lattice result [1721]. In Figure 4, the speed of sound, , is plotted using EoS1 and EoS2 for 3-flavor QGP along with Bannur EoS [17]. Since lattice results are not available for 3 flavors, therefore comparison has not been checked for the above-mentioned flavor. Our flavored results match excellent with the lattice results.

In this paper, we had calculated the dissociation temperatures for the bottomonium states (, , , etc.), by modifying the Cornell potential and incorporating the quasi-particle Debye mass. On that dissociation temperature, we had calculated the screening energy densities, , and the speed of sound which are also listed in Tables 1 and 2 for both EoS1 and EoS2, respectively. We observe from Tables 1 and 2 that the value of is different for different bottomonium states and varies from one EoS to other. If , initial energy density, then there will be no suppression at all, i.e., survival probability, , is equal to 1. With this physical understanding, we analyze our results, , as a function of the number of participants in an expanding QGP.

Here we are using the values as inputs listed in Tables 1 and 2, to calculate for both EOS1 and EOS2, respectively. The experimental data (the nuclear modification factor ) are shown by the squares with error bars whereas circles represent sequential suppression. We had compared our results with the experimental results for the case of for both EoS1 and EoS2 and found good agreement. We observe from Figures 510 that for both the directly and sequentially produced Upsilon () are quite high with the higher values of ’s which is obtained from EOS2 (in Table 2) compared to EOS1 (in Table 1) for both SIQGP and ideal equation of states. We find that the survival probability of sequentially produced is slightly higher compared to the directly produced and is closer to the experimental results. We also observed that sequentially produced nicely matches for the EOS1 compared to the EOS2. The smaller value of screening energy density causes an increase in the screening time and results in more suppression to match with the experimental results.

7. Conclusions

We studied the equation of state for strongly interacting quark-gluon plasma in the framework of strongly coupled plasma with appropriate modifications to take account of color and flavor degrees of freedom and QCD running coupling constant. In addition, we incorporate the nonperturbative effects in terms of nonzero string tension in the deconfined phase, unlike the Coulomb interactions alone in the deconfined phase beyond the critical temperature. Our results on thermodynamic observables, namely, pressure, energy density, and speed of sound, nicely fit the results of lattice equation of state. We had then calculated the dissociation temperatures for the bottomonium states (, , , etc.), by incorporating the quasi-particle Debye mass. On that dissociation temperature, we had calculated the screening energy densities, , and the speed of sound which are listed in Tables 1 and 2 for both EoS1 and EoS2, respectively. By using the above quantities as an input, we have then studied the sequential suppression for bottomonium states at the LHC energy in a longitudinally expanding partonic system, which underwent through the successive preequilibrium and equilibrium phases in the presence of dissipative forces. Bottomonium suppression in nucleus-nucleus collisions compared to - collisions couples the in-medium properties of the bottomonia states with the dynamics of the expanding medium. We have found a good agreement with the experimental data from RHIC 200GeV/nucleon Au-Au collisions, LHC 2.76 TeV/nucleon Pb-Pb, and LHC 5.02 TeV/nucleon Pb-Pb collisions [52, 53]. Here our attempt is to understand suppression systematically in SIQGP in anisotropic medium. It would be of interest to extend the present study by incorporating the contributions of the bulk viscosity. These issues will be taken up separately in the near future.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Vineet Kumar Agotiya acknowledges the UGC-BSR research start up Grant no. F.30-14/2014 (BSR), New Delhi. The authors record their sincere gratitude to the people of India for their generous support for the research in basic sciences.