The present article investigates the effect of multiple slips on axisymmetric magnetohydrodynamics (MHD) buoyant nano-fluid flow over a stretching sheet with radiation and chemical effect. The non-linear partial differential equations were transformed to a non-linear control equation using an appropriate similarity transformation. The governing equations were solved through the finite element method. The influence of physical parameters such as multiple slips, magnetic, thermal radiation, Prandtl number, stretching, Brownian motion, thermophoresis, Schmidt number, Lewis number and chemical reaction on the radial velocity, temperature, solutal concentration and nano-fluid volume fraction profile were investigated. We noted that the boundary layers increases in the presence of multiple slip effects whereas, the effect of thermal slip on Nusselt number increases with the increasing values of magnetic and thermal radiation. To verify the convergence of the numerical solution, the computations were made by reducing the mesh size. Finally, our results are parallel to previous scholarly contributions.