Generating Compatibility Conditions and General Relativity

Abstract

The search for the generating compatibility conditions (CC) of a given operator is a very recent problem met in general relativity in order to study the Killing operator for various standard useful metrics. Accordingly, this paper can be considered as a natural continuation of a previous paper recently published in JMP under the title Minkowski, Schwarschild and Kerr metrics revisited. In particular, we prove that the intrinsic link existing between the lack of formal exactness of an operator sequence on the jet level, the lack of formal exactness of its corresponding symbol sequence and the lack of formal integrability (FI) of the initial operator is of a purely homological nature as it is based on the long exact connecting sequence provided by the so-called snake lemma in homological algebra. It is therefore quite difficult to grasp it in general and even more difficult to use it on explicit examples. It does not seem that any one of the results presented in this paper is known as most of the other authors who studied the above problem of computing the total number of generating CC are confusing this number with the degree of generality introduced by A. Einstein in his 1930 letters to E. Cartan. One of the motivating examples that we provide is so striking that it is even difficult to imagine that such an example could exist. We hope this paper could be used as a source of testing examples for future applications of computer algebra in general relativity and, more generally, in mathematical physics.

Share and Cite:

Pommaret, J. (2019) Generating Compatibility Conditions and General Relativity. Journal of Modern Physics, 10, 371-401. doi: 10.4236/jmp.2019.103025.

1. Introduction

If X is a manifold of dimension n with local coordinates ( x 1 , , x n ) , let us introduce the tangent bundle T = T ( X ) and the cotangent bundle T * = T * ( X ) , the q-symmetric tensor bundle S q T * and the bundle r T * of r-forms. In General Relativity, there may be different solutions of Einstein equations in vacuum like the Minkowski, the Schwarzschild and the Kerr metrics for example. For fixing the notations and with more details, if ω S 2 T * is a nondegenerate metric, that is d e t ( ω ) 0 , and if j q denotes all the derivatives of an object up to order q, we may construct the Christoffel symbols γ through the Levi-Civita isomorphism ( ω , γ ) j 1 ( ω ) and, using the language of jet bundles, ( ω , γ ) is a section of J 1 ( S 2 T * ) that will be simply written ( ω , γ ) J 1 ( S 2 T * ) . Then we can introduce the well-known Riemann tensor ρ = ( ρ l , i j k ) 2 T * T * T with ρ k l , i j + ρ l k , i j = 0 after lowering the upper index by means of ω and δ ρ = 0 where δ : 2 T * T * T 3 T * T is the Spencer δ -map. Introducing the Ricci tensor ρ i j = ρ i , r j r = ρ j i or the Einstein tensor ϵ i j = ρ i j 1 2 ω i j ω r s ρ r s , the 10 non-linear Einstein equations are described by ϵ i j = 0 or, equivalently, by ρ i j = 0 when n = 4 .

Now, if E is a fibered manifold over X with fiber dimension m and local coordinates ( x i , y k ) with i = 1 , , n and k = 1 , , m , we may introduce the tangent bundle T ( E ) over E with local coordinates ( x , y , u , v ) and the vertical bundle V ( E ) with local coordinates ( x , y , u = 0 , v ) = ( x , y , v ) which are both vector bundles over E . We shall denote by the capital letters Ω S 2 T * , Γ S 2 T * T , R 2 T * T * T , the respective linearizations of ω , γ , ρ which are sections of the respective vertical bundles. Introducing the Lie derivative L of geometric objects, it is therefore possible to introduce the corresponding first order Killing operator D : T Ω S 2 T * : ξ L ( ξ ) ω , the first order Christoffel operator Ω Γ S 2 T * T in such a way that C h r i s t o f f e l K i l l i n g : T Γ S 2 T * T : ξ L ( ξ ) γ and the second order Riemann operator Ω R 2 T * T * T in such a way that R i e m a n n K i l l i n g : T R 2 T * T * T : ξ L ( ξ ) ρ both with its contraction Ω S 2 T * called Ricci operator. For example, it is known that 2 ω r k γ i j k = ( i ω r j + j ω i r r ω i j ) that we shall write simply, using formal notations, 2 ω γ = ( ω + ω ω ) and thus 2 ω Γ + 2 γ Ω = ( d Ω + d Ω d Ω ) . We have proved in ( [1] [2] [3] [4] ) that the so-called gravitational waves equations are nothing else than a d ( R i c c i ) by introducing the formal adjoint operator. It is important to notice that the Einstein operator Ω E i j = R i j 1 2 ω i j ω r s R r s is self-adjoint with 6 terms though the Ricci operator is not with only 4 terms. Recently, many physicists (See [5] [6] [7] [8] [9] ) have tried to construct the compatibility conditions (CC) of the Killing operator for various types of background metrics, in particular the three ones already quoted, namely an operator D 1 : S 2 T * F 1 such that D 1 Ω = 0 generates the CC of D ξ = Ω . We have proved in the above references the following crucial results:

・ These CC may contain a certain number of second and third order CC. It is therefore crucial in actual practice to select the successive generating CC of order 1,2,3, till we stop because of noetherian arguments ( [10] ).

・ These CC only depend on the Lie algebra structure (dimension of the solution space and commutation relations) of the corresponding Killing operator, which, even though it is finite dimensional with dimension n ( n + 1 ) / 2 that is 10 obtained for the Minkowski metric, may have dimension 4 for the Schwarzschild metric and dimension 2 for the Kerr metric.

・ The only two canonical sequences that can be constructed from an operator or a system, namely the Janet and Spencer sequences, are structurally quite different. Indeed, the Janet bundles F 0 , , F n appearing in the Janet sequence are concerned wit geometric objects like ω , γ , ρ , while the Spencer bundles C 0 , , C n are far from being related with geometric objects, the simplest example being C 0 = R q J q ( E ) . In the case of Lie equations considered, the central concept is not the system but rather the group as it can be seen at once from the construction of the Vessiot structure equations ( [3] [11] [12] [13] ).

The authors who have studied these questions had in mind that the total number of generating CC could be considered as a kind of “differential transcendence degree”, also called “degree of generality” by A. Einstein in his letters to E. Cartan of 1930 on absolute parallelism ( [14] ), the modern definition being that of the “differential rank” ( [10] [12] [15] [16] ). We must say that Cartan, being unable to explain to Einstein his theory of exterior systems, just copied the work of Janet published in 1920 ( [17] ) in his letters to Einstein, published later on as the only paper he wrote on the PD approach, but without ever quoting Janet who suffered a lot from this behavior and had to turn to mechanics.

Such a result will be obtained in the framework of differential modules as its explanation in the framework of differential systems is much more delicate and technical ( [10] [12] [18] ).

First of all, with our previous assumptions, D = K [ d ] is a noetherian domain and we can restrict our study to finitely generated differential modules which are therefore finitely presented (See [14] for more details). Let thus M be defined by a finite free presentation giving rise to the long exact sequence:

0 L D p D D m p M 0

where the differential operator D is acting on the right by composition with action law ( P , D ) P D , P D , p is the canonical residual projection and L = k e r ( D ) D p . The image i m ( D ) D m is called the differential module of equations and is thus finitely generated because D is a noetherian differential domain.

DEFINITION 1.1: The differential rank r k D ( M ) over D of a differential module M is the differential rank over D of the maximum free differential submodule F of M and we have the short exact sequence 0 F M T 0 where T = M / F is a torsion module over D. In particular, if F D r , then r k D ( M ) = r .

The following useful proposition proves the additivity property of the differential rank and is used in the next two corollaries ( [3] [10] [12] ):

PROPOSITION 1.2: If we have a short exact sequence 0 M f M g M 0 of differential modules, then r k D ( M ) = r k D ( M ) + r k D ( M ) .

COROLLARY 1.3: If D is a linear partial differential operator with coefficients in a differential field K and a d ( D ) is the formal adjoint that can be obtained formally or through an integration by parts, then r k D ( D ) = r k D ( a d ( D ) ) .

COROLLARY 1.4: (Euler-Poincaré characteristic) For any finite free differential resolution of a differential module M, then r k D ( M ) is equal to the alternate sum χ D ( M ) of differential ranks of the free differential modules of the resolution.

We obtain therefore r k D ( L ) p + m r k D ( M ) = 0 and it follows from noetherian arguments that the differential module L D p is finitely generated but not free in general and we may look for a minimum number of generators which may be differentially dependent in general as we shall see in the next examples. It thus remains to provide examples of such computations showing that these two numbers are not related and must therefore be found totally independently in general, apart from the very exceptional situation met when there is only a single generating CC.

In actual practice, working in the system framework, starting with a system R q J q ( E ) of order q on E and introducing the canonical projection Φ : J q ( E ) F 0 = J q ( E ) / R q , we shall construct for each r 0 a family of FI systems B r = i m ( ρ r ( Φ ) ) J r ( F 0 ) such that B r + 1 ρ 1 ( B r ) = J 1 ( B r ) J r + 1 ( F 0 ) J 1 ( J r ( F 0 ) ) projects onto B r , that is B r + 1 is defined by more generating PD equations than the ones defining B r both with its prolongations, and start to get equality when r is large enough in the projective limit B B r + 1 B r B 1 F 0 0 . The striking result is that there may be gaps in the procedure, that is we shall even provide a tricky example where one can have a single generating CC of order 3, then no new generating CC of order 4 and 5, but suddenly a new generating CC of order 6 ending the procedure. We do not believe that such situations were even known to exist.

2. Motivating Examples

We provide below three examples, pointing out that it is quite difficult to exhibit such examples.

EXAMPLE 2.1: With n = 3, m = d i m ( E ) = 2, d i m ( F 0 ) = 5 and K = while keeping an upper index for any unknown, let us consider the following system R 1 J 1 ( E ) with d i m ( R 1 ) = 3 because p a r 1 = { ξ 1 , ξ 2 , ξ 1 2 } and corresponding Janet tabular:

