Third-Order Conditional Lie – Bäcklund Symmetries of Nonlinear Reaction-Diffusion Equations

The third-order conditional Lie–Bäcklund symmetries of nonlinear reaction-diffusion equations are constructed due to themethod of linear determining equations. As a consequence, the exact solutions of the resulting equations are derived due to the compatibility of the governing equations and the admitted differential constraints, which are resting on the characteristic of the admitted conditional Lie–Bäcklund symmetries to be zero.


Introduction
The conditional Lie-Bäcklund symmetry (CLBS) was independently proposed by Zhdanov [1] and Fokas and Liu [2].Exact solutions possessing the form of nonlinear separation of variables [3] and exact solutions defined on the invariant subspaces [4] for nonlinear heat conductivity equations were constructed due to the corresponding symmetry reductions in [1], which cannot be obtained within the framework of Lie's classical symmetry [5], conditional symmetry [6,7], and Lie-Bäcklund symmetry [8].Multishock and multisoliton solutions of nonintegrable equations were obtained by using CLBS in [2], which are regarded as remnants of integrability for the considered nonintegrable equations.
As a consequence, CLBSs were applied to a number of particular subclasses of evolution equations [9,10].It has been noted that CLBS is very effective in studying classifications and reductions of diffusion equations [11][12][13][14][15][16].CLBS can also be used to consider the reductions of initial value problems for nonlinear evolution equations [17][18][19].CLBS for evolution systems were studied in [20,21].
The procedure for determining whether or not a given Lie-Bäcklund symmetry is conditionally invariant of the considered equations is straightforward; however, the determination of the most general CLBS admitted by a given differential equation is a very difficult, if not impossible, problem.CLBSs are generalizations of conditional symmetries, in the same way that Lie-Bäcklund symmetries are generalizations of Lie's classical symmetry.Therefore the procedure for computing CLBSs is about the same as that for the conditional symmetries.We need to determine the form of CLBS presumably.It is proved that CLBSs related to separation of variables [13], sign-invariants [14], and invariant subspaces [15,16] are powerful to study different types of diffusion equations.Are there other choices of CLBS which may yield new interesting results for these considered equations?
Once the CLBS and the governing equation are determined, exact solutions of the discussed equations can be easily constructed due to the compatibility of the considered equations and the additional differential constraints obtained by setting the characteristic of the corresponding CLBS to be zero.Thus, CLBS can be reformulated by using the technicalities of the method of differential constraints.
The key idea of the method of differential constraint was proposed by Yanenko in [22].The general formulation of this method requires that the original system of partial differential equations and the additional differential constraints satisfy some conditions of compatibility.The problem of finding all differential constraints compatible with certain equations can be more complicated than the investigation of solving the original equations.Hence, it is better to content oneself with finding constraints in some fixed classes, and these classes must be chosen using additional considerations.

