“…Non-linearity theories depend heavily on analytical solutions to nonlinear partial differential equations (NPDEs) which can be used to interpret natural phenomena in physical and applied studies such as fluid mechanics, hydrodynamics, mathematical physics, optics, elasticated media, chemical reactions, astrophysics, ecosystems, quantum theory, geology, plasma physics, wave propagation, and shallow water [1][2][3][4][5]. Several methods for obtaining solutions of nonlinear partial differential equations have been improved recently, for instance, the method of the inverse scattering [6], the Bernoulli sub-equation function methods [7,8], and the use of the Laplace transformation for the system which involves the Caputo fractional derivatives [9]. To investigate the ion-acoustic wave constructions in plasma physics, the (3+1)dimensional gKdV-ZK equation is used which is taken in the following configuration [10,11],…”