“…To investigate nonlinear interplay instabilities in the laser-plasma interaction it is desirable to obtain differential equations for the evolution of the macroscopic quantities that characterize the laser beam profile [25,26], including the amplitude, spot size, phase, curvature radius, and centroid. There are various methods for attempting to obtain such envelope equations, including the variation method [27], the moment method [28], and the source dependent expansion (SDE) technique [29]. We can use these mentioned techniques to study the nonlinear relativistic self-focusing phenomena with taking into account the relativistic mass corrections.…”