Advances in Mathematical Physics
In practice, differential constraints related to separation of variables [23] and invariant subspaces [24] are powerful to seek for reductions of different types of evolution equations.Olver and Rosenau [25,26], Olver [27], Kaptsov [28], and Levi and Winternitz [29] showed that many reductions including Bluman and Cole's nonclassical method [6], partial invariance method [30], Clarkson and Kruskal's direct method [31], and Galaktionov's higher-order direct method [3] can be understood within the framework of differential constraints.
B-determining equations [32] and linear determining equations [33,34] were both proposed to construct differential constraints.The method of linear determining equations with free parameters is much easier to illustrate than the one of B-determining equations with undetermined functions.The linear determining equations are more general than the classical determining equations for Lie's generators [30].
In this paper, we will consider the general third-order CLBS with the characteristic of the nonlinear reaction-diffusion equations by using the linear determining equations.Equation ( 2) is widely used as mathematical models of heat conduction and diffusion, filtration of gases and fluids in porous media, process of chemical kinetics, and biological process.It is noted that   =   /  ( = 1, 2, 3) in (1).
There are numerous works aimed to study symmetry classifications and reductions of (2).CLBS was also used to consider the variant forms of (2).Three types of CLBS were applied to study such equations.The CLBS with the characteristic was presented to study functionally separable solutions of (2) in [13].The CLBS with the characteristic was proposed to study the variant forms of (2) in [11,12], which can give symmetry interpretation for first-order signinvariants of the considered equations.The CLBS with the characteristic was proposed to study the various forms of (2) in [15,16].The corresponding linearized form of ( 6) is related to the invariant subspace which is the right solution space of the linear ODE  = 0.
The forms of CLBS in [11][12][13][14][15][16] are all determined in advance.Thus the corresponding results are of course restricted to this ansatz.The third-order differential constraints (1) of the nonlinear diffusion equation   = (    )  are considered in [33].The differential constraints with  = 2, 3 of ( 2) are discussed in [34].However, the study there was worked out for only several special cases.Here we will extend the work there to develop a complete classification for third-order CLBS with the form (1) of ( 2) by applying the method of linear determining equations.New types of CLBS which are not included in the form of ( 3), ( 5), and ( 6) will be derived and therefore a huge number of new exact solutions of (2) will be constructed due to the compatibility of the governing equations and corresponding resulting differential constraints.
Definition 2 (see [1,2]).The evolutionary vector field ( 11) is said to be a CLBS of (10) if the following condition holds, where   denotes the set of all differential consequences of equation  = 0 with respect to .
It is noted that we can exclude all derivatives with respect to  on the set of (10) and thus obtain the vector field (11) with  of the form  =  (, , ,  1 ,  2 , . . .,   ) . ( A direct computation of ( 14) yields On the other hand, we can derive that where   denotes the set of all differential consequences of (10) with respect to .Thus, we have that the condition holds, which is exactly the sufficient condition [28,35] for which  = 0 is a differential constraint of the evolution equation (10).Thus  = 0 defines an invariant manifold of the equation (10), which is right corresponding to the symmetry reductions due to the compatibility of  = 0 and the considered (10).There is an equivalence relation between CLBS and differential constraint for the evolution equation (10).Kaptsov [33,34] proposed a method for finding differential constraints, which uses linear determining equations.The definition of linear determining equation was first given in [33], which is listed as follows.

CLBS of (2)
As stated in [34], the linear determining equation for the differential constraints  = 0 of (2) is right where  1 ,  2 ,  3 ,  4 ∈ R and   is the differential consequences of (2) with respect to .It is easy to see that (2) admits LBVF (11) if the condition (20) holds.
Substituting (1) into (20), the left-hand side and the righthand side are, respectively, simplified as which can be regarded as two polynomials with respect to  5 ,  4 , and  3 .It is noted that the subscripts here denote the differentiation with respect to the indicated variables.Considering the coefficients of  4 and  2 3 , we obtain So  can be represented by  = ℎ (, , ,  1 )  2 +  (, , ,  1 ) .
Here ℎ and  are both functions of , , , and  1 .Collecting the coefficients of  3 , we have Advances in Mathematical Physics Thus, we know that ℎ must satisfy Solving (25), we obtain and (, , ) is the function of , , and .As a consequence, (26) can be simplified as which implies Solving the linear ordinary differential equation about , we get  = (, ) −/2 .The simplification of the above equation about  gives We know that   =  1 must hold and then the above equations about  become Solving this ordinary differential equation, we know that there are three subcases for further discussion, which are listed as (1) (2) (3) For these three different cases, a further identical vanishing of (20) will finally determine the form of  in (1).The detailed procedure is omitted here.We just present the corresponding results in Table 1, where   ( = 1, 2, 3, 4, 5) are arbitrary constants.

Exact Solutions of (2)
In this section, we construct exact solutions of (2) due to the compatibility of differential constraints  = 0 and the governing equation.The equivalence relation between CLBSs and differential constraints for the evolution equation ( 10) is first rigorously proved in [36], where it is also stated that finding of the whole set of CLBSs of an evolution equation is a no-go problem.In fact, the general forms CLBSs  of (2) are completely determined by solving the linear system of determining equations about , which are listed in Table 1.
Solving  = 0, we can obtain  as a function of  with dependent undetermined functions.Substituting the form of  into (2) will yield that the -dependent coefficients in (2) satisfy three-dimensional dynamical system.Here we give several examples to illustrate the reduction procedure.
Example 1. Equation admits the CLBS The corresponding solutions are given as follows.