#### Abstract

We prove that any potential symmetry of a system of evolution equations reduces to a Lie symmetry through a nonlocal transformation of variables. This fact is in the core of our approach to computation of potential and more general nonlocal symmetries of systems of evolution equations having nontrivial Lie symmetry. Several examples are considered.

#### 1. Introduction

The Lie symmetries and their various generalizations have become an inseparable part of the modern physical description of wide range of phenomena of nature from quantum physics to hydrodynamics. Such success of a purely mathematical theory of continuous groups developed by Lie and Engel in 19th century [1] is explained by the remarkable fact that the overwhelming majority of mathematical models of physical, chemical, and biological processes possess nontrivial Lie symmetry.

One can even argue that this very property, invariance under Lie symmetries, distinguishes the popular models of mathematical and theoretical physics from a continuum of possible models in the form of differential or integral equations (see, e.g., [2, 3]). Based on this observation is the symmetry selection principle stating that if an equation describing some physical process contains arbitrary elements, then the latter should be so chosen that the resulting model possesses the highest possible symmetry. In this sense, the Lie theory effectively predicts which equation is the best candidate to serve as a mathematical model of a specific physical, chemical, or biological process.

The procedure of selecting partial differential equations (PDEs) enjoying the highest Lie symmetry from a prescribed class of PDEs is called group classification. In the case when non-Lie symmetries are involved, the more general term, symmetry classification, is used.

In this paper we study symmetries of systems of evolution equations in one spatial variable where , , , . Note that we use the boldface font to denote a multicomponent variable.

The problem of the Lie group classification of PDEs of the form (1.1) has been extensively studied (see, e.g. [4–7] and the references therein). The centerpiece of any approach used in this respect is the classical infinitesimal Lie method. The latter enables to reduce the problem of description of transformation groups admitted by (1.1) to integrating some linear system of PDEs (further details can be found in [8–10]).

However, with all its importance and power, the traditional Lie approach does not provide all the answers to mounting challenges of the modern nonlinear physics. By this very reason there were numerous attempts of generalization of Lie symmetries so that the generalized symmetries retain the most important features of Lie symmetries and allow for a broader scope of applicability. A natural move in this direction would be letting the coefficients of infinitesimal generators of the Lie symmetries to contain not only independent and dependent variables and their derivatives but also integrals of dependent variables, as well. In this way, the so-called nonlocal symmetries have been introduced into mathematical physics.

The concept of nonlocal symmetry of linear PDEs is relatively well understood (see, e.g., [11]). This is not the case for nonlinear differential equations. The problem of developing regular methods for constructing nonlocal symmetries of nonlinear PDEs is still waiting for its Sophus Lie. Still, there are a number of results on nonlocal symmetries for specific equations.

One of the possible approaches to construction of nonlocal symmetries has been suggested by Bluman et al. [12, 13]. They put forward the concept of potential symmetry, which is a special case of nonlocal symmetry. The basic idea of the method for constructing potential symmetries of PDEs can be formulated in the following way. Consider an evolution equation Suppose that it can be rewritten in the form of a conservation law By force of (1.3), we can introduce the new dependent variable and rewrite (1.1) as follows: If the system of two equations (1.4) admits a Lie symmetry such that at least one of the coefficients of its infinitesimal operator depends on , then this symmetry is the nonlocal symmetry for the initial evolution equation (1.2). Here is the inverse of , that is, . This nonlocal symmetry is also called potential symmetry of (1.2).

Pucci and Saccomandi [14] and Saccomandi [15] proved that potential symmetries can be derived using nonclassical symmetries of PDE (1.2). Recently, we established much stronger assertion by associating potential symmetries with classical contact symmetries [16, 17]. More precisely, we proved that any potential symmetry of evolution equation (1.2) can be reduced to contact symmetry by a suitable nonlocal transformation of dependent and independent variables. As a consequence, one can obtain exhaustive description of potential symmetries of (1.2) through classification of contact symmetries of PDEs of the form (1.2).

Some applications of potential symmetries to specific subclasses of the class of PDEs (1.2) can be found in [18–23].

