We build up a model to reproduce the experimentally measured $R_{\mathit{AA}}$ of $J/\psi$ near midrapidty in an $\mathrm{Au}+\mathrm{Au}$ collision at $\sqrt{s_{\mathit{NN}}}=200$ GeV. The model takes into account the $J/\psi$ suppression from the quark-gluon plasma and hadron gas as well as the nuclear absorption of primordial charmonia and the regeneration effects at the hadronization stage and hence is a generalization of the two-component model introduced by Grandchamp and Rapp. The improvements in this work are twofold; the addition of the initial local temperature profile and a consistent use of QCD next-to-leading order (NLO) formula for both the dissociation cross section in the hadron gas and the thermal decay widths in the quark-gluon plasma for the charmonium states. The initial local temperature profile is determined from the assumption that the local entropy density is proportional to a formula involving the number densities of the number of participants and of the binary collisions that reproduces the multiplicities of charged particles at chemical freeze-out. The initial local temperature profile brings about a kink in the $R_{\mathit{AA}}$ curve due to the initial melting of $J/\psi$. The initially formed fireball, composed of weakly interacting quarks and gluons with thermal masses that are extracted from lattice QCD, follows an isentropic expansion with cylindrical symmetry. The fit reproduces well the $\mathrm{Au}+\mathrm{Au}$ as well as the $\mathrm{Cu}+\mathrm{Cu}$ data. The same method is applied to predict the $R_{\mathit{AA}}$ expected from the $\mathrm{Pb}+\mathrm{Pb}$ collision at Large Hadron Collider (LHC) energy.

• Received 10 February 2010

DOI:https://doi.org/10.1103/PhysRevC.81.034914