Moving Frames for Lie Pseudo Groups

Transcription

1 Moving Frames for Lie Pseudo Groups Peter J. Olver University of Minnesota olver Benasque, September, 2009

2 Sur la théorie, si importante sans doute, mais pour nous si obscure, des groupes de Lie infinis, nous ne savons rien que ce qui trouve dans les mémoires de Cartan, première exploration à travers une jungle presque impénétrable; mais celle-ci menace de se refermer sur les sentiers déjà tracés, si l on ne procède bientôt à un indispensable travail de défrichement. André Weil, 1947

3 What s the Deal with Infinite Dimensional Groups? Lie invented Lie groups to study symmetry and solution of differential equations. In Lie s time, there were no abstract Lie groups. All groups were realized by their action on a space. Therefore, Lie saw no essential distinction between finitedimensional and infinite-dimensional group actions. However, with the advent of abstract Lie groups, the two subjects have gone in radically different directions. The general theory of finite-dimensional Lie groups has been rigorously formalized and applied. But there is still no generally accepted abstract object that represents an infinite-dimensional Lie pseudo-group!

4 Ehresmann s Trinity 1953: Lie Pseudo-groups Jets Groupoids

5 Ehresmann s Trinity 1953: Lie Pseudo-groups Jets Groupoids

6 Ehresmann s Trinity 1953: Lie Pseudo-groups Jets Groupoids

7 Ehresmann s Trinity 1953: Lie Pseudo-groups Jets Groupoids

8 Lie Pseudo-groups in Action Lie Medolaghi Vessiot Cartan Ehresmann Kuranishi, Spencer, Goldschmidt, Guillemin, Sternberg, Kumpera,... Relativity Noether s (Second) Theorem

9 Gauge theory and field theories: Maxwell, Yang Mills, conformal, string,... Fluid mechanics, metereology: Euler, Navier Stokes, boundary layer, quasi-geostropic,... Solitons (in dimensions): K P, Davey-Stewartson,... Kac Moody Geometric numerical integration Morphology and shape recognition Linear and linearizable PDEs Lie groups!

10 Moving Frames In collaboration with Juha Pohjanpelto and Jeongoo Cheh, I have recently established a moving frame theory for infinite-dimensional Lie pseudo-groups mimicking the earlier equivariant approach for finite-dimensional Lie groups developed with Mark Fels and others. The finite-dimensional theory and algorithms have had a very wide range of significant applications, including differential geometry, differential equations, calculus of variations, computer vision, Poisson geometry and solitons, numerical methods, relativity, classical invariant theory,...

11 What s New? In the infinite-dimensional case, the moving frame approach provides new constructive algorithms for: Invariant Maurer Cartan forms Structure equations Moving frames Differential invariants Invariant differential operators Basis Theorem Syzygies and recurrence formulae

12 Further applications: = Symmetry groups of differential equations = Vessiot group splitting; explicit solutions = Gauge theories = Calculus of variations = Invariant geometric flows

13 Symmetry Groups Review System of differential equations: ν (x, u (n) ) = 0, ν = 1, 2,..., k By a symmetry, we mean a transformation that maps solutions to solutions. Lie: To find the symmetry group of the differential equations, work infinitesimally. The vector field v = p i=1 ξ i (x, u) x i + q α=1 ϕ α (x, u) u α is an infinitesimal symmetry if its flow exp(t v) is a oneparameter symmetry group of the differential equation.

14 To find the infinitesimal symmetry conditions, we prolong v to the jet space whose coordinates are the derivatives appearing in the differential equation: where v (n) = ϕ J α = D J p i=1 ( ξ i x i + ϕ α p i=1 Infinitesimal invariance criterion: q α=1 u α i ξi ) + n ϕ J α #J=0 u α J p u α J,i ξi i=1 v (n) ( ν ) = 0 whenever = 0. Infinitesimal determining equations: L(x, u; ξ (n), ϕ (n) ) = 0

15 The Heat Equation Symmetry generator: Prolongation: v = τ(t, x, u) t v (2) = v + ϕ t u t = u xx + ξ(t, x, u) x u t + ϕ x u x + ϕ xx ϕ t = ϕ t + u t ϕ u u t τ t u 2 t τ u u x ξ t u t u x ξ u + ϕ(t, x, u) u ϕ x = ϕ x + u x ϕ u u t τ x u t u x τ u u x ξ x u 2 x ξ u u xx + ϕ xx = ϕ xx + u x (2ϕ xu ξ xx ) u t τ xx + u 2 x (ϕ uu 2ξ xu ) 2u x u t τ xu u 3 x ξ uu u2 x u t τ uu + u xx ϕ u u x u xx ξ u u t u xx τ u

16 Infinitesimal invariance: v (3) (u t u xx ) = ϕ t ϕ xx = 0 whenever u t = u xx Determining equations: Coefficient Monomial 0 = 2τ u u x u xt 0 = 2τ x u xt 0 = τ uu u 2 x u xx ξ u = 2τ xu 3ξ u ϕ u τ t = τ xx + ϕ u 2ξ x u x u xx u xx 0 = ξ uu u 3 x 0 = ϕ uu 2ξ xu u 2 x ξ t = 2ϕ xu ξ xx u x ϕ t = ϕ xx 1

17 General solution: ξ = c 1 + c 4 x + 2c 5 t + 4c 6 xt, τ = c 2 + 2c 4 t + 4c 6 t 2, ϕ = (c 3 c 5 x 2c 6 t c 6 x 2 )u + α(x, t), where α t = α xx is an arbitrary solution to the heat equation. Basis for the (infinite-dimensional) symmetry algebra: v 1 = x, v 2 = t, v 3 = u u, v 4 = x x + 2t t, v 5 = 2t x xu u, v 6 = 4xt x + 4t 2 t (x 2 + 2t)u u, v α = α(x, t) u, where α t = α xx. x and t translations, scalings: λ u, and (λ x, λ 2 t), Galilean boosts, inversions, and the addition of solutions stemming from the linearity of the equation.