In the present paper, we generalize the results of [16] for system of evolution equations (1.1) and prove that any potential symmetry of the system in question reduces to classical Lie symmetry under a suitable nonlocal transformation of dependent and independent variables (Sections 2 and 3). Next, we suggest in Section 4 a more general approach to constructing nonlocal symmetries that goes far beyond the concept of potential symmetries. It enables generating systems of evolution equations associated with a given system of the system (1.1), provided the latter admits a nontrivial Lie symmetry. Some applications of the approach in question are presented in Section 4.

#### 2. Conservation Law Representation and Classical Symmetries

*Definition 2.1. *One says that system (1.1) admits complete conservation law representation (CLR) if it can be written in the form
Here , and are -component vectors.

*Definition 2.2. *One says that system (1.1) admits partial CLR if it can be written in the form
Here and are -component and -component vectors, respectively.

Below we present theorems that provide exhaustive characterization of conservation law representability in terms of classical Lie symmetries. We give the detailed proof of the assertion regarding complete CLR; the case of partial CLR is handled in a similar way.

Theorem 2.3. *System (1.1) admits complete CLR if and only if it is invariant under -dimensional commutative Lie algebra , where
**
and besides
*

*Proof. *Suppose that system (1.1) admits CLR (2.1). Introducing new -component function
and eliminating from (2.1), we get
Integrating the obtained system of PDEs with respect to yields
Note that the integration constant is absorbed into the function . Evidently, system (2.7) is invariant under the commutative -dimensional Lie algebra . What is more is that the coefficients of the basis elements of the algebra satisfy condition (2.4).

Let us prove now that the inverse assertion is also true. Suppose that (1.1) admits Lie algebra , whose basis elements have the form (2.3) and satisfy (2.4). Then, there is a change of variables (see, e.g., [8])
reducing basis elements of to the form , . In what follows, we drop the bars.

Now, (1.1) necessarily takes the form
Differentiating (2.9) with respect to and making the (nonlocal) change of dependent variables , we finally get
which completes the proof.

*Note 1. *The fact that symmetry operators are of specific form (2.3) is crucial for the whole procedure of reducing a system of evolution equations to “conserved” form (2.1). If a symmetry group generated by some operator does not preserve the temporal variable (which means that the coefficient of in is nonzero for some ), then this operator cannot be reduced to the canonical form , and the reduction routine does not work.

Theorem 2.4. *System (1.1) admits partial CLR if and only if it is invariant under -dimensional commutative Lie algebra , where
**
and besides
*

#### 3. Potential Symmetries

Potential symmetries of system of evolution equations (1.1) appear in the same way as they do for a single evolution equation. For simplicity, we consider the case of complete CLR. By force of (2.1), we can introduce the new dependent variable , so that Note that is nonlocal variable since .

Suppose now that system (3.1) admits the Lie symmetry such that one of the derivatives does not vanish identically. Rewriting group (3.2) in terms of variables , and and taking into account that yield the nonlocal symmetry of the initial system of evolution equations (1.1). This means, in particular, that the symmetry in question cannot be obtained within the Lie infinitesimal approach. What we are going to prove is that this symmetry can be derived by regular Lie approach if the later is combined with the nonlocal transformation of the dependent variables.

Indeed, let system (1.1) admit complete CLR (2.1). In addition, we suppose that (1.1) possesses potential symmetry. Making the nonlocal change of dependant variables, , we rewrite (2.1) in the form (2.6). As initial system (1.1) admits a potential symmetry, system (3.1) is invariant under the Lie transformation group of the form (3.2).

Integrating (2.6) with respect to yields system of evolution equations Next, we rewrite the Lie symmetry (3.2) by eliminating according to (3.4) which yields By construction, Lie transformation group (3.6) maps the set of solutions of (3.5) into itself. Consequently, (3.6) is the Lie group of contact symmetries of system of evolution equations (3.5).

It is a common knowledge that any contact symmetry of a system of PDEs boils down to the first prolongation of a classical symmetry [24]. Consequently, the derivatives of , and with respect to the third argument vanish identically and we get This group is nothing else than the standard Lie symmetry group of system (3.5).

The same assertion holds true for the case of partial CLR.

