Symmetry Analysis and Conservation Laws of the Zoomeron Equation

In this work, we study the (2 + 1)-dimensional Zoomeron equation which is an extension of the famous (1 + 1)-dimensional Zoomeron equation that has been studied extensively in the literature. Using classical Lie point symmetries admitted by the equation, for the first time we develop an optimal system of one-dimensional subalgebras. Based on this optimal system, we obtain symmetry reductions and new group-invariant solutions. Again for the first time, we construct the conservation laws of the underlying equation using the multiplier method.


Introduction
Many physical phenomena of the real world are governed by nonlinear partial differential equations (NLPDEs).It is therefore absolutely necessary to analyse these equations from the point of view of their integrability and finding exact closed form solutions.Although this is not an easy task, many researchers have developed various methods to find exact solutions of NLPDEs.These methods include the sine-cosine method [1], the extended tanh method [2], the inverse scattering transform method [3], the Hirota's bilinear method [4], the multiple exp-function method [5], the simplest equation method [6,7], non-classical method [8], method of generalized conditional symmetries [9], and the Lie symmetry method [10,11].
This paper aims to study one NLPDE; namely, the (2 + 1)-dimensional Zoomeron equation [12] u xy u tt which has attracted some attention in recent years.Many authors have found closed-form solutions of this equation.For example, the (G /G)−expansion method [12,13], the extended tanh method [14], the tanh-coth method [15], the sine-cosine method [16,17], and the modified simple equation method [18] have been used to find closed-form solutions of (1).The (2 + 1)-dimensional Zoomeron equation with power-law nonlinearity was studied in [19] from a Lie point symmetries point of view and symmetry reductions, and some solutions were obtained.Additionally, in [19], the authors have given a brief history of the (1 + 1)-dimensional Zoomeron equation.See also [20][21][22].
In this paper we first use the classical Lie point symmetries admitted by Equation ( 1) to find an optimal system of one-dimensional subalgebras.These are then used to perform symmetry reductions and determine new group-invariant solutions of (1).It should be noted that such approach was previously used for examination of a wide range of nonlinear PDEs [23][24][25][26][27][28][29][30][31].Furthermore, we derive the conservation laws of (1) using the multiplier method [32,33].
The paper is organized as follows: in Section 2, we compute the Lie point symmetries of (1) and use them to construct the optimal system of one-dimensional subalgebras.These are then used to perform symmetry reductions and determine new group-invariant solutions of (1).In Section 3, we derive conservation laws of (1) by employing the multiplier method.Finally, concluding remarks are presented in Section 4.

Symmetry Reductions and Exact Solutions of (1) Based on Optimal System
In this section, firstly we use the Lie point symmetries admitted by (1) to construct an optimal system of one-dimensional subalgebras.Thereafter, we obtain symmetry reductions and group-invariant solutions based on the optimal system of one-dimensional subalgebras [23,24].

Lie Point Symmetries of (1)
The Lie point symmetries of the Zoomeron Equation (1) are given by [19] which generate a five-dimensional Lie algebra L 5 .

