La description du mouvement des planètes par Johannes Kepler à partir des tables établies par son maître Tycho Brahe est un tournant dans l'histoire des sciences qui prit naissance au début du XVIIe siècle. Malgré le retentissant procès de Galilée (qui doit abjurer sa conception copernicienne de l'univers le 21 juin 1632), le philosophe et mathématicien René Descartes verra, quelques années plus tard, le triomphe de sa méthode, d'après laquelle les sciences de la Nature reposent sur la méthode expérimentale.
Environ cinquante ans plus tard, la théorie sera écrite dans les Principia par Isaac Newton de fin 1684 à 1687, sous l'impulsion de son commanditaire et ami Edmond Halley. Ainsi naquit la mécanique céleste, ainsi que la mécanique rationnelle.
Le XVIIIe sera le siècle où l'on n'aura de cesse de clarifier cette théorie.
ébauche :
Ptolémée et l'équant
Copernic (1543) et la "révolution copernicienne" (cf Koyré)
Tycho Brahe, Uraniborg, et son "équivalence des systèmes"
ébauche :
Newton : les mesures de Richer et Picard (Guyane), de Halley (Sainte-Hélène), de Flamsteed (observatoire de Greenwich, les comètes)
Wallis, Hooke (1679), Halley-Wren-Hooke et le pari.
ébauche :
Galilée et les orbites circulaires, et les satellites médicéens.
Horrocks et le centre de "masse" ; Wallis et le centre de "gravité " du système Terre-Lune.
Depuis la découverte des trois lois de kepler sur le mouvement des planètes, il y a de nombreuses démonstrations relatives au mouvement keplerien dont les plus anciennes virent le jour au XVIIe siècle. En voici quelques-unes.
Laplace la reprendra dans son traité de « Mécanique Céleste ».
Que cela est vite dit dans notre langage moderne ! En réalité, la démonstration géométrique est la remarque classique sur le rôle des podaires dans le cas de champs centraux. Danjon remarque (avec Hamilton) que l'antipodaire de l'inverse d'un cercle est une conique : cela était enseigné encore au baccalauréat des années 1960 (Cf. LEBOSSÉ & EMERY, cours de mathématiques élémentaires).
Quant à Hermann, c'est un tour_de_force :
Il possède trois intégrales premières en coordonnées cartésiennes tirées de et idem en y.
Éliminer la vitesse : on trouve : c'est une ellipse (Cf. conique, Kepler). Mais comment a-t-il trouvé les deux intégrales premières du vecteur excentricité ? par un raisonnement analytico-géométrique horriblement compliqué ! On sait aujourd'hui le faire par la théorie de la représentation linéaire des groupes (Moser et SO(4) : 1968). On trouvera sur ce sujet un article de Guichardet dans la Gazette des mathématiciens (été 2008).
Le problème est plan, si la force est centrale. Le plan de phase est donc ( ). Les deux équations du PFD (principe fondamental de la dynamique) sont :
et la même en y. [Évidemment Ω dépend de r!]. Cette notation est évidemment très réminiscente de celle de Hooke. Mais elle n'a rien à voir, sinon que la symétrie est centrale.
Choisir trois fonctions invariantes par rotation :
Remarquer cette particularité : r² est choisie comme variable, et non r. Et comme J est non-nulle, I va jouer le rôle d'une échelle de temps au moins sur une demi-période, du périgée à l'apogée.
Démontrer que le problème se réduit au système différentiel (S) :
Keill utilise alors l'échelle de temps I ; le système se réduit à :
En éliminant Omega² (et quelle que soit sa valeur ! donc c'est vrai pour toute force centrale !)
.
C'est un vrai tour de force : au début du XVIIIème, on vient de réécrire :
Emmy Noether connaissait-elle cette démonstration due à l'invariance par rotation ?
Puis, l'invariance temporelle donne la conservation de l'énergie :
, où V(I) est l'énergie potentielle relative à la force centrale (=
Ces deux ensembles de surfaces feuillettent l'espace (I, J, K) et leur intersection donne l'orbite du mouvement dans cet espace.
Éliminer K conduit à travailler dans le demi-plan ( ), c'est-à-dire dans un plan de phase presque usuel (on joue avec r² plutôt qu'avec r) :
,
ce qui est l'équation de Leibniz (1689), mais en notation I = r². (Remarquer que tout résulte de cette circonstance, non évidente du temps où les vecteurs n'existaient pas) :
)
et pour finir, dt = dI/2J donne le mouvement sur cette orbite de phase et la primitive de 2J(I) donne l'action S(I) du problème.
Évidemment, actuellement, nous repasserions immédiatement en coordonnées (r et r'). Il n'empêche que voilà décrite la solution de Keill qui témoigne d'une virtuosité tombée dans l'oubli de l'Histoire. La suite est très classique et correspond à différents paramétrages dans le cas de Kepler. L'équation de Leibniz se réécrit dans ce cas :
qui est une conique en J et r, ellipse si H est négatif de grand axe :
Il est usuel alors de paramétrer via l'« anomalie excentrique » :
,
et « miraculeusement » :
,
qui s'intègre en donnant la fameuse équation de Kepler. En contrepartie l'équation en thêta est légèrement plus compliquée à intégrer (primitive de ) d'où :
.
et cette fois le paramétrage adéquat est :
et
ce qui conduit au « miraculeux » dθ = dα ! la trajectoire est donc une ellipse. Mais la deuxième intégration conduit à plus difficile à intégrer (mais tout à fait faisable !)
Partant de l'équation radiale de Leibniz (1689) :
il pose comme nouvelle variable z = C²-r et trouve :
,
identique aux deux équations de départ en x & y : donc il obtient z & x & y linéairement liés, ce qui est la définition d'une ellipse (Cf. conique, discussion). CQFD
Laplace en tire cette fois quatre équations linéaires identiques : d/dt(r^3.Z") = - Z', avec Z = r, x, y, cste. D'où r = a x + by + c.cste : c'est une conique !
|dU/dT|² = 8 GM + 8 E |U|² ; soit par dérivation , avec E négatif.
Donc U décrit une ellipse de Hooke et z =sqrt(U) l'ellipse de Kepler.
On aura reconnu en T(t), l'anomalie excentrique. Ce n'est donc qu'une des méthodes précédentes : mais cette méthode a des prolongements plus importants (Cf. théorème de Bertrand). Voir aussi plus bas.
Note : cette transformation du problème de Kepler en problème de Hooke est assez stupéfiante. Saari (p. 141) s'y attarde un peu plus qu'Arnold (BHHH) ; peut-être est-ce justifié ; voici :
Le problème de régularisation se pose s'il y a collision, c'est-à-dire, C très voisin de zéro. Saari dit : la collision entraîne un changement brutal de 2Pi . Afin de garder la particule sur la droite sans singularité, il "suffit de penser" à garder l'arc -moitié ; soit de changer de jauge (de fonction inconnue): et de variable (transmutation d'échelle de temps) dT = dt/r(t)(ATTENTION au facteur 4!)
La conservation de l'énergie s'écrit 2|U'|²-1 = Eo.r
et l'équation du mouvement : devient :
,
équation LINEAIRE sans le r^3 ! Elle conduit à :
U" -U/r [2|U'|²-1] =0
soit (équation de Hooke).
Le gros avantage de cette solution est qu'elle est stable-numérique : les solutions restent sur la même iso-énergie.
Kustaanheimo (1924-1997) et Stiefel (1909-1978) en 1964 utilisèrent enfin les quaternions pour transformer le problème de Kepler dans R^3 en celui de Hooke dans R^3, via R^4 ! (congrès d'Oberwolfach): ils leur a suffi de prendre la quatrième coordonnées x4 = cste : alors le quaternion U se déplaçait sur la sphère ; ceci mit en exergue la symétrie SO(4) et mieux SO(4,2) qui correspondait à la version spinorielle du problème de Kepler (liée à la solution en coordonnées paraboliques) et mettait en avant le vecteur excentricité. Immédiatement, le traitement des perturbations fut amélioré (Stiefel et Scheifel, 1971), mais aussi la quantification (méthode dite de Pauli (SO(4)), et surtout la quantification lagrangienne SO(4,2), aux orbitales "paraboliques" de Kleinert (1967-1998)(cf Kleinert 2006).
Saari donne des raisons topologiques à l'obstruction du passage de R^2 à R^3 et la nécessité de passer à R^4 (les quaternions) : la relation U^2 = z, ne pouvait se régulariser sur la sphère à cause du célèbre théorème du hérisson de Brouwer-Poincaré. Mais si on ne peut "peigner" S2 , on peut peigner S3 (et même S7 : octonions), ce qui avec les trois vecteurs tangents donne la fameuse matrice 4-4 de la transformation K-S : rappelons que le maître de Stiefel était Hopf lui-même qui dressa la carte de S3 vers S2 : il n'y a pas de hasard, la formation, cela sert ! (cf Oliver (2004)).