This study investigates fluid flow over a rotating disk with a focus on the impact of radiative heat flow of MHD elastico-viscous nanofluid. The analysis accounts for viscous dissipation, variable thermal conductivity, heat source, and chemical reaction. The governing nonlinear coupled partial differential equations were formulated and transformed into a system of nonlinear ordinary differential equations using appropriate similarity variables. The resulting nondimensional models were subjected to the spectral local linearization method (SLLM) for analysis. The numerical solution was validated against existing literature, demonstrating good agreement. The effects of relevant parameters were presented through graphs and tables. It was observed that the elastico-viscous fluid and magnetic parameters hinder fluid profiles, while the opposite effect is observed in temperature and concentration profiles, respectively. Moreover, the thermal boundary layer thickness increases with higher Eckert number, radiation parameter, thermophoresis parameter, Brownian motion parameter, and heat source parameter values.