Prior hybrid nanofluid research is restricted to single or two-phase models, which were dominated by homogeneous nanofluid models. These models adhere a serious flaw with the heat-mass transfer indices. To handle this problem, we develop an entirely novel hybrid nanofluid models based on the fundamental slip mechanisms of nanofluids induced by Brownian and thermophoretic diffusion as well as the chemical reaction of binary species. The laminar flow of MHD rheological model with suspended nanoparticles is striking at the stagnation region laden over a stretching surface. In order to reach the real world scenario, the expression of radiative heat, non-uniform heat source and frictional heating are all included. The coupled Buongiorno's model and Tiwari-Das mixtures models incorporate the effect of particles diffusion and volumetric friction, respectively. It is hypothesized that first order kinetics can describe homogeneous events in the immediate proximity, while isothermal cubic autocatalysis kinetics frequently describe heterogeneous reactions occur on the sheet surface. The restructured mathematical equations are solved with a proficient Runge Kutta Fehlberg technique together with the shooting procedure. According to the outcomes, the presence of higher magnetic field diminished the flow and velocity gradient. The effect of elastic parameters is diminishing frictional coefficient. The effect of first and second Brownian (Nb1, Nb2) and thermophoretic parameters (Nt1, Nt2) are similar for temperature. The mass-heat transfer rate is significant with higher estimations of diffusivity ratio Nbt. The momentum rate of diffusion outweighs the species diffusion rate as the Schmidt number increases, resulting in a decrease in concentration values. Higher intensities of Lorentz's force results in denser streamlines. The streamlines split into an opposite direction because of stagnation point.