Non-Newtonian activity has been found in the manufacturing of coated sheets, foods, fibre optics, drilling muds, and plastics polymers. The Jeffery solvent is one of these fluid models that has really piqued the curiosity of scientists. With such a modelling approach, relaxation and delayed responses are possible. The investigation of boundary layer motions across a stretched cylinder has lots of applications in filament activities, crystallization and sheet fabrication, photonic circuits, recrystallization, and other industries. The convective heat transfer motion of a Jeffery non-Newtonian liquid influenced by a stretchy material with thermal performance is investigated in this paper. Applying transformation variables, the Physical quantities are converted to physical quantitates of the first order. The Runge–Kutta fourth order algorithm, in conjunction with the shooting strategy, is used to provide numerical results for motion and concentration profiles. The impacts of numerous factors on speed and concentration domains, such as the Hartman number, Schmidt number, and chemical reaction parameter, are depicted in graphs. The influence of Deborah numbers and Hartmann numbers on velocities components is the same. As the magnitude of the Hartmann number raised, boundary-layer thickness decreases. The concentration fields shrink as Sc increases and fall for the destructive reaction factor. The surface mass transfer decreases as Hartmann’s number is increased. The heat-transfer rate statistics data are very thoroughly recorded and analysed. The novelty of the present work is to investigate the Jeffery fluid in the presence of a magnetic field along with a chemical reaction over a nonlinear stretching sheet, which has not been investigated yet, so for limiting cases the present work is compared with the published work and good agreement is found.