Algebraic method of group classification
for semi-normalized classes of differential equations
Célestin Kurujyibwami†, Dmytro R. Popovych‡§ and Roman O. Popovych‡◇
† College of Science and Technology, University of Rwanda, P.O. Box: 3900, Kigali, Rwanda
‡ Department of Mathematics and Statistics, Memorial University of Newfoundland,
St. John’s (NL) A1C 5S7, Canada
§ Institute of Mathematics of NAS of Ukraine, 3 Tereshchenkivska Str., 01004 Kyiv, Ukraine
◇ Mathematical Institute, Silesian University in Opava, Na Rybníčku 1, 746 01 Opava, Czechia
E-mails: [email protected], [email protected], [email protected]
We generalize the notion of semi-normalized classes of systems of differential equations, study properties of such classes and extend the algebraic method of group classification to them. In particular, we prove the important theorems on factoring out symmetry groups and invariance algebras of systems from semi-normalized classes and on splitting such groups and algebras within disjointedly semi-normalized classes. Nontrivial particular examples of classes that arise in real-world applications and showcase the relevance of the developed theory are provided. To convincingly illustrate the efficiency of the proposed method, we apply it to the group classification problem for the class of linear Schrödinger equations with complex-valued potentials and the general value of the space dimension. We compute the equivalence groupoid of the class by the direct method and thus show that this class is uniformly semi-normalized with respect to the linear superposition of solutions. This is why the group classification problem reduces to the classification of specific low-dimensional subalgebras of the associated equivalence algebra, which is completely realized for the case of space dimension two. Splitting into different classification cases is based on three integer parameters that are invariant with respect to equivalence transformations. We also single out those of the obtained results that are relevant to linear Schrödinger equations with real-valued potentials.
1 Introduction
Systems of differential equations that appear in real-world applications generally contain constant or functional parameters that are not fixed and should be found using experimental data or additional arguments. An example of such parameters called arbitrary elements is given by potentials in Shrödinger equations. The scope of relevant models is often restricted by certain requirements to their symmetry properties following from the corresponding theory, like certain scale invariance or the Galilean or special relativity principles in physics. In mathematical terminology, this means that modeling differential equations have to possess a certain symmetry group or symmetry groups with certain properties or the widest symmetry group among the possible ones [61]. Thus, to single out such models from a parameterized class of systems of differential equations, one should solve the group classification problem for this class. This is why the theory of group classification plays nowadays a central role in symmetry analysis of differential equations and its application to various sciences.
This theory has its origins in the work of Sophus Lie. His greatest achievement here was solving the group classification problems for three classes differential equations, second-order ordinary differential equations [41], second-order linear partial differential equations in two independent variables [40] and nonlinear Klein–Gordon equations of the form [39]. Both the main approaches to solving group classification problems, which are the direct and the algebraic methods, also originated from these Lie’s studies. The modern development of the theory of group classification was started by Ovsyannikov [60, 61], and this has been the field of intensive research ever since, see [2, 27, 38, 68] and references therein.
The group classification problem for a class of differential equations is complex beyond constructing the Lie symmetry group of a single system of differential equations. It involves classifying solutions of a cumbersome overdetermined system of partial differential equations, where the unknown functions are not only the coefficients of Lie-symmetry vector fields but also the arbitrary elements of the class. This is why the progress in implementing the group classification framework in software is not significant despite there existing many specialized packages for finding Lie symmetries in various symbolic computation systems [3, 14, 15, 16, 43, 72]. Moreover, a considerable fraction of group classification problems presented in the literature is tackled using a brute-force version of the direct method, which is advisable only for some simple classes of differential equations, see a discussion in [58, Section 4]. This definitely inefficient approach requires a rigorous effort to keep track of various cases of integrating the corresponding determining equations for Lie symmetries and furthermore to check inequivalence of these cases. Any lapse in attention leads to missed, overlapping, repeated and/or incorrect classification cases afflicting most papers following this approach. Due to their unhandiness, the corresponding computations are not properly presented if presented at all, which makes the obtained classification results unverifiable and thus not trustworthy.
In the meanwhile, a number of novel techniques for group classification have been developed recently, which include variations to the advanced modification of the direct method called the method of furcate splitting [29, 58, 55] and various flavors of the algebraic method [7, 17, 35, 57, 67, 68]. The former method was suggested and applied for the first time in [55] but was given its name only in [29] and formalized to its full extend in [58]. See also [8, 65, 69, 70] and references therein for examples of its efficient application. It is advisable only for classes of simple structure, which is, at the same time, not appropriate to apply the more powerful algebraic method. In particular, arbitrary elements in such classes should be constants or depend on few arguments, and the reasonable upper bound for the number of arbitrary-element arguments, which was handled in the literature [55], is two. At the same time, the latter method is not so sensitive to the number of arguments of arbitrary elements.
Despite its invention by Lie, the algebraic method of group classification became a common tool of group analysis of differential equations considerably later, only since the 1990s [4, 24, 25, 26, 37, 38, 42, 62, 71, 73]. Moreover, the fact that the class of differential equations under consideration should admit special (normalization) properties for the proper applicability of this method was understood even later [63, 64, 67]. The straightforward application of the algebraic method to non-normalized classes of differential equations results in the so-called preliminary group classification of such classes [1, 17, 28], see the discussion in [7, Section XI]. To extend the scope of applicability of the algebraic method, a number of various techniques and approaches were invented, including splitting the class under consideration into normalized subclasses [36, 67], constructing its normalized superclass or even a hierarchy normalized classes related to it [9, 36, 57, 67], singling out its singular subclasses [7], gauging its arbitrary elements using a wide subgroupoid of its equivalence groupoid [9, 36, 57, 58, 67], mappings between classes that are generated by point transformations [58] and partitioning classification cases into regular and singular ones [68]. One of such approaches is to extend the algebraic method to classes with weaker normalization properties, e.g., to classes that are semi-normalized or normalized in the generalized or the extended sense. The first version of the algebraic method for group classification of semi-normalized classes was suggested in [35] and applied therein for solving the group classification problem for the class of (1+1)-dimensional linear Schrödinger equations with complex-valued potentials. Then this version was used to classify Lie symmetries of (1+1)-dimensional linear evolution equations of arbitrary order greater than two [9]. It was shown in [58] that further development of the algebraic method of group classification requires a deep revision of foundations of the theory of group classification, including the statement of group classification problems.
In the present paper, we enhance and further generalize the above version of this method and apply it to the class of (1+)-dimensional () linear Schrödinger equations with complex-valued potentials
(1) |
where and are the real independent variables, is the complex dependent variable and is an arbitrary smooth complex-valued potential depending on and . Here and in what follows, the subscripts and denote differentiation with respect to and , respectively, the indices , , , and run from to , and we use summation convention over repeated indices. By we denote the subclass of the class that consists of the equations of the form (1) with real-valued potentials .
More specifically, we generalize the notion of semi-normalized classes of differential equations, study properties of such classes and extend the algebraic method of group classification to them. Basic for this extension are the new theorems on factoring out symmetry groups and invariance algebras of systems from semi-normalized classes and on splitting such groups and algebras within disjointedly semi-normalized classes. For the class with the general value of , we compute its equivalence groupoid by the direct method and thus prove its uniform semi-normalization with respect to the linear superposition of solutions. In this way, we reduce the group classification problem for the class to the classification of appropriate low-dimensional subalgebras of its equivalence algebra. This problem is completely solved in the case of . Splitting into different classification cases is based on three integer parameters. These parameters are invariant with respect to equivalence transformations and characterize the dimensions of parts of the corresponding Lie symmetry algebras, which are related to transformation of time, rotations and generalized time-dependent shifts with respect to space variables, respectively. The analogous results for the subclass are singled out from those above. In the study of the class , the present paper succeeds [35], where the case was comprehensively studied, and [34], where the groupoid was computed for the case .
The study of symmetries of Schrödinger equations began in the 1970s with Lie symmetries of linear Schrödinger equations with real-valued time-independent potentials [11, 44, 47, 48, 49, 50], but later the focus shifted to nonlinear Schrödinger equations [5, 12, 19, 20, 21, 22, 23, 24, 36, 62, 66, 67, 74], see also detailed reviews in [67, Section 4] and [35, Section 1]. At the same time, Lie symmetries of linear Schrödinger equations and their generalizations were not properly classified. In particular, this is the case for linear Schrödinger equations with complex-valued potentials although they also arise as models in quantum mechanics, scattering theory, condensed matter physics, quantum field theory and so forth [6, 18, 46]. Moreover, the existing results for the standard case of real-valued time-independent potentials are neither exhaustive nor completely correct, which has been found out in the course of the consideration of more general Schrödinger equations. Thus, in [54], Schrödinger equations with position-dependent mass and stationary real-valued potentials were considered from the Lie-symmetry point of view for arbitrary space dimension, and their Lie symmetries were completely classified for space dimension two. The obtained classification list contains three case families formally intersecting the case of constant mass. As a result, a Lie-symmetry classification case for Schrödinger equations with constant mass and real-valued time-independent potentials that had been missed in [11] was identified; see the conclusion of the present paper for more details. Hence there are a number of open problems even on Lie symmetries of various Schrödinger equations, including linear ones. One can consider wider or new classes of Schrödinger equations [51, 54] or carry out complete, partial or inverse group classification of models with additional properties such as superintegrability [52, 53].
The structure of this paper is as follows. Basic notions related to classes of differential equations are reviewed in Section 2. In Section 3, we extend the notion of a semi-normalized class of differential equations to the semi-normalization with respect to a uniform family of point-symmetry groups of systems from the class under consideration and a subgroup of the corresponding equivalence group, study properties of newly defined semi-normalized classes and prove the theorems on factoring out point-symmetry groups and maximal Lie invariance algebras for systems from such classes. In Section 4, we distinguish between the disjoint and the non-disjoint semi-normalizations, study specific properties of disjointedly semi-normalized classes of differential equations and prove the stronger theorems on splitting point-symmetry groups and maximal Lie invariance algebras within such classes. We also present two important examples, the class that is disjointedly semi-normalized only with respect to a proper subgroup of its equivalence group and a semi-normalized class that is not disjointedly semi-normalized. The additional attribute “uniform” for semi-normalization is introduced in Section 5. Therein we discuss an example of a class that is semi-normalized but not uniformly semi-normalized. More attention is paid to the particularly important case of uniform semi-normalization given by the classes of homogeneous linear systems of differential equations that are uniformly semi-normalized with respect to the linear superposition of solutions. This specific semi-normalization is especially relevant to the symmetry analysis of the linear Schrödinger equations in the second part of the present paper. The algebraic method of group classification is extended to semi-normalized classes in Section 6.
Section 7 is devoted to the preliminary symmetry analysis of the class of (1+)-dimensional linear Schrödinger equations of the form (1) with complex-valued potentials and with an arbitrary value . We compute the equivalence groupoid , the equivalence group and the equivalence algebra of this class. It is proved that the class is uniformly semi-normalized with respect to the linear superposition of solutions. We derive and preliminarily integrate the determining equations for Lie symmetries of equations from the class , which allows us to obtain the general form of Lie-symmetry vector fields, the kernel Lie invariance algebra of the class as well as the classifying condition for such Lie symmetries. Then we analyze properties of subalgebras of (more precisely, of the projection of to the space of independent and dependent variables) that are appropriate as maximal Lie invariance algebras of equations from the class . In particular, we obtain that the least upper bound for the dimensions of these subalgebras equals . As a result, we show that the group classification problem for the class reduces to the classification of appropriate subalgebras of the algebra with respect to the equivalence relation generated by the action of the projected group . The group classification problem for the class in the case is completely solved in Section 8. We use the above results on the class in Section 9 to derive their counterparts for the class .
Lastly, we conclude with a summary of results, their discussion and suggestions for future considerations in Section 10.
Conventions.
Throughout the paper, we assume summation over repeated indices. Subscripts of functions denote differentiation with respect to the corresponding variables. Given a class of differential equations, denotes the natural projection of the space run by the tuple of the independent and dependent variables, relevant jet variables and the arbitrary elements of the class onto the space run by the tuple of the independent and dependent variables only. When talking about pseudogroups, we omit the prefix “pseudo”.111As for many other classes of differential equations, the equivalence group of the class and the point symmetry groups of equations from this class are in fact pseudogroups but not groups. The theory developed should be formulated involving the notions pseudogroups and pseudoalgebras. At the same time, such a modification essentially complicates the consideration and requires revising and enhancing the entire framework of symmetry analysis of differential equations, beginning with the notions of a system of differential equations and of a class of such systems, which we intend to carry out in the future. The notation is used for the identity transformation in the relevant space. In Sections 2–6 with theoretical results on general systems of differential equations and in Sections 7 and 8 with specific results on Schrödinger equations, we use different locally specified notations.
2 Basics notions on classes of differential equations
In order to make the presentation self-contained, we briefly define the notion of a class (of systems) of differential equations and various objects related to point transformations in such classes, mainly following [68, Section 2]. More details can be found in [7, 59, 63, 67, 68] and references therein.
Notation.
Up to Section 6, we use a notation for independent and dependent variables that differs from that of the other sections: denotes the complete -tuple of independent variables and is the -tuple of dependent variables (that is, unknown functions of the independent variables), except particular examples, where natural specific notations of variables are used. The indices and run from to and the indices and run from to .
Consider a system of differential equations : parameterized by a tuple of arbitrary elements where stands for the tuple of derivatives of with respect to up to order , including as the zeroth order derivatives, and denotes the tuple of derivatives of the tuple with respect to its arguments up to order . The arbitrary-element tuple runs through the set of solutions of an auxiliary system of differential equations and differential inequalities of the form , where both and play the role of independent variables, and and are tuples of smooth functions depending on and . Other kinds of inequalities are also possible. By we denote the class of systems with arbitrary elements running through , .
The (point) equivalence groupoid of the class is the small category with as the set of objects and with the set of point transformations222In the case of single dependent variable, one can also consider the contact equivalence groupoid of the class that consists of the admissible contact transformations within this class. of , i.e., of local diffeomorphisms in the space with the coordinates , between pairs of systems from as the set of arrows,
More specifically, elements of are called admissible (point) transformations within the class . The source and target maps are defined by . This map pair inspires the groupoid notation . Admissible transformations and are composable if , and then their composition is , which defines a natural partial multiplication on . For any , the unit at is given by , where is the identity transformation of . This defines the object inclusion map , i.e., the object set can be regarded to coincide with the base groupoid . The inverse of is , where is the inverse of . All required groupoid properties like associativity of the partial multiplication, its consistency with the source and target maps, natural properties of units and inverses are obviously satisfied. The subset of with corresponds to the set of point transformations mapping the system to the system . The vertex group is associated with the point symmetry (pseudo)group of the system , . The fundamental groupoid of the class is the disjoint union of the vertex groups , ,
Since it has the same object set and the same vertex groups as and for any , it is a normal subgroupoid of , which is also called the fundamental subgroupoid of .
Denote by , , the natural projection from the space with the coordinates to the spaces with the coordinates . In this notation, .
The (usual) equivalence group of the class is the (pseudo)group of point transformations in the space with the coordinates that
-
•
are projectable to the spaces with the coordinates for any ,
-
•
are consistent with the contact structure on the space with the coordinates , i.e., is the standard prolongation of to th order jets , and
-
•
map the class onto itself, for any .
If the arbitrary-element tuple depends at most on with , then the group can be considered to act in the space with the coordinates . We call the group the shadow of . Two specific restrictions on equivalence transformations, their projectability and their locality with respect to arbitrary elements, can be weakened to generalize the notion of usual equivalence group in several ways. This results in the notions of generalized equivalence group and extended equivalence group, respectively, or the notions of extended generalized equivalence group if both restrictions are weakened simultaneously [29, 45, 57, 58, 63, 67].
The action groupoid a subgroup of the usual equivalence group of the class ,
is a wide subgroupoid333A subgroupoid of a groupoid is called wide if the sets of objects of and coincide, . of the equivalence groupoid of this class, and . We say that an admissible transformation in the class is generated by an equivalence transformation from if .
The kernel point symmetry group of systems from the class is related to the unfaithful subgroup of under the action on , . Thus, is a normal subgroup of that is isomorphic to , and the isomorphism is established by .
Given a transformation group, it is often required to consider the associated Lie algebra of vector fields, which consists of the generators of one-parameter subgroups of this group. The infinitesimal counterparts of the equivalence group and the kernel point symmetry group of the class and the point symmetry group of a system are called the equivalence algebra and the kernel invariance algebra of the class and the maximal Lie invariance algebra of the system , respectively.
If the equivalence groupoid has specific properties, then the class is more convenient for studying within the framework of symmetry analysis of differential equations [7, 35, 63, 67, 68].
Definition 1.
The class is called normalized if . It is called semi-normalized if .
The Frobenius product is a well-defined subgroupoid of since the fundamental groupoid is a normal subgroupoid of . It is obvious that any normalized class of differential equations is semi-normalized.
When considering different kinds of equivalence groups for the class , the usual, the generalized, the extended or the extended generalized equivalence groups of , we should specify the corresponding kinds of (semi-)normalization, (semi-)normalization in the usual, the generalized, the extended or the extended generalized sense. By default, (semi-)normalization means (semi-)normalization in the usual sense.
More general notion of semi-normalization is presented in the next section.
3 Semi-normalized classes of differential equations
By the definition introduced in [35, Section 3], a class of differential equations is uniformly semi-normalized if there exists a family of point-symmetry groups of systems from this class with certain properties with respect to its equivalence groupoid and its equivalence group, and these point-symmetry groups can be proper subgroups of the corresponding complete point-symmetry groups. At the same time, it was noted therein that a proper subgroup of the equivalence group can be used instead of the entire equivalence group as well. The note was important since already in [35, Section 7] we needed to classify Lie symmetries of a class that can be assumed as uniformly semi-normalized only with respect to a proper subgroup of its equivalence group. Carefully analyzing the above results, we also see that the group classification problem can be solved by the algebraic method for classes that possess at least a part of properties required for uniformly semi-normalized classes in [35]. This is why, following terminology and notation of [59, 68], in the present paper we change the terminology and extend the algebraic method of group classification to more general classes of differential equations with weaker normalization properties.
Definition 2.
Let be the action groupoid of a subgroup of . Suppose that a family of subgroups of the point symmetry groups with the associated subgroups of the vertex groups for each satisfies the property for any . In other words, the family is uniform with respect to the action of on , or -uniform for short. Then the Frobenius product
with is a subgroupoid of , which coincides with the image of under the -action (resp. the -action) of on . If , we call the class semi-normalized with respect to the subgroup of and the family of subgroups of the point symmetry groups, or -semi-normalized for brevity.
Remark 3.
Since is the usual equivalence group of the class , we can say more precisely that this class is -semi-normalized in the usual sense. The other kinds of semi-normalization, in the generalized sense, in the extended sense or in the generalized extended sense, which are associated with equivalence groups of the respective kinds, will be introduced elsewhere.
The case is degenerate; then necessarily for any , and thus and , which is not of interest in the sense of group classification. At the same time, the other improper case, where the corresponding subgroup of coincides with the entire equivalence group , is important. The class that is semi-normalized with respect to the group and the family is genuinely normalized since then . Here and in what follows denotes the identity transformation in the space with the coordinates . If the class is semi-normalized with respect to the group and the family , then it is literally semi-normalized. It is obvious that a normalized class or a class that is semi-normalized with respect to a subgroup of its equivalence group and a family of subgroups of the point symmetry groups of systems from this class is semi-normalized in the above specific sense. More generally, if , then the class is called semi-normalized with respect to the family of subgroups of the point symmetry groups, i.e., we avoid denominating in this specific situation. Moreover, we contract the term “a class that is semi-normalized with respect to a subgroup of its equivalence group and a family of subgroups of the point symmetry groups of systems from this class”, when talking about such classes in general, and just call them semi-normalized.
Theorem 4.
Suppose that the class is semi-normalized with respect to the subgroup of and the symmetry-group family .
(i) Then , is a normal subgroupoid of and for any .
(ii) Furthermore, for each , is a normal subgroup of , is a subgroup of , and the group is the Frobenius product of and , .
(iii) Systems from the class are similar (i.e., they are related by point transformations) if and only if they are -equivalent.
Proof.
(i) The -semi-normalization of the class means that any element can be represented in the form for some and some . Due to the uniformity of , there exists such that . In other words, any admissible transformation can also be represented in the form for some and some .
It is obvious that is a wide subgroupoid of (as well as ), and it coincides with its fundamental subgroupoid. In the notation of the previous paragraph, for any we have
which directly implies that is a normal subgroupoid of . In other words, under the condition , the uniformity of with respect to the actions of on implies the uniformity of with respect to the similar actions of the entire group or even of the entire groupoid .
(ii) We fix an arbitrary and take an arbitrary . Then and, since , the point transformation admits the factorization for some and some . The element of the family is a subgroup of , , and hence the point transformation also belongs to , and consequently to , which is a group as intersection of two groups and is thus a subgroup of . This implies that for any we have the representation , where and , i.e., . It is obvious that as well, which implies .
The fact that is a normal subgroupoid of directly implies that is a normal subgroup of , i.e., is a normal subgroup of . We can also immediately prove the last claim in terms of groups. For an arbitrary and an arbitrary , we consider the composition . As an element of , the transformation admits the factorization with some and some . Since , we have . Due to the uniformity of the family with respect to the action of on , we obtain . Hence the composition belongs to . Thus we have that is a normal subgroup of , .
(iii) If the systems and from the class are related by a point transformation in the space with coordinates , , then is an admissible transformation of this class, . Hence there exist and such that . The existence of means that the systems and are -equivalent. ∎
The member of the family and the subgroup of are called the uniform point-symmetry group and the essential point-symmetry group of the system that are associated with the -semi-normalization of the class , respectively. The knowledge of a family of uniform point-symmetry groups trivializes them in the following sense: since in view of Theorem 4(ii) the group is factorized over for each , then we only need to find the subgroup in order to construct . If either for any or for any , then the factorization is trivial. Theorem 4(i) implies that the family is not only -uniform but also -uniform and even -uniform. We can interpret each of these uniformities as the equivariance of the map under the analogous action on and on .
It also follows from Theorem 4(i) that the condition in Definition 2 can be replaced by the condition or even by the weaker condition since these conditions are equivalent if the family is uniform with respect to the action of on . Roughly speaking, each of these conditions means that the entire equivalence groupoid is generated by distinguished equivalence transformations and transformations from uniform point symmetry groups.
Given a -semi-normalized class , one can extend the subgroup of and elements of in such a way that this class will be semi-normalized with respect to the extended objects.
Theorem 5.
Let the class be semi-normalized with respect to the subgroup of and the symmetry-group family , and let denote the kernel group of . Then the class is - and -semi-normalized, in particular - and -semi-normalized, where with denoting the Frobenius product of and and is any subgroup of that contains , .
Proof.
For each , the uniform symmetry group and the kernel symmetry group are a normal subgroup and a subgroup of in view of Theorem 4 and by the definition of , respectively. Hence the Frobenius product is a subgroup for any .
In view of Theorem 4(i), the symmetry-subgroup family is uniform with respect to any subgroup of . We show that this is also the case for the symmetry-subgroup family .
We take arbitrary and . By the definition of , there exists and such that . The uniformity of with respect to the action of on implies that there exists such that . Since and , we have for some . As a result,
where .
It is obvious that if . As a result, . The two limit cases are and . ∎
We call the claim from Theorem 4(ii) on the factorization the theorem on factoring out point-symmetry groups within semi-normalized classes. The infinitesimal version of this claim can be called the theorem on factoring out maximal Lie invariance algebras within semi-normalized classes. This version follows immediately from Theorem 4 if we replace the groups by the corresponding algebras of generators of the one-parameter subgroups of these groups.
Theorem 6.
Suppose that a class of differential equations is semi-normalized with respect to a subgroup of and a family of symmetry subgroups , and is a subalgebra of that corresponds to the subgroup . Then, for each , the Lie algebras and that are associated with the groups and are a subalgebra and an ideal of the maximal Lie invariance algebra of the system , respectively. Moreover, the algebra is the sum of and , , and .
4 Disjoint and non-disjoint semi-normalization
Definition 7.
Under the conditions of Definition 2, we call the class disjointedly semi-normalized with respect to the subgroup of and the family of subgroups of the point symmetry groups, or disjointedly -semi-normalized for short, if in addition for any .
If the class is disjointedly -semi-normalized, then .444Indeed, let . Then for some , and . Hence , and thus , i.e., . Therefore, , which means that . The inverse implication does not hold in general555Let be the class of the (1+1)-dimensional linear evolution equations , where all variables and parameters are real, the arbitrary-element tuple is with arbitrary constants and . The arbitrary element is artificial and inessential and can be set to be equal to zero by gauge equivalence transformations. For the class , we have and with the action groupoid of the gauge equivalence group of the class . The equivalence group of this class is constituted by the point transformations , , , , , where and are arbitrary constants with , and is an arbitrary sufficiently smooth function of with . The group is the normal subgroup of singled out by the constraints and . For any , the group consists of the point transformations , , , where and are arbitrary constants with , and is an arbitrary solution of . We take the symmetry-subgroup family and the subgroup of constituted by the transformations , , , , . The class is -semi-normalized, and for any . but is obvious if the class does not admit gauge admissible transformations.666Fix an arbitrary . Let , i.e., and for some . Then , , and thus since the class does not admit gauge admissible transformations. Therefore, , which means that . As a result, for any .
In particular, the class that is literally normalized is in fact disjointedly semi-normalized with respect to the group and the family . It is also clear that each disjointedly semi-normalized class is literally semi-normalized.
Disjointedly semi-normalized classes have nicer properties than those presented in Theorem 4 for general semi-normalized classes.
Theorem 8.
If the class is disjointedly semi-normalized with respect to a subgroup of and a symmetry-subgroup family , then
(i) any admissible transformation possesses a unique decomposition of the form with and as well as one of the form with the same and with ;
(ii) for each , the point symmetry group of the system splits over ; more specifically, is a normal subgroup of , is a subgroup of , and the group is the semidirect product of acting on , ;
(iii) for each , the Lie algebras and that are associated with the groups and are a subalgebra and an ideal of the maximal Lie invariance algebra of the system , respectively; moreover, the algebra is the semidirect sum , and , where is a subalgebra of the equivalence algebra of that corresponds to the subgroup .
Proof.
(i) Taking into account the proof of Theorem 4, it is only required to show the uniqueness of the first decomposition. The uniqueness of the second decomposition is proved in the same way. We fix an arbitrary admissible transformation of the class and consider two decompositions of the first kind for , , where and . Then . Since and are subgroupoids of , the left- and right-hand sides of the last equality belong to and , respectively. This means that they both belong to . Hence and .
(ii) In view of the claim on factoring out the groups from Theorem 4, it suffices to show that , but this is a direct consequence of the disjointness property since
Thus, splits over for each .
(iii) This item is the infinitesimal version of item (ii) and follows immediately from (ii) via replacing the groups with the corresponding algebras of generators of the one-parameter subgroups of these groups. It can also be considered as a consequence of Theorem 6. ∎
We would like to emphasize that augmenting the semi-normalization with the disjointedness in Theorem 8 additionally gives the uniqueness of the decompositions in item (i), replacing the Frobenius product by the semidirect product in item (ii) and replacing the sum of vector spaces by the semidirect sum of algebras in item (iii). We call item (ii) the theorem on splitting point-symmetry groups in disjointedly semi-normalized classes and item (iii) the theorem on splitting of maximal Lie invariance algebras in disjointedly semi-normalized classes. In contrast to Theorem 8(ii), the product of and in Theorem 4(ii) is in general a product of subgroups with nontrivial intersection as a product of group subsets but not a semidirect product of subgroups since the condition necessarily holds for any only in disjointedly semi-normalized classes.
Remark 9.
There are classes of differential equations that are disjointedly semi-normalized only with respect to proper subgroups of the corresponding equivalence groups. The first example of such a class was constructed in [35, Section 7] in the course of group classification of (1+1)-dimensional linear Schrödinger equations with complex-valued potentials depending on . Here we reinterpret it in the context of the above claim. More specifically, let in (1). Consider the subclass of the class that consists of the equations of the form (1), where the potential is purely imaginary and (homogeneously) linear with respect to . This subclass can be reparameterized to the class of equations of the form
where is an arbitrary smooth real-valued function of , which is the only arbitrary element.
The equivalence groupoid of the class consists of triples of the form , where is a point transformation in the space of variables with components | |||
(2a) | |||
(2b) | |||
the relation between and is given by | |||
(2c) |
, , , , and are arbitrary real constants with , is a nonzero complex constant, is an arbitrary solution of the initial equation , and . The (usual) equivalence group of the class is constituted by point transformations in the extended space of with components of the form (2), where and . For each , the equation admits the group of point symmetry transformations of the form (2a)–(2b) with and . Elements of the kernel symmetry group of the class additionally satisfy the constrains and , i.e., it consists of the ‘scalings’ of with nonzero complex numbers, : , , . The class is disjointedly semi-normalized with respect to the family and the subgroup of singled out from by the constraint . Admissible transformations with , and have no counterparts among equivalence transformations and thus, for a proper interpretation within the framework of uniform semi-normalization, their transformational parts must and can be included in uniform symmetry groups for the corresponding values of . Each transformation with is the commutator of such transformational parts with zero for any value of and thus necessarily belong to any possible uniform symmetry groups of the class . Moreover, the above transformational parts involve antiderivatives of and of , which are defined up to constant summands, and thus it is also natural to include the transformations with positive real into possible uniform symmetry groups of the class . Therefore, the family is the only natural choice for a family of uniform symmetry groups of the class . Since each of the transformations belongs to the kernel symmetry group of the class , which can be embedded in the equivalence group via trivial prolongation of its transformations to the arbitrary element [17, Proposition 3], the intersection of with uniform symmetry groups cannot coincide with . This means that the class is not disjointedly semi-normalized with respect to the entire equivalence group although it is semi-normalized with respect to and the family . For any we have .
Remark 10.
It is also clear that each disjointedly semi-normalized class is semi-normalized. At the same time, there are semi-normalized classes that are not disjointedly semi-normalized, as can be seen from the example of the class of nonlinear diffusion equations of the form with briefly discussed in a different terminology in [35, Section 3]. Let us properly present this example with extended arguments. The equivalence group of the class consists of the transformations
where , , , , and are arbitrary real constants with , see, e.g., [1, Eq. (3.15)]. For our purpose, it is more convenient to introduce and reparameterize the group as
(3) |
where , , , , and are arbitrary real constants with . The kernel point symmetry group of the class , where denotes the point symmetry group of the equation , is constituted by the point transformations
where , and are arbitrary constants with . The point symmetry group of the nonlinear diffusion equation with consists of by the point transformations
(4) |
where , , , , and are arbitrary constants with and , and, according to [32, 33], we take the natural domain for each transformation of the form (4) and use the modified composition of such transformations that is based on completing by continuity. The group is isomorphic to , where is the subgroup of constituted by all invertible orientation-preserving affine transformations of the affine space , .
Suppose that the class is disjointedly semi-normalized with respect to a subgroup of and a symmetry-subgroup family . The subgroup of should contain a transformation with a nonzero value of the parameter . Since in addition is a normal subgroup of , it then contains the subgroup of elements of with , and . If contains an element of with , then it contains the subgroup of shifts with respect to , with running through . Then the subgroup of singled out be the constraint consists of the transformations , , , where the parameters and freely run through , the parameter either freely runs through or is equal to zero and the parameter runs through a subgroup of the multiplicative group of positive real numbers. Denoting and , we reparameterize elements of the subgroup as , , , where the parameter tuple satisfies the same conditions as the parameter tuple does. In particular, the parameter runs through a subgroup of .
In the parameterization (3) of , it is clear that varying the parameters , and definitely leads to changing general values of the arbitrary element , and varying the other parameters does not affect . This is why the subgroup necessarily contain a subgroup of parameterized by , where the other parameters , and are assumed to be functions of , which are respectively equal to 0, 0, 1 at . Since for any , where denotes the scaling transformation , the elements of can be represented as
the parameter tuple freely runs through , the parameter runs through a subgroup of the multiplicative group of real numbers that is a complement to , , , and, if is free in , are functions of , which are respectively equal to 0, 0, 1 at , and is a free real parameter if takes only the zero values for elements of .
Thus, the ratio runs through the entire set of real numbers if and run through and , respectively. Setting , , , , , if is a free real parameter and otherwise, we obtain nonidentity transformations that are common for and , i.e., .
Theorem 5 implies the following assertions.
Corollary 11.
If the class is -semi-normalized and does not possess inessential equivalence transformations [13, Appendix A] and is a subgroup of that properly contains , then this class is non-disjointedly -semi-normalized.
Proof.
If the class is non-disjointedly -semi-normalized, then the corollary statement is trivial.
Suppose that the class is disjointedly - and -semi-normalized. Let be an element of . For any , consider admissible transformation . Since , we have for some and some . Thus, for some , for some , and . Then , i.e., and . In other words, the equivalence transformations and generate the same subset of and, therefore, should coincide in view of the absence of inessential equivalence transformation for the class , which contradicts the fact that and . ∎
A typical property for the intersections of with elements of the symmetry-subgroup family is given by . Moreover, any semi-normalized class, even a disjointedly semi-normalized one, can be considered as semi-normalized with respect to a wider subgroup of the corresponding equivalence group and a family of wider uniform point symmetry subgroups that satisfy the above property.
Corollary 12.
If a class is semi-normalized with respect to a subgroup of its equivalence group and a symmetry-subgroup family , then it is also semi-normalized with respect to the subgroup of and the symmetry-subgroup family , where for any .
Proof.
Denote by the subgroup of associated with , which is obtained from the group by the trivial (identity) extension of its elements to the arbitrary elements of the class . Since is the normal subgroup of [17, Proposition 3], the Frobenius product is a subgroup of , and then Theorem 5 directly implies the required claim. ∎
In the notation of Corollary 12, If , then the class is non-disjointedly -semi-normalized.
Recall that a normalized class is semi-normalized with respect to its entire equivalence group and the trivial symmetry-subgroup family . Hence Corollary 12 implies, in view of Theorem 4(ii), the following assertion, whose infinitesimal counterpart is Corollary 2 in [17].
Corollary 13.
For any system from a normalized class , the kernel symmetry group of the systems in is a normal subgroup of the point symmetry group of .
5 Uniform semi-normalization and
semi-normalization with respect to linear superposition
In addition to the basic property for any of Definition 2, the attribute “uniform” for the subgroups is also justified by the fact that in most of practically relevant examples of semi-normalized classes, all the subgroups are isomorphic or are at least of similar structure (in particular, of the same dimension). To distinguish such examples, where in addition the subgroups are known a priori, we call the corresponding class uniformly -semi-normalized. If in addition the class is disjointedly -semi-normalized, then we can call it uniformly disjointedly -semi-normalized.
Example 14.
The class defined in Remark 10 is semi-normalized but not uniformly semi-normalized. Indeed, the point symmetry group of the nonlinear diffusion equation with in particular contains the transformations of form (4), which do not belong to , and hence for any choice of a family of uniform symmetry subgroups, such a transformation belongs to . At the same time, the point symmetry group of any equation from the complement of the -orbit of in is contained in . Therefore, the structure of the corresponding element of definitely differs from that of .
A particular but important case of uniform semi-normalized classes, which is also relevant to the linear Schrödinger equations studied in this paper, is given by classes of homogeneous linear systems of differential equations. Below, we follow the presentation in the second part of Section 3 in [35] and modify it according to the framework developed in Sections 3 and 4.
We start with a normalized superclass of (in general) inhomogeneous linear systems of differential equations of the form : with .777We avoid the parameterization of by derivatives of arbitrary elements, considering such involved derivatives as additional arbitrary elements that are related to the original arbitrary elements via additional auxiliary equations as their derivatives. Here we consider the arbitrary-element tuple as consisting of two parts and change its notation from to . The subtuple parameterizes coefficients of the homogeneous linear left-hand side and depends at most on . Each component of the right-hand side runs through the set of sufficiently smooth functions of . Suppose that the class also satisfies the following conditions:
-
1.
Each system from this class is locally solvable.
-
2.
Elements of the pushforward of the equivalence group of the class by are fibre-preserving point transformations of whose components for are affine in .
The second condition means that elements of are of the form
(5) |
The functions and may satisfy additional constraints or even be quite specific but the components of the tuple are arbitrary smooth functions of . The transformations
constitute a normal subgroup of . Furthermore, splits over since , where is the subgroup of the group defined by the constraint . The normalization of the class means that its equivalence groupoid coincides with the action groupoid of .
On the infinitesimal level, the counterpart of the second condition is that the pushforward of the equivalence algebra of the class by consists of the vector fields of the form
where coefficients and may satisfy additional constraints but the are arbitrary smooth functions of . The equivalence algebra is the semidirect sum of the subalgebra and the ideal , which are singled out from by the constraints and , respectively.
Any system from the class is mapped to the associated homogeneous system by the transformation from , where is a solution of . In other words, the class is mapped onto the corresponding class of homogeneous systems by a wide family of admissible transformations from the action groupoid of .
We will use a double interpretation of the class . On the one hand, it can be embedded in the class as the subclass singled out by the constraint , and then all the transformational structures related to are substructures of their counterparts for . In particular, under this interpretation, the equivalence groupoid of is the subgroupoid of that is singled out by the constraints , and for the source and target arbitrary elements and the parameter function tuple . Up to factoring out the insignificant equivalence transformations related to the constraint [13, Appendix A], we identify the equivalence group of the class with its subgroup consisting of the elements of for which the corresponding values of the parameter function tuple belong to the set of common solutions of the systems from . In other words, , where and the normal subgroup of is constituted by the transformations with . Analogously, , where is the equivalence algebra of , is a subalgebra of and is the ideal of consisting of the vector fields with . On the other hand, we can consider as the class of systems of the form , where is the complete tuple of its arbitrary elements that runs through the set . By we denote the natural projection from the space with the coordinates onto the space with the coordinates . Then the corresponding transformational structures are the pushforwards, by , of their counterparts under the former interpretation. An advantage of the latter interpretation is that it avoids the insignificant equivalence transformations in related to the constraint .
The systems and are -equivalent if and only if their homogeneous counterparts and are -equivalent. Here the - and -equivalences can be replaced by the - and -equivalences, respectively. Thus, the group classification of systems from the class reduces to the group classification of systems from the class .
For each we denote by the subgroup of the point symmetry group of that consists of the linear superposition transformations : , , where the tuple runs through the solution set of . The subgroup of and the family satisfy the conditions in Definition 7. Therefore, the class is disjointedly uniformly -semi-normalized. We call this kind of semi-normalization, which is characteristic for classes of homogeneous linear systems of differential equations, uniform semi-normalization with respect to the linear superposition of solutions. If the systems from the class have no common solutions, which is a regular situation, then , and we will not indicate in this case. By Theorem 8(ii), for each the group splits over , and , where . By Theorem 8(iii), the splitting of the point symmetry group induces a splitting of the corresponding maximal Lie invariance algebra, . Here is the essential Lie invariance algebra of , , and the ideal , being the trivial part of , consists of vector fields generating one-parameter symmetry groups of linear superposition of solutions. Thus, the group classification problem for the class reduces to the classification of appropriate subalgebras of the equivalence algebra of this class. The qualification “appropriate” means that the pushforward of such a subalgebra by is the essential Lie invariance algebra of a system from . We can classify appropriate subalgebras of instead of , see Definition 16 below.
Given a class of linear homogeneous systems of differential equations that is uniformly semi-normalized with respect to the linear superposition of solutions, it is not necessary to start by considering the associated normalized superclass of generally inhomogeneous linear systems. The class itself can be the starting point of the analysis. In order to get directly its specific semi-normalization, we need to suppose the following properties of :
-
1.
The transformational part of any admissible transformation in the class is of the form (5). (Then the parameter function tuple necessarily runs through the solution set of the source system.)
-
2.
Any admissible transformation in the class with zero value of belongs to the action groupoid of the equivalence group of this class.
The normal subgroup and the subgroup of and the family are defined as above, and the class is disjointedly uniformly -semi-normalized, i.e., it is uniformly semi-normalized with respect to the linear superposition of solutions. The action groupoid of can be called the essential equivalence groupoid of the class [68, Remark 16].
The class of linear Schrödinger equations fits well into the developed framework, which we use in the present paper to study the group classification problem for this class in general and completely solve it in dimension 1+2 in Sections 7 and 8, respectively.
Remark 15.
There exist classes of homogeneous linear systems of differential equations that are not uniformly semi-normalized with respect to the linear superposition of solutions but disjointedly uniformly semi-normalized with respect to families of point-symmetry subgroups that are wider than the corresponding point-symmetry subgroups associated with the linear superposition of solutions, see Remark 9.
6 Algebraic method of group classification
for semi-normalized classes
The property of -semi-normalization of the class with known and allows one to solve the group classification problem for this class using an algebraic approach.888See [59, Section 2.2] for a modern revisited statement of group classification problems for classes of differential equations.
The subgroups of the groups and the subalgebras of the algebras , where runs through , are in general not known after describing the ingredients of the semi-normalization structure of . According to Theorem 6, we have that for each , the maximal Lie invariance algebra of the system is the sum of and , . Note that in general, the essential Lie invariance algebra nontrivially intersects the uniform Lie invariance algebra , and hence a natural basis of the maximal Lie invariance algebra is constituted by a basis of the intersection and its complements to bases of and of . (In the infinite-dimensional case, bases are replaced by spanning sets of vector fields.) The essential Lie invariance algebras are subalgebras of but usually, only some of the subalgebras of can serve as such algebras.
Definition 16.
We call a subalgebra of (resp. of ) appropriate if it (resp. its pushforward by ) coincides with the essential Lie invariance algebra of a system from the class .
The pushforwards by transformations from (resp. from ) preserve the set of appropriate subalgebras of (resp. of ) since for any and any ,999. The second, the third and the fourth equalities hold since the mapping is a bijection, all elements of and of are projectable to the space with the coordinates and , respectively. We also have that . which generates an equivalence relation within this set. It is natural to call it the -equivalence of appropriate subalgebras. In view of Theorem 4(iii), the -equivalence can be replaced by the - and even the -equivalence.
Making a preliminary analysis of the determining equations for Lie-symmetry vector fields of systems from , one can find constraints that are satisfied by appropriate subalgebras of but not by general ones. These constraints often include upper bounds for the dimensions of the appropriate subalgebras or of their specific parts. This crucially reduces the amount of subalgebras to be classified and makes the classification problem tractable even in the case of infinite-dimensional . Instead of classifying the appropriate subalgebras of , one can classify the appropriate subalgebras of and then push them forward by .
Consequently, we derive the following assertion.
Proposition 17.
Let the class be semi-normalized with respect to a subgroup of and a symmetry-subgroup family that are known, and let be a subalgebra of that corresponds to . Then the solution of the group classification problem for the class reduces to the classification, up to the -equivalence, of appropriate subalgebras of or, equivalently, of the algebra itself.
In view of Theorem 4(iii), Theorem 6 and the above discussion, we suggest the following procedure of solving the group classification problem for a semi-normalized class of systems of differential equations. We formulate this procedure in terms of the algebra .
-
•
When computing the equivalence groupoid and analyzing its structure, construct a subgroup of and a family of uniform point symmetry groups such that the class is -semi-normalized.
-
•
Find the subalgebra of associated with and the uniform Lie invariance algebras associated with .
-
•
Make a preliminary analysis of the system of determining equations for Lie-symmetry vector fields of systems from the class . The system usually splits into two subsystems, the “essential” and the “uniform” ones, which are satisfied by the components of vector fields from the respective parts of the maximal Lie invariance algebras of systems from the class .
-
•
Derive constraints satisfied by the appropriate subalgebras of .
-
•
Classify, modulo the -equivalence, the subalgebras of that satisfy the derived constraints.
-
•
For each subalgebra in the constructed set of -inequivalent families of subalgebras of ,
-
–
successively substitute the components of its basis/spanning elements into ,
-
–
merge all the obtained systems into a single system ,
-
–
solve the system with respect to the arbitrary elements of the class ,
-
–
for each solution of , check whether (it may happen that ).
-
–
-
•
The subalgebra is appropriate if and only if the system is compatible and for a solution of . This allows one to finally select those subalgebras in that are indeed appropriate subalgebras of .
-
•
For each subalgebra in the constructed set of -inequivalent families of appropriate subalgebras of , analyze the equivalence within the set that is induced by the stabilizer of in . If there are -equivalent subalgebras in , one can also consider the equivalence between the corresponding sets . When possible, use these equivalences to reduce the sets and .
-
•
Assemble the set .
The last set constitutes the solution of the group classification problem for the class in terms of essential Lie invariance algebras modulo both the - and the -equivalences [59, Section 2.2]. We call the above procedure the algebraic method of group classification for semi-normalized classes of differential equations. To present the solution of the group classification problem for the class in terms of maximal Lie invariance algebras, one needs to indicate the algebra for each and modify the set of classification cases to .
In fact, the above method is effective for the case of uniform semi-normalization, especially disjointed uniform semi-normalization, including uniform semi-normalized with respect to the linear superposition of solutions. Then all the listed steps of the classification procedure can be properly implemented, and the mere presentation of inequivalent essential Lie-symmetry extensions is definitely sufficient.
It may happen that the same class of systems of differential equations is semi-normalized under different choices of subgroups of its equivalence group and of families of uniform symmetry-subgroups of its systems. The proper selection of these objects helps to solve the group classification problem for such a class in a more efficient way.
Remark 18.
The class from Remark 9 is non-disjointedly semi-normalized with respect to the entire group and the family . At the same time, the interpretation of the class as disjointedly semi-normalized with respect to the family and the proper subgroup of singled out from by the constraint , see Remark 9, is more convenient and effective for group classification of this class. Indeed, then one can apply Theorem 8, which is slightly stronger than Theorems 4 and 6, and, following Proposition 17, handle the smaller group . One can also construct the generalized extended equivalence group of the class or its efficient part [57, 58], considering the reparameterized class for , where additional arbitrary elements and defined by the constraints and , and prove that the reparameterized class is uniformly semi-normalized with respect to the linear superposition of solutions in the generalized sense, i.e., the original class is uniformly semi-normalized with respect to the linear superposition of solutions in the generalized extended sense. Nevertheless, the group classification of based on the last interpretation is essentially less efficient computationally than those based on the above semi-normalizations of in the usual sense, not to mention the absence of the corresponding theoretical background and related techniques for the case of semi-normalization in the generalized extended sense.
Remark 19.
There are examples of the opposite pattern. Given a class of systems of differential equations that is disjointedly uniformly semi-normalized with respect to a subgroup of its equivalence group and a family of uniform symmetry-subgroups of its systems, nevertheless, the group classification problem for this class can be more efficiently solved using its non-disjointed semi-normalization with respect to a narrower subgroup of and a family of wider uniform symmetry-subgroups rather than involving the former (disjointed) semi-normalization. In particular, this is the case for the class of (1+2)-dimensional linear Schrödinger equations of the form (1) with potentials , where is a complex-valued smooth function of and is a real-valued smooth function of , see the consideration related to the equation (20) below. The class is (disjointedly) uniformly semi-normalized with respect to its entire equivalence group and linear superposition of solutions. At the same time, it is more convenient to solve the group classification problem for this class using its non-disjointed semi-normalization with respect to a proper subgroup of and a family of wider uniform symmetry groups than the symmetry groups of linear superposition of solutions. This is the way of treating the case in the proof of Theorem 38 below.
7 Preliminary symmetry analysis of multidimensional
linear Schrödinger equations with complex-valued potentials
To demonstrate the efficiency of the developed theory and the proposed algebraic method of group classification for semi-normalized classes, we study the class of (1+)-dimensional () linear Schrödinger equations with complex-valued potentials, which are of the form (1). First, in this section we carry out a preliminary analysis of transformational properties and Lie symmetries of equations from this class for an arbitrary number of space variables . It turns out that it is possible, without fixing a value of , to find the equivalence groupoid , the equivalence group and the equivalence algebra of the class as well as to preliminarily analyze the determining equations for Lie symmetries of equations from this class and to derive basic properties of such symmetries. In Section 8, we specify these results to the case and use them to carry out the complete group classification of the class with . Recall that the case was exhaustively studied in [35].
Similarly to other Schrödinger equations but in contrast to most systems of differential equations arising in applications, for equations from the class , the independent variables , , …, take values in the set of real numbers, whereas values of the dependent variable and of the arbitrary element are in general complex. To avoid considering complex values, one may replace Schrödinger equations by equivalent systems for two real-valued unknown functions, e.g., the real and the imaginary parts of or the absolute value and the argument of . Nevertheless, such a replacement complicates the study more essentially than the involvement of complex numbers does, especially when using the absolute value and the argument of , which results in nonlinear systems instead of linear ones. This is why we work with the form (1) of equations from the class . In view of the complex nature of and , it is necessary to formally extend the space of variables with and the tuple of the single arbitrary element with . Here and in what follows star denotes the complex conjugation. In particular, we introduce (resp. ) as an additional argument for all functions depending on (resp. ), which includes, e.g., the components of point transformations and of vector fields in the spaces with coordinates or . We thus follow Wirtinger’s approach.101010The variables and with are in fact not independent of each other in the canonical sense since they are uniquely related by the complex conjugation. However, the equalities and with the Wirtinger derivatives and allow one to formally treat and as independent variables and to simultaneously indicate them as arguments of relevant functions instead of . Confining a differential function of , i.e., of in the above interpretation, to the solution set of an equation from the class , we should also take into account its complex conjugate , thus considering the system of two equations instead of the single one, cf. [35, Section 1]. At the same time, there two essential simplifications due to using Wirtinger’s approach. In view of the fact that the invariance (resp. equivalence) conditions hold or do not hold simultaneously for equations from the class and their complex conjugate counterparts, it suffices to check such conditions only for the former equations. Moreover, presenting point transformations in the spaces with coordinates or , we can omit the transformation components for and since they are just the complex conjugates of those for and , respectively.
Notation.
In this and the next sections, and are the real independent variables, is the unknown complex-valued function of and . For a complex value , its conjugate is denoted by , and we define
where is the -component of point transformations in the space with as a coordinate. The indices and run from to , , , and . is the identity matrix. The total derivative operators and are defined as
7.1 Equivalence groupoid
We find the equivalence groupoid of the class using the direct method, see, e.g., [57, Section 2] and [68, Section 2] for a general description of this method and [7, 30, 31] for computational techniques involved. Let denote the equation from with potential . We seek for all point transformations of the form
(6) |
where with , that map a fixed equation from the class to an equation : from the same class.
The following assertion is proved similarly to [67, Lemma 1].
Lemma 20.
Any point transformation connecting two equations from the class satisfies the conditions
Theorem 21.
The equivalence groupoid of the class consists of the triples of the form , where is a point transformation in the space of variables, whose components are
(7a) | |||
(7b) | |||
and the target potential is expressed via the source potential as | |||
(7c) |
Here , , and are arbitrary smooth real-valued functions of with , , is an arbitrary constant orthogonal matrix, and is an arbitrary solution of the initial equation.
Proof.
Let a point transformation connect equations and from the class . Lemma 20 implies that , and with . Acting the total derivative operators and on the equality , we derive
The above equations are equivalent to
(8) |
where , and is the Kronecker delta. In fact, the vector-function is the inverse of the vector-function with respect to , and . We substitute the expressions (6) and (8) for , and and then the expression for , , into the equation . As a result, we derive the equation
Then splitting this equation with respect to various derivatives of and additionally arranging lead to the system
(9) | ||||
(10) | ||||
(11) |
The first equation in (9) together with the condition imply that . Therefore, , i.e., in terms of the tuple . We differentiate this equation for each with respect to and permute the indices , and , which gives the equations
implying for all values of . Since for each the tuples , …, constitute an (orthogonal) basis of , this means that , i.e., the vector-function is affine in with coefficients depending on . The equations constrain the matrix , and thus we have where is a time-dependent orthogonal matrix and is a time-dependent vector.
The general solution of the second equation in (9) is , where and are smooth complex-valued functions of and with .
We substitute the above expression of and into (10) to get
(12) |
Cross-differentiating copies of the equation (12) with different values of and subtracting the results of differentiation from each other give . Hence , or in the matrix terms, where is the transpose of . At the same time, differentiating the orthogonality condition , where is the identity matrix, with respect to results in . Therefore, , from which we have , and thus is a constant orthogonal matrix.
The system consisting of the equations (12) when varying integrates to the following expression for :
(13) |
where and are arbitrary smooth real-valued functions of . Finally, we consider the equation (11), which reduces, under the derived conditions, to the equation
Splitting this equation with respect to in view of the representation for and rearranging, we obtain
(14) | |||
(15) |
We introduce the function , i.e., . The equation (15) is equivalent to the initial linear Schrödinger equation in terms of . After substituting the expression (13) for into (14) then additionally collecting coefficients of , we derive the final expression for . ∎
Corollary 22.
A (1+n)-dimensional linear Schrödinger equation of the form (1) is reduced to the free linear Schrödinger equation by a point transformation if and only if
for some real-valued smooth functions , , and of .
7.2 Equivalence group and equivalence algebra
The concept of equivalence group, which was suggested in [61] and properly formalized in [67, Section 2.3], see also [68] and references therein, plays a central role in the theory of group classification of differential equations. The equivalence algebra and the equivalence group of a class of differential equations can be independently computed using the infinitesimal method and the direct or the algebraic method, respectively. Nevertheless, it is much easier to find them from the equivalence groupoid of the class once this groupoid is known. The following assertions follow from Theorem 21 in a similar way as that we used in [35] to obtain consequences of [35, Theorem 6].
Corollary 23.
The (usual) equivalence group of the class consists of the point transformations in the space with the coordinates whose components are of the form (7) with .111111Strictly following the definition of (usual) equivalence groups, we should assume that elements of the group are point transformations in the space with the coordinates . Here the subscript “” of denotes the collection of jet coordinates corresponding to the derivatives of up to order two, including itself as its zeroth-order derivative, , and similarly for . Nevertheless, the arbitrary element depends only and the relation (7c) between the source and target values of for admissible transformations within the class does not involve nonzero-order derivatives of . Therefore, analogously to [36, footnote 1] and [35, Section 1], we can assume that elements of the group act in the space with the coordinates , thus shrinking the space underlying this group. It suffices to present only the transformation components for . The transformation components for derivatives of are constructed from the -components by the standard prolongation using the chain rule. The transformation components for derivatives of and for are obtained by conjugating their counterparts for derivatives of and for .
Remark 24.
The identity component of consists of the transformations of the form (7), where , and , i.e., . The whole equivalence group is generated by the transformations from its identity component and two discrete transformations, the space reflection , for a fixed and the Wigner time reflection , .
Theorem 25.
The equivalence algebra of the class is the algebra
where , , and run through the set of smooth real-valued functions of . The vector fields spanning are defined by
Proof.
The equivalence algebra of the class is obtained using the knowledge of the identity component of the equivalence group . Searching for infinitesimal generators of one-parameter subgroups of , we represent the parameter function in the form , where is a function of , for better consistency of the group parameterization with the one-parameter subgroup structure of . Then we successively assume one of the transformation parameters , , , and to depend on a continuous parameter and set the other transformation parameters to the trivial values corresponding to the identity transformation, which are for , for and zeroes for , and , where is the identity matrix. The components of the associated infinitesimal generator are given by
The above procedure results in the set of vector fields spanning the algebra . ∎
Corollary 26.
The class is uniformly semi-normalized with respect to the linear superposition of solutions.
The uniform semi-normalization of the class guarantees a specific factorization of point symmetry groups of all equations from this class. For any potential , each element of the point symmetry group of an equation generates an admissible point transformation of . Therefore, the transformation necessarily has the form (7a)–(7b), where the transformation parameters additionally satisfy the equation (7c) with . The symmetry transformations associated with the linear superposition of solutions to the equation which are of the above form with , and , constitute a normal subgroup of the group , which is denoted by and is the trivial part of . In view of the discussion on the uniform semi-normalization with respect to the linear superposition of solutions in [35, Section 3], Corollary 26 implies that the group splits over , , the subgroup of is singled out from by the constraint and will be considered as the only essential part of .
7.3 Analysis of determining equations for Lie symmetries
Using the infinitesimal criterion for Lie symmetries [10, 56, 61], for each potential we can derive the determining equations which are satisfied by the components of vector fields from the maximal Lie invariance algebra of the equation from the class . In this section, we partially integrate the determining equations and preliminary analyze properties of the collection of when varies for .
The infinitesimal criterion states that a vector field , where the components , and are smooth functions of , and is the complex conjugate of , belongs to the algebra if an only if
with being the second prolongation of the vector field , see Section 2. Expanding this expression, we obtain
(16) |
where and . We recall that and are the operators of total derivatives with respect to and , respectively. The substitution of and into (16) and splitting with respect to the various derivatives of and lead to a linear overdetermined system of determining equations for the coefficients of ,
(17a) | |||
(17b) | |||
(17c) |
The subsystem consisting of the equations (17a) and (17b), which do not contain the potential , is integrated immediately. Its general solution is
where , , and are smooth real-valued functions of , is a complex-valued function of and , and is a constant skew-symmetric matrix. Then additionally splitting the last determining equation with respect to , we derive two equations,
(18) | |||
(19) |
Both the equations involve the potential . At the same time, the first equation means that the parameter function is an arbitrary solution of the equation . This is why the second equation is only the truly classifying condition for Lie symmetry generators of equations from the class depending on the potential .
Theorem 27.
The maximal Lie invariance algebra of an equation from is constituted by the vector fields of the form , where
the parameters , , and are arbitrary real-valued smooth functions of and the matrix is an arbitrary constant skew-symmetric matrix that jointly satisfy the classifying condition (19), and the parameter runs through the solution set of the equation .
The kernel invariance algebra of the class , , is obtained by splitting the conditions (18) and (19) with respect to the potential and its derivatives. This results in , , and .
Proposition 28.
The kernel Lie invariance algebra of the class is .
Following [35], denote by the linear span of all vector fields given in Theorem 27 when the potential varies, i.e.,
Here and in what follows the parameters , , and run through the set of real-valued smooth functions of , runs through the set of complex-valued smooth functions of and runs through the solution set of the equation when the potential is fixed. We have since each vector field from either with nonvanishing or or or with jointly vanishing , , , and necessarily belongs to for some . Up to the antisymmetry of the Lie bracket of vector fields, the nonzero commutation relations between vector fields spanning are
The commutation relations between elements of imply that the span is closed with respect to the Lie bracket of vector fields and hence it is a Lie algebra. The algebra can be represented as a semidirect sum of the subalgebra and the ideal , . The kernel invariance algebra is an ideal in and in the entire . The above representation for induces the similar representation for each ,
are a finite-dimensional subalgebra (see Lemma 31 below) and an infinite-dimensional abelian ideal of , respectively. We call the essential Lie invariance algebra of the corresponding equation . The ideal consists of vector fields associated with transformations of linear superposition on the solution set of the equation and therefore it is a trivial part of .
Definition 29.
A subalgebra of is called appropriate if there exists a potential such that .
The relation of the equivalence algebra and is given by , where is the mapping from onto that is induced by the projection of the joint space of the variables and the arbitrary element onto the space of the variables only. Although the projection is not a (local) diffeomorphism, it properly pushes forward vector fields from due to their structure and hence the mapping is well defined. It maps the vector fields , , , , spanning to the vector fields , , , , spanning , respectively. The above relation is stronger than the inclusion implied by the specific uniform semi-normalization of the class .
As the algebra coincides with the set of infinitesimal generators of one-parameter subgroups of the group which is the projection of by , the structure of is agreed with the action of on this algebra. Moreover, both and are invariant with respect to the action of the group . This is why the action of on equations from the class induces the well-defined action of on the essential Lie invariance algebras of these equations, which are subalgebras of . The kernel is obviously an ideal in for any .
Summarizing altogether we have the following proposition:
Proposition 30.
The problem of group classification of (1+n)-dimensional linear Schrödinger equations with complex-valued potentials of the form (1) reduces to the classification of appropriate subalgebras of the algebra with respect to the equivalence relation generated by the action of .
7.4 Further properties of Lie invariance algebras
In order to classify all appropriate subalgebras of , we start with describing the action of the transformations from on the vector fields from . For this purpose, we compute the pushforwards of spanning vector fields of by elementary transformations from .
Given and , the pushforward by is computed as
where into each coefficient of we substitute the expressions of the variables without tildes in terms of the variables with tildes, , and denotes the inverse of .
By , , , and we denote the transformations of the form (7a)–(7b) with and , where the parameters , , , and , successively excluding one of them, are set to the trivial values, which are for , for and zeroes for , and . These transformations, referred to as elementary transformations, generate the entire group . The nonidentity pushforward actions of elementary transformations to the vector fields spanning are exhausted by
where , , , , , , and into each pushforward by we should substitute the expression for given by inverting the relation ; for the other pushforwards. Tildes over vector fields mean that these vector fields are represented in the new variables.
Lemma 31.
for any potential , and this upper bound is least.
Proof.
Fix a potential . Then similarly to the proof of Lemma 1 in [35], the classifying condition (19) implies a system of linear ordinary differential equations of the normal form
where the coefficients and , , , , are functions of . Therefore, the dimension of is bounded by the sum of the number of pairs with (which is associated with rotations and equals ) and the number of arbitrary constants in the general solution of the above system (which equals ). The derived upper bound is least since it coincides with for the zero corresponding to the free Schrödinger equation. ∎
Corollary 32.
.
Proof.
We follow the proof of Lemma 31 by omitting the first equation of the system from this proof and setting and . ∎
Corollary 33.
for any potential .
Proof.
Analogously, we preserve only the last two equations of the system from the proof of Lemma 31 and set and . ∎
Lemma 34.
For all , is a Lie algebra and . Moreover,
where denotes the projection on the space of the variable and .
Proof.
The proof is similar to the one given in [35]. It is based on the famous Lie classification of finite dimensional algebras of vector fields on the real line. ∎
8 Complete group classification of (1+2)-dimensional
linear Schrödinger equations with complex-valued potentials
We complete the group classification of linear Schrödinger equations of the general form (1) only in the case . Even the case is much more cumbersome and requires more study, not to mention the case of general . Specifying results of Section 7.3 to the case , we obtain that for any
For each appropriate subalgebra of , similarly to [35] we introduce five -invariant integers in order to arrange the process of the group classification:
where denotes the projection on the space of the variable . Note that
The parameter tuple is used for the labeling of possible classification cases. In fact, not all values of in are appropriate.
Lemma 35.
.
Proof.
Suppose that and . Then the algebra necessarily contains vector fields of the form and , where . The commutator of these vector fields
should belong to , and hence . Therefore, the case is impossible. ∎
The other constraint for appropriate values of , if , is not obvious at all and can be derived only by analyzing the complete proof of group classification for the class with .
Lemma 36.
Up to the -equivalence, a vector field , where the components of the tuple are linearly independent, reduces to with nonzero .
Proof.
We can represent the tuple as , where and are smooth functions of with . Then we successively apply the transformations and with an appropriate tuple of smooth functions of to the vector field and thus set and , which gives the required form. ∎
Lemma 37.
If , then the vector field with also belongs to .
Proof.
Putting all the above ingredients together, we have the following: any appropriate subalgebra of is spanned by
-
•
the basis vector fields and of the kernel ,
-
•
vector fields , , with linearly independent tuples , …, ,
-
•
vector fields ,
-
•
vector fields , , with linearly independent , …, .
In what follows we use the notation
Theorem 38.
A complete list of inequivalent essential Lie-symmetry extensions in the class of (1+2)-dimensional linear Schrödinger equations with complex-valued potentials is exhausted by the cases presented below. Here is an arbitrary complex-valued smooth function of its arguments or an arbitrary complex constant that satisfies constraints indicated in the corresponding cases, and the other functions and constants take real values.
-
0.
: .
-
1.
: .
-
2.
: .
-
3.
, , , :
-
4.
, : .
-
5.
, and, if , :
. -
6.
: .
-
7.
: .
-
8.
:
. -
9.
, :
. -
10.
, :
. -
11.
, :
. -
12.
, : .
-
13.
, :
. -
14.
, :
. -
15.
: .
-
16.
, :
where is a fundamental set of solutions of the system , and . -
17.
, or :
, where is a fundamental set of solutions of the equation , , , and is a fundamental set of solutions of the equation , , . -
18.
, or :
, where are linearly independent solutions of the system , , and , . -
19.
: .
Remark 39.
The essential Lie invariance algebras listed in Theorem 38 are really maximal for the corresponding potentials if these potentials are -inequivalent to listed potentials with larger essential Lie invariance algebras. In some cases, we have presented simple necessary and sufficient conditions that provide such inequivalence. In other cases, these conditions are not so obvious. For example, in Cases 9–11 the condition of maximal Lie symmetry extension is ( or ) and ( or ), which excludes the values of that are -equivalent to those in Cases 12 and 16–19. Case 8 is already more complicated since for it the similar condition ( or ) and ( or ) should be extended with the condition excluding potentials -equivalent to those in Cases 9–11. Analogously, potentials in Cases 13–15 are -inequivalent to those in Cases 16–19 if and only if or . This condition is necessary and sufficient for the maximality of Lie symmetry extensions given in Cases 14 and 15, and for Case 13 it should be extended to guarantee the exclusion of potentials related to Cases 14 and 15. Schrödinger equations related to Case 16 are not similar to the free Schrödinger equation presented in Case 19 if and only if the parameter functions satisfy at least one of the conditions , , , and . We need additional constraints on these functions to exclude, up to the -equivalence, potentials of Cases 17 and 18.
Remark 40.
In Theorem 38, we in general neglect possible gauges of constant parameters in by discrete and scaling equivalence transformations. For example, a nonzero parameter among and in Case 6, the parameter in Case 9 and the parameter in Case 17 can be set to one. The parameter in Cases 10 and 11 and the imaginary part of in Cases 7 and 12 can be made nonnegative. In Case 18 and, if , in Case 17, we can also assume .
Remark 41.
More generally, the cases that are listed in Theorem 38 are not equivalent to each other. At the same time, each of Cases 0–18 is a family of cases of essential Lie symmetry extensions within the class rather than a single case, and there is an equivalence within each of these case families that is induced by the action of a subgroup of . For example, see below the consideration of the subclass of (1+2)-dimensional linear Schrödinger equations with potentials of the form (20), in particular the formula (21), for the equivalence of potentials in Case 8. Nevertheless, it is usually not convenient or even not feasible to take such an equivalence into account.
Proof.
We single out different cases depending on values of the integers , and . Basis vector fields of a Lie symmetry extension is denoted by
where the range of is equal to , are real constants, and all the other parameters are real-valued functions of .
. The corresponding Lie symmetry extensions are spanned by with linearly independent . By Lemma 34, is an algebra isomorphic to a subalgebra of the algebra . Hence the group classification in this case reduces to the classification of subalgebras of the algebra . Varying we single out the following four subcases.
This is the general case with no extension, i.e., (Case 0).
. The algebra contains a single linearly independent vector field with . Pushing forward by a transformation from , we can set , and , i.e., the vector field reduces to the form . Evaluating the classifying condition (19) at with and gives differential equations in whose general solutions are presented in Cases 1 and 2, respectively.
The Lie symmetry extension is given by vector fields and with linearly independent and . Modulo the -equivalence, we can set and (cf. Lemma 34) and, as in the previous case, . Then, the condition implies that
where and are real constants. This equality splits into the constraints , , , . Taking these constraints into account, we push and forward by transformations from and linearly combine them with and . This allows us to set and when preserving . The integration of the system in obtained by evaluating the classifying condition (19) at and at leads to Case 3. For the extension to be maximal, we need (otherwise we get Case 4 with ) and (otherwise we get Case 7).
The algebra is spanned by the basis elements of the kernel algebra and vector fields with linearly independent coefficients , . In view of Lemma 34, we can assume that , and . The commutation relations of imply . Similarly to the previous case, we can also set , . The commutation relations of with and implies , and , where and are real constants. Up to linearly combining with and , the vector field reduces to . The classifying condition (19) evaluated successively at , and gives the system
whose compatibility with respect to implies and . Thus, we have Case 4, where the extension is maximal if and only if since otherwise it becomes Case 7.
The algebra contains a vector field with and , and additional extensions are realized by vector fields with linearly independent coefficients . The commutation relation for and implies that . Acting by with and on , we set , and thus . Similarly to the case , we partition Lie symmetry extensions into three subcases depending on values of . Moreover, the simplification of the parameter functions , and up to the -equivalence jointly with linearly combining completely coincides with that in the case .
. If the function (resp. ) is constant, then it can be made zero by linearly combining with (resp. ). If one of these functions is not constant, up to the -equivalence it can be set to be equal to . Then the classifying condition (19) implies an equation in whose general solution is presented in Case 5.
The vector field is reduced to up to the -equivalence. As the Lie algebra is closed with respect to Lie bracket of vector fields, we obtain
where and are real constants. Equating the components of the vector fields on both the sides of the last equality and then solving the obtained equations, we derive and , where and are real constants, which can be set to zero due to linearly combining with and . The vector field is reduced to . Successively substituting the components of and into the classifying condition (19) yields two independent equations in : , , whose general solution is presented in Case 6.
The algebra necessarily contains the vector fields and from and vector fields and , , with linearly independent and . Up to the -equivalence the vector fields and are reduced to the form and . As the commutators belong to , we have
where and , , are real constants. The above commutation relations give the system , , , , which implies that and are constants. Therefore, by linearly combining with and , they can be set to be equal to zero. We successively evaluate the classifying condition (19) at , and and obtain the system , , , whose general solution is . At the same time, any such potential admits the Lie symmetry vector field . This means that in fact . For the Lie symmetry extension to be maximal, we require . Therefore, we obtain Case 7.
. We definitely have a vector field with a nonzero tuple in . An additional Lie symmetry extension may be provided by one more vector field of the same form or vector fields , , with linearly independent coefficients . We should suppose separately whether the tuple is either proportional or not proportional to a constant tuple.
1. Let the parameter functions and be linearly dependent. Up to the -equivalence, we can set , , and thus the vector field reduces to . Lemma 37 implies that then the algebra also contains the vector field with . Then substituting the components of the vector fields into the classifying condition (19) yields two dependent equations for , and , whose general solution is
(20) |
where is a complex-valued smooth function of and is a real-valued smooth function of . Thus, any equation from the class with potential of the form (20) is invariant with respect to the vector fields and . The function satisfies the conditions or since otherwise (see the case below).
Consider the subclass of (1+2)-dimensional linear Schrödinger equations with potentials of the form (20) constrained by the above condition for . We can reparameterize the class , assuming the parameter functions and as arbitrary elements instead of . The equivalence groupoid of can be singled out from the equivalence groupoid of the entire class , which is described in Theorem 21. A point transformation connects two equations and from the class if and only if it is of the form (7a)–(7b), where the function is linear fractional in , the matrix is diagonal with , and with arbitrary real constants and . The corresponding arbitrary-element tuples are related in the following way:
(21) |
The equivalence group of the class consists of the point transformations of the form (7a)–(7b) prolonged to the arbitrary elements according to (21), where the parameters satisfy, in addition to all the above conditions for elements of , the constraint . The class is (disjointedly) uniformly semi-normalized with respect to the entire equivalence group and linear superposition of solutions. At the same time, it is already known that each equation , where the potential is of the form (20), admits the one-parameter groups of point symmetry transformations
generated by the vector fields and . This is why it is convenient to include these symmetry transformations in the corresponding uniform symmetry group and to consider the class as (non-disjointedly and uniformly) semi-normalized with respect to a proper subgroup of and a family of wider uniform symmetry groups than the symmetry groups of linear superposition of solutions.
For the group to be closed with respect to the transformation composition, we also need to include the elementary transformations (cf. Section 7.4) with constant values of in since they are the commutators of transformations and . The elementary transformations with constant values of should be in since the expressions for the -components of the transformations and contain antiderivatives of and , which are defined up to constant summands, and there is no canonical choice of these antiderivatives for general values of . Thus, for each equation from the class , we consider the group of its point symmetry transformations of the form (7a)–(7b), where , coincides with the identity matrix, , , , , , and are arbitrary constants, and is an arbitrary solution of . The corresponding Lie algebra is spanned by the vector fields , , , and with running through the solution set of . The suitable subgroup of the group is singled out from this group by the constraint , and the Lie algebra associated with is .
Since all the conditions of Definition 2 are satisfied and the subgroups with running through the set of potentials of the form (20) are known and of the same structure, the class is (uniformly but not disjointedly) semi-normalized with respect to the subgroup of and the family , where . In view of Proposition 17, the group classification of the class reduces to the classification of -inequivalent appropriate subalgebras of the algebra
Denote . The condition and Corollary 33 imply that any Lie symmetry extension in the subclass is spanned by vector fields of the form , where runs through a set of linearly independent quadratic polynomial in . This shows that the computation of inequivalent Lie symmetry extensions reduces to the classification of subalgebras of the algebra . Depending on the values of we obtain the following subcases.
. There is no additional extension, i.e., we have Case 8.
The algebra necessarily contains a vector field with nonzero . Up to the -equivalence, the vector field takes one of the following values: , and . Taking into account the form (20) of , we successively evaluate the classifying condition (19) at each value of and obtain equations for and . The solution of these equations results in Cases 9, 10 and 11, respectively, where we change the notation of . Note that we choose the special form of the last value of for simplifying the representation of the corresponding potentials.
Apart from the vector fields and , the algebra additionally contains at least vector fields , , with linearly independent and . Up to the -equivalence and linearly combining with and , the vector fields reduce to and . Evaluating the classifying condition (19) simultaneously at these values of and for the expression (20) of , we obtain , and , and thus with a nonzero complex constant . At the same time, the equation with this possesses one more Lie symmetry vector field . This shows that the Lie symmetry extension with is not maximal. Therefore, we obtain Case 12.
2. Suppose now that the parameter functions and are linearly independent.
If we have no additional extension, then the algebra is spanned by the vector fields , and .121212In contrast to the previous case, where the parameter functions and are linearly dependent, in the case under consideration the algebra contains, up to linear dependence, only one vector field of the form as . Indeed, suppose that this algebra contains one more vector field , where the tuple is not linearly dependent with . Then the case condition implies that where is a nonconstant function of . Successively evaluating the classifying condition (19) at the vector fields and , we derive two equations for , for which the difference of the second equation and the first equation multiplied by leads, in view of the proportionality of the tuples, to the condition . Splitting it with respect to and integrating the obtained equations, we derive , where are real constants, which contradicts the linear independence of and . Lemma 36 implies that the vector field reduces to the form , where and are smooth functions of with . The substitution of the components of into the classifying condition (19) yields the equation
whose general solution is presented in Case 13.
Otherwise, the algebra also contains a vector field with nonzero , which takes, up to the -equivalence, the form , . The condition means that
where are real constants. Equating the corresponding components of the vector fields on both the sides, we derive the system
(22) |
If , then the parameter functions and are linearly dependent. Hence , i.e., . Integrating the systems (22) and linearly combining with elements from the kernel, we obtain a reduced form of depending on the value of .
If , then with real constants and . Pushing vector fields from forward by , we set . The successive substitution of the components of the vector fields and into the classifying condition (19) gives the system
whose general solution is presented in Case 14.
If , then the solution of the system (22) provides, after linearly combining with and , . In this case no further simplifications are possible. For convenience, we re-denote and . Evaluating the classifying condition (19) at and at yields two equations for ,
Therefore, we obtain Case 15.
Further classification. In view of Lemma 35, the values of not considered so far are and .
Since , the algebra contains at least two vector fields of the form , , where . The substitution of the components of the vector fields into the classifying condition (19) leads to the following system in :
This system can be written as , where the coefficients , and are real-valued functions of that satisfy the conditions
The condition implies that the matrix is symmetric. Thus, the potential is a quadratic polynomial in and with coefficients depending on ,
(23) |
where the functions , and can be set equal to zero up to the -equivalence but we will use this possibility later.
Next consider the subclass of equations from the class with potentials of the form (23). Theorem 21 implies that, similarly to the entire class , the subclass is uniformly semi-normalized with respect to the linear superposition of solutions, and its equivalence group coincides with the equivalence group of . The further group classification of the subclass splits into cases depending on the value of . The conditions and partition the subclass into the subclasses and , respectively. These subclasses are uniformly semi-normalized with respect to the linear superposition of solutions as well, and their equivalence groups also coincide with the equivalence group of the entire class . The classifying condition (19) implies that the subclass is singled out from the class by the conditions , , and , and the subclass is then singled out by requiring that one of these conditions fails to hold. In view of Corollary 22 the above condition means that each equation in the subclass is -equivalent to the (1+2)-dimensional free Schrödinger equation (with ). In other words, the subclass is the -orbit of the free Schrödinger equation, which gives Case 19.
Therefore, the last consideration concerns the subclass , for which and . Up to the -equivalence, the simplified form of potentials for equations from the subclass is
(24) |
Substituting of this form into the classifying condition (19) with and , splitting with respect to and separating real and imaginary parts, we derive the following system for :
(25) |
This system admits a fundamental set of solutions , , and , where the tuples are linearly independent, and . Thus, the algebra contains the four vector fields . The last two solutions correspond to the vector fields and from the kernel algebra . Case 16 presents general equations from the subclass with no additional Lie symmetry extension.
Otherwise, the algebra contains a vector field with nonzero . Up to the -equivalence, the parameter function is reduced to 1, , and the simplified form (24) of is preserved. We substitute the components of jointly with the form (24) of into the classifying condition (19) and split with respect to different powers of . As a result, we obtain that the tuple satisfies the system (25), and the parameter functions and are constrained by the equations
(26) |
Therefore, up to linearly combining with , and , we can set . This reduces to . Consider the cases and separately.
For , the system (26) implies that all and are constants. Up to rotations, we can reduce the matrix to a diagonal matrix . The maximality of the Lie symmetry extension requires or with , which gives Case 17.
If , then, up to translations of time, the general solution of the system (26) is , , , and , where , , and are arbitrary real constants. The Lie symmetry extension is maximal if these constants satisfy the inequalities or . Rearranging the potential in terms of leads to Case 18.
Let us show that the dimension of additional Lie symmetry extension cannot exceed one. Suppose that this is not the case. Then the algebra contains a two-dimensional subalgebra spanned by vector fields and with linearly independent and . Up to the -equivalence, we can assume that and . In a similar way as above, the vector field takes the form . The commutation relation for and implies that . Successively evaluating the classifying condition (19) at and at for the form (24) of and splitting with respect to different powers of provides the equations , which contradicts the condition singling out the subclass from the class . ∎
9 Subclass with real-valued potentials
Consider the subclass of the class , which consists of the equations of the form (1) with real-valued potentials , i.e., it is singled out from the class by the constraint . Since the class is uniformly semi-normalized with respect to the linear superposition of solutions, all results on symmetry analysis of its subclass can be easily derived from the corresponding results for this entire class via singling out objects that are consistent with the constraint . In particular, using Theorem 21, we construct the equivalence groupoid of the class as the maximal subgroupoid of with the set of objects consisting of the real-valued potentials .
Corollary 42.
The equivalence groupoid of the class consists of the triples of the form , where is a point transformation in the space of variables, whose components are of the form (7a)–(7b) and the target potential is expressed via the source potential as in (7c), where in addition we omit the hat accent over and
Corollary 43.
A (1+n)-dimensional linear Schrödinger equation of the form (1) with a real-valued potential is reduced to the free linear Schrödinger equation by a point transformation if and only if
for some real-valued smooth functions , and of .
Corollary 44.
The (usual) equivalence group of the class consists of point transformations in the space with the coordinates whose components are of the form (7), where in addition we omit the hat accent over and
Corollary 45.
The class is uniformly semi-normalized with respect to the linear superposition of solutions.
Remark 46.
The identity component of is constituted by the transformations from this group, where in addition and , i.e., . Similarly to , the entire equivalence group is generated by the transformations from its identity component and two discrete transformations, the space reflection , for a fixed and the Wigner time reflection , .
Corollary 47.
In the notation of Theorem 25, where we omit the -components in all vector fields, the equivalence algebra of the class is the algebra
Lie-symmetry vector fields of equations from the class have the same properties as those within the class that are presented in Sections 7.3 and 7.4 and additionally satisfy the constraint . In particular, the kernel invariance algebra of the class is , and the linear span of essential Lie-symmetry vector of equations from the class is
The problem of group classification of (1+n)-dimensional linear Schrödinger equations with real-valued potentials reduces to the classification of appropriate subalgebras of the algebra with respect to the equivalence relation generated by the action of .
In view of the relation (7c) between source and target potentials for admissible transformations within the class , an equation from this class can be mapped by a point transformation to an equation from the class if and only if the imaginary part of depends at most on . This is why we do not need to solve the group classification problem for the class with from the very beginning, but its solution can be easily derived from that for the class presented in Theorem 38. It suffices to single out, among the cases listed in Theorem 38, those where the imaginary parts of the corresponding potentials are zero. More specifically, we should set in Case 0, in Cases 1–15, in Cases 5 and 13, in Case 6, in Case 8, in Cases 9–11, in Cases 14 and 15, in Case 16 and in Cases 17 and 18.
10 Conclusion
The initial aim of the paper was to study transformational properties of (1+)-dimensional linear Schrödinger equations with complex-valued potentials, which are of the form (1) and constitute the class denoted by , for arbitrary values of and to carry out the complete group of such equations with . It turned out that this effort requires advancing the theoretical background of symmetry analysis in classes of differential equations, which became the main subject of the paper. We have revisited and extended the entire framework of normalized classes of differential equations, including its basic notions and terminology, and proved the theorems on factoring out point-symmetry groups and maximal Lie invariance algebras within semi-normalized classes and the stronger theorems on splitting such groups and algebras within disjointedly semi-normalized classes. We have also constructed several proof-of-concept examples. This has allowed us to extend the algebraic method of group classification to classes that are semi-normalized in a more general sense, whereas previously only the disjointed semi-normalization with respect to the entire corresponding equivalence groups (in the present terminology) was considered.
For the class with an arbitrary , we first computed its equivalence groupoid by the direct method and thus proved that this class is uniformly semi-normalized with respect to the linear superposition of solutions, which motivates the use of the developed version of the algebraic method of group classification. Knowing , we have easily found the equivalence group of the class and constructed the associated equivalence algebra as the set of infinitesimal generators of one-parameter subgroups of the group .
Working within the framework of the algebraic method, we have reduced the group classification of equations from the class to the classification of specific low-dimensional subalgebras of the algebra or, equivalently, . After analyzing the determining equations for Lie symmetries of equations from the class , we have constructed the linear span of the vector fields from the maximal Lie invariance algebras of these equations. This linear span can be represented as the semidirect sum of the so-called essential subalgebra and an ideal related to the linear superposition of solutions, . We have shown that . For each equation from the class , the representation for induces the similar representation for the maximal Lie invariance algebra of , . The analysis of the determining equations has also resulted in the principal constraints for the appropriate subalgebras of the algebra .
The group classification problem of the class with has been completely solved using the results obtained for the general value of . All possible inequivalent families of potentials possessing essential Lie symmetry extensions have been listed in Theorem 38.
In contrast to the case of single space variable studied in [35], the maximal Lie invariance algebras of equations from the class with can be of essentially greater dimensions, and some of these algebras involve the rotation vector field, which cannot appear for . These two facts make the group classification more complex. To come up with this challenge and make the classification rigorous and efficient, we have split the group classification of the class into different cases depending on three integers, , and . These integers characterize the dimensions of parts of the corresponding essential Lie invariance subalgebra that are related to transformations of time, rotations and generalized time-dependent shifts with respect to space variables, respectively, and are invariant under acting by the equivalence group of the class . Not all values of the tuple are appropriate. The constraint for the appropriate values of is quite obvious and is presented in Lemma 35 in the course of the preliminary analysis of Lie symmetries of equations from the class . This is not the case for the other constraint if , which can be derived only from the entire proof of Theorem 38.
We have also comprehensively studied the subclass of , which is constituted by equations of the form (1) with real-valued potentials. One can easily convert the above results for the class into those for the subclass , merely selecting all the transformational objects related to the class that are relevant for real-valued potentials. In this way, we have straightforwardly constructed the equivalence groupoid and the equivalence group of the subclass , proved its uniform semi-normalization with respect to the linear superposition of solutions, described properties of Lie symmetries of its elements and obtained its complete group classification in the case of two space variables.
The similar study for the subclasses and consisting of equations of the form (1) with time-independent complex- and real-valued potentials, respectively, is not straightforward at all. Even singling out the equivalence groupoids of these subclasses as subgroupoids of the groupoids and given in Theorem 21 and Corollary 42 requires significant efforts. The entire group analysis of the subclasses and is complicated due to the fact that these subclasses are not semi-normalized. A particular consequence of this fact is that one should consider two different group classification problems for each of the subclasses and , up to the equivalence generated by the corresponding equivalence groupoid and up to the equivalence generated by the corresponding equivalence group. The group classification lists for and up to the former equivalence can be formally derived from Theorem 38, but they are not too useful if the equivalence groupoids of and are not known. The study of Lie symmetries of Schrödinger equations from the class was initiated by Boyer in [11]. The classification list of potentials was presented for , including the unfounded claim that it can be converted to the lists for and by “straightforward restriction”. It is of common knowledge that Boyer’s classification list has a number of weaknesses, see, e.g., [54]. There are missed and equivalent cases in this list and essential gaps in its derivation. In particular, it is not clear which kind of equivalence had been used. The above weaknesses have not been comprehensively corrected in the literature. Hence the rigorous group analysis of the subclasses and with explicit indication of the used equivalence is still an open problem.
The results on group classification of the class and its subclasses can be used, in particular, to classify Lie reductions of equations from this class and construct their invariant solutions.
Acknowledgements
The authors are grateful to Peter Basarab-Horwath, Michael Kunzinger and Galyna Popovych for helpful discussions and suggestions. The research of C.K. was supported by International Science Programme (ISP) in collaboration with Eastern Africa Universities Mathematics Programme (EAUMP) and the Abdus Salam International Centre for Theoretical Physics (ICTP). The research of D.R.P. was undertaken thanks to funding from the Canada Research Chairs program, the InnovateNL LeverageR&D program and the NSERC Discovery Grant program. The research was also supported in part by the Ministry of Education, Youth and Sports of the Czech Republic (MŠMT ČR) under RVO funding for IČ47813059. R.O.P. expresses his gratitude for the hospitality shown by the University of Vienna during his long staying at the university. D.R.P. and R.O.P. deeply thank to the Armed Forces of Ukraine and the civil Ukrainian people for their bravery and courage in defense of peace and freedom in Europe and in the entire world from russism.
References
- [1] Akhatov I.S., Gazizov R.K. and Ibragimov N.Kh., Nonlocal symmetries. A heuristic approach, J. Soviet Math. 55 (1991), 1401–1450.
- [2] Ames W.F., Anderson R.L., Dorodnitsyn V.A., Ferapontov E.V., Gazizov R.K., Ibragimov N.H. and Svirshchevskii S.R., CRC handbook of Lie group analysis of differential equations. Vol. 1. Symmetries, exact solutions and conservation laws, CRC Press, Boca Raton, FL, 1994.
- [3] Baran H. and Marvan M., Jets. A software for differential calculus on jet spaces and diffieties. Available at http://jets.math.slu.cz.
- [4] Basarab-Horwath P., Lahno V. and Zhdanov R., The structure of Lie algebras and the classification problem for partial differential equations, Acta Appl. Math. 69 (2001), 43–94, arXiv:math-ph/0005013.
- [5] Belmonte-Beitia J., Pérez-García V.M., Vekslerchik V. and Torres P.J., Lie symmetries, qualitative analysis and exact solutions of nonlinear Schrödinger equations with inhomogeneous nonlinearities, Discrete Contin. Dyn. Syst. Ser. B 9 (2008), 221–233, arXiv:0801.1437.
- [6] Bender C.M., Brody D.C. and Jones H.F., Complex extension of quantum mechanics, Phys. Rev. Lett. 89 (2002), 270401, arXiv:quant-ph/0208076.
- [7] Bihlo A., Dos Santos Cardoso-Bihlo E. and Popovych R.O., Complete group classification of a class of nonlinear wave equations, J. Math. Phys. 53 (2012), 123515, arXiv:1106.4801.
- [8] Bihlo A., Poltavets N. and Popovych R.O., Lie symmetries of two-dimensional shallow water equations with variable bottom topography, Chaos 30 (2020), 073132, arXiv:1911.02097.
- [9] Bihlo A. and Popovych R.O., Group classification of linear evolution equations, J. Math. Anal. Appl. 448 (2017), 982–2015, arXiv:1605.09251.
- [10] Bluman G.W. and Kumei S., Symmetries and differential equations, Springer, New York, 1989.
- [11] Boyer C.P., The maximal ‘kinematical’ invariance group for an arbitrary potential, Helv. Phys. Acta 47 (1974), 589–605.
- [12] Boyer C.P., Sharp R.T. and Winternitz P., Symmetry breaking interactions for the time dependent Schrödinger equation, J. Math. Phys., 17 (1976), 1439–1451.
- [13] Boyko V.M., Lokaziuk O.V. and Popovych R.O., Admissible transformations and Lie symmetries of linear systems of second-order ordinary differential equations, J. Math. Anal. Appl. 539 (2024), 128543, arXiv:2105.05139.
- [14] Carminati J. and Vu K., Symbolic computation and differential equations: Lie symmetries, J. Symbolic Comput. 29 (2000), 95–116.
- [15] Cheviakov A.F., GeM software package for computation of symmetries and conservation laws of differential equations, Comput. Phys. Comm. 176 (2007), 48–61.
- [16] Cheviakov A.F., Symbolic computation of local symmetries of nonlinear and linear partial and ordinary differential equations, Math. Comput. Sci. 4 (2010), 203–222.
- [17] Dos Santos Cardoso-Bihlo E., Bihlo A. and Popovych R.O., Enhanced preliminary group classification of a class of generalized diffusion equations, Commun. Nonlinear Sci. Numer. Simul. 16 (2011), 3622–3638, arXiv:1012.0297.
- [18] Fernández F.M., Guardiola R., Ros J. and Znojil M., A family of complex potentials with real spectrum. J. Phys. A 32 (1999), 3105–3116, arXiv:quant-ph/9812026.
- [19] Fushchich W.I. and Moskaliuk S.S., On some exact solutions of the nonlinear Schrödinger equation in three spatial dimensions, Lett. Nuovo Cimento (2) 31 (1981), 571–576.
- [20] Gagnon L. and Winternitz P., Lie symmetries of a generalised nonlinear Schrödinger equation. I. The symmetry group and its subgroups, J. Phys. A 21 (1988), 1493–1511.
- [21] Gagnon L. and Winternitz P., Lie symmetries of a generalised nonlinear Schrödinger equation. II. Exact solutions, J. Phys. A 22 (1989), 469–497.
- [22] Gagnon L., Grammaticos B., Ramani A. and Winternitz P., Lie symmetries of a generalised nonlinear Schrödinger equation. III. Reductions to third-order ordinary differential equations, J. Phys. A 22 (1989), 499–509.
- [23] Gagnon L. and Winternitz P., Exact solutions of the cubic and quintic nonlinear Schrödinger equation for a cylindrical geometry, Phys. Rev. A (3) 39 (1989), 296–306.
- [24] Gagnon L. and Winternitz P., Symmetry classes of variable coefficient nonlinear Schrödinger equations, J. Phys. A 26 (1993), 7061–7076.
- [25] Gazeau J.P. and Winternitz P., Symmetries of variable coefficient Korteweg–de Vries equations, J. Math. Phys. 33 (1992), 4087–4102.
- [26] Güngör F., Lahno V.I. and Zhdanov R.Z., Symmetry classification of KdV-type nonlinear evolution equations, J. Math. Phys. 45 (2004), 2280–2313, arXiv:nlin/0201063.
- [27] Ibragimov N.H., Aksenov A.V., Baikov V.A., Chugunov V.A., Gazizov R.K. and Meshkov A.G., CRC handbook of Lie group analysis of differential equations. Vol. 2. Applications in engineering and physical sciences, CRC Press, Boca Raton, FL, 1995
- [28] Ibragimov N.H., Torrisi M. and Valenti A., Preliminary group classification of equations , J. Math. Phys. 32 (1991), 2988–2995.
- [29] Ivanova N.M., Popovych R.O. and Sophocleous C., Group analysis of variable coefficient diffusion-convection equations. I. Enhanced group classification, Lobachevskii J. Math. 31 (2010), 100–122, arXiv:0710.2731.
- [30] Kingston J.G. and Sophocleous C., On point transformations of a generalized Burgers equation, Phys. Lett. A 155 (1991), 15–19.
- [31] Kingston J.G. and Sophocleous C., On form-preserving point transformations of partial differential equations, J. Phys. A 31 (1998), 1597–1619.
- [32] Koval S.D., Bihlo A. and Popovych R.O., Extended symmetry analysis of remarkable (1+2)-dimensional Fokker–Planck equation, European J. Appl. Math. 34 (2023), 1067–1098, arXiv:2205.13526.
- [33] Koval S.D. and Popovych R.O., Point and generalized symmetries of the heat equation revisited, J. Math. Anal. Appl. 527 (2023), 127430, arXiv:2208.11073.
- [34] Kurujyibwami C., Equivalence groupoid for (1+2)-dimensional linear Schrödinger equations with complex potentials, J. Phys.: Conf. Ser. 621 (2015), 012008.
- [35] Kurujyibwami C., Basarab-Horwath P. and Popovych R.O., Algebraic method for group classification of (1+1)-dimensional linear Schrödinger equations, Acta Appl. Math. 157 (2018), 171–203, arXiv:1607.04118.
- [36] Kurujyibwami C. and Popovych R.O., Equivalence groupoids and group classification of multidimensional nonlinear Schrödinger equations, J. Math. Anal. Appl. 491 (2020), 124271, arXiv:2003.02781.
- [37] Lahno V. and Zhdanov R., Group classification of nonlinear wave equations, J. Math. Phys. 46 (2005), 053301.
- [38] Lahno V., Zhdanov R. and Magda O., Group classification and exact solutions of nonlinear wave equations, Acta Appl. Math. 91 (2006), 253–313, arXiv:nlin.SI/0405069.
- [39] Lie S., Diskussion der Differentialgleichung , Arch. for Math. 6 (1881), 112–124. (Reprinted in: Lie S., Gesammelte Abhandlungen, Vol. 3, B.G. Teubner, Leipzig and H. Aschehoug & Co, Kristiania, 469–478.)
- [40] Lie S., Über die Integration durch bestimmte Integrale von einer Klasse linear partieller Differentialgleichungen, Arch. for Math. 6 (1881), 328–368; translation by N.H. Ibragimov: Lie S., On integration of a class of linear partial differential equations by means of definite integrals, in CRC Handbook of Lie group analysis of differential equations, vol. 2, CRC Press, Boca Raton, FL, 1995, pp. 473–508.
- [41] Lie S., Vorlesungen über continuierliche Gruppen mit geometrischen und anderen Adwendungen, Teubner, Leipzig, 1893.
- [42] Magadeev B.A., Group classification of nonlinear evolution equations, St. Petersburg Math. J. 5 (1994), 345–359.
- [43] Marvan M., Sufficient set of integrability conditions of an orthonomic system, Found. Comp. Math. 9 (2009), 651–674, arXiv:nlin/0605009.
- [44] Miller W., Symmetry and separation of variables, Addison-Wesley Publishing Co., Reading, MA, 1977.
- [45] Meleshko S.V., Generalization of the equivalence transformations, J. Nonlinear. Math. Phys. 3 (1996), 170–174.
- [46] Muga J.G., Palao J.P., Navarro B. and Egusquiza I.L., Complex absorbing potentials, Phys. Rep. 395 (2004), 357–426.
- [47] Niederer U., The maximal kinematical invariance group of the free Schrödinger equation, Helv. Phys. Acta. 45 (1972), 802–810.
- [48] Niederer U., The maximal kinematical invariance group of the harmonic oscillator, Helv. Phys. Acta. 46 (1973), 191–200.
- [49] Niederer U., The group theoretical equivalence of the free particle, the harmonic oscillator and the free fall, in Proceedings of the 2nd International Colloquium on Group Theoretical Methods in Physics, University of Nijmegen, The Netherlands, 1973.
- [50] Niederer U., The maximal kinematical invariance groups of Schrödinger equations with arbitrary potentials, Helv. Phys. Acta. 47 (1974), 167–172.
- [51] Nikitin A.G., Symmetries of Schrödinger equation with scalar and vector potentials, J. Phys. A 53 (2020) 455202, arXiv:2005.10305.
- [52] Nikitin A.G., Superintegrable and scale-invariant quantum systems with position-dependent mass, Ukrainian Math. J. 74 (2022), 405–419, arXiv:2204.09046.
- [53] Nikitin A.G., Superintegrable quantum mechanical systems with position dependent masses invariant with respect to two parametric Lie groups, J. Phys. A 56 (2023), 395203, arXiv:2304.00638.
- [54] Nikitin A.G. and Zasadko T.M., Group classification of Schrödinger equations with position dependent mass, J. Phys. A 49 (2016), 365204, arXiv:1603.00890.
- [55] Nikitin A.G. and Popovych R.O., Group classification of nonlinear Schrödinger equations, Ukrainian Math. J. 53 (2001), 1255–1265, arXiv:math-ph/0301009.
- [56] Olver P.J., Applications of Lie groups to differential equations, Springer, New York, 1993.
- [57] Opanasenko S., Bihlo A. and Popovych R.O., Group analysis of general Burgers–Korteweg–de Vries equations, J. Math. Phys. 58 (2017), 081511, arXiv:1703.06932.
- [58] Opanasenko S., Boyko V. and Popovych R.O., Enhanced group classification of nonlinear diffusion-reaction equations with gradient-dependent diffusion, J. Math. Anal. Appl. 484 (2020), 123739, arXiv:1804.08776.
- [59] Opanasenko S. and Popovych R.O., Mapping method of group classification, J. Math. Anal. Appl. 513 (2022), 126209, arXiv:2109.11490.
- [60] Ovsiannikov L.V. Group properties of nonlinear heat equation Dokl. AN SSSR 125 (1959), 492–495.
- [61] Ovsiannikov L.V., Group analysis of differential equations, Academic Press, New York, 1982.
- [62] Özemir C. and Güngör F., Symmetry classification of variable coefficient cubic-quintic nonlinear Schrödinger equations, J. Math. Phys. 54 (2013), 023502, arXiv:1201.4033.
- [63] Popovych R.O., Classification of admissible transformations of differential equations, Collection of Works of Institute of Mathematics, Kyiv, 3 (2006), 239–254.
- [64] Popovych R.O. and Eshraghi H., Admissible point transformations of nonlinear Schrödinger equations, in Proceedings of 10th International Conference in MOdern GRoup ANalysis (MOGRAN X) (Larnaca, Cyprus, 2004), University of Cyprus, Nicosia, 2005, pp. 167–174.
- [65] Popovych R.O. and Ivanova N.M., New results on group classification of nonlinear diffusion–convection equations, J. Phys. A 37 (2004), 7547–7565, arXiv:math-ph/0306035.
- [66] Popovych R.O., Ivanova N.M. and Eshragi H., Group classfication of (1+1)-dimensional Schrödinger equations with potentials and power nonlinearities, J. Math. Phys. 45 (2004), 3049–3057, arXiv:math-ph/0311039.
- [67] Popovych R.O., Kunzinger M. and Eshraghi H., Admissible transformations and normalized classes of nonlinear Schrödinger equations, Acta Appl. Math. 109 (2010), 315–359, arXiv:math-ph/0611061.
- [68] Vaneeva O.O., Bihlo A. and Popovych R.O., Generalization of the algebraic method of group classification with application to nonlinear wave and elliptic equations, Commun. Nonlinear Sci. Numer. Simul. 91 (2020), 105419, arXiv:2002.08939.
- [69] Vaneeva O.O., Popovych R.O. and Sophocleous C., Extended group analysis of variable coefficient reaction-diffusion equations with exponential nonlinearities, J. Math. Anal. Appl. 396 (2012), 225–242, arXiv:1111.5198.
- [70] Vaneeva O.O., Sophocleous C. and Leach P., Lie symmetries of generalized Burgers equations: application to boundary-value problems, J. Engrg. Math. 91 (2015), 165–176, arXiv:1303.3548.
- [71] Winternitz P. and Gazeau J.P., Allowed transformations and symmetry classes of variable coefficient Korteweg–de Vries equations, Phys. Lett. A 167 (1992), 246–250.
- [72] Wittkopf A.D., Algorithms and Implementations for Differential Elimination, Ph.D. thesis, Simon Fraser University, 2004.
- [73] Zhdanov R.Z. and Lahno V.I., Group classification of heat conductivity equations with a nonlinear source, J. Phys. A 32 (1999), 7405–7418, arXiv:math-ph/9906003.
- [74] Zhdanov R. and Roman O., On preliminary symmetry classification of nonlinear Schrödinger equations with some applications to Doebner–Goldin models, Rep. Math. Phys. 45 (2000), 273–291.