This article investigates the combined effect of second-order velocity slip, Arrhenius activation energy and binary chemical reaction on reactive Casson nanofluid flow in a non-Darcian porous medium. The governing equations of the problem were first non-dimensionalized and later reduced to ordinary nonlinear differential equations by adopting a similarity transformation. The emerging nonlinear boundary value problem was solved by using Galerkin weighted residual method (GWRM). The obtained results were compared with those found in the literature to validate our method. The impact of pertinent parameters on the velocity component, temperature distribution and concentration profile are presented using graphs and were discussed. The computational results show that an increase in second order slip parameter significantly results to an increase in the temperature as well as nanoparticle concentration profiles, while it reduces the velocity profile.