Next Article in Journal
Fuzzy Goal Programming with an Imprecise Intuitionistic Fuzzy Preference Relations
Next Article in Special Issue
Nonlocal Conservation Laws of PDEs Possessing Differential Coverings
Previous Article in Journal
Indirect Contributions to Tumor Dynamics in the First Stage of the Avascular Phase
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geometrical Formulation for Adjoint-Symmetries of Partial Differential Equations

Department of Mathematics and Statistics, Brock University, St. Catharines, ON L2S3A1, Canada
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: Department of Mathematics, Ningbo University, Ningbo 315211, China.
Symmetry 2020, 12(9), 1547; https://doi.org/10.3390/sym12091547
Submission received: 3 September 2020 / Revised: 12 September 2020 / Accepted: 14 September 2020 / Published: 19 September 2020
(This article belongs to the Special Issue Geometric Analysis of Nonlinear Partial Differential Equations)

Abstract

:
A geometrical formulation for adjoint-symmetries as one-forms is studied for general partial differential equations (PDEs), which provides a dual counterpart of the geometrical meaning of symmetries as tangent vector fields on the solution space of a PDE. Two applications of this formulation are presented. Additionally, for systems of evolution equations, adjoint-symmetries are shown to have another geometrical formulation given by one-forms that are invariant under the flow generated by the system on the solution space. This result is generalized to systems of evolution equations with spatial constraints, where adjoint-symmetry one-forms are shown to be invariant up to a functional multiplier of a normal one-form associated with the constraint equations. All of the results are applicable to the PDE systems of interest in applied mathematics and mathematical physics.

1. Introduction

Symmetries are a fundamental coordinate-free structure of a partial differential equation (PDE). In geometrical terms, an infinitesimal symmetry is an evolutionary (vertical) vector field that is tangent to the solution space of a PDE, where the components of the vector field are the solutions of the linearization of the PDE on its solution space (see, e.g., [1,2,3,4]).
Knowledge of the symmetries of a PDE can be used to map given solutions into other solutions, find invariant solutions, detect and find mappings in a target class of PDEs, detect integrability, and find conservation laws through Noether’s theorem when a PDE has a variational (Lagrangian) structure.
Solutions of the adjoint linearization of a PDE on its solution space are known as adjoint-symmetries. This terminology was first introduced and explored for ordinary differential equations (ODEs) in [5,6,7,8] and then generalized to PDEs in [9,10] (see [11] for a recent overview for PDEs). When a PDE lacks a variation structure, then its adjoint-symmetries will differ from its symmetries.
Knowledge of the adjoint-symmetries of a PDE can be used for several purposes just as symmetries can. Specifically, solutions of the PDE can be found analogously to the invariant surface condition associated with a symmetry; mappings into a target class of PDEs can be detected and found analogously to characterizing the symmetry structure of the target class; integrability can be detected analogously to the existence of higher order symmetries; and conservation laws can be determined analogously to symmetries that satisfy a variational condition. In particular, the counterpart of variational symmetries for a general PDE is provided by multipliers, which are well known to be adjoint-symmetries that satisfy a Euler–Lagrange condition.
However, a simple geometrical meaning (apart from abstract formulations) for adjoint-symmetries has yet to be developed in general for PDEs. Several significant new steps toward this goal will be taken in the present paper.
Firstly, for general PDE systems, adjoint-symmetries will be shown to correspond to evolutionary (vertical) one-forms that functionally vanish on the solution space of the system. This formulation has two interesting applications. It will provide a geometrical derivation of a well-known formula that generates a conservation law from a pair consisting of a symmetry and an adjoint-symmetry [9,12]. It also will yield three different actions of symmetries on adjoint-symmetries from Cartan’s formula for the Lie derivative, providing a geometrical formulation of some recent work that used an algebraic viewpoint [13].
Secondly, for evolution systems, these adjoint-symmetry one-forms will be shown to have the structure of a Lie derivative of a simpler underlying one-form, utilizing the flow generated by the system. As a result, adjoint-symmetries of evolution systems will geometrically correspond to one-forms that are invariant under the flow on the solution space of the system. This directly generalizes the geometrical meaning of adjoint-symmetries known for ODEs [8].
Thirdly, a bridge between the preceding results for general PDE systems and evolution systems will be developed by considering evolution systems with spatial constraints. These systems are ubiquitous in applied mathematics and mathematical physics, for example: Maxwell’s equations, incompressible fluid equations, magnetohydrodynamical equations, and Einstein’s equations. For such systems, invariance of the adjoint-symmetry one-form under the constrained flow will be shown to hold up to a functional multiple of the normal one-form associated with the constraint equations.
Throughout, the approach will be concrete, rather than abstract, so that the results can be readily understood and applied to specific PDE systems of interest in applied mathematics and mathematical physics.
The rest of the paper is organized as follows. Section 2 discusses the evolutionary form of vector fields and its counterpart for one-forms in the mathematical framework of calculus in jet space, which will underlie all of the main results. Section 3 reviews the geometrical formulation of symmetries and presents the counterpart geometrical formulation of adjoint-symmetries. In addition, some examples of adjoint-symmetries of physically interesting PDE systems are discussed. Section 4 gives the two applications of adjoint-symmetry one-forms. Section 5 develops the main results for adjoint-symmetries of evolution systems and extends these results to constrained evolution systems. Some concluding remarks are made in Section 6.

2. Vector Fields, One-Form Fields, and Their Evolutionary Form

