Evolution of heavy quark distribution function on quark-gluon plasma: using the iterative Laplace transform Method

The Iterative Laplace Transform Method is used to solve the Fokker-Plank equation for finding the time evolution of the heavy quarks distribution functions such as charm and bottom in quark gluon plasma. The Iterative Laplace Transform Method is a combination of two powerful methods, namely, the Laplace transform method and the Iterative method. This method gives numerical solutions in the form of convergent series with easily computable components. The most outstanding feature of this method is that it provides an analytical solution using the initial conditions only, without any discretization or restrictive assumptions. To verify our solutions, we calculate the nuclear suppression factor, R_AA(p_T), of heavy quarks using simulated data of heavy ion collisions. The results are in good agreement with available experiment data from the PHENIX collaboration