Gravitational Waves and Black Holes: Beyond the Mirror

# Introduction # INTRODUCTION When M. Janet introduced in 1920 the first finite length differential sequence as a footnote of his paper \[1\], he surely did not know about the possibility to use such a sequence in elasticity theory along the way introduced by the brothers E. and F. Cosserat in 1909 \[2\]. Being a visiting student of D. C. Spencer (1912-2001) at Princeton University in 1970, I discovered that he was not even knowing the mathematical foundations of general relativity (GR) studied by his close friend J. A. Wheeler (1911-2008) who was offering 1000 dollars at that time to anybody finding a potential for Einstein equations in vacuum. This possibility is known to exist for Maxwell equations in electromagnetism (EM), usually defined by $d A = F$ while introducing the exterior derivative. I discovered in 1995 the negative solution of this challenge by using (formal) adjoint operators in a systematic way (See ideXlab on the net !), contrary to the general belief of the GR community. As a byproduct, such a result can only be found in books of control theory \[3, 4\]. In 1980, I met Janet who was still alive and he told me about the work of E. Vessiot and the resulting “affair” concerning the Differential Galois Theory \[5, 6\]. In 2015, a few meetings have been organized by the “ Institut Henri Poincaré ” (IHP) in Paris during three months about “ Mathematical General relativity ” . In particular, a Celebration of the 100th Anniversary was held on 16-20 November, largely dedicated to gravitational waves (GW) and followed by two days in honor of A. Lichnerowicz on $1 9 + 2 0$ December. As Lichnerowicz had been my main advisor during more than 20 years and died in 1998, I decided to participate to all these events. The atmosphere was very unpleasant because everybody knew that most sponsors should stop funding. One invited talk ” Are Black Holes Real ” given by S. Klainermann was so badly accepted that I could only answer to my neighbour, a young foreign student, that I was listening to it for the first time. The idea was to distinguish between three kinds of “Reality”, namely “Virtual reality, Physical reality and Mathematical reality”, the latter allowing to write a paper without any mathematical mistake but no comment was done on the mathematical assumptions used at the beginning \[7, 8\]. Less than 6 months later, LIGO announced to have detected GW produced by a couple of merging binary black holes (!) and this event, highly spread in newspapers, has been followed by the diffusion of pictures of black holes \[9\]. Since that time, I started to have doubts and, being specialist of control theory, I decided to use my knowledge for studying the origin of GW. In 2017, I discovered why GW cannot exist because Einstein, copying Beltrami, both ignoring that the Einstein operator, linearization of the Einstein tensor over the M metric, was surprisingly self-adjoint\[10, 11\]. I started to have doubts, not about the proper detection but mainly about the defining equations. Then, I started to have serious doubts when LIGO did stop for 3 years and I don’t speak about the lack of any result for KAGRA after spending 250 millions of dollars. It is at this moment that I decided to care about black holes while taking into account a few recent papers I wrote about the comparison of the M,S and K metrics \[12, 13\] but also as a way to disagree with the approach used by L. Andersson and collaborators met while lecturing at the Albert Einstein Institute (AEI) of Potsdam (October 23-27, 2017) \[14\]. In the Special Relativity paper of Einstein (1905), only a footnote provides a reference to the conformal group of space-time, namely the group of transformations preserving the Minkowski metric $\omega$ up to a function factor, but there is no proof that the conformal factor should be equal to 1. Over a manifold of dimension $n \geq 3$ , this group has $n$ translations, $n ( n - 1 ) / 2$ rotations, 1 dilatation and $n$ non-linear elations introduced by $\mathrm { E }$ . Cartan in 1922, with a total number of $N = ( n + 1 ) ( n + 2 ) / 2$ parameters that is $N = 1 5$ when $n = 4$ \[15, 16\]. However, it is also the number of the Cauchy stress equations (1823), the Cosserat couple-stress equations (1909), the only Clausius virial equation (1870), the Maxwell (1873) and Weyl (1918) equations which are among the most famous partial differential equations that can be found today in any textbook dealing with elasticity theory, continuum mechanics, thermodynamics or electromagnetism. The purpose of this paper is to prove that the form of these equations only depends on the structure of the conformal group for an arbitrary $n \geq 1$ because they are described as a whole by the (formal) adjoint of the first Spencer operator existing in the Spencer differential sequence. Such a group theoretical implication is obtained by applying totally new differential geometric methods in field theory. In particular, when $n = 4$ , the main idea is to enlarge the group from 10 up to 11 or 15 parameters by using the Weyl or conformal group instead of the Poincaré group of space-time. Contrary to the Einstein equations, these equations can be all parametrized by the adjoint of the second Spencer operator through $N = n(n - 1)/2$ potentials. These results bring the need to revisit the mathematical foundations of both General Relativity (GR) and Gauge Theory (GT) according to a clever but rarely quoted paper of H. Poincaré (1901) \[17\]. They strengthen the comments we already made about the dual confusions made by Einstein (1915) while following Beltrami (1892), both using the same operator when $n = 3$ and $n = 4$ but ignoring it is self-adjoint in the framework of differential double duality. They also question the origin and existence of black holes. With some more preliminary details provided in \[18\], we have successively: When $n = 2$ in plane elasticity with Riemann operator $\Omega \to d _ { 2 2 } \Omega _ { 1 1 } + d _ { 1 1 } \Omega _ { 2 2 } - 2 d _ { 1 2 } \Omega _ { 1 2 }$, G.B. Airy found in 1863 the possibility to parametrize the $C a u c h y \ = \ a d ( K i l l i n g )$ operator $\sigma \rightarrow$ $$ \fbox{ $ \begin{array}{rcl} d_1\sigma^{11} + d_2\sigma^{12} & = & 0 \\ d_1\sigma^{21} + d_2\sigma^{22} & = & 0 \end{array} $ } $$ $$ \fbox{ $ \begin{array}{rcccl} & & 1 & & \\ & \nearrow & & \searrow & \\ 2 & \stackrel{\cal{D}}{\longrightarrow} & 3 & \stackrel{{\cal{D}}_1}{\longrightarrow} & 1 \\ & & & & \\ 2 & \stackrel{ad(\cal{D})}{\longleftarrow} & 3 & \stackrel{ad({\cal{D}}_1)}{\longleftarrow} & 1 \end{array} $ } $$ with $\sigma ^ { 1 2 } = \sigma ^ { 2 1 }$ by a single stress function $\phi$ wearing his name through the $A i r y = a d ( R i e m a n n )$ operator $\phi \to ( d _ { 2 2 } \phi , - d _ { 1 2 } \phi , d _ { 1 1 } \phi )$ as follows: $$ \fbox{ $\begin{array}{ccccccccc} & & 2 & \stackrel{D}{\longrightarrow} & 3 & \stackrel{D_1}{\longrightarrow} & 1 & \longrightarrow & 0 \\ & & & & & & & & \\ 0 & \longleftarrow & 2 & \stackrel{ad(D)}{\longleftarrow} & 3 & \stackrel{ad(D_1)}{\longleftarrow} & 1 & & \end{array}$ } $$ Multiplying the Riemann operator on the left by $\phi$ and “integrating by parts”, we obtain: $$ \phi ( d _ { 2 2 } \Omega _ { 1 1 } ) = ( d _ { 2 2 } \phi ) \Omega _ { 1 1 } + d _ { 2 } ( \phi d _ { 2 } \Omega _ { 1 1 } - ( d _ { 2 } \phi ) \Omega _ { 1 1 } ) $$ As for the fact that the Cauchy operator is the adjoint of the Killing operator, this result can be found in any textbook of elasticity. When $n = 3$ in space elasticity, E. Beltrami found in 1892 the possibility to parametrize the $C a u c h y \ = \ a d ( K i l l i n g )$ operator by means of six stress functions $\phi _ { i j } ~ = ~ \phi _ { j i }$ wearing his name through the self-adjoint $B e l t r a m i = a d ( R i e m a n n )$ operator as follows: $$ \fbox{ $\begin{array}{ccccccccccc} & & 3 & \xrightarrow[\;1\;]{\textit{Killing}} & 6 & \xrightarrow[\;2\;]{\textit{Riemann}} & 6 & \xrightarrow[\;1\;]{\textit{Bianchi}} & 3 & \rightarrow & 0 \\ 0 & \leftarrow & 3 & \xleftarrow[\;1\;]{\textit{Cauchy}} & 6 & \xleftarrow[\;2\;]{\textit{Beltrami}} & 6 & \xleftarrow[\;1\;]{\textit{Lanczos}} & 3 & & \end{array}$ } $$ As can be checked, the alternate sum of dimensions in each sequence, called Euler-Poincaré characteristic, does vanish indeed. More generally, studying the Lanczos problems in 2001 for a dimension $n \geq 3$ , $\mathrm { I }$ discovered that the $B e l t r a m i = a d ( R i e m a n n )$ operator can be parametrized by the $L a n c z o s = a d ( B i a n c h i )$ operator in the geometrical and physical adjoint sequences made by operators acting on tensors, giving order of operators and number of components as follows: The geometrical and adjoint physical long exact differential sequences of operators acting on tensors, giving order of operators and number of components as follows: $$ \begin{array}{ccccccccccc} & & n & \xrightarrow[\;1\;]{\textit{Killing}} & \frac{n(n+1)}{2} & \xrightarrow[\;2\;]{\textit{Riemann}} & \frac{n^2(n^2-1)}{12} & \xrightarrow[\;1\;]{\textit{Bianchi}} & \frac{n^2(n^2-1)(n-2)}{24} & \rightarrow & 0 \\ \\ 0 & \leftarrow & n & \xleftarrow[\;1\;]{\textit{Cauchy}} & \frac{n(n+1)}{2} & \xleftarrow[\;2\;]{\textit{Beltrami}} & \frac{n^2(n^2-1)}{12} & \xleftarrow[\;1\;]{\textit{Lanczos}} & \frac{n^2(n^2-1)(n-2)}{24} & & \end{array} $$ When $n = 4$ in space-time, Einstein, probably knowing the work of Beltrami because the comparison needs no comment (\[18\], Proposition 4.1, p 28), made a terrible confusion in 1915 between the $C a u c h y = a d ( K i l l i n g )$ operator and the div operator induced from the Bianchi operator because both have been using the same Einstein operator but ignoring that such an operator is self-adjoint in the framework of differential double duality when $n \geq 3$ (See \[10\] for more details). Looking at \[10\] in order to understand the origin of GW or at section 4, the second order linearized Ricci operator $S _ { 2 } T ^ { * } \to S _ { 2 } T ^ { * } : ( \Omega _ { i j } ) \to ( R _ { i j } )$ with 4 terms is defined by the formula: $$ 2 R _ { i j } = \omega ^ { r s } ( d _ { r s } \Omega _ { i j } + d _ { i j } \Omega _ { r s } - d _ { r i } \Omega _ { s j } - d _ { s j } \Omega _ { r i } ) = 2 R _ { j i } $$ Let us now introduce the linear map $C : S _ { 2 } T ^ { * } \to S _ { 2 } T ^ { *} : \Omega _ { i j } \to \bar { \Omega } _ { i j } = \Omega _ { i j } - \frac { 1 } { 2 } \omega _ { i j } \omega ^ { r s } \Omega _ { r s}$ where $\Omega$ is a perturbation of the Minkowski metric $\omega$ , invertible if and only if $n \geq 3$ . It is well known that the Einstein operator is defined by the same map $C : R _ { i j } \to R _ { i j } - \frac { 1 } { 2 } \omega _ { i j } \omega ^ { r s } R _ { r s } = E _ { i j }$ not depending on any conformal factor. Comparing to \[10\] or to *any other textbook*, the GW are defined by the (strange!) composite operator $\mathcal { X } : \bar { \Omega } \xrightarrow{\;Einstein\;} E$ in such a way that $Einstein : \mathcal { X } \circ C$. Taking the respective adjoint operators, remembering that $Einstein = a d ( Einstein )$ and that we have $a d ( P \circ Q ) = a d ( Q ) \circ a d ( P )$ whenever $P , Q$ are two operators, we obtain: $$ \begin{aligned} Einstein &= a d ( C ) \circ a d ( \mathcal { X } ) = C \circ a d ( \mathcal { X } ) = C \circ Ricci \\ &\Rightarrow a d ( \mathcal { X } ) = Ricci \Rightarrow \mathcal { X } = a d ( Ricci ) \end{aligned} $$ Introducing the test functions $\lambda ^ { i j } = \lambda ^ { j i }$ and setting as usual $\Pi = \omega ^ { i j } d _ { i j }$ , we get: $$ d_{ij} ( {\omega}^{ij} {\lambda}^{rs} + {\omega}^{rs} {\lambda}^{ij} - {\omega}^{sj} {\lambda}^{ri} - {\omega}^{ri} {\lambda}^{sj} ) = {\sigma}^{rs} \Rightarrow d_r({\sigma}^{rs}) = 0 $$ which is exactly the generalized version of the Beltrami parametrization described in the preceding diagram with 10 instead of 20 stress functions when $n = 4$ . However, the above results are showing out the close links existing between group theory and differential sequences. For this reason, we briefly recall the historical framework leading to these new results and the reason for which we have not been able to quote many external references. This paper is also written as a smile to the famous English writer Lewis Carroll who used these words in the novel he wrote in 1871, six years after ” Alice in Wonderland”. The concept of ” group”, introduced in mathematics for the first time by E. Galois (1830), slowly passed from algebra to geometry with the work of S. Lie on Lie groups (1880) and Lie pseudogroups (1890) of transformations. The concept of a finite length differential sequence, now called Janet sequence, has been described for the first time as a footnote by M. Janet (1920). Then, the work of D. C. Spencer (1970) has been the first attempt to use the formal theory of systems of partial differential equations in order to study the formal theory of Lie pseudogroups \[20, 21\]. However, the linear and nonlinear Spencer sequences for Lie pseudogroups, though never used in physics, largely supersede the ” Cartan structure equations ” (1905) and are quite different from the ” Vessiot structure equations ” (1903), introduced for the same purpose but still not known today because they have never been acknowledged by E. Cartan and successors \[5, 6\]. This diagram explains why it has never been possible to connect Gauge Theory (GT) using Maurer-Cartan (MC) equations with torsion $^ +$ curvature with GR using (Riemann curvature alone: $$ \fbox{ $\begin{array}{ccccccc} & & CARTAN & \longrightarrow & SPENCER & \hspace {1cm} & GT \\ & \nearrow & & & & \\ LIE & & \updownarrow & Cauchy, Cosserat, Maxwell & \updownarrow & & \\ & \searrow & & & & \\ & & VESSIOT & \longrightarrow & JANET & \hspace{1cm} & GR \end{array}$ } $$ Example 1.1: In order to explain the difference existing between a Lie group and a Lie pseudogroup of transformations, let us consider the Lie group of projective transformations of the real line and differentiate three times the local action law as follows with $a , b , c , d = c s t$ with $d \neq 0$ : $$ \begin{aligned} y &= \frac { a x + b } { c x + d } \\ \Rightarrow\quad y _ { x } &= \frac { a d - b c } { ( c x + d ) ^ { 2 } } \\ \Rightarrow\quad y _ { x x } &= - 2 \frac { ( a d - b c ) c } { ( c x + d ) ^ { 3 } } \\ \Rightarrow\quad y _ { x x x } &= 6 \frac { ( a d - b c ) c ^ { 2 } } { ( c x + d ) ^ { 4 } } \end{aligned} $$ and get the single Schwarzian third order nonlinear OD equation $( y _ { x x x } / y _ { x } ) - ( 3 / 2 ) ( y _ { x x } / y _ { x } ) ^ { 2 } = 0$ . Setting then $y = x + t \xi + . . .$ and linearizing at the identity $y = x$ when $t 0$ we obtain the infinitesimal Lie equation $\xi _ { x x x } = 0 \Leftrightarrow \partial _ { x x x } \xi = 0$ $$ \begin{array}{ccccccccccccc} & & 0 &\rightarrow & \Theta & \stackrel{j_3}{\rightarrow} & 3 & \stackrel{ d_1}{ \rightarrow} & 3 & \rightarrow & 0 \\ & & & & & & \downarrow & & \parallel & & & \\ & & 0 &\rightarrow & 1 & \stackrel{j_3}{\rightarrow} & 4 & \stackrel{ d_1}{ \rightarrow} & 3 & \rightarrow & 0 \\ & & & & \parallel & & \downarrow & & \downarrow & & & \\ 0 & \rightarrow & \Theta & \rightarrow & 1 & \stackrel{D}{\rightarrow}& 1 & \longrightarrow & 0 & & & \end{array} $$ with a basis of solutions $\{ \theta _ { \tau } \} = \{ \partial _ { x } , x \partial _ { x } , { \textstyle \frac { 1 } { 2 } } x ^ { 2 } \partial _ { x } \}$ generating the space of solutions $\Theta \subset T$ over the constants with $[ \Theta , \Theta ] \subset \Theta$ . We have the three abelian subgroups $\ y = a x , y = x + b , { \frac { 1 } { y } } = { \frac { 1 } { x } } + c$ $n = 1$ with no rotation but $N = 3$ . However, the next examples will prove that no group theoretical technique may be applied when solutions cannot be known. Example 1.2: With $m = n = 2 , q = 1$ and $\omega = ( \alpha \in T ^ { * } , \beta \in \wedge ^ { 2 } T ^ { * } ,$ ), let us consider the Lie operator ${ \mathcal { D } } : T \to \Omega : \xi \to { \mathcal { L } } ( \xi ) \omega = ( A = { \mathcal { L } } ( \xi ) \alpha , B = { \mathcal { L } } ( \xi ) \beta )$ and the first order involutive system: The only Vessiot structure equation is $d \alpha = c \beta$ where now $d$ is the exterior derivative and the only Vessiot structure constant is $c = c s t$ \[5\]. We let the reader check that this system is formally integrable (FI), that is no new PD equation can be obtained by differentiating these three equations in any way, if and only if this condition is satisfied. We have the differential sequence: $$ ~ 0 \to \Theta \to \xi \xrightarrow { \mathcal { D } } \eta \xrightarrow { \mathcal { D } _ { 1 } } \zeta \to 0 $$ with $\eta = ( A _ { 1 } , A _ { 2 } , B )$ and $\zeta \in \wedge ^ { 2 } T ^ { * }$ and we may look for the adjoint differential sequence. Multiplying $( A _ { 1 } , A _ { 2 } , B )$ respectively by $( \mu ^ { \mathrm { 1 } } , \mu ^ { \mathrm { 2 } } , \mu ^ { \mathrm { 3 } } )$ and integrating by parts, we obtain the adjoint operator $a d ( \mathcal { D } )$ in the form: $$ - \alpha _ { 1 } ( \partial _ { 1 } \mu ^ { 1 } + \partial _ { 2 } \mu ^ { 2 } ) - \beta ( \partial _ { 1 } \mu ^ { 3 } - c \mu ^ { 2 } ) = \nu ^ { 1 } , ~ - \alpha _ { 2 } ( \partial _ { 1 } \mu ^ { 1 } + \partial _ { 2 } \mu ^ { 2 } ) - \beta ( \partial _ { 2 } \mu ^ { 3 } + c \mu ^ { 1 } ) = \nu ^ { 2 } $$ Then, multiplying $\zeta = \dot { \partial } _ { 1 } A _ { 2 } - \dot { \partial } _ { 2 } A _ { 1 } - c B$ by $\lambda$ , we obtain $a d ( \mathcal { D } _ { 1 } )$ as: $$ \partial _ { 2 } \lambda = \mu ^ { 1 } , ~ - \partial _ { 1 } \lambda = \mu ^ { 2 } , ~ - c \lambda = \mu ^ { 3 } $$ and obtain the adjoint differential sequence: $\nu \leftarrow \mu \leftarrow \lambda$. in which $a d ( \mathcal { D } )$ may not generate the compatibility conditions (CC) of $a d ( \mathcal { D } _ { 1 } )$ We have therefore to consider the two cases \[4\]: $\bullet$ $c = 0$ : With $\alpha = d x ^ { 1 } , \beta = d x ^ { 1 } \wedge d x ^ { 2 }$ , we have $d \alpha = 0$ and $a d ( \mathcal { D } _ { 1 } )$ is not injective with kernel $\lambda ~ = ~ c s t$ but $a d ( \mathcal { D } )$ is surjective. The Lie pseudogroup $\Gamma$ is made by transfomations $y ^ { 1 } = x ^ { 1 } + a , y ^ { 2 } = x ^ { 2 } + f ( x ^ { 1 } )$ with $a = c s t$ and $f ( x ^ { 1 } )$ arbitrary. $\bullet$ $c \neq 0$ : Like Vessiot himself, we may choose $\alpha = x ^ { 2 } d x ^ { 1 } , \beta : d x ^ { 1 } \wedge d x ^ { 2 }$ and obtain $d \alpha + \beta = 0$ , that is $c = - 1$ . We check that $a d ( \mathcal { D } _ { 1 } )$ is now injective but that $a d ( \mathcal { D } )$ is still surjective. The Lie pseudogroup $\Gamma$ is made by transformations $y ^ { \mathrm { { I } } } = f ( x ^ { \mathrm { { I } } } ) , y ^ { \mathrm { { ' 2 } } } = x ^ { \mathrm { { 2 } } } / ( \partial f / \partial x ^ { \mathrm { { 1 } } } )$ with $f$ invertible. We have thus discovered that the properties of the adjoint sequence largely depend on the Vessiot structure constant $c$ . Finally, if $\alpha = x ^ { 2 } d x ^ { 1 } - x ^ { 1 } d x ^ { 2 }$ , no explicit solution may be provided. Example 1.3: (Contact transformations) With $m = n = 3 , q = 1$ and ground differential field $K = \mathbb { Q } ( x )$ , we may introduce the 1-form $\alpha = d x ^ { 1 } - x ^ { 3 } d x ^ { 2 } \in T ^ { * }$ and consider the system of infinitesimal Lie equations defined by ${ \mathcal { L } } ( \xi ) ( \alpha ) = \rho ( x ) \alpha$ after eliminating the factor $\rho$ . However, we notice that $\alpha$ is not invariant by the contact Lie pseudogroup and cannot be considered as an invariant associated geometric object. In fact, it is known that the corresponding geometric object is a 1-form density $\omega$ leading to the system of infinitesimal Lie equations in Medolaghi form: $$ \Omega _ { i } \equiv ( \mathcal { L } ( \xi ) \omega ) _ { i } \equiv \omega _ { r } \partial _ { i } \xi ^ { r } - \frac { 1 } { 2 } \omega _ { i } \partial _ { r } \xi ^ { r } + \xi ^ { r } \partial _ { r } \omega _ { i } = 0 $$ and to the only Vessiot structure equation \[6\]: $$ \omega _ { 1 } ( \partial _ { 2 } \omega _ { 3 } - \partial _ { 3 } \omega _ { 2 } ) + \omega _ { 2 } ( \partial _ { 3 } \omega _ { 1 } - \partial _ { 1 } \omega _ { 3 } ) + \omega _ { 3 } ( \partial _ { 1 } \omega _ { 2 } - \partial _ { 2 } \omega _ { 1 } ) = c $$ with the only structure constant $c$ . We point out the fact that only the condition of constant Riemannian curvature is known today to be similar and we have explained why such a situation has been produced deliberately by the successors of Cartan. In the present contact situation, we may choose $\omega = ( 1 , - x ^ { 3 } , 0 )$ and get $c = 1$ but we may also choose $\omega = ( 1 , 0 , 0 )$ and get $c = 0$ , these two choices both bringing an involutive system. Our problem will be now to construct the differential sequences: $\xi \xrightarrow { \mathcal { D } } \Omega \xrightarrow { \mathcal { D } _ { 1 } } \zeta$ and its adjoint sequence: $\nu \xleftarrow{ad(\mathcal D)} \mu \xleftarrow{ad(\mathcal D_1)} \lambda$ Linearizing the only Vessiot structure equation, we get the corresponding CC system $\mathcal { D } _ { 1 } \Omega = 0$ $$ \begin{array} { r } { \omega _ { 1 } ( \partial _ { 2 } \Omega _ { 3 } - \partial _ { 3 } \Omega _ { 2 } ) + \omega _ { 2 } ( \partial _ { 3 } \Omega _ { 1 } - \partial _ { 1 } \Omega _ { 3 } ) + \omega _ { 3 } ( \partial _ { 1 } \Omega _ { 2 } - \partial _ { 2 } \Omega _ { 1 } ) } \\ { + ( \partial _ { 2 } \omega _ { 3 } - \partial _ { 3 } \omega _ { 2 } ) \Omega _ { 1 } + ( \partial _ { 3 } \omega _ { 1 } - \partial _ { 1 } \omega _ { 3 } ) \Omega _ { 2 } + ( \partial _ { 1 } \omega _ { 2 } - \partial _ { 2 } \omega _ { 3 } ) \Omega _ { 3 } = 0 } \end{array} $$ Multiplying on the left by a test function $\lambda$ and integrating by parts, we get the operator $a d ( \mathcal { D } _ { 1 } )$ in the form: $$ \left\lbrace \begin{array} { l l l } { \Omega _ { 1 } \to } & { \omega _ { 3 } \partial _ { 2 } \lambda - \omega _ { 2 } \partial _ { 3 } \lambda + 2 ( \partial _ { 2 } \omega _ { 3 } - \partial _ { 3 } \omega _ { 2 } ) \lambda } & { = \mu ^ { 1 } } \\ { \Omega _ { 2 } \to } & { \omega _ { 1 } \partial _ { 3 } \lambda - \omega _ { 3 } \partial _ { 1 } \lambda + 2 ( \partial _ { 3 } \omega _ { 1 } - \partial _ { 1 } \omega _ { 3 } ) \lambda } & { = \mu ^ { 2 } } \\ { \Omega _ { 3 } \to } & { \omega _ { 2 } \partial _ { 1 } \lambda - \omega _ { 1 } \partial _ { 2 } \lambda + 2 ( \partial _ { 1 } \omega _ { 2 } - \partial _ { 2 } \omega _ { 1 } ) \lambda } & { = \mu ^ { 3 } } \end{array} \right. $$ We obtain therefore the crucial formula $2 c \lambda = \omega _ { i } \mu ^ { i }$ showing how the previous sequences are essentially depending on the Vessiot structure constant $c$ . $\bullet$ Indeed, if $c \neq 0$ , then $\mu = 0 \Rightarrow \lambda = 0$ and the operator $a d ( \mathcal { D } _ { 1 } )$ is injective. This is the case when $\omega = ( 1 , - x ^ { 3 } , 0 ) \Rightarrow c = 1 \Rightarrow \lambda = 0$ . $\bullet$ On the contrary, if $c = 0$ , then the operator $a d ( { \mathcal { D } _ { 1 } } )$ may not be injective as can be seen by choosing $\boldsymbol { \omega } = ( 1 , 0 , 0 )$ . Indeed, in this case we get a kernel defined by $\partial _ { 3 } \lambda = 0 , \partial _ { 2 } \lambda = 0$ . Finally, unimodular contact transformations are preserving the 1-form $\alpha = d x ^ { 1 } - x ^ { 3 } d x ^ { 2 }$ , thus also the 2-form $\beta = d \alpha = d x ^ { 2 } \wedge d x ^ { 3 }$ and even the 3-form $\alpha \wedge \beta = d x ^ { 1 } \wedge d x ^ { 2 } \wedge d x ^ { 3 }$ . The Vessiot structure equations for the geometric object $\omega = ( \alpha , \beta )$ are now $d \alpha = c ^ { \prime } \beta , d \beta = c ^ { \prime \prime } \alpha \wedge \beta$ with $0 = d ^ { 2 } \alpha = c ^ { \prime } d \beta = c ^ { \prime } c " { \alpha } \wedge \beta$ and the only Jacobi condition $c ^ { \prime } c ^ { \prime \prime } = 0$ because $\alpha \wedge \beta \neq 0$ \[6\]. Remark 1.4: When $\omega$ is a non-degenerate metric, that is $d e t ( \omega ) \neq 0$ , it is well known since the work of L.P. Eisenhart on Riemannian geometry in 1926 that the first order Killing linear system defined by the $n ( n { + } 1 ) / 2$ PD equations $\Omega _ { i j } \equiv \omega _ { r j } ( x ) \partial _ { i } \xi ^ { r } + \omega _ { i r } ( x ) \partial _ { j } \xi ^ { r } + \xi ^ { r } \partial _ { r } \omega _ { i j } ( x ) = 0$ for a vector field $\xi = \xi ^ { r } \partial _ { r }$ , is formally integrable (FI), that is no new first order PD equation can be obtained iff the metric has a constant Riemannian curvature. We have proved in many books \[20, 21\] and papers that this is the only example of Vessiot equations known today and have explained in \[4\] that the reason for such a poor situation is related to a very unpleasant “Mathematical Affair” having to do with the differential Galois theory and involving the best french mathematicians of the beginning of the last century, namely H. Poincaré, E. Picard and G. Darboux (Original letters have been given directly to me by Janet because of the personal dedication of my first 1978 GB book and can be found in the library of ENS in Paris while a photocopy can be found in my 1988 GB book \[6\]. Meanwhile, mixing differential geometry with homological algebra, M. Kashiwara (1970) has created ” differential homological algebra”, in order to study differential modules by means of double duality and the corresponding extension modules (See \[21\] for references and Zbl 1079.93001). By chance, unexpected arguments have been introduced by the brothers E. and F. Cosserat (1909) in order to revisit elasticity and by H. Weyl (1918) in order to revisit electromagnetism through a unique differential sequence only depending on the structure of the conformal group. However, while the Cosserat brothers were only using (translations + rotations), Weyl has only been dealing with (dilatation $^ +$ elations) as we shall explain \[22, 23\]. After recalling the negative answer we already provided in 1995 \[19\], the main purpose of this paper is to use new techniques of group theory in order to revisit the mathematical foundations of general relativity (GR) and gauge theory (GT) that are leading to gravitational waves. We point out the fact that all the diagrams presented can be obtained by means of computer algebra while using recent packages developed by my former PhD student A. Quadrat and his collaborators \[24\]. The solution of this striking but difficult problem has been announced, as we already said, is a series of lectures given at the Albert Einstein Institute (AEI, Potsdam, October, 23-27, 2017) “General Relativity and Gauge Theory: Beyond the Mirror” (hal-01632085, 09/11/2017). We shall prove in the next section that the study of the Killing operator done in \[14\] by means of purely technical relativistic tools has in fact nothing to do with GR and can be solved only counting with fingers. We also invite the reader to look at the more applied presentation (arXiv: 2302.06585) \[25\]. In this second approach, we advise the reader to have a special look at the photo-elastic beam experiment showing the link that may exist between high level mathematical mathematical tools and their phenomenological approach also done by J.C. Maxwell himself. We finally say that a rough sketch of the Spencer operator for systems with constant coefficients has been provided by F.S. Macaulay in 1916 through ” Inverse Systems” in \[26\]. Before going ahead, let us prove that there may be only two types of differential sequences, the Janet sequence introduced by M. Janet in 1920 and the totally different Spencer sequence introduced by D. C. Spencer in 1970 though both only depend on the Spencer operator \[27, 28\]. Though the mathematical and physical communities still believe that a differential sequence must always be constructed ” step by step”, that is, starting with $\mathcal { D } \xi = \eta$ with generating CC $\mathcal { D } _ { 1 } \eta = 0$ , one may start anew with $\mathcal { D } _ { 1 } \eta = \zeta$ with generating CC $\mathcal { D } _ { 2 } \zeta = 0$ and so on. However, as we shall see, the Poincaré (also called de Rham out of France !) sequence for the exterior derivative may be defined “as a whole”, a fact that led people to believe that the central operator for constructing differential sequences is the exterior derivative, a wrong way indeed (See \[29\] $\mathrm { ~ p ~ } 1 8 5 + 3 9 1$ ). For this, if $E$ is a vector bundle over the base $X$ , we introduce the $q$ jet bundle $J _ { q } ( E )$ with sections $\xi _ { q } : ( x ) \to ( \xi ^ { k } ( x ) , \xi _ { i } ^ { k } ( x ) , \xi _ { i j } ^ { k } ( x ) , \ldots )$ transforming like the sections $j _ { q } ( \xi ) : ( x ) \to$ $( \xi ^ { k } ( x ) , \partial _ { i } \xi ^ { k } ( x ) , \partial _ { i j } \xi _ { i j } ^ { k } ( x ) , \ldots )$ . The Spencer operator $d : J _ { q + 1 } ( T ) \to T ^ { * } \otimes J _ { q } ( T )$ allows to compare these sections by considering the differences ${ } ( \partial _ { i } \xi ^ { k } ( x ) - \xi _ { i } ^ { k } ( x ) , \partial _ { i } \xi _ { j } ^ { k } ( x ) - \xi _ { i j } ^ { k } ( x ) , . . . )$ and so on. For any system $R _ { q } \subset J _ { q } ( E )$ , differentiating once all the given OD or PD equations, we obtain the first prolongation $R _ { q + 1 } \subset J _ { q + 1 } ( E )$ equations and the Spencer operator can be extended to an operator: $$ \displaystyle d:\wedge^{s}T^{*}\otimes R_{q+1}\rightarrow \wedge^{s+1}T^{*}\otimes R_{q}:(\xi^{k}_{\mu,I}(x)\,dx^{I})\rightarrow\big((\partial_{i}\xi^{k}_{\mu,I}(x)-\xi^{k}_{\mu+1_{i},I}(x))\,dx^{i}\wedge dx^{I}\big) $$ We use multi-indices $\boldsymbol { \mu } = ( \mu _ { 1 } , . . . , \mu _ { n } )$ with $\mu + 1 _ { i } = ( \mu _ { 1 } , . . . , \mu _ { i } + 1 , . . . , \mu _ { n } )$ and $\mid \mu \mid = \mu _ { 1 } + \ldots + \mu _ { n }$ both with standard multi-index notation for exterior forms and one can check that $d \circ d = 0$ . The Spencer operator and the exterior derivative are thus interlaced by the above formula. The restriction of $d$ to the terms of upper order defining the symbols $g _ { q + r } = R _ { q + r } \cap S _ { q + r } T ^ { * } \otimes E \subset J _ { q + r } ( E )$ is ” minus” the Spencer map $\delta : \wedge ^ { s } T ^ { * } \otimes g _ { q + 1 } \to \wedge ^ { s + 1 } T ^ { * } \otimes g _ { q }$ with $\delta \circ \delta = 0$ through the formula $\xi _ { \mu , I } ^ { k } d x ^ { I } \rightarrow \xi _ { \mu + 1 _ { i } , I } ^ { k } d x ^ { i } \wedge d x ^ { I }$ because $\xi _ { \mu + 1 _ { i } + 1 _ { j } } ^ { k } d x ^ { i } \wedge d x ^ { j } = 0$ when $\mid \mu \mid = q , \mid \nu \mid = q + 1$ The system $R _ { q }$ is said to be involutive if its symbol $g _ { q }$ is involutive, that is if all the $\delta$ sequences are exact, and the projection $\pi _ { q } ^ { q + 1 } : R _ { q + 1 } \to R _ { q }$ is an epimorphism. For any system $R _ { q } \subset J _ { q } ( E )$ , we may define $R _ { q + r } ^ { ( s ) } \subset R _ { q + r } \subset J _ { q + r } ( E )$ by differentiating $r + s$ times and keeping only the equations of order $q + r$. When $R _ { q }$ is involutive, we may define the Janet bundles $F _ { r }$, for $r = 0 , 1 , \ldots , n$, by the short exact sequences \[20, 29\]: $$ 0 \to \wedge ^ { r } T ^ { * } \otimes R _ { q } + \delta ( \wedge ^ { r - 1 } T ^ { * } \otimes S _ { q + 1 } T ^ { * } \otimes E ) \to \wedge ^ { r } T ^ { * } \times J _ { q } ( E ) \to F _ { r } \to 0 $$ We may pick up a section of $F _ { r }$ , lift it up to a section of $\land ^ { r } T ^ { * } \otimes J _ { q } ( E )$ that we may lift up to a section of $\land ^ { r } T ^ { * } \otimes J _ { q + 1 } ( E )$ and apply $d$ in order to get a section of $\land ^ { r + 1 } T ^ { * } \otimes J _ { q } ( E )$ that we may project onto a section of $F _ { r + 1 }$ in order to construct an operator $\mathcal { D } _ { r + 1 } : F _ { r } \to F _ { r + 1 }$ generating the CC of $\mathcal { D } _ { r }$ in the canonical linear Janet sequence: $$ 0 \longrightarrow \Theta \longrightarrow E \xrightarrow{\mathcal{D}} F_{0} \xrightarrow{\mathcal{D}_{1}} F_{1} \xrightarrow{\mathcal{D}_{2}} \cdots \xrightarrow{\mathcal{D}_{n}} F_{n} \longrightarrow 0 $$ If we have two involutive systems $R _ { q } \subset \hat { R } _ { q } \subset J _ { q } ( E )$ , the Janet sequence for $R _ { q }$ projects onto the Janet sequence for $\hat { R } _ { q }$ and we may define inductively canonical epimorphisms $F _ { r } \to \hat { F } _ { r } \to 0$ for $r = 0 , 1 , . . . , n$ by comparing the previous sequences for $R _ { q }$ and $\hat { R } _ { q }$ . A similar procedure can also be obtained if we define the Spencer bundles $C _ { r }$ for $r = 0 , 1 , . . . , n$ by the short exact sequences \[20, 29\]: $$ 0 \to \delta(\wedge^{r-1}T^{*}\otimes g_{q+1}) \to \wedge^{r}T^{*}\otimes R_{q} \to C_{r} \to 0 $$ We may pick up a section of $C _ { r }$ , lift it to a section of $\wedge ^ { r } T ^ { * } \otimes R _ { q }$ , lift it up to a section of $\land ^ { r } T ^ { * } \otimes R _ { q + 1 }$ and apply $d$ in order to construct a section of $\wedge ^ { r + 1 } \otimes R _ { q }$ that we may project to $C _ { r + 1 }$ in order to construct an operator $D _ { r + 1 } : C _ { r } \to C _ { r + 1 }$ generating the CC of $D _ { r }$ in the canonical linear Spencer sequence which is another completely different resolution of the set $\Theta$ of (formal) solutions of $R _ { q }$ . $$ \fbox{$\displaystyle 0 \longrightarrow \Theta \xrightarrow{j_{q}} C_{0} \xrightarrow{D_{1}} C_{1} \xrightarrow{D_{2}} C_{2} \xrightarrow{D_{3}} \cdots \xrightarrow{D_{n}} C_{n} \longrightarrow 0$} $$ It can be proved that the Spencer sequence for $R _ { q } \subset J _ { q } ( E )$ is the Janet sequence for $R _ { q + 1 } \subset J _ { 1 } ( R _ { q } )$ \[22\]. However, if we have two systems as above, the Spencer sequence for $R _ { q }$ is now contained into the Spencer sequence for $\hat { R } _ { q }$ and we may construct inductively canonical monomorphisms $0 \to C _ { r } \to \hat { C } _ { r }$ for $r = 0 , 1 , . . . , n$ by comparing the previous sequences for $R _ { q }$ and $\hat { R } _ { q }$ . Defining $F _ { 0 }$ and $\Phi = \Phi _ { 0 }$ by the short exact sequence $0 \to R _ { q } \to J _ { q } ( E ) \stackrel { \Phi } { \longrightarrow } F _ { 0 } \longrightarrow 0$ , all the previous results can be combined in the following essential commutative and exact diagram: $$ \fbox{ $ \begin{array}{ccccccc} 0&&0&&0&&\\ \downarrow&&\downarrow&&\downarrow&&\\ \wedge ^ { r - 1 } T ^ { * } \otimes g _ { q + 1 } & \xrightarrow{\delta} & \wedge ^ { r } T ^ { * } \otimes R _ { q } & \longrightarrow & C _ { r } & \longrightarrow & 0\\ \downarrow && \downarrow && \downarrow && \\ \wedge ^ { r - 1 } T ^ { * } \otimes S _ { q + 1 } T ^ { * } \otimes E & \xrightarrow{\delta} & \wedge ^ { r } T ^ { * } \otimes J _ { q } ( E ) & \longrightarrow & C _ { r } ( E ) & \longrightarrow & 0\\ \downarrow && \downarrow\Phi _ { 0 } && \downarrow\Phi _ { r } && \\ \wedge ^ { r - 1 } T ^ { * } \otimes h _ { 1 } & \xrightarrow{\delta} & \wedge ^ { r } T ^ { * } \otimes F _ { 0 } & \longrightarrow & F _ { r } & \longrightarrow & 0\\ \downarrow && \downarrow && \downarrow\\ 0 && 0 && 0 \end{array} $ } $$ in which $h _ { 1 } \subset T ^ { * } \otimes F _ { 0 }$ is defined as the image of $\sigma _ { 1 } ( \Phi )$ in the exact symbol sequence: $$ 0 \longrightarrow g _ { q + 1 } \longrightarrow S _ { q + 1 } T ^ { * } \otimes E \stackrel { \sigma _ { 1 } ( \Phi ) } { \longrightarrow } T ^ { * } \otimes F _ { 0 } $$ It follows that the Janet bundles can be defined ” as a whole” by the short exact sequences: Considering the first order Killing system ${ \mathcal { L } } ( \xi ) \omega = 0$ , adding its first prolongation $\mathcal { L } ( \boldsymbol { \xi } ) \gamma = 0$ while using $\xi _ { 2 }$ instead of $j _ { 2 } ( \xi )$ , we obtain a second order system $R _ { 2 } \subset J _ { 2 } ( T )$ . When $n = 2$ and $\omega$ is the Euclidean or Minkowskian metric, we have a Lie group of isometries with the 3 infinitesimal generators $\{ \partial _ { 1 } , \partial _ { 2 } , x ^ { 1 } \partial _ { 2 } - x ^ { 2 } \partial _ { 1 } \}$ . If we now consider the Weyl group defined by $\mathcal { L } ( \xi ) \omega = 2 A \omega$ and $\mathcal { L } ( \xi ) \gamma = 0$ , we have to add the only dilatation $x ^ { 1 } \partial _ { 1 } + x ^ { 2 } \partial _ { 2 }$ . As for the conformal system $\hat { R } _ { 2 } ~ \subset ~ J _ { 2 } ( T )$ defined by $( \mathcal { L } ( \xi ) \gamma ) _ { i j } ^ { k } = \delta _ { i } ^ { k } A _ { j } + \delta _ { j } ^ { k } A _ { i } - \omega _ { i j } \omega ^ { k r } A _ { r }$ according to \[18\], we have to add the two elations $\theta ^ { 1 } = \textstyle { \frac { 1 } { 2 } } ( ( x ^ { 1 } ) ^ { 2 } + ( x ^ { 2 } ) ^ { 2 } ) \partial _ { 1 } + x ^ { 1 } x ^ { 2 } \partial _ { 2 }$ and $\theta ^ { 2 } = { \textstyle { \frac { 1 } { 2 } } } ( ( x ^ { 1 } ) ^ { 2 } + ( x ^ { 2 } ) ^ { 2 } ) \partial _ { 2 } + x ^ { 1 } x ^ { 2 } \partial _ { 1 }$ . As we have $g _ { 3 } = 0 , \tilde { g } _ { 3 } = 0 , \hat { g } _ { 3 } = 0$ , we have the strict inclusions $R _ { 3 } \subset \tilde { R } _ { 3 } \subset \hat { R } _ { 3 } \subset J _ { 3 } ( T )$ of involutive systems with respective dimensions $3 < 4 < 6 < 2 0$ . Collecting these results, we get the following commutative fundamental diagram $I$ where the upper down arrows are monomorphisms while the lower down arrows are epimorphisms $\Phi _ { 0 } , \Phi _ { 1 } , \Phi _ { 2 }$ obtained by induction \[20, 29\]: $$ \begin{array}{rccccccccccccr} && && && 0 && 0 && 0 && \\ && && && \downarrow && \downarrow && \downarrow && \\ && 0 & \longrightarrow & \hat{\Theta} & \stackrel{j_3}{\longrightarrow} & 6 & \stackrel{D_1}{\longrightarrow} & 12 & \stackrel{D_2}{\longrightarrow} & 6 & \longrightarrow 0 & \\ && && && \downarrow && \downarrow && \downarrow && \\ && 0 & \longrightarrow & \Theta & \stackrel{j_3}{\longrightarrow} & 3 & \stackrel{D_1}{\longrightarrow} & 6 & \stackrel{D_2}{\longrightarrow} & 3 & \longrightarrow 0 & \hspace{3mm}Spencer \\ && && && \downarrow && \downarrow && \downarrow && \\ && 0 & \longrightarrow & 2 & \stackrel{j_3}{\longrightarrow} & 20 & \stackrel{D_1}{\longrightarrow} & 30 & \stackrel{D_2}{\longrightarrow} & 12 & \longrightarrow 0 & \hspace{3mm} hybrid \\ && && \parallel && \hspace{5mm}\downarrow \Phi_0 && \hspace{5mm}\downarrow \Phi_1 && \hspace{5mm}\downarrow \Phi_2 \\ 0 & \longrightarrow & \Theta & \longrightarrow & 2 & \stackrel{\cal{D}}{\longrightarrow} & 17 & \stackrel{{\cal{D}}_1}{\longrightarrow} & 24 & \stackrel{{\cal{D}}_2}{\longrightarrow} & 9 & \longrightarrow 0 & \hspace{3mm} Janet \\ && && && \downarrow && \downarrow && \downarrow && \\ 0 & \longrightarrow & \hat{\Theta} & \longrightarrow & 2 & \stackrel{\hat{\cal{D}}}{\longrightarrow} & 14 & \stackrel{\hat{\cal{D}}_1}{\longrightarrow} & 18 & \stackrel{\hat{\cal{D}}_2}{\longrightarrow} & 6 & \longrightarrow 0 \\ && && && \downarrow && \downarrow && \downarrow \\ && && && 0 && 0 && 0 \end{array} $$ It follows that “Spencer and Janet play at see-saw”, the dimension of each Janet bundle being decreased by the same amount as the dimension of the corresponding Spencer bundle is increased. The Poincaré sequence for the exterior derivative $d$ is $\wedge ^ { 0 } T ^ { * } \stackrel { d } { \longrightarrow } \wedge ^ { 1 } T ^ { * } \stackrel { d } { \longrightarrow } \wedge ^ { 2 } T ^ { * } \longrightarrow 0$ but it is only at the end of the paper that we shall understand the link with Maxwell equations when $n = 4$ . In Special relativity, though surprising it may look like, the above example with $n = 2$ fits with Lorentz transformations if one is using the ” hyperbolic” notations $s ( h ( \phi ) , c h ( \phi ) , t h ( \phi ) = s h ( \phi ) / c h ( \phi )$ with $x ^ { 1 } = x , x ^ { 2 } = c t$ and $d s ^ { 2 } ~ = ~ ( d x ^ { 1 } ) ^ { 2 } - ( d x ^ { 2 } ) ^ { 2 }$ . Indeed, setting $t h ( \phi ) ~ = ~ u / c$ and $t h ( \psi ) ~ = ~ v$ among dimensionless quantities, the Lorentz transformation is now described by the formulas: $\bar { x } ^ { 1 } = c h ( \phi ) x ^ { 1 } - s h ( \phi ) x ^ { 2 } , \bar { x } ^ { 2 } = - s h ( \phi ) x ^ { 1 } + c h ( \phi ) x ^ { 2 }$ . We obtain thus for the composition of speeds $t h ( \phi + \psi ) = ( ( u / c ) + ( v / c ) ) / ( 1 + ( u / c ) ( v / c ) )$ without the need of any ” gedanken experiment” on light signals. A similar result can be obtained with the ordinary ” tangent” for the composition of rotations when using the Euclidean metric $( d s ) ^ { 2 } = ( d x ^ { 1 } ) ^ { 2 } + ( d x ^ { 2 } ) ^ { 2 }$ for the plane $( x ^ { 1 } , x ^ { 2 } )$ . # MOTIVATING EXAMPLES In all the following examples we shall study linear systems of ordinary differential (OD) or partial differential (PD) equations with $m$ dependent variables, $n$ independent variables, order $q$ , Lie groups of dimension $p$ , Lie groups of transformations of a manifold $X$ with tangent bundle $T$ , cotangent bundle $T ^ { * }$ , made by a Lie group $G$ acting on $X$ with a graph $X \times G \to X \times X :$ $( x , a ) \to ( x , y = f ( x , a ) )$ or Lie pseudogroups with $m = n$ and geometric objects $( \omega , \gamma , \rho , \ldots )$ with perturbations $( \Omega , \Gamma , R , \ldots )$ . Vector bundles over $X$ will be denoted by $( E , F , \ldots )$ and their sections will be denoted by $( \xi , \eta , \zeta , \dots )$ . Symmetric covariant tensors will be denoted by $S _ { q } T ^ { * }$ while $r$ -forms will be sections of $\wedge ^ { r } T ^ { * }$ . Whenever needed, we shall introduce the adjoint vector bundle $a d ( E ) = \wedge ^ { n } T ^ { * } \otimes E ^ { * }$ in which $E ^ { * }$ is obtained from $E$ by inverting the local transition maps, exactly like $T ^ { * }$ is obtained from $T$ . The jet bundle of $E$ will be denoted by $J _ { q } ( E )$ as usual with an $q$ injective operator $j _ { q } : E \to J _ { q } ( E ) : \xi \to j _ { q } ( \xi ) = ( \xi , \partial _ { i } \xi , \partial _ { i j } \xi , . . . )$ . Finally, if $K$ is differential field with derivations $\partial _ { i }$ and $d _ { i }$ are commuting formal derivations with $d _ { i } \mid K = \partial _ { i }$ , we shall introduce the non-commutative ring $D = K [ d _ { 1 } , . . . , d _ { n } ]$ of linear differential operators with coefficients in $K$ and we have $d _ { i } a = a d _ { i } + \partial _ { i } a$ in the operator sense. Example 2.1: With $m = 1 , n = 2 , q = 2$ , let us consider the second order system $R _ { 2 } \subset J _ { 2 } ( E )$ written $d _ { 2 2 } \xi - b x ^ { 2 } d _ { 1 } \xi = \eta ^ { 2 }$ , $d _ { 1 2 } \xi - a d _ { 1 1 } \xi = \eta ^ { 1 }$ or $\mathcal { D } \xi = \eta$ with two constant parameter $( a , b )$ and ground differential field $K = \mathbb { Q } ( a , b ) ( x )$ . Three different situations may exist: $\bullet$ If $a \ = \ 0 , b \ = \ 0$ , the reader will check at once the existence of the single first order CC $d _ { 2 } \eta ^ { 1 } - d _ { 1 } \eta _ { 2 } = 0$ as an operator $\mathcal { D } _ { 1 } \eta = 0$ but the second order operator $a d ( \mathcal { D } )$ does not generate the CC of $a d ( \mathcal { D } _ { 1 } )$ which are generated by a single first order CC. $\bullet$ If $a = 1 , b = 0$ , the reader will check at once the existence of the single second order CC $d _ { 2 2 } \eta ^ { 1 } - ( d _ { 1 2 } - d _ { 1 1 } ) \eta ^ { 2 } = 0$ and that $a d ( \mathcal { D } )$ does now generate the CC of $a d ( { \mathcal { D } _ { 1 } } )$ . $\bullet$ It remains to study the case $a = 1 , b = 1$ that will bring surprises. Indeed, differentiating once, we obtain the third order system $R _ { 3 } \subset J _ { 3 } ( E )$ with corresponding Janet tabular: $$ \begin{array} { l c l } { { d _ { 2 2 2 } \xi - x ^ { 2 } d _ { 1 2 } \xi - d _ { 1 } \xi } } & { { = } } & { { d _ { 2 } \eta ^ { 2 } } } \\ { { d _ { 1 2 2 } \xi - x ^ { 2 } d _ { 1 1 } \xi } } & { { = } } & { { d _ { 1 } \eta ^ { 2 } } } \\ { { d _ { 1 1 2 } \xi - x ^ { 2 } d _ { 1 1 } \xi } } & { { = } } & { { d _ { 1 } \eta ^ { 2 } - d _ { 2 } \eta ^ { 1 } } } \\ { { d _ { 1 1 1 } \xi - x ^ { 2 } d _ { 1 1 } \xi } } & { { = } } & { { d _ { 1 } \eta ^ { 2 } - d _ { 2 } \eta ^ { 1 } - d _ { 1 } \eta ^ { 1 } } } \\ { { d _ { 2 2 } \xi - x ^ { 2 } d _ { 1 } \xi } } & { { = } } & { { \eta ^ { 2 } } } \\ { { d _ { 1 2 } \xi - d _ { 1 1 } \xi } } & { { = } } & { { \eta ^ { 1 } } } \end{array} $$ Though the symbol $g _ { 3 } = 0$ defined by $\xi _ { 2 2 2 } = 0 , \xi _ { 1 2 2 } = 0 , \xi _ { 1 1 2 } = 0 , \xi _ { 1 1 1 } = 0 )$ is trivially involutive, this system is not even formally integrable (FI). Indeed, trying all the usual crossed derivatives, we verac $R _ { 2 } ^ { ( 2 ) } \subset R _ { 2 } ^ { ( 1 ) } = R _ { 2 } \subset J _ { 2 } ( E )$ with respective dimension $3 < 4 = 4 < 6$ . After a few tricky substitutions and eliminations, we obtain the totally new second order PD equation: $$ A \equiv d _ { 1 1 } \xi = d _ { 2 2 } \eta ^ { 1 } - d _ { 1 2 } \eta ^ { 2 } + d _ { 1 1 } \eta ^ { 2 } - x ^ { 2 } d _ { 1 } \eta ^ { 1 } \in j _ { 2 } ( \eta ) $$ The hard step is to look for generating CC in the form of an operator $\mathcal { D } _ { 1 } \eta = \zeta$ . We obtain the commutative and exact diagram with $d i m ( R _ { 3 } ) = 4 , d i m ( R _ { 4 } ) = 3$ : $$ \begin{array}{ccccccccccc} & & & & 0 & & 0 & & & & \\ & & & & \downarrow & & \downarrow & & & & \\ & & 0 & \longrightarrow & S_4 T^* \otimes E & \longrightarrow & S_2 T^* \otimes F_0 & \to & h_2 & \to & 0 \\ & & \downarrow & & \downarrow & & \downarrow & & \downarrow & & \\ 0 & \longrightarrow & R_4 & \longrightarrow & J_4(E) & \longrightarrow & J_2(F_0) & \to & Q_2 & \to & 0 \\ & & \downarrow & & \downarrow & & \downarrow & & \downarrow & & \\ 0 & \to & R_3 & \to & J_3(E) & \to & J_1(F_0) & \to & Q_1 & \to & 0 \\ & & & & \downarrow & & \downarrow & & \downarrow & & \\ & & & & 0 & & 0 & & 0 & & \end{array} $$ and the long exact connecting sequence: $0 \longrightarrow R_4 \longrightarrow R_3 \longrightarrow h_2 \longrightarrow Q_2 \longrightarrow Q_1 \longrightarrow 0$. It follows that $Q _ { 1 } \simeq Q _ { 2 } = 0$ because $d i m ( h _ { 2 } ) = 1$ and there cannot exist any first or second order CC. We may start afresh with the new system $R _ { 2 } ^ { \prime } = R _ { 2 } ^ { ( 2 ) } \subset R _ { 2 } \subset J _ { 2 } ( E )$ which is involutive with symbol $g _ { 2 } ^ { \prime } = 0$ . $$ \begin{array}{rcl} d _ { 2 2 } \xi - x ^ { 2 } d _ { 1 } \xi & = & \eta ^ { 2 } \\ d _ { 1 2 } \xi & = & A + \eta ^ { 1 } \\ d _ { 1 1 } \xi & = & A \end{array} \qquad \fbox{$\begin{array}{cc} 1 & 2 \\ 1 & \bullet \\ 1 & \bullet \end{array}$} $$ We obtain therefore at once the Fundamental Diagram $I$ for $R _ { 2 } ^ { \prime }$ : $$ \begin{array}{ccccccccccccc} & & & & & & 0 & & 0 & & 0 & & \\ & & & & & & \downarrow & & \downarrow & & \downarrow & & \\ & & 0 & \rightarrow & \Theta & \stackrel{j_2}{\rightarrow } & 3 & \stackrel{d_1}{\rightarrow} & 6 & \stackrel{d_2}{\rightarrow} & 3 & \rightarrow & 0 \\ & & & & & & \downarrow & & \downarrow & &\parallel & & \\ & & 0 & \rightarrow & 1 & \stackrel{j_2}{\rightarrow} & 6 & \stackrel{d_1}{\rightarrow} &8 & \stackrel{d_2}{\rightarrow} &3 & \rightarrow & 0 \\ & & & & \parallel & & \downarrow & & \downarrow & & \downarrow & & \\ 0 & \rightarrow & \Theta & \rightarrow & 1 & \stackrel{D}{\rightarrow} & 3 & \stackrel{{D}_1}{\rightarrow} &2 & \rightarrow & 0 & & \\ & & & & & & \downarrow & & \downarrow & & & & \\ & & & & & & 0 & & 0 & & & & \end{array} $$ It finally remains to find out the generating CC for the initial second order operator $\mathcal { D } \xi = \eta$ which is neither FI nor involutive. Checking the two dots separately we have the two third order (!) CC: $$ B _ { 1 } \equiv d _ { 1 } A - x ^ { 2 } A - d _ { 1 } \eta ^ { 2 } + d _ { 2 } \eta ^ { 1 } + d _ { 1 } \eta ^ { 1 } = 0 , \qquad B _ { 2 } \equiv d _ { 2 } A - x ^ { 2 } A - d _ { 1 } \eta ^ { 2 } + d _ { 2 } \eta ^ { 1 } = 0 $$ $$ \begin{aligned} d_2 A - d_1 A - d_1 \eta^1 &\equiv d_{222} \eta^1 - d_{122} \eta^2 - d_{122} \eta^1 + 2d_{112} \eta^2 \\ &\quad - d_{111} \eta^2 - x^2 d_{12} \eta^1 - x^2 d_{11} \eta^1 = 0 \end{aligned} $$ Exactly like in \[13\], we now provide the link existing between these two third order CC and the Spencer operator. Indeed, using $R _ { 3 }$ we obtain at once: $$ \begin{array} { r c l c l } { { d _ { 1 } \xi _ { 1 1 } - \xi _ { 1 1 1 } } } & { { = } } & { { d _ { 1 } A - x ^ { 2 } A - d _ { 1 } \eta ^ { 2 } + d _ { 2 } \eta ^ { 1 } + d _ { 1 } \eta ^ { 2 } } } & { { = } } & { { B _ { 1 } } } \\ { { d _ { 2 } \xi _ { 1 1 } - \xi _ { 1 1 2 } } } & { { = } } & { { d _ { 2 } A - x ^ { 2 } A - d _ { 1 } \eta ^ { 2 } + d _ { 2 } \eta ^ { 1 } } } & { { = } } & { { B _ { 2 } } } \end{array} $$ a result that we shall obtain after one more prolongation by the long exact sequence: $$ \begin{array}{c@{\qquad}c} 0 \rightarrow R_5 \rightarrow J_5(E) \rightarrow J_3(F_0) \rightarrow Q_3 \rightarrow 0 & 0 \rightarrow 3 \rightarrow 21 \rightarrow 20 \rightarrow 2 \rightarrow 0 \end{array} $$ We may thus define $F _ { 1 } = Q _ { 3 }$ with $d i m ( F _ { 1 } ) = 2$ and define similarly $F _ { 2 }$ with $d i m ( F _ { 2 } ) = 1$ by the long exact sequence: $$ \begin{aligned} 0 \rightarrow R_6 \rightarrow J_6(E) \rightarrow J_4(F_0) \rightarrow J_1(F_1) \rightarrow F_2 \rightarrow 0 \\ 0 \rightarrow 3 \rightarrow 28 \rightarrow 30 \rightarrow 6 \rightarrow 1 \rightarrow 0 \end{aligned} $$ We have indeed $d _ { 1 } B _ { 2 } - d _ { 2 } B _ { 1 } = 0$ and the exact differential sequence which is not a Janet sequence: $$ \fbox{ $ 0 \longrightarrow \Theta \rightarrow 1 \underset 2{\stackrel{{\cal{D}}}{\longrightarrow}} 2 \underset 3{\stackrel{{\cal{D}}_1}{\longrightarrow}} 2 \underset 1{\stackrel{{\cal{D}}_2}{\longrightarrow}} 1 \longrightarrow 0 $ } $$ More generally, we have the long exact sequences $\forall r \geq 0$ . $$ 0 \to 3 \to J_{r+6}(E) \to J_{r+4}(F_0) \to J_{r+1}(F_1) \to J_r(F_2) \to 0 $$ Using finally the basis $\{ \theta_{\tau}(x) \mid 1 \le \tau \le 3 \} = \{ 1, x^2, x^1 + \frac{1}{6}(x^2)^3 \}$ for the vector space $\nu$ over the constants, the Spencer sequence of the Fundamental Diagram $I$ is the tensor product by $\nu$ of the Poincaré sequence $\wedge ^ { 0 } T ^ { * } \xrightarrow { d } \wedge ^ { 1 } T ^ { * } \xrightarrow { d } \wedge ^ { 2 } T ^ { * } \longrightarrow 0$ for the exterior derivative when $n = 2$ . We shall discover later on that the sitution of the present example with two parameters $( a , b )$ and three different cases is exactly similar to the one that will be provided by the M, S and K metrics with 2 parameters $( m , a )$ and has thus nothing to do with any GR framework. Example 2.2: With again $m = 1 , n = 2 , q = 2 , K = \mathbb { Q } ( x ^ { 2 } )$ , the system $R _ { 2 } \subset J _ { 2 } ( E )$ defined by $d _ { 2 2 } \xi - x ^ { 2 } d _ { 1 } \xi = 0 , d _ { 1 2 } \xi - \xi = 0$ has the only solution $\xi = 0$ because $R _ { 2 } ^ { ( 3 ) } = 0$ with $r = 0 , s = 3$ and only the central hybrid sequence for $j _ { 2 }$ is left: $$ ~ 0 \longrightarrow 1 \xrightarrow { j _ { 2 } } 6 \xrightarrow { D _ { 1 } } 8 \xrightarrow { D _ { 2 } } 3 \xrightarrow { \qquad } 0 $$ Example 2.3: (Macaulay) With $m = 1 , n = 3 , q = 2$ and $P , Q , R \in { \cal D } = \mathbb { Q } [ d _ { 1 } , d _ { 2 } , d _ { 3 } ]$ , let us consider the linear homogeneous second order system of PD equations $R _ { 2 } \subset J _ { 2 } ( E )$ defined by: $$ { \cal P } \xi \equiv d _ { 3 3 } \xi = 0 , ~ Q \xi \equiv d _ { 2 3 } \xi - d _ { 1 1 } \xi = 0 , ~ R \xi \equiv d _ { 2 2 } \xi = 0 $$ but the reader may treat as well the system $( d _ { 3 3 } \xi - d _ { 1 1 } \xi = 0 , d _ { 2 3 } \xi = 0 , d _ { 2 2 } \xi - d _ { 1 1 } \xi = 0 )$ . Of course, this system is FI because it is homogeneous but we let the reader check through the Janet tabular that $g _ { 2 }$ is not involutive though the coordinate system is surely $\delta$ -regular because be have full class 3 and full class 2. All the third order jets vanish but $y _ { 1 2 3 } - y _ { 1 1 1 } = 0$ leading to $d i m ( g _ { 3 } ) = 1 \Rightarrow d i m ( R _ { 3 } ) = 8$ . Finally $g _ { 4 } = 0 \Rightarrow d i m ( R _ { 4 } ) = 8$ and we could believe that we do not need any PP procedure as $R _ { 4 }$ is an ivolutive system because $g _ { 4 } = 0$ is trivially involutive and $R _ { 2 }$ is finite type like the Killing system. It is important to notice that the knowledge of the first second order operator does not provide any way to obtain the third without passing through the second, contrary to the situation existing in the Janet sequence. Such a procedure is rather “experimental” and must be coherent with a theorem saying that the order of generating CC is one plus the number of prolongations needed to reach a 2-acyclic symbol, that is $g _ { 3 }$ must be 2-acyclic \[20\]. Equivalently the $\delta$ -sequence: $$ 0 \to { \textstyle \bigwedge } ^ { 2 } T ^ { * } \otimes g _ { 3 } \stackrel { \delta } { \to } { \textstyle \bigwedge } ^ { 3 } T ^ { * } \otimes g _ { 2 } \to 0 $$ must be exact. We let the reader prove that the corresponding $3 \times 3$ matrix has maximum rank. We recall the dimensions of the following jet bundles: $$ \begin{array} { c c c c c c c c c c c c c c c } { q } & { \to } & { 0 } & { 1 } & { 2 } & { 3 } & { 4 } & { 5 } & { 6 } & { 7 } \\ { S _ { q } T ^ { * } } & { \to } & { 1 } & { 3 } & { 6 } & { 1 0 } & { 1 5 } & { 2 1 } & { 2 8 } & { 3 6 } \\ { J _ { q } ( E ) } & { \to } & { 1 } & { 4 } & { 1 0 } & { 2 0 } & { 3 5 } & { 5 6 } & { 8 4 } & { 1 2 0 } \end{array} $$ and the commutative and exact diagram allowing to construct the Spencer bundles $C _ { r } \subset C _ { r } ( E )$ and the Janet bundles $F _ { r }$ for $r = 0 , 1 , . . . , n$ with $F _ { 0 } = J _ { q } ( E ) / R _ { q }$ , showing that we have indeed with $q = 4$ : $$ C _ { r } = \wedge ^ { r } T ^ { * } \otimes R _ { q } $$ $$ C _ { r } ( E ) = \wedge ^ { r } T ^ { * } \otimes J _ { q } ( E ) / \delta ( \wedge ^ { r - 1 } T ^ { * } \otimes S _ { q + 1 } T ^ { * } \otimes E $$ $$ F _ { r } = \wedge ^ { r } T ^ { * } \otimes J _ { q } ( E ) / ( \wedge ^ { r } T ^ { * } \otimes R _ { q } + \delta ( \wedge ^ { r - 1 } T ^ { * } \otimes S _ { q + 1 } T ^ { \otimes } E ) ) $$ When $R _ { q } \subset J _ { q } ( E )$ is involutive, that is formally integrable (FI) with an involutive symbol $g _ { q }$ , then these three differential sequences are formally exact on the jet level and, in the Spencer sequence: $$ 0 \longrightarrow \Theta \stackrel { j _ { q } } { \longrightarrow } C _ { 0 } \stackrel { D _ { 1 } } { \longrightarrow } C _ { 1 } \stackrel { D _ { 2 } } { \longrightarrow } \dots \stackrel { D _ { n } } { \longrightarrow } C _ { n } \longrightarrow 0 $$ the first order involutive operators $D _ { 1 } , D _ { 2 } , . . . , D _ { n }$ are induced by the Spencer operator $d : R _ { q + 1 } \to$ $T ^ { * } \otimes R _ { q }$ already considered that can be extended to $d : \wedge ^ { r } T ^ { * } \otimes R _ { q + 1 } \to \wedge ^ { r + 1 } \otimes R _ { q }$ . A similar condition is also valid for the Janet sequence: $$ 0 \longrightarrow \Theta \longrightarrow E \stackrel { \mathcal { D } } { \longrightarrow } F _ { 0 } \stackrel { \mathcal { D } _ { 1 } } { \longrightarrow } F _ { 1 } \stackrel { \mathcal { D } _ { 2 } } { \longrightarrow } \dots \stackrel { \mathcal { D } _ { n } } { \longrightarrow } F _ { n } \longrightarrow 0 $$ which can be thus constructed “as a whole” from the previous extension of the Spencer operator (See \[29\], p 185 and 391 for diagrams missing in \[14, 15\] !). However, this result is still not known and not even acknowledged today in mathematicl physics, particularly in general relativity which is never using the Spencer $\delta$ -cohomology in order to define the Riemann or Bianchi operators. The study of the present Macaulay example will be sufficient in order to justify our comment. First of all, as $g _ { 2 }$ is not 2-acyclic and the coeficients are constant, the second order CC are $Q w - R v = 0 , \ R u - P w = 0 , \ P v - Q u = 0$ and the simplest resolution is thus: $$ \begin{array}{r} 0 \longrightarrow \Theta \longrightarrow 1 \underset{2}{\stackrel{\mathcal{D}}{\longrightarrow}} 3 \underset{2}{\stackrel{\mathcal{D}_1}{\longrightarrow}} 3 \underset{2}{\stackrel{\mathcal{D}_2}{\longrightarrow}} 1 \longrightarrow 0 \end{array} $$ Secondly, as the first prolongation of $R _ { 2 }$ becoming involutive is $R _ { 4 }$ because $g _ { 4 } = 0$ , an idea could be to start with the system $R _ { 3 } \subset J _ { 3 } ( E )$ but we have proved in (\[32\], Example 3.14, p 119 to 126) that the simplest formally exact sequence that could be formally exact is quite far from being a Janet sequence as it is: $$ 0 \longrightarrow \Theta \longrightarrow 1 \underset { 3 } { \longrightarrow } 1 2 \underset { 1 } { \longrightarrow } 2 1 \underset { 2 } { \longrightarrow } 4 6 \underset { 1 } { \longrightarrow } 7 2 \underset { 1 } { \longrightarrow } 4 8 \underset { 1 } { \longrightarrow } 1 2 \longrightarrow 0 $$ Indeed, the Euler-Poincaré characteristic is $1 - 1 2 + 2 1 - 4 6 + 7 2 - 4 8 + 1 2 = 0$ but we notice that the orders of the successive operators may vary up and down. It remains to work out the Janet and Spencer sequences in the fundamental diagram I : $$ \begin{array}{rccccccccccccccl} &&&&& 0 &&0&&0& & 0 & & \\ &&&&& \downarrow && \downarrow && \downarrow & & \downarrow & & \\ & 0& \rightarrow & \Theta & \stackrel{j_4}{\longrightarrow} & 8 &\stackrel{d_1}{\longrightarrow} & 24 & \stackrel{d_2}{\longrightarrow} & 24 &\stackrel{d_3}{\longrightarrow} & 8 & \rightarrow & 0 \\ &&&&& \downarrow & & \downarrow & & \downarrow & & \downarrow & & \\ & 0 & \rightarrow & 1 & \stackrel{j_4}{\longrightarrow} & 35 & \stackrel{d_1}{\longrightarrow} & 84 &\stackrel{d_2}{\longrightarrow} & 70 & \stackrel{d_3}{\longrightarrow} & 20 & \rightarrow & 0 \\ & & & \parallel & & \hspace{5mm}\downarrow {\Phi}_0 & &\hspace{5mm} \downarrow {\Phi}_1 & & \hspace{5mm} \downarrow {\Phi}_2& & \hspace{5mm} \downarrow {\Phi}_3 & \\ 0 \rightarrow & \Theta &\rightarrow & 1 & \stackrel{D}{\longrightarrow} & 27 & \stackrel{D_1}{\longrightarrow} & 60 & \stackrel{D_2}{\longrightarrow} & 46 & \stackrel{D_3}{\longrightarrow} & 12 & \rightarrow & 0 \\ &&&&& \downarrow & & \downarrow & & \downarrow & & \downarrow & & \\ &&&&& 0 && 0 && 0 & & 0 & & \end{array} $$ As there is no group background, it is nevertheless not evident that the above Spencer sequence is isomorphic to the tensor product of the Poincaré sequence for the exterior derivative by a vector space $\nu$ of dimension 8. For this, we notice that each solution is a linear combination of polynomials of degree 3 at most in $\mathbb { Q } [ x ^ { 1 } , x ^ { 2 } , x ^ { 3 } ]$ . After tricky substitutions, we just need to choose the basis: $$ \left. \theta _ { \tau } \mid 1 \leq \tau \leq 8 \right. = \{ 1 , x ^ { 1 } , x ^ { 2 } , x ^ { 3 } , x ^ { 1 } x ^ { 2 } , x ^ { 1 } x ^ { 3 } , { \frac { 1 } { 2 } } ( x ^ { 1 } ) ^ { 2 } + x ^ { 2 } x ^ { 3 } , { \frac { 1 } { 6 } } ( x ^ { 1 } ) ^ { 3 } + x ^ { 1 } x ^ { 2 } x ^ { 3 } \} $$ All the operators are of order 1 but $j _ { 4 }$ and $\mathcal { D }$ which are of order 4. Example 2.4:(Contact transformations revisited) Whith $m = 3 , n = 3$ , let $\alpha = d x ^ { 1 } - x ^ { 3 } d x ^ { 2 } \in T ^ { * }$ be the so-called contact 1-form. The Lie pseudogroup of contact transformations $\Gamma = \{ y = f ( x ) \mid$ $d y ^ { 1 } - y ^ { 3 } d y ^ { 2 } = a ( x ) ( d x ^ { 1 } - x ^ { 3 } d x ^ { 2 } ) \}$ preserves $\alpha$ up to a factor $a ( x )$ . Eliminating this factor among the three infinitesimal Lie equations, we obtain two PD equations but this system is neither involutive nor even FI. The system of infinitesimal Lie equations defining the infinitesimal contact transformations $\Theta \subset T$ is obtained by eliminating the factor $\rho ( x )$ in the equations ${ \mathcal { L } } ( \xi ) \alpha = \rho \alpha$ where $\mathcal { L }$ is the standard Lie derivative \[26\]. This system is thus only generated by $\eta ^ { \mathrm { { I } } }$ and $\eta ^ { 2 }$ below but is not involutive and one has to introduce $\eta ^ { 3 }$ in order to obtain the following involutive system $R _ { 1 } \subset J _ { 1 } ( T )$ with two equations of class 3 and one equation of class 2, a result leading to $d i m ( g _ { 1 } ) = 6$ , $d i m ( g _ { 2 } ) = 1 0$ and $d i m ( g _ { 3 } ) = 1 5$ . $$ \fbox{$ \begin{array}{rclc} d_3 \xi^3 + d_2 \xi^2 + 2x^3 d_1 \xi^2 - d_1 \xi^1 & = & \eta^3 & \\ d_3 \xi^1 - x^3 d_3 \xi^2 & = & \eta^2 & \left| \begin{array}{ccc} 1 & 2 & 3 \\ 1 & 2 & 3 \\ 1 & 2 & \bullet \end{array} \right| \\ d_2 \xi^1 - x^3 d_2 \xi^2 + x^3 d_1 \xi^1 - (x^3)^2 d_1 \xi^2 - \xi^3 & = & \eta^1 & \end{array} $} $$ We have one first order CC: $$ d _ { 3 } { \eta } ^ { 1 } - d _ { 2 } { \eta } ^ { 2 } - x ^ { 3 } d _ { 1 } { \eta } ^ { 2 } + { \eta } ^ { 3 } = \zeta $$ showing the origin of $\eta ^ { 3 }$ and finally get the Janet sequence: $$ \begin{array} { r } { \left| \begin{array} { l } { 0 \to \Theta \to 3 \xrightarrow { \mathcal { D } } 3 \xrightarrow { \mathcal { D } _ { 1 } } 1 \to 0 } \end{array} \right| } \end{array} $$ Let us now consider the new pseudogroup $\{ \Gamma ^ { \prime } \subset \Gamma \mid d y ^ { 1 } - y ^ { 3 } d y ^ { 2 } = d x ^ { 1 } - x ^ { 3 } d x ^ { 1 } \}$ preserving the contact form $\alpha$ , thus also $\beta = d \alpha = d x ^ { 2 } \wedge d x ^ { 3 }$ and $\gamma = \alpha \wedge \beta = d x ^ { 1 } \wedge d x ^ { 2 } \wedge d x ^ { 3 }$ which is the so-called volume form. The corresponding infinitesimal Lie equations are successively: $$ d _ { 1 } \xi ^ { 1 } - x ^ { 3 } d _ { 1 } \xi ^ { 2 } 0 , ~ d _ { 2 } \xi ^ { 1 } - x ^ { 3 } d _ { 2 } \xi ^ { 2 } - \xi ^ { 3 } = 0 , ~ d _ { 3 } \xi ^ { 1 } - x ^ { 3 } d _ { 3 } \xi ^ { 2 } = 0 $$ $$ d _ { 1 } \xi ^ { 2 } = 0 , \ d _ { 1 } \xi ^ { 3 } = 0 , \ d _ { 2 } \xi ^ { 2 } + d _ { 3 } \xi ^ { 3 } = 0 $$ $$ d _ { 1 } \xi ^ { 1 } + d _ { 2 } \xi ^ { 2 } + d _ { 3 } \xi ^ { 3 } = 0 $$ Doing the permutation $( 1 2 3 ) ( 3 1 2 )$ , we obtain the first order system $R _ { 1 } ^ { \prime } \subset J _ { 1 } ( T )$ : $$ \begin{array}{rcl} d_3 \xi^3 & = & 0 \\ d_3 \xi^2 & = & 0 \\ d_3 \xi^1 & = & 0 \\ d_2 \xi^3 + d_1 \xi^2 & = & 0 \\ d_2 \xi^1 - x^2 d_2 \xi^2 & = & 0 \\ d_1 \xi^1 - x^2 d_1 \xi^2 - \xi^3 & = & 0 \end{array} \qquad \left| \begin{array}{ccc} 1 & 2 & 3 \\ 1 & 2 & 3 \\ 1 & 2 & 3 \\ 1 & 2 & \bullet \\ 1 & \bullet & \bullet \end{array} \right| $$ and we have the strict inclusions $R _ { 1 } ^ { \prime } \subset R _ { 1 } \subset J _ { 1 } ( T )$ with $6 < 9 < 1 2$ . We have the new Janet sequence with $3 - 6 + 4 - 1 = 0$ : $$ \begin{array} { r } { 0 \to \Theta ^ { \prime } \to 3 \stackrel { \mathcal { D } ^ { \prime } } { \to } 6 \stackrel { \mathcal { D } _ { 1 } ^ { \prime } } { \to } 4 \stackrel { \mathcal { D } _ { 2 } ^ { \prime } } { \to } 1 \to 0 } \end{array} $$ Contrary to other examples, in the present situation it is easy to work out the Janet sequences but much more difficult to work out the corresponding Spencer sequences. In fact, we obtain the following Fundamental diagram $I$ in which the central hybrid sequence, which is at the same time a Janet sequence for $j _ { 1 } i$ and a Spencer sequence for $J _ { 2 } ( T ) \subset J _ { 1 } ( J _ { 1 } ( T ) )$ , is the same for $\Theta ^ { \prime } \subset \Theta$ and all the operators involved are first order: $$ \begin{array}{ccccccccccccccc} & & & & & & 0 & & 0 & & 0 & & 0 & & \\ & & & & & & \downarrow & & \downarrow & & \downarrow & &\downarrow & & \\ & & 0 & \longrightarrow &{ \Theta }'& \stackrel{j_1}{\longrightarrow} & 6 & \stackrel{d_1}{\longrightarrow} & 14 & \stackrel{d_2}{\longrightarrow} & 11 & \stackrel{d_3}{\longrightarrow} & 3 & \longrightarrow & 0 \\ & & & & & & \downarrow & & \downarrow & & \downarrow & & \parallel & & \\ & & 0 & \longrightarrow & 3 & \stackrel{j_1}{\longrightarrow } &12 & \stackrel{d_1 }{\longrightarrow } & 18 & \stackrel{d_2}{\longrightarrow} & 12 & \stackrel{d_3}{\longrightarrow} & 3 & \rightarrow & 0 \\ & & & & \parallel & & \downarrow & & \downarrow & & \downarrow & &\downarrow & & \\ 0 & \longrightarrow & { \Theta}' & \longrightarrow & 3 & \stackrel{ D'}{\longrightarrow} & 6 & \stackrel{D'_1}{\longrightarrow} & 4 & \stackrel{ D'_2}{\longrightarrow} & 1 & \rightarrow & 0 & & \\ & & & & & & \downarrow & &\downarrow & & \downarrow & & & & \\ & & & & & & 0 & & 0 & & 0 & & & & \end{array} $$ In actual practice, when $g _ { q }$ is involutive, its knowledge allows to determines the Hilbert function $d i m ( g _ { q + r } )$ and finally $d i m ( R _ { q + r } )$ . Then $d i m ( C _ { r } )$ may be computed by induction on $r$ . In this example with $q = 1$ , one has the commutative and exact diagram with $d i m ( C _ { 0 } ) = d i m ( R _ { 1 } ) = 9 \Rightarrow$ $d i m ( C _ { 1 } ) = 1 7 \Rightarrow d i m ( C _ { 2 } ) = 1 2$ and so on: $$ \begin{array} { c c c c c c c c c c c c c c c c c } { { 0 } } & { { \to } } & { { R _ { q + 2 } } } & { { \to } } & { { J _ { 2 } ( R _ { q } ) } } & { { \to } } & { { J _ { 1 } ( C _ { 1 } ) } } & { { \to } } & { { C _ { 2 } } } & { { \to } } & { { 0 } } \\ { { } } & { { } } & { { \downarrow } } & { { } } & { { \downarrow } } & { { } } & { { \downarrow } } & { { } } & { { \downarrow } } & { { } } & { { } } \\ { { 0 } } & { { \to } } & { { R _ { q + 1 } } } & { { \to } } & { { J _ { 1 } ( R _ { q } ) } } & { { \to } } & { { C _ { 1 } } } & { { \to } } & { { 0 } } & { { } } & { { } } \end{array} $$ Unhappily, Spencer and collaborators, being unaware on one side of the computer algebra works of Janet and Gröbner or of the work of Vessiot on geometric objects, were unable to compute explicit examples. In particular, anybody can check that the examples presented in the introduction of \[27\] have no relation at all with the remaining of the book. It is thus sometimes easier to exhibit the Janet sequences than the Spencer sequences. By chance the situation will be opposite when we shall deal with general relativity and black holes. Example 2.5: (Bose conjecture) Last but not least, the following example is the most difficult of this paper as it is exhibiting a mixture of second and third order CC \[4\]. With $m = n = 3$ , the trivial differential field $K = \mathbb { Q }$ and the ring $D = K [ d _ { 1 } , d _ { 2 } , d _ { 3 } ]$ of differential operators, let us consider the $D$ module $M$ defined by the second order operator: $$ \mathcal { D } : ( \xi ^ { 1 } , \xi ^ { 2 } , \xi ^ { 3 } ) \to ( \xi _ { 3 3 } ^ { 3 } - \xi _ { 1 } ^ { 1 } = \eta ^ { 1 } , \xi _ { 2 3 } ^ { 3 } - \xi _ { 1 } ^ { 2 } - \xi ^ { 3 } = \eta ^ { 2 } ) $$ A first question is to determine its torsion submodule $t ( M )$ or, equivalently, to know whether $\mathcal { D }$ can be parametrized by a certain operator $\mathcal { D } _ { - 1 }$ through a certain number of arbitrary potential functions $\phi$ . Using double duality, we may construct the adjoint operator $a d ( \mathcal { D } )$ by multiplying these equations respectively by $( \lambda ^ { 1 } , \lambda ^ { 2 } )$ , adding and integrating by parts in order to get successively: $$ \begin{array} { c } { { a d ( \mathcal { D } ) : ( \lambda ^ { 1 } , \lambda ^ { 2 } ) \to ( d _ { 1 } \lambda ^ { 1 } = \mu ^ { 1 } , d _ { 1 } \lambda ^ { 2 } = \mu ^ { 2 } , d _ { 3 3 } \lambda ^ { 1 } + d _ { 2 3 } \lambda ^ { 2 } - \lambda ^ { 2 } = \mu ^ { 3 } ) } } \\ { { a d ( \mathcal { D } _ { - 1 } ) : ( \mu ^ { 1 } , \mu ^ { 2 } , \mu ^ { 3 } ) \to d _ { 3 3 } \mu ^ { 1 } + d _ { 2 3 } \mu ^ { 2 } - d _ { 1 } \mu ^ { 3 } - \mu ^ { 2 } = \nu } } \\ { { { { { \mathcal { D } _ { - 1 } : \phi \to ( d _ { 3 3 } \phi = \xi ^ { 1 } , d _ { 2 3 } \phi - \phi = \xi ^ { 2 } , d _ { 1 } \phi = \xi ^ { 3 } ) } } } \ } } \end{array} $$ There is a new second order CC for $\mathcal { D } _ { - 1 }$ , namely $\eta ^ { 3 } \equiv d _ { 3 3 } \xi ^ { 2 } - ( d _ { 2 3 } - 1 ) \xi ^ { 1 } = 0$ and we check at once that $d _ { 1 } \eta ^ { 3 } = 0$ , that is $t ( M )$ is generated by $\eta ^ { 3 }$ because we have: $$ \begin{array} { r } { ( d _ { 2 3 } - 1 ) \eta ^ { 1 } - d _ { 3 3 } \eta ^ { 2 } = ( \xi _ { 1 3 3 } ^ { 2 } + \xi _ { 3 3 } ^ { 3 } - \xi _ { 1 2 3 } ^ { 1 } ) - ( \xi _ { 3 3 } ^ { 3 } - \xi _ { 1 } ^ { 1 } ) = \xi _ { 1 3 3 } ^ { 2 } - \xi _ { 1 2 3 } ^ { 1 } + \xi _ { 1 } ^ { 1 } = d _ { 1 } \eta ^ { 3 } } \end{array} $$ It follows that $M ^ { \prime } = M / t ( M ) \simeq D$ is a free, thus torsion-free and even projective module defined by $\mathcal { D } ^ { \prime }$ which is parametrized by the injective operator $\mathcal { D } _ { - 1 }$ : $$ \mathcal { D } ^ { \prime } : ( \xi ^ { 1 } , \xi ^ { 2 } , \xi ^ { 3 } ) ( \xi _ { 3 3 } ^ { 3 } - \xi _ { 1 } ^ { 1 } = \eta ^ { 1 } , \xi _ { 2 3 } ^ { 3 } - \xi _ { 1 } ^ { 2 } - \xi ^ { 3 } = \eta ^ { 2 } , \xi _ { 3 3 } ^ { 2 } - \xi _ { 2 3 } ^ { 1 } + \xi ^ { 1 } = \eta ^ { 3 } ) $$ Indeed, we have an isomorphism $M \simeq t ( M ) \oplus M ^ { \prime }$ because the following splitting : $$ d _ { 3 3 } \phi = \xi ^ { 1 } , \ d _ { 2 3 } \phi - \phi = \xi ^ { 2 } \ \Rightarrow \ d _ { 3 } \phi = d _ { 2 } \ \xi ^ { 1 } - d _ { 3 } \xi ^ { 2 } \ \Rightarrow \ \phi = d _ { 2 2 } \xi ^ { 1 } - d _ { 2 3 } \xi ^ { 2 } - \xi ^ { 2 } $$ It is not at all evident that $M ^ { \prime }$ can be defined by only two differentially independent PD equations, the first one for $( \xi ^ { 1 } , \xi ^ { 2 } )$ while the second is providing $\xi ^ { 3 }$ : $$ \xi _ { 3 3 } ^ { 2 } - \xi _ { 2 3 } ^ { 1 } + \xi ^ { 1 } = 0 , \xi _ { 1 2 3 } ^ { 2 } - \xi _ { 1 2 2 } ^ { 1 } + \xi _ { 1 } ^ { 2 } + \xi ^ { 3 } = 0 $$ a result showing that the CC operator $\mathcal { D } ^ { \prime \prime }$ of the parametrizing operator $\mathcal { D } _ { - 1 }$ may be generated by one second order CC and one third order CC, a striking result providing the short exact sequence: $$ \begin{array} { r } { 0 \longrightarrow 1 \xrightarrow { \mathcal { D } _ { - 1 } } 3 \xrightarrow { \mathcal { D } ^ { \prime \prime } } 2 \longrightarrow 0 } \end{array} $$ Now, we may consider the new system: $$ \begin{array}{rcl} \xi _ { 3 3 } ^ { 3 } - \xi _ { 1 } ^ { 1 } & = & \eta ^ { 1 } \\ \xi _ { 3 3 } ^ { 2 } - \xi _ { 2 3 } ^ { 1 } + \xi ^ { 1 } & = & \eta ^ { 3 } \\ \xi _ { 2 3 } ^ { 3 } - \xi _ { 1 } ^ { 2 } - \xi ^ { 3 } & = & \eta ^ { 2 } \\ \xi _ { 1 3 } ^ { 2 } - \xi _ { 1 2 } ^ { 1 } + \xi _ { 3 } ^ { 3 } & = & d _ { 2 } \eta ^ { 1 } - d _ { 3 } \eta ^ { 2 } \end{array} \qquad \left| \begin{array}{ccc} 1 & 2 & 3 \\ 1 & 2 & 3 \\ 1 & 2 & \bullet \\ 1 & \times & \bullet \end{array} \right| $$ The lack of involutivity provides $d _ { 2 2 } \eta ^ { 1 } - d _ { 2 3 } \eta ^ { 2 } - \eta ^ { 2 } = \xi _ { 1 2 3 } ^ { 2 } - \xi _ { 1 2 2 } ^ { 1 } + \xi _ { 1 } ^ { 2 } + \xi ^ { 3 }$ . Conversely,tiv $\xi _ { 1 3 } ^ { 2 } - \xi _ { 1 2 } ^ { 1 } + \xi _ { 3 } ^ { 3 } = 0$ and thus $\xi _ { 1 3 3 } ^ { 2 } - \xi _ { 1 2 3 } ^ { 1 } + \xi _ { 3 3 } ^ { 3 } = 0$ However, we have $d _ { 1 } \eta ^ { 3 } = \xi _ { 1 3 3 } ^ { 2 } - \xi _ { 1 2 3 } ^ { 1 } + \xi _ { 1 } ^ { 1 } = 0$ d $\eta ^ { 1 } = \xi _ { 3 3 } ^ { 3 } - \xi _ { 1 } ^ { 1 } = 0$ . Similarly, we obtain $\xi _ { 1 2 3 } ^ { 2 } - \xi _ { 1 2 2 } ^ { 1 } + \xi _ { 2 3 } ^ { 3 } = 0$ and obtain by subtraction $\eta ^ { 2 } = \xi _ { 2 3 } ^ { 3 } - \xi _ { 1 } ^ { 2 } - \xi ^ { 3 } = 0$ . Remark 2.6: Coming back to the previous motivating examples depending on parameters, we notice that the order of the generating CC may highly depend on the parameters. Such a situation, having nothing to do with physics, is nevertheless quite similar to that of the first order Killing system $R _ { 1 } \subset J _ { 1 } ( T )$ allowing to define the first order Killing operator ${ \mathcal { D } } : T \to S _ { 2 } T ^ { * } : \xi \to { \mathcal { L } } ( \xi ) \omega = \Omega$ through the Lie derivative of a non-degenerate metric $\omega$ , for example in the case of the M, S, K metrics where a prolongation $R _ { q + r } ^ { ( s ) } = \pi _ { q + r } ^ { q + r + s } ( R _ { q + r + s } ) \subset R _ { q + r }$ must be done. Using the Lie algebra $\vec { \mathcal { G } }$ with dimension 10 for $M$ , 4 for $S$ and 2 for $K$ instead of $n$ , the Spencer sequence is always isomorphic to the tensor product of the Poincaré sequence for the exterior derivative by a Lie algebra that may have a very small dimension as we shall see. Accordingly, we claim:\ \ THE IMPORTANT OBJECT IS NOT THE METRIC BUT ITS GROUP OF INVARIANCE In particular, the FACT that third order generating CC for the Killing operator may exist has no physical meaning as nobody is knowing a way to select a best candidate among the possible explicit solutions of Einstein equations in vacuum, a mathematical result questioning the origin and existence of black holes as we shall see !. We also notice the fact that the PP procedure is highly depending on the various parameters involved, namely the only parameter $m$ for the S metric which is reduced to the M metric when $m = 0$ while the K metric depends on the two parameters $( m , a )$ and is reduced to the S metric when $a = 0$ . We study now this comment.  # DIFFERENTIAL TOOLS ## From Group Theory to Differential Sequences Let $G$ be a Lie group with coordinates $( a ^ { \rho } ) = ( a ^ { 1 } , . . . , a ^ { p } )$ acting on a manifold $X$ with a local action map $y = f ( x , a )$ . According to the second fundamental theorem of Lie, if $\theta _ { 1 } , . . . , \theta _ { p }$ are the infinitesimal generators of the effective action of a lie group $G$ on $X$ , then $[ \theta _ { \rho } , \theta _ { \sigma } ] = c _ { \rho \sigma } ^ { \tau } \theta _ { \tau }$ where the $c = ( c _ { \rho \sigma } ^ { \tau } = - c _ { \sigma \rho } ^ { \tau } )$ are the structure constants of a Lie algebra of vector felds which can be identified with $\mathcal { G } = T _ { e } ( G )$ the tangent space to $\mathcal { G }$ at the identity $e \in G$ by using the action. More generally, if $X$ is a manifold and $G$ is a lie group (not acting necessarily on $X$ ), let us consider gauging maps $a : X \to G : ( x ) \to ( a ( x ) )$ . If $x + d x$ is a point of $X$ close to $x$ , then the tangent mapping $T ( a ) : T = T ( X ) \to T ( G ) : d x \to d a = ( \partial a / \partial x ) d x$ will provide a point $a + d a$ close to $a$ on $G$ . We may bring $a$ back to $e$ on $G$ by acting on $a$ with $a ^ { - 1 }$ on the left, that is $b \to a ^ { - 1 } b , \forall b \in G$ . We get therefore a 1-form $a ^ { - 1 } d a = A \in T ^ { * } \otimes { \mathcal { G } }$ and the curvature 2-form $F = ( \partial _ { i } A _ { j } ^ { \tau } ( x ) - \partial _ { j } A _ { i } ^ { \tau } ( x ) - c _ { \rho \sigma } ^ { \tau } A _ { i } ^ { \rho } ( x ) A _ { j } ^ { \sigma } ( x ) = F _ { i j } ^ { \tau } ( x ) ) \in \wedge ^ { 2 } T ^ { * } \otimes \mathcal { G }$ in the nonlinear gauge sequence: $$ \begin{array}{ccccc} X\times G & \longrightarrow & T^*\otimes {\cal{G}} & \longrightarrow & {\wedge}^2T^*\otimes {\cal{G}} \\ a & \longrightarrow & a^{-1}da=A & \longrightarrow & dA-[A,A]=F \end{array} $$ In 1956, at the birth of gauge theory (GT), the above notations were coming from the EM potential $A$ and EM field $d A = F$ of relativistic Maxwell theory. Accordingly, $G = U ( 1 )$ (unit circle in the complex plane) $\longrightarrow \dim ( \mathcal { G } ) = 1$ was the only possibility to get a 1-form $A$ and a 2-form $F$ with vanishing structure constants $c = 0$ . Choosing now $a$ “close” to $e$ , that is $a ( x ) = e + t \lambda ( x ) + . . .$ and linearizing as usual, we obtain the linear operator $d : \wedge ^ { 0 } T ^ { * } \otimes \mathcal { G } \to \wedge ^ { 1 } T ^ { * } \otimes \mathcal { G } : ( \lambda ^ { \tau } ( x ) ) \to ( \partial _ { i } \lambda ^ { \tau } ( x ) )$ and the linear gauge sequence: $$ { \textstyle \wedge } ^ { 0 } T ^ { * } \otimes { \mathcal { G } } \ { \stackrel { d } { \longrightarrow } } \ { \wedge } ^ { 1 } T ^ { * } \otimes { \mathcal { G } } \ { \stackrel { d } { \longrightarrow } } \ { \wedge } ^ { 2 } T ^ { * } \otimes { \mathcal { G } } \ { \stackrel { d } { \longrightarrow } } \ \ldots \ { \stackrel { d } { \longrightarrow } } \ { \wedge } ^ { n } T ^ { * } \otimes { \mathcal { G } } \longrightarrow 0 $$ which is the tensor product by $\mathcal { G }$ of the Poincaré sequence for the exterior derivative. Considering now a Lagrangian on $T ^ { \ast } \otimes \mathcal { G }$ , that is an action $\textstyle W = \int w ( A ) d x$ where $d x = d x ^ { 1 } { \wedge } . . . { \wedge } d x ^ { n }$ , we may vary it. With $A = a ^ { - 1 } d a$ we may introduce $\lambda = a ^ { - 1 } \delta a \in { \mathcal { G } } = \wedge ^ { \cup } T ^ { * } \otimes { \mathcal { G } }$ and get $\delta A _ { i } ^ { \tau } = \partial _ { i } \lambda ^ { \tau } - c _ { \rho \sigma } ^ { \tau } A _ { i } ^ { \rho } \lambda ^ { \sigma }$ (\[20\], p 180-185). Setting $\partial w / \partial A = \mathcal { A } = ( \mathcal { A } _ { \tau } ^ { i } ) \in \wedge ^ { n - 1 } T ^ { * } \otimes \mathcal { G } ^ { * }$ , we obtain the Poincaré equations $\partial _ { i } \mathcal { A } _ { \tau } ^ { i } + c _ { \rho \tau } ^ { \sigma } \mathcal { A } _ { i } ^ { \rho } \mathcal { A } _ { \sigma } ^ { i } = 0$ as the adjoint of the previous operator (up to sign). Setting now $( \delta a ) a ^ { - 1 } = \mu \in { \mathcal G }$ ,we get the adjoint representation $\lambda = a ^ { - 1 } ( ( \delta a ) a ^ { - 1 } ) a = A d ( a ) \mu$ while, introducing $\boldsymbol { B }$ such that $B \mu = A \lambda$ , we get the divergence-like equations $\partial _ { i } B _ { \sigma } ^ { i } = 0$ . In a different setting, if $G$ acts on $X$ , let $\{ \theta _ { \tau } \mid 1 \leq \tau \leq p = d i m ( G ) \}$ be a basis of infinitesimal generators of the action. If $\boldsymbol { \mu } = ( \mu _ { 1 } , . . . , \mu _ { n } )$ is a multi-index of length $| \ \mu \ | = \mu _ { 1 } + \ldots + \mu _ { n }$ and $\mu + 1 _ { i } = ( \mu _ { 1 } , . . , \mu _ { i - 1 } , \mu _ { i } + 1 , \mu _ { i + 1 } , . . . , \mu _ { n } )$ , we may introduce the Lie algebroid $R _ { q } \subset J _ { q } ( T )$ with sections defined by $\xi _ { \mu } ^ { k } ( x ) = \lambda ^ { \tau } ( x ) \partial _ { \mu } \theta _ { \tau } ^ { k } ( x )$ for an arbitrary section $\lambda \in \wedge ^ { 0 } T ^ { * } \otimes { \mathcal G }$ and the trivially involutive operator $j _ { q } : T \to J _ { q } ( T ) : \theta \to ( \partial _ { \mu } \theta , 0 \leq | \mu | \leq q )$ of order $q$ . We finally obtain the Spencer operator $d : R _ { q + 1 } \to T ^ { * } \otimes R _ { q }$ through the chain rule for derivatives \[17\]: $$ ( d \xi _ { q + 1 } ) _ { \mu , i } ^ { k } ( x ) = \partial _ { i } \xi _ { \mu } ^ { k } ( x ) - \xi _ { \mu + 1 _ { i } } ^ { k } ( x ) = \partial _ { i } \lambda ^ { \tau } ( x ) \partial _ { \mu } \theta _ { \tau } ^ { k } ( x ) $$ When $q$ is large enough to have an isomorphism $R _ { q + 1 } \simeq R _ { q } \simeq \Lambda ^ { 0 } T ^ { * } \otimes \mathcal { G }$ and the following linear Spencer sequence in which the operators $D _ { r }$ are induced by $d$ as above: $$ 0 \longrightarrow \Theta \xrightarrow { j _ { q } } R _ { q } \xrightarrow { D _ { 1 } } T ^ { * } \otimes R _ { q } \xrightarrow { D _ { 2 } } \wedge ^ { 2 } T ^ { * } \otimes R _ { q } \xrightarrow { D _ { 3 } } . . . \xrightarrow { D _ { n } } \wedge ^ { n } T ^ { * } \otimes R _ { q } \longrightarrow 0 $$ is isomorphic to the linear gauge sequence but with a completely different meaning because $G$ is now acting on $X$ and $\Theta \subset T$ is such that $[ \Theta , \Theta ] \subset \Theta$ . Surprisingly, these results have NEVER been used in the study of the M, S and K metrics \[13\]. It is not evident, as we already saw, that the projective transformation $y = ( a x + b ) / ( c x + d )$ of the real line with $( a , b , c , d )$ constants is the generic solution of the third order OD equation $( y _ { x x x } / y _ { x } ) - ( 3 / 2 ) ( y _ { x x } / y _ { x } ) ^ { 2 } = 0$ and has three infinitesimal generators $( \partial _ { x } , x \partial _ { x } , \frac{1}{2} x^2 \partial _ { x } )$ providing the generic solution of the linearized OD equation $\xi _ { x x x } = 0$ . The situation of contact transformations met in Section 2 and of Lie pseudogroups in general needs new differential geometric tools as follows, not known by physicists because they necessarily involve the Spencer operator \[20\].  ## LIE ALGEBROIDS If $R _ { q } \subset J _ { q } ( E )$ is a system of order $q$ on $E$ , then $R _ { q + r } = \rho _ { r } ( R _ { q } ) = J _ { r } ( R _ { q } ) \cap J _ { q + r } ( E ) \subset J _ { r } ( J _ { q } ( E ) )$ is called the $r$ -prolongation of $R _ { q }$ and is thus the system obtained after differentiating r times. In actual practice, if the system is defined by PDE $\Phi ^ { \tau } \equiv a _ { k } ^ { \tau \mu } ( x ) \xi _ { \mu } ^ { k } = 0$ the first prolongation is defined by adddng the PDE $d _ { i } \Phi ^ { \tau } \equiv a _ { k } ^ { \tau \mu } ( x ) \xi _ { \mu + 1 _ { i } } ^ { k } + \partial _ { i } a _ { k } ^ { \tau \mu } ( x ) \xi _ { \mu } ^ { k } = 0$ .A Ccordingly, $\xi _ { q } \in R _ { q } \Leftrightarrow a _ { k } ^ { \tau \mu } ( x ) \xi _ { \mu } ^ { k } ( x ) = 0$ and $\xi _ { q + 1 } \in R _ { q + 1 } \Leftrightarrow a _ { k } ^ { \tau \mu } ( x ) \xi _ { \mu + 1 _ { i } } ^ { k } ( x ) + \partial _ { i } a _ { k } ^ { \tau \mu } ( x ) \xi _ { \mu } ^ { k } ( x ) = 0$ as identitties on $X$ or at least over an open subset $U \subset X$ . We finally obtain: $$ \partial _ { i } \Phi ^ { \tau } - d _ { i } \Phi ^ { \tau } \equiv a _ { k } ^ { \tau \mu } ( x ) \big ( \partial _ { i } \xi _ { \mu } ^ { k } ( x ) - \xi _ { \mu + 1 _ { i } } ^ { k } ( x ) \big ) = 0 \Rightarrow d \xi _ { q + 1 } \in T ^ { * } \otimes { \cal R } _ { q } $$ and the Spencer operator restricts to $d : R _ { q + 1 } \to T ^ { * } \otimes R _ { q }$ . Definition 3.B.1: We set $R _ { q + r } ^ { ( s ) } = \pi _ { q + r } ^ { q + r + s } ( R _ { q + r + s } )$ and this system is called the $s$-prolonged system of order $q+r$. For Example, if we consider the second order Killing system $R _ { 1 } \subset J _ { 1 } ( T )$ defined by $x ^ { 2 } \xi _ { 1 } ^ { 1 } + \xi ^ { 2 } = 0 , \xi _ { 2 } ^ { 1 } = 0$ then $R _ { 1 } ^ { ( 1 ) } \subset R _ { 1 }$ is defined by adding $\xi _ { 1 } ^ { 1 } + \xi _ { 2 } ^ { 2 } = 0$ , even though both systems have the same ” solutions”. Definition 3.B.2: The symbol of $R _ { q }$ is the family $g _ { q } = R _ { q } \cap S _ { q } T ^ { * } \otimes E$ of vector spaces over $X$ . The symbol $g _ { q + r }$ of $\mathcal { R } _ { q + r }$ only depends on $g _ { q }$ by a direct prolongation procedure. We may define the vector bundle $F _ { 0 }$ over $\mathcal { R } _ { q }$ by the short exact sequence $0 \to R _ { q } \to J _ { q } ( E ) \to F _ { 0 } \to 0$ and we have the exact induced sequence $0 \rightarrow g_q \rightarrow S_ qT^*\otimes E \rightarrow F_0$ . When $\mid \mu \mid = q$ , we obtain: $$ g _ { q } = \{ v _ { \mu } ^ { k } \in S _ { q } T ^ { * } \otimes E \mid a _ { k } ^ { \tau \mu } ( x ) v _ { \mu } ^ { k } = 0 \} , \mid \mu \mid = q $$ $$ \Rightarrow g _ { q + r } = \rho _ { r } ( g _ { q } ) = \{ v _ { \mu + \nu } ^ { k } \in S _ { q + r } T ^ { * } \otimes E \mid a _ { k } ^ { \tau \mu } ( x ) v _ { \mu + \nu } ^ { k } = 0 \} , $$ $$ \mid \mu \mid = q , \mid \nu \mid = r $$ In general, neither $g _ { q }$ nor $g _ { q + r }$ are vector bundles over $X$ as can be seen in the simple example $x y _ { x } - y = 0 \Rightarrow x y _ { x x } = 0$ . On $\wedge ^ { s } T ^ { * }$ we may introduce the usual bases $\{ d x ^ { I } = d x ^ { i _ { 1 } } \wedge . . . \wedge d x ^ { i _ { s } } \}$ where we have set $I = ( i _ { 1 } < \ldots < i _ { s } )$ . In a purely algebraic setting, one has: Proposition 3.B.3: There exists a map $\delta : \wedge ^ { s } T ^ { * } \otimes S _ { q + 1 } T ^ { * } \otimes E \to \wedge ^ { s + 1 } T ^ { * } \otimes S _ { q } T ^ { * } \otimes E$ which restricts to $\delta : \wedge ^ { s } T ^ { * } \otimes g _ { q + 1 } \to \wedge ^ { s + 1 } T ^ { * } \otimes g _ { q }$ and $\delta ^ { 2 } = \delta \circ \delta = 0$ . Proof: Let us introduce the map: $\omega = \{ \omega _ { \mu } ^ { k } = v _ { \mu , I } ^ { k } d x ^ { I } \}$ $( \delta \omega ) _ { \mu } ^ { k } = d x ^ { i } \wedge \omega _ { \mu + 1 _ { i } } ^ { k }$ $( \delta ^ { 2 } \omega ) _ { \mu } ^ { k } = d x ^ { i } \wedge d x ^ { j } \wedge \omega _ { \mu + 1 _ { i } + 1 _ { j } } ^ { k } = 0$ $a _ { k } ^ { \tau \mu } ( \delta \omega ) _ { \mu } ^ { k } = d x ^ { i } \wedge ( a _ { k } ^ { \tau \mu } \omega _ { \mu + 1 _ { i } } ^ { k } ) = 0$ The kernel of each $\delta$ in the first case is equal to the image of the preceding $\delta$ but this may no longer be true in the restricted case and we set: Definition 3.B.4: Let $B _ { q + r } ^ { s } ( g _ { q } ) \subseteq Z _ { q + r } ^ { s } ( g _ { q } )$ and $H _ { q + r } ^ { s } ( g _ { q } ) = Z _ { q + r } ^ { s } ( g _ { q } ) / B _ { q + r } ^ { s } ( g _ { q } )$ with $H ^ { s } ( g _ { q } ) =$ $H _ { q } ^ { s } ( g _ { q } )$ be the coboundary space $i m ( \delta )$ , cocycle space $k e r ( \delta )$ and cohomology space at $\wedge ^ { s } T ^ { * } \otimes g _ { q + r }$ of the restricted $g _ { q }$ is said to be s-acyclic if $H _ { q + r } ^ { 1 } = \ldots = H _ { q + r } ^ { s } = 0 , \forall r \geq 0$. The symbol $g_q$ is said to be of finite type if $g _ { q + r } = 0$ for $r$ large enough. In particular, if $g _ { q }$ is involutive and finite type, then $g _ { q } = 0$ . Finally, $S _ { q } T ^ { * } \otimes E$ is involutive for any $q \geq 0$ if we set $S _ { 0 } T ^ { \ast } \otimes E = E$ . A first point is provided by the following useful but technical results. As we do not want to provide details about groupoids, we shall introduce a ” copy” $Y$ (target) of $X$ (source) and define simply a Lie pseudogroup $\Gamma \subseteq a u t ( X )$ as a group of transformations solutions of a (in general nonlinear) system $\mathcal { R } _ { q }$ , such that, whenever $y = f ( x ) , z = g ( y ) \in \Gamma$ can be composed, then $z = g \circ f ( x ) \in \Gamma$ , $x = f ^ { - 1 } ( y ) \in \Gamma$ and $y = i d ( x ) = x \in \Gamma$ . Setting $y = x + t \xi ( x ) + . . .$ and passing to the limit when $t \to 0$ , we may linearize the later system and obtain a (linear) system $R _ { q } \subset J _ { q } ( T )$ such that $[ \Theta , \Theta ] \subset \Theta$ . We may use the Frobenius theorem in order to find a generating fundamental set of differential invariants $\left\{ \Phi ^ { \tau } ( y _ { q } ) \right\}$ up to order $q$ which are such that $\Phi ^ { \tau } ( \bar { y } _ { q } ) = \Phi ^ { \tau } ( y _ { q } )$ whenever $\bar { y } = g ( y ) \in \Gamma$ . We obtain the Lie form $\Phi ^ { \tau } ( y _ { q } ) = \Phi _ { \tau } ( i d _ { q } ( x ) ) = \Phi ^ { \tau } ( j _ { q } ( i d ) ( x ) ) = \omega ^ { \tau } ( x )$ of $\mathcal { R } _ { q }$ . Of course, in actual practice one must use sections of $R _ { q }$ instead of solutions and we now prove why the use of the Spencer operator becomes crucial for such a purpose. Indeed, we may define: $$ \{j_{q+1}(\xi), j_{q+1}(\eta)\} = j_q([\xi, \eta]), \forall \xi, \eta \in T \quad \text{(algebraic bracket)} $$ We may obtain by bilinearity a bracket on ${ \cal J } _ { q } ( T )$ extending the bracket on $T$ . $$ [ \xi _ q , \eta _ q ] = \{ \xi _ { q + 1 } , \eta _ { q + 1 } \} + i ( \xi ) d \eta _ { q + 1 } - i ( \eta ) d \xi _ { q + 1 } , \forall \xi _ q , \eta _ q \in J _ q ( T ) \quad \text{(differential bracket)} $$ which does not depend on the respective lifts $\xi _ { q + 1 }$ and $\eta _ { q + 1 }$ of $\xi _ { q }$ and $\eta _ { q }$ in ${ \cal J } _ { q + 1 } ( T )$ . This bracket on sections satisfies the Jacobi identity: $$ [ [ \xi _ q , \eta _ q ] , \zeta _ q ] + [ [ \eta _ q , \zeta _ q ] , \xi _ q ] + [ [ \zeta _ q , \xi _ q ] , \eta _ q ] = 0 , \ \forall \xi _ q , \eta _ q , \zeta _ q \in J _ q ( T ) \quad \text{(Jacobi)} $$ and we set \[20\]: Definition 3.B.5: We say that a vector subbundle $R _ { q } \subset J _ { q } ( T )$ is a system of infinitesimal Lie equations or a Lie algebroid if $[ R _ { q } , R _ { q } ] \subset R _ { q }$ , that is to say $[ \xi _ q , \eta _ q ] \in R _ { q } , \forall \xi _ q , \eta _ q \in R _ { q }$ . Such a definition can be tested by means of computer algebra. We shall also say that $R _ { q }$ is transitive if we have the short exact sequence $0 \to R _ { q } ^ { 0 } \to R _ { q } \stackrel { \pi _ { 0 } ^ { q } } { \to } T \to 0$ . Theorem 3.B.6: The bracket is compatible with prolongations: $$ [ R _ { q } , R _ { q } ] \subset { \cal R } _ { q } \quad \Rightarrow \quad [ R _ { q + r } , R _ { q + r } ] \subset { \cal R } _ { q + r } , \forall r \geq 0 $$ Proof: When $r = 1$ , we have $\rho _ { 1 } ( R _ { q } ) = R _ { q + 1 } = \{ \xi _ { q + 1 } \in J _ { q + 1 } ( T ) \mid \xi _ { q } \in R _ { q } , d \xi _ { q + 1 } \in T ^ { * } \otimes R _ { q } \}$ and we just need to use the following formulas showing how $d$ acts on the various brackets if we set $L ( \xi _ { 1 } ) \zeta = [ \xi , \zeta ] + i ( \zeta ) d \xi _ { 1 }$ (See \[22\] and \[26\] or \[32\] for more details): $$ \begin{aligned} i ( \zeta ) d \{ \xi _ { q + 1 } , \eta _ { q + 1 } \} & = \{ i ( \zeta ) d \xi _ { q + 1 } , \eta _ { q } \} + \{ \xi _ { q } , i ( \zeta ) d \eta _ { q + 1 } \} , \quad \forall \zeta \in T \\ i ( \zeta ) d [ \xi _ { q + 1 } , \eta _ { q + 1 } ] & = [ i ( \zeta ) d \xi _ { q + 1 } , \eta _ { q } ] + [ \xi _ { q } , i ( \zeta ) d \eta _ { q + 1 } ] \\ & \quad + i ( L ( \eta _ { 1 } ) \zeta ) d \xi _ { q + 1 } - i ( L ( \xi _ { 1 } ) \zeta ) d \eta _ { q + 1 } \end{aligned} $$ The right member of the second formula is a section of $R _ { q }$ whenever $\xi _ { q + 1 } , \eta _ { q + 1 } \in R _ { q + 1 }$ . The first formula may be used when $R _ { q }$ is formally integrable. Corollary 3.B.7: The bracket is compatible with the PP procedure: $$ [ R _ { q } , R _ { q } ] \subset { \cal R } _ { q } \quad \Rightarrow \quad [ R _ { q + r } ^ { ( s ) } , R _ { q + r } ^ { ( s ) } ] \subset { \cal R } _ { q + r } ^ { ( s ) } , \forall r , s \geq 0 $$ Example 3.B.8: When $n = 1$ , one has the unusual successive formulas: $$ [ \xi , \eta ] = \xi \partial _ { x } \eta - \eta \partial _ { x } \xi $$ $$ ( [ \xi _ { 1 } , \eta _ { 1 } ] ) _ { x } = \xi \partial _ { x } \eta _ { x } - \eta \partial _ { x } \xi _ { x } $$ $$ ( [ \xi _ { 2 } , \eta _ { 2 } ] ) _ { x x } = \xi _ { x } \eta _ { x x } - \eta _ { x } \xi _ { x x } + \xi \partial _ { x } \eta _ { x x } - \eta \partial _ { x } \xi _ { x x } $$ $$ ( [ \xi _ { 3 } , \eta _ { 3 } ] ) _ { x x x } = 2 \xi _ { x } \eta _ { x x x } - 2 \eta _ { x } \xi _ { x x x } + \xi \partial _ { x } \eta _ { x x x } - \eta \partial _ { x } \xi _ { x x x } $$ They can be used for linear ( $\xi _ { x } = 0$ ), affine ( $\xi _ { x x } = 0$ ) or projective $\xi _ { x x x } = 0$ )transformations. Example 3.B.9: With $n = m = 2$ and $q = 1$ , let us consider the Lie pseudogroup $\Gamma \subset a u t ( X )$ of finite transformations $y = f ( x )$ such that $y ^ { 2 } d y ^ { 1 } = x ^ { 2 } d x ^ { 1 } = \alpha$ . Setting $y = x + t \xi ( x ) + . . .$ and linearizing, we get the Lie operator $\mathcal { D } \xi = \mathcal { L } ( \xi ) \alpha$ where $\mathcal { L }$ is the Lie derivative and the system $R _ { 1 } \subset J _ { 1 } ( T )$ of linear infinitesimal Lie equations already met in Definition 3.B.1. Let us choose the two sections $\xi _ { 1 } , \eta _ { 1 } \in J _ { 1 } ( T )$ : $$ \begin{array} { r l } & { \xi _ { 1 } = \big \{ \xi ^ { 1 } = 0 , \xi ^ { 2 } = - x ^ { 2 } , \xi _ { 1 } ^ { 1 } = 1 , \xi _ { 2 } ^ { 1 } = 0 , \xi _ { 1 } ^ { 2 } = 0 , \xi _ { 2 } ^ { 2 } = 0 \big \} \in R _ { 1 } } \\ & { \eta _ { 1 } = \big \{ \eta ^ { 1 } = x ^ { 2 } , \eta ^ { 2 } = 0 , \eta _ { 1 } ^ { 1 } = 0 , \eta _ { 2 } ^ { 1 } = - x 2 , \eta _ { 1 } ^ { 2 } = 0 , \eta _ { 2 } ^ { 2 } = 1 \big \} \in R _ { 1 } } \end{array} $$ The reader may check that $[ \xi _ { 1 } , \eta _ { 1 } ] \in R _ { 1 }$ . However, we have the strict inclusion $R _ { 1 } ^ { ( 1 ) } \subset R _ { 1 }$ defined by the additional equation $\xi _ { 1 } ^ { 1 } + \xi _ { 2 } ^ { 2 } = 0$ and thus $\xi _ { 1 } , \eta _ { 1 } \notin R _ { 1 } ^ { ( 1 ) }$ though we have $[ R _ { 1 } ^ { ( 1 ) } , R _ { 1 } ^ { ( 1 ) } ] \subset R _ { 1 } ^ { ( 1 ) }$ , a result not evident and lacking in \[14\] because $\xi _ { 1 }$ and $\eta _ { 1 }$ have nothing to do with solutions.  ## JANET VERSUS SPENCER When $R _ { q } \subset J _ { q } ( E )$ is an involutive system on $E$ , there is a Fundamental Diagram $I$ with exact columns, where each operator is induced by the extension $d : \wedge ^ { r } T ^ { * } \otimes J _ { q + 1 } ( E ) \rightarrow \wedge ^ { r + 1 } T ^ { * } \otimes J _ { q } ( E ) :$ $\alpha \otimes \xi _ { q + 1 } \to d \alpha \otimes \xi _ { q } + ( - 1 ) ^ { r } \alpha \wedge d \xi _ { q + 1 }$ of the Spencer operator $d : J _ { q + 1 } ( E ) \to T ^ { * } \otimes J _ { q } ( E ) : \xi _ { q + 1 } \to$ $j _ { 1 } ( \xi _ { q } ) - \xi _ { q + 1 }$ . The upper sequence is the Spencer sequence while the lower sequence is the Janet sequence and the sum $d i m ( C _ { r } ) + d i m ( F _ { r } ) = d i m ( C _ { r } ( E ) )$ does not depend on the system while the epimorphisms $\Phi _ { r }$ are induced inductively by $\Phi = \Phi _ { 0 }$ [20, 29]: $$ \begin{array}{cccccccccccccccc} & & & & & & 0 & & 0 & & 0 & & \\ & & &&&& \downarrow && \downarrow && \downarrow & & & & \\ & & 0& \rightarrow & \Theta &\stackrel{j_q}{\longrightarrow} &C_0 &\stackrel{d_1}{\longrightarrow}& C_1 &\stackrel{d_2}{\longrightarrow} ... \stackrel{d_n}{\rightarrow}& C_n &\rightarrow & 0 \\ &&&&&& \downarrow & & \downarrow & &\downarrow & & \\ & & 0 & \rightarrow & E & \stackrel{j_q}{\longrightarrow} & C_0(E) & \stackrel{d_1}{\longrightarrow} & C_1(E) &\stackrel{d_2}{\longrightarrow} ... \stackrel{d_n}{\longrightarrow} & C_n(E) & \rightarrow & 0 \\ & & & & \parallel && \hspace{5mm}\downarrow {\Phi}_0 & &\hspace{5mm} \downarrow {\Phi}_1 & & \hspace{5mm}\downarrow {\Phi}_n & & \\ 0 & \rightarrow & \Theta &\rightarrow & E & \stackrel{D}{\longrightarrow} & F_0 & \stackrel{D_1}{\longrightarrow} & F_1 & \stackrel{D_2}{\longrightarrow} ... \stackrel{D_n}{\longrightarrow} & F_n & \rightarrow & 0 \\ &&&&& & \downarrow & & \downarrow & &\downarrow & & \\ &&&&& & 0 && 0 &&0 & & \end{array} $$ For later computations, the sequence $J _ { q + 1 } ( E ) \stackrel { d } { \longrightarrow } T ^ { * } \otimes J _ { q } ( E ) \stackrel { d } { \longrightarrow } \wedge { } ^ { 2 } T ^ { * } \otimes J _ { q - 1 } ( E )$ can be described by the images $\partial _ { i } \xi _ { \mu } ^ { k } - \xi _ { \mu + 1 _ { i } } ^ { k } = X _ { \mu , i } ^ { k }$ leading to the identities: $$ \partial _ { i } X _ { \mu , j } ^ { k } - \partial _ { j } X _ { \mu , i } ^ { k } + X _ { \mu + 1 _ { j } , i } ^ { k } - X _ { \mu + 1 _ { i } , j } ^ { k } = 0 $$ We also recall that the linear Spencer sequence for a Lie group of transformations $G \times X \to X$ , which essentially depends on the action because infinitesimal generators are needed, is locally isomorphic to the linear gauge sequence which does not depend on the action any longer as it is the tensor product of the Poincaré sequence by the Lie algebra $\mathcal { G }$ of $G$ . The main idea will be to introduce and compare the three Lie groups of transformations: The Poincaré group of transformations with 10 parameters leading to the Killing system $R _ { 2 }$ : $$ \Omega _ { i j } \equiv ( L ( \xi _ { 1 } ) \omega ) _ { i j } \equiv \omega _ { r j } ( x ) \xi _ { i } ^ { r } + \omega _ { i r } ( x ) \xi _ { j } ^ { r } + \xi ^ { r } \partial _ { r } \omega _ { i j } ( x ) = 0 $$ $$ \Gamma _ { i j } ^ { k } \equiv ( L ( \xi _ { 2 } ) \gamma ) _ { i j } ^ { k } \equiv \xi _ { i j } ^ { k } + \gamma _ { r j } ^ { k } ( x ) \xi _ { i } ^ { r } + \gamma _ { i r } ^ { k } ( x ) \xi _ { j } ^ { r } - \gamma _ { i j } ^ { r } ( x ) \xi _ { r } ^ { k } + \xi ^ { r } \partial _ { r } \gamma _ { i j } ^ { k } ( x ) = 0 $$ The Weyl group of transformations with 11 parameters leading to the Weyl system $\tilde { R } _ { 2 }$ . $$ \Omega _ { i j } \equiv \omega _ { r j } ( x ) \xi _ { i } ^ { r } + \omega _ { i r } ( x ) \xi _ { j } ^ { r } + \xi ^ { r } \partial _ { r } \omega _ { i j } ( x ) = 2 \Lambda ( x ) \omega _ { i j } ( x ) $$ $$ \begin{aligned} \Gamma _ { i j } ^ { k } \equiv & \xi _ { i j } ^ { k } + \gamma _ { r j } ^ { k } ( x ) \xi _ { i } ^ { r } + \gamma _ { i r } ^ { k } ( x ) \xi _ { j } ^ { r } - \gamma _ { i j } ^ { r } ( x ) \xi _ { r } ^ { k } + \xi ^ { r } \partial _ { r } \gamma _ { i j } ^ { k } ( x ) \\ = & 0 \end{aligned} $$ The conformal group of transformations with 15 parameters (4 translations $^ \mathrm { ~ + ~ 6 ~ }$ rotations $+ \ 1$ dilatation $^ { + 4 }$ elations) leading to the conformal Killing system $\hat { R } _ { 2 }$ : $$ \Omega _ { i j } \equiv \omega _ { r j } ( x ) \xi _ { i } ^ { r } + \omega _ { i r } ( x ) \xi _ { j } ^ { r } + \xi ^ { r } \partial _ { r } \omega _ { i j } ( x ) = 2 \Lambda ( x ) \omega _ { i j } ( x ) $$ $$ \begin{aligned} \Gamma _ { i j } ^ { k } \equiv & \xi _ { i j } ^ { k } + \gamma _ { r j } ^ { k } ( x ) \xi _ { i } ^ { r } + \gamma _ { i r } ^ { k } ( x ) \xi _ { j } ^ { r } - \gamma _ { i j } ^ { r } ( x ) \xi _ { r } ^ { k } + \xi ^ { r } \partial _ { r } \gamma _ { i j } ^ { k } ( x ) \\ = & \delta _ { i } ^ { k } \Lambda _ { j } ( x ) + \delta _ { j } ^ { k } \Lambda _ { i } ( x ) - \omega _ { i j } ( x ) \omega ^ { k r } ( x ) \Lambda _ { r } ( x ) \end{aligned} $$ where one has to eliminate the arbitrary function $\Lambda ( x )$ and 1-form $\Lambda _ { i } ( x ) d x ^ { i }$ for finding sections, replacing the ordinary Lie derivative $\mathcal { L } ( \xi )$ by the formal Lie derivative $L ( \xi _ { q } )$ , that is replacing $j _ { q } ( \xi )$ by $\xi _ { q }$ when needed. When $n = 4$ , $\hat { R } _ { 2 }$ is FI but $\hat { g } _ { 2 }$ is only 2-acyclic while $\hat { g } _ { 3 } = 0$ and we have for the involutive $\hat { R } _ { 3 } \simeq \hat { R } _ { 2 }$ (See \[26\] for details and counterexamples): $$ \begin{array}{rccccccccccccccl} &&&&& 0 &&0&&0& &0& & 0 & & \\ &&&&& \downarrow && \downarrow && \downarrow & & \downarrow & & \downarrow & \\ & 0& \rightarrow& \hat{\Theta} &\stackrel{j_3}{\rightarrow}&15 &\stackrel{d_1}{\rightarrow}& 60 &\stackrel{d_2}{\rightarrow} & 90 &\stackrel{d_3}{\rightarrow} & 60 & \stackrel{d_4}{\rightarrow}& 15 &\rightarrow 0 \\ &&&&& \downarrow & & \downarrow & & \downarrow & &\downarrow & & \downarrow & \\ & 0 & \rightarrow & 4 & \underset 3{\stackrel{j_3}{\rightarrow}} & 140 & \stackrel{d_1}{\rightarrow} & 420 &\stackrel{d_2}{\rightarrow} & 504 &\stackrel{d_3}{\rightarrow} & 280 & \stackrel{d_4}{\rightarrow} & 60 & \rightarrow 0 \\ & & & \parallel && \hspace{5mm}\downarrow {\Phi}_0 & &\hspace{5mm} \downarrow {\Phi}_1 & & \hspace{5mm}\downarrow {\Phi}_2 & & \hspace{5mm}\downarrow {\Phi}_3 & & \hspace{5mm} \downarrow {\Phi}_4 & \\ 0 \rightarrow & \hat{\Theta} &\rightarrow & 4 & \underset 3{\stackrel{\hat{D}}{\rightarrow}} & 125 & \underset 1{\stackrel{\hat{D}_1}{\rightarrow}} & 360 & \underset 1{\stackrel{{\hat{D}}_2}{\rightarrow}} & 414 & \underset 1{\stackrel{{\hat{D}}_3}{\rightarrow}} & 220 & \underset 1{\stackrel{{\hat{D}}_4}{\rightarrow} }& 45 & \rightarrow 0 \\ &&&&& \downarrow & & \downarrow & & \downarrow & &\downarrow & & \downarrow & \\ &&&&& 0 && 0 && 0 &&0 & & 0 & \end{array} $$ The top Spencer sequence is the tensor product of the Poincaré sequence by the Lie algebra $\hat { \mathcal G }$ of dimension 15 and we may use the inclusions $R _ { 2 } \subset \bar { R } _ { 2 } \subset \bar { R } _ { 2 } \subset J _ { 2 } ( T )$ with $1 0 < 1 1 < 1 5 < 6 0$ . Working by induction, the minimum formally exact resolution on the jet level is: $$ 0 \longrightarrow \hat { \Theta } \longrightarrow 4 \longrightarrow 9 \longrightarrow 1 0 \longrightarrow 9 \longrightarrow 4 \longrightarrow 0 $$ with "up and down" orders that must be compared to the above canonical Janet sequence. Of course, finding such numbers can be done by means of computer algebra (arXiv: 1603.05030) or through combinatorics (exercise !) but it will never prove that such a sequence is formally exact as it will involve enormous matrices (up to $8 4 0 \textrm { x } 1 1 3 4$ !!!) and cannot be achieved without the help of the Spencer $\delta$ -cohomology, still never introduced in GR or conformal geometry \[20, 26\]. When $\omega$ is the M metric, it follows that $\gamma = 0$ and we obtain therefore: $$ X _ { j , i } ^ { r } - X _ { r i , j } ^ { r } = ( \partial _ { i } \xi _ { r j } ^ { r } - \xi _ { r i j } ^ { r } ) - ( \partial _ { j } \xi _ { r i } ^ { r } - \xi _ { r i j } ^ { r } ) = \partial _ { i } \xi _ { r j } ^ { r } - \partial _ { j } \xi _ { r i } ^ { r } = n \left( \partial _ { i } \Lambda _ { j } - \partial _ { j } \Lambda _ { i } \right) ) $$ Dividing by $n$ , we may thus obtain $( F _ { i j } = \partial _ { i } \Lambda _ { j } - \partial _ { j } \Lambda _ { i } ) \in \wedge ^ { 2 } T ^ { * }$ from $X _ { r j , i } ^ { r } \in T ^ { * } \otimes \hat { g } _ { 2 } \subset C _ { 1 }$ with $d { \boldsymbol { F } } = 0$ because $\hat { g } _ { 3 } = 0$ and thus $\xi _ { r i j } ^ { k } = 0$ in $S _ { 3 } T ^ { * } \otimes T$ . This result is solving the dream of $\mathrm { H }$ . Weyl for exhibiting the conformal origin of electromagnetism in \[22\]. It is however completely contradicting the standard approach of classical gauge theory based on the group $U ( 1 )$ which is not acting on space-time. In addition, the EM field $\mathrm { F }$ is a section of the first Spencer bundle $C _ { 1 }$ in the image of $D _ { 1 }$ because $( \Lambda , \Lambda _ { i } ) \in C _ { 0 } = \hat { R } _ { 3 } \simeq \hat { R } _ { 2 }$ . We apply the linearization procedure to the Riemann tensor: $$ \rho = ( \rho _ { l , i j } ^ { k } = \partial _ { i } \gamma _ { l j } ^ { k } - \partial _ { j } \gamma _ { l i } ^ { k } + \gamma _ { l j } ^ { r } \gamma _ { r i } ^ { k } - \gamma _ { l i } ^ { r } \gamma _ { r j } ^ { k } ) \in \mathbb { \Lambda } ^ { 2 } T ^ { * } \otimes T ^ { * } \otimes T ^ { } $$ Now, as the linearization $\Gamma \in S _ { 2 } T ^ { * } \otimes T$ of $\gamma$ is a tensor, the linearization $R$ of $\rho$ becomes: $$ \begin{array} { r } { R _ { l , i j } ^ { k } = d _ { i } \Gamma _ { l j } ^ { k } - d _ { j } \Gamma _ { l i } ^ { k } + \gamma _ { l j } ^ { r } \Gamma _ { r i } ^ { k } - \gamma _ { l i } ^ { r } \Gamma _ { r j } ^ { k } + \gamma _ { r i } ^ { k } \Gamma _ { l j } ^ { r } - \gamma _ { r j } ^ { k } \Gamma _ { l i } ^ { r } } \end{array} $$ If $\nabla$ is the covariant derivative, we have $\nabla _ { r } \omega _ { i j } = \partial _ { r } \omega _ { i j } - \gamma _ { i r } ^ { s } \omega _ { s j } - \gamma _ { j r } ^ { s } \omega _ { i s } = 0$ and we may move down the index $k$ . Then, using $r$ as a dumb index, we may consider the first order equations: $$ ( L ( \xi _ { 1 } ) \rho ) _ { k l , i j } \equiv R _ { k l , i j } \equiv \rho _ { r l , i j } \xi _ { k } ^ { r } + \rho _ { k r , i j } \xi _ { l } ^ { r } + \rho _ { k l , r j } \xi _ { i } ^ { r } + \rho _ { k l , i r } \xi _ { j } ^ { r } + \xi ^ { r } \partial _ { r } \rho _ { k l , i j } = 0 $$ that can be considered as an infinitesimal first order variation of $\rho$ . As for the Ricci tensor $( \rho _ { i j } ) \in S _ { 2 } T ^ { * }$ , we notice that $\rho _ { i j } = \rho _ { i , r j } ^ { r } = 0 \Rightarrow R _ { i j } \equiv \rho _ { r j } \xi _ { i } ^ { r } + \rho _ { i r } \xi _ { j } ^ { r } + \xi ^ { r } \partial _ { r } \rho _ { i j } = 0$ . Using now a cyclic sums on $( i j r )$ , the Bianchi identities are: $$ \begin{array} { r } { \beta _ { k l , i j r } \equiv \nabla _ { r } \rho _ { k l , i j } + \nabla _ { i } \rho _ { k l , j r } + \nabla _ { j } \rho _ { k l , r i } = 0 \qquad \Leftrightarrow \qquad \beta \equiv \underset { c y c l } { \Sigma } ( \partial \rho - \gamma \rho ) = 0 } \end{array} $$ Their linearizations $B _ { k l , i j r } = 0$ are sections of the vector bundle $F _ { 2 }$ in the short exact sequence: $$ 0 \to F _ { 2 } \to \wedge ^ { 3 } T ^ { * } \otimes g _ { 1 } \overset { \delta } { \longrightarrow } \wedge ^ { 4 } T ^ { * } \otimes T \to 0 $$ $$ \begin{aligned} \dim(F_2) &= \left( \frac{n(n-1)(n-2)}{6} \right)\left( \frac{n(n-1)}{2} \right) - \left( \frac{n(n-1)(n-2)(n-3)}{24} \right)n \\ &= \frac{n^2(n^2-1)(n-2)}{24} \end{aligned} $$ because $d i m ( g _ { 1 } ) = n ( n - 1 ) / 2$ for any nondegenerate metric, that is $2 4 - 4 = 2 0$ when $n = 4$ . Such results cannot be even imagined by somebody not aware of the $\delta$ -cohomology \[30-32\]. We have the linearized cyclic sums of covariant derivatives both with their respective symbolic descriptions, not to be confused with the non-linear corresponding ones: $$ \begin{aligned} B_{kl,rij} &\equiv \nabla_r R_{kl,ij} + \nabla_i R_{kl,jr} + \nabla_j R_{kl,ri} \\ &= 0 \mod(\Gamma) \\ &\Leftrightarrow \frac{\Sigma}{\text{cycl}}(dR - \gamma R - \rho \Gamma) = 0 \end{aligned} $$ We have thus $\omega \longrightarrow \gamma \longrightarrow \rho \longrightarrow \beta$ and the respective linearizations $\Omega \longrightarrow \Gamma \longrightarrow R \longrightarrow B$ . ## Einstein Versus Weyl It remains to prove that, in this new framework, the Ricci tensor only depends on the symbol $\hat { g } _ { 2 } \simeq T ^ { * } \subset S _ { 2 } T ^ { * } \otimes T$ deroonation $\hat { R } _ { 2 } \subset J _ { 2 } ( T )$ $\hat { R } _ { 1 } \subset$ $J _ { 1 } ( T )$ ${ \hat { g } } _ { 1 } \subset T ^ { * } \otimes T$ $\begin{array} { r } { \omega _ { r j } \xi _ { i } ^ { r } + \omega _ { i r } \xi _ { j } ^ { r } - \frac { 2 } { n } \omega _ { i j } \xi _ { r } ^ { r } = 0 } \end{array}$ ing on any conformal factor. The next commutative diagram covers both situations, taking into account that the equations of both the classical and conformal Killing operator are homogeneous. The purely algebraic Spencer map $\delta : g _ { q + 1 } T ^ { * } \otimes g _ { q }$ with symbol $g _ { q } = R _ { q } \cap S _ { q } T ^ { * } \otimes E \subset J _ { q } ( E )$ is induced by $- d$ and all the sequences are exact by definition but perhaps the left column: $$ \begin{array}{rccccccc} & 0 & & 0 & & 0 & \\ & \downarrow & & \downarrow & & \downarrow & \\ 0 \to & g_3 & \to & S_3 T^* \otimes T & \to & S_2 T^* \otimes F_0 & \to F_1 \to 0 \\ & \downarrow \delta & & \downarrow \delta & & \downarrow \delta & \\ 0 \to & T^* \otimes g_2 & \to & T^* \otimes S_2 T^* \otimes T & \to & T^* \otimes T^* \otimes F_0 & \to 0 \\ & \downarrow \delta & & \downarrow \delta & & \downarrow \delta & \\ 0 \to & \wedge^2 T^* \otimes g_1 & \to & \wedge^2 T^* \otimes T^* \otimes T & \to & \wedge^2 T^* \otimes F_0 & \to 0 \\ & \downarrow \delta & & \downarrow \delta & & \downarrow & \\ 0 \to & \wedge^3 T^* \otimes T & = & \wedge^3 T^* \otimes T & \to & 0 & \\ & \downarrow & & \downarrow & & & \\ & 0 & & 0 & & & \end{array} $$ Chasing in this diagram, we discover that $F _ { 1 }$ is just the cohomology $H _ { 1 } ^ { 2 } ( g _ { 1 } )$ of the first vertical column at $\wedge ^ { 2 } T ^ { * } \otimes g _ { 1 }$ that is the quotient of the cocycle bundle $Z _ { 1 } ^ { 2 } ( g _ { 1 } )$ kernel of the map $\delta : \wedge ^ { \cdot 2 } T ^ { * } \otimes g _ { 1 } \wedge ^ { 3 } T ^ { * } \otimes T$ by the coboundary bundle $B _ { 1 } ^ { 2 } ( g _ { 1 } )$ image of the map $\delta : T ^ { * } \otimes g _ { 2 } { \wedge } ^ { 2 } \otimes g _ { 1 }$ . Needless to say that such mathematical methods have never been used in the stuy of classical or conformal Riemannian geometry and the reason for which the reader will never find other references in the mathematical or physical literature (!). THEOREM 3.D.1: Introducing the $\delta$ -cohomology bundles $R i e m a n n = H _ { 1 } ^ { 2 } ( g _ { 1 } )$ at $\land ^ { 2 } T ^ { * } \otimes g _ { 1 }$ and $W e y l = H _ { 1 } ^ { 2 } ( \hat { g } _ { 1 } )$ at $\wedge ^ { 2 } T ^ { * } \otimes \hat { g } _ { 1 }$ while taking into account that $g _ { 1 } \subset { \hat { g } } _ { 1 }$ , $g _ { 2 } = 0$ $\hat { g } _ { 2 } \simeq T ^ { * }$ and $\hat { g } ^ { 3 } = 0$ , we have the commutative and exact " fundamental diagram $I I$ " found in 1983 \[30\]: $$ \begin{array}{ccccccccccc} & & & & & & & & &.& \\ & & & & & & & & 0 & & \\ & & & & & & & & \downarrow & & \\ & & & & & & 0 & & Ricci & & \\ & & & & & & \downarrow & & \downarrow & & \\ & & & & 0 &\longrightarrow & Z^2_1(g_1) & \longrightarrow &Riemann & \longrightarrow & 0 \\ & & & & \downarrow & & \downarrow & & \downarrow & & \\ & & 0 &\longrightarrow & T^*\otimes {\hat{g}}_2 & \stackrel{\delta}{\longrightarrow} & Z^2_1({\hat{g}}_1) & \longrightarrow & Weyl & \longrightarrow & 0 \\ & & & & \downarrow & & \downarrow & & \downarrow & & \\ 0 & \longrightarrow & S_2T^* & \stackrel{\delta}{\longrightarrow}& T^*\otimes T^* &\stackrel{\delta}{\longrightarrow} & {\wedge}^2T^* & \longrightarrow & 0 & & \\ & & & & \downarrow & & \downarrow & & & & \\ & & & & & & & & & \end{array} $$ The splitting sequence $0 \to Ricci \to Riemann \to Weyl \to 0$ provides an unusual interpretation of the successive Ricci, Riemann and Weyl tensor bundles. Similarly, the well known splitting sequence $0 \to S_2 T^* \xrightarrow{\delta} T^* \otimes T^* \xrightarrow{\delta} \wedge^2 T^* \to 0$ provides an unusual conformal interpretation of the EM field $F = (F_{ij} = -F_{ji}) \in \wedge^2 T^*$ in a coherent way with the tentative of H. Weyl in 1918 \[22\]. $$ \begin{aligned} \dim(\mathrm{Riemann}) - \dim(\mathrm{Weyl}) &= \frac{n^2(n^2-1)}{12} - \frac{n(n+1)(n+2)(n-3)}{12} \\ &= \frac{n(n+1)}{2} \\ &= \dim(\mathrm{Ricci}) \end{aligned} $$ and the Weyl operator is of order 3 when $n = 3$ but of order 2 when $n \geq 4$ , a result still neither known but nor even acknowledged today (See arXiv: 1603.05030 for a computer algebra checking by my former PhD student A. Quadrat and \[31\]). We finally point out that the Bianchi-type operator is of order 2 when $n = 4$ but of order 1 when $n = 5$ (See \[26\] for more details). ## Applications Before comparing the Minkowski, Schwarzschild and Kerr metrics as in \[13\] while comparing to \[14\], we point out the purely mathematical fact explaining why gravitational waves and black holes cannot exist, not because of an experimental problem of detection but because of a structural problem of equation. Indeed, introducing the Ricci tensor $\rho _ { i j } = \rho _ { i , r j } ^ { r } = \rho _ { j i }$ as a linear combination of the second order derivatives of the metric $\omega$ , we may linearize it over the locally constant Euclidean or Minkowskian metric in order to obtain a second order linear homogeneous Ricci operator $S_2 T^* \rightarrow S_2 T^* : ( \Omega _ { i j } ) \rightarrow ( R _ { i j } )$ with 4 terms along the formula: $$ 2 R _ { i j } = \omega ^ { r s } ( d _ { r s } \Omega _ { i j } + d _ { i j } \Omega _ { r s } - d _ { r i } \Omega _ { s j } - d _ { s j } \Omega _ { r i } ) = 2 R _ { j i } $$ which is also valid when $n = 2$ with $2 R _ { 1 1 } = 2 R _ { 2 2 } = d _ { 1 1 } \Omega _ { 2 2 } + d _ { 2 2 } \Omega _ { 1 1 } - 2 d _ { 1 2 } \Omega _ { 1 2 }$ and $R _ { 1 2 } = 0$ . Such a result is showing that the important object is not the Riemann tensor or the Einstein tensor but the Ricci tensor according to the previous fundamental diagram II. We obtain therefore: THEOREM 4.1: The so-called defining equations of the gravitational waves that can be found in any textbook of general relativity like \[19\], are nothing else that the adjoint of the Ricci operator which is parametrizing the $C a u c h y = a d ( K i l l i n g )$ operator but is going " backwards", that is from right to left as we saw in the Introduction. Proof. Introducing the test functions $\lambda ^ { i j } = \lambda ^ { j i }$ , we get: $$ \lambda ^ { i j } R _ { i j } = \omega ^ { r s } \lambda ^ { i j } ( d _ { r s } \Omega _ { i j } + d _ { i j } \Omega _ { r s } - d _ { r i } \Omega _ { s j } - d _ { s j } \Omega _ { r i } ) $$ Integrating by parts while setting as usual $\Pi = \omega ^ { i j } d _ { i j }$ for the Dalembert operator and exchanging the dumb indices, we get for the adjoint operator: $$ {\omega}^{ij} d_{ij} {\lambda}^{rs}+ {\omega}^{rs}d_{ij}{\lambda}^{ij}-{\omega}^{sj}d_{ij}{\lambda}^{ri}- {\omega}^{ri}d_{ij}{\lambda}^{sj} = {\sigma}^{rs} $$ that is exactly the equations of the gravitational waves leading to the Cauchy identities: $$ d _ { r } \sigma ^ { r s } \equiv \omega ^ { i j } d _ { r i j } \lambda ^ { r s } + \omega ^ { r s } d _ { r i j } \lambda ^ { i j } - \omega ^ { s j } d _ { r i j } \lambda ^ { r i } - \omega ^ { r i } d _ { r i j } \lambda ^ { s j } = 0 $$ There is absolutely no need to set $d _ { i } \lambda ^ { i j } = 0$ and we obtain the adjoint sequences when $n = 4$ . $$ \begin{array}{ccccccc} & & 4 & \xrightarrow{Killing} & 10 & \xrightarrow{Ricci} & 10 \\ \\ 0 & \leftarrow & 4 & \xleftarrow{Cauchy} & 10 & \xleftarrow{ad(Ricci)} & 10 \end{array} $$ without any reference to the Bianchi operator and the induced div operator or even to the Einstein operator. Hence, gravitational waves cannot exist, not for a problem of detection but for a problem of equation, as we have only obtained " $a$ " parametrization of the Cauchy operator which is not even minimum but we can find a minimum one by keeping only $\lambda _ { i j }$ with $\textit { i } < \textit { j }$ as in \[18\]. Indeed, we have already obtained another parametrization with 20 test functions by using the $B e l t r a m i = a d ( R i e m a n n )$ operator in the Introduction with $n = 4$ , a result proving that the Cauchy operator can be parametrized. For the advanced reader, this intrinsic structural fact means that the differential module defined by the Cauchy operator is torsion-free \[4\]. We finally point out that the Airy, Beltrami or Einstein parametrizations are not responsible for earthquakes as there is no relation between the test functions $\lambda$ and the deformation $\Omega$ of the metric $\omega$ , even if $n = 2$ as we saw with $\lambda ^ { i j } R _ { i j } = \phi ( d _ { 1 1 } \Omega _ { 2 2 } + d _ { 2 2 } \Omega _ { 1 1 } - 2 d _ { 1 2 } \Omega _ { 1 2 } )$ and Airy function $\phi$ . We also point out that the differential module defined by the Riemann or Ricci operator has a nonzero torsion submodule generated by the 10 components of the Weyl tensor which are separately killed by the Dalembert operator. In the case of the Killing operator for the Minkowski metric, according to H. Poincaré, the geometrical and adjoint physical long exact differential sequences of operators acting on tensors, giving order of operators and number of components, are exactly the ones we have presented in the Introduction for various dimensions $n$ . The main problem is that the corresponding Killing operators and systems for the M, S or (worst !) K metrics are not involutive and not even formally integrable, a key problem because the Janet sequences cannot be easily constructed. Hence, our main target will be to apply the PP procedure in each case, a difficult task indeed. As we know that the corresponding Spencer sequence is isomorphic to the tensor product of the Poincaré differential sequence for the exterior derivative by a Lie algebra of dimension $n ( n + 1 ) / 2$ , its adjoint sequence is also exact, the reason for which the physical sequence is also exact, that is each operator is generating the CC of the previous one. In particular, the $C a u c h y = a d ( K i l l i n g )$ operator generates the CC of the $B e l t r a m i = a d ( R i e m a n n )$ operator which generates the CC of the $L a n c z o s = a d ( B i a n c h i )$ operator. We have proved in many recent publications that this result cannot be extended easily to the conformal situation because the acyclicity properties of the symbols highly depend on the dimension $n$ and order $q$ \[31\]. It thus remains to study the cases of the Schwarzschild and Kerr metrics, a much more delicate problem. EXAMPLE 4.2: (Kerr metric) We now write the Kerr metric in Boyer-Lindquist coordinates $( t , r , \theta , \phi )$ . $$ \omega = \frac{\rho^2 - mr}{\rho^2}\,dt^2 - \frac{\rho^2}{\Delta}\,dr^2 - \rho^2\,d\theta^2 - \frac{2amr\sin^2\theta}{\rho^2}\,dt\,d\phi - \left(r^2 + a^2 + \frac{mra^2\sin^2\theta}{\rho^2}\right)\sin^2\theta\,d\phi^2 $$ where we have set $\Delta = r ^ { 2 } - m r + a ^ { 2}$ and $\rho ^ { 2 } = r ^ { 2 } + a ^ { 2 } c o s ^ { 2 } ( \theta )$ as usual and we recover the S-metric when $a = 0$ with $\begin{array} { r } { A ( r ) = 1 - \frac { m } { r } } \end{array}$: $$ \omega = A ( r ) d t ^ { 2 } - ( 1 / A ( r ) ) d r ^ { 2 } - r ^ { 2 } d \theta ^ { 2 } - r ^ { 2 } s i n ^ { 2 } ( \theta ) d \phi $$ We notice that $t$ or $\phi$ do not appear in the coefficients of the metric. We shall change the coordinate system in order to confirm these results by using computer algebra. The idea is to use the so-called " rational polynomial " coefficients as follows with $c = c o s ( \theta ) \Rightarrow d c = - s i n ( \theta ) d \theta \Rightarrow$ $d c ^ { 2 } = ( 1 - c ^ { 2 } ) d \theta ^ { 2 }$ and set $x ^ { 0 } = t$ , $x ^ { 1 } = r$ , $x ^ { 2 } = c = c o s ( \theta )$ , $x ^ { 3 } = \phi$ ). We obtain over the differential field $K = \mathbb { Q } ( a , m ) ( t , r , c , \phi )$ : $$ \omega = \frac{\rho^2 - mr}{\rho^2}\,dt^2 - \frac{\rho^2}{\Delta}\,dr^2 - \frac{\rho^2}{1-c^2}\,dc^2 - \frac{2amr(1-c^2)}{\rho^2}\,dt\,d\phi - (1-c^2)\left(r^2 + a^2 + \frac{ma^2r(1-c^2)}{\rho^2}\right)\,d\phi^2 $$ with now $\Delta = r ^ { 2 } - m r + a ^ { 2 }$ and $\rho ^ { 2 } = r ^ { 2 } + a ^ { 2 } c ^ { 2 }$ and we have $d e t ( \omega ) = - ( r ^ { 2 } + a ^ { 2 } c ^ { 2 } ) ^ { 2 }$ in a coherent way with the fact that the S metric that can be written in the new system of coordinates: $$ \omega = A ( r ) d t ^ { 2 } - \frac { 1 } { A ( r ) } d r ^ { 2 } + \frac { r ^ { 2 } } { ( 1 - c ^ { 2 } ) } d c ^ { 2 } - r ^ { 2 } ( 1 - c ^ { 2 } ) d \phi ^ { 2 } $$ Looking at the symbol $g _ { 1 } \subset T ^ { * } \otimes T$ , elementary linear combinatorics allow to prove \[18, 22\]: $$ \omega _ { 3 3 } \xi _ { 3 } ^ { 3 } + \omega _ { 0 3 } \xi _ { 3 } ^ { 0 } = 0 \ m o d ( \xi ) , \omega _ { 3 3 } \xi _ { 0 } ^ { 3 } + \omega _ { 0 0 } \xi _ { 3 } ^ { 0 } = 0 \ m o d ( \xi ) , \omega _ { 3 3 } \xi _ { 0 } ^ { 0 } - \omega _ { 0 3 } \xi _ { 3 } ^ { 0 } = 0 \ m o d ( \xi ) , $$ Then, multiplying $\Omega _ { 2 2 }$ by $\omega _ { 1 1 }$ , $\Omega _ { 1 1 }$ by $\omega _ { 2 2 }$ and adding, we finally obtain: $$ 2 ( \omega _ { 1 1 } \omega _ { 2 2 } ) ( \xi _ { 1 } ^ { 1 } + \xi _ { 2 } ^ { 2 } ) + \xi \partial ( \omega _ { 1 1 } \omega _ { 2 2 } ) = 0 $$ However, we have also successively: $$ \left\{ \begin{array} { c c l } { { R _ { 0 3 , 0 3 } } } & { { \equiv } } & { { 2 \rho _ { 0 3 , 0 3 } ( \xi _ { 0 } ^ { 0 } + \xi _ { 3 } ^ { 3 } ) + \xi \partial \rho _ { 0 3 , 0 3 } = 0 } } \\ { { R _ { 1 2 , 1 2 } } } & { { \equiv } } & { { 2 \rho _ { 1 2 , 1 2 } ( \xi _ { 1 } ^ { 1 } + \xi _ { 2 } ^ { 2 } ) + \xi \partial \rho _ { 1 2 , 1 2 } = 0 } } \\ { { R _ { 0 1 , 2 3 } } } & { { \equiv } } & { { \rho _ { 0 1 , 2 3 } ( \xi _ { 0 } ^ { 0 } + \xi _ { 1 } ^ { 1 } + \xi _ { 2 } ^ { 2 } + \xi _ { 3 } ^ { 3 } ) + \xi \partial \rho _ { 0 1 , 2 3 } = 0 } } \end{array} \right. $$ Now, the coefficients of the metric are rational functions in $K$ and the various geometric objects appearing in $\omega \gamma \rho$ can be obtained through the rules of differential algebra. It is thus possible to obtain the 13 non-zero components of the Riemann tensor for the Kerr metric according to K. R. Koehler in (http://kias.dyndns.org/crg/blackhole.html) by adding factorizations as follows: We shall distinguish the 6 non-vanishing components: $$ \left\{ \begin{array}{rcl} \rho_{01,01} & = & -\frac{mr(2(r^2-mr+a^2)+a^2(1-c^2))(r^2-3a^2c^2)}{2(r^2+a^2c^2)^3(r^2-mr+a^2)} \\ \rho_{02,02} & = & \frac{mr(r^2-mr+a^2+2a^2(1-c^2))(r^2-3a^2c^2)}{2(1-c^2)(r^2+a^2c^2)^3} \\ \rho_{03,03} & = & \frac{mr(1-c^2)(r^2-mr+a^2)(r^2-3a^2c^2)}{2(r^2+a^2c^2)^3} \\ \rho_{12,12} & = & -\frac{mr(r^2-3a^2c^2)}{2(1-c^2)(r^2+a^2c^2)(r^2-mr+a^2)} \\ \rho_{13,13} & = & -\frac{(1-c^2)mr(r^4-2a^2c^2r^2+4a^2r^2-2a^4c^2+3a^4-2a^2mr(1-c^2))(r^2-3a^2c^2)}{2(r^2+a^2c^2)^3(r^2-mr+a^2)} \\ \rho_{23,23} & = & \frac{mr(2r^4-a^2c^2r^2+5a^2r^2-a^4c^2+3a^4-a^2mr(1-c^2))(r^2-3a^2c^2)}{2(r^2+a^2c^2)^3} \end{array} \right. $$ from the 7 components that are vanishing when $a = 0$ . $$ \left\{ \begin{array}{rcl} \rho_{01,23} & = & \frac{amc(2r^2-a^2c^2+3a^2)(3r^2-a^2c^2)}{2(r^2+a^2c^2)^3} \\ \rho_{02,31} & = & -\frac{amc(r^2-2a^2c^2+3a^2)(3r^2-a^2c^2)}{2(r^2+a^2c^2)^3} \\ \rho_{03,12} & = & -\frac{amc(3r^2-a^2c^2)}{2(r^2+a^2c^2)^2} \\ \rho_{02,10} & = & \frac{3a^2mc(3r^2-a^2c^2)}{2(r^2+a^2c^2)^3} \\ \rho_{02,32} & = & \frac{amr(3r^2-mr+3a^2)(r^2-3a^2c^2)}{2(r^2+a^2c^2)^3} \\ \rho_{13,23} & = & -\frac{3a^2mc(1-c^2)(r^2+a^2)(3r^2-a^2c^2)}{2(r^2+a^2c^2)^3} \\ \rho_{01,13} & = & \frac{amr(1-c^2)(3r^2+3a^2-2mr)(r^2-3a^2c^2)}{2(r^2+a^2c^2)^3(r^2-mr+a^2)} \end{array} \right. $$ We have finally to add the 8 vanishing components: $$ \begin{aligned} \rho_{01,03} &= 0, & \rho_{01,12} &= 0, & \rho_{02,03} &= 0, & \rho_{02,12} &= 0, \\ \rho_{03,13} &= 0, & \rho_{03,23} &= 0, & \rho_{12,13} &= 0, & \rho_{12,23} &= 0. \end{aligned} $$ In fact, as the Riemann tensor has $n ^ { 2 } ( n ^ { 2 } - 1 ) / 1 2$ components, that is 20 when $n = 4$ , we have to take into account the only identity: $$ \rho _ { 0 1 , 2 3 } + \rho _ { 0 2 , 3 1 } + \rho _ { 0 3 , 1 2 } = 0 \Rightarrow R _ { 0 1 , 2 3 } + R _ { 0 2 , 3 1 } + R _ { 0 3 , 1 2 } = 0 $$ We obtain therefore $\xi \partial ( \rho _ { 1 2 , 1 2 } / ( \omega _ { 1 1 } \omega _ { 2 2 } ) ) = 0$ but we have also $\xi \partial ( \rho _ { 0 3 , 0 3 } \rho _ { 1 2 , 1 2 } / d e t ( \omega ) ) = 0$ . The following invariants are obtained successively in a coherent way: $$ \omega _ { 1 1 } \omega _ { 2 2 } = { \frac { ( r ^ { 2 } + a ^ { 2 } c ^ { 2 } ) ^ { 2 } } { ( 1 - c ^ { 2 } ) ( r ^ { 2 } - m r + a ^ { 2 } ) } } \Rightarrow \mid \rho _ { 1 2 , 1 2 } \mid / ( \omega _ { 1 1 } \omega _ { 2 2 } ) = { \frac { m r ( r ^ { 2 } - 3 a ^ { 2 } c ^ { 2 } ) } { 2 ( r ^ { 2 } + a ^ { 2 } c ^ { 2 } ) ^ { 3 } } } $$ Also, as $a \in K$ , then $\rho _ { 0 1 , 2 3 }$ and $\rho _ { 0 2 , 1 3 }$ can be both divided by $a$ and we get the new invariant: $$ \textstyle \rho _ { 0 1 , 2 3 } / \rho _ { 0 3 , 1 2 } = \frac { 2 r ^ { 2 } - a ^ { 2 } c ^ { 2 } + 3 a ^ { 2 } } { r ^ { 2 } + a ^ { 2 } c ^ { 2 } } $$ These results are leading to $\xi ^ { 1 } = 0 , \xi ^ { 2 } = 0$ , thus to $\xi _ { 1 } ^ { 1 } = 0 , \xi _ { 2 } ^ { 2 } = 0$ and $\xi _ { 0 } ^ { 0 } + \xi _ { 3 } ^ { 3 } = 0$ after substitution in the equations defining the first order symbol $g _ { 1 }$ of $R _ { 1 }$ . In the case of the S-metric with $a = 0$ , the previous division has no meaning and we have only $\xi _ { 1 } = 0$ as the only equation of zero order. Let us now introduce the new equation: $$ R _ { 0 1 , 1 3 } \equiv \rho _ { 0 1 , 1 3 } ( \xi _ { 0 } ^ { 0 } + 2 \xi _ { 1 } ^ { 1 } + \xi _ { 3 } ^ { 3 } ) - \rho _ { 1 3 , 1 3 } \xi _ { 0 } ^ { 3 } - \rho _ { 0 1 , 0 1 } \xi _ { 3 } ^ { 0 } + ( \rho _ { 0 1 , 2 3 } + \rho _ { 0 2 , 1 3 } ) \xi _ { 1 } ^ { 2 } = 0 $$ As we have $\xi _ { 0 } ^ { 0 } + \xi _ { 3 } ^ { 3 } = 0$ and $\xi _ { 1 } ^ { 1 } = 0$ , we obtain therefore a linear equation of the form: $$ ~ \rho _ { 1 3 , 1 3 } \xi _ { 0 } ^ { 3 } + \rho _ { 0 1 , 0 1 } \xi _ { 3 } ^ { 0 } - ( \rho _ { 0 1 , 2 3 } + \rho _ { 0 2 , 1 3 } ) \xi _ { 1 } ^ { 2 } = 0 ~ $$ Similarly, we have also: $$ R _ { 0 1 , 0 2 } \equiv \rho _ { 0 1 , 0 2 } ( 2 \xi _ { 0 } ^ { 0 } + \xi _ { 1 } ^ { 1 } + \xi _ { 2 } ^ { 2 } ) - ( \rho _ { 0 1 , 2 3 } + \rho _ { 0 2 , 1 3 } ) \xi _ { 0 } ^ { 3 } + \rho _ { 0 1 , 0 1 } \xi _ { 2 } ^ { 1 } + \rho _ { 0 2 , 0 2 } \xi _ { 1 } ^ { 2 } = 0 $$ and we obtain therefore a linear equation of the form: $$ 2 \rho _ { 0 1 , 0 2 } \xi _ { 0 } ^ { 0 } - ( \rho _ { 0 1 , 2 3 } + \rho _ { 0 2 , 1 3 } ) \xi _ { 0 } ^ { 3 } + \rho _ { 0 1 , 0 1 } \xi _ { 2 } ^ { 1 } + \rho _ { 0 2 , 0 2 } \xi _ { 1 } ^ { 2 } = 0 $$ In the case of the S-metric, that is when $a = 0$ , we obtain respectively $\xi _ { 3 } ^ { 0 } = 0$ and $\xi _ { 2 } ^ { 1 } = 0$ as in \[18\] because $\xi _ { 0 } ^ { 0 } \simeq \xi _ { 3 } ^ { 0 }$ The previous linear system has thus a rank equal to 2 and we obtain therefore because $\xi _ { 0 } ^ { 3 } \simeq \xi _ { 3 } ^ { 0 } , \xi _ { 1 } ^ { 2 } \simeq \xi _ { 2 } ^ { 1 }$ . $$ \xi _ { 3 } ^ { 0 } = 0 , \xi _ { 2 } ^ { 1 } = 0 , \Leftrightarrow \xi _ { 0 } ^ { 3 } = 0 , \xi _ { 1 } ^ { 2 } = 0 , \xi _ { 0 } ^ { 0 } = 0 , \xi _ { 3 } ^ { 3 } = 0 $$ It remains to study the following 4 linear equations, namely \[13\]: $$ R _ { 0 1 , 0 3 } = 0 , \ R _ { 0 3 , 2 3 } = 0 , \ R _ { 0 3 , 1 3 } = 0 , \ R _ { 0 2 , 0 3 } = 0 $$ The rank of the previous system with respect to the 4 jet coordinates $( \xi _ { 0 } ^ { 1 } , \xi _ { 0 } ^ { 2 } , \xi _ { 3 } ^ { 1 } , \xi _ { 3 } ^ { 2 } )$ is equal to 2, for both the S and K-metrics thanks to the two striking identities: $$ \begin{array} { r } { R _ { 0 3 , 1 3 } + a ( 1 - c ^ { 2 } ) R _ { 0 1 , 0 3 } = 0 , \quad R _ { 0 2 , 0 3 } + \frac { a } { ( r ^ { 2 } + a ^ { 2 } ) } R _ { 0 3 , 2 3 } = 0 } \end{array} $$ Two prolongations only provide 6 additional equations of order one that we provide in the following list which is obtained $m o d ( j _ { 2 } ( \Omega ) )$ , namely: $$ \xi ^ { 1 } = 0 , \xi ^ { 2 } = 0 , \xi _ { 2 } ^ { 1 } = 0 , \xi _ { 3 } ^ { 0 } = 0 , \xi _ { 3 } ^ { 1 } + l i n ( \xi _ { 0 } ^ { 1 } , \xi _ { 0 } ^ { 2 } ) = 0 , \xi _ { 3 } ^ { 2 } + l i n ( \xi _ { 0 } ^ { 1 } , \xi _ { 0 } ^ { 2 } ) $$ We have therefore obtained the inclusion of Lie algebroids $R _ { 1 } ^ { ( 2 ) } \subset R _ { 1 } ^ { ( 1 ) } = R _ { 1 } \subset J _ { 1 } ( T )$ with respective dimensions $4 < 1 0 = 1 0 < 2 0$ . Using the standard $k e r - c o k e r$ long exaxt sequence of prolongations: $$ 0 \to R_3 \to J_3(T) \to J_3(S_2 T^*) \to Q_2 \to 0, \quad 0 \to 4 \to 140 \to 150 \to 14 \to 0 $$ we discover that the initial Killing system for the Kerr metric has 14 compatibility conditions of second order contrary to the 20 existing for the Minkowski metric. Such a result has been obtained totally independently of any specific GR technical object like the Teukolski scalars or the Killing-Yano tensors introduced in \[14, 15\]. However, this system is not involutive because its symbol is finite type but non-zero. Using one more prolongation, all the sections (care again) vanish but $\xi ^ { 0 }$ and $\xi ^ { 3 }$ , a result leading to $d i m ( R _ { 1 } ^ { ( 3 ) } ) = 2$ in a coherent way with the only nonzero Killing vectors $\{ \partial _ { t } , \partial _ { \phi } \}$ . We have indeed: $$ \xi _ { 0 } ^ { 1 } = 0 ~ \xi _ { 0 } ^ { 2 } = 0 ~ \Big | \ \Leftrightarrow ~ \xi _ { 3 } ^ { 1 } = 0 , ~ \xi _ { 3 } ^ { 2 } = 0 ~ \Rightarrow ~ \xi _ { 1 } ^ { 0 } = 0 , ~ \xi _ { 1 } ^ { 0 } = 0 , ~ \xi _ { 2 } ^ { 0 } = 0 , ~ \xi _ { 2 } ^ { 3 } = 0 $$ Taking therefore into account that the metric only depends on ( $x ^ { 1 } = r , x ^ { 2 } = c o s ( \theta ) \rangle$ we obtain after three prolongations the inclusions of first order systems: $$ R _ { 1 } ^ { ( 3 ) } \subset R _ { 1 } ^ { ( 2 ) } \subset R _ { 1 } ^ { ( 1 ) } = R _ { 1 } \subset J _ { 1 } ( T ) \Leftrightarrow 2 < 4 < 1 0 = 1 0 < 2 0 $$ Surprisingly and contrary to the situation found for the $S$ metric, we have now an involutive first order system with only solutions $( \xi ^ { 0 } = c s t , \xi ^ { 1 } = 0 , \xi ^ { 2 } = 0 , \xi ^ { 3 } = c s t$ ) and notice that $R _ { 1 } ^ { ( 3 ) }$ does not depend any longer on the parameters $( m , a ) \in K$ . The difficulty is to know what second members must be used along the procedure met for all the motivating examples, in particular we have again identities to zero like $d _ { 0 } \xi ^ { 1 } - \xi _ { 0 } ^ { 1 } = 0 , d _ { 0 } \xi ^ { 2 } - \xi _ { 0 } ^ { 2 } = 0$ . We finally obtain 14 second order generating CC and their prolongations as we already said but also 6 third order CC coming from the 6 following components of the Spencer operator, namely: $$ \begin{aligned} d_1 \xi^1 - \xi_1^1 &= 0, & d_2 \xi^1 - \xi_2^1 &= 0, & d_3 \xi^1 - \xi_3^1 &= 0, \\ d_1 \xi^2 - \xi_1^2 &= 0, & d_2 \xi^2 - \xi_2^2 &= 0, & d_3 \xi^2 - \xi_3^2 &= 0. \end{aligned} $$ a result that cannot be even imagined from \[14\]. Of course, proceeding like in the motivating examples, we must substitute in the right members the values obtained from $j _ { 2 } ( \Omega )$ and set for example $\xi _ { 1 } ^ { 1 } = - { \textstyle \frac { 1 } { 2 \omega _ { 1 1 } } } \xi ^ { r } \partial _ { r } \omega _ { 1 1 }$ while replacing $\xi ^ { 1 }$ and $\xi ^ { 2 }$ by the corresponding linear combinations of the Riemann tensor already obtained for the right members of the two zero order equations. The corresponding Fundamental Diagram $I$ is no longer depending on $( m , a )$ as follows: $$ \begin{array}{rcccccccccccccl} &&&&& 0 &&0&&0&& 0 & &0& \\ &&&&& \downarrow && \downarrow && \downarrow & & \downarrow & & \downarrow & \\ & 0& \rightarrow& \Theta &\stackrel{j_1}{\rightarrow}& 2 &\stackrel{D_1}{\rightarrow}& 8 &\stackrel{D_2}{\rightarrow} & 12 &\stackrel{D_3}{\rightarrow}& 8 &\stackrel{D_4}{\rightarrow}& 2 &\rightarrow 0 \\ &&&&& \downarrow & & \downarrow & & \downarrow & & \downarrow &&\downarrow & \\ & 0 & \rightarrow & 4 & \stackrel{j_1}{\rightarrow} & 20 & \stackrel{D_1}{\rightarrow} & 40 &\stackrel{D_2}{\rightarrow} & 40 &\stackrel{D_3}{\rightarrow} & 20 &\stackrel{D_4}{\rightarrow} & 4 & \rightarrow 0 \\ & & & \parallel && \downarrow & & \downarrow & & \downarrow & & \downarrow & & \downarrow & \\ 0 \rightarrow & \Theta &\rightarrow & 4 & \stackrel{\cal{D}}{\rightarrow} & 18 & \stackrel{{\cal{D}}_1}{\rightarrow} & 32 & \stackrel{{\cal{D}}_2}{\rightarrow} & 28 & \stackrel{{\cal{D}}_3}{\rightarrow} & 12 &\stackrel{{\cal{D}}_4}{\rightarrow} & 2 & \rightarrow 0 \\ &&&&& \downarrow & & \downarrow & & \downarrow & & \downarrow & &\downarrow & \\ &&&&& 0 && 0 && 0 &&0 &&0 & \end{array} $$ with the Euler-Poincaré characteristic $4 - 1 8 + 3 2 - 2 8 + 1 2 - 2 = 0$ . However, the only intrinsic concepts associated with a differential sequence are the ” extension modules ” that only depend on the Kerr differential module but not on the differential sequence and we repeat once more that: THE ONLY IMPORTANT CONCEPT IS THE GROUP INVOLVED, NOT THE METRIC. Needless to say that the group involved in this case has no physical usefulness. EXAMPLE 4.3: (Schwarzschild metric) According to what we said, the situation of the Schwarzschild metric is much simpler when $a = 0$ with only six non-zero components of the Riemannn tensor: $$ \begin{array}{rcl} {\rho}_{01,01} & = & - \frac{m}{r^2(r-m)} \\ {\rho}_{02,02} & = & \frac{m}{2r} \\ {\rho}_{03,03} & = & \frac{m(1-c^2)}{2r} \\ {\rho}_{12,12} & = & - \frac{m}{2(r-m)} \\ {\rho}_{13,13} & = & - \frac{m(1-c^2)}{2(r-m)} \\ {\rho}_{23,23} & = & mr(1-c^2) \end{array} $$ We have now the inclusions of algebroids: $$ R _ { 1 } ^ { ( 3 ) } \subset R _ { 1 } ^ { ( 2 ) } \subset R _ { 1 } ^ { ( 1 ) } = R _ { 1 } \subset J _ { 1 } ( T ) 4 < 5 < 1 0 = 1 0 < 2 0 $$ The subsystem $R _ { 1 } ^ { ( 2 ) }$ is def by g the 5 nest rd eutios: $$ \mid \xi ^ { 1 } = 0 , \xi _ { 2 } ^ { 1 } = 0 , \xi _ { 3 } ^ { 1 } = 0 , \xi _ { 2 } ^ { 0 } = 0 , \xi _ { 3 } ^ { 0 } = 0 $$ while $R _ { 1 } ^ { ( 3 ) }$ is obtained by adding aain $\xi _ { 0 } ^ { \perp } = 0$ . We have 15 generating second order CC and 3 third order CC provided by the Spencer operator: $$ d _ { 1 } \xi ^ { 1 } - \xi _ { 1 } ^ { 1 } = 0 , d _ { 2 } \xi ^ { 1 } - \xi _ { 2 } ^ { 1 } = 0 , d _ { 3 } \xi ^ { 1 } - \xi _ { 3 } ^ { 1 } = 0 $$ in which we have to substitute $\xi ^ { 1 } \in j _ { 2 } ( \Omega )$ $, \xi _ { 1 } ^ { 1 } \in j _ { 2 } ( \Omega ) , \xi _ { 2 } ^ { 1 } \in j _ { 2 } ( \Omega ) , \xi _ { 3 } ^ { 1 } \in j _ { 2 } ( \Omega ) .$ It thus follows from the previous results of this section, which have been obtained without the need of any purely relativistic tool, that black holes cannot exist as they are contradicting \[14,15\] while showing that the underlying mathematical problem is a purely formal one with no link with GR along the given motivating example.  When a linear partial differential operator $\mathcal { D } \xi = \eta$ is given, a direct problem is to look for the generating compatibility conditions $\mathcal { D } _ { 1 } \eta = 0$ that must be satisfied by $\eta$ . Similarly, if $\mathcal { D } _ { 1 } \eta = \zeta$ is given, one may look for CC of the form $\mathcal { D } _ { 2 } \zeta = 0$ and so on. The mathematical community (and we do not speak about the physical community !) is of course aware of such a ” step by step” way but is not at all aware of the existence of another “as a whole” procedure allowing to define the various differential operators of the differential sequence thus obtained apart from the very specific situation of the Poincaré (in France !) sequence for the exterior derivative that admits a unique defining formula for each operator separately. The best known case is that of Riemannian geometry and its application to general relativity with the successive Killing, Riemann and Bianchi operators of first, second and first order respectively. In particular, we may ask ” Who knows about the Spencer operator and the corresponding Spencer sequence $\ell ^ { \ n }$ at the heart of this paper. In the Introduction, we have explained and illustrated through many motivating examples that, when a second order differential operator $\mathcal { D }$ is depending on constant or variable coefficients, its generating compatibility conditions (CC) may be of first, second, third and even sixth or higher order, a result largely depending on the parameters. In the meantime, we have shown that the solution of this problem for a system of order $q$ cannot be obtained without bringing such a system to an involutive form or at least to a formally integrable form of order $q + r$ after differentiating $r + s$ times the equations while keeping only the equations left up to order $q + r$ in such a way that the order of the CC is at most $r + s + 1$ . From a very different point of view, the Spencer differential sequence is obtained by bringing any involutive system $R _ { q } \subset J _ { q } ( E )$ to a first order involutive system $R _ { q + 1 } \subset J _ { 1 } ( R _ { q } )$ having an isomorphic space of solutions or, with a more precise language, allowing to define a differential module isomorphic to the differential module $M$ defined by the initial system. The quotient of the Spencer sequence for the first order trivially involutive first order system $J _ { q + 1 } ( E ) \subset J _ { 1 } ( J _ { q } ( E ))$ by the previous Spence sequence which is induced by the inclusion $R _ { q } \subset J _ { q } ( E )$ is the well defined finite length differential Janet sequence introduced by M. Janet as a footnote in 1920 which is thus providing another resolution of the same space of solutions or of the differential module $M$ already defined. According to a very difficult theorem of (differential) homological algebra, the only objects that do not depend of the resolution used are the (differential) extension modules that are measuring the fact that the corresponding dual sequence made by the respective formal adjoint of the operators involved and going thus ” backwards” (that is from right to left) may not be exact, that is each operator may not generate the CC of the preceding one. it thus follows that the Spencer and Janet sequences will bring the same formal informations as a whole, even though, in actual practice, we proved that they can be completely different. It may happen, for example with the Schwarzschild and Kerr metrics, with $q = 1 , r = 0 , s = 3$ , that the final corresponding FI systems will not depend any longer on the parameters involved initially. Accordingly, the only important object to consider is not the metric but its group $G$ of invariance which is used through the fact that the Spencer sequence is the tensor product of the Poincaré sequence by its Lie algebra $\vec { \mathcal { G } }$ , the main formal reason for which black holes cannot exist. We have thus finally proved that the main idea, along the tentative of H. Weyl in 1918, is not to shrink the dimension of this group from 10 down to to 4 or 2 parameters by using the S or K metrics instead of the M metric but, on the contrary, to enlarge the group from 10 up to 11 or 15 parameters by using the Weyl or conformal group instead of the Poincar’e group of space-time while using the adjoint of the respective Spencer sequences \[17\]. It will follow that the first set of Maxwell equations is obtained by a projection of the second Spencer operator $D _ { 2 }$ that can be parametrized by a projection of the first Spencer operator $D _ { 1 }$ , a result contradicting the basic assumption of classical gauge theory in which they are induced by $D _ { 3 }$ . The main problem today is that, in the minimum resolution of the conformal Killing operator, the generating CC of the second order Weyl operator are also made by a second order operator, a result confirmed by my PhD student A. Quadrat (INRIA) in 2016 (See \[24, 26\] or arXiv: 1603.05030) but still not acknowledged and showing that conformal geometry but me revisited by using the Spencer $\delta$ -cohomology \[31\]. Finally, in a more general framework, when a Lie group $G$ is acting on a manifold $X$ of dimension $n$ , the Spencer sequence is always locally and formally exact, being isomorphic to the tensor product of the Poincaré sequence by the Lie algebra of $G$ . On the contrary, the corresponding formally exact Janet sequence may have Janet bundles of high dimensions. The last operator $\mathcal { D } _ { n }$ is always surjective while its adjoint is always injective. However, $a d ( \mathcal { D } _ { n - 1 } )$ may not define all the CC of $a d ( \mathcal { D } _ { n } )$ because $a d ( D _ { n } )$ may fail to be injective. Applying these methods to the conformal group of transformations when $n = 4$ , we discovered in the very recent \[17\] the common conformal origin of the Cauchy, Cosserat, Clausius/Poisson and Maxwell equations. For example the EM field $F$ comes from the composition of epimorphisms ($60 \to 16 \to 6 \to 0$): $$ \hat{C}_1 \to \hat{C}_1 / \tilde{C}_1 = (T^* \otimes \hat{R}_2) / (T^* \otimes \tilde{R}_2) \simeq T^* \otimes (\hat{R}_2 / \tilde{R}_2) \simeq T^* \otimes \hat{g}_2 \simeq T^* \otimes T^* \xrightarrow{\delta} \wedge^2 T^* $$ while the EM potential comes from the composition of epimorphisms ( $1 5 4 0$ ): $$ { \hat { C } } _ { 0 } { \hat { C } } _ { 0 } / { \tilde { C } } _ { 0 } \simeq { \hat { R } } _ { 2 } / { \tilde { R } } _ { 2 } \simeq { \hat { g } } _ { 2 } \simeq T ^ { * } $$ and the parametrization $d A = F$ is induced by $D _ { 1 }$ while the Maxwell equation $d { \boldsymbol { F } } = 0$ is induced by $D _ { 2 }$ . Using the short exact splitting sequence $0 \to S _ { 2 } T ^ { * } { \overset { \delta } { \longrightarrow } } T ^ { * } \otimes T ^ { * } { \overset { \delta } { \longrightarrow } } \wedge ^ { 2 } T ^ { * } \to 0$ leading to the isomorphism $T ^ { * } \otimes T ^ { * } \simeq S _ { 2 } T ^ { * } \oplus { \wedge } ^ { 2 } T ^ { * } \simeq ( R _ { i j } ) \oplus ( F _ { i j } )$ with $1 6 = 1 0 + 6$ , that only depends on the elations of the conformal group, along the Fundamental Diagram $I I$ showing the common conformal origin of electromagnetism and gravitation, already published in 1983 \[30\]. More generally, taking the quotient of the Spencer sequence for $\hat { C } _ { r }$ by the Spencer sequence for $\tilde { \Theta }$ , we obtain a formally exact sequence with bundles $\hat { C } _ { r } / \tilde { C } _ { r }$ which is isomorphic to the tensor product of the Poincar\ ’e sequence by a vector space of dimension $d i m ( \hat { g } _ { 2 } ) = d i m ( T ^ { * } ) = 4$ that we can project onto the formally exact Poincaré sequence, with a shift by one step contradicting the use of $U ( 1 )$ in GT. We obtain therefore the following commutative and exact diagram in which we notice that the first vertical short exact sequence on the left is nothing else than the bottom short exact sequence of the previous fundamental diagram II: $$ \begin{array}{rcccccccccccccl} & & & & 0 & & 0 & & 0 & & 0 & & & & \\ & & & & \downarrow & & \downarrow & & \downarrow & & \downarrow & & & & \\ & & 1 & \rightarrow & 10 & \rightarrow & 20 & \rightarrow & 15 & \rightarrow & 4 & \rightarrow & 0 & gravitation \\ & & & & \downarrow & & \downarrow & & \downarrow & & \parallel & & & & \\ & & 4 & \stackrel{d}{ \rightarrow} & 16 & \stackrel{d}{\rightarrow} & 24 & \stackrel{d}{\rightarrow} & 16 & \stackrel{d}{\rightarrow} & 4 & \rightarrow & 0 & & \\ & & \parallel && \downarrow & & \downarrow & & \downarrow & & \downarrow & & & & \\ 1 & \stackrel{d}{\rightarrow} & 4 &\stackrel{d}{\rightarrow} & 6 & \stackrel{d}{\rightarrow} & 4 & \stackrel{d}{\rightarrow} & 1 & \rightarrow & 0 &\hspace{1cm} & & electromagnetism & \\ & & & & \downarrow & & \downarrow & & \downarrow & & & & & & \\ & & & & 0 & & 0 & & 0 & & & & & & \end{array} $$ A snake chase in this diagram proves that gravitation has to do with the conformal factor $\phi$ when $F = 0$ as the first operator on the upper left is just the composition $D _ { 1 } \circ d: \phi \to {\partial}_i \phi \to {\partial}_{ij} \phi$ according to the fundamental diagram II. Also, using the dual summation $\rho ( x ) ( \partial _ { 4 } \xi ^ { 4 } - \xi _ { 4 } ^ { 4 } ) + g ^ { i } ( x ) ( \partial _ { i } A -$ $A _ { i } ) + g ^ { i j } ( x ) ( \partial _ { i } A _ { j } - 0 )$ with $\xi _ { i j r } ^ { k } = 0$ and integrating by parts, we get the induction equations: $$ \xi ^ { 4 } \partial _ { 4 } \rho = 0 , \xi _ { 4 } ^ { 4 } = A \partial _ { i } g ^ { i } + \rho = 0 , \xi _ { 4 i } ^ { 4 } = A _ { i } \rightarrow \partial _ { r } g ^ { r i } + g ^ { i } = 0 $$ leading to the Poisson equation for $\phi$ when $g ^ { i j } = \phi \omega ^ { i j }$ and thus to a new conformal scheme for electromagnetism and gravitation. It is also a reason for which elations are also sometimes called and confused with “accelerations” because and accelerometer on an inertial platform in a rocket can only measure the difference $\partial _ { 4 } \xi _ { 4 } ^ { k } - \xi _ { 4 4 } ^ { k }$ with $k = 1 , 2 , 3$ between the acceleration (time derivative of the speed as a rotation of space-time) and the gravitational field, while one has $\xi _ { 4 4 } ^ { k } = \omega ^ { k r } A _ { r }$ and the ditive ruleof secn rderts $$ \bar{f} = g \circ f \Rightarrow {\bar{f}}^k_{44} = f^k_{44} + g ^k_{44} $$ holds whenever $f _ { 1 } = g _ { 1 } = i d _ { 1 } = j _ { 1 } ( i d )$ for the one jet of the identity transformation $y = x$ (Compare \[25\] to \[33\]). However, no one of these results could have been even imagined by Weyl because the Riemann operator must be replaced by the Weyl operator while the Bianchi operator must be replaced by a second order CC operator, a result still not known today in the conformal framework \[33, 34\]. In the spirit of the Cosserat brothers on elasticity in 1909 and H. Weyl on electromagnetism in 1918, we have tried to explain and illustrate why the structure of physical theories must depend on a Lie group of transformations considered as a Lie pseudogroup. The main idea is then to introduce the corresponding Janet (1920) or Spencer (1970) differential sequences both with their respective adjoint differential sequences needed for variational calculus \[28, 35\]. In a coherent way with the dream of Weyl, the most striking result is the “shift” that must be used in the physical interpretation of the sequences obtained, in the sense that the EM field as a 2-form is coming from the first Spencer bundle and not from the second (See \[16\] for other examples). These two sequences may be sometimes quite different though always homologically equivalent and we have tried to illustrate them through the study of gravitational waves and black holes. As both of them are largely unknown, no other reference can be found and only the future will decide about the best one that must be used !. After 1950 the introduction of homological algebra in pure mathematics has been more than an evolution but rather quite a revolution in algebraic geometry. Following the wish of one of the reviewers that we thank, our aim in this appendix is to explain that differential homological algebra, namely the way to add the word ” differential ” in front of ALL the concepts already introduced (module, resolution, duality, extension modules, ..) will bring quite a revolution in physics. The main problem is that a full account should need a whole paper \[35\] or even a whole book \[21\] and that we have seen in the core of this paper that the differential geometric framework is already difficult while the geometric part is somehow even more difficult. By chance, in a more general framework, we have proved that the controllability of a control system is a structural property of the control system, not depending on the choice of inputs and outputs among the control variables and can be described in terms of vanishing of extension modules \[4\]. Roughly, it means that a control system is controllable if and only if the system can be parametrized or, equivalently, if there does not exist any observable quantity which is solution of an autonomous system for itself, a more physical way to say that the differential module defined by the system is torsion-free \[35\]. Accordingly we shall explain through two examples how to pass ” beyond the mirror”. First of all, we sketch the idea hidden beyond ” differential extension modules” by using only the language of operators along the content of this paper. For this, we notice that, if a given linear differential operator $\mathcal { D } _ { 1 } \eta = \zeta$ is give, a ” direct problem ” is to look for compatibility conditions (CC) in the form of a linear operator $\mathcal { D } _ { 2 } \zeta = 0$ and so on. We hope that the reader is now convinced that it may already be a difficult problem to provide ” a prior” any idea of the order of the CC operator as we have provided examples in which the order is ranging from 1 up to 6. However, conversely, if a system described by $\mathcal { D } _ { 1 } \eta = 0$ , a much more difficult ” inverse problem ” is to find a parametrization, that is an operator $\mathcal { D } \xi = \eta$ having exactly such CC. Now, using the well known properties of the adjoint procedure $\mathcal { D } a d ( \mathcal { D } )$ with $a d ( a d ( \mathscr { D } ) ) = \mathscr { D }$ , if $\mathcal { D } \xi = \eta$ has generating CC $\mathcal { D } _ { 1 } \eta = 0$ , then $\mathcal { D } _ { 1 } \circ \mathcal { D } = 0 \Rightarrow a d ( \mathcal { D } ) \circ a d ( \mathcal { D } _ { 1 } ) = 0$ but $a d ( \mathcal { D } )$ may not generate all the CC of $a d ( \mathcal { D } _ { 1 } )$ while $a d ( D _ { 1 } )$ may not generate all the CC of $a d ( D _ { 2 } )$ and $a d ( \mathcal { D } _ { 2 } )$ may not generate all the CC of $a d ( \mathcal { D } _ { 3 } )$ and so on, as can be seen in the motivating examples. The introduction of extension modules is a way to measure such a gap. Starting with the system, we sketch on an example the five steps of the differential double duality test must be used in order to know whether this system can be parametrized or not: STEP 1: Write down the system $\mathcal { D } _ { 1 } \eta = \zeta$ : $d _ { 1 2 } \eta ^ { 1 } + d _ { 2 2 } \eta ^ { 2 } = \zeta$ . STEP 2: Write down the adjoint operator $\cdots$ backwards: multiplying on the left by a test function $\lambda$ and integrating by parts, we get $$ {\eta}^1 \rightarrow d_{12} \lambda = {\mu}^1, {\eta}^2 \rightarrow d_{22} \lambda = {\mu}^2 $$ STEP 3: As any operator is the adjoint of its adjoint, find its CC as an operator $a d ( \mathcal { D } ) \mu = \nu$ $d _ { 2 } \mu ^ { 1 } - d _ { 1 } \mu ^ { 2 } = \nu$ . STEP 4: Write down its adjoint as an operator $\mathcal { D } \xi = \eta$ . $\mu ^ { 1 } \to - d _ { 2 } \xi = \eta ^ { 1 } , \mu ^ { 2 } \to d _ { 1 } \xi = \eta ^ { 2 }$ Step 5: Find its CC as an operator $\mathcal { D } _ { 1 } ^ { \prime } \eta = 0$ and compare to $\mathcal { D } _ { 1 }$ .\ CONCLUDE: If $\mathcal { D } _ { 1 }$ and $\mathcal { D } _ { 1 } ^ { \prime }$ are identical or have the same space of solutions, that is produce the same involutive system, then $\mathcal { D } _ { 1 }$ is parametrized by $\mathcal { D }$ . Otherwise, if the space of solutions of $\mathcal { D } _ { 1 } ^ { \prime }$ is smaller than that of $\mathcal { D } _ { 1 }$ , then $\mathcal { D } _ { 1 }$ cannot be parametrized. It means that there is at least one new CC which is an autonomous element which is satisfying at least one OD or PD equation for itself, a fact showing that the initial system, considered as a control system, cannot be controllable (See \[21, 35\] or Zbl 179.93001 for more details and examples, in particular the study of RCL electrical circuits). We notice that we don’t need to separate $( \eta ^ { 1 } , \eta ^ { 2 } )$ into input and output. END: In the present example, $\mathcal { D } _ { 1 } ^ { \prime }$ is defined by $d _ { 1 } \eta ^ { 1 } + d _ { 2 } \eta ^ { 2 } = \zeta ^ { \prime }$ and we have $d _ { 2 } \zeta ^ { \prime } = \zeta$ , that is $d _ { 2 } \zeta ^ { \prime } = 0$ whenever $\zeta = 0$ . It follows that $\mathcal { D }$ CANNOT be parametrized. Ofc course, such a result could have been found directly in this elementary example but such a test is unavoidable in general. In order to recapitulate, our purpose in this appendix is to explain how differential homological algebra is able to combine these two procedures. For this we ask the reader to spend a few dollars in order to realize the next experiment and try to understand why such a fact is clearly showing that gravitational waves cannot exist as we have explained in the Introduction and in theorem 4.1 or in \[36\].\ Example 5.1: (Double pendulum) Let us consider a thin rigid bar of length $L$ able to slide horizontally with a position $x$ function of time $t$ and time derivative $d$ . We may attach a pendulum of length $l _ { 1 }$ , mass $m _ { 1 }$ , (small) angle $\theta _ { 1 }$ from the vertical at one end and a pendulum of length $l _ { 2 }$ , mass $m _ { 2 }$ , (small) angle $\theta _ { 2 }$ from the vertical at the other end. At equilibrium, the two pendulums are fixed. If one of the pendulums is slightly moved and $l _ { 1 } \neq l _ { 2 }$ , moving the bar with enough skill, one can bring the full system to equilibrium, that is $x = c s t , \theta _ { 1 } = \theta _ { 2 } = 0$ and we say that the full system is controllable by using the movement of the bar. Now, if $l _ { 1 } = l _ { 2 } = l$ and one pendulum is moved while the other is untouched, then any way to stop the first wll bring the other to have the same movement and the system is not controllable that is, in particular, there is no way to stop both pendulums as before. In order to write down the two OD equations of the system, one has to consider the Newton law applied to the inertial forces $d ^ { 2 } ( x + l \theta )$ , the gravity force $m g$ and the tension of the thread. Projecting on the perpendicular to each pendulum to eliminate the last one, we may divide by the respective mass in order to get the two OD equations in the form of a surjective operator $\mathcal { D } _ { 1 } \eta = \zeta$ with $\eta = ( \eta ^ { 1 } = \theta ^ { 1 } , \eta ^ { 2 } = \theta ^ { 2 } , \eta ^ { 3 } = x$ ): $$ d ^ { 2 } x + l _ { 1 } d ^ { 2 } \theta ^ { 1 } + g \theta ^ { 1 } = \zeta ^ { 1 } , \qquad d ^ { 2 } x + l _ { 2 } d ^ { 2 } \theta ^ { 2 } + g \theta ^ { 2 } = \zeta ^ { 2 } $$ Multiplying the firstOD equation by the test function $\lambda ^ { 1 }$ , the seond by the test fncin $\lambda ^ { 2 }$ , adding and integrating by parts, we obtain the adjoint operator in the form $a d ( \mathcal { D } _ { 1 } ) \lambda = \mu$ , namely: $$ \theta ^ { 2 } l _ { 2 } d ^ { 2 } \lambda ^ { 2 } + g \lambda ^ { 2 } = \mu ^ { 2 } , x d ^ { 2 } \lambda ^ { 1 } + d ^ { 2 } \lambda ^ { 2 } = \mu ^ { 3 } $$ First of all, if $\mu ^ { 1 } = \mu ^ { 2 } = \mu ^ { 3 } = 0$ , we obtain two zero order OD equations where the second one is obtain by differentiating twice the first and substituting: $$ l _ { 2 } \lambda ^ { 1 } + l _ { 1 } \lambda ^ { 2 } = 0 , \qquad \quad ( l _ { 2 } / l _ { 1 } ) \lambda ^ { 1 } + ( l _ { 1 } / l _ { 2 } ) \lambda ^ { 2 } = 0 $$ As the determinant of this linear system is equal to $l _ { 1 } - l _ { 2 }$ , it follows that $a d ( \mathcal { D } _ { 1 } )$ is injective if and only if $l _ { 1 } \neq l _ { 2 }$ . Let us now prove that $\mathcal { D } _ { 1 }$ can be parametrized if and only if such a condition is satified. First of all we notice that $\lambda ^ { 1 } , \lambda ^ { 2 } \in j _ { 2 } ( \mu )$ and, using a tricky PP procedure or computer algebra, we discover that $a d ( { \mathcal { D } _ { 1 } } )$ has only one CC $a d ( \mathcal { D } ) \mu = \nu$ of order 4. Taking the adjoint, we obtain a well defined fourth order injective parametrization in the form $\mathcal { D } \xi = \eta$ as follows: $$ - l_1 l_2 d^4\xi - g(l_1+l_2)d^2 \xi - g^2 \xi = x, \, \, l_2 d^4\xi +g d^2 \xi = {\theta}^1, \, \, l_1 d^4 \xi + gd^2 \xi = {\theta}^2 $$ This parametrization is injective and we have the short exact adjoint sequences: $$ \begin{array}{ccccccccc} 0 & \longrightarrow & 1 & \stackrel{{\cal{D}}}{\longrightarrow} & 3 & \stackrel{{\cal{D}}_1}{\longrightarrow} & 2 & \longrightarrow & 0 \\ & & & & & & & & \\ 0 & \longleftarrow & 1 & \stackrel{ad({\cal{D}})}{\longleftarrow} & 3 & \stackrel{ad({\cal{D}}_1)}{\longleftarrow} & 2 & \longleftarrow & 0 \end{array} $$ At no moment we have ever been speaking about inputs, outputs and Kalman test !. It remains to study the case $l _ { 1 } = l _ { 2 } = l$ . Subtracting the second equation from the first while setting $\theta = \theta ^ { 1 } - \theta ^ { 2 }$ , we get $l d ^ { 2 } \theta = g \theta = 0$ , that is the standard Od equation of a single pendulum. Of course, we notice at once that the observable $\theta$ is an autonomous element and its behaviour cannot be controlled. Equivalently, $\theta$ is generating the torsion submodule $t ( M )$ of the differential module $M$ defined by $\mathcal { D } _ { 1 }$ . The main fact is that now $\mathcal { D } _ { 1 }$ cannot any longer be parametrized?. For this, we may multiply each OD equation by two test functions as before and integrate by parts in order to obtain: $$ l d ^ { 2 } \lambda ^ { 1 } + g \lambda ^ { 1 } = \mu ^ { 1 } , l d ^ { 2 } \lambda ^ { 2 } + g \lambda ^ { 2 } = \mu ^ { 2 } , d ^ { 2 } ( \lambda ^ { 1 } + \lambda ^ { 2 } ) = \mu ^ { 3 } $$ Summing the two first OD equations while setting $\lambda = \lambda ^ { 1 } + \lambda ^ { 2 }$ , we get: $$ d ^ { 2 } \lambda = \mu ^ { 3 } , \ l d ^ { 2 } \lambda + g \lambda = \mu ^ { 1 } + \mu ^ { 2 } $$ It follows that $a d ( \mathcal { D } _ { 1 } )$ is no longer injective with kernel $\lambda = 0 \Leftrightarrow \lambda ^ { 1 } + \lambda ^ { 2 } = 0$ . We also obtain $g \lambda = \mu ^ { 1 } + \mu ^ { 2 } - l \mu ^ { 3 }$ and, substituting, we obtain one single CC of order 2 for $a d ( \mathcal { D } _ { 1 } )$ , namely $a d ( \mathcal { D } )$ as follows: $$ d ^ { 2 } ( \mu ^ { 1 } + \mu ^ { 2 } - l \mu ^ { 3 } ) - g \mu ^ { 3 } = \nu $$ Multiplying by the test function $\xi$ and integrating by parts, we get $\mathcal { D }$ : $$ \mu ^ { 1 } d ^ { 2 } \xi = \theta ^ { 1 } , \mu ^ { 2 } d ^ { 2 } \xi = \theta ^ { 2 } , \mu ^ { 3 } - l d ^ { 2 } \xi - g \xi = x $$ The 3 CC are described by $\mathcal { D } _ { 1 } ^ { \prime } \neq \mathcal { D } _ { 1 }$ as follows: $$ d ^ { 2 } x + l d ^ { 2 } \theta ^ { 1 } + g \theta ^ { 1 } = 0 , d ^ { 2 } x + l d ^ { 2 } \theta ^ { 2 } + g \theta ^ { 2 } = 0 , \theta ^ { 1 } - \theta ^ { 2 } = 0 $$ We recover the fact that the new CC is introducing $\theta = \theta ^ { 1 } - \theta ^ { 2 }$ which is indeed a torsion element with autonomous OD equation $d ^ { 2 } \theta + g \theta = 0$ . This operator $\mathcal { D }$ of order 2 is injective and we have the short exact adjoint sequences: $$ \begin{array}{ccccccccc} & & & & & & 3 & & \\ & & & & & {\cal{D}}'_1 \nearrow \hspace{5mm} & & & \\ 0 & \longrightarrow & 1 & \stackrel{{\cal{D}}}{\longrightarrow} & 3 & \stackrel{{\cal{D}}_1}{\longrightarrow} & 2 & \longrightarrow & 0 \\ & & & & & & & & \\ 0 & \longleftarrow & 1 & \stackrel{ad({\cal{D}})}{\longleftarrow} & 3 & \stackrel{ad({\cal{D}}_1)}{\longleftarrow} & 2 & & \end{array} $$ It is absolutely impossible to recover all these results and understand their meaning without differential double duality and differential homological algebra. Example 5.2: (General relativity) We now prove that the above situation is .. exactly similar to the one we have described in GR in order to explain why GW cannot exist. When $n = 4$ , let us prove that the 10 Einstein equations cannot admit a potential or, equivalently, that the Ricci operator cannot be parametrized and that the Einstein operator is therefore useless as it cannot have a mathematical origin contrary to the Ricci operator, according to the fundamental diagram II. We shall discover that the existence of autonomous elements (torsion elements in the module framework) … are nothing else than the $2 0 - 1 0 = 1 0$ components of the Weyl bundle already introduced in the fundamental diagram II, a result highly not evident at first sight !. First of all, according to the parametrization test, to the result of the Introduction and to Theorem 4.1, we have at once the negative test procedure: $$ \begin{array}{ccccccc} & & & & & \stackrel{Riemann }{ } & 20 \\ & & & & & \nearrow & \\ & & 4 & \stackrel{Killing}{\longrightarrow} & 10 & \stackrel{Ricci}{\longrightarrow} & 10 \\ & & & & & & \\ 0 & \longleftarrow & 4 & \stackrel{Cauchy}{\longleftarrow} & 10 & \stackrel{ad(Ricci)}{\longleftarrow} & 10 \end{array} $$ The 5 steps are indicated as follows while taking into account that $C a u c h y = a d ( K i l l i n g )$ , a result that can be found in any textbook of continuum mechanics or elasticity: $R i c c i \longrightarrow a d ( R i c c i ) \longrightarrow C a u c h y \longrightarrow K i l l i n g \longrightarrow R i e m a n n$ Accordingly and in a coherent way with the previous example, the Ricci operator and thus the Einstein operator cannot be parametrized and we have $2 0 - 1 0 = 1 0$ generating torsion elements, that is to say elements satisfying at least one PD equation for each one. Of course, such an approach has never been followed because the fundamental diagram II is not known though it provides the splitting $R i e m a n n = R i c c i \oplus W e y l$ coming from the splitting $T ^ { * } \otimes T ^ { * } = S _ { 2 } T ^ { * } \oplus \wedge ^ { 2 } T ^ { * }$ . My chance has been to know very well the results of A. Lichnerowicz who has been my advisor during 25 years from 1973 to his death in 1998. Indeed, I got in mind the so-called (in France !) ” Lichnerowicz waves ” and I am now able to give a short proof of the linearized case \[18\]:\ Theorem 5.3: The Dalembertian of each component of the Weyl tensor is a linear differential consequence of the 10 Einstein equations $R _ { i j } = 0$ . Proof. Linearizing over the Euclide or Minkowski metric $\omega$ while using the Bianchi identities and taking into account that $R _ { k l , i j } = R _ { i j , k l }$ , we get with $d ^ { r } = \omega ^ { r i } d _ { i }$ : $$ d _ { r } R _ { k l , i j } + d _ { i } R _ { k l , j r } + d _ { j } R _ { k l , r i } = 0 \Rightarrow d ^ { r } R _ { r l , i j } = d ^ { r } R _ { i j , r l } = d _ { i } R _ { l j } - d _ { j } R _ { l i } $$ $$ d_{rs} R_{kl,ij} + d_{is} R_{kl,jr} + d_{js} R_{kl,ri} = 0 \Rightarrow \Box R_{kl,ij} + \omega^{rs} d_{is} R_{kl,jr} + \omega^{rs} d_{js} R_{kl,ri} = 0 $$ $$ \Box R_{kl,ij} = d_i (d_k R_{lj} - d_l R_{kj}) - d_j (d_k R_{li} - d_l R_{ki}) $$ Finally, using the splitting formula for defining the components $\sigma _ { l , i j } ^ { k } = \rho _ { l , i j } ^ { k } - ( \sum \rho _ { r s } )$ of the Weyl tensor from the components of the Riemann tensor with now $\sigma _ { l , r j } ^ { r } = 0$ , we obtain by linearization $\begin{array} { r } { \Sigma _ { l , i j } ^ { k } = R _ { l , i j } ^ { k } - ( \sum R _ { r s } ) } \end{array}$ It is absolutely impossible to extend the above result to the conformal Riemannian geometry because the analogue of the Bianchi operator is now of order 2, a result still neither known and nor even acknowledged today in the long exact sequence that we have already provided at the end of Section 3C, where $\mathcal { D }$ is the conformal Killing operator and $\mathcal { D } _ { 1 }$ is the Weyl operator \[15, 26, 31\]. This is the main reason for which the Spencer sequences and their adjoint sequences must be used !. 

Generating HTML Viewer...