To begin, some essential tools [3,11,14] from calculus in jet space will be reviewed. This will set the stage for a discussion of the evolutionary form of vector fields and its counterpart for one-forms, as needed for the main results in the subsequent sections.
Independent variables are denoted x i , i = 1 , , n , and dependent variables are denoted u α , α = 1 , , m . Derivative variables are indicated by subscripts employing a multi-index notation: I = { i 1 , , i N } , u I α = u i 1 i N α : = x i 1 x i N u α , | I | = N ; I = , u I α : = u α , | I | = 0 . Some useful notation is as follows: k u will denote the set { u I α } | I | = k of all derivative variables of order k 0 ; u ( k ) will denote the set { u I α } 0 | I | k of all derivative variables of all orders up to k 0 . The summation convention of summing over any repeated (multi-)index in an expression is used throughout.
Jet space is the coordinate space J = ( x i , u α , u j α , ) . A smooth function u α = ϕ α ( x ) : R n R m determines a point in J : at any x i = ( x 0 ) i ; the values ( u 0 ) α : = ϕ α ( x 0 ) and the derivative values ( u 0 ) J α : = j 1 j N ϕ α ( x 0 ) for all orders N 1 give a map,
u α = ϕ α ( x ) x 0 ( ( x 0 ) i , ( u 0 ) α , ( u 0 ) j α , ) J .
In jet space, the primitive geometric objects consist of partial derivatives x i , u J α , and differentials d x i , d u J α . They are related by duality (hooking) relations:
x i d x j = δ i j ,
u I α d u J β = δ α β δ J I .
It will be useful to also introduce the geometric contact one-forms:
Θ I α = d u I α u I i α d x i .
Under the evaluation map (1), the pull back of a contact one-form vanishes.
Total derivatives are given by D i = x i + u i J α u J α , which corresponds to the chain rule under the evaluation map (1). Higher total derivatives are defined by D J = D j 1 D j N , J = { j 1 , , j N } , | J | = N . For J = , D = id is the identity operator, where | | = 0 . In particular, D J u α = u J α , and D J d u α = d u J α .
A differential function is a function f ( x , u ( k ) ) defined on a finite jet space J ( k ) = ( x i , u α , u j α , , u j 1 j k α ) of order k 0 . The Frechet derivative of a differential function f is given by
f = f u I α D I
which acts on (differential) functions F α . The adjoint-Frechet derivative of a differential function f is given by
( f * ) α = ( 1 ) | I | D I f u I α
which acts on (differential) functions F, where the right-hand side is viewed as a composition of operators.
The Frechet second-derivative is given by
f ( F 1 , F 2 ) = f u I α u J β ( D I F 1 α ) ( D J F 2 β ) .
This expression is symmetric in the pair of functions ( F 1 α , F 2 α ) .
The commutator of two differential functions f 1 and f 2 is given by [ f 1 , f 2 ] = f 2 ( f 1 ) f 1 ( f 2 ) .
The Euler operator (variational derivative) is given by
E u α = ( 1 ) | I | D I u I α .
It characterizes total divergence expressions: E u α ( f ) = 0 holds identically iff f = D i F i for some differential vector function F i ( x , u ( k ) ) . The product rule takes the form:
E u α ( f 1 f 2 ) = f 1 * ( f 2 ) α + f 2 * ( f 1 ) α .
The higher Euler operators
E u α I = I J ( 1 ) | J | D J u I J α
characterize higher order total derivative expressions: E u α I ( f ) = 0 holds identically iff f = D i 1 D i | I | F i 1 i | I | for some differential tensor function F i 1 i | I | ( x , u ( k ) ) .
The Frechet derivative is related to the Euler operator by:
f ( F ) = F α E u α ( f ) + D i Γ i ( F ; f ) , Γ i ( F ; f ) = ( D J F α ) E u i J α ( f ) .
The Frechet derivative and its adjoint are related by
F 2 f ( F 1 ) F 1 α f * ( F 2 ) α = D i Ψ i ( F 1 , F 2 ; f ) , Ψ i ( F 1 , F 2 ; f ) = ( D K F 2 ) ( D J F 1 α ) E u i J α K ( f ) .

Evolutionary Vector Fields and One-Form Fields

A vector field in jet space is defined as the geometric object,
P i x i + P I α u I α
whose components are differential functions. Similarly, a one-form field in jet space is defined as the geometric object,
Q i d x i + Q α I d u I α
whose components are differential functions. Total derivatives D i = x i + u i I α u I α represent trivial vector fields that annihilate contact one-forms: D i Θ J α = 0 .
Geometric counterparts of partial derivatives u J α are evolutionary (vertical) differentials d u J α , where d is the evolutionary version of d: d 2 = 0 , d x i = 0 . They satisfy the duality (hooking) relation:
u I α d u J β = δ α β δ J I .
An evolutionary (vertical) vector field is the geometric object
P I α u I α
whose components are differential functions. Every vector field X = P i x i + P I α u I α has a unique evolutionary form X ^ = X P i D i = P ^ I α u I α given by the components P ^ I α = P I α P i u i I α . Its dual counterpart is an evolutionary (vertical) one-form field,
Q α I d u I α
whose components are differential functions.
For later developments, it will be useful to define the functional pairing relation,
P I α u I α , Q α I d u I α = P I α Q α I d x
between evolutionary vector fields and evolutionary one-form fields. In the local form, this pairing is given by the expression:
P I α Q α I mod total D .
Two evolutionary one-forms will be considered functionally equivalent iff their pairings with an arbitrary evolutionary vector field agree,
P I α u I α , Q 1 β J d u J β = P I α u I α , Q 2 β J d u J β ,
or in the local form,
P I α ( Q 1 α I Q 2 α I ) = 0 mod total D .
The functional equivalence of one-forms is closely related to the notion of functional one-forms in the variational bi-complex. See [3] for details.

3. Geometric Formulation of Symmetries and Adjoint-Symmetries

