LettreMartial Simonin à Henri Poincaré, 12 avril 1896

Nice le 12 Avril 1896

Cher Monsieur,

Je me permets de vous adresser les résultats que j’ai obtenus au sujet de l’orbite d’Hécube (108).

D’après les indications que vous avez bien voulu me donner, j’ai introduit une constante de plus dans mes formules qui n’en contenaient que 3 au lieu de 4. De plus j’ai tenu compte, pour les termes principaux, du carré de l’excentricité ; j’ai ainsi obtenu, dans les expressions du périhélie et de l’excentricité, le terme en \(\begin{array}{l}\sin\\\cos\end{array}(J - \overline{3\theta_0 - g_0}t)\) que j’avais introduit empiriquement au mois d’octobre dernier pour diminuer les résidus.1

J’ai résolu à nouveau les équations qu’on obtient par la méthode d’identification de Delaunay; pour cela je me suis servi de ce que j’avais une valeur approché de chacune de ces inconnues pour négliger certains termes et en conserver d’autres, sans tenir compte de l’exposant de la masse perturbatrice.

En outre dans l’expression du demi grand-axe, j’ai introduit un terme à longue période contenant le carré de la masse. Enfin j’ai déterminé les perturbations périodiques les plus importantes du demi grand-axe, de la longitude moyenne, du périhélie et de l’excentricité. En se servant des formules de Le Verrier, on trouve les coefficients suivants exprimés en secondes d’arc :

pour la longitude:

