<p>The paper focuses on the effects of temperature, concentration, and velocity slip on the steady, laminar boundary layer flow and heat transfer dynamics of Sisko nanofluids over a nonlinearly stretching sheet. By applying appropriate similarity variables and transformations, the governing partial differential equations, formulated using the Buongiorno model, are reduced to a system of nonlinear ordinary differential equations, which are then solved numerically using the shooting method. The study analyses the effects of parameters such as the stretching parameter, thermophoresis number, concentration slip parameter, Brownian motion number, Lewis number, and thermal slip parameter on temperature, velocity, and nanoparticle concentration profiles. Additionally, the impact on mass and heat transfer rates is assessed. The results show that increasing the velocity slip parameter decreases the fluid velocity but increases both temperature and concentration near the surface. Shear-thinning fluids exhibit thicker boundary layers, while shear-thickening fluids result in sharper gradients and thinner layers. Thermal and concentration relaxation times reduce the intensity of heat and mass transfer within the boundary layer. Excellent agreement is found when comparing the current findings with numerical solutions that are already in use and results that have been published. These findings provide a significant understanding of the behavior of Sisko nanofluids over nonlinearly stretched surfaces and emphasize the crucial role that slip effects play in maximizing heat transfer efficiency. CFD simulations have played a vital role in investigating nanofluid behavior, optimizing system configurations, and improving heat transfer performance.</p>