On the incompleteness of Ibragimov's conservation law theorem and its equivalence to a standard formula using symmetries and adjoint-symmetries

A conservation law theorem stated by N. Ibragimov along with its subsequent extensions are shown to be a special case of a standard formula that uses a pair consisting of a symmetry and an adjoint-symmetry to produce a conservation law through a well-known Fr\'echet derivative identity. Also, the connection of this formula (and of Ibragimov's theorem) to the standard action of symmetries on conservation laws is explained, which accounts for a number of major drawbacks that have appeared in recent work using the formula to generate conservation laws. In particular, the formula can generate trivial conservation laws and does not always yield all non-trivial conservation laws unless the symmetry action on the set of these conservation laws is transitive. It is emphasized that all local conservation laws for any given system of differential equations can be found instead by a general method using adjoint-symmetries. This general method is a kind of adjoint version of the standard Lie method to find all local symmetries and is completely algorithmic. The relationship between this method, Noether's theorem, and the symmetry/adjoint-symmetry formula is discussed.


Introduction
The most well-known method for finding conservation laws of differential equations (DEs) is Noether's theorem [1] which is applicable to any system of one or more DEs admitting a variational formulation in terms of a Lagrangian. Noether's theorem shows that every local symmetry preserving the variational principle of a given variational system yields a non-trivial local conservation law. Moreover, for variational systems that do not possess any differential identities, every non-trivial local conservation law arises from some local symmetry that preserves the variational principle.
However, there are many physically and mathematically interesting DEs that are not variational systems, and this situation has motivated a lot of work in the past few decades to look for some generalization of Noether's theorem which could be applied to non-variational DEs. One direction of work has been to replace the need for a variational principle by introducing some other structure, but still making use of the local symmetries of a given DE system to produce local conservation laws. In fact, a general formula is available that yields local conservation laws from local symmetries combined with solutions of the adjoint of the symmetry determining equations. This formula first appears (to the knowledge of the author) in a 1986 paper by G. Caviglia [2] and was later derived independently in a 1990 Russian paper by F. Lunev [3] as well as in a 1997 paper by the author and G. Bluman [4]. In the latter paper, solutions of the adjoint of the symmetry determining equations were called adjoint-symmetries; these solutions are also known as cosymmetries in the literature on integrable systems [5]. Essentially the same formula appears in a more abstract form in the cohomological framework for finding conservation laws, summarized in Ref. [6,7,8].
In recent years, a similar conservation law formula has been popularized by N. Ibragimov [9,10,11,12], and subsequently extended by others [13,15,16,17], where a "nonlinear selfadjointness" condition is required to hold for the given DE system. However, in several papers [18,19,17], this formula sometimes is seen to produce only trivial conservation laws, and sometimes the formula does not produce all admitted conservation laws. Also, in a number of papers [20,21,13,22,15,16,17], the use of translation symmetries is mysteriously avoided and other more complicated symmetries are used instead, without explanation.
The purpose of the present paper is to make several relevant remarks: (1) Ibragimov's conservation law formula is a simple re-writing of a special case of the earlier formula using symmetries and adjoint-symmetries; (2) Ibragimov's "nonlinear self-adjointness" condition in its most general form is equivalent to the existence of an adjoint-symmetry for a general DE system, and reduces to the existence of a symmetry in the case of a variational DE system; (3) this formula does not always yield all admitted local conservation laws and it produces trivial conservation laws whenever the symmetry is a translation and the adjoint-symmetry is translation-invariant; (4) the computation to find adjoint-symmetries (and hence to apply the formula) is just as algorithmic as the computation of local symmetries; (5) most importantly, if all adjoint-symmetries are known for a given DE system (whether or not it has a variational formulation) then they can be used directly to obtain all local conservation laws, providing a kind of generalization of Noether's theorem to general DE systems.
All of these remarks have been pointed out briefly in Ref. [23], and remark (2) has been discussed in Ref. [16,17], but it seems worthwhile to give a comprehensive discussion for all of the remarks (1)- (5), with examples, as the formula continues to be used in recent papers when a complete, general method for finding all local conservation laws is available instead. In particular, for any given DE system, a full generalization of the content of Noether's theorem is provided by a direct method using adjoint-symmetries, based on the framework shown in Ref. [24,25] and presented in an algorithmic fashion in Ref. [4,26,27,28]. In the case when a DE system is variational, adjoint-symmetries reduce to symmetries, and the direct method reproduces the relationship between symmetries and conservation laws in Noether's theorem but without the need for a Lagrangian. A detailed review and further development of this general method appears in Ref. [29]. Consequently, there is no need for any kind of special methods to find local conservation laws, just as there is no need to use special methods to find local symmetries, because the relevant determining equations can be solved in a direct algorithmic manner.
The remainder of the present paper is organized as follows. Remarks (1) and (2) will be demonstrated in section 2. Remark (3), along with some further consequences and properties related to the action of symmetries, will be explained in section 3. Remarks (4) and (5) will be briefly discussed in section 4. Throughout, the class of DEs u tt − u xx + a(u)(u 2 t − u 2 x ) + b(u)u t + c(u)u x + m(u) = 0 with be used as a running example to illustrate the main points, and the notation in Ibragimov's work will be used to allow the simplest possible comparison of the results. Some concluding remarks are made in section 5. Many examples of conservation laws of wave equations and other evolution equations can be found in Ref. [30,24,31] and references therein.
2. Symmetries, adjoint-symmetries, and "nonlinear self-adjointness" As preliminaries, a few basic tools from variational calculus will be reviewed. Let x = (x 1 , . . . , x n ) be n ≥ 1 independent variables, u = (u 1 , . . . , u m ) be m ≥ 1 dependent variables, and let ∂ k u denote all kth order partial derivatives of u with respect to x. Introduce an index notation for the components of x and u: x i , i = 1, . . . , n; and u α , α = 1, . . . , m. In this notation, the components of ∂ k u are given by u α i 1 ···i k , α = 1, . . . , m, i q = 1, . . . , n, with q = 1, . . . , k. Summation is assumed over each pair of repeated indices in any expression. The coordinate space J = (x, u, ∂u, ∂ 2 u, . . .) is called the jet space associated with the variables x, u. A differential function is a locally smooth function of finitely many variables in J. Total derivatives with respect to x applied to differential functions are denoted The necessary tools that will now be introduced are the Fréchet derivative and its adjoint derivative, the Euler operator and its product rule, and the Helmholtz conditions. Given a set of M ≥ 1 differential functions f a (x, u, ∂u, . . . , ∂ N u), a = 1, . . . , M, with differential order N ≥ 1, the Fréchet derivative is the linearization of the functions as defined by This linearization can be viewed as a local directional derivative in jet space, corresponding to the action of a generatorX = w α ∂ u α in characteristic form,X(f ) = δ w f , where w = (w 1 (x, u, ∂u, . . . , ∂ k u), . . . , w m (x, u, ∂u, . . . , ∂ k u)) is a set of m arbitrary differential functions.
It is useful also to view the Fréchet derivative as a linear differential operator acting on w. Then integration by parts defines the Fréchet adjoint derivative which is a linear differential operator acting on a set of M ≥ 1 arbitrary differential functions v = (v 1 (x, u, ∂u, . . . , ∂ k u), . . . , v M (x, u, ∂u, . . . , ∂ k u)). These two derivatives (2.1) and (2.2) where the associated vector Ψ i (v, w; f ) is given by the explicit formula The Euler operator E u α , or variational derivative, is defined in terms of the Fréchet derivative through the variational relation which is obtained from integration by parts, yielding Here, for simplicity, f (x, u, ∂u, . . . , ∂ N u) is a single differential function. In particular, an explicit formula for Φ i (w; f ) is given by The Euler operator (2.6) has the following three important properties: First, it obeys the product rule Second, its kernel is given by total divergences holding for some differential vector function F i . Third, its image consists of differential functions characterized by the Helmholtz conditions where w α is a set of arbitrary differential functions. There are several common alternative notations for the Fréchet derivative and its adjoint: (v) appears in the literature on integrable systems and in Ref. [29]; [24]; δ w f = L[u]w is used in the early work of Anco and Bluman [4,26,27] and in the book [31]. In contrast, Ibragimov [9,12] uses

