- Vecteur de Runge-Lenz
-
- Dans cet article les vecteurs et leurs normes sont indiqués respectivement en gras et italique. Par exemple : .
En mécanique classique, le vecteur de Runge-Lenz ou invariant de Runge-Lenz est un vecteur utilisé principalement pour décrire la forme et l'orientation de l'orbite d'un corps astronomique autour d'un autre, comme dans le cas d'une planète autour d'une étoile. Pour deux corps en interaction gravitationnelle, le vecteur de Runge-Lenz est une constante du mouvement (en), ce qui signifie qu'il prend la même valeur en n'importe quel point de l'orbite[1] ; de manière équivalente on dit que le vecteur de Runge-Lenz se conserve. Plus généralement le vecteur de Runge-Lenz est conservé pour n'importe quel problème à deux corps interagissant par le biais d'une force centrale variant comme l'inverse du carré de la distance entre eux. De tels problèmes sont appelés « problèmes de Kepler »[2].
L'atome d'hydrogène est un problème de Kepler puisqu'il comprend deux charges en interaction électrostatique, une autre force centrale en carré inverse de la distance. Le vecteur de Runge-Lenz fut essentiel dans les premières descriptions quantiques du spectre d'émission de l'atome d'hydrogène[3] après le développement de l'équation de Schrödinger. Cependant cette approche est aujourd'hui très peu utilisée. En mécaniques classique et quantique, les grandeurs conservées correspondent généralement à une symétrie du problème. La conservation du vecteur de Runge-Lenz est associée à une symétrie inhabituelle : le problème de Kepler est mathématiquement équivalent à une particule se déplaçant librement sur une 3-sphère[4], ce qui implique que le problème est symétrique pour certaines rotations dans un espace à quatre dimensions[5]. Cette symétrie supérieure résulte de deux propriétés du problème de Kepler : le vecteur vitesse se déplace toujours dans un cercle parfait et, pour une énergie totale donnée, tous les cercles de vitesse s'interceptent en deux mêmes points[6].
Le vecteur de Runge–Lenz est nommé d'après Carl Runge et Wilhelm Lenz. Il est également connu sous le nom de vecteur de Laplace (d'après Pierre-Simon de Laplace) bien qu'aucun de ces scientifiques ne l'ait découvert. Le vecteur de Runge-Lenz a en réalité été redécouvert plusieurs fois[7] et il est équivalent au vecteur excentricité (en) de la mécanique céleste[8]. Plusieurs généralisations du vecteur de Runge-Lenz ont été définies pour tenir compte de la relativité générale, du champ électromagnétique et des différents types de forces centrales.
Contexte
Une particule seule se déplaçant sous l'effet d'un force centrale conservative possède au minimum quatre constantes du mouvement : son énergie totale, et les trois composante du vecteur moment cinétique L. La trajectoire de la particule est contenue dans un plan défini par sa quantité de mouvement initiale p (ou de manière équivalente par sa vitesse v) et par le rayon-vecteur r entre la particule et le centre de force (voir ci-dessous Figure 1).
Comme défini ci-dessous (voir Définitions mathématiques), le vecteur de Runge-Lenz A est toujours contenu dans le plan de la trajectoire pour une force centrale. Cependant A n'est constant que pour une force centrale en carré inverse[1]. Pour de nombreuses forces ce vecteur A n'est pas constant mais change en norme et en direction ; si le champ de force centrale suit approximativement une loi en carré inverse, le vecteur A est approximativement constant en norme mais tourne lentement en direction. Un vecteur de Runge-Lenz généralisé se conservant peut être défini pour toutes les forces centrales mais le vecteur est alors une fonction compliquée de la position et n'est généralement pas exprimable de façon analytique[9],[10].
Le plan de la trajectoire est perpendiculaire au vecteur moment cinétique L qui est constant ; ceci peut être exprimé par le produit scalaire r·L = 0 ; de la même façon, A étant contenu dans ce plan : A·L = 0.
Histoire des redécouvertes
Le vecteur de Runge-Lenz A est une constante du mouvement du problème de Kepler, utile pour décrire les orbites astronomiques, telles que le mouvement des planètes. Néanmoins son usage n'a jamais été très répandu parmi les physiciens peut-être parce qu'il est moins intuitif que la quantité de mouvement ou le moment cinétique. C'est la raison pour laquelle il a été redécouvert plusieurs fois au cours des trois derniers siècles[7]. Jakob Hermann est le premier à montrer que A est conservé dans le cas particulier des forces centrales en carré inverse[11], et travailla sur ses connexions avec l'excentricité orbitale dans le cas des orbitales elliptiques. Les travaux de Hermann sont généralisés dans leur forme moderne par Jean Bernoulli en 1710[12]. À la fin du XVIIIe siècle, Pierre-Simon de Laplace redécouvre la conservation de A, l'obtenant de manière analytique et non plus de façon géométrique[13]. Vers le milieu du XIXe siècle, William Rowan Hamilton détermine le vecteur excentricité, qui est équivalent[8], et l'utilise pour montrer que le vecteur quantité de mouvement p décrit un cercle pour un mouvement effectué sous une force centrale en carré inverse (Figure 3)[6]. Au début du XXe siècle, Josiah Willard Gibbs trouve le même vecteur par analyse vectorielle[14]. La méthode de Gibbs est citée en exemple par Carl Runge dans un livre allemand sur les vecteurs[15] qui fut référencé par Wilhelm Lenz dans son article sur le traitement quantique de l'atome d'hydrogène[16]. En 1926, le vecteur est utilisé par Wolfgang Pauli pour déterminer le spectre de l'hydrogène en utilisant la mécanique matricielle et non l'équation de Schrödinger[3] ; après l'article de Pauli il commence à être connu sous le nom de vecteur de Runge–Lenz.
Définitions mathématiques
Pour une particule seule mise en mouvement par une force centrale en carré inverse de la distance décrite par l'équation , le vecteur de Runge-Lenz A est défini mathématiquement par la formule[1]
où
- est le masse de la particule ponctuelle se déplaçant sous l'effet de la force centrale,
- est le vecteur quantité de mouvement,
- est son vecteur moment cinétique,
- est un paramètre décrivant l'intensité de la force,
- est le vecteur position de la particule (Figure 1), et
- est le vecteur unitaire correspondant, i.e., où r est la norme de r.
Étant donné que la force est supposée conservative, l'énergie totale E est une constante du mouvement (en)
De plus la force étant centrale, le moment cinétique L est également constant et définit le plan dans lequel la particule se déplace. Le vecteur de Runge-Lenz est perpendiculaire à L car p ∧ L et r sont perpendiculaires à L. Il s'ensuit que A est contenu dans le plan de l'orbite.
Cette définition de A est applicable pour une particule seule et ponctuelle de masse m se déplaçant sous l'effet d'une force donnée. Cependant la même définition peut être étendue au problème à deux corps, tel que le problème de Kepler, en prenant pour m la masse réduite des deux corps et pour r le vecteur entre les deux corps.
Une alternative pour la même constante du mouvement est également possible. La plus commune est de diviser par mk pour définir le vecteur excentricité (en)
Obtention de l'orbite de Kepler
La forme et l'orientation des orbites dans le problème à deux corps de Kepler peuvent être déterminées à partir du vecteur de Runge-Lenz comme suit[1]. Le produit scalaire de A par le vecteur position r donne l'équation
où θ est l'angle entre r et A (Figure 2). En permutant le produit mixte
et un réarrangement conduit à la formule de définition d'une conique
et de paramètre p
Le demi grand-axe a de la conique peut être obtenu à partir du paramètre et l'excentricité
où le signe moins désigne une ellipse et le signe plus une hyperbole.
La norme de A mène à une équation faisant intervenir l'énergie E
ce qui peut être réécrit en termes d'excentricité
Donc, si l'énergie est négative (état lié) l'excentricité est inférieure à un et l'orbite est une ellipse. Si l'énergie est positive (état libre ou état de diffusion), l'excentricité est supérieure à un est l'orbite est une hyperbole. Pour finir, si l'énergie est exactement égale à zéro l'orbite est une parabole (excentricité égale à 1). Dans tous les cas A est colinéaire à l'axe de symétrie de la conique et est dirigé du centre de force vers le periapside, point de plus courte approche.
Hodographes circulaires de la quantité de mouvement
La conservation du vecteur de Runge-Lenz A et du vecteur moment cinétique L est utile pour montrer que le vecteur p se déplace sur un cercle dans un mouvement à force centrale suivant une loi en carré inverse[6],[7]. Le produit vectoriel de A et L mène à une équation pour p
En prenant L le long de l'axe z et le grand-axe suivant l'axe x on obtient l'équation :
En d'autres termes, le vecteur quantité de mouvement p décrit un cercle de rayon mk/L centré en (0, A/L). L'excentricité e correspond au cosinus de l'angle η vu à la Figure 3. Il est également intéressant d'introduire la variable . Cet hodographe circulaire permet d'illustrer la symétrie du problème de Kepler.
Constantes du mouvement et superintégrabilité
Les sept grandeurs scalaires E, A et L (étant de vecteurs, les deux dernières comptent comme trois grandeurs conservées chacune) sont reliées par deux équations A · L = 0 et A2 = m2k2 + 2 m E L2, ce qui laisse 5 constantes du mouvement indépendantes. Ceci est cohérent avec les six conditions initiales du système (la position initiale et le vecteur vitesse initiale ayant chacune trois composantes) qui détermine l'orbite de la particule puisque la date initiale ne peut pas être déterminée par une constante du mouvement. Étant donné que la norme de A (et donc l'excentricité e de l'orbite) peut être déterminée par le moment cinétique L et l'énergie E, seule la direction de A est conservée indépendamment ; en outre, comme A doit être perpendiculaire à L, il ne fournit en fait qu'une seule constante du mouvement.
Un système mécanique à d degrés de liberté peut avoir au maximum 2d - 1 constantes du mouvement indépendantes, puisqu'il y a 2d conditions initiales et que la date initiale ne peut être déterminée par une constante du mouvement[17]. Un système qui possède plus de d constantes du mouvement est appelé superintégrable et un système à 2d − 1 constantes est dit maximalement superintégrable[18].
Puisque l'équation de Hamilton-Jacobi dans un système de coordonnées ne mène qu'à d constantes du mouvement, les systèmes intégrables sont séparables dans plus d'un système de coordonnées[19] Le problème de Kepler est maximalement superintégrable puisqu'il possède 3 degrés de liberté et 5 constantes du mouvement indépendantes ; son équation de Hamilton-Jacobi est séparable à la fois en coordonnées sphériques et en coordonnées paraboliques (en)[17]. Les systèmes maximalement superintégrables mènent à des orbites fermées à une dimension dans l'espace des phases, puisque l'orbite est l'intersection des iso-surfaces de leurs constantes du mouvement dans l'espace des phases. Ce systèmes peuvent être quantifiés canoniquement en utilisant uniquement des relations de commutation comme expliqué plus bas[20].
Évolution dans des potentiels perturbés
Le vecteur de Runge-Lenz A est conservé uniquement pour une force centrale suivant loi strictement en carré inverse. Dans la plupart des problèmes pratiques tels que les mouvement de planètes, l'énergie potentielle d'interaction entre deux corps ne découle pas rigoureusement d'une telle loi et peut inclure un terme additif, appelé potentiel de perturbation h(r). Dans de tels cas le vecteur de Runge-Lenz tourne lentement dans le plan de l'orbite ce qui correspond à une lente précession du périastre de l'orbite. Par hypothèse le potentiel de perturbation h(r) est conservatif et central, ce qui implique que l'énergie totale E et le vecteur moment cinétique L sont conservés. En conséquence le mouvement reste contenu dans un plan perpendiculaire à L et la norme de A se conserve selon l'équation A2 = m2k2+2mEL2. Le potentiel de perturbation h(r) peut être n'importe quelle fonction mais doit être significativement plus petit que le potentiel découlant de la force en carré inverse entre les deux corps.
La vitesse à laquelle le vecteur de Runge-Lenz tourne donne des informations sur le potentiel de perturbation h(r). En utilisant la théorie des perturbations canoniques et les coordonnées action-angle, on obtient immédiatement[1] que A tourne à la vitesse
où T est la période orbitale et dans laquelle l'identité L dt = m r2 dθ a été utilisée pour convertir l'intégrale temporelle en intégrale angulaire (Figure 5). L'expression entre crochets, 〈h(r)〉, représente la moyenne sur une période du potentiel de perturbation, c'est-à-dire la moyenne pour la particule décrivant complètement une fois son orbite. Mathématiquement la moyenne temporelle, qui est la grandeur entre accolades, permet d'éliminer les fluctuations dans le taux de rotation.
Cette approche a été utilisée pour permettre une vérification de la théorie d'Einstein de la relativité générale, qui ajoute un terme cubique de perturbation au potentiel gravitationnel newtonien[21]
En injectant cette fonction dans l'intégrale et en utilisant l'équation
pour exprimer r en fonction de θ, la précession du périastre due à des perturbations non-newtonniennes est[21]
ce qui correspond étroitement à l'anomalie observée de la précession de Mercure[22] et des pulsars binaires[23]. Cet accord avec l'expérience est considéré comme une preuve importante de la validité de la relativité générale[24],[25].
Crochets de Poisson
Les trois composantes Li du vecteur moment cinétique L ont les crochets de Poisson[1]
où i=1,2,3 et εijs est le symbole de Levi-Civita ; l'indice de sommation s est utilisé ici pour éviter l'ambiguïté avec le paramètre k de la force vue plus haut.
Comme expliqué plus bas une autre échelle du vecteur de Runge-Lenz, D, peut être définie dans la même unité que le moment cinétique en divisant A par p0. Le crochet de Poisson de D avec le vecteur moment cinétique L peut être sous une forme similaire[26]
Le crochet de Poisson de D avec lui-même dépend du signe de E. Pour des énergies négatives (systèmes liés de trajectoire fermée, elliptique si la force suit une loi en carré inverse) le crochet de Poisson est :
tandis que pour une énergie positive (systèmes libres de trajectoire ouverte, hyperbolique si la force est en carré inverse de la distance) le crochet a le signe opposé
Les opérateurs de Casimir pour des énergies négatives sont définis par
et ont des crochets de Poisson nuls avec toutes les composantes de D et L
C2 est nul de manière évidente, puisque les deux vecteurs sont toujours perpendiculaires. Cependant l'autre invariant C1 n'est pas trivial et dépend seulement de m, k et E. Cet invariant permet d'obtenir les niveaux d'énergie de l'atome d'hydrogène en utilisant uniquement les relations de commutation de la mécanique quantique à la place de la méthode plus commune s'appuyant sur l'équation de Schrödinger.
Mécanique quantique de l'atome d'hydrogène
Les crochets de Poisson fournissent une méthode simple de quantification canonique d'un système mécanique ; la relation de commutation de deux opérateurs quantiques est égal au crochet de Poisson des variables classiques correspondantes multiplié par [27]. En menant le calcul de cette quantification et en déterminant les valeurs propres de l'opérateur de Casimir C1 du problème de Kepler, Wolfgang Pauli a obtenu spectre énergétique des atomes hydrogénoïdes (Figure 6) et donc leur spectre d'émission atomique[3] Ce résultat élégant a été obtenu avant le développement de la mécanique ondulatoire[28].
Une subtilité de l'opérateur quantique associé au vecteur de Runge-Lenz A est que les opérateurs impulsion et moment cinétique ne commutent pas ; le produit vectoriel de p et L doit alors être défini avec précaution[26]. Généralement, les opérateurs pour les composantes cartésiennes As sont définis en utilisant un produit symétrique :
à partir desquels les opérateurs d'échelle correspondant peuvent être définis
Un premier opérateur de Casimir normalisé peut de la même manière être défini
où H−1 est l'inverse de l'opérateur énergie hamiltonien et I est l'opérateur identité. En appliquant cet opérateur aux états propres des opérateurs moment cinétique total, moment cinétique azimutal et énergie, les valeurs propres du premier opérateur de Casimir C1 sont n2 - 1 ; elles sont indépendantes des nombres quantiques l et m, rendant les états d'énergies dégénérés[26]. Les états d'énergie sont donnés par
ce qui est équivalent à la formule de Rydberg pour les atomes hydrogénoïdes (Figure 6).
Conservation et symétrie
La conservation du vecteur de Runge-Lenz correspond à une symétrie subtile du système. En mécanique classique les symétries sont des opérations continues qui permettent de passer d'une orbite à une autre sans changer l'énergie du système ; en mécanique quantique les symétries sont des opérations continues qui mélangent les orbitales atomiques de même énergie c'est-à-dire qui provoquent une dégénérescence des niveaux d'énergie. Une grandeur se conservant est généralement associée à de telles symétries[1]. Par exemple toute force centrale est symétrique par le groupe de rotations SO(3) ce qui mène à la conservation du moment cinétique L. Classiquement une rotation complète du système ne modifie pas l'énergie de l'orbite ; quantiquement, les rotations combinent les harmoniques sphériques de même nombre quantique l sans changer l'énergie.
La symétrie d'une force centrale suivant une loi en carré inverse est plus élevée et plus compliquée. La symétrie particulière du problème de Kepler résulte de la conservation simultanée du moment cinétique L et du vecteur de Runge-Lenz A (tel que défini plus haut) et, en mécanique quantique, assure que les niveaux d'énergie de l'atome d'hydrogène ne dépendent pas des nombres quantiques l et m. La symétrie est cependant plus subtile car elle prend place dans un espace de dimension supérieure ; de telles symétries sont parfois appelées « symétries cachées »[29].
Classiquement la haute symétrie du problème de Kepler permet une perturbation continue des orbites conservant l'énergie mais pas le moment cinétique ; en d'autres termes, des orbites de même énergie mais de moments cinétiques différents (donc d'excentricités différentes) peuvent être transformées continûment les unes en les autres. De manière quantique cela consiste à faire une combinaison d'orbitales différant par leurs nombres quantiques l et m telles que des orbitales atomiques s (l=0) et p (l=1) par exemple. De telles combinaisons ne peuvent être faites en considérant les rotations et translations habituelles à trois dimensions, mais correspondent à une rotation dans une espace de plus grande dimension.
Pour des énergies négatives, c'est-à-dire es systèmes liés, le groupe de plus haute symétrie du système est SO(4), qui conserve la longueur dans un espace à quatre dimensions
En 1935, Vladimir Fock montre que le problème de Kepler d'un système quantique lié est équivalent au problème d'une particule libre se déplaçant sur une 3-sphère dans un espace à quatre dimensions[4]. Plus précisément Fock montre que l'équation de Schrödinger des fonctions d'onde dans l'espace des moments est la projection stéréographique des harmoniques sphériques sur la sphère.
Notes et références
Références
- H. Goldstein (en), Classical Mechanics, Addison Wesley, 1980, 2e éd., p. 102–105,421–422
- (en) V. Arnold, Mathematical Methods of Classical Mechanics, 2nd ed., New York, Springer-Verlag, 1989, 2e éd. (ISBN 978-0-387-96890-2), p. 38 (traduit du russe)
- (de) W. Pauli, « Über das Wasserstoffspektrum vom Standpunkt der neuen Quantenmechanik », dans Zeitschrift für Physik, vol. 36, 1926, p. 336–363
- (de) V. Fock, « Zur Theorie des Wasserstoffatoms », dans Zeitschrift für Physik, vol. 98, 1935, p. 145–154
- (de) V. Bargmann (en), « Zur Theorie des Wasserstoffatoms: Bemerkungen zur gleichnamigen Arbeit von V. Fock », dans Zeitschrift für Physik, vol. 99, 1936, p. 576–582
- (en) W. R. Hamilton, « The hodograph or a new method of expressing in symbolic language the Newtonian law of attraction », dans Proceedings of the Royal Irish Academy, vol. 3, 1847, p. 344ff
- (en) H. Goldstein, « Prehistory of the Runge–Lenz vector », dans American Journal of Physics (en), vol. 43, 1975, p. 735–738
(en) H. Goldstein, « More on the prehistory of the Runge–Lenz vector », dans American Journal of Physics, vol. 44, 1976, p. 1123–1124 - (en) W. R. Hamilton, « Applications of Quaternions to Some Dynamical Questions », dans Proceedings of the Royal Irish Academy, vol. 3, 1847, p. Appendix III
- (en) D. M. Fradkin, « Existence of the Dynamic Symmetries O4 and SU3 for All Classical Central Potential Problems », dans Progress of Theoretical Physics, vol. 37, 1967, p. 798–812
- (en) T. Yoshida, « Two methods of generalisation of the Laplace–Runge–Lenz vector », dans European Journal of Physics, vol. 8, 1987, p. 258–259
- (it) J. Hermann, « Metodo d'investigare l'Orbite de' Pianeti, nell'ipotesi che le forze centrali o pure le gravità … », dans Giornale de Letterati D'Italia, vol. 2, 1710, p. 447–467
J. Hermann, « Extrait d'une lettre de M. Herman à M. Bernoulli datée de Padoüe le 12. Juillet 1710 », dans Histoire de l'académie royale des sciences, vol. 1732, 1710, p. 519–521 - J. Bernoulli, « Extrait de la Réponse de M. Bernoulli à M. Herman datée de Basle le 7. Octobre 1710 », dans Histoire de l'académie royale des sciences, vol. 1732, 1710, p. 521–544
- P.-S. Laplace, Traité de mécanique celeste, 1799, « Tome I, Premiere Partie, Livre II », p. 165 et suivantes
- (en) J. W. Gibbs, Wilson EB, Vector Analysis, New York, Scribners, 1901, p. p. 135
- C. Runge, Vektoranalysis, Leipzig, Hirzel, 1919, p. Volume I
- (de) W. Lenz, « Über den Bewegungsverlauf und Quantenzustände der gestörten Keplerbewegung », dans Zeitschrift für Physik, vol. 24, 1924, p. 197–207
- L. Landau et E. Lifchitz, Physique théorique, vol. 1 : Mécanique, Moscou, MIR, 1982 [détail de l’édition]
- (en) N. W. Evans, « Superintegrability in classical mechanics », dans Physical Review A, vol. 41, 1990, p. 5666–5676
- (en) A. Sommerfeld, Atomic Structure and Spectral Lines, London, Methuen, 1923, p. 118
- (en) N. W. Evans, « Group theory of the Smorodinsky–Winternitz system », dans Journal of Mathematical Physics, vol. 32, 1991, p. 3369–3375
- (de) A. Einstein, « Erklärung der Perihelbewegung des Merkur aus der allgemeinen Relativitätstheorie », dans Sitzungsberichte der Preussischen Akademie der Wissenschaften, vol. 1915, 1915, p. 831–839
- U. Le Verrier, « Lettre de M. Le Verrier à M. Faye sur la Théorie de Mercure et sur le Mouvement du Périhélie de cette Planète », dans Comptes Rendus de l'Academie de Sciences (Paris), vol. 49, 1859, p. 379–383
- C. M. Will, General Relativity, an Einstein Century Survey, Cambridge, Cambridge University Press, 1979, SW Hawking and W Israel, eds.e éd., p. Chapter 2
- (en) A. Pais , Subtle is the Lord: The Science and the Life of Albert Einstein, Oxford University Press, 1982
- N. T. Roseveare, Mercury's Perihelion from Le Verrier to Einstein, Oxford University Press, 1982
- (en) A. Bohm, Quantum Mechanics: Foundations and Applications, Springer Verlag, 1986, 2e éd., p. 208–222
- (en) P. Dirac, Principles of Quantum Mechanics, 4th revised edition, Oxford University Press, 1958
- (de) E. Schrödinger, « Quantisierung als Eigenwertproblem », dans Annalen der Physik, vol. 384, 1926, p. 361–376
- (en) G. E. Prince et C. J. Eliezer, « On the Lie symmetries of the classical Kepler problem », dans Journal of Physics A: Mathematical and General, vol. 14, 1981, p. 587–596
Notes
- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Laplace-Runge-Lenz vector » (voir la liste des auteurs)
Voir aussi
Article connexe
Lien externe
Démonstration des lois de Kepler et propriétés d'une ellipse, cours de mécanique par Bernard Gisin (site personnel)
Wikimedia Foundation. 2010.