\(366''\sin(\lambda'-\lambda)\), \(37''\sin(\lambda'-\varpi)\), \(-34''\sin(\lambda'-\varpi')\), \(77''\sin(5\lambda'-3\lambda-2\varpi)\)
\(441\sin2(\lambda'-\lambda)\), \(52\sin(\lambda'-\varpi)\), \(228 \sin(3\lambda'-2\lambda-\varpi')\) \(-112\sin(5\lambda'-3\lambda-\varpi-\varpi')\)
\(182\sin3(\lambda'-\lambda)\), \(-344\sin(3\lambda'-2\lambda-\varpi)\) \(83\sin(4\lambda'-3\lambda-\varpi')\)  
\(86\sin4(\lambda'-\lambda)\), \(-123\sin(4\lambda'-3\lambda-\varpi)\) \(46\sin(5\lambda'-4\lambda-\varpi')\)  
\(44\sin5(\lambda'-\lambda)\), \(-61\sin(5\lambda'-4\lambda-\varpi)\)    

et pour l’excentricité :

\(-39''\cos(\lambda - \varpi)\), \(82''\cos(3\lambda'-\lambda-2\varpi)\)
\(-83\cos(\lambda'-\varpi)\), \(-59\cos(5\lambda'-3\lambda-2\varpi)\)
\(189\cos(3\lambda'-2\lambda-\varpi)\), \(-53\cos(3\lambda'-\lambda-\varpi-\varpi')\)
\(70\cos(4\lambda'-3\lambda-\varpi)\), \(43\cos(5\lambda'-3\lambda-\varpi-\varpi')\)
\(32\cos(5\lambda'-4\lambda-\varpi)\),  

Le terme à longue période contenant le carré de l’excentricité de l’orbite de Jupiter a toujours été négligé; par suite je n’ai pas tenu compte des perturbations à courte période contenant ce carré ; le plus fort des termes ainsi négligés a pour coefficient 36” de longitude.

Les formules qui m’ont servi sont, en ne transcrivant pas les perturbations à courte période et en posant : \(\mathcal{L} = \sqrt{a}\), \(\lambda=\) longitude moyenne, \(\xi = \sqrt{L}e\cos\varpi\), \(\eta = \sqrt{L}e\sin\varpi\),2 \[\begin{aligned} \mathcal{L}&=G_0 + A - \frac{G_0}{2}B^2 - \frac{E_1^2}{2}\cos^2\theta_0t - BE_1G_0\cos\theta_0t + B'E_1G_0\sin\theta_0t\\ &\quad + A_1\cos\left(J - \overline{\theta_0 - g_0}t - g' - h'\right),\\ \lambda&=J + nt + A\frac{\partial(\theta_0 - g_0)}{\partial G_0}t + \frac{3}{2G_0^3}B^2t + B(\theta_1 - g_1)\sin\theta_0t + B'(\theta_1 - g_1)\cos\theta_0t\\ &\quad + B^2(\theta_2 - g_2)\sin 2\theta_0t - C + B_1\sin(J - \overline{\theta_0 - g_0}t - g' - h'),\\ \xi&=\left[E_1\sqrt{G_0} + A\frac{\partial(E_1\sqrt{G_0})}{\partial G_0} + \sqrt{G_0}E'B^2\right]\cos(J - \overline{\theta_0 - g_0}t)\\ &\quad + \left[AE_1\sqrt{G_0}\frac{\partial(\theta_0 - g_0)}{\partial G_0} + \frac{3e_1\sqrt{G_0}}{2G_0^3}B^2\right]t\sin(J - \overline{\theta_0 - g_0}t)\\ &\quad + B\sqrt{G_0}\cos(J + g_0t) - B'\sqrt{G_0}\sin(J + g_0t) + BE_2\sqrt{G_0}\cos(J - \overline{2\theta_0 - g_0}t)\\ &\quad + B^2E_3\sqrt{G_0}\cos(J - \overline{3\theta_0 - g_0}t) + C_1\cos(g' + h'),\\ \eta&=\left[E_1\sqrt{G_0} + A\frac{\partial(E_1\sqrt{G_0})}{\partial G_0} + \sqrt{G_0}E'B^2\right]\sin(J - \overline{\theta_0 - g_0}t)\\ &\quad - \left[AE_1\sqrt{G_0}\frac{\partial(\theta_0 - g_0)}{\partial G_0} + \frac{3e_1\sqrt{G_0}}{2G_0^3}B^2\right]t\cos(J - \overline{\theta_0 - g_0}t)\\ &\quad + B\sqrt{G_0}\sin(J + g_0t) + B'\sqrt{G_0}\cos(J + g_0t) + BE_2\sqrt{G_0}\sin(J - \overline{2\theta_0 - g_0}t)\\ &\quad + B^2E_3\sqrt{G_0}\sin(J - \overline{3\theta_0 - g_0}t) + C_1\sin(g' + h').\\\end{aligned}\] Vu la difficulté du calcul des coefficients \(\theta_1\) \(g_1\) \(\theta_2\) \(g_2\) etc. …, j’ai conservé ces notations malgré leur complication. Les valeurs numériques des constantes de ces formules sont :

\(n=613'',576\) \(J=166^{\circ},59'\) \(\theta_0=16'',176\) \(A_1=[\overline{3}.167997]_n\)
\(G_0=[0,842200]\) \(E_1=[\overline{2}.363093]\) \(g_0=0'',856\) \(B_1=[\overline{2}.394700]_n\)
\(A=[\overline{2}.527508]\) \(E'=[0,008790]_n\) \(\theta_1-g_1=[0,439290]\) \(C_1=[\overline{2}.662960]_n\)
\(B=[\overline{2}.988653]\) \(E_2=[\overline{3}.080174]_n\) \(\theta_2-g_2=[\overline{1}.744784]_n\)  
\(B'=[\overline{3}.888869]\) \(E_3=[\overline{1}.373974]_n\) \(\frac{\partial(\theta_0-g_0)}{\partial G_0}=[\overline{3}.093583]_n\)  
\(C=[\overline{2}.211722]\)      

J’ai mis entre crochets [ ], au lieu des nombres, leurs logarithmes. On obtient ainsi les éléments d’Hécube rapportés à l’écliptique et à l’équinoxe de 1850,0; l’époque est 1897 septembre 23,5, temps moyen de Berlin.3

Je rappelle rapidement qu’on a posé \(\theta=\lambda-2\lambda'\), puis \[\begin{aligned} \theta&=&\theta_0(t+c) + \theta_1\sin\theta_0(t+c) + \theta_2\sin 2\theta_0(t+c)\\ g&=&g_m + g_0(t+c) + g_1\sin\theta_0(t+c) + g_2\sin 2\theta_0(t+c)\\ e&=&e_0 + e_1\sin\theta_0(t+c) + e_2\sin 2\theta_0(t+c) + e_3\sin 3\theta_0(t+c)\\\end{aligned}\] \(e_1\), \(e_2\), \(e_3\) peuvent se mettre sous la forme \(\begin{array}{|rcl}e_1&=&E_1 + E'e_0^2\\ e_2&=&E_2e_0 + \frac{E_2'}{e_0}\\ e_3&=&E_3e_0^2 + \cdots\\ \end{array}\)

Les formules qui donnent \(\Omega\), longitude du nœud ascendent et \(i\), inclinaison de l’orbite, sont: \[\begin{gathered} \operatorname{tg}(\Omega - J + \overline{\theta_0 - g_0}t) =[0,001810]\operatorname{tg}\left[ [\overline{1}.996483](353^{\circ}15'9'' - J +\overline{\theta_0 - g_0}t)\right]\\ i=[\overline{2}.886269]\left\{[\overline{1}.996486] - [\overline{3}.616297]\cos[0,297543]\{353^{\circ}15'9'' - J + \overline{\theta_0 - g_0}t\}\right\}^{\frac{1}{2}}\end{gathered}\]

En comparant les ascensions droites et les déclinaisons données par les observations et par ces formules, on obtient les différences: \(o-c\)4

Année \(\alpha\) \(\delta\)
1869 \(-7^s\) \(+68''\)
71 \((-20)\) \(-165\)
74 \(+9\) \(+54\)
75 \(+8\) \(-67\)
76 \(+4\) \(+6\)
77 \(+3\) \(+38\)
78 \(-9\) \(-18\)
80 \(-6\) \(+43\)
81 \(+5\) \(-48\)
86 \(-1\) \(+8\)
88 \(+6\) \(+77\)
89 \(-9\) \(-43\)
92 \((-18)\) \(+120\)
1894 \(+7\) \(+71\)

Etant donnés ces résidus, au lieu de calculer des lieux normaux, j’ai conservé, pour chaque opposition, une observation fictive, moyenne de plusieurs observations faites le même soir par divers observateurs.

On peut remarquer que les observations de 1871 et de 1892 sont distantes d’environ 21 ans, et que les arguments \(3\lambda'-2\lambda\), \(3\lambda'-\lambda\), \(5\lambda'-3\lambda\) et \(\lambda'\) ont pour période environ 11 ans. Je crois devoir attribuer ces deux résidus, non pas aux termes périodiques que j’ai négligés, mais plutôt à la suppression de tous les termes de \(e^3\). D’ailleurs monsieur Perrotin5 a bien voulu me faire remarquer que dans le \(n^\circ\) 368 de l’astronomical Journal p. 62, M Hill a trouvé pour la longitude de moyenne de Cérés des résidus de \(-40''\) en 1857 et \(+40''\) en 1866, quoique les éléments de cette planète ne présentent rien de particulier.6

Comme je ne pouvais diminuer mes résidus en changeant le moyen mouvement, l’époque ou la masse, et que j’avais rencontré de grandes difficultés dans les calculs des divers coefficients numériques, j’ai cherché des coefficients empiriques; le premier terme de \(\xi\) et de \(\eta\) difficile à calculer donnerait de meilleurs résultats si on lui ajoutait le coefficient empirique: \([\overline{3}.231551]_n\). On peut remarquer aussi que le terme \(B^2E_3\sqrt{G_0}\cos(J - \overline{3\theta_0 - g_0}t)\) est plus important que le terme \(BE_2\sqrt{G_0}\cos(J - \overline{2\theta_0 - g_0}t)\). En outre on voit aisément que les observations de 1871–74–75 seraient mieux représentées si on diminuait la longitude du périhélie, et celles de 1892–94, si on l’augmentait. L’introduction du terme à longue période contenant le carré de l’excentricité de l’orbite de Jupiter ne donnerait pas de meilleurs résultats.

En résumé je ne vois pas comment avec les formules données plus haut, on peut obtenir des résidus inférieurs à ceux que j’ai transcrits ci-dessus. Je recours donc encore une fois à vos bienveillants conseils, trop heureux si vous êtes un peu satisfait des efforts que j’ai faits et des résultats obtenus.

Je vous prie de vouloir bien excuser la longueur de cette lettre.

Daignez agréer, cher Monsieur, l’expression de toute ma gratitude et de mon entier dévouement,

M. Simonin


Apparat critique 

  1. Dans un premier temps, Simonin se sert des variables classiques : \[L = \sqrt{a}, \qquad G = \sqrt{a(1-e^2)},\qquad H = G \cos i, \qquad \gamma = \sin \frac{i}{2},\]\(a\) est le demi grand-axe, \(e\) l’excentricité, \(i\) l’inclinaison. Après avoir exprimé les équations du système avec ces variables et obtenu une série de formules permettant de calculer \(e\) et une expression \(\theta\) égale à \[l+2g+2h-2l'-2g'-2h'\]\(l\) est l’anomalie moyenne, \(g\) la distance du nœud ascendant au périhélie et \(h\) la longitude du nœud, Simonin effectue un changement de variables qui offre, entre autres avantages, celui de faire disparaître \(e_0\) du dénominateur Simonin (1897, p. 17) : \[\left\{ \begin{array}{r@{\; = \;}l} L&L,\\ \lambda & l+g+h,\\ \eta & \sqrt{L} \, e \sin(g+h),\\ \xi & \sqrt{L} \, e \cos(g+h), \end{array} \right. .\] Simonin signale que les formules obtenues pour ces quatres variables donnent pour employer le langage de M. Poincaré, toutes les solutions périodiques du problème restreint. Ces formules dépendent de quatre constantes \(G_0\) une constante arbitraire qui apparaît lorsque l’on considère la première du système \(2L-G=G_0\), \(e_0\), \(g_m\) et \(c\) des constantes d’intégration. En posant \(J = g_m + h\) et \(n = \theta_0-g_0+2n'\), Simonin obtient une solution périodique de la première sorte (au sens de Poincaré) : \[\left\{ \begin{array}{r@{\; = \;}l} L&G_0,\\ \lambda & J + nt,\\ \eta & e_{11}\sqrt{G_0} \sin [J-(\theta_0 - g_0)t],\\ \xi & e_{11}\sqrt{G_0} \cos [J-(\theta_0 - g_0)t], \end{array} \right. .\] Simonin poursuit son calcul en suivant la stratégie de Poincaré :

    Dès qu’on connaît une solution périodique du problème, toutes celles qui en diffèrent peu s’obtiennent par l’intégration d’un système d’équations d’un système d’équations linéaires et homogènes ; ce sont ces équations que M. Poincaré a appelées équations aux variations. Simonin (1897, p. 20)

    Le traitement des équations aux variations fait apparaître des termes en \(\begin{array}{l}\sin\\\cos\end{array}[J - (\theta_0 - g_0)t]\), \(\begin{array}{l}\sin\\\cos\end{array}[J - (2\theta_0 - g_0)t]\) et \(\begin{array}{l}\sin\\\cos\end{array}[J - (3\theta_0 - g_0)t]\) dans les séries de \(\eta\) et \(\xi\) (Voir Simonin (1897, p. 24).↩︎

  2. Ces formules sont celles que Simonin appelle les expressions définitives de \(L\), \(\lambda\), \(\eta\) et \(\xi\) Simonin (1897, p. 30). Comme le rappelle Simonin, \(A,B,B',C\) sont des constantes d’intégration qu’il faut déterminer à partir des observations et les autres coefficients sont des fonctions de \(G_0\).↩︎

  3. Les valeurs numériques obtenues dans la thèse Simonin (1897, p. 53) sont sensiblement différentes.↩︎

  4. \(o-c\) désigne les différences entre les valeurs observées et celles calculées.↩︎

  5. Joseph Perrotin était depuis 1884 le directeur de l’observatoire de Nice dans lequel travaillait Perrotin.↩︎

  6. Hill (1896.)↩︎


Références

Hill, G. W. 1896 Jupiter-perturbations of Ceres of the first order and the derivation of the mean elements. Astronomical Journal 16, pp. 57–62. Lien externe. ↩︎

Simonin, M 1897 Sur l’orbite de (108) Hécube. Ph.D. Thesis, Faculté des sciences de Paris, Paris. Lien externe. ↩︎

Titre (dcterms:title)

Martial Simonin à Henri Poincaré, 12 avril 1896

Incipit (ahpo:incipit)

Je me permets de vous adresser les résultats que j'ai obtenus ...

Date (ahpo:writingDate)

1896-04-12

Expéditeur (ahpo:sentBy)

Destinataire (ahpo:sentTo)

Adresse (ahpo:destinationAddress)

Lieu (ahpo:writtenAt)

Chapitre (ahpo:publishedIn)

Lieu d’archivage (ahpo:archivedAt)

Type (ahpo:documentType)

(fr) Lettre autographe signée

Section (dans le livre) (ahpo:sectionNumber)

2

Droits (ahpo:rightsHolder)

Archives Henri Poincaré

Nombre de pages (ahpo:numberOfPages)

7

Noms cités dans l'apparat (ahpo:citeApparatName)

Langue (ahpo:language)

fr

Éditeur (dcterms:publisher)

Licence (dcterms:license)

« Martial Simonin à Henri Poincaré, 12 Avril 1896 ». La Correspondance Entre Henri Poincaré, Les Astronomes Et Les géodésiens. Archives Henri Poincaré, s. d., Archives Henri Poincaré, s. d, La correspondance d'Henri Poincaré, accessed 9 August 2020, http://henripoincare.fr/s/correspondance/item/8373

Scan