Conservation laws and symmetries.
Consider an Nth-order system of M ≥ 1 DEs F = (F 1 (x, u, ∂u, . . . , ∂ N u), . . . , F M (x, u, ∂u, . . . , ∂ N u)) = 0. (2.14) The space of solutions u(x) of the system will be denoted E. When the number of independent variables x is n = 1, each DE is an ordinary differential equation (ODE), whereas when the number of independent variables x is n ≥ 2, each DE is a partial differential equation (PDE). The number, m, of dependent variables u need not be the same as the number, M, of DEs in the system. A local infinitesimal symmetry [24,31,32] of a given DE system (2.14) is a generator whose prolongation leaves invariant the DE system, which holds on the whole solution space E of the system. (In this determining equation, the notation E means that the given DE system as well as its differential consequences are to be used.) The differential functions ξ i and η α in the symmetry generator are called the symmetry characteristic functions. When acting on the solution space E, an infinitesimal symmetry generator can be formally exponentiated to produce a one-parameter group of transformations exp(ǫprX), with parameter ǫ, where the infinitesimal transformation is given by Two infinitesimal symmetries are equivalent if they have the same action (2.17) on the solution space E of a given DE system. An infinitesimal symmetry is thereby called trivial if it leaves all solutions u(x) unchanged. This occurs iff its characteristic functions satisfy the relation 18) The corresponding generator (2.15) of a trivial symmetry is thus given by which has the prolongation prX triv. | E = ξ i D i . Conversely, any generator of this form on the solution space E represents a trivial symmetry. Thus, any two generators that differ by a trivial symmetry are equivalent. The differential order of an infinitesimal symmetry is defined to be the smallest differential order among all equivalent generators. Any symmetry generator is equivalent to a generator given bŷ holding on the whole solution space E of the given system. Note that the action of prX is the same as a Fréchet derivative (2.1), and hence an equivalent, modern formulation [24,31,29] of this invariance (2.22) is given by the symmetry determining equation (2.23) (Recall, the notation E means that the given DE system as well as its differential consequences are to be used in these determining equations.) In jet space J, a group of transformations exp(ǫprX) with a non-trivial generator X in general will not act in a closed form on x, u, and derivatives ∂ k u up to a finite order, except [24,31] for point transformations acting on (x, u), and contact transformations acting on (x, u, ∂u). Moreover, a contact transformation is a prolonged point transformation when the number of dependent variables is m > 1 [24,31]. A point symmetry is defined as a symmetry transformation group on (x, u), whose generator is given by characteristic functions of the form corresponding to the infinitesimal point transformation Likewise, a contact symmetry is defined as a symmetry transformation group on (x, u, ∂u) whose generator corresponds to an infinitesimal transformation that preserves the contact relations u α i = D i u α . The set of all admitted point symmetries and contact symmetries for a given DE system comprises its group of Lie symmetries. The corresponding generators of this group comprise a Lie algebra [24,32,31].
A local conservation law of a given DE system (2.14) is a divergence equation which holds on the whole solution space E of the system, where C = (C 1 (x, u, ∂u, . . . , ∂ r u), . . . , C n (x, u, ∂u, . . . , ∂ r u)) (2.27) is the conserved current vector. In the case when one of the independent variables represents a time coordinate, and the remaining n−1 independent variables represent space coordinates, namely x = (t, x 1 , . . . , x n−1 ), then C 1 = T is a conserved density and (C 2 , . . . , C n ) = X is a spatial flux vector, while the conservation law has the form of a local continuity equation (D t T + Div X)| E = 0. (Similarly to the symmetry determining equation, the notation E here means that the given DE system as well as its differential consequences are to be used.) A conservation law (2.26) is locally trivial if holds for some differential antisymmetric tensor function Θ ij (x, u, ∂u, . . . , ∂ r−1 u) on E, since any total curl is identically divergence free, D i (D j Θ ij ) = D i D j Θ ij = 0 due to commutativity of total derivatives. Two conservation laws are said to be locally equivalent if, on the solution space E, their conserved currents differ by a locally trivial current (2.28). The differential order of a conservation law is defined to be the smallest differential order among all locally equivalent conserved currents. (Sometimes a local conservation law is itself defined as the equivalence class of locally equivalent conserved currents.) For a given DE system (2.14), the set of all non-trivial local conservation laws (up to local equivalence) forms a vector space on which the local symmetries of the system have a natural action [33,24,31]. In particular, the infinitesimal action of a symmetry (2.15) on a conserved current (2.27) is given by [24] When the symmetry is expressed in characteristic form (2.20), its action has the simple form The conserved currents C i X and C iX are locally equivalent, 32) which follows from the relation (2.21).
A DE system is variational if it arises as the Euler-Lagrange equations of a local Lagrangian. This requires that the number of equations in the system is the same as the number of dependent variables, M = m, and that the differential order N of the system is even, in which case the system is given by where the Lagrangian is a differential function L(x, u, ∂u, . . . , ∂ N/2 u). (2.34) The necessary and sufficient conditions [24,32,31,29] for a given DE system (2.14) to be variational consist of the Helmholtz conditions (2.13), which are given by where w α is a set of arbitrary differential functions. Note these conditions (2.35) are required to hold identically in jet space J (and not just on the solution space E of the DE system).

