The flow and heat transfer of non-Newtonian nanofluids has an extensive range of applications in oceanography, the cooling of metallic plates, melt-spinning, the movement of biological fluids, heat exchangers technology, coating and suspensions. In view of these applications, we studied the steady Marangoni driven boundary layer flow, heat and mass transfer characteristics of a nanofluid. A non-Newtonian second-grade liquid model is used to deliberate the effect of activation energy on the chemically reactive non-Newtonian nanofluid. By applying suitable similarity transformations, the system of governing equations is transformed into a set of ordinary differential equations. These reduced equations are tackled numerically using the Runge–Kutta–Fehlberg fourth-fifth order (RKF-45) method. The velocity, concentration, thermal fields and rate of heat transfer are explored for the embedded non-dimensional parameters graphically. Our results revealed that the escalating values of the Marangoni number improve the velocity gradient and reduce the heat transfer. As the values of the porosity parameter increase, the velocity gradient is reduced and the heat transfer is improved. Finally, the Nusselt number is found to decline as the porosity parameter increases.