Analysis has been conducted to analyze the effects of second order slip flow and heat transfer of Jeffrey nanofluid over a stretching sheet with non linear thermal radiation and chemical reaction. The effects of Brownian motion and thermophoresis occur in the transport equations. The velocity, temperature and nanoparticle concentration profiles are analyzed with respect to the involved parameters of interest namely Brownian motion parameters, thermophoresis parameter, magnetic parameter, radiation parameter, Prandtl number, Lewis number, chemical reaction parameter, and Deborah number, Convergence of the derived solutions was checked and the influence of embedded parameters was analyzed by plotting graphs. It was noticed that the velocity increases with an increase in the Deborah number. We further found that for fixed values of other parameters, numerical values of the skin friction coefficient, local Nusselt numbers and Sherwood numbers were computed and examined. A comparative study between the previous published and present results in a limiting sense is found in an excellent agreement.