This analysis intends to address the coupled effect of phase change heat transfer, thermal radiation, and viscous heating on the MHD flow of an incompressible chemically reactive nanofluid in the vicinity of the stagnation point toward the stretching surface, taking a Jeffrey fluid as the base fluid. Convergent analytical solutions for the nonlinear boundary layer equations are obtained by the successive application of scaling variables and the highly efficacious homotopy analysis method. Error analysis is implemented to endorse the convergence of the solutions. Through parametric examination, influence of various physical parameters occurring in analysis of the profiles of velocity, temperature, and nanoparticle concentration, coefficient of surface drag, rates of mass and heat transfer is explored pictorially. The Deborah number and the melting parameter are found to enhance velocity, and the associated momentum boundary layers are thicker, whereas the magnetic field depreciates the flow rate. Temperature is observed to enhance with the thermophoresis parameter, Prandtl number and Eckert number, whereas a reduction is seen with the thermal radiation parameter and Brownian motion parameter. Nanoparticle concentration is depleted by the chemical reaction parameter, the thermophoresis parameter, and the Lewis number.