18 Symmetry generator: Prolongation: where The Korteweg devries equation v = τ(t, x, u) t v (3) = v + ϕ t u t + u xxx + uu x = 0 + ξ(t, x, u) x u t + ϕ x + ϕ(t, x, u) u u x + + ϕ xxx u xxx ϕ t = ϕ t + u t ϕ u u t τ t u 2 t τ u u x ξ t u t u x ξ u ϕ x = ϕ x + u x ϕ u u t τ x u t u x τ u u x ξ x u 2 x ξ u ϕ xxx = ϕ xxx + 3u x ϕ u +

19 Infinitesimal invariance: v (3) (u t + u xxx + uu x ) = ϕ t + ϕ xxx + u ϕ x + u x ϕ = 0 on solutions Infinitesimal determining equations: τ x = τ u = ξ u = ϕ t = ϕ x = 0 ϕ = ξ t 2 3 uτ t ϕ u = 2 3 τ t = 2 ξ x τ tt = τ tx = τ xx = = ϕ uu = 0 General solution: τ = c 1 + 3c 4 t, ξ = c 2 + c 3 t + c 4 x, ϕ = c 3 2c 4 u.

20 Basis for symmetry algebra g KdV : v 1 = t, v 2 = x, v 3 = t x + u, v 4 = 3t t + x x 2u u. The symmetry group G KdV is four-dimensional (x, t, u) (λ 3 t + a, λ x + c t + b, λ 2 u + c)

21 v 1 = t, v 2 = x, v 3 = t x + u, v 4 = 3t t + x x 2u u. Commutator table: v 1 v 2 v 3 v 4 v v 1 v v 1 3 v 2 v 3 0 v v 3 v 4 v 1 3 v 2 2 v 3 0 Entries: [v i, v j ] = k C k ij v k. Ck ij structure constants of g

22 Navier Stokes Equations u t + u u = p + ν u, u = 0. Symmetry generators: v α = α(t) x + α (t) u α (t) x p v 0 = t s = x x + 2t t u u 2 p p r = x x + u u w h = h(t) p

23 Kadomtsev Petviashvili (KP) Equation Symmetry generators: ( u t u u x u xxx ) x ± 3 4 u yy = 0 v f = f(t) t y f (t) y + ( 1 3 x f (t) 2 9 y2 f (t) ) x + ( 2 3 u f (t) x f (t) 4 27 y2 f (t) ) u, w g = g(t) y 2 3 y g (t) x 4 9 y g (t) u, z h = h(t) x h (t) u. = Kac Moody loop algebra A (1) 4

24 Main Goals Given a system of partial differential equations: Find the structure of its symmetry (pseudo-) group G directly from the determining equations. Find and classify its differential invariants. Use symmetry reduction or group splitting to construct explicit solutions both invariant and non-invariant.

25 Main Goals Given a system of partial differential equations: Find the structure of its symmetry (pseudo-) group G directly from the determining equations. Find and classify its differential invariants. Use symmetry reduction or group splitting to construct explicit solutions.

26 Main Goals Given a system of partial differential equations: Find the structure of its symmetry (pseudo-) group G directly from the determining equations. Find and classify its differential invariants. Use symmetry reduction or group splitting to construct explicit solutions.

27 Main Goals Given a system of partial differential equations: Find the structure of its symmetry (pseudo-) group G directly from the determining equations. Find and classify its differential invariants. Use symmetry reduction or group splitting to construct explicit solutions.

28 Pseudo-groups M smooth (analytic) manifold Definition. A pseudo-group is a collection of local diffeomorphisms ϕ : M M such that Identity: 1 M G, Inverses: ϕ 1 G, Restriction: U dom ϕ = ϕ U G, Continuation: dom ϕ = [ U κ and ϕ U κ G = ϕ G, Composition: im ϕ dom ψ = ψ ϕ G.

29 Lie Pseudo-groups Definition. A Lie pseudo-group G is a pseudo-group whose transformations are the solutions to an involutive system of partial differential equations: F (z, ϕ (n) ) = 0. called the nonlinear determining equations. = analytic (Cartan Kähler) Key complication: abstract object G

30 A Non-Lie Pseudo-group Acting on M = R 2 : X = ϕ(x) Y = ϕ(y) where ϕ D(R) is any local diffeomorphism. Cannot be characterized by a system of partial differential equations (x, y, X (n), Y (n) ) = 0

31 Theorem. Any regular non-lie pseudo-group can be completed to a Lie pseudo-group with the same differential invariants. Completion of previous example: X = ϕ(x), Y = ψ(y) where ϕ, ψ D(R).

32 Infinitesimal Generators g Lie algebra of infinitesimal generators of the pseudo-group G z = (x, u) local coordinates on M Vector field: v = m a = 1 ζ a (z) z a = p i=1 ξ i x i + q ϕ α α=1 u α Vector field jet: j n v ζ (n) = (... ζ b A... ) ζ b A = #A ζ b z A = k ζ b z a 1 z a k

33 The infinitesimal generators of G are the solutions to the Infinitesimal (Linearized) Determining Equations L(z, ζ (n) ) = 0 ( ) Remark: If G is the symmetry group of a system of differential equations (x, u (n) ) = 0, then ( ) is the (involutive completion of) the usual Lie determining equations for the symmetry group.

34 The Diffeomorphism Pseudo-group M smooth m-dimensional manifold D = D(M) pseudo-group of all local diffeomorphisms Z = ϕ(z) z = (z 1,..., z m ) source coordinates Z = (Z 1,..., Z m ) target coordinates

35 Jets For 0 n : Given a smooth map ϕ : M M, written in local coordinates as Z = ϕ(z), let j n ϕ z denote its n jet at z M, i.e., its n th order Taylor polynomial or series based at z. J n (M, M) is the n th order jet bundle, whose points are the jets. Local coordinates on J n (M, M): (z, Z (n) ) = (... z a... ZA b... ), Zb A = k Z b z a 1 z a k

36 Diffeomorphism Jets The n th order diffeomorphism jet bundle is the subbundle D (n) = D (n) (M) J n (M, M) consisting of n th order jets of local diffeomorphisms ϕ : M M. The Inverse Function Theorem tells us it is defined by the non-vanishing of the Jacobian determinant: det( Z a b ) = det( Za / z b ) 0

37 Pseudo-group Jets A Lie pseudo-group G D defines the subbundle G (n) = {F (z, Z (n) ) = 0} D (n) consisting of the jets of pseudo-group diffeomorphisms, and therefore characterized by the pseudo-group s nonlinear determining equations. Local coordinates on G (n), e.g., the restricted diffeomorphism jet coordinates z c, ZB a, are viewed as the pseudo-group parameters, playing the same role as the local coordinates on a Lie group G. The pseudo-group jet bundle G (n) does not form a group, but rather a groupoid under composition of Taylor polynomials/series.

38 Double fibration: Groupoid Structure σ (n) M G (n) τ (n) M σ (n) (z, Z (n) ) = z source map τ (n) (z, Z (n) ) = Z target map You are only allowed to multiply h (n) g (n) if σ (n) (h (n) ) = τ (n) (g (n) ) Composition of Taylor polynomials/series is well-defined only when the source of the second matches the target of the first.

39 One-dimensional case: M = R Source coordinate: x Target coordinate: X Local coordinates on D (n) (R) g (n) = (x, X, X x, X xx, X xxx,..., X n ) Jet: X [ h ] = X + X x h X xx h X xxx h3 + = Taylor polynomial/series at a source point x

40 Groupoid multiplication of diffeomorphism jets: (X, X, X X, X XX,... ) (x, X, X x, X xx,... ) = (x, X, X X X x, X X X xx + X XX X 2 x,... ) = Composition of Taylor polynomials/series The higher order terms are expressed in terms of Bell polynomials according to the general Fàa di Bruno formula. The groupoid multiplication (or Taylor composition) is only defined when the source coordinate X of the first multiplicand matches the target coordinate X of the second.

41 Structure of Lie Pseudo-groups The structure of a finite-dimensional Lie group G is specified by its Maurer Cartan forms a basis µ 1,..., µ r for the right-invariant one-forms: dµ k = i<j C k ij µi µ j

42 What should be the Maurer Cartan forms of a Lie pseudo-group? Cartan: Use exterior differential systems and prolongation to determine the structure equations. I propose a direct approach based on the following observation: The Maurer Cartan forms for a pseudo-group can be identified with the right-invariant one-forms on the jet groupoid G ( ). The structure equations can be determined immediately from the infinitesimal determining equations.

43 What should be the Maurer Cartan forms of a Lie pseudo-group? Cartan: Use exterior differential systems and prolongation to determine the structure equations. I propose a direct approach based on the following observation: The Maurer Cartan forms for a pseudo-group can be identified with the right-invariant one-forms on the jet groupoid G ( ). The structure equations can be determined immediately from the infinitesimal determining equations.

44 What should be the Maurer Cartan forms of a Lie pseudo-group? Cartan: Use exterior differential systems and prolongation to determine the structure equations. I propose a direct approach based on the following observation: The Maurer Cartan forms for a pseudo-group can be identified with the right-invariant one-forms on the jet groupoid G ( ). The structure equations can be determined immediately from the infinitesimal determining equations.

45 What should be the Maurer Cartan forms of a Lie pseudo-group? Cartan: Use exterior differential systems and prolongation to determine the structure equations. I propose a direct approach based on the following observation: The Maurer Cartan forms for a pseudo-group can be identified with the right-invariant one-forms on the jet groupoid G ( ). The structure equations can be determined immediately from the infinitesimal determining equations.

46 What should be the Maurer Cartan forms of a Lie pseudo-group? Cartan: Use exterior differential systems and prolongation to determine the structure equations. I propose a direct approach based on the following observation: The Maurer Cartan forms for a pseudo-group can be identified with the right-invariant one-forms on the jet groupoid G ( ). The structure equations can be determined immediately from the infinitesimal determining equations.

47 The Variational Bicomplex The differential one-forms on an infinite jet bundle split into two types: horizontal forms contact forms Consequently, the exterior derivative d = d M + d G on D ( ) splits into horizontal (manifold) and contact (group) components, leading to the variational bicomplex structure on the algebra of differential forms on D ( ).

48 For the diffeomorphism jet bundle Local coordinates: D ( ) J (M, M) z 1,... z m, Z 1,... Z m,... ZA b,... }{{}}{{} source target jet }{{} Horizontal forms: dz 1,..., dz m Basis contact forms: Θ b A = d G Zb A = dzb A m a = 1 Z a A,a dza

49 One-dimensional case: M = R Local coordinates on D ( ) (R) (x, X, X x, X xx, X xxx,..., X n,... ) Horizontal form: Contact forms: dx Θ = dx X x dx Θ x = dx x X xx dx Θ xx = dx xx X xxx dx.

50 Maurer Cartan Forms The Maurer Cartan forms for the diffeomorphism pseudo-group are the right-invariant one-forms on the diffeomorphism jet groupoid D ( ). Key observation: The target coordinate functions Z a are right-invariant. Thus, when we decompose dz a = σ a + µ a horizontal contact the two constituents are also right-invariant.

51 µ b A = DA Z µb = D Z a 1 D Z anµ b b = 1,..., m, #A 0 Invariant horizontal forms: σ a = d M Z a = m b = 1 Z a b dzb Invariant total differentiation (dual operators): D Z a = m b = 1 ( Z a b ) 1 D z b Thus, the invariant contact forms are obtained by invariant differentiation of the order zero contact forms: µ b = d G Z b = Θ b = dz b m a = 1 Z b a dza

52 One-dimensional case: M = R Contact forms: Θ = dx X x dx Θ x = D x Θ = dx x X xx dx Θ xx = D 2 x Θ = dx xx X xxx dx Right-invariant horizontal form: σ = d M X = X x dx Invariant differentiation: D X = 1 X x D x

53 Invariant contact forms: µ = Θ = dx X x dx µ X = D X µ = Θ x X x = dx x X xx dx X x µ XX = D 2 X µ = X x Θ xx X xx Θ x X 3 x µ n = D n X µ = X x dx xx X xx dx x + (X2 xx X x X xxx ) dx X 3 x.

54 The Structure Equations for the Diffeomorphism Pseudo group Maurer Cartan series: dµ b A = C b,b,c A,c,d µ b [ H ] = A µc B µd C 1 A! µb A HA H = (H 1,..., H m ) formal parameters dµ[ H ] = µ[ H ] ( µ[ H ] dz ) dσ = dµ[ 0 ] = µ[ 0 ] σ

55 The Structure Equations for the Diffeomorphism Pseudo group Maurer Cartan series: dµ b A = C b,b,c A,c,d µ b [ H ] = A µc B µd C 1 A! µb A HA H = (H 1,..., H m ) formal parameters dµ[ H ] = µ[ H ] ( µ[ H ] dz ) dσ = dµ[ 0 ] = µ[ 0 ] σ

56 Structure equations: One-dimensional case: M = R dσ = µ X σ dµ[ H ] = dµ dh [ H ] (µ[ H ] dz) where σ = X x dx = dx µ µ[ H ] = µ + µ X H µ XX H2 + µ[ H ] dz = σ + µ X H µ XX H2 + dµ[ H ] dh = µ X + µ XX H µ XXX H2 +

57 In components: dσ = µ 1 σ dµ n = µ n+1 σ + = σ µ n+1 n 1 i = 0 [ n+1 2 ] j = 1 ( ) n i n 2j + 1 n + 1 µ i+1 µ n i ( n + 1 j ) µ j µ n+1 j. = Cartan

58 The Maurer Cartan Forms for a Lie Pseudo-group The Maurer Cartan forms for G are obtained by restricting the diffeomorphism Maurer Cartan forms σ a, µ b A to G( ) D ( ). The resulting one-forms are no longer linearly independent.

59 Theorem. The Maurer Cartan forms on G ( ) satisfy the invariant infinitesimal determining equations L(... Z a... µ b A... ) = 0 ( ) obtained from the infinitesimal determining equations by replacing L(... z a... ζ b A... ) = 0 ( ) source variables z a by target variables Z a derivatives of vector field coefficients ζ b A by right-invariant Maurer Cartan forms µ b A

60 The Structure Equations for a Lie Pseudo-group Theorem. The structure equations for the pseudo-group G are obtained by restricting the universal diffeomorphism structure equations dµ[ H ] = µ[ H ] ( µ[ H ] dz ) to the solution space of the linearized involutive system L(... Z a,... µ b A,... ) = 0.

61 The Korteweg devries Equation u t + u xxx + uu x = 0 Diffeomorphism Maurer Cartan forms: µ t, µ x, µ u, µ t T, µt X, µt U, µx T,..., µu U, µt T T, µt T X,...

62 Infinitesimal determining equations: τ x = τ u = ξ u = ϕ t = ϕ x = 0 ϕ = ξ t 2 3 uτ t ϕ u = 2 3 τ t = 2 ξ x τ tt = τ tx = τ xx = = ϕ uu = 0 Maurer Cartan determining equations: µ t X = µt U = µx U = µu T = µu X = 0, µ u = µ x T 2 3 Uµt T, µu U = 2 3 µt T = 2 µx X, µ t T T = µt T X = µt XX = = µu UU =... = 0.

63 Basis (dim G KdV = 4): µ 1 = µ t, µ 2 = µ x, µ 3 = µ u, µ 4 = µ t T. Substituting into the full diffeomorphism structure equations yields the structure equations for g KdV : dµ 1 = µ 1 µ 4, dµ 2 = µ 1 µ U µ1 µ µ2 µ 4, dµ 3 = 2 3 µ3 µ 4, dµ 4 = 0. dµ i = C i jk µj µ k

64 dµ 1 = µ 1 µ 4, dµ 2 = µ 1 µ U µ1 µ µ2 µ 4, dµ 3 = 2 3 µ3 µ 4, dµ 4 = 0. In general, the pseudo-group structure equations live on the principal bundle G ( ) ; if G is a finite-dimensional Lie group, then G ( ) M G, and the usual Lie group structure equations are found by restriction to the target fibers {Z = c} G. Note that the constructed basis µ 1,..., µ r of g might vary from fiber to fiber.

65 Lie Kumpera Example Linearized determining system X = f(x) U = u f (x) ξ x = ϕ u ξ u = 0 ϕ u = ϕ u

66 Maurer Cartan forms: σ = u U dx = f x dx, τ = U x dx + U u du = u f xx dx + f x du f x 2 µ = dx U u dx = df f x dx, ν = du U x dx U u du = u f x 2 ( df x f xx dx ) µ X = du u du U x dx U = df x f xx dx f x, µ U = 0 ν X = U u (du x U xx dx) U x u (du U x dx) ν U = du u + du U x dx U = u f x 3 ( df xx f xxx dx ) + u f xx f x 4 ( df x f xx dx ) = df x f xx dx f x

67 First order linearized determining equations: ξ x = ϕ u ξ u = 0 ϕ u = ϕ u First order Maurer Cartan determining equations: µ X = ν U µ U = 0 ν U = ν U First order structure equations: dµ = dσ = ν σ U, dν = ν X σ ν τ U dν X = ν XX σ ν X (τ + 2 ν) U

68 Action of Pseudo-groups on Submanifolds a.k.a. Solutions of Differential Equations G Lie pseudo-group acting on p-dimensional submanifolds: N = {u = f(x)} M For example, G may be the symmetry group of a system of differential equations (x, u (n) ) = 0 and the submanifolds the graphs of solutions u = f(x).

69 Prolongation J n = J n (M, p) n th order submanifold jet bundle Local coordinates : z (n) = (x, u (n) ) = (... x i... u α J... ) Prolonged action of G (n) on submanifolds: (x, u (n) ) (X, Û (n) ) Coordinate formulae: Û α J = F J α (x, u(n), g (n) ) = Implicit differentiation.

70 Differential Invariants A differential invariant is an invariant function I : J n R for the prolonged pseudo-group action Invariant differential operators: I(g (n) (x, u (n) )) = I(x, u (n) ) = curvature, torsion,... D 1,..., D p = arc length derivative I(G) the algebra of differential invariants

71 The Basis Theorem Theorem. The differential invariant algebra I(G) is locally generated by a finite number of differential invariants I 1,..., I l and p = dim S invariant differential operators D 1,..., D p meaning that every differential invariant can be locally expressed as a function of the generating invariants and their invariant derivatives: D J I κ = D j1 D j2 D jn I κ. = Lie groups: Lie, Ovsiannikov = Lie pseudo-groups: Tresse, Kumpera, Kruglikov Lychagin, Muñoz Muriel Rodríguez, Pohjanpelto O

72 Key Issues Minimal basis of generating invariants: I 1,..., I l Commutation formulae for the invariant differential operators: [ D j, D k ] = p i=1 Syzygies (functional relations) among Y i jk D i = Non-commutative differential algebra the differentiated invariants: Φ(... D J I κ... ) 0 = Codazzi relations

73 Computing Differential Invariants The infinitesimal method: v(i) = 0 for every infinitesimal generator v g = Requires solving differential equations. Moving frames. (Cartan; PJO Fels Pohjanpelto ) Completely algebraic. Can be adapted to arbitrary group and pseudo-group actions. Describes the complete structure of the differential invariant algebra I(G) using only linear algebra & differentiation! Prescribes differential invariant signatures for equivalence and symmetry detection.

74 Moving Frames In the finite-dimensional Lie group case, a moving frame is defined as an equivariant map ρ (n) : J n G = All classical moving frames can be thus interpreted.

75 However, we do not have an appropriate abstract object to represent our pseudo-group G. Consequently, the moving frame will be an equivariant section ρ (n) : J n H (n) of the pulled-back pseudo-group jet groupoid: G (n) H (n) M J n.

76 Moving Frames for Pseudo Groups Definition. A (right) moving frame of order n is a rightequivariant section ρ (n) : V n H (n) defined on an open subset V n J n. = Groupoid action. Proposition. A moving frame of order n exists if and only if G (n) acts freely and regularly.

77 Freeness For Lie group actions, freeness means no isotropy. For infinitedimensional pseudo-groups, this definition cannot work, and one must restrict to the transformation jets of order n, using the n th order isotropy subgroup: Definition. G (n) = { g (n) G (n) z (n) z g (n) z (n) = z (n) } At a jet z (n) J n, the pseudo-group G acts freely if G (n) z (n) = { 1 (n) z } locally freely if G (n) z (n) is a discrete subgroup of G (n) z the orbits have dim = r n = dim G (n) z

78 Persistence of Freeness Theorem. If n 1 and G (n) acts locally freely at z (n) J n, then it acts locally freely at any z (k) J k with π n k(z(k) ) = z (n) for all k > n.

79 The Normalization Algorithm To construct a moving frame : I. Compute the prolonged pseudo-group action u α K U α K = F α K (x, u(n), g (n) ) by implicit differentiation. II. Choose a cross-section to the pseudo-group orbits: u α κ J κ = c κ, κ = 1,..., r n = fiber dim G (n)

80 III. Solve the normalization equations U α κ J κ = F α κ J κ (x, u (n), g (n) ) = c κ for the n th order pseudo-group parameters g (n) = ρ (n) (x, u (n) ) IV. Substitute the moving frame formulas into the unnormalized jet coordinates u α K = F K α(x, u(n), g (n) ). The resulting functions form a complete system of n th order differential invariants IK α (x, u(n) ) = FK α (x, u(n), ρ (n) (x, u (n) ))

81 Invariantization A moving frame induces an invariantization process, denoted ι, that projects functions to invariants, differential operators to invariant differential operators; differential forms to invariant differential forms, etc. Geometrically, the invariantization of an object is the unique invariant version that has the same cross-section values. Algebraically, invariantization amounts to replacing the group parameters in the transformed object by their moving frame formulas.

82 Invariantization In particular, invariantization of the jet coordinates leads to a complete system of functionally independent differential invariants: ι(x i ) = H i Phantom differential invariants: ι(u α J ) = Iα J I α κ J κ = c κ The non-constant invariants form a functionally independent generating set for the differential invariant algebra I(G)

83 Replacement Theorem I(... x i... u α J... ) = ι( I(... xi... u α J... ) ) = I(... H i... I α J... ) Differential forms = invariant differential forms ι( dx i ) = ω i i = 1,..., p Differential operators = invariant differential operators ι ( D x i ) = D i i = 1,..., p

84 Recurrence Formulae Invariantization and differentiation do not commute The recurrence formulae connect the differentiated invariants with their invariantized counterparts: D i I α J = Iα J,i + M α J,i = M α J,i correction terms

85 Once established, the recurrence formulae completely prescribe the structure of the differential invariant algebra I(G) thanks to the functional independence of the non-phantom normalized differential invariants. The recurrence formulae can be explicitly determined using only the infinitesimal generators and linear differential algebra!

86 Korteweg devries Equation Prolonged Symmetry Group Action: T = e 3λ 4 (t + λ 1 ) X = e λ 4 (λ 3 t + x + λ 1 λ 3 + λ 2 ) U = e 2λ 4 (u + λ 3 ) U T = e 5λ 4 (u t λ 3 u x ) U X = e 3λ 4 u x U T T = e 8λ 4 (u tt 2λ 3 u tx + λ 3 2 u xx ) U T X = D X D T U = e 6λ 4 (u tx λ 3 u xx ) U XX = e 4λ 4 u xx.

87 Cross Section: T = e 3λ 4 (t + λ 1 ) = 0 X = e λ 4 (λ 3 t + x + λ 1 λ 3 + λ 2 ) = 0 U = e 2λ 4 (u + λ 3 ) = 0 U T = e 5λ 4 (u t λ 3 u x ) = 1 Moving Frame: λ 1 = t, λ 2 = x, λ 3 = u, λ 4 = 1 5 log(u t + uu x )

88 Moving Frame: λ 1 = t, λ 2 = x, λ 3 = u, λ 4 = 1 5 log(u t + uu x ) Invariantization: ι(u K ) = U K λ1 = t,λ 2 = x,λ 3 = u,λ 4 =log(u t +uu x )/5 Phantom Invariants: H 1 = ι(t) = 0 H 2 = ι(x) = 0 I 00 = ι(u) = 0 I 10 = ι(u t ) = 1

89 Normalized differential invariants: u I 01 = ι(u x ) = x (u t + uu x ) 3/5 I 20 = ι(u tt ) = u tt + 2uu tx + u2 u xx (u t + uu x ) 8/5 I 11 = ι(u tx ) = I 02 = ι(u xx ) = u tx + uu xx (u t + uu x ) 6/5 u xx (u t + uu x ) 4/5 I 03 = ι(u xxx ) =. u xxx u t + uu x

90 Invariantization: ι( F (t, x, u, u t, u x, u tt, u tx, u xx,... ) ) = F (ι(t), ι(x), ι(u), ι(u t ), ι(u x ), ι(u tt ), ι(u tx ), ι(u xx ),... ) = F (H 1, H 2, I 00, I 10, I 01, I 20, I 11, I 02,... ) = F (0, 0, 0, 1, I 01, I 20, I 11, I 02,... ) Replacement Theorem: 0 = ι(u t + u u x + u xxx ) = 1 + I 03 = u t + uu x + u xxx u t + uu x. Invariant horizontal one-forms: ω 1 = ι(dt) = (u t + uu x ) 3/5 dt, ω 2 = ι(dx) = u(u t + uu x ) 1/5 dt + (u t + uu x ) 1/5 dx.

91 Invariant differential operators: D 1 = ι(d t ) = (u t + uu x ) 3/5 D t + u(u t + uu x ) 3/5 D x, D 2 = ι(d x ) = (u t + uu x ) 1/5 D x. Commutation formula: [ D 1, D 2 ] = I 01 D 1 Recurrence formulae: D 1 I 01 = I I I 01 I 20, D 2 I 01 = I I I 01 I 11, D 1 I 20 = I I I 01 I I2 20, D 2 I 20 = I I 01 I I2 01 I I 11 I 20, D 1 I 11 = I 21 + I I 01 I I 11 I 20, D 2 I 11 = I 12 + I 01 I I2 01 I I2 11, D 1 I 02 = I I 01 I I 02 I 20, D 2 I 02 = I I2 01 I I 02 I 11,..

92 Generating differential invariants: I 01 = ι(u x ) = Fundamental syzygy: u x (u t + uu x ) 3/5, I 20 = ι(u tt ) = u tt + 2uu tx + u2 u xx (u t + uu x ) 8/5. D 2 1 I I 01 D 1 I 20 D 2 I 20 + ( 1 5 I I 01) D1 I 01 D 2 I I 01 I I2 01 I I3 01 = 0.

93 Horizontal coframe Lie Tresse Kumpera Example X = f(x), Y = y, U = u f (x) d H X = f x dx, d H Y = dy, Implicit differentiations D X = 1 f x D x, D Y = D y.

94 Prolonged pseudo-group transformations on surfaces S R 3 X = f Y = y U = u f x U X = u x f 2 x U XX = u xx f 3 x u f xx f 3 x 3u x f xx f 4 x U Y = u y f x u f xxx f 4 x + 3u f 2 xx f 5 x U XY = u xy f 2 x Coordinate cross-section u y f xx f 3 x U Y Y = u yy f x = action is free at every order. X = f = 0, U = u f x = 1, U X = u x f 2 x u f xx f 3 x = 0, U XX = = 0.

95 Moving frame f = 0, f x = u, f xx = u x, f xxx = u xx. Differential invariants U Y J = u y u U XY J 1 = uu xy u x u y u 3 U Y Y J 2 = u yy u Invariant horizontal forms d H X = f x dx u dx, d H Y = dy dy, Invariant differentiations D 1 = 1 u D x D 2 = D y

96 Higher order differential invariants: D m 1 Dn 2 J Recurrence formulae: J,1 = D 1 J = uu xy u x u y u 3 = J 1, J,2 = D 2 J = uu yy u2 y u 2 = J 2 J 2. D 1 J = J 1, D 2 J = J 2 J 2, D 1 J 1 = J 3, D 2 J 1 = J 4 3 J J 1, D 1 J 2 = J 4, D 2 J 2 = J 5 J J 2,

97 The Master Recurrence Formula p d H IJ α = ( D i IJ α ) ωi = i=1 p i=1 I α J,i ωi + ψ α J where ψ α J = ι( ϕ α J ) = Φα J (... Hi... IJ α... ;... γb A... ) are the invariantized prolonged vector field coefficients, which are particular linear combinations of γa b = ι(ζb A ) invariantized Maurer Cartan forms prescribed by the invariantized prolongation map. The invariantized Maurer Cartan forms are subject to the invariantized determining equations: L(H 1,..., H p, I 1,..., I q,..., γ b A,... ) = 0

98 d H I α J = p i=1 I α J,i ωi + ψ α J (... γb A... ) Step 1: Solve the phantom recurrence formulas 0 = d H I α J = p i=1 I α J,i ωi + ψ α J (... γb A... ) for the invariantized Maurer Cartan forms: γ b A = p i=1 J b A,i ωi ( ) Step 2: Substitute ( ) into the non-phantom recurrence formulae to obtain the explicit correction terms.

99 Only uses linear differential algebra based on the specification of cross-section. Does not require explicit formulas for the moving frame, the differential invariants, the invariant differential operators, or even the Maurer Cartan forms!

100 The Korteweg devries Equation (continued) Recurrence formula: di jk = I j+1,k ω 1 + I j,k+1 ω 2 + ι(ϕ jk ) Invariantized Maurer Cartan forms: ι(τ) = λ, ι(ξ) = µ, ι(ϕ) = ψ = ν, ι(τ t ) = ψ t = λ t,... Invariantized determining equations: λ x = λ u = µ u = ν t = ν x = 0 ν = µ t ν u = 2 µ x = 2 3 λ t λ tt = λ tx = λ xx = = ν uu = = 0 Invariantizations of prolonged vector field coefficients: ι(τ) = λ, ι(ξ) = µ, ι(ϕ) = ν, ι(ϕ t ) = I 01 ν 5 3 λ t, ι(ϕ x ) = I 01 λ t, ι(ϕ tt ) = 2I 11 ν 8 3 I 20 λ t,...

101 Phantom recurrence formulae: 0 = d H H 1 = ω 1 + λ, 0 = d H H 2 = ω 2 + µ, 0 = d H I 00 = I 10 ω 1 + I 01 ω 2 + ψ = ω 1 + I 01 ω 2 + ν, 0 = d H I 10 = I 20 ω 1 + I 11 ω 2 + ψ t = I 20 ω 1 + I 11 ω 2 I 01 ν 5 3 λ t, = Solve for λ = ω 1, µ = ω 2, ν = ω 1 I 01 ω 2, λ t = 3 5 (I 20 + I 01 )ω (I 11 + I2 01 )ω2. Non-phantom recurrence formulae: d H I 01 = I 11 ω 1 + I 02 ω 2 I 01 λ t, d H I 20 = I 30 ω 1 + I 21 ω 2 2I 11 ν 8 3 I 20 λ t, d H I 11 = I 21 ω 1 + I 12 ω 2 I 02 ν 2I 11 λ t, d H I 02 = I 12 ω 1 + I 03 ω I 02 λ t,.

102 D 1 I 01 = I I I 01 I 20, D 2 I 01 = I I I 01 I 11, D 1 I 20 = I I I 01 I I2 20, D 2 I 20 = I I 01 I I2 01 I I 11 I 20, D 1 I 11 = I 21 + I I 01 I I 11 I 20, D 2 I 11 = I 12 + I 01 I I2 01 I I2 11, D 1 I 02 = I I 01 I I 02 I 20, D 2 I 02 = I I2 01 I I 02 I 11,..

103 Lie Tresse Kumpera Example (continued) X = f(x), Y = y, U = u f (x) Phantom recurrence formulae: 0 = dh = ϖ 1 + γ, 0 = di 10 = J 1 ϖ 2 + ϑ 1 γ 2, 0 = di 00 = J ϖ 2 + ϑ γ 1, 0 = di 20 = J 3 ϖ 2 + ϑ 3 γ 3, Solve for pulled-back Maurer Cartan forms: γ = ϖ 1, γ 2 = J 1 ϖ 2 + ϑ 1, Recurrence formulae: dy = ϖ 2 γ 1 = J ϖ 2 + ϑ, γ 3 = J 3 ϖ 2 + ϑ 3, dj = J 1 ϖ 1 + (J 2 J 2 ) ϖ 2 + ϑ 2 J ϑ, dj 1 = J 3 ϖ 1 + (J 4 3 J J 1 ) ϖ 2 + ϑ 4 J ϑ 1 J 1 ϑ, dj 2 = J 4 ϖ 1 + (J 5 J J 2 ) ϖ 2 + ϑ 5 J 2 ϑ,

104 Gröbner Basis Approach Identify the cross-section variables with the complementary monomials to a certain algebraic module J, which is the pull-back of the symbol module of the pseudo-group under a certain explicit linear map. = Compatible term ordering. = Algebraic specification of compatible moving frames of all orders n > n.

105 Theorem. Suppose G acts freely at order n. Then a system of generating differential invariants is contained in the non-phantom normalized differential invariants of order n and those differential invariants corresponding to a Gröbner basis for the module J >n.

106 The Symbol Module Linearized determining equations L(z, ζ (n) ) = 0 t = (t 1,..., t m ), T = (T 1,..., T m ) { T = P (t, T ) = m a = 1 P a (t) T a } R[t] R m R[t, T ] I T symbol module s = (s 1,..., s p ), S = (S 1,..., S q ), { Ŝ = T (s, S) = q α = 1 T α (s) S α } R[s] R q R[s, S]

107 Define the linear map s i = β i (t) = t i + q α = 1 u α i t p+α, i = 1,..., p, S α = B α (T ) = T p+α p i=1 u α i T i, α = 1,..., q. Prolonged symbol module: J = (β ) 1 (I) N leading monomials s J S α = normalized differential invariants I α J K complementary monomials s K S β = phantom differential invariants I β K

108 The Symbol Module Vector field: Vector field jet: v = m a = 1 ζ b (z) z b j v ζ ( ) = (... ζ b A... ) ζ b A = #A ζ b z A = k ζ b z a 1 z a k Determining Equations for v g L(z;... ζ b A... ) = 0 ( )

109 Duality t = (t 1,..., t m ) T = (T 1,..., T m ) Polynomial module: { T = P (t, T ) = Dual pairing: m a = 1 T P a (t) T a } R[t] R m R[t, T ] (J T M z ) j v ; t A T b = ζ b A.

110 Each polynomial τ(z; t, T ) = m b = 1 #A n h b A (z) ta T b induces a linear partial differential equation L(z, ζ (n) ) = j v ; τ(z; t, T ) T = m b = 1 #A n h b A (z) ζb A = 0

111 The Linear Determining Equations Annihilator: L = (J g) Determining Equations j v ; τ = 0 for all η L v g Symbol = highest degree terms: Σ[L(z, ζ (n) )] = Λ[τ(z; t, T )] = m b = 1 #A=n h b A (z) ta T b.

112 Symbol submodule: I = Λ(L) = Formal integrability (involutivity)

113 Prolonged Duality Prolonged vector field: v ( ) = p i=1 ξ i (x, u) x i + α,j ϕ α J (x, u(k) ) u α J s = ( s 1,..., s p ), s = (s 1,..., s p ), S = (S 1,..., S q ) Prolonged polynomial module: { Ŝ = σ(s, S, s) = p c i s i + q i=1 α = 1 σ α (s) S α } R p (R[s] R q )

114 Ŝ T J z ( ) Dual pairing: v ( ) ; s i = ξ i v ( ) ; S α = Q α = ϕ α p i=1 u α i ξi v ( ) ; s J S α = ϕ α J = Φα J (u(n) ; ζ (n) )

115 Algebraic Prolongation Prolongation of vector fields: p : J g g ( ) j v v ( ) Dual prolongation map: p : S T j v ; p (σ) = p(j v) ; σ = v ( ) ; σ On the symbol level, p is algebraic

116 Prolongation Symbols Define the linear map β : R 2m R m s i = β i (t) = t i + q α = 1 u α i t p+α, i = 1,..., p, S α = B α (T ) = T p+α p i=1 u α i T i, α = 1,..., q. Pull-back map β [ σ(s 1,..., s p, S 1,..., S q ) ] = σ( β 1 (t),..., β p (t), B 1 (T )..., B q (T ) )

117 Lemma. are The symbols of the prolonged vector field coefficients Σ(ξ i ) = T i Σ( ϕ α ) = T α+p Σ(Q α ) = β (S α ) = B α (T ) Σ( ϕ α J ) =β (s J S α ) = β (s j1 s jn S α ) = β j1 (t) β jn (t) B α (T )

118 Prolonged annihilator: Z = (p ) 1 L = (g ( ) ) v ( ) ; σ = 0 for all v g σ Z Prolonged symbol subbundle: U = Λ(Z) J (M, p) S Prolonged symbol module: J = (β ) 1 (I) Warning: : But U J U n = J n when n > n n order of freeness.

119 Algebraic Recurrence Polynomial: σ(i (k) ; s, S) = α,j h a J (I(k) ) s J S α Ŝ Differential invariant: I σ = α,j h a J (I(k) ) I α J

120 Recurrence: D i I σ = I Di σ I s i σ + R i,σ order I σ = n σ J n, n > n = order I Di σ = n + 1 order R i,σ n

121 I symbol module Algebra = Invariants determining equations for g M T / I complementary monomials t A T b pseudo-group parameters Maurer Cartan forms N leading monomials s J S α normalized differential invariants I α J

122 K = S / N complementary monomials s K S β cross-section coordinates u β K = cβ K phantom differential invariants I β K J = (β ) 1 (I) Freeness: β : K M

123 Generating Differential Invariants Theorem. The differential invariant algebra is generated by differential invariants that are in one-to-one correspondence with the Gröbner basis elements of the prolonged symbol module plus, possibly, a finite number of differential invariants of order n.

124 Syzygies Theorem. Every differential syzygy among the generating differential invariants is either a syzygy among those of order n, or arises from an algebraic syzygy among the Gröbner basis polynomials in J.