Consider a general PDE system of order N consisting of M equations,
G A ( x , u ( N ) ) = 0 , A = 1 , , M
where x i , i = 1 , , n , are the independent variables and u α , α = 1 , , m , are the dependent variables. The space of formal solutions u α ( x ) of the PDE system will be denoted E .
There are many equivalent starting points for the formulation of infinitesimal symmetries. For the present purpose, the most useful one is given by evolutionary vector fields and utilizes only the Frechet derivative. A symmetry is a vector field,
X P = P α ( x , u ( k ) ) u α
whose component functions P α ( x , u ( k ) ) are non-singular on E and satisfy the linearization of the PDE system on E ,
( pr X P G A ) | E = G ( P ) A | E = 0 .
This is the symmetry determining equation, and the functions P α are called the characteristic of the symmetry.
In this setting, an adjoint-symmetry consists of functions Q A ( x , u ( l ) ) that are non-singular on E and that satisfy the adjoint linearization of the PDE system on E ,
G * ( Q ) α | E = 0 .
This is the adjoint-symmetry determining equation.
In particular, the two determining equations (24) and (25) are formal adjoints of each other. They coincide only in two cases: either G = G * , which is the necessary and sufficient condition for a PDE system to be a Euler–Lagrange equation (namely, possess a variational structure) [1,3,11]; or G = G * , which is the necessary and sufficient condition for a PDE system to be a linear, constant-coefficient system of odd order [10].
Since P α has the geometrical status as the components of the vector field (23), a natural question is whether Q A has any status given by the components of some other geometrical object [11,12].
It will be useful to work with a coordinate-free description of the PDE system (22) in jet space. Such a system of equations ( G 1 ( x , u ( N ) ) , , G M ( x , u ( N ) ) ) = 0 describes a set of M surfaces in the finite space J ( N ) ( x , u , u , , N u ) . Total derivatives of these equations, ( D I G 1 ( x , u ( N ) ) , , D I G M ( x , u ( N ) ) ) = 0 , correspondingly describe sets of surfaces in the higher derivative finite spaces J ( N + | I | ) ( x , u , u , , N + | I | u ) . Altogether, the set comprised by the equations and the derivative equations for all orders | I | 0 corresponds to an infinite set of surfaces in jet space, which can be identified with the solution space E .
As is well known, symmetry vector fields geometrically describe tangent vector fields with respect to E . To see this explicitly, first consider the identities:
d G A = ( G A ) u I α d u I α ,
G ( P ) A = pr X P G A = pr X P d G A .
Now, observe that d G A is the normal one-form to the surfaces G A = 0 . The symmetry determining equation (24) then shows that the prolonged vector field pr X P is annihilated by the normal one-form and hence is tangent to these surfaces iff X P is a symmetry of the PDE system.
This normal one-form (26) provides a natural way to associate a one-form to an adjoint-symmetry via:
ϖ Q = Q A ( x , u ( l ) ) d G A .
A functionally equivalent one-form is obtained through integration by parts:
Q A d G A = Q A ( G A ) ( d u ) = G * ( Q ) α d u α mod total D .
Evaluating this one-form on the solution space E then gives
ϖ Q | E = 0 mod total D .
Thus, a one-form ϖ Q functionally vanishes on the surfaces E iff its components Q A are an adjoint-symmetry.
This establishes a main geometrical result.
Theorem 1.
Adjoint-symmetries describe evolutionary one-forms Q A d G A that functionally vanish on the solution space E of a PDE system (22).
These developments have used evolutionary (vertical) vector fields and evolutionary one-forms. It is straightforward to reformulate everything in terms of full vector fields and full one-forms.
First, consider the normal one-form
d G A = ( G A ) x i d x i + ( G A ) ( d u ) = ( G A ) ( Θ ) + ( ( G A ) x i + ( G A ) ( u i ) ) d x i = ( G A ) ( Θ ) + D i G A d x i
which yields the relation
d G A | E = ( G A ) ( Θ ) | E .
Then, observe:
Q A d G A | E = Q A ( G A ) ( Θ ) | E = ( G A ) * ( Q A ) α | E Θ α mod total D .
As a consequence, Q A d G A | E vanishes mod total D iff Q A satisfies the adjoint-symmetry determining Equation (25). Moreover, the determining equation itself can be expressed directly in terms of the one-form Q A d G A | E by E Θ α ( Q A d G A ) | E = ( G A ) * ( Q A ) | E = 0 .
Proposition 1.
The adjoint-symmetry determining Equation (25) can be expressed geometrically as:
E Θ α ( Q A d G A ) | E = 0 .

Examples of Adjoint-Symmetries

