_{1}

^{*}

In 1916, F.S. Macaulay developed specific localization techniques for dealing with “unmixed polynomial ideals” in commutative algebra, transforming them into what he called “inverse systems” of partial differential equations. In 1970, D.C. Spencer and coworkers studied the formal theory of such systems, using methods of homological algebra that were giving rise to “differential homological algebra”, replacing unmixed polynomial ideals by “pure differential modules”. The use of “differential extension modules” and “differential double duality” is essential for such a purpose. In particular, 0-pure differential modules are torsion-free and admit an “absolute parametrization” by means of arbitrary potential like functions. In 2012, we have been able to extend this result to arbitrary pure differential modules, introducing a “relative parametrization” where the potentials should satisfy compatible “differential constraints”. We recently noticed that General Relativity is just a way to parametrize the Cauchy stress equations by means of the formal adjoint of the Ricci operator in order to obtain a “minimum parametrization” by adding sufficiently many compatible differential constraints, exactly like the Lorenz condition in electromagnetism. In order to make this difficult paper rather self-contained, these unusual purely mathematical results are illustrated by many explicit examples, two of them dealing with contact transformations, and even strengthening the comments we recently provided on the mathematical foundations of General Relativity and Gauge Theory. They also bring additional doubts on the origin and existence of gravitational waves.

The main purpose of this paper is to prove how apparently totally abstract mathematical tools, ranging among the most difficult ones existing in differential geometry and homological algebra today, can also become useful and enlighten many engineering or physical concepts (see the review Zbl 1079.93001 for the only application to control theory).