{ Φ 5 ξ 3 2 = 0 Φ 4 ξ 3 1 = 0 Φ 3 ξ 2 2 = 0 Φ 2 ξ 2 1 + ξ 1 2 = 0 Φ 1 ξ 1 1 = 0 1 2 3 1 2 3 1 2 1 2 1 ×

It is easy to check that all the second order jets vanish and that the general solution { ξ 1 = a x 2 + b , ξ 2 = c x 1 + d | a + c = 0 } depends on 3 arbitrary constants. As the non-multiplicative variable written with the sign ×cannot be used, the symbol g 1 is not involutive because it is finite type with g 2 = 0 . This system is trivially FI because it is made by homogeneous PD equations. We have the following commutative diagrams:

0 0 0 0 S 3 T * E S 2 T * F 0 h 2 0 0 R 3 J 3 ( E ) J 2 ( F 0 ) Q 2 0 π 2 3 π 1 2 0 R 2 J 2 ( E ) J 1 ( F 0 ) Q 1 0 0 0 0 0

0 0 0 0 20 30 10 0 0 3 40 50 13 0 π 2 3 π 1 2 0 3 20 20 3 0 0 0 0 0

0 0 0 S 3 T * E S 2 T * F 0 h 2 0 δ δ 0 T * S 2 T * E T * T * F 0 T * h 1 0 δ δ 0 2 T * g 1 2 T * T * E 2 T * F 0 0 δ δ 0 3 T * E = 3 T * E 0 0 0

0 0 0 20 30 10 0 δ δ 0 36 45 9 0 δ δ 0 3 18 15 0 δ δ 0 2 = 2 0 0 0

The next result points out the importance of the Spencer δ -cohomology. Indeed, we shall prove that the last symbol diagram is commutative and exact. In particular, the lower left map δ is surjective and thus the upper right induced map h 2 T * Q 1 is also surjective while these two maps have isomorphic kernels.

For this, we notice that the 3 components of 2 T * g 1 are { v 1,12 2 , v 1,13 2 , v 1,23 2 } and the map δ is described by the two linear equations:

w 123 1 v 1 , 23 1 + v 2 , 31 1 + v 3 , 12 1 = v 1 , 13 2 = 0 , w 123 2 v 1 , 23 2 + v 2 , 31 2 + v 3 , 12 2 = v 1 , 23 2 = 0

that is to say by two linearly independent equations. Accordingly, in the left column we have:

d i m ( H 2 ( g 1 ) ) = d i m ( Z 2 ( g 1 ) ) = d i m ( k e r ( δ ) ) = 1

An unusual snake-type diagonal chase left to the reader as an exercise proves that the induced map h 2 T * Q 1 is surjective with a kernel isomorphic to H 2 ( g 1 ) . This is indeed a crucial result because it also proves that the additional CC has only to do with the single second order component of the Riemann tensor in dimension 2, a striking result that could not even be imagined by standard methods. Moreover, we know that if a system R q J q ( E ) is FI, for example when it is homogeneous like in this case, and its symbol g q S q T * E is such that s is the smallest integer such that g q + s becomes 2-acyclic (or involutive), then the generating CC are of order at most s + 1 ( [3] [10] [12] ).

Collecting the above results, we find the 3 first order differentially independent generating CC coming from the Janet tabular and the additional single second order generating CC describing the 2-dimensional Riemann operator, that is the linearized Riemann tensor in the space ( x 1 , x 2 ) :

{ Ψ 4 d 22 Φ 1 + d 11 Φ 3 d 12 Φ 2 = 0 Ψ 3 d 3 Φ 3 d 2 Φ 5 = 0 Ψ 2 d 3 Φ 2 d 2 Φ 4 d 1 Φ 5 = 0 Ψ 1 d 3 Φ 1 d 1 Φ 4 = 0

An elementary computation provides the second order CC:

d 22 Ψ 1 + d 11 Ψ 3 d 12 Ψ 2 d 3 Ψ 4 = 0

The corresponding differential sequence written with differential modules over the ring D = K [ d 1 , d 2 , d 3 ] = K [ d ] is:

0 D 2 D 4 2 D 5 1 D 2 p M 0

where p is the canonical (residual) projection. We check indeed that 1 4 + 5 2 = 0 but this sequence is quite far from being even strictly exact. Of course, as R 2 is involutive, we may set C r = r T * R 2 and obtain the corresponding canonical second Spencer sequence which is induced by the Spencer operator:

0 Θ j 2 C 0 D 1 C 1 D 2 C 2 D 3 C 3 0

with dimensions:

0 Θ j 2 3 1 D 1 9 1 D 2 9 1 D 3 3 0

Proceeding inductively as we did for finding the second order CC, we may obtain by combinatorics the following formally exact sequence:

0 Θ 2 1 5 2 13 1 19 1 12 1 3 0

with Euler-Poincaré characteristic 2 5 + 13 19 + 12 3 = 0 but, as before, there is a matrix 260 × 280 at least and we doubt about the use of computer algebra, even on such an elementary example. With F 0 = J 1 ( E ) / R 1 , the starting long exact sequence used as a middle row of the first diagram with dimensions:

0 3 40 50 13 0

and we have 13 = ( 3 + 3 × 3 ) + 1 , that is three generating first order CC which are differentially independent, plus their 9 prolongations, plus one second order CC which is nevertheless not differentially independent. Hence we have a total number of 3 + 1 = 4 generating CC but this number has nothing to do with any differential transcendence degree because Ψ 4 is differentially algebraic over { Ψ 1 , Ψ 2 , Ψ 3 } .

We finally compute the corresponding (canonical) Janet sequence by quotient. For this, we must use the trivial second Spencer sequence:

0 E j 2 C 0 ( E ) D 1 C 1 ( E ) D 2 C 2 ( E ) D 3 C 3 ( E ) 0

namely:

0 2 2 20 1 40 1 30 1 8 0

with 2 20 + 40 30 + 8 = 0 . The (canonical) Janet sequence is thus:

0 Θ E 2 F 0 1 F 1 1 F 2 1 F 3 0

with now F 0 = J 2 ( E ) / R 2 and F r = C r ( E ) / C r , r = 0 , 1 , 2 , 3 and dimensions:

0 Θ 2 2 17 1 31 1 21 1 5 0

so that we have again 2 17 + 31 21 + 5 = 0 in a coherent way with the fact that r k D ( M ) = 0 .

EXAMPLE 2.2: With n = 3 , m = 1 , q = 2 and K = ( x 1 , x 2 , x 3 ) = ( x ) , let us consider the following linear inhomogeneous system:

y 33 x 2 y 1 = v , y 12 = u

・ Step 1: The symbol g 2 is defined by v 33 = 0 , v 12 = 0 may not be involutive or the coordinate system may not be δ -regular. However, changing linearly the local coordinates with x 1 x 1 , x 2 x 2 + x 1 , x 3 x 3 , we obtain the Janet tabular for g 2 :

{ v 33 = 0 v 22 + v 12 = 0 1 2 3 1 2

and thus the Janet tabular for g 3 :

{ v 333 = 0 v 233 = 0 v 223 + v 123 = 0 v 222 + v 122 = 0 v 133 = 0 v 122 + v 112 = 0 1 2 3 1 2 1 2 1 2 1 1

We let the reader check as an exercise that g 2 is not 2-acyclic by counting the dimensions in the long sequence:

0 g 4 δ T * g 3 δ 2 T * g 2 δ 3 T * T *

and that g 3 is involutive, thus 2-acyclic, with characters ( 0,0,4 ) . It follows that d i m ( g 2 ) = 6 2 = 4 , d i m ( g 3 ) = 0 + 0 + 4 = 4 , d i m ( g 4 ) = 4, . We obtain from the main theorem ρ r ( R 2 ( 1 ) ) = R 2 + r ( 1 ) . It is easy to check that R 2 ( 1 ) = R 2 with d i m ( R 2 ) = 8 , d i m ( R 3 ) = 8 + 4 = 12 , d i m ( R 4 ) = 11 + 4 = 15 , d i m ( R 5 ) = 56 39 = 13 + 4 = 17 but things are changing after that. As such a property is intrinsic, coming back to the original system of coordinates, we have after one more prolongation:

{ y 1233 = u 33 y 1233 + x 2 y 112 + y 11 = v 12 x 2 y 112 = x 2 u 1

and thus y 11 = u 33 v 12 x 2 u 1 . We may thus consider the new second order system R 2 = R 2 ( 2 ) R 2 with a strict inclusion and d i m ( R 2 ( 2 ) ) = 7 :

y 33 x 2 y 1 = v , y 12 = u , y 11 = w , w = u 33 v 12 x 2 u 1

We may start again with R 2 and study its symbol g 2 defined by the 3 linear equations with the following Janet tabular obtained after doing the same change of local coordinates as before:

{ v 33 = 0 v 22 + v 12 = 0 v 12 + v 11 = 0 1 2 3 1 2 1

This symbol is neither 2-acyclic nor involutive but its prolongation g 3 , defined by the 8 equations:

{ v 333 = 0 v 233 = 0 v 223 v 113 = 0 v 222 + v 111 = 0 v 133 = 0 v 123 + v 113 = 0 v 122 v 111 = 0 v 112 + v 111 = 0 1 2 3 1 2 1 2 1 2 1 1 1 1

is involutive with characters ( 0,0,2 ) and we may consider again the system:

{ y 33 x 2 y 1 = v y 12 = u y 11 = w = u 33 v 12 x 2 u 1

Instead of doing the same change of variables, writing out the system R 3 and study its formal inegrability with corresponding 9 + 11 = 20 CC for ( u , v , w ) , an elementary but tedious computation, we shall use a trick, knowing in advance that the generating CC must be of order 1 + 1 = 2 because g 2 had to get one prolongation in order to become involutive and thus 2-acyclic.

・ Step 2: It thus remains to find out the CC for ( u , v ) in the initial inhomogeneous system. As we have used two prolongations in order to exhibit R 2 , we have second order formal derivatives of u and v in the right members. Now, from the above argument, we have second order CC for the new right members and could hope therefore for a fourth order generating CC. The trick is to use the three different brackets of operators that can be obtained. We have in a formal way:

[ d 33 x 2 d 1 , d 12 ] y = y 1233 x 2 y 112 y 1233 + x 2 y 112 + y 11 = y 11 = u 33 x 2 u 1 v 12

[ d 2 , d 11 ] y = 0 = d 2 ( d 11 y ) d 1 ( d 12 y ) = ( u 233 v 122 x 2 u 12 u 1 ) u 1

brings the third order CC:

A u 233 v 122 x 2 u 12 2 u 1 = 0

[ d 33 x 2 d 1 , d 11 ] y = y 1133 x 2 y 111 y 1133 + x 2 y 111 = 0 = ( d 33 x 2 d 1 ) w v 11 = u 3333 v 1233 2 x 2 u 133 + x 2 v 112 v 11 + ( x 2 ) 2 u 11

brings the fourth order CC:

B u 3333 v 1233 2 x 2 u 133 + x 2 v 112 v 11 + ( x 2 ) 2 u 11 = 0

We have indeed the identity A 33 x 2 A 1 B 2 = 0 and thus ( A , B ) are differentially dependent, that is B is a new generating fourth order CC which is not a consequence of the prolongations of A. Again, the total number of generating CC, that is 1 + 1 = 2 , has nothing to do with the differential transcendence degree of the CC differential module which is d i m ( F 0 ) d i m ( E ) = 2 1 = 1 .

EXAMLE 2.3: With the same n = 3 , m = 1 , q = 2 and K = ( x 1 , x 2 , x 3 ) = ( x ) , we now prove that a slight change of the equations may provide quite important changes in the number and order of the CC. Such an example is the only one that we could have found in more than 40 years of computing CC in mathematics and applications. For this, let us consider the new system:

y 33 x 2 y 1 = v , y 22 = u

Before starting, we first notice that it is a prioiri not evident to discover that R = R is a finite dimensional vector space over K with d i m K ( R ) = 6 . However such a result can be obtained by direct integration (Compare to the Janet example treated in the introduction of [12] ).

・ Step 1: The symbol g 2 is defined by v 33 = 0 , v 22 = 0 may not be involutive or the coordinate system may not be δ -regular. However, we obtain the Janet tabular for g 2 :

{ v 33 = 0 v 22 = 0 [ 1 2 3 1 2 ]

and thus the Janet tabular for g 3 :

{ v 333 = 0 v 233 = 0 v 223 = 0 v 222 = 0 v 133 = 0 v 122 = 0 [ 1 2 3 1 2 1 2 1 2 1 1 ]

We let the reader check as an exercise that g 2 is not 2-acyclic by counting the dimensions in the long sequence:

0 g 4 δ T * g 3 δ 2 T * g 2 δ 3 T * T *

and that g 3 is involutive, thus 2-acyclic, with characters ( 0,0,4 ) as in the previous example. It follows that d i m ( g 2 ) = 6 2 = 4 , d i m ( g 3 ) = 0 + 0 + 4 = 4 , d i m ( g 4 ) = 4, . We obtain from the main theorem ρ r ( R 2 ( 1 ) ) = R 2 + r ( 1 ) . It is easy to chek that R 2 ( 1 ) = R 2 with d i m ( R 2 ) = 8 , d i m ( R 3 ) = 8 + 4 = 12 , d i m ( R 4 ) = 11 + 4 = 15 , d i m ( R 5 ) = 56 39 = 13 + 4 = 17 , . We have after two prolongation:

{ y 2233 = u 33 y 2233 + x 2 y 122 + 2 y 12 = v 22 x 2 y 122 = x 2 u 1

and thus 2 y 12 = u 33 v 22 x 2 u 1 . We may thus consider the new second order system R 2 = R 2 ( 2 ) R 2 with a strict inclusion and d i m ( R 2 ) = 7 :

y 33 x 2 y 1 = v , y 22 = u , y 12 = w , 2 w = u 33 v 22 x 2 u 1

We may start again with R 2 and study its symbol g 2 defined by the 3 linear equations with the following Janet tabular obtained after doing the same change of local coordinates as before:

{ v 33 = 0 v 22 = 0 v 12 = 0 [ 1 2 3 1 2 1 ]

This symbol is not involutive but its prolongation g 3 , defined by the 8 equations:

{ v 333 = 0 v 233 = 0 v 223 = 0 v 222 = 0 v 133 = 0 v 123 = 0 v 122 = 0 v 112 = 0 [ 1 2 3 1 2 1 2 1 2 1 1 1 1 ]

is involutive with characters ( 0,0,2 ) and we may consider again the system:

{ y 33 x 2 y 1 = v y 22 = u y 12 = w

with 2 w = u 33 v 22 x 2 u 1 . As before, instead of writing out the system R 3 and studying its formal inegrability by an elementary but tedious computation, we shall use a trick, knowing in advance that the generating CC must be of order at least 1 + 1 = 2 because g 2 had to get one prolongation in order to become involutive and thus 2-acyclic.

・ Step 2: It thus remains to find out the CC for ( u , v ) in the initial inhomogeneous system. As we have used two prolongations in order to exhibit R 2 , we have second order formal derivatives of u and v in the right members. Now, from the above argument, we have second order CC for the new right members and could hope therefore for a fourth order generating CC. The trick is to use the three different brackets of operators that can be obtained. We obtain in a formal way:

[ d 33 x 2 d 1 , d 22 ] y = y 2233 x 2 y 122 y 2233 + x 2 y 122 + 2 y 12 = 2 y 12 = u 33 x 2 u 1 v 22

Then:

2 [ d 2 , d 12 ] y = 0 = d 2 ( 2 d 12 y ) 2 d 1 ( d 22 y ) = ( u 233 v 222 x 2 u 12 u 1 ) 2 u 1

brings the third order CC:

A u 233 v 222 x 2 u 12 3 u 1 = 0

2 [ d 33 x 2 d 1 , d 12 ] y = 2 y 1233 2 x 2 y 112 2 y 1233 + 2 x 2 y 112 + 2 y 11 = 2 y 11 = 2 ( d 33 x 2 d 1 ) w 2 v 12 = u 3333 v 2233 2 x 2 u 133 + x 2 v 122 2 v 12 + ( x 2 ) 2 u 11

brings the new first order equation:

2 y 11 = 2 ω = u 3333 v 2233 2 x 2 u 133 + x 2 v 122 2 v 12 + ( x 2 ) 2 u 11 = 0

Accordingly, we may start afresh with the new system R 2 = R 2 ( 4 ) R 2 R 2 which is surprisingly of finite type with d i m ( R 2 ) = 10 4 = 6 , g 3 = 0 and defined by the 4 second order PD equations:

y 33 x 2 y 1 = v , y 22 = u , y 12 = w , y 11 = ω , w j 2 ( u , v ) , ω j 4 ( u , v )

We obtain therefore:

2 y 112 = u 23333 v 22233 2 x 2 u 1233 2 u 133 + x 2 v 1222 v 122 + ( x 2 ) 2 u 112 + 2 x 2 u 11 2 y 112 = u 133 + v 122 + x 2 u 11

and thus a CC of order 5, namely:

B u 23333 v 22233 2 x 2 u 1233 3 u 133 + x 2 v 1222 + ( x 2 ) 2 u 112 + 3 x 2 u 11 = 0

However, we have indeed the identity A 33 x 2 A 1 B = 0 and thus B j 2 ( A ) , that is B is a not a new generating fifth order CC as it is only a consequence of the prolongations of A. Using now the bracket of operators [ d 11 , d 33 ] = 0 that has not been already used, we get:

2 y 1133 = u 333333 v 223333 2 x 2 u 13333 + x 2 v 12233 2 v 1233 + ( x 2 ) 2 u 1133 2 y 1133 + 2 x 2 y 111 = 2 v 11 2 x 2 y 111 = x 2 u 13333 + x 2 v 12233 + 2 ( x 2 ) 2 u 1133 ( x 2 ) 2 v 1122 + 2 x 2 v 112 ( x 2 ) 3 u 111

We obtain therefore a new sixth order CC:

C u 333333 v 223333 3 x 2 u 13333 + 2 x 2 v 12233 2 v 1233 + 3 ( x 2 ) 2 u 1133 ( x 2 ) 2 v 1122 + 2 x 2 v 112 ( x 2 ) 3 u 111 2 v 11 = 0

which cannot be a differential consequence of A. After tedious computations, one can find the differential identity:

A 3333 2 x 2 A 133 + ( x 2 ) 2 A 11 C 2 = 0

The corresponding simplest free resolution, written with differential modules, is thus:

0 D 4 D 2 6 D 2 2 D p M 0

Again, the total number of generating CC, that is 1 + 1 = 2 , has nothing to do with the differential transcendence degree of the CC differential module which is still d i m ( F 0 ) d i m ( E ) = 2 1 = 1 because r k D ( M ) = 0 .

3. Mathematical Tools

Instead of starting with a linear system R q J q ( E ) of order q on E, let us start with a bundle map Φ : J q ( E ) F 0 with R q = k e r ( Φ ) and let us consider the linear PD operator D : E j q J q ( E ) Φ F 0 . Using the canonical inclusion J q + r ( E ) J r ( J q ( E ) ) , let us define the r- prolongation ρ r ( Φ ) : J q + r ( E ) J r ( J q ( E ) ) J r ( Φ ) J r ( F 0 ) . The general case of the successive prolongations with r 0 is described by the following commutative and exact diagram:

0 0 0 0 g q + r + 1 S q + r + 1 T * E σ r + 1 ( Φ ) S r + 1 T * F 0 h r + 1 0 0 R q + r + 1 J q + r + 1 ( E ) ρ r + 1 ( Φ ) J r + 1 ( F 0 ) Q r + 1 0 π q + r q + r + 1 π r r + 1 0 R q + r J q + r ( E ) ρ r ( Φ ) J r ( F 0 ) Q r 0 0 0 0

with symbol-map induced in the upper symbol sequence ( [19] [20] [21] ).

Chasing in this diagram while applying the “snake” lemma ( [10] [22] [23] ), we obtain the long exact connecting sequence:

0 g q + r + 1 R q + r + 1 R q + r h r + 1 Q r + 1 Q r 0

which is thus connecting in a tricky way FI (lower left) with CC (upper right). Needless to say that absolutely no classical procedure can produce such a result which is thus totally absent from the GR papers already quoted.

Setting H ( R q + r ) = R q + r / π q + r q + r + 1 ( R q + r + 1 ) , we have equivalently the shorter long exact sequence:

0 H ( R q + r ) h r + 1 Q r + 1 Q r 0

As a possible interpretation, d i m ( Q r ) is the total number of CC of order 0,1, up to r included. However, the problem to solve is to study the structure of the projective limit of vector bundles made by the induced epimorphisms Q r + 1 Q r . Of course, as it is mostly realized in the examples, we have to suppose that R q is sufficiently regular in such a way that the R q + r are vector bundles r 0 and that the R q + r ( s ) = π q + r q + r + s ( R q + r + s ) are also vector bundles, such a situation being in particular always realized when R q J q ( E ) or D are defined over a differential field K. In this case, introducing the filtered noetherian ring D = K [ d 1 , , d n ] = K [ d ] of differential operators with coefficients in K, we may introduce a differential module M with induced filtration 0 = M 0 M 1 M q m = M in such a way that the system R = R = h o m K ( M , K ) associated with M with R q = h o m K ( M q , K ) is of course automatically FI (care). Following Macaulay in ( [24] ), we have already proved in many places ( [3] [10] ) that R is a differential module for the Spencer operator d : R T * R : f d x i d i f with d i : R R : R q + 1 R q defined by the explicit formula:

( d i f ) μ k = i f μ k f μ + 1 i k K

It is important to notice that such an operator/system is far from being formally integrable because:

i ( j f μ k f μ + 1 j k ) j ( i f μ k f μ + 1 i k ) = j f μ + 1 i k i f μ + 1 j k

As can be seen from the examples previously presented, starting with Ψ r for a given r, the main problem is to compare the epimorphism Ψ r + 1 : J r + 1 ( F 0 ) Q r + 1 with the morphism ρ 1 ( Ψ r ) : J r + 1 ( F 0 ) J 1 ( Q r ) in the following commutative diagram which may not be exact:

0 0 0 0 0 g q + r + 1 S q + r + 1 T * E σ r + 1 ( Φ ) S r + 1 T * F 0 σ 1 ( Ψ r ) T * Q r 0 R q + r + 1 J q + r + 1 ( E ) ρ r + 1 ( Φ ) J r + 1 ( F 0 ) ρ 1 ( Ψ r ) J 1 ( Q r ) π q + r q + r + 1 π r r + 1 0 R q + r J q + r ( E ) ρ r ( Φ ) J r ( F 0 ) Ψ r Q r 0 0 0 0

where the central row is induced from the long exact sequence:

0 J 1 ( R q + r ) J 1 ( J q + r ( E ) ) J 1 ( J r ( F 0 ) ) J 1 ( Q r ) 0

and may not be exact.

PROPOSITION 3.1: We have only in general:

B r + 1 = i m ( ρ r + 1 ( Φ ) ) = k e r ( Ψ r + 1 ) k e r ( ρ 1 ( Ψ r ) ) = ρ 1 (Br)

Proof: Denoting the Spencer operator by d in place of the standard notation D of the literature that could be confused with the ring D of differential operators, we have the following commutative diagram:

0 R q + r + 1 J q + r + 1 ( E ) ρ r + 1 ( Φ ) B r + 1 J r + 1 ( F 0 ) d d d d 0 T * R q + r T * J q + r ( E ) ρ r ( Φ ) T * B r T * J r (F0)

As B r + 1 projects onto B r and d i B r + 1 B r , it follows from ( [12] , Propositions 10, p 83) or ( [10] , Remark 2.9, p 315) that B r + 1 ρ 1 ( B r ) . We have thus a projective limit of systems, each one being defined by more equations than the preceding one and such a procedure must finish with a FI system that can even be prolonged, as we shall see in the examples, in order to obtain an involutive system that may be used to start a Janet sequence. The decision to stop is provided by the maximum order of the CC obtained, namely of order bounded by r + s + 1 = t if the system R q + r ( s ) is involutive or at least with a 2-acyclic symbol. ,

The idea is to use the composite morphism S r T * F 0 Q r while chasing in order to prove that any element of J r + 1 ( F 0 ) killed by ρ 1 ( Ψ r ) can be decomposed into the sum of an element in i m ( ρ r ( Φ ) ) plus an element in S r + 1 T * F 0 killed by σ 1 ( Ψ r ) . With more details, setting for simplicity S r T * F 0 = S r ( F 0 ) , introducing the coboundary bundle B ( S r + 1 ( F 0 ) ) = i m ( σ r + 1 ( Φ ) ) and the cocycle bundle Z ( S r + 1 ( F 0 ) ) = k e r ( σ 1 ( Ψ r ) ) , we may define the corresponding cohomology bundle H ( S r + 1 ( F 0 ) ) = Z ( S r + 1 ( F 0 ) ) / B ( S r + 1 ( F 0 ) ) . We may also define similarly H ( J r + 1 ( F 0 ) ) = k e r ( ρ 1 ( Ψ r ) ) / i m ( ρ r + 1 ( Φ ) ) and we obtain the following crucial proposition (See [4] , Example 2.A.9) through a chase left to the reader as an exercise:

PROPOSITION 3.2: There exists a short exact sequence:

0 H ( R q + r ) H ( S r + 1 ( F 0 ) ) H ( J r + 1 ( F 0 ) ) 0

Let us now deal with the symbol cohomology by chasing in the following commutative diagram:

0 0 0 0 0 g q + r + 1 S q + r + 1 T * E S r + 1 T * F 0 T * Q r δ δ δ 0 T * g q + r T * S q + r T * E T * S r T * F 0 T * Q r δ δ δ 0 2 T * g q + r 1 2 T * S q + r 1 T * E 2 T * S r 1 T * F 0 0 δ δ 0 3 T * g q + r 2 = 3 T * S q + r 2 T * E

where neither the first nor the second upper columns may be exact and where the left column may not be exact, unless gq is involutive or 2-acyclic. Chasing with the same notations, we obtain:

PROPOSITION 3.3: There exists an exact sequence:

0 H 2 ( g q + r 1 ) H ( S r + 1 ( F 0 ) ) T * H ( S r ( F 0 ) )

The upper left arrows are not in general epimorphisms and it may be sometimes useful to consider h r as a kind of symbol in the more abstract diagram:

0 0 0 0 g q + r + 1 S q + r + 1 T * E S r + 1 T * F 0 h r + 1 0 δ δ δ δ 0 T * g q + r T * S q + r T * E T * S r T * F 0 T * h r 0 δ δ δ δ 0 2 T * g q + r 1 2 T * S q + r 1 T * E 2 T * S r 1 T * F 0 2 T * h r 1 0 δ δ δ 0 3 T * g q + r 2 3 T * S q + r 2 T * E 3 T * S r 2 T * F 0 δ δ 0 4 T * g q + r 3 4 T * S q + r 3 T * E

where the rows are now exact. However, understanding the meaning of h r as a kind of new symbol may not be possible unless h r + 1 T * h r is a monomorphism, that is when g q is 2-acyclic and h r is 1-acyclic, that is when g q is also 3-acyclic (or involutive). Once more, we understand the crucial importance of 2-acyclicity but we recall that the only symbol known to be 2-acyclic without being involutive is the symbol of the conformal Killing system whenever n 4 , which is also 3-acyclic whenever n 5 ( [3] [11] [12] [13] ).

4. Applications

MACAULAY EXAMPLE REVISITED:

With m = 1 , n = 3 , q = 2 , K = , let us introduce two operators P , Q D = K [ d 1 , d 2 , d 3 ] and consider the second order system R 2 J 2 ( E ) used by Macaulay as early as in 1916 ( [4] [24] ):

Q y y 33 = v , P y y 13 y 2 = u

We have the strict inclusions:

R 2 ( 2 ) R 2 ( 1 ) R 2 J 2 ( E ) 6 < 7 < 8 < 10

As g 2 ( 2 ) g 2 ( 1 ) g 2 are involutive, we obtain ρ r ( R 2 ( 2 ) ) = R r + 2 ( 2 ) by using the Prolongation/Projection (PP) procedure. We exhibit the parametric jets of the bundles that will be used in the following diagrams:

p a r 2 = { y , y 1 , y 2 , y 3 , y 11 , y 12 , y 22 , y 23 }

p a r 3 = { y , y 1 , y 2 , y 3 , y 11 , y 12 , y 22 , y 111 , y 112 , y 122 , y 222 , y 223 }

p a r 4 = { y , y 1 , y 2 , y 3 , y 11 , y 12 , y 111 , y 112 , y 122 , y 222 , y 1111 , y 1112 , y 1122 , y 1222 , y 2222 , y 2223 }

and thus d i m ( R 2 ) = 8, d i m ( R 3 ) = 12, d i m ( R 4 ) = 16 . More generally, we let the reader prove that d i m ( R r + 2 ) = 4 r + 8, r 0 , thus d i m ( R r + 4 ) = 4 r + 16, r 0 and d i m ( g r + 4 ) = r + 6, r 0 .

We have the Janet tabular for R 2 ( 2 ) :

{ y 33 = v y 23 = v 1 u 3 y 22 = v 11 u 13 u 2 y 13 y 2 = u [ 1 2 3 1 2 1 2 1 ]

The two CC are:

A v 13 v 2 u 33 = 0 , A 1 v 113 v 12 u 133 = 0

while the other ones are what we called identity to zero like:

d 2 ( y 13 y 2 ) d 1 ( y 23 ) + d 2 y 2 = 0 = u 2 ( v 11 u 13 ) + ( v 11 u 13 u 2 ) = 0

There is thus only one generating CC of order 2, namely A = 0 , given by the commutation relation P Q Q P 0 and the corresponding operator D 1 is thus surely formally subjective. Setting F 1 = Q 2 , we obtain the following diagram with exact central and lower rows whenever r 1 .

0 0 0 0 0 g r + 4 S r + 4 T * E S r + 2 T * F 0 S r T * F 1 0 0 R r + 4 J r + 4 ( E ) J r + 2 ( F 0 ) J r ( F 1 ) 0 0 R r + 3 J r + 3 ( E ) J r + 1 ( F 0 ) J r 1 ( F 1 ) 0 0 0 0

0 0 0 0 0 r + 6 ( r + 5 ) ( r + 6 ) / 2 ( r + 3 ) ( r + 4 ) ( r + 1 ) ( r + 2 ) / 2 0 0 4 r + 16 ( r + 5 ) ( r + 6 ) ( r + 7 ) / 6 ( r + 3 ) ( r + 4 ) ( r + 5 ) / 3 ( r + 1 ) ( r + 2 ) ( r + 3 ) / 6 0 0 4 r + 12 ( r + 4 ) ( r + 5 ) ( r + 6 ) / 6 ( r + 2 ) ( r + 3 ) ( r + 4 ) / 3 r ( r + 1 ) ( r + 2 ) / 6 0 r + 2 0 0 0 0

0 0 0 0 0 g r + 4 S r + 4 T * E S r + 2 T * F 0 S r T * F 1 0 δ δ δ δ 0 T * g r + 3 T * S r + 3 T * E T * S r + 1 T * F 0 T * S r 1 T * F 1 0 δ δ δ 0 2 T * g r + 2 2 T * S r + 2 T * E 2 T * S r T * F 0 δ δ 0 3 T * g r + 1 = 3 T * S r + 1 T * E 0 0

where S 4 T * E S 4 T * and F 1 Q 2 while Q 1 = 0 as there is no CC of order 1. From the snake lemma and a chase, we obtain the long exact connecting sequence when r = 0 :

0 g 4 R 4 R 3 h 2 F 1 0

0 6 16 12 3 1 0

relating FI (lower left) to CC (upper right). By composing the epimorphism S 2 T * F 0 h 2 with the epimorphism h 2 F 1 , we obtain an epimorphism S 2 T * F 0 F 1 and the long exact sequence:

0 g 4 S 4 T * E S 2 T * F 0 F 1 0

which is nevertheless not a long ker/coker exact sequence by counting the dimensions as we have 6 15 + 12 1 = 2 0 .

The above diagrams illustrate perfectly the three propositions of Section 2. We have in particular:

H ( J r + 2 ( F 0 ) ) = 0 H ( J r + 1 ( F 0 ) ) = 0 , H ( S r + 2 ( F 0 ) ) = H ( R r + 3 ) 0

and the formally exact sequence, which is nevertheless not strictly exact though 1 2 + 1 = 0 :

0 Θ E 2 D F 0 2 D 1 F 1 0

0 Θ 1 2 1 0

We remind the reader that, contrary to the situation met with FI systems where the exactness on the jet level is obtained inductively from the exactness on the symbol level, here we discover that we may have the exactness on the jet level without having exactness on the symbol level.

EXAMPLE 2.1 REVISITED:

First of all, let us compute the dimensions and the parametric jets that will be used in the following diagrams.

p a r 1 = p a r 2 = { y , y 1 , y 2 , y 3 , y 1 2 } ,

n = d i m ( X ) = 3, m = d i m ( E ) = 2, d i m ( R 1 ) = d i m ( R 2 ) = 3, d i m ( g 1 ) = 1, g 2 = 0 g 3 = 0

0 0 0 0 S 3 T * E S 2 T * F 0 T * Q 1 0 0 R 3 J 3 ( E ) J 2 ( F 0 ) J 1 ( Q 1 ) 0 π 4 3 π 1 2 0 R 2 J 2 ( E ) J 1 ( F 0 ) Q 1 0 0 0 0 0

0 0 0 0 20 30 9 0 0 3 40 50 12 0 π 2 3 π 1 2 0 3 20 20 3 0 0 0 0 0

0 0 0 0 S 3 T * E S 2 T * F 0 T * Q 1 0 δ δ 0 T * S 2 T * E T * T * F 0 T * Q 1 0 δ δ 0 2 T * g 1 2 T * T * E 2 T * F 0 0 δ δ 0 3 T * E = 3 T * E 0 0 0

0 0 0 0 20 30 9 0 δ δ 0 36 45 9 0 δ δ 0 3 18 15 0 δ δ 0 2 = 2 0 0 0

It is not at all evident to study these diagrams. We have d i m ( B ( S 2 ( F 0 ) ) ) = 20 < d i m ( Z ( S 2 ( F 0 ) ) ) = 30 9 = 21 . We have already proved that d i m ( H ( S 2 ( F 0 ) ) ) = 21 20 = 1 = d i m ( H 2 ( g 1 ) ) = 3 2 = 1 , a result not evident at first sight explaining why the only second order additional generating CC is nothing else than the Riemann tensor in dimension equal to 2.

We have explained in ( [4] ) that such a system has its origin in the study of the integration of the Killing system for the Schwarzschild metric, which is not FI. With more details, let us use the Boyer-Lindquist coordinates ( t , r , θ , ϕ ) = ( x 0 , x 1 , x 2 , x 3 ) instead of the Cartesian coordinates ( t , x , y , z ) and consider the Schwarzschild metric ω = A ( r ) d t 2 ( 1 / A ( r ) ) d r 2 r 2 d θ 2 r 2 sin 2 ( θ ) d ϕ 2 and ξ = ξ i d i T , let us introduce ξ i = ω r i ξ r with the 4 formal derivatives ( d 0 = d t , d 1 = d r , d 2 = d θ , d 3 = d ϕ ) . With speed of light c = 1 and A = 1 m r where m is a constant, the metric can be written in the diagonal form:

( A 0 0 0 0 1 / A 0 0 0 0 r 2 0 0 0 0 r 2 sin 2 (θ))

Using the notations that can be found in the theory of differential modules, let us consider the Killing equations:

Ω L ( ξ ) ω = 0 Ω i j d i ξ j + d j ξ i 2 γ i j r ξ r = 0

where we have introduced the Christoffel symbols γ while setting A = r A in the differential field K of coefficients. As in the previous Macaulay example and in order to avoid any further confusion between sections and derivatives, we shall use the sectional point of view and rewrite the previous equations in the symbolic form L ( ξ 1 ) ω = Ω S 2 T * where L is the formal Lie derivative:

{ ξ 3 , 3 + sin ( θ ) cos ( θ ) ξ 2 + r A sin 2 ( θ ) ξ 1 = 1 2 Ω 33 ξ 2 , 3 + ξ 3 , 2 2 cot ( θ ) ξ 3 = Ω 23 ξ 1 , 3 + ξ 3 , 1 2 r ξ 3 = Ω 13 ξ 0 , 3 + ξ 3 , 0 = Ω 03 ξ 2 , 2 + r A ξ 1 = 1 2 Ω 22 ξ 1 , 2 + ξ 2 , 1 2 r ξ 2 = Ω 12 ξ 0 , 2 + ξ 2 , 0 = Ω 02 ξ 1 , 1 + A 2 A ξ 1 = 1 2 Ω 11 ξ 0 , 1 + ξ 1 , 0 A A ξ 0 = Ω 01 ξ 0 , 0 A A 2 ξ 1 = 1 2 Ω 00

This system R 1 J 1 ( T ) is far from being involutive because it is finite type with second symbol g 2 = 0 defined by the 40 equations v i j k = 0 in the initial coordinates. From the symmetry, it is clear that such a system has at least 4 solutions, namely the time translation t ξ 0 = 1 ξ 0 = A and, using Cartesian coordinates ( t , x , y , z ) , the 3 space rotations y z z y , z x x z , x y y x .

These results are brought by the formal Lie derivative of the Weyl tensor because the Ricci tensor vanishes by assumption and we have the splitting R i e m a n n R i c c i W e y l according to the fundamental diagram II that we discovered as early as in 1988 ( [25] ), still not acknowledged though it can be found in ( [1] [2] [3] ). In particular, as the Ricci part is vanishing by assumption, we may identify the Riemann part with the Weyl splitting part as tensors ( [3] ) and it is possible to prove (using a tedious direct computation or computer algebra) that there are only 6 non-zero components. It is important to notice that this result, bringing a strong condition on the zero jets because of the Lie derivative of the Weyl tensor and thus on the first jets, involves indeed the first derivative of the Weyl tensor because we have a term in ( A ) . When Ω = 0 , we obtain after 2 prolongations the additional 5 new first order PD equations:

ξ 1 = 0 , ξ 1 , 2 = 0 , ξ 1 , 3 = 0 , ξ 0 , 2 = 0 , ξ 0 , 3 = 0

As we are dealing with sections, ξ 1 = 0 does imply ξ 1 , 1 = 0 and ξ 0 , 0 = 0 but does not imply ξ 1 , 0 = 0 , these later condition being only brought by one additional prolongation and we have the strict inclusions R 1 ( 3 ) R 1 ( 2 ) R 1 ( 1 ) = R 1 . Hence, it remains to determine the dimensions of the subsystems R 1 = R 1 ( 2 ) and R 1 = R 1 ( 3 ) with the strict inclusion R 1 R 1 , exactly again like in the Macaulay example. Knowing that d i m ( R 1 ) = d i m ( R 2 ) = 10 , d i m ( R 3 ) = 5 , d i m ( R 4 ) = 4 , we have thus obtained the 15 equations defining R 1 with d i m ( R 1 ) = 20 15 = 5 and the 16 equations defining R 1 with d i m ( R 1 ) = 20 16 = 4 , namely:

ξ 3 , 3 + sin ( θ ) cos ( θ ) ξ 2 = 0 ξ 2 , 3 + ξ 3 , 2 2 cot ( θ ) ξ 3 = 0 ξ 1 , 3 = 0 ξ 0 , 3 = 0 ξ 2 , 2 = 0

ξ 1 , 2 = 0 ξ 0 , 2 = 0 ξ 3 , 1 2 r ξ 3 = 0 ξ 2 , 1 2 r ξ 2 = 0 ξ 1 , 1 = 0

ξ 0 , 1 A A ξ 0 = 0 ξ 3 , 0 = 0 ξ 2 , 0 = 0

ξ 0 , 1 A A ξ 0 = 0 ξ 3 , 0 = 0 ξ 2 , 0 = 0

Setting now in an intrinsic way ξ 0 = A ξ 0 , ξ 1 = 1 A ξ 1 , ξ 2 = r 2 ξ 2 and in a non-intrinsic way (care) ξ 3 = r 2 ξ 3 , we may even simplify these equations and get a system not depending on A anymore:

{ ξ 3 3 + sin ( θ ) cos ( θ ) ξ 2 = 0 ξ 3 2 + ξ 2 3 2 cot ( θ ) ξ 3 = 0 ξ 3 1 = 0 ξ 3 0 = 0 ξ 1 3 = 0 ξ 1 2 = 0 ξ 1 1 = 0 ξ 1 0 = 0 ξ 0 3 = 0 ξ 0 2 = 0 ξ 0 1 = 0 ξ 0 0 = 0 ξ 2 2 = 0 ξ 2 1 = 0 ξ 2 0 = 0 ξ 1 = 0

It is easy to check that R 1 ( 3 ) , having minimum dimension equal to 4, is formally integrable, though not involutive as it is finite type, and to exhibit 4 solutions linearly independent over the constants. Indeed, we must have ξ 0 = c where c is a constant and we may drop the time variable not appearing elsewhere while using the equation ξ 1 = 0 . It follows that ξ 2 = f ( θ , ϕ ) , ξ 3 = g ( θ , ϕ ) while f , g are solutions of the first, second and fifth equations of Killing type wih a general solution depending on 3 constants, a result leading to an elementary probem of 2-dimensional elasticity left to the reader as an exercise. The system R 1 ( 3 ) is formally integrable while the system R 2 ( 2 ) is involutive. Having in mind the PP procedure, it follows that the CC could be of order 2, 3 and even 4. Equivalently, we may cut the integration of this system into three systems:

1) First of all, we have ξ 1 = 0 and thus ξ 0 1 = 0 , ξ 1 1 = 0 , ξ 2 1 = 0 , ξ 3 1 = 0 .

2) Then, we may consider ξ 0 0 = 0 , ξ 1 0 = 0 , ξ 2 0 = 0 , ξ 3 0 = 0 ξ 0 = c .

3) Finally, we arrive to the FI system with the same properties as the ones found for Example 2.1:

{ ξ 3 3 + sin ( θ ) cos ( θ ) ξ 2 = 0 ξ 3 2 + ξ 2 3 2 cot ( θ ) ξ 3 = 0 ξ 1 3 = 0 ξ 1 2 = 0 ξ 2 2 = 0

that is with 3 generating first order CC and 1 additional second order generating CC.

Proceeding like in the motivating examples, we may introduce the inhomogeneous systems:

{ ξ 1 = U , ξ 1 , 2 = V 2 , ξ 1 , 3 = V 3 , ξ 0 , 2 = W 2 , ξ 0 , 3 = W 3 } j 2 (Ω)

and we finally obtain 16 PD equations, namely ξ 1 = U plus the 15 PD equations:

ξ 0 , 0 = 1 2 Ω 00 + A A 2 U , ξ 0 , 1 A A ξ 0 = Ω 01 d 0 U j 3 ( Ω ) , ξ 0 , 2 = W 2 , ξ 0 , 3 = W 3

ξ 1 , 0 = d 0 U j 3 ( Ω ) , ξ 1 , 1 = 1 2 Ω 11 A 2 A U , ξ 1 , 2 = V 2 , ξ 1 , 3 = V 3

ξ 2 , 0 = Ω 02 W 2 , ξ 2 , 1 2 r ξ 2 = Ω 12 V 2 , ξ 2 , 2 = 1 2 Ω 22 r A U

ξ 3 , 0 = Ω 03 W 3 , ξ 3 , 1 2 r ξ 3 = Ω 13 V 3 , ξ 3 , 2 + ξ 2 , 3 2 cot ( θ ) ξ 3 = Ω 23 ,

ξ 3 , 3 + sin ( θ ) cos ( θ ) ξ 2 = 1 2 Ω 33 r A sin 2 ( θ ) U

As a byproduct, we have d i m ( R 3 ) = 5, d i m ( R 4 ) = 4 and we obtain 15 second order CC in j 2 ( Ω ) along the ker/coker exact sequence:

0 R 3 J 3 ( T ) J 2 ( S 2 T * ) Q 2 0 B 2 0 0

0 5 140 150 15 0 135

Then, we have identities to zero like d 0 ξ 1 ξ 0 1 = 0 but we have also surely the three third order CC like d 1 ξ 1 ξ 1 1 = 0 , d 2 ξ 1 ξ 2 1 , d 3 ξ 1 ξ 3 1 , then perhaps the other third order CC d 2 ξ 3 0 d 3 ξ 2 0 = 0 and perhaps even fourth order CC like d 0 ξ 1 0 d 1 ξ 0 0 = 0 which is containing the leading term d 00 U after substitution. However, we have the linearization formulas:

ρ k l , i j = ω k r ρ l , i j r R k l , i j = ω k r R l , i j r + ρ l , i j r Ω k r

ρ i j = ρ i , r j r ω r s R r i , s j = R i j + ω r s ρ i , r j t Ω s t R i j

and obtain therefore the formulas:

2 r 3 m sin 2 ( θ ) R 23 , 30 + 2 3 Ω 02 = ξ 0 , 2 = 2 r 3 A 3 m R 01 , 12 + 1 3 Ω 02 R 02 = R 0 , 12 1 + R 3 , 32 3 = 0 2 r 3 m R 23 , 02 + 2 3 Ω 03 = ξ 0 , 3 = 2 r 3 A 3 m R 01 , 13 + 1 3 Ω 03 R 03 = R 0 , 13 1 + R 0 , 23 2 = 0

with two similar ones for ξ 1,2 and ξ 1,3 showing the unexpected partition of the Ricci tensor:

{ R i j } = { R 00 , R 11 , R 22 , R 33 } + { R 12 , R 13 , R 02 , R 03 } + { R 01 , R 23 }

determined by the 15 = 10 + 5 = 4 + 4 + 2 second order CC that we have exhibited.

Now, after one prolongation, we get:

ξ 1 , 00 + ( A A 2 ( A ) 2 4 ) ξ 1 + 1 2 d 1 Ω 00 d 0 Ω 01 A 2 A Ω 00 + A A 24 Ω 11 = 0

and thus d 1 ξ 0 , 0 d 0 ξ 0 , 1 = 0 . Similarly, we have:

ξ 1 , 01 + A 2 A ξ 1 , 0 1 2 d 0 Ω 11 d 0 ( d 1 ξ 1 + A 2 A ξ 1 Ω 11 ) = 0

and thus d 1 ξ 1 , 0 d 0 ξ 1 , 1 = 0 . It follows that d 1 U + A 2 A 1 2 Ω 11 = 0 is a generating CC of order 3 but d 01 U A 2 A d 0 U d 0 Ω 11 = 0 is not a generating CC of order 4.

In order to proceed further on, we notice that the generating CC of order 3 already found can be written as:

d 1 U + A 2 A U 1 2 Ω 11 = 0 , d 2 U V 2 = 0 , d 3 U V 3 = 0

Using crossed derivatives, we get:

d 1 V 2 + A 2 A d 2 U 1 2 d 2 Ω 11 = 0 , d 1 V 3 + A 2 A d 3 U 1 2 d 3 Ω 11 = 0 , d 2 V 3 d 3 V 2 = 0

and thus d 1 ξ 1 , 2 d 2 ξ 1 , 1 = 0 , d 1 ξ 1 , 3 d 3 ξ 1 , 1 = 0 , d 2 ξ 1 , 3 d 3 ξ 1 , 2 = 0 .

However, in order to prove that d 2 ξ 0 , 3 d 3 ξ 0 , 2 = 0 or equivalently that d 2 W 3 d 3 W 2 = 0 , the previous procedure cannot work but we must never forget that U , V 2 , V 3 , W 2 , W 3 both belong to j 2 ( Ω ) . Introducing the formal Lie derivative R = L ( ξ 1 ) ρ , we recall that:

W 2 = 2 r 3 A 3 m R 01 , 12 + 1 3 Ω 02 , W 3 = 2 r 3 A 3 m R 01 , 13 + 1 3 Ω 03

Hence, linearizing the Bianchi identity:

1 ρ 01 , 23 + 2 ρ 01 , 31 + 3 ρ 01 , 12 = 0 d 1 R 01 , 23 + d 2 R 01 , 31 + d 3 R 01 , 12 + = 0

we have proved in ( [4] ) that the third order CC d 2 W 3 d 3 W 2 = 0 is not a generating one because it is just a differential consequence of the second order CC R 01 , 23 = 0 .

Finally, as already noticed, the symbol g 1 g 1 T * T is not involutive and even 2-acyclic because otherwise there should only be first order CC for the right members defining the system R 1 J 1 ( T ) . As a byproduct, we have, at least on the symbol level, the second order CC:

d 22 ξ 3 , 3 + d 33 ξ 2 , 2 d 23 ( ξ 3 , 2 + ξ 2 , 3 ) = 0

and thus:

d 22 ( 1 2 Ω 33 r A sin 2 ( θ ) ξ 1 sin ( θ ) cos ( θ ) ξ 2 ) + d 33 ( 1 2 Ω 22 r A ξ 1 ) d 23 ( Ω 23 + 2 cot ( θ ) ξ 3 ) = 0

containing surely d 22 ξ 1 , d 22 ξ 2 , d 33 ξ 1 , d 23 ξ 3 and thus surely d 2 U , d 2 V 2 , d 3 V 3 , producing therefore a third order CC that cannot be reduced by means of any Bianchi identity, that is we finally have 15 generating second order CC and 4 new generating third order CC, in a manner absolutely similar to that of all the motivating examples of this paper.

As shown in ( [4] ), the study of the Killing system for the Kerr metric is even more difficult because the space of solutions is reduced from 4 already given to the 2 infinitesimal generators { t , ϕ } only. Accordingly, we discover that the Schwarzschild and the Kerr metrics do behave quite differently and there is thus no hope at all for selecting specific solutions of the Einstein equations in vacuum. We consider this result as a key challenge when questioning the origin and existence of gravitational waves in general relativity and believe this problem has never been pointed out clearly for the very simple reason that the underlying mathematics are not known by physicists.

EXAMPLE 2.2 REVISITED:

Coming back to the system R 2 = R 2 ( 2 ) R 2 with a strict inclusion and second members ( u , v , w = u 33 v 12 x 2 u 1 ) , let us exchange x 1 with x 2 in order to have an involutive third order symbol g 3 in δ -regular coordinates and consider the system R 3 with now w = u 33 v 12 x 1 u 2 in the new coordinates:

{ y 333 x 1 y 23 = v 3 y 233 x 1 y 22 = v 2 y 223 = w 3 y 222 = w 2 y 133 x 1 y 12 = v 1 y 123 = u 3 y 122 = u 2 y 112 = u 1 y 33 x 1 y 2 = v y 22 = w y 12 = u 1 2 3 1 2 1 2 1 2 1 1 1 1

This new system is easily seen to be involutive and we have 3 + 8 + 9 = 20 first order CC if we consider the second members as just simple notations. Substituting and taking now into account that we have in fact u 2 = d 2 u formally and so on, all these CC reduce to identities to zero of the form 0 = 0 but, using again the original coordinates, A w 2 u 1 = 0 , A 2 = 0 , A 3 = 0 , B w 33 v 11 x 2 w 1 = 0 , a system which is not FI. Accordingly, the generating CC are described by A of order 3 and B of order 4 with A 33 x 2 A 1 B 2 = 0 . I remain to check that this result is coherent with the diagrams of the previous section.

For this, we let the reader compute by hands or with computer algebra the following dimensions d i m ( R 2 ) = 1 + 3 + 4 = 8 , d i m ( R 3 ) = d i m ( J 3 ( E ) ) d i m ( J 1 ( F 0 ) ) = 12 because there is no CC of order 1, d i m ( R 4 ) = 15 d i m ( Q 2 ) = d i m ( R 4 ) d i m ( J 4 ( E ) ) + d i m ( J 2 ( F 0 ) ) = 15 35 + 20 = 0 because there is no CC of order 2, d i m ( R 5 ) = 17 d i m ( Q 3 ) = d i m ( R 5 ) d i m ( J 5 ( E ) ) + d i m ( J 3 ( F 0 ) ) = 17 56 + 40 = 1 because there is only 1 CC of order 3. We have therefore the long sequence:

0 R 6 J 6 ( E ) J 4 ( F 0 ) J 1 ( Q 2 ) 0

0 19 84 70 4 0

and obtain d i m ( H ( J 4 ( F 0 ) ) ) = ( 70 4 ) ( 84 19 ) = 66 65 = 1 both with d i m ( Q 4 ) = 5 , in a coherent way with the only CC A of order 3. We let the reader prove that we have similarly d i m ( Q 5 ) = 13 by taking into account the fact that B 2 = A 33 x 2 A 1 . In order to take into account the existence of a new generating CC of order 4, we let the reader check that d i m ( Q 4 ) = 5 and set F 1 = Q 4 in order to define a fourth order operator D 1 : F 0 F 1 by the involutive system:

{ B u 3333 + = 0 A 3 u 2333 + = 0 A 2 u 2233 + = 0 A 1 u 1233 + = 0 A u 233 + = 0 1 2 3 1 2 1 2 1

Starting anew from this operator, we obtain the first order involutive system:

{ d 3 A 3 x 2 A 1 B 2 = 0 d 3 A 2 d 2 A 3 = 0 d 3 A 1 d 1 A 3 = 0 d 3 A A 3 = 0 d 2 A 1 d 1 A 2 = 0 d 2 A A 2 = 0 d 1 A A 1 = 0 1 2 3 1 2 3 1 2 3 1 2 3 1 2 1 2 1 1 2 3 1 2 3 1 2 3 1 2 1 2 3

where each Janet tabular is induced from the preceding one till the end of the procedure as in ( [12] , p 153, 154 for details). We also notice that this system brings automatically the Spencer operator.

We obtain therefore the following differential sequence:

0 Θ E 2 D F 0 4 D 1 F 1 1 D 2 F 2 1 D 3 F 3 1 D 4 F 4 0

0 Θ 1 2 2 4 5 1 7 1 4 1 1 0

which is formally exact on the jet level, even if it is not strictly exact because the first operator is not FI, and we check that 1 2 + 5 7 + 4 1 = 0 . We notice that the part between F 0 and F 4 is typically a Janet sequence for D 1 .

It follows that we have the following long exact sequence on the level of jets, r 5 :

0 R r + 9 J r + 9 ( E ) J r + 7 ( F 0 ) J r + 3 ( F 1 ) J r + 2 ( F 2 ) J r + 1 ( F 3 ) J r ( F 4 ) 0

a result leading to:

d i m ( R r + 9 ) = 1 ( r + 10 ) ( r + 11 ) ( r + 12 ) / 6 2 ( r + 8 ) ( r + 9 ) ( r + 10 ) / 6 + 5 ( r + 4 ) ( r + 5 ) ( r + 6 ) / 6 7 ( r + 3 ) ( r + 4 ) ( r + 5 ) / 6 + 4 ( r + 2 ) ( r + 3 ) ( r + 4 ) / 6 1 ( r + 1 ) ( r + 2 ) ( r + 3 ) / 6 = 2 r + 25

and thus to d i m ( R r + 4 ) = 2 r + 15, r 0 , a result not evident to grasp at first sight because it comes from the lack of formal integrability of R 2 and the strict inclusion R 2 ( 2 ) R 2 .

EXAMPLE 2.3 REVISITED:

Coming back to the systems R 2 with second members ( u , v ) and R 2 = R 2 ( 2 ) R 2 with a strict inclusion and second members ( u , v , 2 w = u 33 v 22 x 2 u 1 ) , let us exchange x 1 with x 2 in order to have an involutive third order symbol g 3 in δ -regular coordinates but the system R 3 , with now w = u 33 v 12 x 1 u 2 in the new coordinates, is not FI. Hence, we must start anew with the system R 2 = R 2 ( 4 ) R 2 with a strict inclusion, described by the 4 PD equations:

{ y 33 x 2 y 1 = v y 22 = u y 12 = w y 11 = ω

where 2 ω = u 3333 v 2233 2 x 2 u 133 + x 2 v 122 2 v 12 + ( x 2 ) 2 u 11 . Using one prolongation, we get the third order PD equations:

{ y 333 x 2 y 13 = v 3 y 233 y 1 = v 2 + x 2 w y 223 = u 3 y 222 = u 2 y 133 = v 1 + x 2 ω y 123 = w 3 y 122 = u 1 = w 2 A y 113 = ω 3 y 112 = w 1 = ω 2 B y 111 = ω 1 1 2 3 1 2 1 2 1 2 1 1 1 1 1 1

and we discover that the symbol g 2 is finite type because g 3 = 0 . As we had to use one prolongation in order to get a 2-acyclic symbol, we obtain sixth order CC ( A , B , C ) j 2 ( u , v , w , ω ) . We refer the reader to ( [10] , p 315) or ( [13] , p 83) for more details on this delicate result.

Using the notations of the last section, we now provide the systems B 2 , B 3 , B 4 , B 5 , B 6 together and we notice the following striking results:

B 2 = J 2 ( F 0 ) , B 3 ρ 1 ( B 2 ) J 3 ( F 0 ) , B 4 = ρ 1 ( B 3 ) J 4 ( F 0 ) ,

B 5 = ρ 1 ( B 4 ) J 5 ( F 0 ) , B 6 ρ 1 ( B 5 ) J 6 (F0)

reaching therefore the following involutive system of order 6 where we did not quote B = 0 because we already proved that B = A 33 x 2 A 1 :

{ C u 333333 + = 0 A 333 u 233333 + = 0 A 233 u 223333 + = 0 A 223 u 222333 + = 0 A 222 u 222233 + = 0 A 133 u 123333 + = 0 A 123 u 122333 + = 0 A 122 u 122233 + = 0 A 113 u 112333 + = 0 A 112 u 112233 + = 0 A 111 u 111233 + = 0 A 33 u 23333 + = 0 A 23 u 22333 + = 0 A 22 u 22233 + = 0 A 13 u 12333 + = 0 A 12 u 12233 + = 0 A 11 u 11233 + = 0 A 3 u 2333 + = 0 A 2 u 2233 + = 0 A 1 u 1233 + = 0 A u 233 + = 0 1 2 3 1 2 1 2 1 2 1 2 1 1 1 1 1 1

Starting anew from this operator providing 46 CC, we obtain the first order involutive system:

{ d 3 A 333 2 x 2 A 133 + ( x 2 ) 2 A 11 C 2 = 0 d 3 A 233 d 2 A 333 = 0 d 3 A 223 d 2 A 233 = 0 d 3 A 222 d 2 A 223 = 0 d 3 A 1 d 1 A 3 = 0 d 3 A A 3 = 0 d 2 A 133 d 1 A 233 = 0 d 2 A 1 d 1 A 2 = 0 d 2 A A 2 = 0 d 1 A 33 A 133 = 0 d 1 A 1 A 11 = 0 d 1 A A 1 = 0 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 1 2 1 2 1 1 1

with 20 equations of class 3, 16 equations of class 2 and 10 equations of class 1. There are 36 CC providing an involutive system with 26 equations of class 3, 10 equations of classs 2 but no equation of class 1. We get a final system of 10 CC of class 3 without any CC. Like in the preceding application, we have thus obtained the following formally exact sequence:

0 Θ E 2 D F 0 6 D 1 F 1 1 D 2 F 2 1 D 3 F 3 1 D 4 F 4 0

0 Θ 1 2 2 6 21 1 46 1 36 1 10 0

with 1 2 + 21 46 + 36 10 = 0 , a part of it being a Janet sequence as before. Similarly, we get:

d i m ( R r + 11 ) = 1 ( r + 12 ) ( r + 13 ) ( r + 14 ) / 6 2 ( r + 10 ) ( r + 11 ) ( r + 12 ) / 6 + 21 ( r + 4 ) ( r + 5 ) ( r + 6 ) / 6 46 ( r + 3 ) ( r + 4 ) ( r + 5 ) / 6 + 36 ( r + 2 ) ( r + 5 ) ( r + 6 ) / 6 10 ( r + 1 ) ( r + 2 ) ( r + 3 ) / 6 = 18

or even d i m ( R r + 6 ) = 18, r 0 as a striking result indeed that can be checked directly through the exact sequences:

0 R 8 J 8 ( E ) J 6 ( F 0 ) F 1 0

0 18 165 168 21 0

0 R 6 J 6 ( E ) J 4 ( F 0 ) Q 4 0

0 18 84 70 4 0

0 R 5 J 5 ( E ) J 3 ( F 0 ) Q 3 0

0 17 56 40 1 0

and comes from the fact that d i m K ( R ) = 6 < or, equivalently, d i m K ( M ) = 6 < .

The reader not familiar with the formal theory of differential systems or modules may be surprised by the fact the two dimensions just found do not coincide at all because 6 < 18 . However, we have indeed d i m ( R 3 ) = d i m ( R 2 ) = d i m ( R 2 ( 4 ) ) = 6 and the exact sequence:

0 g 4 R 4 R 3 R 3 / R 3 ( 1 ) 0

0 4 15 12 1 0 d i m ( R 3 ( 1 ) ) = 12 1 = 11

showing that d i m ( R 3 ( 1 ) ) < d i m ( R 3 ) < d i m ( R 4 ) with 11 < 12 < 15 .However, we have the general Theorem 2.A.7 in ( [20] ) providing the useful prolongation/ projection (PP) procedure, namely that we have ρ r ( R q ( 1 ) ) = R q + r ( 1 ) , r 0 whenever the symbol g q of R q is 2-acyclic. In the present case, we have indeed ρ r ( R 3 ( 1 ) ) = R r + 3 ( 1 ) , r 0 because g 3 is known to be involutive, and the final system R 3 ( 4 ) is involutive with zero symbol, providing R 2 ( 4 ) which is only FI but with dimension 6. This situation is quite tricky indeed because prolongations are filling up successively the PD equations of order 2, then 3 and so on, adding therefore:

{ y 12 = 0 } , { y 123 = 0 , y 122 = 0 , y 112 = 0 } , { y 11 = 0 } , { y 113 = 0 , y 112 = 0 , y 113 = 0 } ,

5. Conclusions

When a differential operator D of order q is given, the problem of finding its compatibility conditions (CC) is to look for a new operator D 1 of a certain order s such that D 1 η = 0 must be satisfied in order to be able to solve the inhomogeneous system D ξ = η . This is an old problem first solved as a footnote by M. Janet in 1920 ( [10] [11] [12] [17] ) and finally studied by D.C. Spencer with collaborators around 1970 ( [19] [20] [21] ). The main idea is to construct a finite length differential sequence by repeating this procedure anew with D 1 and so on till one eventually ends with D n according to Janet when n is the number of independent variables. It soon became clear that constructing D 1 is largely depending on various intrinsic properties of D .

・ If D is involutive, then D 1 , , D n are first order involutive operators in the corresponding Janet sequence that can be constructed “step by step” as above but also “as a whole” like in the Poincaré sequence for the exterior derivative.

・ If D is only formally integrable (FI), that is all the equations of order q + r of the corresponding homogeneous system can be obtained by only r prolongations, then the order of D 1 is s + 1 when s is the smallest integer such that the symbol of order q + s becomes 2-acyclic. Such a result is still not acknowledged today by physicists even though it is essential for studying the conformal Killing system of space-time in general relativity.

・ If D is not even FI, not only the construction of D 1 may become very difficult but also a strange phenomenon may appear, namely one can start to find CC of order s 1 , then no new CC other than the ones generated by these CC up to order s 2 > s 1 when suddenly new generating CC may appear, generating all the CC up to order s 3 > s 2 and so on till the procedure ends.

This delicate question has been recently raised by physicists in order to study the system of Killing equations for certain useful metrics solutions of the Einstein equations in vacuum (Minkowski gives s = 2 while Schwarzschild gives s 1 = 2 , s 2 = 3 ). Needless to say that computer algebra is of quite a poor help in this case because the dimensions of the jet spaces and the size of the matrices involved (up to 420 × 460 for r = 0 in the last example) may increase drastically ( [26] ).

The aim of this paper has been first to provide illustrating examples of the above situations and one of them with s 1 = 3 , s 2 = 6 seems to be the only one known in the literature today. In addition, we have solved the (general) generating problem by using new differential homological algebraic methods, with the hope that computer algebra will soon become of some help in a near future.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Pommaret, J.-F. (2013) Journal of Modern Physics, 4, 223-239.
https://doi.org/10.4236/jmp.2013.48A022
[2] Pommaret, J.-F. (2017) Journal of Modern Physics, 8, 2122-2158.
https://doi.org/10.4236/jmp.2017.813130
[3] Pommaret, J.-F. (2018) New Mathematical Methods for Physics. NOVA Science Publisher, New York.
[4] Pommaret, J.-F. (2018) Journal of Modern Physics, 9, 1970-2007.
https://doi.org/10.4236/jmp.2018.910125
https://arxiv.org/abs/1805.11958v2
[5] Andersson, L. (2015) Spin Geometry and Conservation Laws in the Kerr Spacetime.
https://arxiv.org/abs/1504.02069
[6] Shah, A.G., Whitting, B.F., Aksteiner, S., Andersson, L. and Bäckdahl, T. (2016) Gauge Invariant perturbation of Schwarzschild Spacetime.
https://arxiv.org/abs/1611.08291
[7] Khavkine, I. (2017) Journal of Geometry and Physics, 113, 131-169.
https://arxiv.org/abs/1409.7212
https://doi.org/10.1016/j.geomphys.2016.06.009
[8] Aksteiner, S. and Bäckdahl, T. (2018) All Local Gauge Invariants for Perturbations of the Kerr Spacetime.
https://arxiv.org/abs/1803.05341
https://doi.org/10.1103/PhysRevLett.121.051104
[9] Khavkine, I. (2018) Compatibility Complexes of Overdetermined PDEs of Finite Type, with Applications to the Killing Equations.
https://arxiv.org/abs/1805.03751
[10] Pommaret, J.-F. (2001) Partial Differential Control Theory, Kluwer, Dordrecht (1000 p).
[11] Pommaret, J.-F. (1978) Systems of Partial Differential Equations and Lie Pseudogroups, Gordon and Breach, New York. Russian Translation: MIR, Moscow, 1983.
[12] Pommaret, J.-F. (1994) Partial Differential Equations and Group Theory. Kluwer.
https://doi.org/10.1007/978-94-017-2539-2
[13] Pommaret, J.-F. (2016) Deformation Theory of Algebraic and Geometric Structures. Lambert Academic Publisher (LAP), Saarbrucken, Germany.
https://arxiv.org/abs/1207.1964
[14] Cartan, E. and Einstein, A. (1979) Letters on Absolute Parallelism. Princeton University Press, Princeton.
[15] Kashiwara, M. (1995) Algebraic Study of Systems of Partial Differential Equations, Mémoires de la Société Mathématique de France 63 (Transl. from Japanese of His 1970 Master’s Thesis).
[16] Schneiders, J.-P. (1994) Bulletin de la Societe Royale des Sciences de Liege, 63, 223-295.
[17] Janet, M. (1920) Journal de Mathématiques, 8, 65-151.
[18] Pommaret, J.-F. (2013) Multidimensional Systems and Signal Processing, 26, 405-437.
https://doi.org/10.1007/s11045-013-0265-0
[19] Spencer, D.C. (1965) Bulletin of the American Mathematical Society, 75, 1-114.
[20] Goldschmidt, H. (1968) Annales scientifiques de l’école normale supérieure, 4, 417-444.
https://doi.org/10.24033/asens.1168
[21] Goldschmidt, H. (1968) Annales scientifiques de l’école normale supérieure, 4, 617-625.
https://doi.org/10.24033/asens.1173
[22] Northcott, D.G. (1966) An Introduction to Homological Algebra. Cambridge University Press, Cambridge.
[23] Rotman, J.J. (1979) An Introduction to Homological Algebra, Pure and Applied Mathematics. Academic Press, Cambridge.
[24] Macaulay, F.S. (1916) The Algebraic Theory of Modular Systems, Cambridge.
[25] Pommaret, J.-F. (1988) Lie Pseudogroups and Mechanics. Gordon and Breach, New York.
[26] Quadrat, A. and Robertz, R. (2014) Acta Applicandae Mathematicae, 133, 187-234.
https://hal-supelec.archives-ouvertes.fr/hal-00925533
https://doi.org/10.1007/s10440-013-9864-x

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.