“…[2,9] Previously, the propagation dynamics of beams were well discussed through numerical simulation [10] or inverse scattering modeling and so on. [11] In recent years, people often use Jacobi elliptic function, [12] variable separation method, [13] variational method (VM), [14][15][16] and other integration algorithms to solve relevant problems. Moreover, some properties of beams in propagation process have been verified in many aspects, such as nonlocalization in dispersion shock wave, [17] Anderson localization, [18] gradient mutation, guide wave distribution coupling, and photon condensation.…”