To illustrate the results, some examples of PDEs that possess non-trivial adjoint-symmetries will be given.
The Korteweg–de Vries (KdV) equation
u t + u u x + u x x x = 0
for shallow water waves is an example of an evolutionary wave equation. Its symmetries X = P u are the solutions of the determining equation
G ( P ) | E = ( D t P + D x ( u P ) + D x 3 P ) | E = 0 ,
with G = D t + D x u + D x 3 being the Frechet derivative of the KdV equation, where P is a non-singular function of t, x, u, and derivatives of u on the space of KdV solutions E . The determining equation for adjoint-symmetries ϖ = Q G ( d u ) is the adjoint equation
G * ( Q ) | E = ( D t Q u D x Q D x 3 Q ) | E = 0 ,
where Q is a non-singular function of t, x, u, and derivatives of u on E .
KdV adjoint-symmetries up to first-order Q ( t , x , u , u t , u x ) are given by [9] the span of
Q ( 1 ) = 1 , Q ( 2 ) = u , Q ( 3 ) = t u x .
The first two are part of a hierarchy of higher order adjoint-symmetries generated by a recursion operator R = D x 2 + 1 3 u + 1 3 D x 1 u D x applied to Q = u . The third one along with all of the ones in the hierarchy are related to symmetries of the KdV equation through the Hamiltonian operator H = D x . If a linear combination of the lowest order adjoint-symmetries is used like an invariant surface condition, c 1 + c 2 ( t u x ) + c 3 u = 0 , then this yields u = ( c 2 x c 1 ) / ( c 2 t + c 3 ) , which is a similarity solution of the KdV equation.
An example of a non-evolutionary equation is,
Δ ϕ t + ϕ x Δ ϕ y ϕ y Δ ϕ x = 0
which governs the vorticity Ω = Δ ϕ for incompressible inviscid fluid flow in two spatial dimensions, where the fluid velocity has the components v = ( ϕ y , ϕ x ) . The symmetries X = P ϕ of this equation are the solutions of the determining equation,
G ( P ) | E = ( D t Δ P + ϕ x D y Δ P + Δ ϕ y D x P ϕ y D x Δ P Δ ϕ x D y P ) | E = 0 ,
where P is a non-singular function of t, x, y, ϕ , and derivatives of ϕ on the space of vorticity solutions E , with G = D t Δ + ϕ x D y Δ + Δ ϕ y D x ϕ y D x Δ Δ ϕ x D y being the Frechet derivative of the vorticity equation given in terms of the total Laplacian operator Δ = D x 2 + D y 2 . The determining equation for adjoint-symmetries ϖ = Q G ( d ϕ ) is the adjoint equation,
G * ( Q ) | E = ( D t Δ Q + D y Δ ( ϕ x Q ) + D x ( Δ ϕ y Q ) D x Δ ( ϕ y Q ) D y ( Δ ϕ x Q ) | E = 0 ,
where Q is a non-singular function of t, x, y, ϕ , and derivatives of ϕ on E .
The first-order adjoint-symmetries Q ( t , x , y , ϕ , ϕ t , ϕ x , ϕ y ) are given by [13] the span of,
Q ( 1 ) = x 2 + y 2 , Q ( 2 ) = ϕ , Q ( 3 ) = f ( t ) , Q ( 4 ) = x f ( t ) , Q ( 5 ) = y f ( t ) ,
where f ( t ) is an arbitrary smooth function. If a linear combination of these adjoint-symmetries is used like an invariant surface condition, c 1 ( x 2 + y 2 ) + c 2 ϕ + c 3 f ( t ) + c 4 x f ( t ) + c 5 y f ( t ) = 0 , then taking c 2 = 1 gives ϕ = c 1 ( x 2 + y 2 ) + ( c 3 + c 4 x + c 5 y ) f ( t ) , which is a constant vorticity solution, with Ω = 2 c 1 and v = ( 2 c 1 y + c 5 f ( t ) , 2 c 1 x + c 4 f ( t ) ) .
Maxwell’s equations in free space are an example of an evolution system with spatial constraints:
E t × B = 0 , B t + × E = 0 , · E = · B = 0
(in relativistic units with the speed of light set to one). The symmetries X = P E · E + P B · B of this system are the solutions of the determining equations
G P E P B | E = ( D t P E × P B ) | E ( D t P B + × P E ) | E ( · P E ) | E ( · P B ) | E = 0 ,
where P E and P B are non-singular vector functions of t, x, y, z, E , B , and derivatives of E , B on the space of Maxwell solutions E , with G = D t × × D t · 0 0 · being the Frechet derivative of the system in terms of the total derivative operator = ( D x , D y , D z ) . The determining equation for adjoint-symmetries ϖ = Q E Q B Q E Q B G d E d B is the adjoint equation
G * Q E Q B Q E Q B | E = ( D t Q E + × Q B Q E ) | E ( D t Q B × Q E Q B ) | E = 0 ,
where the vectors Q E , Q B , and the scalars Q E , Q B , are non-singular functions of t, x, y, z, E , B , and derivatives of E , B on E . Note that the adjoint * here includes a matrix transpose applied to the row matrix comprising the adjoint-symmetry vector and scalar functions.
Because Maxwell’s equations are a linear system and contain constraints, it possesses three types of adjoint-symmetries [15,16]: elementary adjoint-symmetries such that Q E , Q B , Q E , Q B are functions only of t, x, y, z; gauge adjoint-symmetries given by Q E = χ E , Q B = χ B , Q E = D t χ E , Q B = D t χ B in terms of scalars χ E and χ B that are arbitrary non-singular functions of t, x, y, z, E , B , and derivatives of E , B on E ; and a hierarchy of linear adjoint-symmetries. The linear adjoint-symmetries of zeroth order are given by the span of
Q E = ξ × B + ζ E , Q B = ξ × E + ζ B , Q E = ξ · E , Q B = ξ · B
and
Q E = ξ × E ζ B , Q B = ξ × B + ζ E , Q E = ξ · B , Q B = ξ · E
where
ξ = a 0 + a 1 × x + a 2 t + a 3 x + a 4 t x + ( a 5 · x ) x 1 2 a 5 ( x · x + t 2 ) , ζ = a 0 + a 2 · x + a 3 t + 1 2 a 4 ( x · x + t 2 ) + ( a 5 · x ) t ,
in terms of arbitrary constant scalars a 0 , a 3 , a 4 and arbitrary constant vectors a 0 , a 1 , a 2 , a 5 , with x = ( x , y , z ) . The pair ( ξ , ζ ) represents a conformal Killing vector in Minkowski space R 3 , 1 .
These two zeroth-order adjoint-symmetries are related by the duality symmetry ( E , B ) ( B , E ) . The linear first-order adjoint-symmetries are more complicated and involve conformal Killing–Yano tensors. All higher order adjoint-symmetries can be obtained from the zeroth and first-order adjoint-symmetries by taking Lie derivatives with respect to conformal Killing vectors. Their explicit description can be found in [15,16]. An unexplored question is whether the lowest order adjoint-symmetries can be used like an invariant surface condition to produce solutions of Maxwell’s equations.

4. Some Applications

Two geometrical applications of Theorem 1 will be presented. The first application is a geometrical derivation of a well-known formula that generates a conservation law from a pair consisting of a symmetry and an adjoint-symmetry. This derivation will use the functional pairing (18). The second application is a geometrical derivation of three actions of symmetries on adjoint-symmetries. These symmetry actions have been obtained in recent work using an algebraic point of view [13]. They will be shown here to arise from Cartan’s formula for the Lie derivative of an adjoint-symmetry one-form (28).
It will be useful to work with the determining equations for symmetries and adjoint-symmetries off of the solution space E of a given PDE system (22). More precisely, the determining equations will be expressed in the full jet space containing E .
Remark 1.
A PDE system (22) will be assumed to be regular [11], so that Hadamard’s lemma holds: a differential function f satisfies f | E = 0 iff f = R f ( G ) , where R f is a linear differential operator whose coefficients are non-singular on E .
Consequently, for symmetries, G ( P ) A | E = 0 holds iff
G ( P ) A = R P ( G ) A ,
and likewise for adjoint-symmetries, G * ( Q ) α | E = 0 holds iff
G * ( Q ) α = R Q ( G ) α ,
where R P and R Q are linear differential operators whose coefficients are non-singular on E .

4.1. Conservation Laws from Symmetries and Adjoint-Symmetries

The functional pairing (18) between a symmetry vector field (23) and an adjoint-symmetry one-form (28) is given by,
pr X P , ϖ Q = pr P α u α , Q A d G A = Q A G ( P ) A d x
from identity (27). This pairing in local form (19) is the expression,
Q A G ( P ) A mod total D .
There are two different ways to evaluate it.
First, since X P is a symmetry, Q A G ( P ) A = Q A R P ( G ) A . Second, since ϖ Q is an adjoint-symmetry, Q A G ( P ) A = G * ( Q ) α P α + D i Ψ i ( P , Q ) G = P α R Q ( G ) α + D i Ψ i ( P , Q ; G ) , where
Ψ i ( P , Q ; G ) = ( D K Q A ) ( D J P α ) E u i J α K ( G A ) .
Hence, on E , Q A G ( P ) A | E = D i Ψ i ( P , Q ) G | E = 0 , which is equivalent to pr X P , ϖ Q | E = 0 . This establishes the following conservation law.
Theorem 2.
Vanishing of the functional pairing (51) for any symmetry (23) and any adjoint-symmetry (28) corresponds to a conservation law
D i Ψ i ( P , Q ; G ) | E = 0
holding for the PDE system G A = 0 , where the conserved current Ψ i ( P , Q ; G ) is given by expression (53).

4.2. Action of Symmetries on Adjoint-Symmetries

For any PDE system (22), its set of adjoint-symmetries is a linear space, and as shown in [13], symmetries of the PDE system have three different actions on this space.
The primary symmetry action can be derived from the Lie derivative of an adjoint-symmetry one-form with respect to a symmetry vector field.
Proposition 2.
If ϖ Q is an adjoint-symmetry one-form (28), namely ϖ Q | E = 0 ( mod total D ) , then its Lie derivative with respect to any symmetry vector X P = P α u α yields an adjoint-symmetry one-form,
L X P ϖ Q | E = ϖ S P ( Q ) | E = 0 ( mod total D )
where
S P ( Q ) A = Q ( P ) A + R P * ( Q ) A
are its components.
Here and throughout, R P and R Q are the linear differential operators determined by Equations (49) and (50). The adjoints of these operators are denoted R P * and R Q * .
Proof. 
Recall that the Lie derivative has the following properties: it acts as a derivation; it commutes with the differential d; it reduces to the Frechet derivative when acting on a differential function.
By the use of these properties,
L X P ϖ Q = L X P ( Q A d G A ) = ( L X P Q A ) d G A + Q A L X P ( d G A ) = Q ( P ) A d G A + Q A d ( G ( P ) A ) = Q ( P ) A d G A + Q A d ( R P ( G ) A ) .
The last term can be simplified on E : Q A d ( R P ( G ) A ) | E = Q A R P ( d G ) A | E = R P * ( Q ) A d G A ( mod total D ) . This yields
L X P ϖ Q | E = ( ( Q ( P ) A + R P * ( Q ) A ) d G A ) | E ( mod total D ) ,
completing the derivation. □
There is an elegant formula, due to Cartan, for the Lie derivative in terms of the operations d and . This formula gives rise to two additional symmetry actions.
Theorem 3.
The terms in Cartan’s formula
L X P ϖ Q = d ( pr X P ϖ Q ) + pr X P ( d ϖ Q )
evaluated on E each yield an action of symmetries on adjoint symmetries. The action produced by the Lie derivative term has the components (56), and the actions produced by the differential term and the hook term respectively have the components
S 1 P ( Q ) = R P * ( Q ) A R Q * ( P ) A ,
S 2 P ( Q ) = Q ( P ) A + R Q * ( P ) A .
Proof. 
Consider the first term on right-hand side in the formula (59). It can be evaluated in two different ways. Firstly, pr X P ( Q A d G A ) = Q A G ( P ) A = Q A R P ( G ) A yields
d ( pr X P ( Q A d G A ) ) | E = d ( Q A R P ( G ) A ) | E = ( Q A R P ( d G A ) ) | E = ( R P * ( Q ) A d G A ) | E .
Secondly, Q A d G A = R Q ( G ) α Θ α + Q A ( D i G A ) d x i ( mod total D ) gives pr X P ( Q A d G A ) = pr X P ( R Q ( G ) α Θ α + Q A ( D i G A ) d x i ( mod total D ) ) = R Q ( G ) α P α ( mod total D ) . This yields
d ( pr X P ( Q A d G A ) ) | E = d ( R Q ( G ) α P α ( mod total D ) ) | E = ( R Q ( d G ) α P α ( mod total D ) ) | E = ( R Q * ( P ) A d G A ( mod total D ) ) | E .
Then, equating expressions (62) and (63) leads to the result:
( ( R P * ( Q ) A R Q * ( P ) A ) d G A ) | E = 0 ( mod total D ) | E .
This equation shows that the symmetry action (60) produces an adjoint-symmetry.
Now, consider the second term on the right-hand side in formula (59). Similarly to the first term, it can be evaluated in two different ways. Firstly, d ϖ Q = d Q A d G A yields
pr X P ( d Q A d G A ) = Q ( P ) A d G A G ( P ) A d Q A = Q ( P ) A d G A R P ( G ) A d Q A .
Hence, on E ,
( pr X P ( d Q A d G A ) ) | E = ( Q ( P ) A d G A ) | E .
Secondly, d ϖ Q = d ( R Q ( G ) α Θ α + Q A ( D i G A ) d x i ) ( mod total D ) gives
d ϖ Q | E = ( R Q ( d G ) α Θ α + Q A ( D i d G A ) d x i ) | E ( mod total D ) .
This yields
( pr X P ( R Q ( d G ) α Θ α + Q A ( D i d G A ) d x i ) ) | E = ( R Q ( G ( P ) ) α Θ α P α R Q ( d G ) α + Q A ( D i G ( P ) A ) d x i ) | E = ( R Q * ( P ) A d G A ) | E ( mod total D ) .
Equating expressions (66) and (68) then gives the equation
( ( Q ( P ) A + R Q * ( P ) A ) d G A ) | E = 0 ( mod total D ) | E ,
showing that the symmetry action (61) produces an adjoint-symmetry. □
Observe that the three actions (56), (60) and (61) are related by:
S 1 P ( Q ) + S 2 P ( Q ) = S P ( Q ) .
Each action is a mapping on the linear space of adjoint-symmetries Q A . The algebraic properties of these actions can be found in [13].

5. Geometrical Adjoint-Symmetries of Evolution Equations

A general system of evolution equations of order N has the form
u t α = g α ( x , u , x u , , x N u )
where t is the time variable, x i , i = 1 , , n , are now the space variables, and u α , α = 1 , , m , are the dependent variables. The space of formal solutions u α ( t , x ) of the system will be denoted E .
The developments for general PDE systems can be specialized to evolution systems, with G α = u t α g α via identifying the indices A = α ( M = m ). On E , since u t α can be eliminated through the evolution equations, the components of symmetries and adjoint-symmetries can be assumed to contain only u α and its spatial derivatives in addition to t and x i . Hereafter, multi-indices will refer to spatial derivatives.
A symmetry is thereby an evolutionary vector field,
X P = P α ( t , x , x u , , x k u ) u α
satisfying the linearization of the evolution system on E :
( pr X P ( u t α g α ) ) | E = ( D t P α g ( P ) α ) | E = 0 .
Off of E , D t P α = ( P t + P ( g ) ) α + P ( G ) α , whereby R P = P . Consequently, the symmetry determining equation (73) can be expressed simply as:
( P t + [ g , P ] ) α = 0 .
The determining equation for adjoint-symmetries Q α ( t , x , x u , , x l u ) is given by the adjoint linearization of the evolution system on E :
( D t Q g * ( Q ) ) α | E = 0 .
Similar to the symmetry case, here, R Q = Q off of E , and the adjoint-symmetry determining equation simply becomes
( Q t + Q ( g ) + g * ( Q ) ) α = 0 .
These two determining equations have a geometrical formulation given by a Lie derivative defined in terms of a flow arising from the evolution system, similar to the situation for ODEs [8]. Specifically, observe that D t u α | E = g α , and hence, D t f | E = f t + f ( g ) for any differential function f. This motivates introducing the flow vector field,
Y = t + g α u α
which is related to the total time derivative by prolongation,
pr Y = D t | E = t + ( D I g α ) u I α .
Associated with this flow vector field is the Lie derivative
L t : = L pr Y
which acts on differential functions by L t f = pr Y ( f ) = D t f | E . On evolutionary vector fields (72), this Lie derivative acts in the standard way as a commutator:
L t pr X P = pr ( ( pr Y ( P ) pr X P ( g ) ) α u α ) = pr ( ( P t + P ( g ) g ( P ) ) α u α ) = pr ( ( P t + [ g , P ] ) α u α ) .
Thus, the symmetry determining equation (74) can be formulated as the vanishing of the Lie derivative expression (80). This establishes the following well-known geometrical result.
Proposition 3.
A symmetry of an evolution system (71) is an evolutionary vector field (72) that is invariant under the associated flow (79).
In particular, the resulting Lie-derivative vector field
L t X P = ( P t + [ g , P ] ) α u α
vanishes iff the functions P α are the components of a symmetry.
A similar characterization will now be given for adjoint-symmetries, based on viewing the adjoint relation between the determining equations (74) and (76) as a duality relation between vectors and one-forms.
Introduce the evolutionary one-form:
ω Q = Q α ( t , x , x u , , x l u ) d u α .
Its Lie derivative is given by
L t ω Q = ( L t Q α ) d u α + Q α L t ( d u α ) = ( Q t + Q ( g ) ) α d u α + Q α d ( L t u α ) = ( Q t + Q ( g ) ) α d u α + Q α d g α = ( Q t + Q ( g ) + g * ( Q ) ) α d u α ( mod total D ) .
This shows that the adjoint-symmetry determining equation (76) can be formulated as the functional vanishing of the Lie derivative expression (83).
Theorem 4.
An adjoint-symmetry of an evolution system (71) is an evolutionary one-form (82) that is functionally invariant under the associated flow (79).
In particular, the resulting Lie-derivative one-form
L t ω Q = ( Q t + Q ( g ) + g * ( Q ) ) α d u α ( mod total D )
functionally vanishes iff the functions Q α are the components of an adjoint-symmetry. This one-form (84) is functionally equivalent to the adjoint-symmetry one-form (28) introduced for a general PDE system. To see the relationship in detail, observe that:
ϖ Q = Q α d G α = Q α d ( u t α g α ) = Q α ( D t ( d u α ) g ( d u ) α ) = ( D t Q α + g * ( Q ) α ) d u α ( mod total D ) = L t ω Q ( mod total D ) .
An interesting question is how to extend this relationship to more general PDE systems.

Evolution Equations with Spatial Constraints

A wide generalization of evolution systems occurring in applied mathematics and mathematical physics is given by systems comprised of evolution equations with spatial constraints. Some notable examples are Maxwell’s equations, incompressible fluid equations, magnetohydrodynamical equations, and Einstein’s equations.
The constraints in such systems in general consist of spatial equations
C Υ ( x , u , x u , , x N u ) = 0 , Υ = 1 , , M
that are compatible with the evolution equation (71). Compatibility means that the time derivative of the constraints vanishes on the solution space E of the whole system, ( D t C Υ ) | E = 0 . For systems that are regular [11], Hadamard’s lemma implies that the system obeys a differential identity,
D t C Υ = C ( G ) Υ + D ( C ) Υ
where G α = u t α g α denotes the evolution equation (71), and where D is a linear differential spatial operator whose coefficients are non-singular on E . Equivalently, the constraints must obey the identity C ( g ) Υ = D ( C ) Υ . A comparison of the differential order of each side of this identity shows that D is of the same order N as the evolution equations, namely:
D = 0 | I | N R I Λ Υ D I .
The full system consists of n + M equations G α = 0 , C Υ = 0 . Note that, in the previous notation (22), ( G α , C Υ ) = ( G A ) with A = ( α , Υ ) .
The symmetry determining equation is given by the linearization of the full system on E , which is comprised by the evolution part (73) and the constraint part
( pr X P C Υ ) | E = C ( P ) Υ | E = 0 .
Off of E , C ( P ) Υ = R C ( C ) Υ , where R C is a linear differential spatial operator whose coefficients are non-singular on E . Hence, the determining equations (73) and (89) can be stated as:
( P t + [ g , P ] ) α | E C = 0 , C ( P ) Υ | E C = 0
where E C denotes the solution space of the spatial constraint equation (86).
The adjoint-symmetry determining equation is given by the adjoint linearization of the full system on E , which comprises evolution terms and additional constraint terms:
( D t Q g * ( Q ) + C * ( q ) ) α | E = 0 .
Here, the components of an adjoint-symmetry consist of
( Q α ( t , x , x u , , x l u ) , q Υ ( t , x , x u , , x l u ) )
with Q α being associated with the evolution equations as before, while q Υ is associated with the constraint equations. Similar to the symmetry case, the determining equation can be stated as:
( Q t + Q ( g ) + g * ( Q ) C * ( q ) ) α | E C = 0 .
These determining equations for symmetries and adjoint-symmetries have a geometrical formulation in terms of a constrained flow (77), generalizing the previous formulation for evolution systems as follows.
Theorem 5.
A symmetry of a constrained evolution system (71) and (86) is an evolutionary vector field (72) that is invariant under the associated constrained flow (79) and that preserves the constraints.
The proof of this result is simply the observation that, first, the determining Equation (89) corresponds to the constraints being preserved, and second, the Lie derivative of the symmetry vector field (81) along the flow vanishes on the constraint solution space.
Theorem 6.
An adjoint-symmetry of a constrained evolution system (71) and (86) is an evolutionary one-form (82) that is functionally invariant under the associated constrained flow (79), up to a functional multiple of the normal one-form d C Υ arising from the constraints.
The proof is given by the earlier computation (84) for the Lie derivative of the adjoint-symmetry one-form. This computation shows that the adjoint-symmetry determining Equation (93) now can be expressed as:
L t ω Q | E C = ( C * ( q ) α d u α ) | E C = ( q Υ d C Υ ) | E C ( mod total D )
where d C Υ is the normal one-form given by the constraints viewed as surfaces in jet space.
The Lie-derivative one-form (94) is functionally equivalent to the adjoint-symmetry one-form (28) introduced for a general PDE system. In the present notation, the full system of evolution and constraint equations (71) and (86) consists of ( G α , C Υ ) = 0 , and the corresponding one-form associated with this system is given by ϖ Q , q = Q α d G α + q Υ d C Υ . Now, using the relation (85), observe that:
ϖ Q , q = q Υ d C Υ L t ω Q ( mod total D ) .
There is a class of adjoint-symmetries arising from the summed product of arbitrary functions χ Υ ( t , x ) and the components of the the differential identity (87). This yields, after integration by parts,
0 = χ Υ ( D t C Υ C ( G ) Υ D ( C ) Υ ) = D t ( χ Υ C Υ ) + D i Ψ i ( χ , G ; C ) D i Φ i ( χ , C ; R ) ( D t χ + D * ( χ ) ) Υ C Υ C * ( χ ) α G α
where Φ i ( χ , C ; R ) = 0 | I | N 1 ( 1 ) | J | D J ( χ Υ R i I Λ Υ ) D I / J C Λ from expression (88). Hence,
D t ( χ Υ C Υ ) + D i ( Ψ i ( χ , G ; C ) Φ i ( χ , C ; R ) ) = C * ( χ ) α G α + ( D t χ + D * ( χ ) ) Υ C Υ
has the form of a conservation law off E , with ( C * ( χ ) α , ( D t χ + D * ( χ ) ) Υ ) being the multiplier. As is well known, every multiplier for a regular PDE system is an adjoint-symmetry [1,3,11,17,18]. This can be proven here by applying the Euler operator E u α and using its product rule. Consequently,
Q α = C * ( χ ) α , q Υ = ( D t χ + D * ( χ ) ) Υ
are components of an adjoint-symmetry, involving the arbitrary functions χ Υ ( t , x ) . Such adjoint-symmetries are a counterpart of gauge symmetries, and accordingly are called gauge adjoint-symmetries [11].
The corresponding gauge adjoint-symmetry one-form is given by
ω χ = C * ( χ ) α d u α = χ Υ d C Υ ( mod total D )
and satisfies the geometrical relation
L t ω χ | E C = ( ( D t χ + D * ( χ ) ) Υ d C Υ ) | E C ( mod total D ) .
This establishes the following geometrical result.
Theorem 7.
A gauge adjoint-symmetry (98) is functionally equivalent to a normal one-form ω χ associated with the constraint equation (86). Under the evolution flow, it is mapped into another normal one-form.
The preceding developments for general systems of evolution equations with spatial constraints have used the classical notion of symmetries and adjoint-symmetries. It would be interesting to extend the formulation and the results by considering a notion of conditional symmetries and corresponding conditional adjoint-symmetries based on the spatial constraints.
Specifically, on the solution space of the full system, consider a symmetry given by an evolutionary vector field (72) that satisfies
( P t + [ g , P ] ) α | E C = 0
where E C denotes the solution space of the spatial constraint Equation (86). Such conditional symmetries (101) differ from classical symmetries (90) by relaxing the condition that the constraints are preserved. Their natural adjoint counterpart is given by an evolutionary one-form (82) satisfying
( Q t + Q ( g ) + g * ( Q ) ) α | E C = 0 .
which is the adjoint of the determining Equation (101). Such conditional adjoint-symmetries (102) differ from classical adjoint-symmetries (93) by excluding the terms arising from the spatial constraints.
This notion of conditional symmetries and adjoint-symmetries is more general than the classical notion because the conditional determining equations hold on E C instead of the whole jet space.

6. Concluding Remarks

The main results showing how adjoint-symmetries correspond to evolutionary one-forms with certain geometrical properties provides a first step towards giving a fully geometrical interpretation for adjoint-symmetries. In particular, for systems of evolution equations, adjoint-symmetries can be geometrically described as one-forms that are invariant under the flow generated by the system on the solution space. This interesting result has a straightforward generalization to systems of evolution equations with spatial constraints. Consequently, the results presented here are applicable to all PDE systems of interest in applied mathematics and mathematical physics.
One direction for future work will be to translate and generalize these results into the abstract geometrical setting of secondary calculus [2,19] developed by Vinogradov and Krasil’shchik and their co-workers.
It will also be interesting to fully develop the use of adjoint-symmetries in the study of specific PDE systems, as outlined in the Introduction: finding exact solutions, detecting and finding mappings into a target class of PDEs, and detecting integrability, which are the counterparts of some important uses of symmetries. Another use of adjoint-symmetries, which has been introduced very recently [20], is for finding pre-symplectic operators.

Author Contributions

Conceptualization, S.C.A. and B.W.; methodology, S.C.A. and B.W.; writing, original draft preparation, S.C.A.; writing, review and editing, S.C.A. and B.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

S.C.A. is supported by an NSERC Discovery grant. W.B. thanks the Department of Mathematics & Statistics, Brock University, for support during a research visit when this work was completed.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bluman, G.W.; Cheviakov, A.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2009. [Google Scholar]
  2. Krasil’shchik, I.S.; Vinogradov, A.M. (Eds.) Symmetries and Conservation Laws for Differential Equations of Mathematical Physics; Translations of Math. Monographs 182; American Mathematical Society: Providence, RI, USA, 1999. [Google Scholar]
  3. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  4. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  5. Bluman, G.W.; Anco, S.C. Symmetry and Integration Methods for Differential Equations; Springer: New York, NY, USA, 2002. [Google Scholar]
  6. Sarlet, W.; Cantrijn, F.; Crampin, M. Pseudo-symmetries, Noether’s theorem and the adjoint equation. J. Phys. A Math. Gen. 1987, 20, 1365–1376. [Google Scholar] [CrossRef]
  7. Sarlet, W.; Bonne, J.V. REDUCE procedures for the study of adjoint symmetries of second-order differential equations. J. Symb. Comput. 1992, 13, 683–693. [Google Scholar] [CrossRef] [Green Version]
  8. Sarlet, W. Construction of adjoint symmetries for systems of second-order and mixed first- and second-order ordinary differential equations. Math. Comput. Model. 1997, 25, 39–49. [Google Scholar] [CrossRef]
  9. Anco, S.C.; Bluman, G. Direct construction of conservation laws from field equations. Phys. Rev. Lett. 1997, 78, 2869–2873. [Google Scholar] [CrossRef]
  10. Anco, S.C.; Bluman, G. Direct construction method for conservation laws of partial differential equations Part II: General treatment. Eur. J. Appl. Math. 2002, 41, 567–585. [Google Scholar] [CrossRef] [Green Version]
  11. Anco, S.C. Generalization of Noether’s theorem in modern form to non-variational partial differential equations. In Recent Progress and Modern Challenges in Applied Mathematics, Modeling and Computational Science; Springer: New York, NY, USA, 2017; Volume 79, pp. 119–182. [Google Scholar]
  12. Anco, S.C. On the incompleteness of Ibragimov’s conservation law theorem and its equivalence to a standard formula using symmetries and adjoint-symmetries. Symmetry 2017, 9, 33. [Google Scholar] [CrossRef] [Green Version]
  13. Anco, S.C.; Wang, B. Algebraic structures for adjoint-symmetries and symmetries of partial differential equations. arXiv 2020, arXiv:2008.07476. [Google Scholar]
  14. Nestruev, J. Smooth Manifolds and Observables; Graduate Texts in Mathematics 220; Springer: Berlin, Germany, 2002. [Google Scholar]
  15. Anco, S.C.; Pohjanpelto, J. Classification of local conservation laws of Maxwell’s equations. Acta Appl. Math. 2001, 69, 285–327. [Google Scholar] [CrossRef]
  16. Anco, S.C.; Pohjanpelto, J. Symmetries and currents of massless neutrino fields, electromagnetic and graviton fields. In CRM Proceedings and Lecture Notes (Workshop on Symmetry in Physics); American Mathematical Society: Providence, RI, USA, 2004; Volume 34, pp. 1–12. [Google Scholar]
  17. Vinogradov, A.M. The C-spectral sequence, Lagrangian formalism, and conservation laws I. The linear theory. J. Math. Anal. Appl. 1984, 100, 1–40. [Google Scholar] [CrossRef] [Green Version]
  18. Vinogradov, A.M. The C-spectral sequence, Lagrangian formalism, and conservation laws II. The nonlinear theory. J. Math. Anal. Appl. 1984, 100, 41–129. [Google Scholar] [CrossRef] [Green Version]
  19. Vinogradov, A.M. Introduction to Secondary Calculus. In Proceedings of the Conference Secondary Calculus and Cohomology Physics; Henneaux, M., Krasil’shchik, I.S., Vinogradov, A.M., Eds.; Contemporary Mathematics; American Mathematical Society: Providence, RI, USA, 1998. [Google Scholar]
  20. Anco, S.C.; Wang, B. A formula for symmetry recursion operators from non-variational symmetries of partial differential equations. arXiv 2020, arXiv:2004.03743. [Google Scholar]

Share and Cite

MDPI and ACS Style

Anco, S.C.; Wang, B. Geometrical Formulation for Adjoint-Symmetries of Partial Differential Equations. Symmetry 2020, 12, 1547. https://doi.org/10.3390/sym12091547

AMA Style

Anco SC, Wang B. Geometrical Formulation for Adjoint-Symmetries of Partial Differential Equations. Symmetry. 2020; 12(9):1547. https://doi.org/10.3390/sym12091547

Chicago/Turabian Style

Anco, Stephen C., and Bao Wang. 2020. "Geometrical Formulation for Adjoint-Symmetries of Partial Differential Equations" Symmetry 12, no. 9: 1547. https://doi.org/10.3390/sym12091547

APA Style

Anco, S. C., & Wang, B. (2020). Geometrical Formulation for Adjoint-Symmetries of Partial Differential Equations. Symmetry, 12(9), 1547. https://doi.org/10.3390/sym12091547

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop