This paper investigates the non-linear coaxial (or coupled mode) propagation of Laguerre–Gaussian (LG) (in particular L01 mode) and Gaussian electromagnetic (em) beams in a homogeneous plasma characterized by ponderomotive and relativistic non-linearities. The formulation is based on numerical solution of non-linear Schrödinger wave equation under Jeffreys–Wentzel–Kramers–Brillouin approximation, followed by paraxial approach applicable in the vicinity of intensity maximum of the beams. A set of coupled differential equations for spot size (beam width) and phase evolution with space corresponding to coupled mode has been derived and numerically solved to determine the propagation dynamics. Using focusing equation a critical condition describing the self-trapped (i.e., spatial soliton) mode of laser beam propagation in the plasma has been discussed; as a consequence oscillatory focusing/defocusing of the beams in coupled mode propagation have been analyzed and presented graphically. As an important outcome, significant enhancement in the intensity of LG beam is noticed when it is coupled with the Gaussian mode.