Lie -Poisson structures over differential algebras
Zharinov V.V
Steklov Mathematical InstituteE-mail: [email protected]
Abstract
In this paper we use key elements of the Olver’s approach to Hamiltonian evolution
equations in partial derivatives and propose an algebraic construction appropriate
for Hamiltonian evolution systems with constraints.
Keywords: differential algebra, differential bicomplex, Lie-Poisson structure,
Hamiltonian mapping, Hamiltonian evolution system of PDEs, constraint.
1 Introduction.
In his book [1], Chapter VII, P.J. Olver has developed a Hamiltonian approach
to evolution systems of differential equations in partial derivatives.
In essence, the main component of his technics is the construction of a Lie-Poisson structure
over the differential algebra of the differential functions (see Subsection 4.2, below),
associated with evolution systems without constraints.
In this paper we use key elements of the Olver’s scheme and present a general algebraic
construction appropriate, in particular, for Hamiltonian evolution systems of partial differential equations with constraints.
Namely, in a pure algebraic manner we describe Lie-Poisson structures over
arbitrary differential algebras. The main result is Theorem 1,
the algebraic analog of the Theorem 7.8 from the book [1],
characterizing Hamiltonian operators over differential algebra under study.
Our abstract results can be used as follows.
For a given system of constraints (a system of partial differential equations in the space
variables) one should define an appropriate differential algebra and then search for
Hamiltonian operators compatible with the given evolution system.
Note, both these problems are rather difficult.
Moreover, the second one is solvable extremely rare,
each solvable case is a great success and an important step in the study
of the given evolution system with constraints.
We use the following general notations:
F=R,C, N={1,2,3,…}⊂Z+={0,1,2,…};
M={1,…,m}, m∈N;
I=Z+M={i=(i1,…,im)∣iμ∈Z+, μ∈M}.
2 Algebraic essentials.
2.1 The basic algebra.
Let F be an associative commutative unital algebra over the number field F.
We denote by D=D(F) the set of all differentiations of the algebra F,
i.e., all F-linear mappings X:F→F (i.e., X∈EndF(F)),
satisfying the Leibniz rule
[TABLE]
The set D has two matching structures, namely,
the structure of a Lie algebra with the commutator X,Y↦[X,Y]=X∘Y−Y∘X
as the Lie bracket,
the structure of a F-module with (f⋅X)g=f⋅(Xg),
and the matching condition [X,g⋅Y]=(Xg)⋅Y+g⋅[X,Y]
for all f,g∈F, X,Y∈D.
For any index sets A,B the F-module
[TABLE]
is defined, where the multiplication f⋅η, f∈F, η∈FBA,
is defined component-wise,
i.e., (f⋅η)βα=f⋅ηβα, α∈A, β∈B.
For any F-module FBA and any differentiation X∈D the linear mapping
X:FBA→FBA is defined component-wise,
η=(ηβα)↦Xη=(Xηβα).
In particular, the Leibniz rule takes the form:
X(f⋅η)=(Xf)⋅η+f⋅(Xη), for all X∈D, f∈F, η∈FBA.
We denote by F∘BA the F-module of all elements
η=(ηβα)∈FBA, s.t. for any upper index α∈A
only a finite number of components ηβα=0.
In particular, F∘A=FA for any index set A,
while F∘B consists of all finite elements η=(ηβ)∈FB,
and F∘B=FB iff (i.e., if and only if) the index set B is finite.
Elements ϕ=(ϕα)∈FA with the upper indices (superscripts)
may be considered as “vectors”, elements ξ=(ξα)∈F∘A with the lower indices
(subscripts) may be considered as “covectors”, and elements with multiple indices
may be considered as “tensors”.
The summation over repeated upper and lower indices is assumed
and it may be considered as “pairing” (contraction).
The index sets may be infinite,
but we shall ensure that all formally infinite sums in fact will be finite.
2.2 The de Rham complex.
The F-modules Ωq=Ωq(F), q∈Z+, of q-forms over
the algebra F with values in F are defined as follows:
[TABLE]
The direct sum Ω=⊕q∈Z+Ωq (here ⊕=⊕F)
has the natural structure of an exterior algebra.
The exterior differential d∈EndF(Ω) is defined by the Cartan formula:
[TABLE]
for all q∈Z+, ω∈Ωq, X0,…,Xq∈D
(here and hereafter, the “checked” arguments are understood to be omitted).
Note, dq=d∣Ωq:Ωq→Ωq+1, q∈Z+, and d∘d=0.
Thus, the de Rham complex {Ωq;dq} of the algebra F is defined
(see, e.g., [3]).
The cohomology spaces of this complex are H^{q}=H^{q}({\cal F})=\operatorname{Ker}d^{q}\big{/}\operatorname{Im}d^{q-1},
where q∈Z+, Kerdq={ω∈Ωq∣dω=0},
Imdq−1={ω=dχ∣χ∈Ωq−1}=dΩq−1 (Ω−1=0).
2.3 The differential algebra. The basic assumptions.
The algebra F is called differential if a subalgebra DH=DH(F)
of the Lie algebra D is specified with a F-basis D={Dμ∣μ∈M},
where commutators [Dμ,Dν]=0, μ,ν∈M.
Note, a differential algebra F is in fact the pair (F,D)
(see, e.g., [4],[5]).
Let (F,D) be a differential algebra under study.
Assumption 1**.**
We assume that the F-module D is splitted into a direct sum
D=DV⊕FDH
(here and hereafter, V stands for the vertical component
while H stands for the horizontal component),
where DV=DV(F) is a subalgebra of D,
such that [D,DV]⊂DV
(i.e. [Dμ,V]∈DV for all V∈DV and μ∈M).
Moreover, we assume that the Lie algebra DV has a formal F-basis
{∂a∣a∈A}, where [∂a,∂b]=0, a,b∈A,
A is an index set.
In this case, [Dμ,∂a]=Γμab⋅∂b,
where the symbol Γ=(Γμab)∈F∘MAA, a,b∈A, μ∈M.
Note, that here the set M is finite, while the set A is as a rule infinite,
and according to the above definition of the F-module F∘MAA,
for any index b∈A only a finite number of coefficients Γμab=0.
Assumption 2**.**
We assume that for any f∈F the action ∂af=0
only for a finite number of indices a∈A (i.e. ∂f=(∂af)∈F∘A).
In particular, DV={V=ϕa⋅∂a∣ϕa∈F}, and the action
of the formally infinite sum ϕa⋅∂a on F is well defined.
Assumption 3**.**
We assume that the differential algebra (F,D) is of the du Bois-Reymond type,
i.e., it has the property:
the equality ϕ⋅ψ=Dμχμ is valid for a fixed ϕ∈F and all
ψ∈F with some χμ=χμ(ϕ,ψ)∈F, μ∈M, iff ϕ=0.
Let us denote by
DE=DE(F)={V∈DV∣[Dμ,V]=0, μ∈M} the set of all evolutionary
differentiations of the algebra F.
The set DE is a subalgebra of the Lie algebra DV, because
[DE,DE]⊂DE due to the Jacobi identity for commutators,
but it is not a submodule of the F-module DV (see, e.g., [2]).
Let us consider the set E=Ker∇={ϕ=(ϕa)∈FA∣∇ϕ=0},
where
[TABLE]
Clear, the set E has the structure of a linear subspace
of the linear space FA.
Proposition 1**.**
The mapping ev:E→DE, ϕ=(ϕa)↦U=evϕ=ϕa⋅∂a,
is an isomorphism of linear spaces.
Moreover, the structure of the Lie algebra on DE
defines the isomorphic structure on E by the rule:
[TABLE]
where ξa=Uψa−Vϕa, a∈A.
2.4 The differential bicomplex.
The splitting D=DV⊕FDH of differentiations
leads to the splitting of the differential forms,
Ω=⊕p,q∈Z+Ωpq, where
[TABLE]
in particular, Ω00=ΩV0=ΩH0=F (here, ⊕=⊕F).
In the same way, the exterior differential d:Ω→Ω
splits in the two components: d=dV+dH, where
[TABLE]
The identity d∘d=0 implies the identities:
[TABLE]
Further, the F-module ΩV1=HomF(DV;F)=DV∗
has the dual F-basis {ρa∣a∈A}, ρa(∂b)=δba, a,b∈A,
so ΩV1={ωa⋅ρa∣ω=(ωa)∈F∘A}.
In general, for any p∈N the F-module
[TABLE]
where only a finite number of coefficients ωa1…ap=0.
In the same way, the F-module ΩH1=HomF(DH;F)=DH∗
has the dual F-basis {ϑμ∣μ∈M}, ϑμ(Dν)=δνμ,
so ΩH1={ωμ⋅ϑμ∣ωμ∈F}.
For any q∈N we have ΩHq=∧FqΩH1
(in fact, ΩHq=0 for q>m).
The equality d(ω∧χ)=(dω)∧χ+(−1)ordωω∧(dχ),
ω,χ∈Ω, implies the analogous equalities for dV and dH.
It is easy to verify that
dVL=∂aL⋅ρa, dHL=DμL⋅ϑμ,
L∈Ω0=F;
dVρa=0, dHρa=Γμba⋅ρb∧ϑμ,
a∈A, μ∈M;
dVϑμ=0, dHϑμ=0, μ∈M.
Hence, for example, let
ω=ωa⋅ρa∈ΩV1, χ=χμ⋅ϑμ∈ΩH1,
then
dVω=dVωa∧ρa+ωa⋅dVρa=∂[aωb]⋅ρa∧ρb,
where \partial_{[a}\omega_{b]}=\frac{1}{2}\big{(}\partial_{a}\omega_{b}-\partial_{b}\omega_{a}\big{)};
dHω=dHωa∧ρa+ωa⋅dHρa=(−Dμωa+Γμab⋅ωb)⋅ρa∧ϑμ;
dVχ=dVχμ∧ϑμ+χμ⋅dVϑμ=∂aχμ⋅ρa∧ϑμ ;
dHχ=dHχμ∧ϑμ+χμ⋅dHϑμ=∂[μχν]⋅ϑμ∧ϑν,
where \partial_{[\mu}\chi_{\nu]}=\frac{1}{2}\big{(}\partial_{\mu}\chi_{\nu}-\partial_{\nu}\chi_{\mu}\big{)}.
In particular, the bicomplex {Ωpq;dVpq,dHpq} is defined
with the cohomology spaces
H^{pq}_{V}=\operatorname{Ker}d^{pq}_{V}\big{/}\operatorname{Im}d^{p-1,q}_{V} and
H^{pq}_{H}=\operatorname{Ker}d^{pq}_{H}\big{/}\operatorname{Im}d^{p,q-1}_{H}, p,q∈Z+.
2.5 The basic ingredients.
We need the bottom-right corner of the above bicomplex:
{diagram}
Let us set
ϑ=ϑ1∧⋯∧ϑm,
ϑμ=(−1)μ−1ϑ1∧…ϑμˇ⋯∧ϑm,
so ϑμ∧ϑν=δνμϑ, μ,ν∈M.
With these notations we have:
Ω0,m−1={ω=ψμ⋅ϑμ∣ψμ∈F}≃FM;
Ω0m={ω=L⋅ϑ∣L∈F}≃F;
\Omega^{1,m-1}=\big{\{}\omega=\chi^{\mu}_{a}\cdot\rho^{a}\wedge\boldsymbol{\vartheta}_{\mu}\ \big{|}\ \boldsymbol{\chi}=(\chi^{\mu}_{a})\in\overset{\,\circ}{{\cal F}}{}^{\mathrm{M}}_{\mathbf{A}}\big{\}}\simeq\overset{\,\circ}{{\cal F}}{}^{\mathrm{M}}_{\mathbf{A}};
\Omega^{1m}=\big{\{}\omega=f_{a}\cdot\rho^{a}\wedge\boldsymbol{\vartheta}\ \big{|}\ \mathbf{f}=(f_{a})\in\overset{\,\circ}{{\cal F}}{}_{\mathbf{A}}\big{\}}\simeq\overset{\,\circ}{{\cal F}}{}_{\mathbf{A}}.
Further,
dH0,m−1:Ω0,m−1→Ω0m,
ω=ψμ⋅ϑμ↦dHω=Divψ⋅ϑ;
dH1,m−1:Ω1,m−1→Ω1m,
ω=χaμ⋅ρa∧ϑμ↦dHω=(∇∗χ)a⋅ρa∧ϑ;
dV0m:Ω0m→Ω1m,
ω=L⋅ϑ↦dVω=(∂L)a⋅ρa∧ϑ;
where
Div=−D∗:FM→F, ψ=(ψμ)↦Divψ=Dμψμ
(D:F→FM, L↦ζ=(ζμ),ζμ=DμL);
∇∗:F∘AM→F∘A, χ=(χaμ)↦∇∗χ,
(∇∗χ)a=−Dμχaμ+Γμab⋅χbμ;
∂:F→F∘A, L↦∂L, (∂L)a=∂aL.
Thus,
H^{0m}_{H}=\Omega^{0m}\big{/}d^{0,m-1}_{H}\Omega^{0,m-1}={\cal F}\big{/}\operatorname{Div}{\cal F}^{\mathrm{M}},
DivFM=ImDiv;
H^{1m}_{H}=\Omega^{1m}\big{/}d^{1,m-1}_{H}\Omega^{1,m-1}=\overset{\,\circ}{{\cal F}}{}_{\mathbf{A}}\big{/}\nabla^{*}\overset{\,\circ}{{\cal F}}{}^{\mathrm{M}}_{\mathbf{A}},
∇∗F∘AM=Im∇∗;
[dV0m]:HH0m→HH1m, [ω]↦[dV][ω]=[dVω],
i.e. [d^{0m}_{V}]:{\cal F}\big{/}\operatorname{Div}{\cal F}^{\mathrm{M}}\to\overset{\,\circ}{{\cal F}}{}_{\mathbf{A}}\big{/}\nabla^{*}\overset{\,\circ}{{\cal F}}{}^{\mathrm{M}}_{\mathbf{A}}, [L]↦[∂L].
We shall use the notation
\mathfrak{F}=H^{0m}_{H}={\cal F}\big{/}\operatorname{Div}{\cal F}^{\mathrm{M}}=\big{\{}\textstyle{\int\!}L=[L]=L+\operatorname{Div}{\cal F}^{\mathrm{M}}\ \big{|}\ L\in{\cal F}\big{\}};
{\cal S}=H^{1m}_{H}=\overset{\,\circ}{{\cal F}}{}_{\mathbf{A}}\big{/}\nabla^{*}\overset{\,\circ}{{\cal F}}{}^{\mathrm{M}}_{\mathbf{A}}=\big{\{}[\mathbf{f}]=\mathbf{f}+\nabla^{*}\overset{\,\circ}{{\cal F}}{}^{\mathrm{M}}_{\mathbf{A}}\ \big{|}\ \mathbf{f}\in\overset{\,\circ}{{\cal F}}{}_{\mathbf{A}}\};
δ=[dV0m]:F→S, ∫L↦δ∫L=[∂L]=∂L+∇∗F∘AM.
There are the natural pairings
⟨⋅,⋅⟩:F∘A×FA→F, f=(fa), ϕ=(ϕa)↦⟨f,ϕ⟩=fa⋅ϕa;
⟨⋅,⋅⟩:F∘AM×FMA→F, χ=(χaμ), η=(ημa)↦⟨χ,η⟩=χaμ⋅ημa;
and the linear mapping (the natural projection in the exact sequence of the F-linear spaces
0→DivFM→F→F→0)
∫:F→F,L↦∫L=L+DivFM.
It is to verify the following two statements.
Proposition 2**.**
For any ϕ∈E the image evϕDivFM⊂DivFM,
in particular the mapping evϕ:F→F is defined by the rule:
evϕ∫K=∫evϕK for any K∈F.
Proposition 3**.**
The mapping ∇∗:F∘AM→F∘A is the Lagrange dual
for the mapping ∇:FA→FMA,
i.e., for any χ∈F∘AM and ϕ∈FA
the Green formula ⟨χ,∇ϕ⟩−⟨∇∗χ,ϕ⟩=Divψ
holds, where ψ=(ψμ)∈FM, ψμ=χaμ⋅ϕa
(in other words, ∫(⟨χ,∇ϕ⟩−⟨∇∗χ,ϕ⟩)=0).
We shall need the dual F-linear space
[TABLE]
the last equality due to the du Bois-Reymond property of the differential algebra (F,D)
(remind, E=Ker∇).
Note, that for any ϕ∈FA the F-linear mapping ϕ:S→F
is defined by the rule:
[f]↦∫⟨f,ϕ⟩ for all f∈F∘A.
In particular, there is defined the pairing
[TABLE]
3 The Lie -Poisson structure.
3.1 Definition.
The Lie-Poisson structure over the differential algebra (F,D)
is a bilinear mapping (Lie-Poisson bracket)
[TABLE]
with the properties:
\big{\{}\textstyle{\int\!}K,\textstyle{\int\!}L\big{\}}+\big{\{}\textstyle{\int\!}L,\textstyle{\int\!}K\big{\}}=0;
(skew-symmetry)
\mathbf{J}\mathbf{I}\big{(}\textstyle{\int\!}K,\textstyle{\int\!}L,\textstyle{\int\!}M\big{)}=\big{\{}\textstyle{\int\!}K,\big{\{}\textstyle{\int\!}L,\textstyle{\int\!}M\big{\}}\big{\}}+\text{c.p.}=0;
(Jacobi identity)
where the abbreviation “c.p.” stands for the cyclic permutation
of arguments ∫K,∫L,∫M∈F.
In this case, the pair (F,{⋅,⋅}) is a Lie algebra.
Remark 1**.**
The Lie-Poisson bracket is not a Poisson bracket in the proper sense,
it lacks the product rule property:
[TABLE]
the product ∫K⋅∫L is not even defined in the F-linear space F.
Definition 1**.**
We define the bracket {⋅,⋅}:F×F→F as follows:
[TABLE]
where the F-linear mapping
Λ:F∘A→FA, f↦ϕ=Λf, has the properties:
Λ∘∇∗=0 and ∇∘Λ=0.
The definition is correct due to the last two equalities, because in this case
[TABLE]
Proposition 4**.**
For all K,L∈F the representation
[TABLE]
holds.
The skew-symmetry and the Jacobi identity impose additional restrictions on the mapping Λ.
The skew-symmetry reduces to the property
[TABLE]
it surely will be performed if the mapping Λ is the Lagrange skew-adjoint, i.e. if
∫(⟨f,Λg⟩+⟨g,Λf⟩)=0 for all f,g∈F∘A.
The last condition is also necessary if the image Im∂=∂F⊂F∘A
is of the du Bois-Reymond type in the proper sense.
The Jacobi identity is much more complicate.
First let us recall one simple fact from the cohomologies of Lie algebras.
Let V be a F-linear space and
[TABLE]
be the graded linear space of forms over V with values in V.
Suppose that a skew-symmetric bilinear operation
[⋅,⋅]:V×V→V is defined.
Let us define the graded linear endomorphism
[TABLE]
by the Cartan formula:
[TABLE]
where u0,…,uq∈V.
Then the endomorphism d is an exterior differential, i.e., d∘d=0, iff the Jacobi identity
[TABLE]
holds (i.e., V has the structure of a Lie algebra).
In particular, let u∈V=Ω0(V) then
[TABLE]
In our case, V=F and [⋅,⋅]={⋅,⋅}.
Thus, the following statement is valid.
Proposition 5**.**
The Jacoby identity for the skew-symmetric bilinear mapping
{⋅,⋅}:F×F→F, \textstyle{\int\!}K,\textstyle{\int\!}L\mapsto\big{\{}\textstyle{\int\!}K,\textstyle{\int\!}L\big{\}}
holds iff the corresponding linear mapping d:Ω(F)→Ω(F)
is the exterior differential.
To get any practical results concerning the Jacobi identity
we need additional technical instruments and assumptions.
3.2 Linear differential operators.
For a multi-index i=(i1,…,im)∈I we set Di=(D1)i1…(Dm)im.
We denote by
[TABLE]
the unital associative algebra of all polynomials (differential polynomials)
in the indeterminates D=(D1,…,Dm) with coefficients in F.
Every polynomial P(D)∈F[D] defines the linear mapping
(linear differential operator)
[TABLE]
*The multiplication in F[D] is defined by the composition rule
of the differential operators. *
In the same way, for index sets A,B we denote by
[TABLE]
the F[D]-module of all polynomials in the indeterminates
D=(D1,…,Dm) with coefficients in FBA.
Every polynomial P(D)∈FBA[D] defines the linear mappings:
[TABLE]
For example,
D∈FM[D],D:F→FM,K↦DK,(DK)μ=DμK,μ∈M;
Div∈FM[D],Div:FM→F,ψ=(ψμ)↦Divψ=Dμψμ;
∇∈FMAA,∇:FA→FMA,ϕ=(ϕa)↦η=(ημa),ημa=Dμϕa+Γμba⋅ϕb;
∇∗∈FMAA,∇∗:F∘AM→F∘A,χ=(χaμ)↦f=(fa),fa=−Dμχaμ+Γμab⋅χbμ.
Every differential polynomial
[TABLE]
has the Lagrange dual polynomial
[TABLE]
and the Green’s formula
[TABLE]
holds, where the “current” ψ=(ψμ)∈FM, ψμ=ψμ(P,K,L),
can be calculated explicitly using the integration by parts method.
In the general case, when the polynomial P(D)∈FBA[D] the situation is similar.
For example, let
[TABLE]
where ϕa=(P(D)g)a=Piab⋅Digb (here, a,b∈A).
Then the Lagrange dual polynomial
[TABLE]
where \xi^{a}=(P^{*}(D)\mathbf{f})^{a}=(-D)^{i}\big{(}P^{ba}_{i}\cdot f_{b}\big{)},
and the Green’s formula takes the form
[TABLE]
the “current” ψ∈FM is calculated explicitly.
In particular, ∇∗ is the Lagrange dual for ∇,
while Div is the Lagrange dual for −D.
Proposition 6**.**
For every polynomial P(D)=\big{(}P^{a}_{ib}\cdot D^{i}\big{)}:{\cal F}^{\mathbf{B}}\to{\cal F}^{\mathbf{A}},
its Lagrange dual P^{*}(D)=\big{(}(-D)^{i}\operatorname{{\scriptstyle\circ}\,}P^{a}_{ib}\big{)}:\overset{\,\circ}{{\cal F}}{}_{\mathbf{A}}\to\overset{\,\circ}{{\cal F}}{}_{\mathbf{B}}
and any ϕ∈E the following statements hold:
[\operatorname{ev}_{\boldsymbol{\phi}},P(D)]=\operatorname{ev}_{\boldsymbol{\phi}}P(D)=\big{(}(\operatorname{ev}_{\boldsymbol{\phi}}P^{a}_{ib})\operatorname{{\scriptstyle\circ}\,}D^{i}\big{)};
[\operatorname{ev}_{\boldsymbol{\phi}},P^{*}(D)]=\operatorname{ev}_{\boldsymbol{\phi}}P^{*}(D)=\big{(}(-D)^{i}\operatorname{{\scriptstyle\circ}\,}(\operatorname{ev}_{\boldsymbol{\phi}}P^{a}_{ib})\big{)}=[\operatorname{ev}_{\boldsymbol{\phi}},P(D)]^{*},
i.e., the evolutionary differentiation evϕ acts coefficient-wise.
Proof.
The proof based on the characteristic property of the evolutionary differentiations:
[Dμ,evϕ]=0 for all μ∈M and ϕ∈E.
∎
3.3 The additional assumptions.
To move further we need the following two additional assumptions.
Assumption 4**.**
We assume that there exists a linear differential operator j=j(D):FA→FA,
A is an index set, s.t.
the composition ∇∘j=0, i.e., Imj⊂Ker∇;
the commutator [evϕ,j]=0 for any ϕ∈E.
In more detail, j(D)=(jiαa⋅Di)∈FAA[D],
ϕ=(ϕα)↦jϕ=ϕ=(ϕa), ϕa=jiαa⋅Diϕα.
We define E=FA and E∗=F∘A.
Proposition 7**.**
The following statements hold:
the Lagrange dual polynomial \mathbf{j}^{*}=\big{(}(-D)^{i}\operatorname{{\scriptstyle\circ}\,}j^{a}_{i\alpha}\big{)}\in{\cal F}^{\mathbf{A}}_{\mathrm{A}}[D]
is defined, j∗:F∘A→E∗,
f=(fa)↦f=(fα), f_{\alpha}=(-D)^{i}\big{(}j^{a}_{i\alpha}\cdot f_{a}\big{)};
the composition j∗∘∇∗=0, i.e., Im∇∗⊂Kerj∗;
the commutator [evϕ,j∗]=[evϕ,j]∗=0 for any ϕ∈E;
*the linear mapping δ=j∗∘δ=j∗∘∂:F→E∗
acts by the rule: *
∫L↦δL=(δαL), \delta_{\alpha}L=(-D)^{i}\big{(}j^{a}_{i\alpha}\cdot\partial_{a}L\big{)}, α∈A;
∫⟨δL,ϕ⟩=∫L∗ϕ=∫evjϕL*, L∈F, ϕ∈E,
where *
L∗:E→F,
ϕ=(ϕα)↦L∗ϕ=∂aL⋅jiαa⋅Diϕα.
Proof.
The proof is based on the definitions, Proposition 6 and the du Bois-Reymond
property of the algebra F.
∎
Assumption 5**.**
We assume that the mapping Λ:F∘A→FA has the form
Λ=j∘Λ∘j∗, where the skew-adjoint polynomial
Λ=(Λiαβ⋅Di)∈FAA[D], while the Lagrange dual polynomial
Λ∗=((−D)i∘Λiαβ)=−Λ∈FAA[D], and
Λ,Λ∗:E∗=F∘A→E=FA.
Proposition 8**.**
Under the above assumptions Λ∈FAA[D]
is the skew-adjoint polynomial, and compositions ∇∘Λ=0, Λ∘∇∗=0.
Proposition 9**.**
For every ϕ∈E the equality [evϕ,Λ]∗=−[evϕ,Λ] holds.
Proof.
See Proposition 6.
∎
3.4 The Hamiltonian mappings.
We assume that all five assumptions listed above hold.
In this case the Lie-Poisson bracket takes the form:
[TABLE]
By the construction this bracket is skew-symmetric, and our aim here to find
a possibly simple and effective test for the skew-adjoint differential polynomial
Λ:E∗→E to be Hamiltonian, i.e., to satisfy the Jacobi identity.
For every R∈F we denote ϕ(R)=ΛδR∈E and ϕ(R)=jϕ(R)∈E.
Lemma 1**.**
For all K,L,M∈F the following equality
[TABLE]
holds.
Proof.
Indeed,
[TABLE]
∎
Theorem 1**.**
The Jacoby identity has the following representation:
[TABLE]
Proof.
Indeed,
[TABLE]
∎
Corollary 1**.**
If a differential polynomial Λ:E∗→E is skew-adjoint,
and the commutator [evϕ,Λ]=0 for all ϕ∈E,
then the polynomial Λ is Hamiltonian, i.e., the Jacobi identity holds.
4 The main example – evolution system without constraints.
4.1 The algebra.
Here (cf., [1], Chapter VII), the algebra F=Cfin∞(XU), where
X=RM={x=(xμ)∣xμ∈R, μ∈M}
is the space of independent variables;
U=RA={u=(uα)∣uα∈R, α∈A}
is the space of dependent variables;
U=RIA={u=(uiα)∣uiα∈R, α∈A, i∈I}
is the space of differential variables, uα=u0α;
Cfin∞(XU) is the algebra of all smooth functions on
the infinite dimensional space XU=X×U,
depending on a finite number of the arguments xμ,uiα.
4.2 The differential algebra.
The Lie subalgebra DH has the F-basis D={Dμ∣μ∈M},
where the total partial derivatives
[TABLE]
are characterized by the chain rule property:
\partial_{x^{\mu}}\big{(}F(x,\mathbf{u})\big{|}_{\mathbf{u}=\mathbf{j}\phi(x)}\big{)}=\big{(}D_{\mu}F(x,\mathbf{u})\big{)}\big{|}_{\mathbf{u}=\mathbf{j}\phi(x)},\quad F\in{\cal F},
ϕ=(ϕα(x))=FA, jϕ=ϕ=(ϕiα(x))∈FIA,
ϕiα(x)=Diϕα(x).
The Lie subalgebra DV has the F-basis {∂uiα∣α∈A, i∈I},
thus, here, a=(iα), A=(IA),
and the commutator
[TABLE]
4.3 The key ingredients.
There are defined the mappings:
∇:FIA→FMIA,ϕ=(ϕiα)↦η=(ημiα),ημiα=Dμϕiα−ϕi+(μ)α;
j:FA→FIA,ϕ=(ϕα)↦ϕ=(ϕiα),ϕiα=Diϕα=δikδβα⋅Dkϕβ.
One can verify by the induction that
[TABLE]
Let us set E=FA, then we get the exact sequence
[TABLE]
Moreover, the commutator [evϕ,j]=0 for all ϕ∈E,
because the differential polynomial
j=j(D)∈FIAA[D] has constant coefficients
jikβα=δikδβα.
Further, the F-module ΩH1 has the dual F-basis {dxμ∣μ∈M},
and for the F-module ΩV1 it is convenient to chose the dual F-basis
[TABLE]
In particular, dVδuiα=0, dHδuiα=−δui+(μ)α∧dxμ,
α∈A , i∈I.
There are defined the Lagrange dual mappings:
∇∗:F∘AMI→F∘AI,χ=(χαμi)↦f=(fαi),fαi=Dμχαμi+δj+(μ)iχαμj;
j∗:F∘AI→FA,f=(fαi)↦j∗f=f=(fα),fα=(−D)iδikfαk.
Let us set E∗=FA.
Proposition 10**.**
Let f=(fαi),g=(gαi)∈F∘AI, gαi=δ0ifα,
fα=(j∗f)α∈E,
then
[TABLE]
the summation is over all k,j∈I,
remind that only a finite number of components fαj=0.
Proof.
The proof is done by the direct check.
∎
Corollary 2**.**
The linear space F∘AI splits into the direct sum of the linear spaces,
F∘AI=E∗⊕F∇∗F∘AMI,
where the injection mapping E∗→F∘AI acts by the rule:
f=(fα)↦f=(fαi), fαi=δ0ifα.
Corollary 3**.**
The following sequence of the linear spaces
[TABLE]
is exact.
Corollary 4**.**
There is the representation {\cal S}=\overset{\,\circ}{{\cal F}}{}^{\mathbb{I}}_{\mathrm{A}}\big{/}\nabla^{*}\overset{\,\circ}{{\cal F}}{}^{\mathbb{I}\mathrm{M}}_{\mathrm{A}}={\cal E}^{*}.
4.4 The verification of the assumptions.
Let us check that here all five assumptions are fulfilled.
The Lie algebra D=D(F) is splitted
into the direct sum of vertical and horizontal subalgebras,
D=DV⊕FDH.
For every f∈F the derivative ∂uiαf=0 only for a finite number
of indices a=(iα)∈A=(IA).
The differential algebra (F,D), F=Cfin∞(XU),
is of the du Bois-Reymond type (see, e.g., [1]).
The mapping j:E→E here enjoys the properties:
Imj=Ker∇ and [evϕ,j]=0 for all ϕ∈E,
so Assumption 4 is fulfilled. Moreover, here Im∇∗=Kerj∗, also.
The Assumption 5, in fact, is a choice of the representation for the mapping Λ,
it doesn’t lead to any additional restrictions.
5 Conclusion.
Let us summarize.
It was shown that the Lie-Poisson structures admit a pure algebraic formulation
as Hamiltonian mappings over differential algebras.
The main example, in particular, proves that the five assumptions, listed above,
are compatible and up to date.
The crucial component of the scheme is the symbol
Γ=(Γμab)∈FMAA,
determining all further constructions.
The main (and the most difficult) problem at this stage is to find
a mapping j:FA→FA with the necessary properties
(i.e., a jet mapping).
What are possible applications of our construction?
Suppose that an evolution system of partial differential equations with constraints
is given and one wants to try to represent it as a Hamiltonian system.
To apply the above scheme in this situation one should to specify an appropriate
differential algebra compatible with constraints and satisfying to the three basic assumptions.
The natural way to do this provides the algebraic approach to partial differential equations
(see, e.g., [2]). After that one should act as it is described in the introduction.
Note, that applications of our construction are not restricted to evolution systems
of partial differential equations. In particular, this method can be applied,
e.g., in the following researches: [6], [7], [8], [9],
[10], [11], [12], [13],[14], [15], [16], etc..