In the second section, we first sketch and then recall, with more specific references, the main (difficult) mathematical results on differential extension modules and differential double duality that are absolutely needed in order to understand the purity concept and, in particular, the so-called purity filtration of a differential module ( [

The third section will present for the first time a few explicit motivating academic examples in order to illustrate the above mathematical results, in particular the unexpected striking situations met in the study of contact and unimodular contact structures.

In the fourth section, we finally provide examples of applications, studying the mathematical foundations of OD/PD control theory (CT) ( [

Let D = K [ d 1 , ⋯ , d n ] = K [ d ] be the ring of differential operators with coefficients in a differential field K of characteristic zero, that is such that ℚ ⊂ K , with n commuting derivations ∂ 1 , ⋯ , ∂ n and commutation relations d i a = a d i + ∂ i a , ∀ a ∈ K . If y 1 , ⋯ , y m are m differential indeterminates, we may identify D y 1 + ⋯ + D y m = D y with D m and consider the finitely presented left differential module M = M D with presentation D p → D m → M → 0 determined by a given linear multidimensional system with n independent variables, m unknowns and p equations. Applying the functor h o m D ( • , D ) , we get the exact sequence 0 → h o m D ( M , D ) → D m → D p → N D → 0 of right differential modules that can be transformed by a side-changing functor to an exact sequence of finitely generated left differential modules. This new presentation corresponds to the formal adjoint a d ( D ) of the linear differential operator D determined by the initial presentation but now with p unknowns and m equations, obtaining therefore a new finitely generated left differential module N = N D and we may consider h o m D ( M , D ) as the module of equations of the compatibility conditions (CC) of a d ( D ) , a result which is not evident at first sight (see [

It is however essential to notice that such an approach is leading to a “vicious circle” because the only constructive way to check whether M is torsion-free or not is to use the differential double duality. For this reason, we briefly recall the five steps of the corresponding test allowing to know whether a given differential system or operator can be parametrized or not:

• STEP 1: Write down the system in the form of a differential operator D 1 .

• STEP 2: Construct its formal adjoint a d ( D 1 ) .

• STEP 3: Construct the generating CC of such an operator as a new operator and use the fact that a d ( a d ( P ) ) = P , ∀ P ∈ D in order to denote it by a d ( D ) .

• STEP 4: Write down D = a d ( a d ( D ) ) .

• STEP 5: As a d ( D ) ∘ a d ( D 1 ) = a d ( D 1 ∘ D ) = a d ( 0 ) = 0 ⇔ D 1 ∘ D = 0 , we just need to check whether D 1 generates the CC of D or not.

If “yes”, we shall say that D 1 is parametrized by D .

If “no”, we shall say that D 1 cannot be parametrized.

The purpose of this paper is to extend such a result to a much more general situation, that is when M is not torsion-free, by using unexpected results first found by F.S. Macaulay in 1916 ( [

Introducing t r ( M ) = { m ∈ M | c d ( D m ) > r } where the codimension of Dm is n minus the dimension of the characteristic variety determined by m in the corresponding system for one unknown, we may define the purity filtration as in ( [

0 = t n ( M ) ⊆ t n − 1 ( M ) ⊆ ⋯ ⊆ t 1 ( M ) ⊆ t 0 ( M ) = t ( M ) ⊆ M

The module M is said to be r-pure if t r ( M ) = 0 , t r − 1 ( M ) = M or, equivalently, if c d ( M ) = c d ( N ) = r , ∀ N ⊂ M and a torsion-free module is a 0-pure module. Moreover, when K = k = c s t ( K ) is a field of constants and m = 1 , a pure module is unmixed in the sense of Macaulay, that is defined by an ideal having an equidimensional primary decomposition.

Example 2.1: As an elementary example with K = k = ℚ , m = 1 , n = 2 , p = 2 , the differential module defined by d 22 y = 0 , d 12 y = 0 is not pure because z ′ = d 2 y satisfies d 2 z ′ = 0 , d 1 z ′ = 0 while z ″ = d 1 y only satisfies d 2 z ″ = 0 and ( ( χ 2 ) 2 , χ 1 χ 2 ) = ( χ 2 ) ∩ ( χ 1 , χ 2 ) 2 . We obtain therefore the purity filtration 0 = t 2 ( M ) ⊂ t 1 ( M ) ⊂ t 0 ( M ) = t ( M ) = M with strict inclusions as 0 ≠ z ′ ∈ t 1 ( M ) while z ″ ∈ t 0 ( M ) but z ″ ∉ t 1 ( M ) .

From the few (difficult) references ( [

0 → e x t 1 ( N ) → M → e x t 0 ( e x t 0 ( M ) ) → e x t 2 ( N ) → 0

allows to test the torsion-free property of M in actual practice by using the double-duality formula t ( M ) = e x t 1 ( N ) as in ( [

Independently of the previous results, the following procedure, where one may have to change linearly the independent variables if necessary, is the heart towards the next effective definition of involution. It is intrinsic even though it must be checked in a particular coordinate system called δ -regular ( [

• Equations of class n: Solve the maximum number β q n of equations with respect to the jets of order q and class n. Then call ( x 1 , ⋯ , x n ) multiplicative variables.

• Equations of class i ≥ 1 : Solve the maximum number β q i of remaining equations with respect to the jets of order q and class i. Then call ( x 1 , ⋯ , x i ) multiplicative variables and ( x i + 1 , ⋯ , x n ) non-multiplicative variables.

• Remaining equations of order ≤ q − 1 : Call ( x 1 , ⋯ , x n ) non-multiplicative variables.

In actual practice, we shall use a Janet tabular where the multiplicative “variables” are represented by their index in upper left position while the non-multiplicative variables are represented by dots in lower right position ( [

DEFINITION 2.2: A system of PD equations is said to be involutive if its first prolongation can be achieved by prolonging its equations only with respect to the corresponding multiplicative variables. In that case, we may introduce the

Cartan characters α q i = m ( q + n − i − 1 ) ! ( q − 1 ) ! ( n − i ) ! − β q i for i = 1 , ⋯ , n and we have

d i m ( g q ) = ∑ α q = α q 1 + ⋯ + α q n and d i m ( g q + 1 ) = ∑ i α q i = 1 α q 1 + ⋯ + n α q n . Moreover, one can exhibit the Hilbert polynomial d i m ( R q + r ) in r with leading term ( α / d ! ) r d with d ≤ n when α is the smallest non-zero character in the case of an involutive symbol. Such a prolongation allows to compute in a unique way the principal ( p r i ) jets from the parametric ( p a r ) other ones. This definition may also be applied to nonlinear systems as well.

REMARK 2.3: For an involutive system with β = β q n < m , then ( y β + 1 , ⋯ , y m ) can be given arbitrarily and may constitute the input variables in control theory, though it is not necessary to make such a choice. In this case, the intrinsic number α = α q n = m − β > 0 is called the n-character and is the system counterpart of the so-called “differential transcendence degree” in differential algebra and the “rank” in module theory. As we shall see in the next Section, the smallest non-zero character and the number of zero characters are intrinsic numbers that can most easily be known by bringing the system to involution and we have α q 1 ≥ ⋯ ≥ α q n ≥ 0 .

In the situation of the last remark, the following procedure will generalize for PD control systems the well known first order Kalman form of OD control systems where the derivatives of the input do not appear ( [

• Look at the equations of class n solved with respect to y n 1 , ⋯ , y n β .

• Use integrations by parts like:

y n 1 − a ( x ) y n β + 1 = d n ( y 1 − a ( x ) y β + 1 ) + ∂ n a ( x ) y β + 1 = y ¯ n 1 + ∂ n a ( x ) y β + 1

• Modify y 1 , ⋯ , y β to y ¯ 1 , ⋯ , y ¯ β in order to “absorb” the various y n β + 1 , ⋯ , y n m only appearing in the equations of class n.

We have the following unexpected result providing what we shall call a reduced Spencer form:

THEOREM 2.4: The new equations of class n contain y 1 , ⋯ , y β and their jets but only contain y i β + 1 , ⋯ , y i m with 0 ≤ i ≤ n − 1 while the equations of class 1, ⋯ , n − 1 no longer contain y β + 1 , ⋯ , y m and their jets. Accordingly, as we shall see in the next section, any torsion element, if it exists, only depends on y ¯ 1 , ⋯ , y ¯ β .

If χ 1 , ⋯ , χ n are n algebraic indeterminates or, in a more intrinsic way, if χ = χ i d x i ∈ T * is a covector and D : E → F : ξ → a k τ μ ( x ) ∂ μ ξ k ( x ) is a linear involutive operator of order q, we may introduce the characteristic matrix a ( x , χ ) = ( a k τ μ ( x ) χ μ , | μ | = μ 1 + ⋯ + μ n = q ) and the resulting map σ χ ( D ) : E → F is called the symbol of D at χ . Then there are two possibilities:

• If m a x χ r k ( σ χ ( D ) ) < m ⇔ α q n > 0 : the characteristic matrix fails to be injective for any covector.

• If m a x χ r k ( σ χ ( D ) ) = m ⇔ α q n = 0 : the characteristic matrix fails to be injective if and only if all the determinants of the m × m submatrices vanish. However, one can prove that this algebraic ideal a ∈ K [ χ ] is not intrinsically defined and must be replaced by its radical r a d ( a ) made by all polynomials having a power in a . This radical ideal is called the characteristic ideal of the operator.

DEFINITION 2.5: For each x ∈ X , the algebraic set defined by the characteristic ideal is called the characteristic set of D at x and V = ∪ x ∈ X V x is called the characteristic set of D while we keep the word “variety” for an irreducible algebraic set defined by a prime ideal.

One has the following important theorem ( [

THEOREM 2.6: (Hilbert-Serre) The dimension d ( V ) of the characteristic set, that is the maximum dimension of the irreducible components, is equal to the number of non-zero characters while the codimension c d ( V ) = n − d ( V ) is equal to the number of zero characters, that is to the number of “full” classes in the Janet tabular of an involutive system.

If P = a μ d μ ∈ D = K [ d ] with implicit summation on the multi-index, the highest value of | μ | with a μ ≠ 0 is called the order of the operator P and the ring D with multiplication ( P , Q ) → P ∘ Q = P Q is filtred by the order q of the operators. We have the filtration 0 ⊂ K = D 0 ⊂ D 1 ⊂ ⋯ ⊂ D q ⊂ ⋯ ⊂ D ∞ = D . Moreover, it is clear that D, as an algebra, is generated by K = D 0 and T = D 1 / D 0 with D 1 = K ⊕ T if we identify an element ξ = ξ i d i ∈ T with the vector field ξ = ξ i ( x ) ∂ i of differential geometry, but with ξ i ∈ K now. It follows that D = D D D is a bimodule over itself, being at the same time a left D-module by the composition P → Q P and a right D-module by the composition P → P Q . We define the adjoint map a d : D → D o p : P = a μ d μ → a d ( P ) = ( − 1 ) | μ | d μ a μ and we have a d ( a d ( P ) ) = P . It is easy to check that a d ( P Q ) = a d ( Q ) a d ( P ) , ∀ P , Q ∈ D . Such a definition can also be extended to any matrix of operators by using the transposed matrix of adjoint operators (see [

Accordingly, if y = ( y 1 , ⋯ , y m ) are differential indeterminates, then D acts on y k by setting d μ y k = y μ k with d i y μ k = y μ + 1 i k and y 0 k = y k . We may therefore use the jet coordinates in a formal way as in the previous section. Therefore, if a system of OD/PD equations is written in the form:

Φ τ ≡ a k τ μ y μ k = 0

with coefficients a k τ μ ∈ K , we may introduce the free differential module D y = D y 1 + ⋯ + D y m ≃ D m and consider the differential submodule I = D Φ ⊂ D y which is usually called the module of equations, both with the differential module M = D y / D Φ or D-module and we may set M = M D if we want to specify the ring of differential operators. The work of Macaulay only covers the case m = 1 with K replaced by k ⊆ c s t ( K ) . Again, we may introduce the formal prolongation with respect to d i by setting:

d i Φ τ ≡ a k τ μ y μ + 1 i k + ( ∂ i a k τ μ ) y μ k

in order to induce maps d i : M → M : y ¯ μ k → y ¯ μ + 1 i k if we use to denote the residue D y → M : y k → y ¯ k by a bar as in algebraic geometry. However, for simplicity, we shall not write down the bar when the background will indicate clearly if we are in Dy or in M.

As a byproduct, the differential modules we shall consider will always be finitely generated ( k = 1 , ⋯ , m < ∞ ) and finitely presented ( τ = 1 , ⋯ , p < ∞ ). Equivalently, introducing the matrix of operators D = ( a k τ μ d μ ) with m columns and p rows, we may introduce the morphism D p → D D m : ( P τ ) → ( P τ Φ τ ) : P → P Φ = P D over D by acting with D on the left of these row vectors while acting with D on the right of these row vectors and the presentation of M is defined by the exact cokernel sequence D p → D m → M → 0 . It is essential to notice that the presentation only depends on K , D and Φ or D , that is to say never refers to the concept of (explicit or formal) solutions. It is at this moment that we have to take into account the results of the previous section in order to understand that certain presentations will be much better than others, in particular to establish a link with formal integrability and involution.

DEFINITION 2.7: It follows from its definition that M can be endowed with a quotient filtration obtained from that of D m which is defined by the order of the jet coordinates y q in D q y . We have therefore the inductive limit 0 ⊆ M 0 ⊆ M 1 ⊆ ⋯ ⊆ M q ⊆ ⋯ ⊆ M ∞ = M with d i M q ⊆ M q + 1 and M = D M q for q ≫ 0 with prolongations D r M q ⊆ M q + r , ∀ q , r ≥ 0 . We shall set g r ( M q ) = G q = M q / M q − 1 and g r ( M ) = G = ⊕ q G q .

Having in mind that K is a left D-module for the action ( D , K ) → K : ( d i , a ) → ∂ i a and that D is a bimodule over itself, we have only two possible constructions:

DEFINITION 2.8: We define the system R = h o m K ( M , K ) = M * and set R q = h o m K ( M q , K ) = M q * as the system of order q. We have the projective limit R = R ∞ → ⋯ → R q → ⋯ → R 1 → R 0 . It follows that f q ∈ R q : y μ k → f μ k ∈ K with a k τ μ f μ k = 0 defines a section at order q and we may set f ∞ = f ∈ R for a section of R. For a ground field of constants k, this definition has of course to do with the concept of a formal power series solution. However, for an arbitrary differential field K, the main novelty of this new approach is that such a definition has nothing to do with the concept of a formal power series solution (care) as illustrated in ( [

DEFINITION 2.9: We may define the right differential module e x t 0 ( M ) = h o m D ( M , D ) .

PROPOSITION 2.10: When M is a left D-module, then R is also a left D-module.

Proof: As D is generated by K and T as we already said, let us define:

( a f ) ( m ) = a f ( m ) , ∀ a ∈ K , ∀ m ∈ M

( ξ f ) ( m ) = ξ f ( m ) − f ( ξ m ) , ∀ ξ = a i d i ∈ T , ∀ m ∈ M

In the operator sense, it is easy to check that d i a = a d i + ∂ i a and that ξ η − η ξ = [ ξ , η ] is the standard bracket of vector fields. We finally get ( d i f ) μ k = ( d i f ) ( y μ k ) = ∂ i f μ k − f μ + 1 i k and thus recover exactly the Spencer operator though this is not evident at all. We also get ( d i d j f ) μ k = ∂ i j f μ k − ∂ i f μ + 1 j k − ∂ j f μ + 1 i k + f μ + 1 i + 1 j k ⇒ d i d j = d j d i , ∀ i , j = 1, ⋯ , n and thus d i R q + 1 ⊆ R q ⇒ d i R ⊂ R induces a well defined operator R → T * ⊗ R : f → d x i ⊗ d i f . This result has been discovered (up to sign) by Macaulay in 1916 ( [

Q.E.D.

DEFINITION 2.11: t r ( M ) is the greatest differential submodule of M having codimension > r.

PROPOSITION 2.12: c d ( M ) = c d ( V ) = r ⇔ α q n − r ≠ 0 , α q n − r + 1 = ⋯ = α q n = 0 ⇔ t r ( M ) ≠ M , t r − 1 ( M ) = ⋯ = t 0 ( M ) = t ( M ) = M and this intrinsic result can be most easily checked by using the standard or reduced Spencer form of the system defining M.

We are now in a good position for defining and studying purity for differential modules.

DEFINITION 2.13: M is r-pure ⇔ t r ( M ) = 0 , t r − 1 ( M ) = M ⇔ c d ( D m ) = r , ∀ m ∈ M . More generally, M is pure if it is r-pure for a certain 0 ≤ r ≤ n . In particular, M is 0-pure if t ( M ) = 0 and, if c d ( M ) = r but M is not r-pure, we may call M / t r ( M ) the pure part of M. It follows that t r − 1 ( M ) / t r ( M ) is equal to zero or is r-pure (see the picture in [

0 = t n ( M ) ⊆ t n − 1 ( M ) ⊆ ⋯ ⊆ t 1 ( M ) ⊆ t 0 ( M ) = t ( M ) ⊆ M

PROPOSITION 2.14: t r ( M ) does not depend on the presentation or on the filtration of M.

EXAMPLE 2.15: If K = ℚ and M is defined by the involutive system y 33 = 0 , y 23 = 0 , y 13 = 0 , then z = y 3 satisfies d 3 z = 0 , d 2 z = 0 , d 1 z = 0 and c d ( D z ) = 3 while z ′ = y 2 only satisfies d 3 z ′ = 0 and c d ( D z ′ ) = 1 . We have the purity filtration 0 = t 3 ( M ) ⊂ t 2 ( M ) = t 1 ( M ) ⊂ t 0 ( M ) = t ( M ) = M with one gap and two strict inclusions.

We now recall the definition of the extension modules e x t D i ( M , D ) that we shall simply denote by e x t i ( M ) and the way to use their dimension or codimension. We point out once more that these numbers can be most easily obtained by bringing the underlying systems to involution in order to get informations on M from informations on G. We divide the procedure into four steps that can be achieved by means of computer algebra ( [

• STEP 1: Construct a free resolution of M, say:

⋯ → F i → ⋯ → F 1 → F 0 → M → 0

• STEP 2: Suppress M in order to obtain the deleted sequence:

⋯ → F i → ⋯ → F 1 → F 0 → 0

• STEP 3: Apply h o m D ( • , D ) in order to obtain the dual sequence heading backwards:

⋯ ← h o m D ( F i , D ) ← ⋯ ← h o m D ( F 1 , D ) ← h o m D ( F 0 , D ) ← 0

• STEP 4: Define e x t i ( M ) to be the cohomology at h o m D ( F i , D ) of the dual sequence in such a way that e x t 0 ( M ) = h o m D ( M , D ) .

The following nested chain of difficult propositions and theorems can be obtained, even in the non-commutative case, by combining the use of extension modules and bidualizing complexes in the framework of algebraic analysis. The main difficulty is to obtain first these results for the graded module G = g r ( M ) by using techniques from commutative algebra before extending them to the filtered module M as in ( [

THEOREM 2.16: The extension modules do not depend on the resolution of M used.

PROPOSITION 2.17: Applying h o m D ( • , D ) provides right D-modules that can be transformed to left D-modules by means of the side changing functor and vice-versa. Namely, if N D is a right D-module, then D N = ∧ n T ⊗ K N D is the converted left D-module while, if N = N D is a left D-module, then N D = ∧ n T * ⊗ K N is the converted right D-module.

PROPOSITION 2.18: Instead of using h o m D ( • , D ) and the side changing functor in the module framework, we may use a d in the operator framework. Namely, to any operator D : E → F we may associate the formal adjoint a d ( D ) : ∧ n T * ⊗ F * → ∧ n T * ⊗ E * with the useful though striking relation r k D ( a d ( D ) ) = r k D ( D ) .

PROPOSITION 2.19: e x t i ( M ) is a torsion module ∀ 1 ≤ i ≤ n but e x t 0 ( M ) = h o m D ( M , D ) may not be a torsion module.

EXAMPLE 2.20: When M is a torsion module, we have h o m D ( M , D ) = 0 . Indeed, if m ∈ M , we may find 0 ≠ P ∈ D such that P m = 0 . Hence, if f ∈ h o m D ( M , D ) , we have P f ( m ) = f ( P m ) = f ( 0 ) = 0 in D and thus f ( m ) = 0 , ∀ m ∈ M , that is to say f = 0 because D is an integral domain. When n = 3 and the torsion-free module M is defined by the formally surjective d i v operator, the formal adjoint of d i v is − g r a d which defines a torsion module. Also, when n = 1 as in classical control theory, a controllable system with coefficients in a differential field allows to define a torsion-free module M which is free in that case because a finitely generated module over a principal ideal domain is free if and only if it is torsion-free and h o m D ( M , D ) is thus also a free module.

THEOREM 2.21: e x t i ( M ) = 0 , ∀ i < c d ( M ) and ∀ i ≥ n + 1 .

THEOREM 2.22: c d ( e x t i ( M ) ) ≥ i .

THEOREM 2.23: c d ( M ) ≥ r ⇔ e x t i ( M ) = 0 , ∀ i < r .

PROPOSITION 2.24: c d ( M ) = r ⇒ c d ( e x t r ( M ) ) = r and e x t r ( M ) is r-pure.

PROPOSITION 2.25: e x t r ( e x t r ( M ) ) is equal to 0 or is r-pure, ∀ 0 ≤ r ≤ n .

PROPOSITION 2.26: If we set t − 1 ( M ) = M , there are exact sequences ∀ 0 ≤ r ≤ n :

0 → t r ( M ) → t r − 1 ( M ) → e x t r ( e x t r (M) )

THEOREM 2.27: If c d ( M ) = r , then M is r-pure if and only if there is a monomorphism 0 → M → e x t r ( e x t r ( M ) ) of left differential modules.

THEOREM 2.28: M is pure ⇔ e x t s ( e x t s ( M ) ) = 0, ∀ s ≠ c d ( M ) .

COROLLARY 2.29: If M is r-pure with r ≥ 1 , then it can be embedded into a differential module L having a free resolution with only r operators.

The previous theorems are known to characterize purity but it is however evident that they are not very useful in actual practice. For more details on these two results which are absolutely out of the scope of this paper, see ( [

THEOREM 2.30: When M is r-pure, the characteristic ideal is thus unmixed, that is a finite intersection of prime ideals having the same codimension r and the characteristic set is equidimensional, that is the union of irreducible algebraic varieties having the same codimension r.

In 2012 we have provided a new effective test for checking purity while using the involutivity of the Spencer form with four steps as follows ( [

• STEP 1: Compute the involutive Spencer form of the system and the number r of full classes.

• STEP 2: Select only the equations of class 1 to d ( M ) = n − r of this Spencer form which are making an involutive system over K [ d 1 , ⋯ , d ( n − r ) ] .

• STEP 3: Using differential biduality for such a system, check if it defines a torsion-free module M ( n − r ) and work out a parametrization.

• STEP 4: Substitute the above parametrization in the remaning equations of class n − r + 1, ⋯ , n of the Spencer form in order to get a system of PD equations which provides the parametrizing module L in such a way that M ⊆ L and L has a resolution with r operators.

THEOREM 2.31: As purity is an intrinsic property, we may work with an involutive Spencer form and M is r-pure if the classes n − r + 1, ⋯ , n are full and the module M ( n − r ) defined by the equations of class 1 + ⋯ + class ( n − r ) is torsion-free. Hence M is 0-pure if it is torsion-free.

We shall now illustrate and apply this new procedure in the next two sections.

EXAMPLE 3.1: With n = 3 , m = 1 and K = ℚ , let us consider the following polynomial ideal:

a = ( ( χ 3 ) 2 , χ 2 χ 3 − χ 1 χ 3 , ( χ 2 ) 2 − χ 1 χ 2 ) ⊂ K [ χ 1 , χ 2 , χ 3 ] = K [ χ ]

We shall discover that it is not evident to prove that it is an unmixed polynomial ideal and that the corresponding differential module is 1-pure.

The first result is provided by the existence of the primary decomposition obtained from the two existing factorizations ( [

a = ( ( χ 3 ) 2 , χ 2 − χ 1 ) ∩ ( χ 3 , χ 2 ) = q ′ ∩ q ″

Taking the respective radical ideals, we get the prime decomposition:

r a d ( a ) = ( χ 3 , χ 2 − χ 1 ) ∩ ( χ 3 , χ 2 ) = p ′ ∩ p ″ = r a d ( q ′ ) ∩ r a d (q″)

The corresponding involutive system is:

{ y 33 = 0 y 23 − y 13 = 0 y 22 − y 12 = 0 1 2 3 1 2 • 1 2 •

with characters α 2 3 = 1 − 1 = 0 , α 2 2 = 2 − 2 = 0 , α 2 1 = 3 − 0 = 3 and d i m ( g 2 ) = ∑ α = 3 .

Setting ( z 1 = y , z 2 = y 1 , z 3 = y 2 , z 4 = y 3 ), we obtain the involutive first order Spencer form:

{ z 3 4 = 0 , z 3 3 − z 1 4 = 0 , z 3 2 − z 1 4 = 0 , z 3 1 − z 4 = 0 z 2 4 − z 1 4 = 0 , z 2 3 − z 1 3 = 0 , z 2 2 − z 1 3 = 0 , z 2 1 − z 3 = 0 z 1 1 − z 2 = 0 1 2 3 1 2 • 1 • •

with new characters α 1 3 = 4 − 4 = 0 , α 1 3 = 4 − 4 = 0 , α 1 1 = 4 − 1 = 3 and similarly d i m ( g 1 ) = ∑ α = 3 . Both class 3 and class 2 are full while class 1 is defining a torsion-free module M ( 1 ) over K [ d 1 ] by means of a trivially involutive system of class 1. Hence the differential module M is such that c d ( M ) = 2 and is 1-pure because it is 1-pure in this presentation.

Suppressing the bar for the various residues, we are ready to exhibit the relative parametrization defining the parametrization module L because we may choose the 3 potentials ( z 1 = y , z 3 , z 4 ) while taking into account that z 2 = y 1 = d 1 y :

{ z 3 4 = 0 z 3 3 − z 1 4 = 0 z 3 1 − z 4 = 0 z 2 4 − z 1 4 = 0 z 2 3 − z 1 3 = 0 z 2 1 − z 3 = 0 1 2 3 1 2 3 1 2 3 1 2 • 1 2 • 1 2 •

Both ( y , z 3 , z 4 ) are torsion elements and we can eliminate ( z 3 , z 4 ) in order to find the desired system that must be satisfied by y which is showing the inclusion M ⊂ L but we have indeed M = L because z 3 = y 2 , z 4 = y 3 . It follows that M admits a free resolution with only 2 operators, a result following at once from the last Janet tabular, contrary to the previous one.

The reader may treat similarly the example a = ( χ 1 , χ 2 ) ∩ ( χ 3 , χ 4 ) and look at ( [

EXAMPLE 3.2: With n = 3 , m = 1 , q = 2 , K = ℚ , D = K [ d 1 , d 2 , d 3 ] , let us consider the differential module M defined by the second order system P y ≡ y 33 = 0 , Q y ≡ y 13 − y 2 = 0 first considered by Macaulay in 1916 ( [

First of all, in order to find out the codimension c d ( M ) = 2 , we have to consider the equivalent involutive system:

{ Φ 4 ≡ y 33 = u Φ 3 ≡ y 23 = d 1 u − d 3 v Φ 2 ≡ y 22 = d 11 u − d 13 v − d 2 v Φ 1 ≡ y 13 − y 2 = v 1 2 3 1 2 • 1 2 • 1 • •

The Janet tabular on the right allows at once to compute the characters α 2 3 = 0 , α 2 2 = 0 , α 2 1 = α = 3 − 1 = 2 and to construct the following strictly exact sequence of differential modules:

0 → D → D 4 → D 4 → D → p M → 0

Also, we have

r a d ( a ) = r a d ( ( χ 3 ) 2 , χ 2 χ 3 , ( χ 2 ) 2 , χ 1 χ 3 ) = ( χ 3 , χ 2 ) = p ⇒ d i m ( V ) = 1 .

As the classes 3 and 2 are full, it follows that d ( M ) = d ( D y ) = 1 ⇒ c d ( M ) = n − 1 = 2 if we denote simply by y the canonical residue y ¯ of y after identifying D with D y . We have constructed explicitly in ( [

0 → D → 2 D 1 D 2 → D D → p M → 0

Indeed, introducing differential duality through the functor h o m D ( • , D ) and the respective adjoint operators, we may define the torsion left differential module N by the long exact sequence:

0 ← N ← q D ← a d ( D 1 ) D 2 → a d ( D ) D ← 0

showing that r k D ( M ) = 0 ⇒ r k D ( N ) = 1 − 2 + 1 = 0 because of the additivity property of the differential rank and the vanishing of the Euler-Poincaré characteristic of the full sequence. It follows that M = e x t 2 ( N ) = e x t 2 ( e x t 2 ( M ) ) .

Similarly, using certain parametric jet variables as new unknowns, we may set z 1 = y , z 2 = y 1 , z 3 = y 2 , z 4 = y 3 in order to obtain the following involutive first order system with no zero order equation:

{ class 3 d 3 z 1 − z 4 = 0 , d 3 z 2 − z 3 = 0 , d 3 z 3 = 0 , d 3 z 4 = 0 class 2 d 2 z 1 − z 3 = 0 , d 2 z 2 − d 1 z 3 = 0 , d 2 z 3 = 0 , d 2 z 4 = 0 class 1 d 1 z 1 − z 2 = 0 , d 1 z 4 − z 3 = 0 1 2 3 1 2 • 1 • •

where we have separated the classes while using standard computer algebra notations this time instead of the jet notations used in the previous example. Contrary to what could be believed, this operator does not describe the Spencer sequence that could be obtained from the previous Janet sequence but we can use it exactly like a Janet sequence or exactly like a Spencer sequence. We obtain therefore a long strictly exact sequence of differential modules with only first order operators while replacing D y by D z = D z 1 + D z 2 + D z 3 + D z 4 as follows:

0 → D 2 → 1 D 8 → 1 D 10 → 1 D 4 → M → 0

and we still have the vanishing Euler-Poincaré characteristic 2 − 8 + 10 − 4 = 0 .

The differential module M 1 is defined over K [ d 1 ] by the two PD equations of class 1 and is easily seen to be torsion-free with the two potentials ( z 1 = y , z 4 ) . Substituting into the PD equations of class 2 and 3, we obtain the generating differential constraints:

{ d 3 z 1 − z 4 = 0 d 3 z 4 = 0 d 2 z 1 − d 1 z 4 = 0 d 2 z 4 = 0 1 2 3 1 2 3 1 2 • 1 2 •

They define the parametrization module L and the inclusion M ⊆ L is obtained by eliminating z 4 but we have indeed M = L because z 4 = d 3 y .

EXAMPLE 3.3: We have provided in ( [

{ y 44 = 0 y 34 = 0 y 33 = 0 y 24 − y 13 = 0 1 2 3 4 1 2 3 • 1 2 3 • 1 2 • •

From the Janet tabular we may construct at once the Janet sequence:

0 → Θ → 1 → D 4 → D 1 4 → D 2 1 → 0

where D 1 is defined by the involutive system:

{ d 4 ( y 34 ) − d 3 ( y 44 ) = 0 d 4 ( y 33 ) − d 3 ( y 34 ) = 0 d 4 ( y 24 − y 13 ) − d 2 ( y 44 ) + d 1 ( y 34 ) = 0 d 3 ( y 24 − y 13 ) − d 2 ( y 34 ) + d 1 ( y 33 ) = 0 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 •

and so on. We have therefore a free resolution of M with 3 operators:

0 → D → D 4 → D 4 → D → M → 0

and thus discover that p d ( M ) ≤ 3 .

However, we have

r a d ( ( χ 4 ) 2 , χ 3 χ 4 , ( χ 3 ) 2 , χ 2 χ 4 − χ 1 χ 3 ) = ( χ 4 , χ 3 ) = p ⇒ d i m ( V ) = 1 .

Let us transform the initial second order involutive system for y into a first order involutive system for ( z 1 = y , z 2 = y 1 , z 3 = y 2 , z 4 = y 3 , z 5 = y 4 ) as follows:

{ d 4 z 1 − z 5 = 0 , d 4 z 2 − d 1 z 5 = 0 , d 4 z 3 − d 1 z 4 = 0 , d 4 z 4 = 0 , d 4 z 5 = 0 d 3 z 1 − z 4 = 0 , d 3 z 2 − d 1 z 4 = 0 , d 3 z 3 − d 2 z 4 = 0 , d 3 z 4 = 0 , d 3 z 5 = 0 d 2 z 1 − z 3 = 0 , d 2 z 2 − d 1 z 3 = 0 , d 2 z 5 − d 1 z 4 = 0 d 1 z 1 − z 2 = 0 1 2 3 4 1 2 3 • 1 2 • • 1 • • •

with five equations of full class 4, five equations of full class 3, three equations of class 2 and finally one equation of class 1. The equations of classes 2 and 1 are providing an involutive system over ℚ [ d 1 , d 2 ] defining a torsion-free module M ( 2 ) that can be parametrized by setting z 1 = y , z 2 = d 1 y , z 3 = d 2 y , z 4 = d 2 z , z 5 = d 1 z with only 2 arbitrary potentials ( y , z ) . Substituting in the other equations of classes 3 and 4, we finally discover that L is defined by the involutive system describing the relative parametrization:

{ d 4 y − d 1 z = 0 d 4 z = 0 d 3 y − d 2 z = 0 d 3 z = 0 1 2 3 4 1 2 3 4 1 2 3 • 1 2 3 •

We have the strict inclusion M ⊂ L obtained by eliminating z because now z ∉ D y if we take the residue or, equivalently, the residue of z does not belong to M. The differential module L defined by the above system is therefore 2-pure with a strict inclusion M ⊂ L and admits a free resolution with only 2 operators according to its Janet tabular.

EXAMPLE 3.4: (Contact structure) With n = m = 3 and K = ℚ ( x 1 , x 2 , x 3 ) let us introduce the so-called contact 1-form α = d x 1 − x 3 d x 2 and consider the first order system of infinitesimal Lie equations obtained by eliminating the contact factor ρ from the equations L ( ξ ) α = ρ α . We let the reader check that he will obtain only the two equations Φ 1 = 0 , Φ 2 = 0 which is nevertheless neither formally integrable nor even involutive. Using crossed derivatives one obtains the involutive system:

{ Φ 3 ≡ ∂ 3 ξ 3 + ∂ 2 ξ 2 + 2 x 3 ∂ 1 ξ 2 = 0 Φ 2 ≡ ∂ 3 ξ 1 − x 3 ∂ 3 ξ 2 = 0 Φ 1 ≡ ∂ 2 ξ 1 − x 3 ∂ 2 ξ 2 + x 3 ∂ 1 ξ 1 − ( x 3 ) 2 ∂ 1 ξ 2 − ξ 3 = 0 1 2 3 1 2 3 1 2 •

with the unique CC Ψ ≡ ∂ 3 Φ 1 − ∂ 2 Φ 2 − x 3 ∂ 1 Φ 2 + Φ 3 = 0 . The following injective absolute parametrization is well known and we let the reader find it by using differential double duality:

ϕ − x 3 ∂ 3 ϕ = ξ 1 , − ∂ 3 ϕ = ξ 2 , ∂ 2 ϕ + x 3 ∂ 1 ϕ = ξ 3 ⇒ ξ 1 − x 3 ξ 2 = ϕ

We obtain the Janet sequence

0 → 1 → D − 1 3 → D 3 → D 1 1 → 0 ϕ ξ Φ Ψ

with formally exact adjoint sequence:

0 ← 1 ← a d ( D − 1 ) 3 ← a d ( D ) 3 ← a d ( D 1 ) 1 → 0 θ ν μ λ

and the resolution of the trivially torsion-free module M ≃ D :

0 → D → D 3 → D 3 → M → 0

which splits totally because it is made with free and thus projective modules.

EXAMPLE 3.5: (Unimodular contact structure) With n = m = 3 and K = ℚ ( x 1 , x 2 , x 3 ) let us introduce the 1-form ω = d x 1 − x 3 d x 2 used as a geometric object and consider the first order system of infinitesimal Lie equations from the equations L ( ξ ) ω = 0 . One obtains the system using jet notations:

ξ 3 1 − x 3 ξ 3 2 = 0 , ξ 2 1 − x 3 ξ 2 2 − ξ 3 = 0 , ξ 1 1 − x 3 ξ 1 2 = 0

We let the reader prove that these three PD equations are differentially independent and we obtain the free resolution of M:

0 → D 3 → D D 3 → M → 0

and its adjoint sequence is:

0 ← N ← D 3 ← a d ( D ) D 3 ← 0

because r k D ( M ) = r k D ( N ) = 3 − 3 = 0 , that is both M and N are torsion modules with N = e x t 1 ( M ) ⇒ M = e x t 1 ( N ) = e x t 1 ( e x t 1 ( M ) ) and M is surely 1-pure. However, this system is not formally integrable, as it can be checked directly through crossed derivatives or by noticing that L ( ξ ) d ω = 0 with d ω = d x 2 ∧ d x 3 and L ( ξ ) ( ω ∧ d ω ) = 0 with ω ∧ d ω = d x 1 ∧ d x 2 ∧ d x 3 . Hence, we have to add the 3 first order equations:

ξ 2 2 + ξ 3 3 = 0 , ξ 1 3 = 0 , ξ 1 2 = 0 ⇒ ξ 1 1 = 0

Exchanging x 1 and x 3 , we obtain the equivalent involutive system in δ -regular coordinates:

{ ξ 3 3 = 0 ξ 3 2 = 0 ξ 3 1 = 0 ξ 2 2 + ξ 1 3 = 0 ξ 2 1 + x 1 ξ 1 3 − ξ 3 = 0 ξ 1 1 − x 1 ξ 1 2 = 0 1 2 3 1 2 3 1 2 3 1 2 • 1 2 • 1 • •

The differential module M ( 2 ) over K [ d 1 , d 2 ] is defined by the three bottom equations. Setting now ϕ = ξ 1 − x 1 ξ 2 , we deduce from the last bottom equation that ξ 2 = − d 1 ϕ and thus ξ 1 = ϕ − x 1 d 1 ϕ . Finally, substituting in the equation before the last, we get ξ 3 = d 2 ϕ . We have thus obtained an injective parametrization of M ( 2 ) which is therefore torsion-free and M is 2-pure in a coherent way. Substituting into the three upper equations, we obtain the desired relative parametrization by adding the differential constraint d 3 ϕ = 0 . Coming back to the original coordinates, we obtain the relative parametrization:

ϕ − x 3 d 3 ϕ = ξ 1 , − d 3 ϕ = ξ 2 , d 2 ϕ = ξ 3 with d 1 ϕ = 0 ⇒ d 2 ϕ + x 3 d 1 ϕ = ξ 3

which is thus strikingly obtained from the previous contact parametrization by adding the only differential constraint d 1 ϕ = 0 obtained by substituting it in the new system of Lie equations.

Before studying applications to mathematical physics, we shall start with an example describing in an explicit way the Janet and Spencer sequences used thereafter, both with their link, namely the relations existing between the dimensions of the respective Janet and Spencer bundles.

EXAMPLE 4.1: When n = m = 2 , q = 2 , ω is the Euclidean metric of X = ℝ 2 with Christoffel symbols γ and metric density ω ˜ = ω / d e t ( ω ) , we consider the two involutive systems of linear infinitesimal Lie equations R 2 ⊂ R ˜ 2 ⊂ J 2 ( T ) respectively defined by { L ( ξ ) ω = 0, L ( ξ ) γ = 0 } and { L ( ξ ) ω ˜ = 0, L ( ξ ) γ = 0 } . We have g 2 = g ˜ 2 = 0 and construct the following successive commutative and exact diagrams followed by the corresponding dimensional diagrams that are used in order to construct effectively the respective Janet and Spencer differential sequences while comparing them.

0 0 0 ↘ ↓ ↓ 0 → S 2 T * ⊗ T → F ″ 0 → 0 ↓ ↓ ↘ ↓ 0 → R 2 → J 2 ( T ) → F 0 → 0 ↓ ↓ ↓ 0 → R 1 → J 1 ( T ) → F ′ 0 → 0 ↓ ↓ ↓ 0 0 0

0 0 0 ↘ ↓ ↓ 0 → 6 → 6 → 0 ↓ ↓ ↘ ↓ 0 → 3 → 12 → 9 → 0 ↓ ↓ ↓ 0 → 3 → 6 → 3 → 0 ↓ ↓ ↓ 0 0 0

In the present situation we notice that R 2 = ρ 1 ( R 1 ) = J 1 ( R 1 ) ∩ J 2 ( T ) ⊂ J 1 ( J 1 ( T ) ) and thus F 0 ≃ J 1 ( F ′ 0 ) with F ″ 0 ≃ T * ⊗ F ′ 0 ≃ S 2 T * ⊗ T by counting the dimensions because we have surely F 0 ⊂ J 1 ( F ′ 0 ) with g 2 = 0 .

0 0 0 ↓ ↓ ↓ 0 → S 3 T * ⊗ T → T * ⊗ F 0 → F 1 → 0 ↓ ↓ ↓ ∥ 0 → R 3 → J 3 ( T ) → J 1 ( F 0 ) → F 1 → 0 ↓ ↓ ↓ ↓ 0 → R 2 → J 2 ( T ) → F 0 → 0 ↓ ↓ ↓ 0 0 0

0 0 0 ↓ ↓ ↓ 0 → 8 → 18 → 10 → 0 ↓ ↓ ↓ ∥ 0 → 3 → 20 → 27 → 10 → 0 ∥ ↓ ↓ ↓ 0 → 3 → 12 → 9 → 0 ↓ ↓ ↓ 0 0 0

SPENCER

0 0 0 ↓ ↓ ↓ 0 → Θ → j q C 0 → D 1 C 1 → D 2 C 2 → 0 ↓ ↓ ↓ 0 → T → j q C 0 ( T ) → D 1 C 1 ( T ) → D 2 C 2 ( T ) → 0 ∥ ↓ Φ 0 ↓ Φ 1 ↓ Φ 2 0 → Θ → T → D F 0 → D 1 F 1 → D 2 F 2 → 0 ↓ ↓ ↓ 0 0 0

JANET

0 → Θ ˜ → j 2 4 → D 1 8 → D 2 4 → 0 0 → Θ → j 2 3 → D 1 6 → D 2 3 → 0 SPENCER ↓ ↓ ↓ 0 → 2 → j 2 12 → D 1 16 → D 2 6 → 0 HYBRID ∥ ↓ Φ 0 ↓ Φ 1 ↓ Φ 2 0 → Θ → 2 → D 9 → D 1 10 → D 2 3 → 0 JANET 0 → Θ ˜ → 2 → D 8 → D 1 8 → D 2 2 → 0

0 0 ↓ ↓ 0 → T * ⊗ ( C ˜ 0 / C 0 ) = C ˜ 1 / C 1 → 0 ↓ ↓ ↓ 0 → S 3 T * ⊗ T → T * ⊗ F 0 → F 1 → 0 ∥ ↓ ↓ 0 → S 3 T * ⊗ T → T * ⊗ F ˜ 0 → F ˜ 1 → 0 ↓ ↓ ↓ 0 0 0

0 0 ↓ ↓ 0 → 2 → 2 → 0 ↓ ↓ ↓ 0 → 8 → 18 → 10 → 0 ∥ ↓ ↓ 0 → 8 → 16 → 8 → 0 ↓ ↓ ↓ 0 0 0

C ˜ r / C r ≃ k e r ( F r → F ˜ r ) ≃ ∧ r T * ⊗ ( C ˜ 0 / C 0 ) ⇒ d i m ( C ˜ r ) − d i m ( C r ) = d i m ( F r ) − d i m ( F ˜ r ) ⇒ d i m ( C r ) + d i m ( F r ) = d i m ( C ˜ r ) + d i m ( F ˜ r )

In this new situation, we now notice that

R ˜ 2 ⊊ ρ 1 ( R ˜ 1 ) = J 1 ( R ˜ 1 ) ∩ J 2 ( T ) ⊂ J 1 ( J 1 (T) )

and the induced morphism F ˜ 0 → J 1 ( F ˜ ′ 0 ) is thus no longer a monomorphism though we still have an isomorphism F ˜ ″ 0 ≃ S 2 T * ⊗ T because g ˜ 2 = 0 again. Finally, we may extend such a procedure to the conformal group of space-time by considering the system of infinitesimal conformal transformations of the Minkowski metric defined by the first order system R ^ 1 ⊂ J 1 ( T ) in such a way that we have the strict inclusions R 2 ⊂ R ˜ 2 ⊂ R ^ 2 ⊂ J 2 ( T ) with d i m ( g ^ 2 ) = n = 4 .

For this, we just need to introduce the metric density ω ^ = ω ( | d e t ( ω ) | ) − 1 n and consider the system L ( ξ ) ω ^ = 0 ( [

Ω ^ i j ≡ ω ^ r j ∂ i ξ r + ω ^ i r ∂ j ξ r − 2 n ω ^ i j ∂ r ξ r + ξ r ∂ r ω ^ i j = 0

A) Control Theory:

EXAMPLE 4.A.1: (OD control theory) In classical control theory we have n = 1 and the only independent variable is the time, simply denoted by x but we may choose any ground differential field like K = ℚ ( x ) . In that case, we shall refer to ( [

{ y 1 → λ x x = μ 1 y 2 → λ x + a λ = μ 2 1 •

This system has a trivially involutive zero symbol but is not even formally integrable and we have to consider:

{ λ x x = μ 1 λ x + a λ = μ 2 ( ∂ x a − a 2 ) λ = μ x 2 − μ 1 − a μ 2 1 • •

We have thus two possibilities:

• We have a x − a 2 ≠ 0 and the adjoint system has the only zero solution, that is the adjoint operator is injective. In this case N = 0 and thus t ( M ) = e x t 1 ( N ) = 0 that is M is torsion-free. However, as n = 1 it follows that D = K [ d ] is a principal ideal ring which is therefore free and thus projective ( [

• The Riccati equation a x − a 2 = 0 is satisfied, for example if a = − 1 / x , and we get the CC μ x 2 − μ 1 − a μ 2 = 0 . Multiplying by a test function ξ and integrating by parts, we get the adjoint operator:

{ μ 1 → − ξ = η 1 μ 2 → − ξ x − a ξ = η 2

with only one first order generating CC, namely ∂ x η 1 − η 2 + a η 1 = 0 . It follows that N ≠ 0 ⇒ e x t 1 ( N ) ≠ 0 is a torsion module generated by the residue of z = y x 1 − y 2 + a y 1 , even though y 1 and y 2 are free separately and M is not pure. We obtain indeed a torsion element as we can check at once that z x − a z = 0 and wish good luck for control people to recover simply this result even on such an elementary example because the Kalman criterion is only working for systems with constant coefficients (compare [

EXAMPLE 4.A.2: (PD control theory) With n = 2 , let us consider the (trivially involutive) inhomogeneous single first order PD equations with two independent variables ( x 1 , x 2 ) , two unknown functions ( η 1 , η 2 ) and a second member ζ :

∂ 2 η 1 − ∂ 1 η 2 + x 2 η 2 = ζ ⇔ D 1 η = ζ

The ring of differential operators is D = K [ d 1 , d 2 ] with K = ℚ ( x 1 , x 2 ) . Multiplying on the left by a test function λ and integrating by parts, the corresponding adjoint operator is described by:

{ η 1 → − ∂ 2 λ = μ 1 η 2 → ∂ 1 λ + x 2 λ = μ 2 ⇔ a d ( D 1 ) λ = μ

Using crossed derivatives, this operator is injective because λ = ∂ 2 μ 2 + ∂ 1 μ 1 + x 2 μ 1 and we even obtain a lift λ → μ → λ . Substituting, we get the two CC:

{ ∂ 22 μ 2 + ∂ 12 μ 1 + x 2 ∂ 2 μ 1 + 2 μ 1 = ν 1 ∂ 12 μ 2 + ∂ 11 μ 1 + 2 x 2 ∂ 1 μ 1 + x 2 ∂ 2 μ 2 + ( x 2 ) 2 μ 1 − μ 2 = ν 2 1 2 1 •

This system is involutive and the corresponding generating CC for the second member ( ν 1 , ν 2 ) is:

∂ 2 ν 2 − ∂ 1 ν 1 − x 2 ν 1 = 0

Therefore ν 2 is differentially dependent on ν 1 but ν 1 is also differentially dependent on ν 2 . Multiplying on the left by a test function θ and integrating by parts, the corresponding adjoint system of PD equations is:

{ ν 1 → ∂ 1 θ − x 2 θ = ξ 1 ν 2 → − ∂ 2 θ = ξ 2 ⇔ a d ( D − 1 ) θ = ξ

Multiplying now the first equation by the test function ξ 1 , the second equation by the test function ξ 2 , adding and integrating by parts, we get the canonical parametrization D ξ = η :

{ μ 2 → ∂ 22 ξ 1 + ∂ 12 ξ 2 − x 2 ∂ 2 ξ 2 − 2 ξ 2 = η 2 μ 1 → ∂ 12 ξ 1 − x 2 ∂ 2 ξ 1 + ξ 1 + ∂ 11 ξ 2 − 2 x 2 ∂ 1 ξ 2 + ( x 2 ) 2 ξ 2 = η 1 1 2 1 •

of the initial system with zero second member. This system is involutive and the kernel of this parametrization has differential rank equal to 1 because ξ 1 or ξ 2 can be given arbitrarily.

Keeping now ξ 1 = ξ while setting ξ 2 = 0 , we get the first second order minimal parametrization ξ → ( η 1 , η 2 ) :

{ ∂ 22 ξ = η 2 ∂ 12 ξ − x 2 ∂ 2 ξ + ξ = η 1 1 2 1 •

This system is again involutive and the parametrization is minimal because the kernel of this parametrization has differential rank equal to 0. With a similar comment, setting now ξ 1 = 0 while keeping ξ 2 = ξ ′ , we get the second order minimal parametrization ξ ′ → ( η 1 , η 2 ) :

{ ∂ 11 ξ ′ − 2 x 2 ∂ 1 ξ ′ + ( x 2 ) 2 ξ ′ = η 1 ∂ 12 ξ ′ − x 2 ∂ 2 ξ ′ − 2 ξ ′ = η 2

which is again easily seen to be involutive by exchanging x 1 with x 2 .

With again a similar comment, setting now ξ 1 = ∂ 1 ϕ , ξ 2 = − ∂ 2 ϕ in the canonical parametrization, we obtain the third different second order minimal parametrization:

{ x 2 ∂ 22 ϕ + 2 ∂ 2 ϕ = η 2 x 2 ∂ 12 ϕ − ( x 2 ) 2 ∂ 2 ϕ + ∂ 1 ϕ = η 1 1 2 1 •

We are now ready for understanding the meaning and usefulness of what we have called “relative parametrization” in ( [

{ ∂ 2 ξ 1 + ∂ 1 ξ 2 = 0 − x 2 ∂ 2 ξ 2 − 2 ξ 2 = η 2 − x 2 ∂ 1 ξ 2 + ( x 2 ) 2 ξ 2 + ξ 1 = η 1 1 2 1 2 1 •

with 2 equations of class 2 (thus with class 2 full) and only 1 equation of class 1.

In a different way, we may add the differential constraint ∂ 1 ξ 1 + ∂ 2 ξ 2 = 0 but we have to check that it is compatible with the previous parametrization. For this, we have to consider the following second order system which is easily seen to be involutive with 2 second order equations of (full) class 2, (only) 2 second order equations of class 1 and 1 equation of order 1:

{ ∂ 22 ξ 2 + ∂ 12 ξ 1 = 0 ∂ 22 ξ 1 + ∂ 12 ξ 2 − x 2 ∂ 2 ξ 2 − 2 ξ 2 = η 2 ∂ 12 ξ 2 + ∂ 11 ξ 1 = 0 ∂ 12 ξ 1 − x 2 ∂ 2 ξ 1 + ξ 1 + ∂ 11 ξ 2 − 2 x 2 ∂ 1 ξ 2 + ( x 2 ) 2 ξ 2 = η 1 ∂ 2 ξ 2 + ∂ 1 ξ 1 = 0 1 2 1 2 1 • 1 • • •

The 4 generating CC only produce the desired system for ( η 1 , η 2 ) as we wished.

We cannot impose the condition D − 1 θ = ξ already found as it should give the identity 0 = η .

It is however also important to notice that the strictly exact long exact sequence:

0 → D → D 1 D 2 → D D 2 → D − 1 D → 0

splits because we have a lift ζ → η → ζ , namely:

ζ → ( − ∂ 1 ζ + x 2 ζ = η 1 , − ∂ 2 η 2 = η 2 ) → ∂ 2 η 1 − ∂ 1 η 2 + x 2 η 2 = ζ

We have thus an isomorphism D 2 ≃ D ⊕ M in the resolution 0 → D → D 1 D 2 → p M → 0 and all the differential modules defined from the operators involved are projective, thus torsion-free or 0-pure with vanishing r-extension modules e x t r ( ) = 0, ∀ r ≥ 1 .

As an exercise, we finally invite the reader to study the situation met with the system ∂ 2 η 1 − ∂ 1 η 2 + a ( x ) η 2 whenever a ∈ K (Hint: The controllability condition is now ∂ 1 a ≠ 0 ). The comparison with the previous OD case needs no comment.

B) Electromagnetism:

Most physicists know the Maxwell equations in vacuum, eventually in dielectrics and magnets, but are largely unaware of the more delicate constitutive laws involved in field-matter couplings like piezzoelectricity, photoelasticity or streaming birefringence. In particular they do not know that the phenomenological laws of these phenomena have been given ... by Maxwell ( [

Though it is quite useful in actual practice, the following approach to Maxwell equations cannot be found in any textbook. Namely, avoiding any variational calculus based on given Minkowski constitutive laws F ∼ F between field F and induction F for dielectric or magnets, let us use differential duality and define the first set M 1 of Maxwell equations by d : ∧ 2 T * → ∧ 3 T * while the second set M 2 will be defined by a d ( d ) : ∧ 4 T * ⊗ ∧ 2 T → ∧ 4 T * ⊗ T with d : T * → ∧ 2 T * , in a totally independent and intrinsic manner, using now contravariant tensor densities in place of covariant tensors. As we have already proved since a long time in ( [

potential → d field → induction → a d ( d ) current A → F → F → J

which is responsible for EM waves, though it is equivalent to the composition in the right square:

pseudopotential → a d ( d ) induction → field → d ∧ 3 T *

The main difference is that we need to set J = 0 in the first approach because of M 2 while we get automatically such a vanishing assumption in the second approach because of M 1 , avoiding therefore the Lorenz condition as in ( [

potential = ( A i ) → d ( ∂ i A j − ∂ j A i = F i j ) = field → M 1 ( ∂ i F j k + ∂ j F k i + ∂ k F i j = 0 ) current = ( ∂ i F i j = J j ) ← M 2 ( F i j ) = induction ← a d ( d ) pseudopotential

A F T * → d ∧ 2 T * → d = M 1 ∧ 3 T * ⋮ ↕ ⋮ ∧ 4 T * ⊗ T ← a d ( d ) = M 2 ∧ 4 T * ⊗ ∧ 2 T ← a d ( d ) ∧ 4 T * ⊗ ∧ 3 T ↕ ↕ ↕ ∧ 3 T * ← d ∧ 2 T * ← d T * J F

Using symbolic notations with an Euclidean metric instead of the Minkowski one because they are both locally constant while using the constitutive law F = F for simplicity in vacuum while raising or lowering the indices by means of the metric, we have the parametrization d i A j − d j A i = F i j and obtain by composition in the left upper square:

d i ( d i A j − d j A i ) = d i d i A j − d j ( d i A i ) = J j ⇒ d j ( d i d i A j − d j d i A i ) = d j J j = 0

with implicit summations on i and j. However, such a result does not prove at all that the density of current does not satisfy other CC. Nevertheless, we have:

LEMMA 4.B.1: The system d i d i A j − d j ( d i A i ) = J j is involutive whenever d j J j = 0 but the differential module defined by the corresponding homogeneous system is not torsion-free.

Proof: When J = 0 , this second order system is formally integrable because it is homogeneous. However, even if we know a priori that necessarily d j J j = 0 , it is not evident that such a condition is also sufficient, contrary to what is claimed in the literature. When n = 4 , using the Euclidean metric for simplicity, one can rewrite the system in the form:

{ d 44 A 3 + d 33 A 3 + d 22 A 3 + d 11 A 3 − d 3 ( d 4 A 4 + d 3 A 3 + d 2 A 2 + d 1 A 1 ) = J 3 d 44 A 2 + d 33 A 2 + d 22 A 2 + d 11 A 3 − d 2 ( d 4 A 4 + d 3 A 3 + d 2 A 2 + d 1 A 1 ) = J 2 d 44 A 1 + d 33 A 1 + d 22 A 1 + d 11 A 3 − d 1 ( d 4 A 4 + d 3 A 3 + d 2 A 2 + d 1 A 1 ) = J 1 d 34 A 3 + d 24 A 2 + d 14 A 1 − d 33 A 4 − d 22 A 4 − d 11 A 4 = − J 4 1 2 3 4 1 2 3 4 1 2 3 4 1 2 3 •

Let us check that the second order symbol is involutive with three equations of class 4 and only one equation of class 3. Indeed, we have successively:

d 344 A 3 + d 244 A 2 + d 144 A 1 − d 334 A 4 − d 224 A 4 − d 114 A 4 = − d 4 J 4 − d 344 A 3 − d 333 A 3 − d 223 A 3 − d 113 A 3 + d 334 A 4 + d 333 A 3 + d 233 A 2 + d 133 A 1 = − d 3 J 3 − d 244 A 2 − d 233 A 2 − d 222 A 2 − d 112 A 2 + d 224 A 4 + d 223 A 3 + d 222 A 2 + d 122 A 1 = − d 2 J 2 − d 144 A 1 − d 133 A 1 − d 122 A 1 − d 111 A 1 + d 114 A 4 + d 113 A 3 + d 112 A 2 + d 111 A 1 = − d 1 J 1

Summing the four equations, the left member is vanishing and we get the desired only CC. The four characters are ( 1 < 7 < 12 < 16 ) and we check that d i m ( g 2 ) = 1 + 7 + 12 + 16 = 36 = ( 10 × 4 ) − 4 . This procedure cannot be avoided though it cannot be found in the literature.

In addition, rewriting the homogeneous system in the form d i d i A j − d i d j A i = 0 with implicit summation on i, we obtain successively ( [

d i d i F j k = d i d i ( d j A k − d k A j ) = d j ( d i d i A k ) − d k ( d i d i A j ) = d j ( d i d k A i ) − d k ( d i d j A i ) = d j d k ( d i A i ) − d k d j ( d i A i ) = 0

It follows that each component of the field is a torsion element of the corresponding differential module M, which is killed by the Dalembertian. Accordingly, M is not torsion-free and thus not pure because we have just proved that r k D ( M ) = 1 . In particular, each component of the potential is free by itself. Such a situation is similar to the one met in the study of Single Input/Single Output (SISO) ordinary differential control systems. For example, if we have y x x − u x = 0 , then z = y x − u is a torsion element with z x = 0 even though both the input u and the output y are free separately. Bringing the system to first order by choosing the new unknowns ( z 1 = y , z 2 = y x , z 3 = u ) , we obtain the Spencer form over K = ℚ :

z x 1 − z 2 = 0 , z x 2 − z x 3 = 0

Setting a new z ¯ 1 = z 1 , z ¯ 2 = z 2 − z 3 , z ¯ 3 = z 3 , we obtain the reduced Spencer form:

z ¯ x 1 − z ¯ 2 − z ¯ 3 = 0 , z ¯ x 2 = 0 ⇔ z ¯ = A z + B u ⇒ r k K ( B , A B ) = 1 < 2

which is a Kalman form because no jet of the input z ¯ 3 = z 3 = u is indeed appearing.

Finally, the character α 2 n is obtained by considering d n n A j − d j n A n for the equation giving J j . For J n we get d n n A n − d n n A n = 0 and thus α 2 n = n − ( n − 1 ) = 1 a result showing that the corresponding differential module has rank 1 and there is thus only one CC, namely d j J j = 0 with implicit summation on j.

Q.E.D.

We now prove that we may add the Lorenz condition d i A i = 0 to bring the rank to zero. Indeed, we have now the inhomogeneous system d i d i A j = J j and the differential constraint thus brought is compatible with the conservation of current. The corresponding homogeneous system obtained by adding the Lorenz constraint has second order symbol obtained by considering both d i d i A j = 0 and d i j A i = 0 or d i j A j = 0 . We obtain therefore d n n A j = 0 , d n n A n = 0 showing that we have now α 2 n = 0 and a torsion differential module. As a more important and effective result that does not seem to be known, we have:

PROPOSITION 4.B.2: When n = 4 , the system:

{ Ψ j ≡ d 44 A j + ⋯ + d 11 A j = J j d 4 Φ − Ψ 4 ≡ d 34 A 3 + d 24 A 2 + d 14 A 1 − d 33 A 4 − d 22 A 2 − d 11 A 4 = − J 4 d 3 Φ ≡ d 34 A 4 + ⋯ + d 13 A 1 = 0 d 2 Φ ≡ d 24 A 4 + ⋯ + d 12 A 1 = 0 d 1 Φ ≡ d 14 A 4 + ⋯ + d 11 A 1 = 0 Φ ≡ d 4 A 4 + ⋯ + d 1 A 1 = 0 1 2 3 4 1 2 3 • 1 2 3 • 1 2 • • 1 • • • • • • •

is involutive with four equations of class 4, two equations of class 3, one equation of class 2 and one equation of class 1. The 11 resulting CC only provide the conservation of current.

Proof: Using the corresponding Janet tabular on the right, one can check at once that the 4 CC brought by the only first order equation Φ = 0 do not bring anything new, as they amount to crossed derivatives, and that we are only left with the 4 upper dots on the right side. However, for i = 1 , 2 , 3 , we have d 4 ( d i Φ ) = d i ( d 4 Φ − Ψ 4 ) + d i Ψ 4 and we are thus only left with a single CC, getting successively:

d 4 ( d 4 Φ − Ψ 4 ) ≡ d 344 A 3 + d 244 A 2 + d 144 A 1 − d 334 A 4 − d 224 A 2 − d 114 A 1 − d 3 Ψ 3 ≡ − d 344 A 3 − d 333 A 3 − d 223 A 3 − d 113 A 3 − d 2 Ψ 2 ≡ − d 244 A 2 − d 233 A 2 − d 222 A 2 − d 112 A 2 − d 1 Ψ 1 ≡ − d 144 A 1 − d 133 A 1 − d 122 A 1 − d 111 A 1

d 3 ( d 3 Φ ) ≡ d 334 A 4 + d 333 A 3 + d 233 A 2 + d 133 A 1 d 2 ( d 2 Φ ) ≡ d 224 A 4 + d 223 A 3 + d 222 A 2 + d 122 A 1 d 1 ( d 1 Φ ) ≡ d 114 A 4 + d 113 A 3 + d 112 A 2 + d 111 A 1

Summing these 7 equations, we are left with the identity − ( d 44 Φ + ⋯ + d 11 Φ ) + d j Ψ j = d j J j = 0 . It is important to notice that no other procedure can prove that we have an involutive symbol in δ -regular coordinates and this is the only way to compute effectively all the four characters ( 0 < 6 < 11 < 15 ) with 6 + 11 + 15 = 32 = ( 4 × 10 ) − ( 4 + 4 ) for the dimension of the symbol of order 2, a result not evident at first sight. Accordingly, the so-called Lorenz gauge condition is only a pure “artifact” amounting to a relative minimum parametrization with no important physical meaning as it can be avoided by using only the EM field F ( [

Q.E.D.

REMARK 4.B.3: Using the exactness of the Euler-Poincaré characteristic for the Poincaré sequence and its adjoint sequence, we have:

0 = 1 − n + n ( n − 1 ) / 2 − r k D ( d = M 1 ) ⇒ r k D ( d = M 1 ) = ( n − 1 ) ( n − 2 ) / 2

0 = 1 − n + n ( n − 1 ) / 2 − r k D ( d = M 2 ) ⇒ r k D ( a d ( d ) = M 2 ) = ( n − 1 ) ( n − 2 ) / 2

by counting the ranks from left to right in the upper row and then from right to left in lower row of the previous diagram. We obtain therefore:

r k D ( d ) + r k D ( d = M 1 ) = r k D ( a d ( d ) = M 2 ) + r k D ( a d ( d ) ) = 3 + 3 = d i m ( ∧ 2 T * ) = 6

when n = 4 by using successively the exactness of the upper and lower rows of the diagram. Hence, the image of the dotted arrow on the left is equal to the image of a d ( d ) = M 2 , a result explaining why both operators have the same CC, namely the conservation of current obtained by extending the previous diagram one step more to the left. As the linear Spencer sequence for a Lie group action (isometries, isometries plus dilatation, conformal isometries in Example 4.1) is (locally) isomorphic to the tensor product of the Poincaré sequence by a finite dimensional Lie algebra ( 10 < 11 < 15 when n = 4 in Example 4.1), the same comment remains valid. This result justifies by itself the specific importance of the linear Spencer sequence for infinitesimal Lie equations.

Such a new approach to a classical result is nevertheless bringing a totally unsatisfactory consequence. Using the well known correspondence between electromagnetism (EM) and elasticity (EL) used for all engineering computations with finite elements:

E M p o t e n t i a l ↔ E L d i s p l a c e m e n t , E M f i e l d ↔ E L s t r a i n , E M i n d u c t i o n ↔ E L s t r e s s

where EL means elasticity, and instead of the left upper square in the diagram, ... we have to consider the right upper square.

We finally prove that the use of the linear and nonlinear Spencer operators drastically changes the previous standard procedure in a way that could not even be imagined with classical methods. For such a purpose, we make a few comments on the implicit summation appearing in differential duality. For example, we have, up to a divergence:

X k , r X , r k = X k , r ( ∂ r ξ k − ξ r k ) = − ∂ r ( X k , r ) ξ k − X k , r ξ r k + ⋯

In the conformal situation, we have ξ 1 1 = ξ 2 2 = ⋯ = ξ n n = 1 n ξ r r and obtain therefore, as factor of the firs jets:

X 1 ,1 ξ 1 1 + X 2 ,2 ξ 2 2 + ⋯ + X n , n ξ n n = ( X r , r ) 1 n ξ r r = ( X r , r ) ξ 1 1

Going to the next order, we get as in ( [

X 1 1, r ∂ r ξ 1 1 = − ( ∂ r X 1 1, r ) ξ 1 1 + ⋯

Collecting the results and changing the sign, we obtain for the first time the Cosserat equation for the dilatation, namely the so-called virial equation that we provided in 2016 ( [

∂ r X 1 1, r + X r , r = 0

generalizing the well known Cosserat equations for the rotations provided in 1909 ( [

∂ r X i j , r + X i , j − X j , i = 0

As for EM, substituting ∂ i ξ r j r − ∂ j ξ r i r in the dual sum F i < j F i < j = 1 2 F i j F i j and integrating by parts, we get a part of the Cosserat equations for the elations, namely:

∂ r F i r − J i = 0 ⇒ ∂ i J i = 0 ⇒ X r , r = 0

saying that the trace of the EM impulsion-energy tensor must vanish ( [

REMARK 4.B.4: When g q + 1 = 0 , introducing the linear Spencer sequence with Spencer bundles C r = ∧ r T * ⊗ R ^ q , we obtain the following diagram describing the Cosserat procedure ( [

potential field C 0 → D 1 C 1 → D 2 C 2 ⋮ ↕ ⋮ ∧ n T * ⊗ C 0 * ← a d ( D 1 ) ∧ n T * ⊗ C 1 * ← a d ( D 2 ) ∧ n T * ⊗ C 2 * ↕ ↕ ↕ ∧ n T * ⊗ M q ← a d ( D 1 ) ∧ n − 1 T * ⊗ M q ← a d ( D 2 ) ∧ n − 2 T * ⊗ M q induction pseudopotential

where we have used the isomorphism R q * ≃ M q . It just remains to consider the various Spencer bundles C r ⊂ C ˜ r ⊂ C ^ r that we have already considered with q = 2 (see [

We sum up all these results in the following tabular only depending on the Spencer operator:

We notice that F = ( F i j ) ∈ ∧ 2 T * , R = ( R i j ) ∈ S 2 T * and the part of the first Spencer bundle made by the 1-forms with value in the n = 4 elations provides the splitting:

( F , R ) ∈ ∧ 2 T * ⊕ S 2 T * ≃ T * ⊗ T *

because of the well known exactness of the Spencer δ -sequence:

0 → S 2 T * → δ T * ⊗ T * → δ ∧ 2 T * → 0

As a byproduct, the surprising shift of exterior powers is simply obtained by applying the Spencer δ -map to C ^ r / C ˜ r = ∧ r T * ⊗ g ^ 2 ≃ ∧ r T * ⊗ T * because δ ( ∧ r T * ⊗ T * ) ⊂ ∧ ( r + 1 ) T * .

C) General Relativity:

Roughly speaking, we shall say that a parametrization of an operator is minimal if its corresponding operator defines a torsion module or, equivalently, if the kernel of the parametrizing operator has differential rank equal to 0. It is not so well known even today that, up to an isomorphism, the Cauchy stress operator essentially admits only one parametrization in dimension n = 2 which is minimum but the situation is quite different in dimension n = 3 . Indeed, the parametrization found by E. Beltrami in 1892 with 6 potentials ( [

With standard notations, denoting by Ω ∈ S 2 T * a perturbation of the non-degenerate metric ω , it is well known (see [

2 R i j = ω r s ( d i j Ω r s + d r s Ω i j − d r i Ω s j − d s j Ω r i ) = 2 R j i

Multiplying by test functions ( λ i j ) ∈ ∧ 4 T * ⊗ S 2 T and integrating by parts on space-time, we obtain the following four terms describing the so-called gravitational waves equations:

( □ λ r s + ω r s d i j λ i j − ω s j d i j λ r i − ω r i d i j λ s j ) Ω r s = σ r s Ω r s

where □ is the standard Dalembertian. Accordingly, we have:

d r σ r s = ω i j d r i j λ r s + ω r s d r i j λ i j − ω s j d r i j λ r i − ω r i d r i j λ s j = 0

The basic idea used in GR has been to simplify these equations by adding the differential constraints d r λ r s = 0 in order to find only □ λ r s = σ r s , exactly like in the Lorenz condition for EM. Before going ahead, it is important to notice that when n = 2 , the only Lagrange multiplier λ is just the Airy function ϕ and, using an integration by parts, we have the identity:

ϕ ( d 11 Ω 22 + d 22 Ω 11 − 2 d 12 Ω 12 ) = d 22 ϕ Ω 11 − 2 d 12 ϕ Ω 12 + d 11 ϕ Ω 22 + d i v ( )

providing the Airy parametrization of the Cauchy stress equations:

σ 11 = d 22 ϕ , σ 12 = σ 21 = − d 12 ϕ , σ 22 = d 11 ϕ

where the Airy function has, of course, nothing to do with any perturbation of the metric.

However, even if it is clear that the constraints are compatible with the Cauchy equations, we do believe that the following result is not known as it does

not contain any reference to the usual Einstein tensor E i j = R i j − 1 2 ω i j t r (R)

where t r ( R ) = ω r s R r s , which is therefore useless because it contains 6 terms instead of 4 terms only, even though the corresponding operator is self-adjoint.

PROPOSITION 4.C.1: The system made by □ λ r s = σ r s and d r λ r s = 0 is a relative minimum involutive parametrization of the Cauchy equations describing the formal adjoint of the Killing operator, that is Cauchy = a d ( Killing ) as operators.

Proof: For each given s = 1 , 2 , 3 , 4 the system under study is exactly the system used for studying the Lorenz condition in Proposition 4.B.1. Accordingly, nothing has to be changed in the proof of this proposition and we get an involutive second order system with d r σ r s = 0 as only CC in place of the conservation of current. Needless to say that this result has nothing to do with any concept of gauge theory as it is sometimes claimed ( [

Q.E.D.

In 1916, F.S. Macaulay used a new localization technique for studying unmixed polynomial ideals. In 2012, we have generalized this procedure in order to study pure differential modules, obtaining therefore a relative parametrization in place of the absolute parametrization already known for torsion-free modules and equivalent to controllability in the study of OD or PD control systems. Such a result is showing that controllability does not depend on the choice of the control variables, despite what engineers still believe. Meanwhile, we have pointed out the existence of minimum parametrizations obtained by adding, in a convenient but generally not intrinsic way, certain compatible differential constraints on the potentials. We have proved that this is exactly the kind of situation met in control theory, in EM with the Lorenz condition and in GR with gravitational waves. However, the systematic use of adjoint operators and differential duality is proving that the physical meaning of the potentials involved has absolutely nothing to do with the one usually adopted in these domains. Therefore, these results bring the need to revisit the mathematical foundations of Electromagnetism and Gravitation, thus of Gauge Theory and General Relativity, in particular Maxwell and Einstein equations, even if they seem apparently well established.

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

Pommaret, J.-F. (2019) Differential Homological Algebra and General Relativity. Journal of Modern Physics, 10, 1454-1486. https://doi.org/10.4236/jmp.2019.1012097

K differential field containing ℚ with n commuting derivations ∂ 1 , ⋯ , ∂ n .

d 1 , ⋯ , d n formal derivatives acting on the m differential indeterminates y 1 , ⋯ , y m .

D = K [ d 1 , ⋯ , d n ] = K [ d ] ring of differential operators P , Q with coefficients in K.

L , M , N filtered differential modules over D.

μ = ( μ 1 , ⋯ , μ n ) multi-index with length | μ | = μ 1 + ⋯ + μ n and μ + 1 i = ( μ 1 , ⋯ , μ i + 1 , ⋯ , μ n ) .

X manifold with tangent, cotangent, symmetric and exterior bundles T , T * , S q T * , ∧ r T * .

y q = ( y k , y i k , y i j k , ⋯ ) = ( y μ k ) , 0 ≤ | μ | ≤ q jet coordinates up to order q.

a d ( D ) formal adjoint defined by 〈 λ , D ξ 〉 = 〈 a d ( D ) λ , ξ 〉 + d i v ( ) for any test row vector λ .

E , F vector bundles over X or free modules over D.

V characteristic algebraic variety with dimension d and codimension c d = n − d over K.

0 = t n ( M ) ⊆ ⋯ ⊆ t 0 ( M ) ⊆ M purity filtration with t r ( M ) = { m ∈ M | c d ( D m ) > r } .

t 0 ( M ) = t ( M ) = { m ∈ M | ∃ 0 ≠ P ∈ D , P m = 0 } the torsion submodule of M.

R = h o m K ( M , K ) differential module over D associated with M, also called inverse system.

d : R → T * ⊗ K R = f → d x i ⊗ d i f Spencer operator with ( d i f ) μ k = ∂ i f μ k − f μ + 1 i k .

T * ⊗ X T * ≃ S 2 T * ⊕ ∧ 2 T * with R = ( R i j ) ∈ S 2 T * Ricci tensor and F = ( F i j ) ∈ ∧ 2 T * EM field.