arxiv:math-ph/ v1 1 Jun 2005

Transcription

1 PRE-MULTISYMPLECTIC CONSTRAINT ALGORITHM FOR FIELD THEORIES arxiv:math-ph/ v1 1 Jun 2005 Manuel de León Instituto de Matemáticas y Física Fundamental, CSIC C/ Serrano 123. E Madrid. Spain Jesús Marín-Solano, Departamento de Matemática Económica, Financiera y Actuarial, UB Av. Diagonal 690. E Barcelona. Spain Juan Carlos Marrero, Departamento de Matemática Fundamental, Fac. Matemáticas, U. La Laguna La Laguna, Tenerife. Spain Miguel C. Muñoz-Lecanda, Narciso Román-Roy, Departamento de Matemática Aplicada IV Edificio C-3, Campus Norte UPC C/ Jordi Girona 1. E Barcelona. Spain September 27, 2017 Abstract We present a geometric algorithm for obtaining consistent solutions to systems of partial differential equations, mainly arising from singular covariant first-order classical field theories. This algorithm gives an intrinsic description of all the constraint submanifolds. The field equations are stated geometrically, either representing their solutions by integrable connections or, what is equivalent, by certain kinds of integrable m-vector fields. First, we consider the problem of finding connections or multivector fields solutions to the field equations in a general framework: a pre-multisymplectic fibre bundle (which will be identified with the first-order jet bundle and the multimomentum bundle when Lagrangian and Hamiltonian field theories are considered. Then, the problem is stated and solved in a linear context, and a pointwise application of the results leads to the algorithm for the general case. In a second step, the integrability of the solutions is also studied. Finally, the method is applied to Lagrangian and Hamiltonian field theories and, for the former, the problem of finding holonomic solutions is also analized. Key words: Fibre bundles, connections, multisymplectic manifolds, constraints, field theories. AMS s. c.(2000: 35Q99, 37J05, 53C05, 53D99, 55R10, 58A20, 70S05. PACS (1999: Vh, Ef. mdeleon@imaff.cfmacc.csic.es jmarin@ub.edu jcmarrer@ull.es matmcml@ma4.upc.edu matnrr@ma4.upc.edu

2 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 2 Contents 1 Introduction 2 2 Linear theory Statement of the problem. Equivalences Maps induced by a section Characterization of solutions The general multisymplectic case Statement of the problem Conditions for the existence of solutions on a submanifold of the total space The pre-multisymplectic constraint algorithm The integrability algorithm Application to Lagrangian and Hamiltonian field theories Lagrangian and Hamiltonian field theories Lagrangian and Hamiltonian algorithms Almost-regular Lagrangians and integrable Euler-Lagrange m-vector fields An example: affine Lagrangian densities 24 1 Introduction Systems of singular differential equations have been a matter of increasing interest, especially during the last 30 years, and they have been studied separately in theoretical physics and in some technical areas such as engineering of electric networks or control theory. The fundamental characteristic of these kinds of systems is that the existence and uniqueness of solutions are not assured. In particular, this situation arises in mechanics when dynamical systems described by singular Lagrangians are considered. Furthermore, these systems do not have a nice Hamiltonian description, since not all the momenta are available, and there is a submanifold of the momentum phase space where, in general, the dynamical equations have no solution everywhere. The same problems arise when considering systems of PDE s associated with field theories described by singular Lagrangians (indeed, many field theories are singular, for instance electromagnetism, as well as in some other applications related with optimal control theories. Dirac [7] was pioneering in solving the problem for the Hamiltonian formalism of singular me-

3 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 3 chanical systems, by developing a constraint algorithm which gives, in the favourable cases, a final constraint submanifold where admissible solutions to the dynamics exist (in the sense that the dynamical evolution remains on this manifold. Dirac s main aim was to apply this procedure to field theories. After Dirac, a lot of work was done in order to geometrize his algorithm. The first important step was the work by Gotay et al [15], and its application to the Lagrangian formalism [16, 17]. Other algorithms were given later, in order to find consistent solutions of the dynamical equations in the Lagrangian formalism of singular systems (including the problem of finding holonomic solutions [1, 24, 38], and afterwards, new geometric algorithms were developed to be applied both in the Hamiltonian and the Lagrangian formalisms [18, 20, 23, 35, 37, 42]. The Lagrangian and Hamiltonian descriptions of field theories, termed the multisymplectic approach, is the natural extension of time-dependent mechanics. Therefore, in order to understand the constraint algorithm for field theories in a covariant formalism, the first step was to develop the algorithmic procedures for time-dependent systems. This work was provided in [4, 5, 14, 19, 22, 27, 28, 29, 31, 33, 36, 43]. A basic geometric study of these systems can be found in [6]. Furthermore, a qualitative description of constraint algorithms for field theories was made in [9, 11]. Working within the framework of the multisymplectic description for these theories, we present in this paper a geometric algorithm for finding the maximal submanifold where there are consistent solutions to the field equations of singular theories. This algorithm gives an intrinsic description of all the constraint submanifolds. The problem is stated in a generic pre-multisymplectic fibre bundle,inordertogiveasolutiontobothlagrangianandhamiltonianfieldtheories, aswellasother possible kinds of systems of partial differential equations. In this framework, the solutions to these equations are given geometrically by integrable connections or, what is equivalent, by integrable locally decomposable m-vector fields which are transverse to the fibre projection. The key point consists in using an auxiliar connection for constructing different geometrical structures needed to develop the algorithm, by following the same methods introduced in [28] for time-dependent singular systems. This technique (the use of a connection was used for the first time in [3], in order to obtain (global Hamiltonian functions, and afterwards applied both in the Lagrangian and Hamiltonian formalisms for this and other purposes (see [8, 9, 14, 33, 34, 40]. An exhaustive use of this technique in mechanics and field theory can be found in [31, 32, 39]. First, the problem is reduced to another in the realm of linear algebra, and solved in this context, and then the results are applied to the general pre-multisymplectic framework. In this way, a constraint algorithm can be developed giving a sequence of submanifolds which, in the best case, ends in some final constraint submanifold where field equations have consistent solutions (connections or m-vector fields, although not necessarily integrable. The problem of integrability is considered and solved separately. Finally, Lagrangian and Hamiltonian field theories are particular cases where the above results are applied straightforwardly, although in the Lagrangian case the problem of finding holonomic solutions must be also analized. The paper is organized as follows: First, in Section 2, we state and solve the algebraic version of the problem. Then, in Section 3, we pose the general problem in the context of a pre-multisymplectic fiber bundle and, applying the results obtained in the previous Section, the solution is achieved after studying the additional problem of integrability. After this, Section 4 is devoted to giving the application to Lagrangian and Hamiltonian field theories, including the problem of finding holonomic solutions in the Lagrangian formalism. Finally, as a classical example, field theories described by affine Lagrangians are analyzed in Section 5. An Appendix about multivector fields and connections is also included, in order to make the paper more self-contained and readable.

4 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 4 Manifolds are real, paracompact, connected and C. Maps are C. Sum over crossed repeated indices is understood. 2 Linear theory 2.1 Statement of the problem. Equivalences The problem we want to solve can be first posed and solved in a linear algebraic way. In fact, let W and E be R-vector spaces (although, instead of R, another field of characteristic different from 2 can be used, with dim E = m, and dim W = m+n. Let σ: W E be a surjective morphism, and denote V(σ = ker σ, and by j: V(σ W the natural injection. Consider the exact sequence 0 V(σ j W σ E 0 (1 Supposethat η Λ m E is a volume element; denote ω = σ η, andassumethat a formω Λ m+1 W and a subspace C of W are given. We denote this collection of data as (σ;η,ω;c. Next we consider the following problems in (σ;η,ω;c: Statement 1 To find a m-vector X Λ m C satisfying that: 1. X is decomposable. 2. i(xω = i(xω = 0. Statement 2 To find a subspace H C satisfying that: 1. dim H = dim E = m. 2. σ H : H E is an isomorphism. 3. [i(wω] H = 0, w W. Observe that condition 2 is equivalent to W = H V(σ. Statement 3 To find a linear map h: E C W satisfying that: 1. σ h = Id E. 2. [i(wω] Imh = 0, for every w W. Proposition 1 Statements 1, 2, and 3 are equivalent, that is, from every solution to some of these problems we can obtain a solution to the others. ( Proof (1 = 2 Let X Λ m C be a solution to the problem 1. As a consequence of the first condition, we have X = w 1... w m, with w α C. If e α = σ(w α E, for every α = 1,...,m, by the second condition we have η(e 1,...,e m = 1, and hence {e α } is a basis of E. Consider the subspace H = w 1,...,w m. We have obviously that dim H = m and that the restriction σ H : H E is an isomorphism. Furthermore, [i(wω](w 1,...,w m = 0, for every w W. Thus H is a solution to problem 2. (2 = 3 Let H be a solution to problem 2. So σ H is an isomorphism. If j H : H C is the natural injection, let h: E C be defined as h := j H (σ H 1. This map is a solution to problem 3 because the first condition holds straightforwardly and, as Im h = H, the second condition holds.

5 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 5 (3 = 1 Let h be a solution to problem 3. If {e 1,...,e m } is a basis of E satisfying that η(e 1,...,e m = 1, let w α = h(e α, and X = w 1... w m. Then X Λ m C is a solution to problem 1, because it is decomposable, and Furthermore, if w W, i(xω = ω(w 1,...,w m = η(e 1,...,e m = 1. [i(xω](w = Ω(w 1,...,w m,w = ( 1 m [i(wω](w 1,...,w m = 0 since h is a solution to problem 3, and w 1,...,w m Imh. 2.2 Maps induced by a section Consider the exact sequence (1, and let : E W be a section of σ. Denote H( := Im. We have the splitting W = H( V(σ H( is called the horizontal subspace of, and V(σ is the vertical subspace of σ. Note that σ H( is an isomorphism. The above splitting induces the natural projections σ H : W H( W ; σ V : W V(σ W with σ H + σv = Id W; and, for every w W, we write w = w H + wv, where wh H( and w V V(σ are called the horizontal and vertical components of w induced by. In the same way we have the induced splitting W = H ( V (σ whereh ( is identifiedwith theset {β W ;β σ V = 0}, andv (σ with{β W ;β σ H = 0}, in a natural way. This splitting of W induced by gives rise to a bigradation in Λ k W given by Λ k W = (Λ p H ( Λ q V (σ p,q=0,...,k; p+q=k Now, let Z Λ m E such that η(z = 1. With this condition, Z is unique and decomposable, since dim E = m. Consider Y η = Λm (Z Λ m W, which verifies the following properties: 1. Y η is decomposable, because if Z = e 1... e m, then Y η = (e 1... (e m. 2. ω(y η = 1, since ω(y η = σ η(λ m (Z = η[(λ m σ Λ m (Z] = η[λ m (σ (Z] = η(z = 1. Y η is said to be the m-vector associated to and η, and it generates Λm H(. The bigradation in Λ k W induces a splitting of Ω as follows: Ω = Ω (m,1 + Ω, Ω (m,1 being a (m + 1-form of bidegree (m,1, and Ω a (m + 1-form that includes the rest of components. Moreover, we have: Proposition 2 Ω (m,1 = ω γ η, where γ η := i(y η Ω. Then Ω = Ω +ω γ η.

6 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 6 ( Proof As Y η generates Λ m H(, it suffices to prove that Ω (m,1 and ω γ η coincide when acting on Y η v, for every v V(σ. Thus, as γ η vanishes on H(, we obtain Ω (m,1 (Y η v = Ω(Y η v = [i(y η Ω](v = γ η (v (ω γ η (Y η v = ω(y η γ η (v = γ η (v Finally, if h: E C is a linear map, induces a splitting h = h H +hv, where hh = σh h, and h V = σv h. Then, we introduce the map (endomorphism of W h V = hv σ = σv h σ: W V(σ W 2.3 Characterization of solutions In what follows, we assume that: Assumption 1 The (m+1-form Ω is of bidegree (m 1,2. Hence Ω = Ω (m,1 +Ω (m 1,2 (2 This is equivalent to demanding that i(v 1 i(v 2 i(v 3 Ω = 0, for every v 1,v 2,v 3 V(σ. Note that if U and V are real vector spaces of finite dimension then U V = {h : U V h is linear}. Thus, the auxiliar section induces the R-bilinear map Ω : E C (E H( V (σ h (h H,i(i([ h V ] Y η Ω V(σ (3 where i([ h V ] Yη is the m-vector on W defined as follows: for every β1,...,β m W, m (i([ h V ] Yη (β1,...,β m := Yη (β1,...,[ h V ]t (β α,...,β m Observe that, if Yη = w 1... w m, with w α W, then m i([ h V ] Yη = w 1... h V (w α... w m α=1 α=1 Theorem 1 The necessary and sufficient condition for a linear map h: E C to be a solution to the problem posed in Statement 3 is that Ω (h = (j H( (σ H( 1, (γ η V(σ (4 where j H( : H( W denotes the natural injection, and j H( (σ H( 1 is the horizontal lift associated with.

7 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 7 ( Proof (= Suppose that the linear map h: E C is a solution to the problem posed in Statement 3. Consider the linear map ϕ: E W defined by We have that and therefore, ϕ := h j H( (σ H( 1 : E W. σ ϕ = σ h σ j H( (σ H( 1 = Id Id = 0 h H = j H( (σ H( 1, h V = ϕ Now, suppose that e 1,...,e m E such that η(e 1,...,e m = 1, and let w α = (e α, for α = 1,...,m; thus Y η = w 1... w m. We obtain that h(e α = h H (e α +h V (e α = w α +h V (e α As h is a solution to the problem, using the splitting (2, for every v W, we have and the result follows. 0 = Ω(h(e 1,...,h(e m,v = Ω(w 1 +h V (e 1,...,w m +h V (e m,v m = Ω (m,1 (w 1,...,w m,v+ Ω (m 1,2 (w 1,...,h V (e α,...,w m,v α=1 = γ η (v+(i(i([ h V ] Y η Ω(v ( = Suppose that there exists a linear map h: E C such that (4 holds; that is, h H = σ H h = j H( (σ H( 1, (i(i([ h V ] Y η Ω V(σ = γ η V(σ = i(y η Ω V(σ. First we prove that h is a section of σ. In fact, σ h = σ (h H +hv = σ hh = σ j H( (σ H( 1 = Id E. Furthermore, let e 1,...,e m E, with η(e 1,...,e m = 1, and let w α = (e α, for α = 1,...,m. We have Y η = w 1... w m. Now h(e α = (j H( (σ H( 1 (e α +h V (e α = w α +h V (e α and we must prove that, if w W, then Ω(h(e 1,...,h(e m,w = 0. Note that, as h is a section of σ, it induces a splitting W = h(e V(σ, and hence w = wh H + wv h, where wh h Imh and wh V V(σ. Then Ω(h(e 1,...,h(e m,wh H = 0, and it suffices to prove that Ω(h(e 1,...,h(e m,v = 0, for every v V(σ. In fact, Ω(h(e 1,...,h(e m,v = Ω(w 1 +h V (e 1,...,w m +h V (e m,v m = Ω (m,1 (w 1,...,w m,v+ Ω (m 1,2 (w 1,...,h V (e α,...,w m,v = (i(y η Ω(v+ α=1 m Ω (m 1,2 (w 1,..., h V (w α,...,w m,v α=1 = γ η (v+(i(i([ h V ] Y η Ω(v = 0 Now, from Theorem 1, we deduce that:

8 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 8 Corollary 1 A linear map h: E C is a solution to the problem posed in Statement 3 if, and only if h H = j H( (σ H( 1, [i(i([ h V ]t Y η Ω] V(σ = γ η V(σ Let V(σ 0 W be the annihilator of V(σ. It is clear that the vector spaces H ( and V(σ 0 are isomorphic. The orthogonal complement of C with respect to Ω and is the subspace (C Ω of (E V(σ 0 V(σ defined by (C Ω := (Im Ω 0. (5 Then, from Theorem 1, we obtain: Theorem 2 There exists a solution to the problem posed in Statement 3 if, and only if, h (j H( (σ H( 1 γ η (Z = 0, for every (h,z (C Ω. (6 Note that, if C = W and (W Ω following Proposition: = {0}, then it is clear that (6 holds. This is the case in the Proposition 3 If the (m+1-form Ω Λ m+1 W given by Ω = Ω ω γ η (7 is 1-nondegenerate (that is, the map Ω : W Λ m W, defined by Ω (v = i(vω, for every v W, is injective, then (W Ω = {0}. ( Proof Let (h,z (W Ω. From the definitions of Ω and (W Ω (eqs. (3 and (5, we obtain that h (h H +i(zi(i([ h V ] Yη Ω = 0 (8 for every h Lin(E,W. In particular, this implies that h (h = 0, for every h E H(, and hence h = 0. Therefore, using (8, we deduce that i(zi(i([ h V ] Y η Ω = 0, for every h E W. As a consequence, from (7 and from assumption 1, it follows that i(zω = 0 and, since Ω is 1-nondegenerate, we have that Z = 0. 3 The general multisymplectic case 3.1 Statement of the problem The problem we wish to solve arises from the Lagrangian and Hamiltonian formalisms in field theories, although other kinds of systems can also be stated in this way. The general geometrical setting for these kinds of systems consists in giving a fibred manifold κ: F M (which in what follows is assumed to be a fibre bundle, where dim M = m > 1 and dim F = n + m, and M is an orientable manifold with volume form η Ω m (M. We denote

9 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 9 ω = κ η. We write (U;x µ,y j, µ = 1,...,m, j = 1,...,n, for local charts of coordinates in F adapted to the fibred structure, and such that ω = dx 1... dx m d m x. Let Ω Ω m+1 (F be a closed form, and consider the triad (F,Ω,ω. The form Ω is said to be a multisymplectic form if it is 1-nondegenerate, that is, if the map Ω : TF Λ m T F, defined by Ω (v = i(vω, for every v W, is injective. In this case, the system described by the above triad is called a multisymplectic system. Otherwise, the form is said to be a pre-multisymplectic form, and the system is pre-multisymplectic. The problem is stated as follows: Statement 4 Given apre-multisymplectic system (F,Ω,ω, we want to find asubmanifold j C : C F, and a κ-transverse, locally decomposable and integrable m-vector field X C along C, in the fibration κ: F M, such that i(x C (yω(y = 0, for every y C. (9 First we obviate the integrability condition. Hence the problem consists in finding a submanifold C F and a locally decomposable m-vector field X C X m (F along C such that i(x C (yω(y = 1, i(x C (yω(y = 0, for every y C. (10 (Note that the first equation implies that X C is κ-transverse. Taking into account Remark 5 in the Appendix and Proposition 1, we have: Proposition 4 If C is a submanifold of F, then there exists a solution to the problem stated in Statement 4 if, and only if, at every point y C, there is h y T κ(y M T yc = Lin(T κ(y M,T y C such that 1. h y is κ-transverse (that is, it is a connection along C: T y κ TyC h y = Id. (11 2. For every (X 1 κ(y,...,(x m κ(y T κ(y M, and Y y T y F, Ω(y(h y ((X 1 κ(y,...,h y ((X m κ(y,y y = 0. (12 In order to solve this problem, the use of an arbitrary connection in the fibration κ: F M is required. Thus, let be a connection in κ: F M, and Yη the corresponding locally decomposable m-vector field on F such that i(yη ω = 1. As is well-known (see Appendix and Section 2.2, the connection induces a splitting Λ k T F = (Λ p H ( Λ q V (κ p,q=0,...,k; p+q=k where H( F is the horizontal subbundle associated with the connection and V(κ F is the vertical subbundle of the fibration κ : F M. Thus, we have that Ω = Ω (m,1 +Ω, Ω (m,1 beinga(m+1-formofbidegree(m,1andω a(m+1-form. Moreover, asastraightforward consequence of Proposition 2, we have that:

10 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 10 Proposition 5 Ω (m,1 = ω γ η, where γ η := i(y η Ω. Hence Ω = Ω +ω γ η. In what follows, we assume that the following condition holds: Assumption 2 The (m+1-form Ω is of bidegree (m 1,2. By Proposition 5, this is equivalent to demanding that i(z 1 i(z 2 i(z 3 Ω = 0,for every Z 1,Z 2,Z 3 X V(κ (F. Remark 1 The above assumption is justified because this is the situation in the Lagrangian and Hamiltonian formalism of field theories (see Propositions 7 and Conditions for the existence of solutions on a submanifold of the total space Taking into account the above considerations, the necessary and sufficient condition for the existence of solutions to the problem posed in the Statement 10 arises from the results obtained in Sections 2.2 and 2.3. The key consists in working at every point of the manifolds involved in this problem. Thus, if y C, the following identifications can be made: E T κ(y M, W T y F, C T y C, V(σ V y (κ Then we may consider the R-linear map defined by Ω (y : T κ(y M T yc (T κ(y M H y( V y (κ Ω (y(h y = ((h y H,i(i([ (h y V ]t (Y η (y(ω(y V y(κ. (13 Therefore, Theorem 1 and Corollary 1 lead to the following results: Theorem 3 Let y C. Then, there exists a linear map h y T κ(y M T yc such that (11 and (12 hold if, and only if, Ω(y(h y = ((T y κ H( 1, γ η (y Vy(κ where (T y κ H( 1 : T κ(y M H y ( is the horizontal lift at y associated with the connection. (Observe that (T y κ H( 1 T κ(y M H y(. Corollary 2 If y C, and h y T κ(y M T yc, then (11 and (12 hold if, and only if, (h y H = (T y κ H( 1, [i(i([ (h y V ]t (Y η (y(ω(y] Vy(κ = (γ η (y Vy(κ. Remark 2 If y C, let V y (κ 0 T yw be the annihilator of the vertical subspace V y (κ at the point y. Then we have that [Lin(T κ(y M,H y ( ] = Tκ(y M H y( = T κ(y M V y (κ 0 = Lin(T κ(y M,V y (κ 0.

11 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 11 If h y : T κ(y M H y ( and h y: T κ(y M V y(κ 0 are linear maps, {(X 1 κ(y,...,(x m κ(y } is a basis of T κ(y M such that {α 1 κ(y,...,αm κ(y } is the dual basis of T κ(y M, and η(κ(y = (α 1 κ(y (α m κ(y, X η (κ(y = (X 1 κ(y (X m κ(y, then, taking (Y i y = (T y κ H( ((X i κ(y and β i y = αi κ(y T yκ, for all i {1,...,m}, we deduce that {(Y 1 y,...,(y m y } and {βy 1,...,βm y } are a basis of H y( and V y (κ 0, respectively. Moreover, if h y ((X i κ(y = (h y j i (Y j y, h y (αi κ(y = (h y i j βj y, for every i {1,...,m} it follows that < h y,h y >= h y (h y = (h y j i (h y i j. Now, if y C, the orthogonal complement (T y C Ω (T κ(y M V y (κ 0 V y (κ defined by with respect to Ω and is the subspace of (T y C Ω = (Im Ω (y0. (14 As in Theorem 2, from Theorem 3 we obtain Theorem 4 Let y C. Then, there exists a linear map h y T κ(y M T yc such that if, and only if, (h y H = (T y κ H( 1, [i(i([ (h y V ]t (Y η (y(ω(y] Vy(κ = (γ η (y Vy(κ h y (T yκ H( 1 γ η (y(z y = 0, for every (h y,z y (T y C Ω. (15 Note that if (T y C Ω = {0} then it is clear that (15 holds. Thus, from Proposition 3, we have: Proposition 6 If the (m+1-form Ω on F given by Ω = Ω ω γη is, the map Ω : TF Λ m T F is injective, then is 1-nondegenerate, that (T y F Ω = {0}, for every y F. 3.3 The pre-multisymplectic constraint algorithm Now we apply the above results in order to solve the problem stated in Section 3.1. The procedure is algorithmic, and gives a sequence of subsets {C i } of F. Then, we assume that: Assumption 3 Every subset C i of this sequence is a regular submanifold of F, and its natural injection is an embedding. Thus, we consider the submanifold C 1 F where a solution exists, that is, C 1 = {y F h y Lin(T κ(y M,T y F such that (h y H = (T yκ H( 1, [i(i([ (h y V ] (Y η (y(ω(y] V y(κ = (γ η (y V y(κ}.

12 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 12 Then, using the results of Section 3.2, we deduce that there is a locally decomposable section X 1 of the vector bundle Λ m T C1 F C 1 such that (i(x 1 ω C1 = 1 and (i(x 1 Ω C1 = 0. However, in general, h y (T κ(y M is not a subspace of T y C 1 and then X 1 is not tangent to C 1 or, in other words, in general, X 1 is not a connection in the fibration κ: F M along C 1. Therefore, we consider the submanifold C 2 = {y 1 C 1 h y1 Lin(T κ(y1 M,T y1 C 1 such that (h y1 H = (T y 1 κ H( 1, [i(i([ (h y1 V ] (Y η (y 1(Ω(y 1 ] Vy1 (κ = (γ η (y 1 Vy1 (κ}. Then, there is a locally decomposable section X 2 of the vector bundle Λ m T C2 C 1 C 2 such that (i(x 2 ω C2 = 1 and (i(x 2 Ω C2 = 0. However, in general, X 2 is not a connection in the fibration κ: F M along C 2. Following this process, we obtain a sequence of constraint submanifolds ji i+1 For every i 1, C i is called the ith constraint submanifold. j i 1 i C i j1 2 j 1 C 1 C 0 F. (16 This procedure is called the pe-multisymplectic constraint algorithm. We have two possibilities: There exists an integer k > 0 such that dimc k m 1. This means that the equations have no solution on a submanifold of F. There exists an integer k > 0 such that C k+1 = C k C f. In such a case, there exists a connection X f in the fibration κ: F M along C f such that i(x f (y f (Ω(y f = 0, for every y f C f. In this case, C f is called the final constraint submanifold. This is the situation which is interesting to us. Note that the existence of a connection in the fibration κ: F M along C f implies that κ(c f is an open subset of M and that κ Cf : C f κ(c f is a fibration (see Remark 5 in the Appendix. In particular, dimc f m. Next we give an intrinsic characterization of the constraints which define the constraint submanifolds C i. For this purpose, we consider the vector bundle over F, whose fiber over the point y F is W(κ, = (κ (T M H( F V (κ W y (κ, = (T κ(y M H y( V y(κ = Lin(T κ(y M,H y ( V y(κ. Thehorizontalliftassociatedwiththeconnection andthe1-formγ η induceasection((tκ H( 1, (γ η V(κ of this vector bundle given by ((Tκ H( 1, (γ η V(κ(y = ((T y κ H( 1, (γ η (y V y(κ, for every y F. Furthermore, let W Ci (κ, be the vector bundle over the submanifold C i whose fiber at the point y i C i is W yi (κ,. Moreover, we may consider the orthogonal complement (T y i C i Ω of T y i C i with respect to Ω and given by (see (14 (T y i C i Ω = {(h y i,z yi Lin(T κ(y i M,V y i (κ 0 V yi (κ h y i ((h yi H + i(z yi i(i([ (h yi V ]t (Y η (y(ω(y i = 0, for every h yi Lin(T κ(yi M,T yi C i }.

13 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 13 Note that (T y i C i Ω W y i (κ,. Furthermore, if (T C i Ω is a vector subbundle of rank r of WC i (κ, (that is, the dimension of (T y i C i Ω is r, for every y i C i then one may choose a set of r local sections {(h 1,Z 1,...,(h r,z r } of the vector bundle W (κ, F such that {(h 1,Z 1 Ci,...,(h r,z r Ci } is a local basis of the space Γ((T C i Ω of sections of the vector subbundle (T C i Ω C i. In addition, using Theorem 4, we deduce Theorem 5 Every submanifold C i (i 1 in the sequence (16 may be defined as C i = {y i 1 C i 1 ((Tκ H( 1, (γ η V(κ (y i 1,(T y i 1 C i 1 Ω = 0}. Therefore, if (T C i 1 Ω is a vector subbundle of rank r of W C i 1 (κ, and {(h 1,Z 1 (i 1,..., (h r,z r (i 1 } is a set of sections of the vector bundle W (κ, F spanning locally the space Γ((T C i 1 Ω, then C i, is defined locally, as a submanifold of C i 1, as the zero set of the functions ξ (i j C (F given by ξ (i j = ((Tκ H( 1, (γ η V(κ ((h j,z j (i 1. These functions are called ith-generation constraints. 3.4 The integrability algorithm Suppose that after applying the premultisymplectic constraint algorithm we have a final constraint submanifold C f F and a connection defined by the multivector field X f in the fibration κ: F M along C f such that (9 holds on C f, that is, i(x f (yω(y = 0, for every y C f. (17 However, X f is not, in general, a flat connection. Nevertheless, in many cases, one may find a submanifold I f of C f such that (X f If is a flat connection in the fibration κ: F M along I f and (9 holds for (X f If. Next we present an algorithm which enables us to find this submanifold (which is an adapted version of that given in [10]. This is a local algorithm, that is, we are in fact working on suitable open sets in C f. Hence, let X f m µ=1 X µ be a solution to (17. Integrability condition: The condition that X f is flat is equivalent to demanding that the distribution spanned by X 1,...,X m is involutive. Then, if c f = dimc f, let Z 1,...,Z n m X(F, such that {X 1,...,X m,z 1,...,Z cf m} is a local basis of the module of vector fields on C f. Therefore, for every pair X µ,x ν (1 µ,ν m we have [X µ,x ν ] = f ρ µνx ρ +ζ l µνz l for some functions f ρ µν,ζ l µν. Consider the system ζl µν = 0 and let I 1 = {y C f ; ζ l µν(y = 0, µ,ν,l}. We have three options: 1. I 1 = C f. Then the distribution spanned by X 1,...,X m is involutive, and (X f Cf is a flat connection in the fibration κ: F M along C f.

14 M. de León et al, Pre-multisymplectic constraint algorithm for field theories I 1 =. Then the distribution spanned by X 1,...,X m is not involutive at any point in C f, and hence the m-vector field X f is not integrable. 3. I 1 is a proper subset of C f. In this case we assume that I 1 is a closed submanifold of C f and the functions ζ l µν are the constraints locally defining I 1. The distribution spanned by X 1,...,X m is involutive on I 1 ; that is, the m-vector field X f is integrable on I 1. If X f is tangent to I 1, then (X f I1 defines a flat connection in κ: F M along I 1 and (9 holds on I 1 which implies that the problem is solved. Nevertheless, this is not the case in general, so we need the following: Tangency condition: Consider the set I 2 := {y I 1 ; X f (y Λ m T y I 1 } For I 2 we have the same problem, so we define inductively, for i > 1, I i := {y I i 1 ; X f (y Λ m T y I i 1 } and assume that we obtain a sequence... I i... I 1 C f such that I i is a non-empty (closed submanifold of F, for all i, or I i =, for some i. Observe that the locally decomposable m-vector field X f = X 1... X m is tangent to I i (with I i if, and only if, X µ is tangent to I i, for every µ. Thus, using the constraints, we have that, if {ζ α (i i } is a basis of constraints defining locally I i in I i 1, the tangency condition is 0 = X µ (ζ α (i i (for every µ,α i, that is, we have I i I i+1 := {y I i ; X µ (ζ (i α i (y = 0, µ,α i }, for every i 1. The above algorithm ends at step f in one of the following two options: 1. dimi f m 1. In such a case, we deduce that it is not possible to find a submanifold I of C f such that (X f I is a flat connection in the fibration κ: F M along I. Therefore, we must consider (if it exists another connection X f along C f such that i(x f (yω(y = 0, for every y C f, and then we must repeat the above procedure. 2. I f+1 = I f. In this case I f is a submanifold of F and we deduce that (X f If is a flat connection in the fibration κ: F M along I f such that i(x f (y(ω(y = 0, for every y I f. Thus, the problem is solved. As in Section 3.3, we remark that the existence of a connection in the fibration κ: F M along I f implies that κ(i f is an open subset of M and that κ If : I f κ(i f is a fibration. In particular, dimi f m. We will call this procedure the integrability algorithm for decomposable m-vector fields. 4 Application to Lagrangian and Hamiltonian field theories 4.1 Lagrangian and Hamiltonian field theories (For details on the construction of the Lagrangian and Hamiltonian formalisms of field theories, see for instance, [2], [3], [9], [10], [11], [12], [13], [21], [25], [26], [39], [41].

15 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 15 A first-order classical field theory is described by its configuration fibre bundle π: E M and a Lagrangian density which is a π 1 -semibasic m-form, L, on J 1 π (the first-order jet bundle of π: E M. L is usually written as L = L( π 1 η Lω, where L C (J 1 π is the Lagrangian function associated with L and ω, and π 1 : J 1 π E and π 1 := π π: J 1 π M are the natural projections. The Poincaré-Cartan m and (m + 1-forms associated with the Lagrangian density L are defined using the vertical endomorphism V of the bundle J 1 π Θ L := i(vl+l Ω m (J 1 π ; Ω L := dθ L Ω m+1 (J 1 π Then a Lagrangian system is a couple (J 1 π,ω L. The Lagrangian system is regular if Ω L is 1- nondegenerate. Elsewhere it is called singular. In a natural chart of coordinates (x α,y A,v A α in J 1 π (adapted to the bundle structure, and such that ω = dx 1... dx m dx m we have dvν B dya d m 1 x α 2 L y B vα A dy B dy A d m 1 x α ( + 2 L vν B vα A vαdv A ν B d m 2 L x+ y B vα A vα A L y B + 2 L x α vα B Ω L = 2 L v B ν v A α ( (where d m x = dx 1... dx m and d m 1 x α i ( x α 2 L equivalent to det vα v A ν B (ȳ 0, for every ȳ J 1 π. dy B d m x (18 d m x. Locally, the regularity condition is The Lagrangian problem associated with a Lagrangian system (J 1 π,ω L consists in finding sections φ Γ(M,E (where Γ(M,E denotes the set of sections of π, such that (j 1 φ i(xω L = 0, for every X X(J 1 π In natural coordinates this is equivalent to demanding that φ satisfies the Euler-Lagrange equations. The problem of finding these sections can be formulated equivalently as follows: to find the integral sections of a class of holonomic m-vector fields {X L } X m (J 1 π, such that i(x L Ω L = 0, for every X L {X L } (HolonomicmeansthatX L isintegrableanditsintegral sectionsareholonomic. Thisisequivalent to demanding that X L is integrable and semi-holonomic, that is, it satisfies the condition i(x L V = 0. Semi-holonomic (not necessarily integrable locally decomposable m-vector fields which are solution to these equations are called Euler-Lagrange m-vector fields for (J 1 π,ω L. For the Hamiltonian formalism of field theories, we take as the multimomentum bundle the manifold J 1 π Λ m 2 T E/π Λ m T M, where Λ m 2 T E Mπ is the bundle of m-forms on E vanishing by the action of two π-vertical vector fields. It is a bundle τ 1 = π τ 1 : J 1 π M, where τ 1 : J 1 π E is the natural projection. Natural charts of coordinates in Mπ and J 1 π (adapted to the bundle structure, and such that ω τ 1 η = dx 1... dx m dx m are denoted by (x α,y A,p α A,p and (xα,y A,p α A, respectively. As Mπ is a subbundle of Λ m T E (the multicotangent bundle of E of order m, then Mπ is endowed with canonical forms: the tautological form Θ Ω m (Mπ, and the multisymplectic form Ω := dθ Ω m+1 (Mπ. They are known as the multimomentum Liouville m and (m+1- forms. Their local expressions are Θ = p α Ady A d m 1 x α +pd m x, Ω = dp α A dy A d m 1 x α dp d m x Now, if (J 1 π,ω L is a Lagrangian system, the extended Legendre map associated with L, FL: J 1 π Mπ, is defined by: ( FLȳ(Z 1,...,Z m := (Θ L ȳ( Z 1,..., Z m, for ȳ J 1 π, where

16 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 16 Z 1,...,Z m T π 1 (ȳe, and Z 1,..., Z m TȳJ 1 π are such that Tȳπ 1 Zα = Z α. Then, using the natural projection µ: Mπ J 1 π, we define the restricted Legendre map associated with L as FL := µ FL. Their local expressions are FL x ν = x ν, FL y A = y A, FL p ν A = L vν A FL x ν = x ν, FL y A = y A, FL p ν A = L vν A, FL p = L v A ν L v A ν We have that FL Θ = Θ L, and FL Ω = Ω L. (J 1 π,ω L is a regular Lagrangian system if FL is a local diffeomorphism (this definition is equivalent to that given above. Elsewhere (J 1 π,ω L is a singular Lagrangian system. As a particular case, (J 1 π,ω L is a hyper-regular Lagrangian system if FL is a global diffeomorphism. A singular Lagrangian system (J 1 π,ω L is almost-regular if P := FL(J 1 π is a closed submanifold of J 1 π, FL is a submersion onto its image, and for every ȳ J 1 π, the fibres FL 1 (FL(ȳ are connected submanifolds of J 1 π. If (J 1 π,ω L is an almost-regular Lagrangian system then P is a fibre bundle over E and M (the natural projections are denoted by τ0 1: P E and τ1 0 := π τ1 0 : P M and the µ-transverse submanifold P = FL(J 1 π Mπ is diffeomorphic to P (and we denote by j 0 : P Mπ the natural imbedding. This diffeomorphism is denoted µ: P P, and it is just the restriction of the projection µ to P. Then, taking h := µ 1, we define the Hamilton-Cartan (m + 1-form Ω 0 h = ( j 0 h Ω, which verifies that FL 0 Ω0 h = Ω L (where FL 0 is the restriction map of FL onto P. Then h is called a Hamiltonian section, and (P,Ω 0 h is the Hamiltonian system associated with the almost-regular Lagrangian system (J 1 π,ω L (see [26]. If (J 1 π,ω L is a hyper-regular Lagrangian system, then P = J 1 π, and the construction is the same. In addition, FL(J 1 π is a 1-codimensional embedded submanifold of Mπ, which is transverse to the projection µ, and is diffeomorphic to J 1 π. This diffeomorphism is µ 1, when µ is restricted to FL(J 1 π, and coincides with the map h := FL FL 1, when it is restricted onto its image. h is the Hamiltonian section in this case, and the associated Hamiltonian system is denoted by (J 1 π,ω h, where Ω h = h Ω. In a local chart of natural coordinates, the Hamiltonian section is specified by a local Hamiltonian function H C (U, U J 1 π, such that h(x α,y A,p α A (x α,y A,p α A,p = H, where H(x α,y A,p α A = (FL 1 (v A α and Ω h = dp α A dya d m 1 x α +dh d m x. L v A α L = p α A (FL 1 vα A (FL 1 L The Hamiltonian problem associated with the Hamiltonian system (P,Ω 0 h (for (J1 π,ω h is analogous, consists in finding sections ψ o Γ(M,P such that ψ o i(x 0 Ω 0 h = 0, for every X 0 X(P As in the Lagrangian case, these sections are the integral sections of a class of integrable and τ 1 0 -transverse m-vector fields {X H o } X m (P satisfying that i(x Ho Ω 0 h = 0, for every X H o {X Ho }. m-vector fields satisfying these conditions (but not necessarily integrable are called Hamilton-De Donder-Weyl m-vector fields for (P,Ω 0 h.

17 M. de León et al, Pre-multisymplectic constraint algorithm for field theories Lagrangian and Hamiltonian algorithms Let (J 1 π,ω L bealagrangian system. If is an Ehresmannconnection in the fibration π 1 : J 1 π M, let Y η be the corresponding m-vector field on J 1 π. Then, we have: Proposition 7 The Poincaré-Cartan (m + 1-form may be written as Ω L = ω (γ L η +Ω L, where (γ L η = i(y η Ω L Ω 1 (J 1 π, and Ω L is a (m+1-form on J1 π of bidegree (m 1,2 with respect to the connection. ( Proof If y J 1 π and v 1,v 2,v 3 V y ( π 1 (V( π 1 being the vertical bundle of π 1 then, from (18 we have that i(v 1 v 2 v 3 Ω L (y = 0. Thus, the result follows from Propositions 2 and 5. If (J 1 π,ω L is a hyperregular Lagrangian system (the regular case is analogous FL is a global diffeomorphism. Moreover, if Ω h is the Hamilton-Cartan (m+1-form on J 1 π, then FL Ω h = Ω L. (19 Furthermore, as FL is a global diffeomorphism, the connection induces a connection in the fibration τ 1 : J 1 π M in such a way that FL Y η = Y η, (20 where Yη is the m-vector field on J 1 π associated with and the volume form η. Thus, from (19, (20 and Proposition 7, we obtain: Proposition 8 The Hamilton-Cartan (m + 1-form may be written as Ω h = Ω h +ω (γ h η, where (γ h η = i(y η Ω h, and Ω h is a (m+1-form on J 1 π of bidegree (m 1,2 with respect to the connection. Furthermore, we may prove the following result: Proposition 9 If (J 1 π,ω L is a regular Lagrangian system, then the (m+1-forms Ω L are 1-nondegenerate. and Ω h ( Proof As FL is a diffeomorphism and FL Ω h nondegenerate. The local expression of Ω h is Ω ( θ being a 1-form such that θ x α i( x β Ω h = A,α;α β h = dp α A dya d m 1 x α +θ d m x, = 0, for every α. As a consequence, dp α A dya d m 2 x αβ θ = Ω L, it suffices to prove that Ω h is 1- ( y A dy A d m 1 x β θ ( p α A dp α A dm 1 x β (21

18 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 18 ( i y A ( i p α A Ω h = α Ω h = dy A d m 1 x α +θ ( dp α A dm 1 x α +θ y A d m x, (22 ( p α A d m x. (23 Thus, if X = λ β +µ A x β y A +να A p α is a local vector field such that i(xω h = 0 then, from A (21, it follows that λ β = 0, for every β, which implies that (see (22 and (23 µ A dp α A d m 1 x α +µ A θ( y Adm x νady α A d m 1 x α +ναθ( A p α d m x = 0. A Therefore, µ A = 0 and νa α = 0, for every A and α, that is, X = 0. If the Lagrangian is regular, then from Propositions 6 and 9, we obtain that (T y J1 π Ω L = {0}, for every y J 1 π. Thus, there exist locally decomposable m-vector fields X L on J 1 π such that Moreover, we have i(x L ω = 1, i(x L Ω L = 0. Proposition 10 If (J 1 π,ω L is a regular Lagrangian system and X L is a locally decomposable m- vector field on J 1 π such that i(x L ω = 1 and i(x L Ω L = 0 then X L is an Euler-Lagrange m-vector field for L. ( Proof We must prove that X L is semi-holonomic, that is, i(x L V = 0. For this purpose, we consider local fibred coordinates (x α,y A,vα A on J1 π. Then, since i(x L ω = 1, it follows that ( X L = Λ m α=1 x α +ΓA α y A +ΓA αβ with Γ A α and Γ A αβ local real functions on J1 π. Furthermore, from (18, we deduce that 0 = (i(x L Ω L ( v B ν v A β = ( 1 m i(x L (i( v B ν Ω L = ( 1 m+1 (Γ A α vα A 2 L vα A, for all B and ν. vb ν Therefore, using the fact that L is regular, we conclude that which implies that X L is semi-holonomic. Γ A α = v A α, for all A and α, Hence, if (J 1 π,ω L is a regular Lagrangian system, then the existence of classes of Euler- Lagrange m-vector fields for L is assured in J 1 π. In the same way, for the Hamiltonian formalism, the existence of Hamilton-De Donder-Weyl m-vector fields is assured everywhere in J 1 π (note that if X L is an Euler-Lagrange m-vector field for L then (FL X L is a Hamilton-De Donder-Weyl m-vector field on J 1 π. In both cases, the solution is not unique. For singular (almost-regular Lagrangian systems, the existence of Euler-Lagrange m-vector fields is not assured except perhaps on some submanifold S f J 1 π, where the solution is not

19 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 19 unique. In order to find this submanifold we apply the algorithm developed in Section 3.3 to the system (J 1 π,ω L, by doing the identifications κ: F M with π 1 : J 1 π M, and Ω with Ω L. Thus we obtain obtain a sequence ji i+1 j i 1 i N i j1 2 j 1 N 1 N 0 J 1 π. (24 which, in the best of cases stabilizes in the final constraint submanifold N f where there exist m-vector fields X N f on N f, solution to the equations (i(x N f Ω L Nf = 0, (i(x N f ω Nf = 1. (25 But X N f will not be, in general, an Euler-Lagrange m-vector field on N f (that is, it is not semiholonomic, and, in addition, X N f will not in general be an integrable m-vector field. The problem of finding integrable Euler-Lagrange m-vector fields (i.e., holonomic is discussed and solved in the next Section. Now, we consider the Hamiltonian system (P,Ω 0 h. Let 0 be a connection in the bundle τ 0 1: P M and denote by Y 0 η the corresponding m-vector field on P associated with 0 and η. Then, we have: Proposition 11 The Hamilton-Cartan (m + 1-form may be written as Ω 0 h = Ω 0 h +ω 0 (γ h 0 η, where ω0 = τ1 0 η, (γ h 0 η = i(y 0 η Ω h, and Ω 0 h is a (m+1-form on P of bidegree (m 1,2 with respect to the connection 0. ( Proof If ȳ = FL 0 (y P, with y J 1 π, and v 1, v 2, v 3 V y ( τ 1 0 then, since FL 0 : J 1 π P is a submersion and τ 1 0 FL 0 = π 1, it follows that there exist v 1,v 2,v 3 V y ( π 1 such that (T y FL 0 (v i = v i, for i {1,2,3}. Thus, using that (FL 0 Ω 0 h = Ω L, we deduce that This proves the result. i( v 1 v 2 v 3 Ω 0 h (ȳ = 0. Hamilton-De Donder-Weyl m-vector fields do not exist, in general, in P, and then we must apply the algorithmic procedure developed in Section 3.3 to the system (P,Ω 0 h, by doing the identifications κ: F M with τ 1 P : P M, and Ω with Ω 0 h. Thus we obtain a sequence ji i+1 j i 1 i P i j1 2 j 1 P 1 P 0 P. (26 which, in the best of cases stabilizes in the final constraint submanifold P f of P where there exist m-vector fields X P f on P f, solution to the equations Of course the solution X P f is not unique. (i(x P f Ω 0 h P f = 0, (i(x P f ω 0 P f = 1. (27

20 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 20 Remark 3 The Lagrangian and Hamiltonian pre-multisymplectic algorithms are equivalent in the following sense: at every level j of the Lagrangian and Hamiltonian algorithms, the submanifolds of the sequences (24 and (26 are FL-related, that is, FL(N j = P j and FL j = FL Nj : N j P j is a submersion such that FL 1 j (FL j (x j = FL 1 0 (FL 0(x j, for x j N j. Moreover, if N f is the final constraint submanifold (in the Lagrangian level and X N f is a locally decomposable m- vector field on N f such that equations (25 hold and, in addition, X N f is FL f -projectable to an m-vector field X P f on P f then X P f is locally decomposable and equations (27 hold. Conversely, if X P f is a locally decomposable m-vector field on P f satisfying equations (27 and X N f is a locally decomposable m-vector field on N f which is FL f -projectable on X P f then X N f satisfies equations (25 (see [26, 30] for a detailed discussion on this topic. Finally, the Hamilton-De Donder-Weyl m-vector fields X P f are not integrable, in general. In fact, if we have that X P f = X P f 1... X P f m, where X P f α are (local vector fields on P f, for all α, and {X P f 1,...,XP f m, Z 1,..., Z p } is a local basis of the vector bundle TP f P f then [X P f α,x P f β ] = f γ αβ XP f γ + ζ l αβ Z l for some functions f γ αβ and ζ l αβ on P f. Therefore, we must apply the integrability algorithm of Section 3.4, and we obtain a sequence... J i... J 1 P f, such that J i is a non-empty (closed submanifold of S f, with J 1 = {y P f ζ αβ l (y = 0} J i = {y J i 1 X P f (y Λ m T y J i 1 }, for i 2. In the best cases, there exists an integer i such that J i+1 = J i. Then, J f = J i+1 = J i is a submanifold of P f, and X J f = (X P f Jf is an integrable Hamilton-De Donder-Weyl m-vector field in J f. 4.3 Almost-regular Lagrangians and integrable Euler-Lagrange m-vector fields Let (J 1 π,ω L be an almost-regular Lagrangian system, and N f the final constraint submanifold (in the Lagrangian setting. Then, there exists a locally decomposable m-vector field X N f on N f such that (i(x N f Ω L Nf = 0, (i(x N f ω Nf = 1. But, in general, X N f is not an Euler-Lagrange m-vector field on N f and, in addition, X N f will not in general be an integrable m-vector field. In order to solve these problems, first we construct a submanifold S f of N f where there exists a locally decomposable m-vector field X S f such that (i(x S f Ω L Sf = 0, (i(x S f ω Sf = 1, (i(x S f V Sf = 0. In fact, from the above discussion we know that we can choose the m-vector field X N f on N f such that it projects via FL f (the restriction of FL to N f onto an m-vector field X P f on P f. Then, we consider the subset S f of N f defined by S f = {x N f /(i(x N f V(x = 0}. (28

21 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 21 In [30] (see also [26], it was proved that (i(x N f V(x KerT x (FL f = KerT x (FL 0 (29 and that for every x N f, S f FL 1 f (FL f(x = S f FL 1 0 (FL 0(x is a single point in S f. The above result allows us to introduce a well-defined map s f : P f N f such that S f = s f (P f, FL f s f = Id. Thus, s f : P f N f is a global section of the submersion FL f : N f P f and, therefore, S f is an embedded submanifold of N f and the map s f : P f S f is a diffeomorphism (for more details, see [26, 30]. Now, defining the m-vector field X S f on S f by X S f = (Λ m Ts f X P f, then we have [30]: Theorem 6 X S f is an Euler-Lagrange m-vector field on S f for the Lagrangian L, that is, X S f is a locally decomposable m-vector field on S f and (i(x S f Ω L Sf = 0, (i(x S f ω Sf = 1, (i(x S f V Sf = 0. Next, we give a local description of the submanifold S f and of the Euler-Lagrange m-vector field X S f ( on S f. Since L is almost-regular, ( it follows that the rank of the partial Hessian matrix 2 L 2 L vα A is constant. Let rank vb β vα A = pm + q, with 0 p n 1 and 0 q m, vb β and assume that the first pm + q rows of this matrix are independent. Denote by Ṽ(π1 J 1 π the vector subbundle of the vertical bundle V(π 1 J 1 π of π 1 : J 1 π M generated by the local vector fields { } vα A,,..., vq p+1, for 1 A p and 1 α m. v p+1 1 Then, there exist sections {X p+1 q+1,...,xp+1 m,xα A }, with p+2 A n and 1 α m, of the vector bundle Ṽ(π1 J 1 π such that {W p+1 q+1,...,wp+1 m,wα A }, with p+2 A n and 1 α m, is a local basis of Ker(T(FL 0, where W p+1 β = v p+1 +X p+1 β, Wα A = v A β α Now, suppose that X N f = X N f 1... X N f m, with ( X N f α = x α +ΓA α +Xα A, (p+2 A n, 1 α m, q+1 β m. (30 y A +ΓA αβ v A β N f, for α {1,...,m}. Then, using that X N f is FL f -projectable, it follows that the functions Γ A α are constant on the fibers of FL f : N f P f. But, as FL 1 f (FL f(x = FL 1 0 (FL 0(x, for every x N f (see Remark 3, we obtain that W p+1 W p+1+i γ (Γ A α = 0, γ {q +1,...,m} γ (Γ A α = 0, i {1,...,n p 1}, γ {1,...,m}. (31 Furthermore i(x N f V = (Γ A α vα A vα A.

22 M. de León et al, Pre-multisymplectic constraint algorithm for field theories. 22 Thus, from (29 and (30, we have that Note that the functions i(x N f V = (Γ p+1 γ vγ p+1 Wγ p+1 +(Γ p+1+i γ vγ p+1+i Wγ p+1+i. (32 ζγ p+1 = Γ p+1 γ vγ p+1, ζ p+1+i γ = Γ p+1+i γ γ {q +1,...,m} v γ p+1+i, i {1,...,n p 1}, γ {1,...,m}. (33 are independent on N f. In fact (see (30, (31 and (33 Wγ p+1 (ζ p+1 γ = δ γγ, Wγ p+1 (ζ γ p+1+i = 0, W γ p+1+i (ζγ p+1 = 0, Wᾱ p+1+i (ζ γ p+1+j = δ ij δᾱ γ. Moreover, using(28 and(32, weconcludethat{ζγ p+1,ζ γ p+1+i }, withγ {1,...,m}, i {1,...,n p 1} and γ {1,...,m}, is a set of local independent constraint functions defining S f as a submanifold of N f, that is, S f = {x N f /(Γ p+1 γ vγ p+1 (x = 0, (Γ p+1+i γ v γ p+1+i (x = 0}. Finally, a direct calculation proves that the Euler-Lagrange m-vector field X S f on S f is given by X S f = X S f 1... X S f m, with X S f α = (X N f α +X N f α (ζγ p+1 Wγ p+1 +X N f α (ζ γ p+1+i W p+1+i γ Sf, for every α. X S f is not, in general, integrable. In fact, if {X S f 1,...,XS f m,z 1,...,Z s } is a local basis of the vector bundle TS f S f then we have that for some functions f γ αβ and ζl αβ. [X S f α,x S f β ] = fγ αβ XS f γ +ζ l αβ Z l, Therefore, we must apply the integrability algorithm of Section 3.4. Then, we obtain a sequence... I i... I 1 S f, such that I i is a non-empty (closed submanifold of S f, with I 1 = {x S f ζ γ αβ (x = 0} I i = {x I i 1 X S f (x Λ m T x I i 1 }, for i 2 In the best cases, there exists an integer i such that I i+1 = I i. Then, I f = I i+1 = I i is a submanifold of S f and X I f = (X S f If is an integrable Euler-Lagrange m-vector field on I f, and hence it is holonomic. In fact: Theorem 7 If U is an open subset of M and s : U M I f is an integral section of X I f then there exists a section φ : U M E of the projection π : E M such that s = j 1 φ and φ is a solution to the Euler-Lagrange equations for L. ( Proof We have that (i(x f Ω L If = 0, (34 (i(x f V If = 0. (35