Theorem 3.1. *Let system of evolution equations (1.1) admit complete or partial CLR and be invariant under a potential symmetry. Then, there exists a (nonlocal) change of variables mapping (1.1) into another system of the form (1.1) so that the potential symmetry of (1.1) becomes the standard Lie symmetry of the transformed system.*

This assertion is, in fact, the no-go theorem for potential symmetries of systems of evolution equations. It states that the concept of potential symmetry does not produce essentially new symmetries. The system admitting potential symmetry is equivalent to the one admitting the standard Lie symmetry, which is the image of the potential symmetry in question.

However, there is more to it. Theorem 3.1 implies the regular algorithm for group classification system of nonlinear evolution equations admitting nonlocal symmetries. Again, for the sake of simplicity, we consider the case of complete CLR.

Indeed, let system of evolution equations (1.1) be invariant under -dimensional Lie algebra . Here are commuting operators of the form (2.3) and their coefficients satisfy constraint (2.4). Basis operator is of the generic form

Making an appropriate change of variables, we can reduce the operators to the canonical forms, namely, , . Then, system (1.1) necessarily takes the form (3.5).

Let (3.7) be the Lie transformation group generated by the symmetry operator . Calculating the first prolongation of formulas (3.7) we get the transformation rule for the first derivatives of :

Now, we differentiate (2.6) with respect to and make the following change of dependent variables: which yields Formulas (3.7), (3.9) provide the image of the transformation group (3.7) under the mapping (3.9), so that Here .

Consequently, if one of the derivatives, , , , does not vanish identically, then (3.12) is the nonlocal symmetry group of system of evolution equations (3.11).

The same line of reasoning applies to the case when system (1.1) admits partial CLR.

We summarize the above speculations in the form of the procedure for computation of nonlocal symmetries of systems of evolution equations associated with a given system of the form (1.1).

Let system of evolution equations (1.1) be invariant under -dimensional Lie symmetry algebra . For simplicity, we consider the case of complete CLR.

*Procedure. **Classification of Potential Symmetries of (1.1)*(1)Calculate inequivalent subalgebras of the algebra . (2)Select those subalgebras , which contain commutative subalgebras of operators of the form (2.3). (3)For each commutative subalgebra perform change of variables reducing its basis elements to the canonical forms and transform the initial system (1.1) accordingly. (4)Perform nonlocal transformation (3.10). (5)Eliminate “old” dependent variables from (3.7) in order to derive symmetry group (3.12) of the transformed system of evolution equations (3.11). (6)Verify that there is, at least, one derivative from the list , , that does not vanish identically. If this is the case, then (3.12) is the nonlocal (potential) symmetry of (3.11).

The steps needed to implement the above procedure for the case of system of evolution equations admitting partial CLR are the same, the only difference is that intermediate formulas (3.7)–(3.12) are more cumbersome, since we need to distinguish between two sets of dependent variables and (see, (2.2)).

Note that by force of Theorems 2.3 and 2.4, any potential symmetry of equations of the form (1.1) can be obtained in the above-described manner.

As an example, we consider the Galilei-invariant nonlinear Schrödinger equation introduced in [25] where , , is an arbitrary real constant, and is an arbitrary complex-valued function. Equation (3.13) admits the Lie algebra of the Galilei group having the following basis operators [25]:

Operators , commute and the rank of the matrix of coefficients of operators , , , and is equal to 2. Consequently, there is a change of variables that reduces , to canonical forms , . Indeed, making the change of variables transforms , to become , . So we can apply Procedure 1 to (3.13) transformed according to (3.15). As a result, the transformed operator becomes the potential symmetry of the transformed nonlinear system of two evolution equations.

#### 4. Some Generalizations

Denote the class of partial differential equations of the form (1.1) as . Then any system of the form belongs to , and (ii) its image under nonlocal transformation also belongs to . Existence of such nonlocal transformation is in the core of our approach to classifying nonlocal symmetries of systems of evolution equations.

