A significance of Marangoni-driven convective magnetohydrodynamics (MHD) was investigated in a planned study. The present study examines the impact of radiation and chemical reactions on the two dimensional boundary layer flow of bi-viscous Bingham fluid on a thermosolutal Marangoni boundary with magnetic field and heat source/sink. The physical flow problem is mathematically modelled into Navier–Stokes equations. These nonlinear partial differential equations are then mapped into a set of nonlinear ordinary differential equations using similarity transformation. The elegant analytical method is used to obtain analytic approximations for the resulting system of nonlinear differential equations. The analytical solution for the temperature and concentration profiles is expressed in terms of Kummer's confluent hypergeometric function. The features of flow characteristics such as velocity, temperature, and concentration profiles in response to the variations of the emerging parameters are simulated and examined with a physical explanation through graphs. Surface tension is considered to be corresponding to heat and mass. The results reveal that the velocity profile decreases with increasing the magnetic field, while the Marangoni number value increases with increasing the velocity profile. Moreover, the temperature profile rises with rising the radiation parameter and whereas the value chemical reaction parameter upsurges with decline concentration profile. The present problem has great interest due to its applications in industrial applications such as drying of silicon wafers, thin layers of paint, glues, in heat exchangers, crystal growth in space, biological fluids, blood, synovial liquids, plasma, lubricants, many pharmaceutical products, proteins, sewage sludge, china clay, and many more.