Optimal System of One-Dimensional Subalgebras
In this subsection, we use the Lie point symmetries of (1) to compute an optimal system of one-dimensional subalgebras.We employ the method given in [23,24], which takes a general element from the Lie algebra and reduces it to its simplest equivalent form by using the chosen adjoint transformations where ε is a real number, and [X i , X j ] denotes the commutator defined by The table of commutators of the Lie point symmetries of Equation ( 1) and the adjoint representations of the symmetry group of (1) on its Lie algebra are given in Tables 1 and 2, respectively.Then, Tables 1 and 2 are used to construct the optimal system of one-dimensional subalgebras for Equation (1).
Using Tables 1 and 2, we can construct an optimal system of one-dimensional subalgebras, which is given by {X 3 , X 4 , X 5 ,

Symmetry Reductions
In this subsection, we use the optimal system of one-dimensional subalgebras computed in the previous subsection, and present symmetry reductions of (1) to two-dimensional partial differential equations.
For the first operator X 3 of the optimal system, we have the three invariants s = t, r = x, f = u, and using these invariants, (1) reduces to Likewise for X 4 , the invariants s = ty, r = xy, f = u transforms (1) to Similarly, the invariants s = t, r = y − x, f = u of X 2 + X 3 reduces (1) to The symmetry X 1 + X 5 has invariants s = x, r = ye −2t , f = uy 1/2 , and these reduce (1) to Using the invariants s = x/t, r = y/t, f = tu of X 4 + X 5 , (1) reduces to The operator X 1 + X 2 + X 3 has invariants s = x − t, r = y − t, f = u, and with the use of these invariants, (1) reduces to Finally, X 1 + X 2 + X 5 has invariants s = x − t, r = ye −2t , f = ue t , and their use reduces (1) to

Group-Invariant Solutions
We now obtain group-invariant solutions based on the optimal system of one-dimensional subalgebras.However, in this paper we are looking only at some interesting cases.

Case 1. X 5 = 2y∂/∂y − u∂/∂u
The associated Lagrange system to the operator X 4 yields three invariants s = t, r = x, u = y −1/2 U(r, s), which give group-invariant solution u = y −1/2 U(s, r) and transforms (1) to This equation has three Lie point symmetries, viz., The symmetry Γ 1 − νΓ 2 gives the two invariants z = r + νs and F = U.Using these invariants, (2) transforms to the nonlinear third-order ordinary differential equation Integrating (3) twice with respect to z, we obtain where k 1 and k 2 are constants of integration.The solutions of this equation are given by where k 3 is a constant of integration and erfi(z) is the imaginary error function [34].Thus, solutions of (1) are The associated Lagrange system to this operator yields the three invariants ( The Lie point symmetries of the above equation are The symmetry Γ 2 gives the two invariants z = s and U = r −1/2 F, and using these invariants, (2) transforms to the nonlinear third-order ordinary differential equation Integrating (6) twice with respect to z, we obtain where k 1 and k 2 are constants of integration.The solution of this equation is given by where k 3 is a constant of integration.Thus, a solution of ( 1) is which is a steady-state solution.
The associated Lagrange system to this symmetry operator gives three invariants, viz., which give group-invariant solution u = U(s, r) and reduces (1) to The Lie point symmetries of the above equation are The symmetry Γ 1 − νΓ 2 gives the two invariants z = r + νs and F = U.Using these invariants, (8) transforms to the nonlinear fourth-order ordinary differential equation Integrating ( 9) twice with respect to z, we obtain where k 1 and k 2 are constants of integration.This equation can not be integrated in the closed form.However, by taking k 1 = 0, one can obtain its solution in the closed form in the following manner.Multiplying (10) with k 1 = 0 by F and integrating, we obtain where k 3 is a constant of integration.The solution of this equation is given by 0 and sn is the Jacobi elliptic sine function [35].Thus, a solution of (1) is For k 3 = 0 we have the solution given by and when C = 0 we have Likewise, one may obtain more group-invariant solutions using the other symmetry operators of the optimal system of one-dimensional subalgebras.For example, the symmetry operator X 2 + X 3 of the optimal system gives us the group-invariant solution (2.9) of [19] in terms of the Airy functions.

Conservation Laws of (1)
Conservation laws describe physical conserved quantities, such as mass, energy, momentum and angular momentum, electric charge, and other constants of motion [32].They are very important in the study of differential equations.Conservation laws can be used in investigating the existence, uniqueness, and stability of the solutions of nonlinear partial differential equations.They have also been used in the development of numerical methods and in obtaining exact solutions for some partial differential equations.
A local conservation law for the (2 + 1)-dimensional Zoomeron Equation ( 1) is a continuity equation holding for all solutions of Equation ( 1), where the conserved density T and the spatial fluxes X and Y are functions of t, x, y, u.The results in [11] show that all non-trivial conservation laws arise from multipliers.Specifically, when we move off of the set of solutions of Equation ( 1), every non-trivial local conservation law ( 12) is equivalent to one that can be expressed in the characteristic form holding off of the set of solutions of Equation (1) where Q(x, y, t, u . ..) is the multiplier, and where ( T, X, Ỹ) differs from (T, X, Y) by a trivial conserved current.On the set of solutions u(x, y, t) of Equation ( 1), the characteristic form (13) reduces to the conservation law (12).In general, a function Q(x, t, u . ..) is a multiplier if it is non-singular on the set of solutions u(x, y, t) of Equation ( 1), and if its product with Equation ( 1) is a divergence expression with respect to t, x, y.There is a one-to-one correspondence between non-trivial multipliers and non-trivial conservation laws in characteristic form.
The determining equation to obtain all multipliers is δ δu where δ/δu is the Euler-Lagrange operator given by Equation ( 14) must hold off of the set of solutions of Equation (1).Once the multipliers are found, the corresponding non-trivial conservation laws are obtained by integrating the characteristic Equation ( 13) [11].
We will now find all multipliers Q(x, y, t, u) and obtain corresponding non-trivial (new) conservation laws.
The determining Equation ( 14) splits with respect to the variables u t , u x , u y , u tt , u tx , u ty , u xy , u yy , u ttt , u ttx , u txy , u tyy , u xyy , u tttx , u txyy .This yields a linear determining system for Q(x, y, t, u) which can be solved by the same algorithmic method used to solve the determining equation for infinitesimal symmetries.By applying this method, for Equation (1), we obtain the following linear determining equations for the multipliers: Q ty (t, x, y, u) = 0, ( Q yyy (t, x, y, u) = 0, ( Q tt (t, x, y, u) − Q yy (t, x, y, u) = 0.
It is straightforward using Maple to set up and solve this determining system ( 15)-( 18), and we get the four multipliers given by For each solution Q, a corresponding conserved density and flux can be derived (up to local equivalence) by integration of the divergence identity (13) [11,36].We obtain the following results.
, which give group-invariant solution u = e −t U(s, r) and transforms (1) to U U