2.2.
Ibragimov's conservation law formula. The starting point is the well-known observation [24] that any Nth-order system of M ≥ 1 DEs (2.14) can be embedded into a larger system by appending an "adjoint variable" for each DE in the system, where this set of M ≥ 1 variables v = (v 1 , . . . , v M ) is taken to satisfy the adjoint of the linearization of the original DE system. Specifically, the enlarged DE system is given by for u α (x) and v a (x), in Ibragimov's notation. This system (2.36)-(2.37) comprises the Euler-Lagrange equations of the Lagrangian function This embedding relationship can used to show that every symmetry of the original DE system (2.36) can be extended to a variational symmetry of the Euler-Lagrange system (2.39). The proof is simplest when the symmetries are formulated in characteristic form (2.20). LetX = P α (x, u, ∂u, . . . , ∂ r u)∂/∂u α (2.40) be any local symmetry generator (in characteristic form) admitted by the DE system (2.36). Under some mild regularity conditions [29] on the form of these DEs, the symmetry determining equation (2.23) implies that the characteristic functions P α satisfy 42) is some linear differential operator whose coefficients R P b a , R P bi a , . . . , R P bi 1 ···ir a are differential functions that are non-singular on solution space E of the DE system (2.14). Now consider the action of this symmetry generator (2.40) on the Lagrangian (2.38). From the operator relation (2.41) followed by integration by parts, the symmetry action is given by ) is the adjoint of the operator (2.42), with the non-singular coefficients Although the Lagrangian is not preserved, the expression (2.43) for the symmetry action shows that if the symmetry is extended to act on v viâ then under this extended symmetry the Lagrangian will be invariant up to a total divergence, This completes the proof. A useful remark is that the vectorΘ i in the total divergence (2.47) is a linear expression in terms of F a (and total derivatives of F a ), and hence this vector vanishes whenever u(x) is a solution of the DE system (2.36). Consequently,Θ i is a trivial current for the Euler-Lagrange system (2.39).
Some minor remarks are that the proof given by Ibragimov [9] does not take advantage of the simplicity of working with symmetries in characteristic form and also glosses over the need for some regularity conditions on the DE system so that the symmetry operator relation (2.41) will hold. Moreover, that proof is stated only for DE systems in which the number of equations is the same as the number of dependent variables, M = m. Now, since the extended symmetry (2.46) is variational, Noether's theorem can be applied to obtain a corresponding conservation law for the Euler-Lagrange system (2.39), without the need for any additional conditions. The formula in Noether's theorem comes from applying the variational identity (2.5) to the Lagrangian (2.38), which yields for any generatorX The total divergence term D i Φ i (P ; L) is given by the formula (2.8) derived using the Euler operator (2.6). This yields When this variational identity (2.48) is combined with the action (2.47) of the variational symmetry (2.46) on the Lagrangian, the following Noether relation is obtained: ) is any solution of the Euler-Lagrange system (2.39), the Noether relation (2.51) yields a local conservation law where E(u, v) denotes the solution space of the system (2.39) (including its differential consequences). This conservation law is locally equivalent to the conservation law formula underlying Ibragimov's work [9,12], which is given by A side remark is that the locally equivalent conservation law (2.53) also can be derived from Noether's theorem if the extended symmetry (2.46) is expressed in canonical form as obtained from relations (2.20)-(2.21). In particular, the corresponding form of the variational identity (2.48) becomes while the action of the symmetry (2.54) on the Lagrangian is given by is a solution of the DE system (2.36). Hence, the Noether relation obtained from combining equations (2.55)-(2.56) yields the conserved current C i = Φ i (P, v; F ) − ξ i L modulo the locally trivial current Θ i . If the original symmetry (2.40) being used is a point symmetry, then this trivial current Θ i can be shown to vanish identically, which is the situation considered in Ibragimov's papers [9,12] and in nearly all subsequent applications in the literature.
2.3. "Nonlinear self-adjointness". The conservation law (2.52) holds for all solutions (u(x), v(x)) of the Euler-Lagrange system (2.39). It seems natural to restrict this to solutions of the original DE system (2.36) for u(x) by putting v = 0. However, the resulting conserved current is trivial, Φ i (P ; L)| v=0 = Φ i (P ; 0) = 0, because L is a linear expression in terms of v. Consequently, some other way must be sought to project the solution space E(u, v) of the Euler-Lagrange system onto the solution space E of the original DE system (2.36).
Ibragimov's first paper [9] proposes to put v = u, which is clearly a significant restriction on the form of the original DE system (2.36). In particular, this requires that F * α | v=u = F a hold identically, where the DE system is assumed to have the same number of equations as the number of dependent variables, M = m, which allows the indices a = α to be identified. He calls such a DE system F a = 0 "strictly self-adjoint". This definition is motivated by the case of a linear DE system, since linearity implies (δ u F ) a = F a and (δ * u F ) α = F * α | v=u are identities, whereby a linear DE system with M = m is "strictly self-adjoint" iff it satisfies (δF ) α = (δ * F ) α which is the condition for self-adjointness of a linear system. However, for nonlinear DE systems, the definition of "strictly self-adjoint" conflicts with the standard of definition [24,6] in variational calculus that a general DE system F a = 0 is self-adjoint iff its associated Fréchet derivative operator is self-adjoint, (δF ) a = (δ * F ) α , which requires M = m.
Ibragimov subsequently [10] proposed to have v = φ(u), which he called "quasi-selfadjointness". A more general proposal v = φ(x, u) was then introduced first in Ref. [13] and shortly later appears in Ibragimov's next paper [12], with the condition that F * α | v=φ(x,u) = λ α β F β must hold for some coefficients λ α β , again with M = m. This condition is called "weak self-adjointness" in Ref. [13] and "nonlinear self-adjointness" in Ref. [12]. Ibragimov also mentions an extension of this definition to v = φ(x, u, ∂u, . . . , ∂ s u), but does not pursue it. Later he applies this definition in Ref. [14] to a specific PDE, where λ α β is extended to be a linear differential operator. However, unlike in the previous papers, no conservation laws are found from using this extension. A subsequent paper [15] then uses this extension, which is called "nonlinear self-adjointness through a differential substitution", to obtain conservation laws for several similar PDEs. Finally, the same definition is stated more generally in Ref. [17] for DE systems with M = m: where the coefficients λ α β , λ α βi , . . . , λ α βi 1 ···ir are differential functions. These developments lead to the following conservation law theorem, which is a generalization of Ibragimov's main theorem [9,12] to arbitrary DE systems (not restricted by M = m), combined with the use of a differential substitution [12,15,17].
Theorem 1. Suppose a system of DEs (2.14) satisfies . . , ∂ s u) that are nonsingular on the solution space E of the DE system, where F * α is the adjoint linearization (2.2) of the system. Then any local symmetry admitted by the DE system yields a local conservation law (2.26) given in an explicit form by the conserved current (2.50) with v a = φ a and P α = η α − ξ i u α i . An important remark is that all of the functions φ a , λ α a , λ α ai , . . . , λ α ai 1 ···ip must be nonsingular on E, as otherwise the condition (2.58) can be satisfied in a trivial way. This point is not mentioned in any of the previous work [9,12,13,15,17].
The "nonlinear self-adjointness" condition (2.58) turns out to have a simple connection to the determining equations for symmetries. This connection is somewhat obscured by the unfortunate use of non-standard definitions and non-standard notation in Ref. [9,12]. Nevertheless, it is straightforward to show that equation (2.58) is precisely the adjoint of the determining equation (2.23) for symmetries formulated as an operator equation (2.41).
2.4. Adjoint-symmetries and a formula for generating conservation laws. For any given DE system (2.14), the adjoint of the symmetry determining equation (2.23) is given by for a set of differential functions Q a (x, u, ∂u, . . . , ∂ r u). (Similarly to the symmetry determining equation, the notation E here means that the given DE system as well as its differential consequences are to be used.) These differential functions are called an adjoint-symmetry [4], in analogy to the characteristic functions of a symmetry (2.40), and so the equation (2.60) is called the adjoint-symmetry determining equation. As shown in Ref. [29], this analogy has a concrete geometrical meaning in the case when a DE system is an evolutionary system , where t is a time coordinate and x i , i = 1, . . . , n−1, are space coordinates. In this case, Q α can be viewed as the coefficients of a 1-form or a covector Q α du α , in analogy to P α being the coefficients of a vector P α ∂/∂u α . The condition for P α ∂/∂u α to be a symmetry can be formulated as L f P α ∂/∂u α | E = 0 where L f denotes the Lie derivative [24,29] with respect to the time evolution vectorX = f α ∂/∂u α . Then the condition for Q α du α to be an adjoint-symmetry is equivalent to L f Q α du α | E = 0. (Note the awkwardness in the index positions here comes from Ibragimov's choice of index placement F α for a DE system with M = m. A better notation would be F α , and F a when M = m, which is used in Ref. [4,26,27,31].) In the case when a DE system is variational (2.33), the symmetry determining equation Then the adjoint-symmetry determining equation (2.60) reduces to the symmetry determining equation (2.23), with Q a (x, u, ∂u, . . . , ∂ r u) = P α (x, u, ∂u, . . . , ∂ r u), where the indices a = α can be identified, due to M = m. Consequently, adjoint-symmetries of any variational DE system are the same as symmetries.
Other aspects of adjoint-symmetries and their connection to symmetries are discussed in Ref. [34].
Now, under some mild regularity conditions [29] on the form of a general DE system (2.14), the adjoint-symmetry determining equation (2.60) implies that the functions Q a satisfy is some linear differential operator whose coefficients R Q b α , R Q bi α , . . . , R Q bi 1 ···ir α are differential functions that are non-singular on solution space E of the DE system (2.14). In Ibragimov Therefore, the following equivalence has been established.
Proposition 2. For a general DE system (2.14), the condition (2.58) of "nonlinear selfadjointness" coincides with the condition of existence of an adjoint-symmetry (2.60). When a DE system is variational (2.33), these conditions reduce to the condition of existence of a symmetry.
One remark is that the formulation of "nonlinear self-adjointness" given here is more general than what appears in Ref. [12,15,17] since those formulations assume that the DE system has the same number of equations as the number of dependent variables, M = m. Another remark is that the meaning of "nonlinear self-adjointness" shown here in the case of variational DE systems has not previously appeared in the literature.
Example: Consider the class of semilinear wave equations , and a mass-type coefficient m(u). In Ref. [20], conditions under which a slightly more general family of wave equations is "nonlinearly self-adjoint" (2.58) are stated for v = φ(u). These results will be generalized here by considering v = φ(t, x, u). A first observation is that this class of wave equations admits an equivalence transformation u →ũ = f (u), with f ′ = 0, which can be used to put a = 0 by f (u) = exp(A(u))du where A ′ = a. (Equivalence transformations were not considered in Ref. [20], and so their results are considerably more complicated than is necessary.) This transformation gives (2.63) In Ibragimov's notation, the condition of "nonlinear self-adjointness" with v = φ(t, x, u) is For comparison, the determining equation (2.23) for local symmetriesX = P (t, x, u, u t , u x , . . .)∂/∂u (in characteristic form) is given by Its adjoint is obtained by multiplying by Q(t, x, u, u t , u x , . . .) and integrating by parts, which yields After the D x terms are expanded out, this gives the determining equation (2.60) for local adjoint-symmetries which coincides with the "nonlinear self-adjointness" condition (2.65) extended to differential substitutions [13,12,16] given by v = Q(t, x, u, u t , u x , . . .). All adjoint-symmetries of lowestorder form Q(t, x, u) can be found in a straightforward way. After Q(t, x, u) is substituted into the determining equation (2.67), and u tt is eliminated through the wave equation (2.63), the determining equation splits with respect to the variables u t and u x , yielding a linear overdetermined system of four equations (after some simplifications): It is straightforward to derive and solve this determining system by Maple. Hereafter, the conditions b ′ = 0, c ′ = 0, m ′′ = 0, m(0) = 0 (2.69) will be imposed, which corresponds to studying wave equations (2.63) whose lower-order terms are nonlinear and homogeneous. The general solution of the determining system (2.68) then comprises three distinct cases (as obtained using the Maple package 'rifsimp'), after merging. This leads to the following complete classification of solution cases shown in table 1. The table is organized by listing each solution Q and the conditions on b, c, m for which it exists. (From these conditions, a classification of maximal linear spaces of multipliers can be easily derived.) Note that if the transformation u →ũ = exp(A(u))du is inverted, then Q transforms toQ = exp(A(u))Q. (Also note that, under the restriction Q = φ(u) considered in Ref. [20], the classification reduces to just the first case with m = const. and Q = 1.) The Fréchet derivative operator in the symmetry determining equation (2.23) and the adjoint of this operator in the adjoint-symmetry determining equation (2.60) are related by the integration-by-parts formula (2.3). For a general DE system (2.14), this formula is given by where the vector Ψ i (P, Q; F ) is given by the explicit expression (2.4) with v = Q, w = P , and f = F . As shown in Ref. [2,3,4], this vector Ψ i (P, Q; F ) will be a conserved current whenever the differential functions P α and Q a respectively satisfy the symmetry and adjointsymmetry determining equations. Moreover, it is straightforward to see which implies Ψ i (P, v; F ) = Φ i (P ; L) holds (up to the possible addition of a total curl). When the relation (2.72) is combined with Propositions 1 and 2, the following main result is obtained.
Theorem 2. For any DE system (2.14) admitting an adjoint-symmetry (2.60) (namely, a "nonlinearly self-adjoint system" in the general sense), the conserved current (2.50) derived from applying Noether's theorem to the extended Euler-Lagrange system (2.39) using any given symmetry (2.46) is equivalent to the conserved current obtained using the adjointsymmetry/symmetry formula (2.70).
This theorem shows that the "nonlinear self-adjointness" method based on Ibragimov's theorem as developed in papers [9,12,13,15,17] for DE systems with M = m is just a special case of the adjoint-symmetry/symmetry formula (2.70) introduced for general DE systems in prior papers Ref. [2,3,4] which were never cited. Moreover, the adjoint-symmetry/symmetry formula (2.70) has the advantage that there is no need to extend the given DE system by artificially adjoining variables to get an Euler-Lagrange system.
Another major advantage of the adjoint-symmetry/symmetry formula is that it can be used to show how the resulting local conservation laws are, in general, not necessarily nontrivial and comprise only a subset of all of the non-trivial local conservation laws admitted by a given DE system. In particular, in many applications of Theorem 1, it is found that some non-trivial symmetries, particularly translation symmetries, only yield trivial conservation laws [18,19,17], and that some local conservation laws are not produced even when all admitted symmetries are used. These observations turn out to have a simple explanation through the equivalence of Theorem 1 and the adjoint-symmetry/symmetry formula (2.70), as explained in the next section.
Example: For the semilinear wave equation (2.63), the extended Euler-Lagrange system in Ibragimov's notation consists of where F * is defined by the adjoint-symmetry equation ( where R * P is the adjoint of the operator R P defined by relation (2.41) for the point symmetry holding off of the solution space of the wave equation (2.74a). In particular, R P can be obtained by a straightforward computation of δ P F = R P (F ), where the terms in δ P F are simplified by using the equations τ u = ξ u = 0, τ t = ξ x , and τ x = ξ t that arise from splitting the determining equation (2.66). This yields and thus R * P = τ D t + ξD x + η u . (2.77) Hence, the variational symmetry is simply X ext. = τ ∂/∂t + ξ∂/∂x + η∂/∂u − η u v∂/∂v (2.78) which is a point symmetry.
The action of this variational symmetry on the Lagrangian L = vF is given by F . This symmetry action then can be combined with the variational identity (2.55) to get the Noether relation are obtained from formula (2.50). This yields a conservation law on the solution space E(u, v) of the Euler-Lagrange system F = 0, F * = 0. Since L| E(u,v) = 0, this conservation law is locally equivalent to the conservation law (2.52) which is given by Moreover, from the identity (2.70) relating the symmetry equation (2.66) and the adjointsymmetry equation (2.67), the conserved current (Ĉ t ,Ĉ x ) in the conservation law (2.83) is the same as the conserved current (Ψ t , Ψ x ) in the adjoint-symmetry/symmetry formula In Ref. [20], the conservation law formula (2.82) is used to obtain a single local conservation law for a special case of the wave equation (2.63) given by b = −c = − ln(u) and d = 0, corresponding toũ tt −ũ xx − (ũ 2 t −ũ 2 x ) +ũ(ũ t −ũ x ) = 0 after an equivalence transformation u →ũ = e −u is made. The formula is applied to the adjoint-symmetryQ = e −ũ and the point symmetry X = e (t+x)/2 ∂/∂ũ with characteristicP = e (t+x)/2 , which respectively correspond to Q = 1 and X = e (t+x)/2 u∂/∂u with P = e (t+x)/2 u. The likely reason why the obvious translation symmetries X = ∂/∂t and X = ∂/∂x were not considered in Ref. [20] is that these symmetries lead to locally trivial conservation laws whenQ = e −ũ is used.
To illustrate the situation, consider the translation symmetries Notice that for Q = const. the conserved currents (Ψ t , Ψ x ) obtained from the two translation symmetries vanish. This implies that Ibragimov's theorem (2.82) yields just trivial conserved currents (Ψ t , Ψ x ) for some cases of the wave equation (2.63) when a non-trivial conserved current exists. A full explanation of why this occurs will be given in the next section.

Properties of conservation laws generated by the adjoint-symmetry/symmetry formula and Ibragimov's theorem
To determine when a conserved current is locally trivial, or when two conserved currents are locally equivalent, it is useful to have a characteristic (canonical) form for local conservation laws, in analogy to the characteristic form for local symmetries.
Any local conservation law (2.26) can be expressed as a divergence identity [24] D i C i = R C a F a + R C ai D i F a + · · · + R C ai 1 ···ir D i 1 · · · D ir F a (3.1) by moving off of the solution space E of the system, where R C a , R C ai , . . . , R C ai 1 ···ir a are some differential functions that are non-singular on E, under some mild regularity conditions [29] on the form of the DEs (2.14). Integration by parts on the terms on the right-hand side in this identity (3.1) then yields ) reduces to the conserved vector in the given conservation law (2.26). Hence, is a locally equivalent conservation law. The identity (3.2) is called the characteristic equation [24] for the conservation law (2.26), and the set of differential functions (3.3) is called the conservation law multiplier [24]. In general a set of functions f a (t, x, u, ∂u, ∂ 2 u, . . . ∂ s u) will be a multiplier iff it is non-singular on E and its summed product with the DEs F a in the system has the form of a total divergence. For a given local conservation law, the multiplier arising from the integration by parts formula (3.3) will be unique iff the coefficient functions in the characteristic equation (3.1) are uniquely determined by the conserved vector C i . This uniqueness holds straightforwardly for any DE system consisting of a single equation that can be expressed in a solved form for a leading derivative [29]. For DE systems containing more than one equation, some additional technical requirements are necessary [24]. In particular, it is necessary that a DE system have no differential identities [24], and it is sufficient that a DE system have a generalized Cauchy-Kovalevskaya form [24,25,26,27]. A concrete necessary and sufficient condition, which leads to the following uniqueness result, is stated in Ref. [29].
Proposition 3. For any closed DE system (2.14) having a solved form in terms of leading derivatives and having no differential identities, a conserved current is locally trivial (2.28) iff its corresponding multiplier (3.3) vanishes when evaluated on the solution space of the system. This class of DE systems includes nearly all systems of physical interest, apart from systems such as the Maxwell equations and the incompressible fluid equations, which possess differential identities. Often the distinction between systems with and without differential identities is overlooked in the literature on conservation law multipliers.
The importance of Proposition 3 is that in a wide class of DE systems it establishes that a unique characteristic form for locally equivalent conservation laws is provided by multipliers. From this result, it is now straightforward to derive a simple condition to detect when a local conservation law given by the adjoint-symmetry/symmetry formula (2.70) is locally trivial (2.28).
Let P α (x, u, ∂u, . . . , ∂ r u) be the characteristic functions defining a symmetry (2.23), and let Q a (x, u, ∂u, . . . , ∂ s u) be a set of differential functions defining an adjoint-symmetry (2.60). Then the adjoint-symmetry/symmetry formula (2.70) yields a local conservation law (2.71). The characteristic equation of this conservation law is given by substituting the symmetry identity (2.41) and the adjoint-symmetry identity (2.61) into the formula (2.70) to get Integration by parts gives (3.8) Hence, the conservation law multiplier is given by [36] This yields the following result.
Proposition 4. The adjoint-symmetry/symmetry formula (2.70) for a given DE system (2.14) produces a locally trivial conservation law if the condition holds for the given symmetry and adjoint-symmetry pair, where P α (x, u, ∂u, . . . , ∂ r u) is the set of characteristic functions of the symmetry (2.23) and Q a (x, u, ∂u, . . . , ∂ s u) is the set of functions defining the adjoint-symmetry (2.60). This condition (3.10) is also sufficient whenever the DE system (2.14) belongs to the class stated in Proposition 3.
Through the equivalence stated in Theorem 2, which relates the adjointsymmetry/symmetry formula (2.70) and the generalized version of Ibragimov's conservation law formula in Theorem 1, it follows that "nonlinear self-adjointness" through a differential substitution with v = Q produces a conservation law (2.50) that is locally trivial when Q a (x, u, ∂u, . . . , ∂ s u) and P α (x, u, ∂u, . . . , ∂ r u) satisfy condition (3.10).
A useful remark is that the triviality condition (3.10) can be checked directly, without the need to derive the local conservation law itself.
Example: For the semilinear wave equation (2.63), consider the conserved currents obtained in table 2, which are generated from the three adjoint-symmetries Q 1 = e m 2 t+m 3 x , Q 2 = e αx+βt , Q 3 = e γx q(x ∓ t), (3.11) and the two translation symmetries (2.86). The operators δ P F = R P (F ) associated to the characteristics P 1 = −u t , P 2 = −u x (3.12) of these two symmetries (2.86) are given by the formula (2.76), which yields For adjoint-symmetries of the form Q(t, x, u), the operator δ * Q F = R Q (F ) is easily found to be R Q = Q u . (3.14) Hence, the operators associated to the three adjoint-symmetries (3.11) are simply The triviality condition (3.10) is then given by 16) This shows that the two conserved currents obtained from Q 1 will be trivial when m 2 = 0 and m 3 = 0 hold, respectively, and that likewise the two conserved currents obtained from Q 2 will be trivial when β = 0 and α = 0 hold, respectively. Similarly, for Q 3 , the first conserved current will be trivial when q ′ = 0 holds, while the second conserved current will be trivial when q ′ + γq = 0 holds, corresponding to q = e −γ(x±t) . These trivial cases can be seen to occur directly from the explicit expressions for the conserved currents (Ψ t , Ψ x ) in table 2.
In general, while the adjoint-symmetry/symmetry formula (2.70) (and hence Ibragimov's theorem) looks very appealing, it has major drawbacks that in many examples [20,21,13,22,15,16,17] the selection of a symmetry must be fitted to the form of the adjoint-symmetry to produce a non-trivial conservation law, and that in other examples [18,19,17] no nontrivial conservation laws are produced when only translation symmetries are available. More importantly, it is not (as is sometimes claimed) a generalization of Noether's theorem to non-variational DE systems.
As a reinforcement of these statements, consider the situation of variational DE systems, where adjoint-symmetries coincide with symmetries. Then the adjoint-symmetry/symmetry formula (2.70) produces a conserved current directly from any pair of symmetries admitted by a given variational DE system. But, from Noether's theorem, this conserved current must also arise directly from some variational symmetry of the system. Moreover, if the pair of symmetries being used are variational symmetries that happen to commute with each other, then the resulting conserved current turns out to be trivial as shown in Ref. [35].
To understand these aspects and other properties of the formula, the determining equations for multipliers are needed.
3.1. Multiplier determining equations. All conservation law multipliers for any given DE system can be determined from the property (2.10)-(2.11) that a differential function is a total divergence iff it is annihilated by the Euler operator (2.6). Specifically, when this property is applied directly to the characteristic equation (3.2) for local conservation laws, it yields the determining equations E u α (Q a C F a ) = 0, α = 1, . . . , m (3.17) which are necessary and sufficient [24] for a set of differential functions Q a C (x, u, ∂u, . . . , ∂ s u) to be a multiplier for a local conservation law (2.26). Note these equations (3.17) must hold identically in jet space (and not just on the solution space E of the DE system).
The multiplier determining equations (3.17) have a close connection to the determining equation (2.60) for adjoint-symmetries. This can be immediately seen from the product rule (2.9) obeyed by the Euler operator, which gives holding identically in jet space J(x, u, ∂u, ∂ 2 u, . . .). Notice if this equation (3.18) is restricted to the solution space E ⊂ J of the given DE system (2.14), then it coincides with the adjoint-symmetry determining equation (2.60). Hence, every conservation law multiplier is an adjoint-symmetry. This is a well-known result [24,4,26,27,31,29]. What is not so well-known are the other conditions [4,29] that an adjoint-symmetry must satisfy to be a conservation law multiplier. These conditions arise from splitting the determining equation (3.18) with respect to F a and its total derivatives. As shown in Ref. [29], the splitting can be derived by using the adjoint-symmetry identity (2.61) combined with the expression Then, in the determining equation (3.18), the coefficients of F a , D i F a , and so on yield the system of equations [29] (δ * and where R Q a α and R Q ai 1 ···iq α are the coefficient functions of the linear differential operator (2.62) determined by equation (3.20), and where E u α is the Euler operator (2.6) and E (i 1 ···iq) u α is a higher-order Euler operator defined by [24,29]  for an arbitrary differential function f (x, u, ∂u, . . . , ∂ s u). This system (3.20)-(3.21) constitutes a determining system for conservation law multipliers. Its derivation requires the same technical conditions on the form of the DE system (2.14) as stated in Proposition 3. This well-know result [4,26,27,29] gives a precise relationship between adjoint-symmetries and multipliers, or equivalently between "nonlinear self-adjointness" and multipliers. In particular, it provides necessary and sufficient conditions for an adjoint-symmetry to be a multiplier. The simplest situation is when adjoint-symmetries of the lowest-order form Q α (x, u) are considered, which corresponds to "nonlinear self-adjointness" without differential substitutions. In this case, the only condition is equation (3.21a), which reduces to R Q a α + ∂Q a ∂u α = 0. This condition is, in general, non-trivial. (Unfortunately, some recent work [16] incorrectly asserts that, for any DE system, every adjoint-symmetry of the form Q α (x, u) is a multiplier.) When Theorem 3 is applied to variational DE systems, it yields the following well-known connection [4,26,27,29] with Noether's theorem.
Corollary 1. For a variational DE system (2.33), the multiplier determining system (3.20)-(3.21) reduces to a determining system for variational symmetries. In particular, the determining equation for adjoint-symmetries (3.20) coincides with the determining equation for symmetries (2.23), and the Helmholtz-type conditions (3.21) coincide with the necessary and sufficient conditions for a symmetry to be variational (namely, that prX(L) = D i Γ i holds for some differential vector function Γ i , where L is the Lagrangian (2.34)).
Note that, in this modern formulation of Noether's theorem, the use of a Lagrangian is completely by-passed through the Helmholtz-type conditions (3.21).
Example: For the semilinear wave equation (2.63), the determining equation for multipliers of lowest-order form Q C (t, x, u) is given by Since Q C (t, x, u) does not depend on derivatives of u, this determining equation splits with respect to the variables u t , u x , u tt , u xx , giving an overdetermined linear system which can be derived and solved directly by Maple. This provides the simplest computational route to finding all multipliers of lowest-order form. The connection between multipliers and adjointsymmetries arises when the determining equation (3.23) is instead split into the two terms δ * Q C F and δ * F Q C , which are given by where the operator R Q C is obtained from expression (3.14). Hence, on the solution space E of the wave equation (2.63), the multiplier determining equation reduces to the adjointsymmetry equation (2.67). Off of the solution space E, the multiplier determining equation which splits with respect to F , yielding ∂Q C ∂u = 0. The Helmholtz-type equation (3.27) directly shows that all adjoint-symmetries of the form Q(t, x) are conservation laws multipliers Q C (t, x), and so the three adjoint-symmetries (3.11) each determine a non-trivial conserved current through the characteristic equation These conserved currents (Ĉ t ,Ĉ x ) can be derived in terms of the multipliers Q C (t, x, u) in several different ways. One simple way is by applying integration by parts to the terms in Q C F to get a total time derivative D tĈ t plus a total space derivative D xĈ x , which yields (Ĉ t ,Ĉ x ). Another way is by takingĈ t (t, x, u, u t , u x ) andĈ x (t, x, u, u t , u x ) as unknowns and splitting the characteristic equation with respect to u tt , u tx , u xx to get a linear system of determining equations that can be integrated. The resulting conserved currents are shown in table 3. The specific relationship between these conserved currents and the conserved currents derived in table 2 will be explained in the next subsection.
3.2. Conservation laws produced by a multiplier/symmetry pair. From Theorem 3 and Proposition 3, every multiplier admitted by a given DE system determines, up to local equivalence, a conserved current for the system. Since multipliers are adjoint-symmetries, the adjoint-symmetry/symmetry formula (2.70) can be applied by using any multiplier (3.3) together with any symmetry (2.40). The resulting conserved current produced this way is given by where P α (x, u, ∂u, . . . , ∂ r u) is a given symmetry characteristic and Q a C (x, u, ∂u, . . . , ∂ s u) is a given multiplier. The following result characterizing these conserved currents will now be established for DE systems in the class stated in Proposition 3. The case of DE systems consisting of a single DE has appeared previously in Ref. [23].
Theorem 4. Let Ψ i (P, Q C ; F ) be the conserved current produced from the adjointsymmetry/symmetry formula (2.70) by using any multiplier Q a C (x, u, ∂u, . . . , ∂ s u) together with any symmetry characteristic P α (x, u, ∂u, . . . , ∂ r u). This conserved current Ψ i (P, Q C ; F ) is locally equivalent to a conserved current (2.30) that is given by the infinitesimal action of the symmetryX P = P α ∂/∂u α applied to the conserved current C i determined by the multiplier Q a C . In particular, Ψ i (P, Q C ; F ) and C i are related by for some differential antisymmetric tensor function Θ ij (x, u, ∂u, . . . , ∂ k u).
The proof consists of showing that both conserved currents Ψ i (P, Q C ; F ) and prX P (C i ) have the same multiplier. Consider the local conservation law determined by the multiplier Q a C . The symmetryX P applied to the characteristic equation (3.2) of this conservation law yields (3.31) The second term in this equation can be expressed as using the symmetry identity (2.41) combined with integration by parts, where Γ i (Q C , F ; P )| E = 0. Next, the first term in equation (3.31) can be expressed as through the Fréchet derivative identity (2.3) combined with the multiplier determining equation (3.18) and the adjoint-symmetry identity (2.61). Integration by parts then yields where Γ i (P, F ; Q C )| E = 0. Substitution of expressions (3.34) and (3.32) into equation (3.31) gives (3.35) Finally, since prX P commutes with total derivatives [24,29], this yields whereΓ i | E = 0 is a locally trivial conserved current, and where is the multiplier (3.9) of the local conservation law (3.7) from the adjointsymmetry/symmetry formula (2.70) with Q a = Q a C . This completes the proof. Theorem 4 is a generalization of a similar result [2,3,35] for variational DE systems, where the adjoint-symmetry/symmetry formula (2.70) reduces to a formula using any pair of symmetries.
Corollary 2. For a variational DE system, let Ψ i (P, Q C ; F ) be the conserved current produced from the adjoint-symmetry/symmetry formula (2.70) by using any symmetry characteristic P α (x, u, ∂u, . . . , ∂ r u) together with any multiplier Q α C (x, u, ∂u, . . . , ∂ s u) given by a variational symmetry characteristic. The conserved current Ψ i (P, Q C ; F ) is locally equivalent to a conserved current (2.30) that is given by the infinitesimal action of the symmetrŷ X P = P α ∂/∂u α applied to the conserved current C i determined by the multiplier Q α C . Moreover, through Noether's theorem, the multiplier of this conserved current Ψ i (P, Q C ; F ) is the characteristic of a variational symmetry given by the commutator of the symmetrieŝ X P = P α ∂/∂u α andX Q C = Q α C ∂/∂u α . Several basic properties of the adjoint-symmetry/symmetry formula (2.70) can be deduced from Theorem 4, as first shown in Ref. [23] for DE systems consisting of a single DE.
Theorem 5. (i) For a given DE system (2.14), let Q a C be the multiplier for a local conservation law in which the components of the conserved current C i have no explicit dependence on x. Then, using any translation symmetry X = a i ∂/∂x i , with characteristic P α = −a i u α i where a i is a constant vector, the conserved current Ψ i (P, Q C ; F ) is locally trivial. (ii) For a given DE system (2.14) that possesses a scaling symmetry X = a (i) where a (i) , b (α) are constants, let Q a C be the multiplier for a local conservation law in which the components of the conserved current C i are scaling homogeneous. Then, using the charac- i of the scaling symmetry, the conserved current Ψ i (P, Q C ; F ) is locally equivalent to a multiple w of the conserved current C i determined by Q a C . This multiple, w = const., is the scaling weight of the conserved integral given by ∂Ω C i dS i where Ω is any closed domain in R n and ∂Ω is its boundary surface.
The proof is a straightforward extension of the proof in Ref. [23] and will be omitted. Part (i) of this theorem explains the observations made in many recent papers in which Ibragimov's theorem gave only trivial local conservation laws. This will happen whenever the only local symmetries admitted by a DE system are translations and the only admitted adjoint-symmetries have no dependence on x.
Part (ii) of the theorem first appeared in Ref. [28]. It shows that the local conservation laws admitted by any DE system with a scaling symmetry can be obtained from an algebraic formula using the conservation law multipliers. This explains why in many recent papers the use of scaling symmetries in Ibragimov's theorem has produced non-trivial local conservation laws.
A more important point comes from putting together Theorem 4 and Theorem 3. Together, these two theorems show that the adjoint-symmetry/symmetry formula (2.70) cannot produce any "new" local conservation laws, since any local conservation law admitted by a given DE system must already arise directly from a multiplier. Moreover, for this formula to generate all of the local conservation laws for a given DE system, it seems plausible that the set of admitted symmetries needs to act transitively on set of admitted local conservation laws, so then every multiplier arises from some symmetry applied to some multiplier. The need for a transitive action is especially clear from Corollary 2, since if a pair of commuting variational symmetries is used in the formula, then the resulting local conservation law will have a trivial multiplier and hence will be a locally trivial conservation law. These significant deficiencies should discourage the unnecessary use of the adjointsymmetry/symmetry formula (2.70) -and consequently the unnecessary use of Ibragimov's theorem -when local conservation laws are being sought for a given DE system. It is much simpler and more direct to find all multipliers and then derive the conserved currents determined by these multipliers, as will be explained further in the next section.
Example: For the semilinear wave equation (2.63), table 2 shows the conserved currents obtained from the adjoint-symmetry/symmetry formula (2.70). Each of these conserved currents can be checked to satisfy the characteristic equation (3.28) withĈ t = Ψ t and C x = Ψ x , where the resulting multipliers Q Ψ are shown in table 4. There is a simple relationship (3.37) between each multiplier Q Ψ and the adjoint-symmetry/symmetry pair Q, P used to generate the conserved current (Ψ t , Ψ x ). In particular, from expressions (3.11)-(3.15) for the symmetry characteristics, adjoint-symmetries, and their associated operators R P and R Q , the relationship (3.37) yields in accordance with table 4. Table 4. Multipliers from the adjoint-symmetry/symmetry formula −u x e γx q e γx (±(q ′′ + γq ′ )u + (q ′ + γq)(u t + b du)), In particular, consider the case when m(u) is zero and both b(u), c(u) are arbitrary, so then the only admitted multiplier of lowest-order form Q(t, x, u) is Q = 1 (up to a multiplicative constant), as shown by table 1. In this case, it is straightforward to show that (by solving the relevant determining equations) there are no first-order multipliers and that the only admitted point symmetries are generated by the translations (2.86). Consequently, when the set of multipliers Q(t, x, u, u t , u x ) is considered, a single non-trivial conservation law C t = u t + ∫ b(u) du, C x = −u x + ∫ c(u) du is admitted by the wave equation u tt − u xx + b(u)u t + c(u)u x = 0 with b(u) and c(u) arbitrary, whereas all of the conserved currents (Ψ t , Ψ x ) obtained from Ibragimov's theorem (2.82), or from the simpler equivalent adjoint-symmetry/symmetry formula (2.83), are trivial! Note that, correspondingly, the symmetry action on the set of non-trivial conservation laws given by the set of multipliers Q(t, x, u, u t , u x ) is not transitive. This example succinctly illustrates the incompleteness of these formulas for generating conservation laws.

A direct construction method to find all local conservation laws
The results stated in Proposition 3 and Theorems 3, 4, and 5 have been developed in Ref. [4,26,27,28,23] and extended in Ref. [36,29]. This collective work provides a simple, algorithmic method to find all local conservation laws for any given system of DEs. The method is based on the general result that all local conservation laws arise from multipliers as given by the solutions of a linear system of determining equations, where the multipliers are simply adjoint-symmetries subject to certain Helmholtz-type conditions.
Consequently, all multipliers can be found by either of the two following methods [4]: (1) directly solve the full determining system for multipliers; or (2) first, solve the determining equation for adjoint-symmetries, and next, check which of the adjoint-symmetries satisfy the Helmholtz-type conditions. The adjoint-symmetry determining equation is simply the adjoint of the symmetry determining equation, and hence it can be solved by the standard algorithmic procedure used for solving the symmetry determining equation [24,32,31]. Likewise, the same procedure works equally well for solving the multiplier determining system.
A natural question is, in practice, at which differential orders s ≥ 0 will multipliers or adjoint-symmetries Q a (x, u, ∂u, . . . , ∂ s u) be found?
One answer is that the same situation arises for symmetries. Normally, point symmetries are sought first, since many DE systems admit point symmetries, and since relatively fewer DE systems admit contact symmetries or higher-order symmetries. Indeed, the existence of a sufficiently high-order symmetry is one main definition of an integrable system [37], as this can indicate the existence of an infinite hierarchy of successively higher-order symmetries. For multipliers, the most physically important conserved currents always have a low differential order. Based on numerous examples, a concrete definition of a low-order multiplier that seems to characterize these physically important conserved currents, and distinguishes them from higher-order conserved currents arising for integrable systems, has been introduced in recent work [36,29].
Another answer is that it is straightforward just to find all multipliers or adjointsymmetries with a specified differential order s = 0, 1, 2, . . ., going up to any desired maximum finite order. Moreover, in some situations a standard descent/induction argument [38,39,40] can be used to find the multipliers or adjoint-symmetries to all orders s ≥ 0.
Once a set of multipliers has been found for a given DE system, the corresponding conserved currents are straightforward to find in an explicit form. Several different methods are available.
One algorithmic method is the direct integration of the characteristic equation [31,41] defining the conserved current. Another algorithmic method is the use of a homotopy integral formula. This method has several versions [24,4,26,27,29], all of which involve trade-offs between the simplicity of the integration versus the flexibility of avoiding singularities (if any) in the integrand.
However, purely algebraic methods for the construction of conserved currents from multipliers are known. One algebraic method is the use of a scaling formula [28,23,31] which is given by the adjoint-symmetry/symmetry formula. This applies only to DE systems that admit a scaling symmetry, but it has recently been extended to general DE systems by incorporating a dimensional analysis method as shown in Ref. [29]. In particular, with the use of this dimensional-scaling method, the construction of conserved currents becomes completely algebraic.
Therefore, the general method just outlined provides a completely algorithmic computational way to derive all local conservation laws for any given DE system. In particular, there is no need to resort to any special methods or ansatzes (such as the "abc" technique [42], partial Lagrangians [43], "nonlinear self-adjointness" [9,12,13,15,17], undetermined coefficients [44]) which at best just yield a subset of all of the local conservation laws admitted by a DE system or just apply to restricted classes of DE systems.
Example: The semilinear wave equation (2.63) can be expected to admit conserved currents that depend nonlinearly on u t and u x , in addition to the previous conserved currents in table 3, all of which have linear dependence on u t and u x . The multipliers (3.11) for the latter conserved currents have the form Q C (t, x). Conserved currents that depend nonlinearly on u t and u x will arise from multipliers Q C (t, x, u, u t , u x ) that have explicit dependence on u t and u x . It is straightforward to find all such multipliers by using Maple to set up and solve the multiplier equation (3.23), which splits respect to the variables u tx , u tt , u xx , u txx , u ttx , u txx , u ttt , u xxx , giving an overdetermined linear system. Alternatively, the multiplier equation (3.23) can be split instead into the two terms δ * Q C F and δ * F Q C , which provides a direct connection between multipliers and adjoint-symmetries. In particular, the first term in the multiplier equation (3.23) consists of where the operator R Q C is found to be given by for multipliers Q C (t, x, u, u t , u x ), through u tt = u xx −b(u)u t −c(u)u x −m(u) and its differential consequences. The second term in the multiplier equation (3.23) is given by for multipliers Q C (t, x, u, u t , u x ). On the solution space E of the wave equation (2.63), the terms (4.3) vanish, while the other terms (4.1) reduce to the adjoint-symmetry equation (2.67). Off of the solution space E, the these terms (4.1) and (4.3) become a linear combination of F , D t F , D x F , whose coefficients must vanish separately. This splitting is found to yield a single Helmholtz-type equation Taken together, this Helmholtz-type equation (4.5) and the adjoint-symmetry equation (2.67) constitute the determining system (3.20)-(3.21) for finding all first-order multipliers Q C (t, x, u, u t , u x ) admitted by the wave equation (2.63).
The most computationally effective way to solve equations (4.5) and (2.67) in the determining system is by changing variables from t, x, u, u t , u x to µ = 1 2 (t + x), ν = 1 2 (t − x), u, u µ = u t + u x , u ν = u t − u x , based on null coordinates for the wave equation (2.63). In these new variables, the general solution of the determining system consists of three distinct cases (as obtained using the Maple package 'rifsimp'), after the nonlinearity and homogeneity conditions (2.69) are imposed on b(u), c(u), d(u). The resulting multipliers, after merging cases, are shown in table 5. Each multiplier determines a non-trivial conserved current through the characteristic equation (3.28). These conserved currents (Ĉ t ,Ĉ x ) can be derived in terms of the multipliers Q C (t, x, u, u t , u x ) in the same way discussed previously for lowest-order multipliers. The results are shown in table 6. Table 5. First-order multipliers Table 6. First-order conserved currents × γ (b + c) du + 2(u t + u x ) + m 1

Concluding remarks
The conservation law theorem stated by Ibragimov in Ref. [9,12] for "nonlinear selfadjoint" DEs and subsequent extensions of this theorem in Ref. [13,15,17] are not new. In its most general form, this theorem is simply a re-writing of a standard formula [2,3,4] that uses a pair consisting of a symmetry and an adjoint-symmetry to produce a conservation law through a well-known Fréchet derivative identity [2,3,24,31,29]. Unfortunately, no references to prior literature are provided in Ibragimov's papers, which may give the impression that the results are original. One aspect that is novel is the derivation of the formula by using an auxiliary Lagrangian, although it does not in any way simplify either the formula or its content. Moreover, the condition of "nonlinear self-adjointness" is nothing but a re-writing of the condition that a DE system admits an adjoint-symmetry [4,29], and this condition automatically holds for any DE system that admits a local conservation law.
The present paper shows how the symmetry/adjoint-symmetry formula is directly connected to the action of symmetries on conservation laws, which explains a number of major drawbacks in trying to use the formula -and hence in applying Ibragimov's theoremas a method to generate conservation laws. In particular, the formula can generate trivial conservation laws and does not always yield all non-trivial conservation laws unless the symmetry action on the set of these conservation laws is transitive, which cannot be known until all conservation laws have been found.
A broader point, which is more important, is that there is a completely general method [31,29] using adjoint-symmetries [2,3,4,26,27] to find all local conservation laws for any given DE system. This method is a kind of adjoint version of the standard Lie method to find all local symmetries. The method is algorithmic [29] and the required computations are no more difficult than the computations used to find local symmetries.