A mathematical model is envisioned to discourse the impact of Thompson and Troian slip boundary in the carbon nanotubes suspended nanofluid flow near a stagnation point along an expanding/contracting surface. The water is considered as a base fluid and both types of carbon nanotubes i.e., single-wall (SWCNTs) and multi-wall (MWCNTs) are considered. The flow is taken in a Dacry-Forchheimer porous media amalgamated with quartic autocatalysis chemical reaction. Additional impacts added to the novelty of the mathematical model are the heat generation/absorption and buoyancy effect. The dimensionless variables led the envisaged mathematical model to a physical problem. The numerical solution is then found by engaging MATLAB built-in bvp4c function for non-dimensional velocity, temperature, and homogeneous-heterogeneous reactions. The validation of the proposed mathematical model is ascertained by comparing it with a published article in limiting case. An excellent consensus is accomplished in this regard. The behavior of numerous dimensionless flow variables including solid volume fraction, inertia coefficient, velocity ratio parameter, porosity parameter, slip velocity parameter, magnetic parameter, Schmidt number, and strength of homogeneous/heterogeneous reaction parameters are portrayed via graphical illustrations. Computational iterations for surface drag force are tabulated to analyze the impacts at the stretched surface. It is witnessed that the slip velocity parameter enhances the fluid stream velocity and diminishes the surface drag force. Furthermore, the concentration of the nanofluid flow is augmented for higher estimates of quartic autocatalysis chemical.