It is not but natural to ask whether there are other types of nonlocal transformations of the class that can be utilized to generate nonlocal symmetries. Remarkably, such nonlocal transformations do exist. Sokolov [26] put forward the idea of group approach to generating such transformations for a single evolution equation. It is straightforward to modify his approach to handle systems of evolution equations, as well. As an illustration, we consider system (4.1). It is invariant under the -dimensional Lie algebra . The simplest set of functionally-independent invariants of the algebra can be chosen as follows: , , , , . Now, we define the transformation so that , , are invariants of the symmetry group of the system under study. In the case under consideration, we have , , and . As we established in Section 2, applying this transformation to any equation of the form (4.1) yields system of evolution equations that belongs to . What is more, is that Lie symmetry group of (4.1) is mapped into symmetry group of the transformed system and some of the basis operators of the latter become nonlocal ones.

Consider as the next example system of evolution equations This system is invariant under the -dimensional Lie algebra . The simplest set of functionally independent first integrals reads as , , , . Consequently, change of variables (4.2) takes the form Note that we dropped the bars and replaced with .

Transforming (4.3) according to (4.4) we get or, equivalently, Integrating twice yields Note that integration constants are absorbed by the function .

So that nonlocal transformation (4.4) maps a subset of equations from into . Consequently, it can be used to generate nonlocal symmetries of the initial system (4.3).

Let system (4.3) be invariant under the Lie transformation group Computing the second prolongation of the above formulas, we get the transformation law for the functions , Combining (4.8) and (4.9) yields the symmetry group of system of evolution equations (4.7), where are nonlocal variables. Now, if one of the derivatives does not vanish identically, then (4.10) is the nonlocal symmetry group of system of evolution equations (4.7).

It is important to emphasize that the symmetry algebra is not obliged to be commuting. The necessary condition is that the corresponding transformation group has to preserve the temporal variable, , that is, basis elements of have to be of the form

As an illustration, consider the following system of second-order evolution equations: This system is invariant under the -dimensional Lie algebra . Note that the algebra is not commutative. The set of invariants of the algebra can be chosen as follows: Making the change of variables we rewrite (4.13) in the form Taking into account that the operators and commute, differentiating (4.16) with respect to , and replacing with , we finally get The above system is obtained from the initial one through the change of dependent variables , . Consequently, if system (4.13) admits symmetry (4.8), then system (4.17) admits the following transformation group: with . Again, if one of the derivatives does not vanish identically, then (4.18) is the nonlocal invariance group of system of evolution equations (4.16).

The procedure for calculation of nonlocal symmetries of system (1.1) suggested in the previous section yields those nonlocal symmetries which are potential, since the nonlocal transformation was chosen *a priori*. Allowing for a nonlocal transformation to be determined by symmetry group of the system under study yields a more general algorithm for constructing nonlocal symmetries.

Let system of evolution equations (1.1) be invariant under -dimensional Lie symmetry algebra . Then the following procedure can be used to construct nonlocal symmetries of (1.1).

*Procedure. **Classification of Nonlocal Symmetries of (1.1)*(1)Calculate inequivalent subalgebras of the algebra . (2)Select those subalgebras , which contain basis elements of the form (4.12).(3) For each construct functionally independent invariants , and make change of variables
(4)Eliminate “old” dependent variables from (4.20) in order to derive symmetry group of the transformed system of evolution equations. (5)Verify that there is, at least, one function from the list that depends on for some . If this is the case, then is the nonlocal symmetry of (3.11).

#### 5. Conclusion

One of the principal results of the paper is Theorem 3.1 stating that any potential symmetry of system of evolution equations (1.1) reduces to a Lie symmetry by an appropriate nonlocal transformation of dependent and independent variables. The nonlocal transformation in question is a superposition of the local change of variables and of the nonlocal change of dependent variables The explicit form of transformations (5.1) is defined by the Lie symmetry admitted by the corresponding system (1.1).

We obtain as a by-product exhaustive characterization of systems (1.1) that can be represented in the form of conservation law(s), in terms of the Lie symmetries preserving the temporal variable, , (see Theorems 2.3 and 2.4).

In Section 4, we generalize the above reasoning in order to obtain nonlocal symmetries which are not potential. The basic idea is replacing (5.2) with a more general nonlocal transformation. This transformation is determined by invariants of the Lie symmetry algebra of the system under study.

We intend to devote one of our future publications to systematic study of nonlocal symmetries of systems of nonlinear evolution equations (1.1) within the framework of the approach developed in Section 4.