Principe de moindre action
331228
31405371
2008-07-08T12:56:05Z
LyricV
210569
/* Résumé anthropomorphique : référence plus proche de ce qui est écrit*/
[[Image:Joseph Louis Lagrange.jpg|thumb|right|[[Joseph-Louis Lagrange]] est le premier a définir mathématiquement le principe de moindre action.]] {{Sciences}}
En [[physique]], le '''principe de moindre action''' est l'hypothèse selon laquelle la dynamique d'une quantité physique (la position, la vitesse et l'accélération d'une particule, ou les valeurs d'un [[champ (physique)|champ]] en tout point de l'espace, et leurs variations) peut se déduire à partir d'une unique grandeur appelée [[action (physique)|action]], dépendant de la quantité physique considérée, et en supposant que les valeurs dynamiques permettent à l'action d'avoir une valeur minimale entre deux instants donnés.
La plupart des équations fondamentales de la physique peuvent être formulées à partir du principe de moindre action. C'est notamment le cas en [[Mécanique newtonienne|mécanique classique]], en [[électromagnétisme]], en [[relativité générale]] et en [[théorie quantique des champs]].
== Formulation historique ==
Dans ''Principe de la moindre quantité d'action pour la mécanique'' ([[1744]]), [[Pierre Louis Moreau de Maupertuis|Maupertuis]] définit l'action comme suit :
{{Début citation}}L'Action est proportionnelle au produit de la masse par la vitesse et par l'espace. Maintenant, voici ce principe, si sage, si digne de l'Être suprême : lorsqu'il arrive quelque changement dans la Nature, la quantité d'Action employée pour ce changement est toujours la plus petite qu'il soit possible.{{Fin citation}}
[[Pierre de Fermat|Fermat]], [[Samuel König|König]] et [[Gottfried Wilhelm von Leibniz|Leibniz]] avaient avancé le même principe sous le nom de « principe d’économie naturelle »<ref>Voir par exemple http://serge.mehl.free.fr/chrono/Maupertuis.html</ref>; lequel deviendra le principe de conservation de l’énergie avec les travaux d’[[Leonhard Euler|Euler]], de [[Joseph-Louis Lagrange|Lagrange]], de [[Charles Gustave Jacob Jacobi|Jacobi]] et de [[Hermann Ludwig von Helmholtz|Helmholtz]].
==Vulgarisation et interprétations==
===Résumé anthropomorphique===
Le principe de moindre action dit qu'en mécanique un corps prend la direction qui lui permet de dépenser le moins d'énergie dans l'immédiat (ou d'acquérir le plus d'énergie dans l'immédiat), en tenant compte qu'il doit y avoir continuité du mouvement (positions et vitesses) s'il y a continuité des conditions physiques<ref>On peut voir à ce sujet le chapitre 3, p117-134, et les exercices ''E3.4'' et ''E3.5'', p144-145, du livre ''Mécanique : de la formulation lagrangienne au chaos hamiltonien'' de Claude Gignoux et Bernard Silvestre-Brac, éditeur EDP-Sciences, 2002, ISBN 2868835848</ref>. <!--<ref>Principe de Maupertuis dans : §44 de ''Mécanique'',Landau-Lifchitz, traducteur Claude Ligny, 1960, Éditeur MIR ou p173 à 175 de ''Calcul des variations'' par Mr Pierre Bérest, 1997, Éditeur ellipse, ISBN 2 7298 9704 5</ref>-->
Il est à remarquer qu'en reliant deux points, la trajectoire prise par le corps n'est pas toujours celle qui lui fait dépenser ''globalement'' le moins d'énergie car c'est la dépense ''immédiate'' d'énergie qui est minimisée (comme si le corps ne percevait que les conditions de son environnement immédiat) et si le chemin parcouru est long, un chemin plus court avec une dépense d'énergie immédiate plus élevée peut permettre une dépense globale inférieure. Une analogie avec la consommation en carburant d'une voiture peut être faite.
Note : dans ce « résumé », "''énergie''" signifie ''[[énergie cinétique]]'', et une "''dépense d'énergie''" signifie que de l'énergie cinétique se transforme en [[énergie potentielle]].
===Méthode et interprétation en physique classique===
L'action se présente comme la sommation, le long du parcours du système, de la différence entre l'[[énergie cinétique]] et l'[[énergie potentielle mécanique|énergie potentielle]]. La minimisation se fait par une méthode variationnelle : à points extremum fixés, temps de trajet fixé, et trajet variable, on cherche les conditions imposées au trajet pour qu'il minimise l'action.
On peut interpréter cela comme équivalent aux deux conditions suivantes:
* La trajectoire que suit un corps est celle qui permet la transformation instantanée de l'[[énergie cinétique]] en [[énergie potentielle mécanique|énergie potentielle]] la plus petite possible (donc aussi la plus lente sur la trajectoire), ou la transformation immédiate dans le sens inverse la plus grande possible (donc la plus rapide possible sur la trajectoire).
* La transformation (et donc la trajectoire) est déterminée par les conditions initiales (position et vitesse) et les conditions de l'environnement physique : il doit y avoir continuité de la trajectoire s'il y a continuité du milieu physique.
Il y a parfois un échange cyclique entre ces deux énergies (balancier sans frottement, satellite à orbite elliptique,...) ou une stabilisation provisoire (bille immobile ou posée au fond d'un trou, satellite à orbite circulaire,...).
La chute libre d'un corps est l'exemple type de la transformation de l'énergie potentielle (gravitationnelle) en énergie cinétique. Le ralentissement et l'arrêt (avant sa chute) d'un corps lancé verticalement est un exemple de la transformation inverse.
Les frottements imposent une transformation plus compliquée car ils engendrent de la [[Énergie thermique|chaleur]], qui est l'énergie cinétique des molécules des matériaux, mais en négligeant cette forme d'énergie, on peut utiliser le Principe de moindre action en considérant que de l'énergie cinétique se perd (sort du système étudié).
===Un problème métaphysique ?===
Le principe de moindre action utilise l'hypothèse de deux points fixes sur le parcours du mobile : un point de départ, mais aussi un point d'arrivée. Cela a souvent été critiqué comme étant l'utilisation dans le raisonnement d'une « cause finale », ce qui est contraire à la causalité qui suit la flèche du temps en physique.
En fait, si le point de départ est doté de conditions initiales (coordonnées et vitesse), le point d'arrivée n'a pas de coordonnées précises ni de vitesse imposée : il existe, c'est tout. L'existence du point final dans le raisonnement permet d'émettre l'hypothèse de l'existence d'un ''trajet'' à partir de l'état initial et de déterminer ses conditions (équations d'Euler-Lagrange), mais n'impose aucune autre condition en dehors de la continuité indiquée plus haut (ce travail peut même montrer que seul un trajet de longueur nulle est possible dans les cas de stabilité du mobile).
===Un principe démontré===
Avant Lagrange ce principe se concevait à partir de considérations métaphysiques, indépendamment de tout autre principe physique.
Lagrange, en 1756, fut celui qui donna au Principe de moindre action son expression mathématique efficace qui est toujours d'actualité. Il fut aussi celui qui développa la [[mécanique analytique]] et démontra, dans son ouvrage de 1788, ce principe à partir de la conservation de l'énergie et du ''principe des vitesses virtuelles'' (nommé aussi [[principe de d'Alembert]]). Le ''principe des vitesses virtuelles'' est le [[lois du mouvement de Newton|principe fondamental de la dynamique de Newton]] exprimé dans le langage de l'[[analyse (mathématiques)|analyse mathématique]], alors balbutiante au regard de ses développements ultérieurs.
Cette démonstration met un point final aux interrogations métaphysiques sur le principe de moindre action : le principe est équivalent à un principe physique de Newton, non sujet aux critiques métaphysiques. Mais cette conclusion est souvent oubliée, et dans ce cas l'artifice mathématique de la « cause finale » est discuté comme étant problématique.
===De l'action classique à l'action relativiste===
Suivant le système étudié, et le cadre théorique dans lequel on le considère, l'expression mathématique du principe de moindre action change légèrement de forme.
C'est un des rares principes ayant survécu aux multiples mutations de la physique, mais il a rarement été à l'origine d'une découverte : il est plutôt utilisé pour reformuler ou redémontrer des lois trouvées par d'autres biais. Sa plus grande contribution a sans doute été de mettre [[William Rowan Hamilton|W. R. Hamilton]] sur la voie de ses travaux théoriques (voir: [[Mécanique hamiltonienne]]).
En physique relativiste, les équations d'Euler-Lagrange restent inchangées, mais le lagrangien n'est plus égal à la différence entre l'énergie cinétique et l'énergie potentielle. En fait, à partir de la relativité il est apparu que le principe de moindre action se base sur l'existence d'une trajectoire continue, paramétrée par le temps, qui minimise une fonction ou la différence entre des fonctions du système étudié, déterminées à partir de principes généraux, tels que par exemples :
* Comme la trajectoire dans l'espace-temps ne dépend pas du repère d'où on l'observe, l'action qui la détermine, ainsi que les fonctions qui composent l'action, sont invariantes par changement de repère.
* L'indépendance de corps implique l'additivité de leurs actions et de leurs lagrangiens, pour que les trajectoires puissent être déterminées séparément en appliquant la méthode variationnelle.
Il se trouve qu'en physique classique, ces fonctions du système sont les énergies cinétiques et potentielles, ce n'est plus le cas en relativité.
En physique relativiste, et en l'absence de champ électromagnétique, on montre que la fonction du corps qui est minimisée dans le principe est particulièrement simple : il s'agit de <math>-mc\tau</math>, où <math>\tau</math> est « [[temps propre]] » du trajet, qui est à la fois le temps s'écoulant dans le référentiel du corps au cours du trajet et la longueur de la trajectoire mesurée par la métrique de l'espace : cela revient à maximiser le « [[temps propre]] », du fait du signe <math>-</math> et de la constance de la masse <math>m</math> et de la vitesse de la lumière <math>c</math>.<br />
Un champ électromagnétique amène des différences de parcours entre les corps, suivant leurs charges et leurs répartitions.<br />
Et comme en physique classique, toutes les équations peuvent être obtenues sans le principe de moindre action.
===Sa formulation en mécanique quantique===
Dans le but de trouver une formulation plus simple de l'[[électrodynamique quantique]], vers 1940, [[Richard Feynman]] chercha une formulation du principe de moindre action en [[mécanique quantique]]. La solution lui vint d'une idée que [[Paul Dirac]] avait exprimée dans un article<ref>Florence Martin-Robine ; ''Histoire du principe de moindre action'', Vuibert (2006), p 206-209; et R.Feynmann, Conférence Nobel, dans ''La nature de la physique'', Seuil (1980).</ref>.<br />
Le principe a ainsi permis une reformulation de cette branche de la physique sous forme d'[[intégrale de chemin]] qui s'est révélé, en effet, plus simple que la formulation hamiltonnienne pour l'[[électrodynamique quantique]].<br />
Cette formulation a donné lieu à des interprétations telle que « la particule teste tous les chemins possibles avec des probabilités différentes ».<ref>Florence Martin-Robine ; ''Histoire du principe de moindre action'', Vuibert (2006), p 209.</ref>
Comme on peut s'y attendre, la formulation quantique permet de retrouver, à la limite classique, la formulation habituelle et le chemin qui rend extrémale l'action classique est un col de l'intégrale : seul celui ci contribue de manière significative dans l'intégrale.<br /><br />
{{Loupe|Intégrale de chemin}}
===Court historique<ref>''Histoire du principe de moindre action'' par Florence Martin-Robine, chez Vuibert, 2006</ref>===
[[Image:Hilbert.jpg|thumb|right|En 1916, [[David Hilbert]] a redémontré les équations de la gravitation de la [[relativité générale]] à l'aide du principe de moindre action.<br>(Photo prise en 1912.)]]
L'idée que la trajectoire minimise une durée ou une longueur est d'abord née chez [[Pierre de Fermat]] vers 1655 pendant son étude de l'optique (voir [[Principe de Fermat]]). Même si elle a interessé [[Gottfried Wilhelm von Leibniz|Leibnitz]] et [[Isaac Newton|Newton]], c'est [[Pierre Louis Moreau de Maupertuis|Maupertuis]], vers 1740, qui fera progresser la formulation verbale et mathématique d'un « principe de moindre action » pour la mécanique. [[Leonhard Euler|Euler]], en développant l'[[analyse (mathématiques)|analyse mathématique]], commença à reformuler ce principe, mais c'est [[Joseph-Louis Lagrange|Lagrange]] qui lui donnera sa méthode et sa forme définitive en 1755, pour ensuite l'inclure comme une simple conséquence de sa [[mécanique analytique]].
En 1827, [[William Rowan Hamilton|Hamilton]], en cherchant à appliquer ce principe à l'optique, développa une nouvelle approche basée sur l'étude de l'énergie par la méthode analytique : la [[mécanique hamiltonienne]], que [[Charles Gustave Jacob Jacobi|Jacobi]] peaufinera vers 1840.
Depuis sa formulation, ce principe a guidé de nombreux scientifiques dans leurs recherches, notamment [[Louis de Broglie|de Broglie]] vers 1920 dans son travail sur la [[théorie des quanta]]. En 1916, [[David Hilbert|Hilbert]] a redémontré les équations de la gravitation de la [[relativité générale]] à l'aide du principe, et [[Richard Feynman]], en 1942, a proposé une nouvelle formulation du principe dans sa thèse de doctorat intitulée ''Le Principe de moindre action en mécanique quantique'', permettant une réécriture de la [[mécanique quantique]].
==En mécanique classique==
===Une alternative de présentation===
* Soit on expose l'action et le lagrangien habituels de la physique classique (non relativiste), puis on détermine les équations d'Euler-Lagrange.
* Soit on définit abstraitement l'action et le lagrangien (''à la manière de [[Lev Landau|Landau]] et [[Evguéni Lifchitz|Lifchitz]]''<ref>L. Landau et E. Lifchitz, trad. Claude Ligny,''Physique théorique'', tome 1 « Mécanique », éditions MIR, Moscou, 1982, {{ISBN|5-03-000198-0|5-03-000197-2}}</ref>), et on détermine leurs formes et leurs propriétés qu'imposent les principes de la physique, ainsi que les équations d'Euler-Lagrange.
Dans cet article, seule la première présentation sera donnée.
===Définition du lagrangien et de l'action classiques===
Considérons pour simplifier un point matériel décrit, '''dans un repère galiléen''', par ''un seul [[degré de liberté]]'', noté <math>q(t)</math> à l'instant <math>t</math> <ref>La généralisation à un nombre quelconque de degrés de liberté ne pose pas de problème de principe.</ref>.
Le lagrangien est la différence entre l'énergie cinétique (du point matériel) et l'énergie potentielle (due à l'environnement physique) :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|<math> L(q,\dot{q},t) \ = \ \frac{1}{2} \, m \, \dot{q}^2 \ - \ V(q,\dot{q},t)\ </math>
}}
où <math>\frac{1}{2} \, m \, \dot{q}^2</math> est l'énergie cinétique du système et <math> \ V(q,\dot{q},t)</math> est l'énergie potentielle, qui en général ne dépend pas de <math>\dot{q}</math>.
L'action de la trajectoire, étant la somme totale de la différence entre l'énergie cinétique et l'énergie potentielle durant la trajectoire, est définie par:
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>S \ = \ \int_{t_i}^{t_f} L(q(t),\dot{q}(t),t)\ dt</math>
}}
où <math>t_i</math> et <math>t_f</math> désignent respectivement l'instant initial et l'instant final.
=== Les équations d'Euler-Lagrange ===
* Dire que le chemin « minimise localement l'action » signifie que pour tout autre chemin ayant les mêmes conditions initiales et finales, et suffisamment proche du chemin minimisant, la valeur de l'action est plus grande.
* Avec certaines conditions initiales et/ou finales, le chemin minimisant localement l'action peut ne pas exister, mais s'il existe, il est unique (à cause des conditions initiales et de la continuité du mouvement).
Le lagrangien n'est pas défini de manière unique : l'ajout au lagrangien d'une fonction <math>f(q,\dot{q},t) = \frac{dF}{dt}</math> ajoute à l'action une fonction <math>\ F(t_{final}) - F(t_{initial})</math> qui ne dépend que des extrémités et qui s'annule quand on fait varier l'action par rapport au chemin.
Le chemin <math>q(t)</math> effectivement suivi par le point matériel entre les instants <math>t_i</math> et <math>t_f</math> fixés est un ''extremum'' de l'action (car il lui fait atteindre sa valeur minimale), donc en faisant une variation du chemin, on a :
<center>
{| border=1 cellpadding=10 style="border-collapse:collapse"
|-----
| |<math> \delta S[q] \ = \ 0 </math>
|}</center>
Les équations (d'Euler-Lagrange) que l'on en déduit sont:
{{Bloc emphase|fond=#fff8ff;border:4px double #888| <math>\frac{d~~ }{dt} \frac{\partial L}{\partial \dot{q}}(q, \dot{q},t) \ - \ \frac{\partial L}{\partial q}(q, \dot{q},t) \ = \ 0</math>
}}
{{Boîte déroulante|titre=Démonstration des équations d'Euler-Lagrange|contenu=
<center>
Soit <math>q_0(t)</math> le chemin réellement suivi entre les instants <math>t_i</math> et <math>t_f</math>, vérifiant les conditions aux limites :
<math> q_0(t_i) \ = \ q_i \quad ; \quad q_0(t_f) \ = \ q_f</math>
Soit alors un chemin quelconque <math>q(t)</math> vérifiant les ''mêmes'' conditions aux limites :
<math> q(t_i) \ = \ q_i \quad ; \quad q(t_f) \ = \ q_f</math>
On peut toujours décomposer ce chemin quelconque <math>q(t)</math> sous la forme :
<math> q(t) \ = \ q_0(t) \ + \ \epsilon \, f(t) </math>
où <math>\epsilon</math> est un nombre et <math>f(t)</math> une fonction continue qui doit vérifier les conditions aux limites :
<math> f(t_i) \ = \ f(t_f) \ = \ 0</math>
Le principe de moindre action nous dit alors que, lorsque <math> \epsilon </math> non nul tend vers zéro :
<math> \delta S[q_0] \ = \ S[q_0 + \epsilon f] \ - \ S[q_0] \ = \ 0 \ + \ O(\epsilon^2)</math>
Calculons l'action <math>S[q_0 + \epsilon f]</math> au premier ordre :
<math>S[q_0 + \epsilon f] \ = \ \int_{t_i}^{t_f} L(q_0(t) + \epsilon f(t), \dot{q}_0(t) + \epsilon \dot{f}(t) \, ) \ dt</math>
Un développement limité au premier ordre du lagrangien donne, en omettant la dépendance explicite en temps pour simplifier l'écriture :
<math>L(q_0 + \epsilon f, \dot{q}_0 + \epsilon \dot{f} \, ) \ = \ L(q_0, \dot{q}_0\, ) \ + \epsilon \ \left[ \, \frac{\partial L(q_0, \dot{q}_0)}{\partial q} \ f \ + \ \frac{\partial L(q_0, \dot{q}_0)}{\partial \dot{q}} \ \dot{f} \, \right] \ + \ O(\epsilon^2)</math>
La variation de l'action est donc égale à :
<math>\delta S[q_0] \ = \ \epsilon \ \int_{t_i}^{t_f} \left[ \, \frac{\partial L(q_0, \dot{q}_0)}{\partial q} \ f(t) \ + \ \frac{\partial L(q_0, \dot{q}_0)}{\partial \dot{q}} \ \dot{f}(t) \, \right] \ dt \ + \ O(\epsilon^2)</math>
On peut réécrire le second terme en utilisant une intégration par partie :
<math>\int_{t_i}^{t_f} \frac{\partial L(q_0, \dot{q}_0)}{\partial \dot{q}} \ \dot{f}(t) \ dt \ = \ \left[ \frac{\partial L(q_0, \dot{q}_0)}{\partial \dot{q}} \ f(t)\right]_{t_i}^{t_f} \ - \ \int_{t_i}^{t_f} \frac{d~~ }{dt} \left(\frac{\partial L(q_0, \dot{q}_0)}{\partial \dot{q}}\right) \ f(t) \ dt \ </math>
Le terme entre crochets est nul en raison des conditions aux limites imposées à la fonction <math>f(t)</math>. On en déduit que la variation de l'action s'écrit au premier ordre :
<math>\delta S[q_0] \ = \ \epsilon \ \int_{t_i}^{t_f} f(t) \ \left[ \, \frac{\partial L(q_0, \dot{q}_0)}{\partial q} \ - \ \frac{d~~ }{dt} \left(\frac{\partial L(q_0, \dot{q}_0)}{\partial \dot{q}}\right) \, \right] \ dt \ + \ O(\epsilon^2)</math>
Cette variation doit être ''nulle'' d'après le principe de moindre action, ceci pour tout <math>\epsilon \ne 0</math> et toute fonction <math>f(t) \ne 0</math> vérifiant les conditions aux limites, donc l'''intégrand doit être nul'', et on obtient les équations d'Euler-Lagrange :
<math>\frac{d~~ }{dt} \frac{\partial L}{\partial \dot{q}}(q_0, \dot{q}_0,t) \ - \ \frac{\partial L}{\partial q}(q_0, \dot{q}_0,t) \ = \ 0</math>
</center>
}}
===L'équation fondamentale de la dynamique newtonienne===
Supposons que <math>\ V(q,\dot{q},t) = V(q,t) </math> .
Avec l'expression du lagrangien classique, on obtient: <math>\frac{d~~ }{dt} \frac{\partial L}{\partial \dot{q}}(q, \dot{q},t) \ = \ m.\ddot{q} </math>
Les équations d'Euler-lagrange donnent:
{| align="center" border="0"
|<math> m.\ddot{q} = - \frac{\partial V}{\partial q}(q_0,t) </math>
|}
Ce qui exprime les [[lois du mouvement de Newton]] avec:
{| align="center" border="0"
|(Somme des forces extérieures) = <math> - \frac{\partial V}{\partial q}(q_0,t) </math>
|}
Si toutes les forces en jeu dérivent d'un potentiel, le principe de moindre action peut être considéré comme une réécriture des lois du mouvement de Newton.
===Cas où il y a des forces de [[frottement]]s===
On suppose que le système évolue dans un milieu homogène (mêmes propriétés à tous les endroits) et isotrope (mêmes propriétés dans toutes les directions), et dont la viscosité engendre des frottements :
* S'il n'y a qu'un seul degré de liberté <math> \ q </math> du système :
:En première approximation, les forces de frottements peuvent se modéliser par <math> \ f = -\ a.\dot{q} </math> où <math> \ a </math> est une constante positive dépendant des caractéristiques géométriques du corps et de la viscosité du milieu (air, eau, ...), le signe <math>\ - </math> indiquant que les frottements sont orientés dans le sens inverse du mouvement.
:On peut donc utiliser le potentiel <math> \ V_{\ fr} = \frac{a}{2}.(\dot{q})^2 </math> dans le lagrangien.
* S'il y a plusieurs degrés de liberté, le lagrangien s'écrit <math> L(q,\dot{q},t) \ = \frac{1}{2} \Sigma_i \, m_i \, \dot{q}_{i}^2 \ - \ V(q,\dot{q},t)\ </math>
:Et la force de frottement dans la direction de la j-ième coordonnée est <math> \ f_{j} = -\Sigma_i \ a_{i,j}.\dot{q}_i </math>
:Ce qui est une égalité scalaire, où les <math> \ a_{i,j} </math> sont des constantes, ainsi que dans le cas précédent.
:Il n'y a de potentiel que si <math> \ a_{i,j} = a_{j,i} </math> , ce qui est toujours vrai dans un fluide homogène et isotrope.
:Dans le cas où il y a un potentiel (appelé ''fonction de dissipation''), il a la forme d'une fonction quadratique : <math> \ F = \frac{1}{2} \Sigma_{ij} \ a_{i,j}.\dot{q}_i \dot{q}_j </math>
:On vérifie facilement que : <math> \ f_{j} = - \frac{\partial F}{\partial \dot{q}_j} </math>
:Dans tous les cas, pour préserver l'utilisation du lagrangien, et pour respecter l'équation fondamentale de la dynamique, on écrit :
:<math>\frac{d~~ }{dt} \frac{\partial L}{\partial \dot{q}_j}(q, \dot{q},t) \ - \ \frac{\partial L}{\partial q_j}(q, \dot{q},t) \ =\ f_{j} = -\Sigma_i \ a_{i,j}.\dot{q}_i </math>
Si la fonction de dissipation <math> \ F </math> est utilisable, on démontre que <math> \frac{dE}{dt} = -2F </math> (voir plus bas pour la définition de l'énergie ''E''), ainsi cette fonction quantifie la dissipation d'énergie du système au cours du temps.
Pour une vue plus exhaustive sur les frottements, consulter le [[b:Tribologie|wikilivre ''Tribologie'']]
===Particule chargée dans un [[champ électromagnétique]]===
Ici <math>\; e</math> est la charge de la particule.
<math> L = L(\vec{q},\dot{\vec{q}},t) = \frac{1}{2}m \ \dot{\vec{q}}^2 \ - \ V(q,\dot{\vec{q}},t)\ </math>
avec : <math> \ V(q,\dot{\vec{q}},t)= e\phi(\vec{q},t) - e \dot{\vec{q}}. \vec{A}(\vec{q},t)</math> ; la donnée <math> \left(\phi(\vec{q},t), \vec{A}(\vec{q},t)\right)</math> est appelée « potentiel électromagnétique ».
Donc :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> L = {1 \over 2} m \dot{\vec{q}}^2 - e\phi + e. \dot{\vec{q}} \cdot \vec{A}</math>
}}
En posant <math> \vec{E} = - \frac{\partial \vec{A}}{\partial t} - \vec{grad} \phi </math> le champ électrique et <math> \vec{B} = \vec{rot} \vec{A} </math> le champ magnétique.
Les équations d'Euler-Lagrange donnent :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> m \ddot{\vec{q}~}= e \left(\vec{E} + \dot{\vec{q}~} \wedge \vec{B}\right)</math>
}}
Où <math> \vec{F} = e \left(\vec{E} + \dot{\vec{q}} \wedge \vec{B}\right)</math> est appelé [[force de Lorentz]].
Historiquement, la force de Lorentz a été trouvée avant l'idée du potentiel électromagnétique.
===L'impulsion===
L'[[Moment linéaire|impulsion]] est la variable conjuguée de la vitesse dans la [[transformation de Legendre|transformée de Legendre]] du lagrangien.
Elle est définie par : <math> p = \frac{\partial L}{\partial \dot{q}}(q, \dot{q},t) </math>
* Si <math> \ V = V(q,t) </math> alors <math> p = m \dot{q} </math>
:Donc <math> \dot{q} = \frac{1}{m} p </math> , d'où : <math> L = L(q,p,t) = \frac{1}{2m} \ p^2 \ - \ V(q,t)\ </math>
* Si <math> \ V = V(q,\dot{q},t) </math> est ''affine'' par rapport à <math> \dot{q} </math> (sinon il s'agit d'une force de frottement) alors <math> p = m \dot{q} - \frac{\partial V}{\partial \dot{q}}</math>
::où <math>\frac{\partial V}{\partial \dot{q}}</math> est indépendant de <math> \dot{q} </math>.
:L'impulsion est utilisable comme variable :
:En remarquant que: <math> \dot{q} = \frac{1}{m} \left( p + \frac{\partial V}{\partial \dot{q}}\right) </math>, on a <math> L = L(q,p,t) = \frac{1}{2m} \ \left( p + \frac{\partial V}{\partial \dot{q}}\right)^2 \ - \ V(q,p,t)\ </math>
:En développant le carré, on obtient:
::<math> L = L(q,p,t) = \frac{1}{2m} \ p^2 \ - \ V'(q,t)\ </math>
Dans le lagrangien <math> \ L(q,p,t) </math>, le potentiel est indépendant de <math> \ p </math>.
Avec l'impulsion comme variable, les équations d'Euler-Lagrange ne changent pas de forme : <math> \ p </math> y prend la place de <math> \dot{q} </math>.
===Invariances et constantes du mouvement===
Le [[Théorème de Noether (physique)|théorème de Noether]] montre qu'une invariance du lagrangien par une transformation impose une grandeur invariante du système. On peut utiliser l'énergie du système pour arriver aux mêmes invariants (mis à part le premier cas).
====Indépendance par rapport au temps ====
Si les forces en présence sont indépendantes du temps ou si le système est fermé, alors le lagrangien ne dépend pas explicitement du temps : <math> \frac{\partial L}{\partial t } = 0 </math> ou encore <math> L = L(\vec{q}, \dot{\vec{q}~}) </math>
En dérivant le lagrangien par rapport au temps, on démontre que l'énergie totale
{| align="center" border="0"
|<math> E = \frac{\partial L}{\partial \dot{\vec{q}~}} .\dot{\vec{q}~} - L = \frac{1}{2} \, m \, \dot{\vec{q}~}^2 + V\left(\vec{q}\right) </math>
|}
{| align="center" border="0"
|est constante par rapport au temps.
|}
:C’est-à-dire : <math>\frac{dE }{dt} = 0 </math>
Ceci peut se démontrer directement à partir de l'équation fondamentale de la dynamique.
L'énergie du système peut aussi être définie comme la [[Transformation de Legendre|transformée de Legendre]] du lagrangien.
Dans le cas particulier où il s'agit d'un potentiel électromagnétique, et en utilisant les notations vues ci-dessus pour ce cas, on montre que : <math>E = \frac{1}{2m}.\left(\vec{p} - e.\vec{A}\right)^2 + e.\phi </math>
====Translation dans l'espace====
Pour un système fermé, et du fait de l'homogénéité (mêmes propriétés à des endroits différents) de l'espace, une translation du système de vecteur constant n'en change pas les propriétés, et donc ne change pas le lagrangien.
On en tire la conclusion de l'invariance de l'[[Moment linéaire|impulsion]] du système. Cette impulsion est la somme des impulsions des éléments du système.
====Rotation dans l'espace====
Pour un système fermé, et du fait de l'isotropie de l'espace (mêmes propriétés dans des directions différentes) de l'espace, une rotation fixée du système n'en change pas les propriétés, et donc ne change pas le lagrangien.
On en tire la conclusion de l'invariance du [[Moment angulaire|moment cinétique]] du système. Ce moment cinétique est la somme des moments cinétiques des éléments du système.
====Changement de repère galiléen====
Un changement de repère galiléen ne change pas les propriétés du système, par contre il permet de montrer qu'il existe un centre de masse (ou centre d'inertie) du système : particule virtuelle dont la masse est la masse totale du système, et animée d'un mouvement rectiligne uniforme en cas d'absence d'interaction avec l'extérieur.
Dans ce dernier cas, en posant
* <math> \mu \ </math> = ''masse totale''
* <math> \ V </math> = ''vitesse du centre de masse''
* ''énergie interne'' : <math> \ E_{int} = \Sigma_{i} \ m_{i} \ v_{i}^2 + U </math>
* ''lagrangien interne'' : <math> \ L_{int} = \Sigma_{i} \ m_{i} \ v_{i}^2 - U </math>
:avec <math> \ v_{i} </math>= vitesse relative du ième corps par rapport au centre de masse
:et <math> \ U </math> = énergie potentielle d'interactions entre les corps du système.
Le lagrangien du système peut s'écrire : <math> \ L = \frac{\mu \ V^2}{2} + \ L_{int} </math>
et l'énergie <math> \ E = \frac{\mu \ V^2}{2} + \ E_{int} </math>
Le référentiel à privilégier pour faciliter les calculs est le référentiel du centre de masse dont l'impulsion et le moment cinétique sont ceux vus ci-dessus.
====Changement de jauge====
La [[jauge de Lorenz|jauge]] est l'ensemble des unités de mesures utilisées pour mesurer les différentes caractéristiques du système. Un changement de jauge qui ne change le lagrangien que d'un facteur multiplicatif constant permet de montrer facilement certaines propriétés quand l'énergie potentielle est une fonction homogène des coordonnées (ce qui est souvent le cas).
Plus précisément : supposons que <math> \ U(a.q_1,a.q_2,....,a.q_n) = a^k.U(q_1,q_2,...,q_n) </math>
Un changement global de jauge est un changement de mesures : <math> \ q \to a.q</math> et <math> \ t \to b.t </math>
Ce qui signifie que : <math> \frac{q}{q'} = a </math> et <math> \frac{t}{t'} = b </math> où <math>\ q </math>, <math>\ q'</math> d'une part, et <math>\ t</math> et <math>\ t'</math> de l'autre, sont des mesures différentes des mêmes distances physiques ou temporelles.
Alors <math> v = \frac{dq}{dt} </math> subit le changement de mesure <math>\ v \to \frac{a}{b}.v </math> , donc l'énergie cinétique subit un changement de mesures de facteur <math> \frac{a^2}{b^2} </math>
Pour que le lagrangien ne soit que multiplié par un nombre constant, il faut que : <math> \frac{a^2}{b^2} = a^k </math> , c'est-à-dire : <math> \ b = a^{1-\frac{k}{2}} </math>
* Par exemple avec le [[Loi universelle de la gravitation|potentiel gravitationnel newtonien]], on a : <math> \ k = -1 </math>, donc <math> \ b = a^{\frac{3}{2}} </math> , c'est-à-dire : <math> \ q^3 </math> proportionnel à <math>\ t^2 </math> , ce qui correspond à la [[Lois de Kepler|3{{e}} loi de Kepler]].
==En [[relativité restreinte]]==
Le lecteur doit prendre garde que dans cette partie, et celle concernant la relativité générale, on n'étudie que le potentiel du champ électromagnétique, et la lettre <math>\ V</math> désigne une vitesse, ainsi qu'indiqué ci-dessous.
===Avec ou sans quadri-écriture===
* En relativité restreinte, les corps évoluent dans l'espace-temps de Minkowski où chaque référentiel galiléen a ses coordonnées d'espace <math>\ x_1;x_2;x_3</math> et sa coordonnée de temps <math>\ x_0 = c.t</math>, subissant toutes une modification en cas de changement de référentiel galiléen. Il n'y a donc plus de temps absolu, pourtant le temps d'un référentiel quelconque, ''galiléen ou non'', permet toujours de paramétrer l'évolution d'un système physique.
:En choisissant de repérer le système dans un référentiel galiléen quelconque, donc avec les coordonnées <math>\ (x_0;x_1;x_2;x_3)</math>, on peut choisir un temps <math>t_0</math> d'un autre référentiel quelconque, ''galiléen ou non'', pour paramétrer son évolution.
:Le lagrangien <math>\ L = L(x_0;x_1;x_2;x_3;V_0;V_1;V_2;V_3) = L(x_i;V_i)</math> exprimé à l'aide des coordonnées et de la vitesse peut donc s'écrire <math>\ L(x_i;V_i) = L(x_i(t_0);V_i(t_0))</math> , avec <math>V_i = \frac{dx_i}{dt_0}</math>.
* Si on choisit <math>t_0 = \frac{x_0}{c}~~</math> le temps du référentiel des coordonnées, le lagrangien et les équations qui en sont tirées donnent, à l'approximation aux petites vitesses devant <math>c~</math>, le lagrangien et les propriétés de la mécanique classique. On dira alors travailler sans la quadri-écriture car seules les coordonnées spatiales apparaissent en général.
* Si on choisit <math>t_0 =~~</math> le temps propre, avec <math>dt_0 = \frac{ds}{c}~~</math>, on obtient des résultats équivalents mais dont l'écriture est jugée plus élégante et s'approche de celle de la relativité générale. On remarquera qu'un repère propre n'est galiléen que si le corps est libre. Avec le temps propre, on dira travailler en quadri-écriture car les quatre coordonnées du référentiel apparaissent dans les calculs.
* Si on choisit <math>t_0 </math> un temps autre quelconque, on peut travailler plus facilement avec les dérivées partielles qu'en utilisant les deux autres temps précédents; mais les résultats, bien qu'équivalents, ont une écriture moins maniable et moins élégante. Dans ce cas, on dira aussi travailler en quadri-écriture, pour la même raison.
; Avec la quadri-écriture
Par commodité, nous adopterons la [[convention de sommation d'Einstein]] dans l'espace de Minkowski : pour deux quadri-vecteurs <math>\ (V_0;V_1;V_2;V_3)\ </math> et <math>\ (U_0;U_1;U_2;U_3)</math>, on définit le produit scalaire <math>\ V^iU_i</math> par <math> V^iU_i = \Sigma_{i=0}^3V^i.U_i = V_0.U_0-V_1.U_1-V_2.U_2-V_3.U_3 = V_iU^i </math> , avec <math>\ V^0 = V_0</math> et pour i=1;2;3 <math>\ V^i = -V_i</math>
On a alors : <math>\ (V_0)^2 - (V_1)^2 - (V_2)^2 - (V_3)^2 = V^iV_i</math>
On montre que : <math> \frac{\partial V_jV^j}{\partial V_i} = 2V^i</math>
De manière similaire, on écrira : <math>\partial^i = \frac{\partial}{\partial x_i}</math> et <math>\partial_i = \frac{\partial}{\partial x^i}</math>
En utilisant un temps quelconque indéterminé <math>\ t_0</math>, l'action <math>\ S = \int_{t_{0i}}^{t_{0f}}L_0 dt_0 </math> permet d'obtenir les équations d'Euler-Lagrange, relativistes mais obtenues de la même manière que dans le cas classique, avec une coordonnées de plus :
:<math>\frac{d~~ }{dt_0} \frac{\partial L_0}{\partial V_j } \ - \ \frac{\partial L_0}{\partial x_j} \ = \ 0~~</math> pour j=0;1;2;3
Il est important de remarquer que comme dans le cas classique, l'action et le lagrangien ne sont pas définis de manière unique : l'action est définie à l'addition près d'une fonction des extrémités du trajet et du temps, et le lagrangien est défini à l'addition près de la dérivée d'une fonction du temps (qui une fois intégrée donne une fonction des extrémités et du temps).
===Cas d'un corps libre===
====Lagrangien d'un corps libre====
Déterminons l'action et le lagrangien relativiste d'un corps libre.
Dans aucun référentiel galiléen la quadri-vitesse n'est nulle car le corps avance au moins dans la dimension temporelle.
Le lagrangien relativiste d'un corps libre doit, aux petites vitesses et en première approximation, être égal (peut-être à une constante additive près : l'ajout d'une constante ne change pas les équations d'Euler-Lagrange) au lagrangien classique.
Dans l'espace-temps de Minkowski, l'action détermine la trajectoire, et celle-ci ne dépend pas du référentiel d'où on l'observe. Donc l'action ne dépend pas des coordonnées, et, pour un corps libre, dépend seulement de la vitesse et est invariante par les transformations de Lorentz :
:<math>\ S = \int_{t_i}^{t_f}L dt </math> est invariant par les transformations de Lorentz
Dans le référentiel propre du corps, <math>\ dt = dt_0</math> est la variation du temps propre du corps ; <math>\ L = L_0</math> et la vitesse spatiale du corps est nulle. Dans un référentiel galiléen, et avec l'hypothèse que le corps est libre, la quadri-vitesse est constante dans le temps (et n'est jamais nulle) donc le lagrangien aussi car il dépend de la seule vitesse.
Ainsi, dans le référentiel propre du corps le lagrangien propre, <math>\ L_0 </math>, est une constante dans le temps.
Vu depuis un autre référentiel galiléen, se déplaçant par rapport au référentiel propre à la vitesse spatiale <math>\ v </math> constante, on a : <math>\ dt_0 = \sqrt{1-\frac{v^2}{c^2}}.dt </math>
Donc : <math>\ S = \int_{t_i}^{t_f} L_0.\sqrt{1-\frac{v^2}{c^2}} dt </math>
où <math>\ v = </math> vitesse spatiale relative entre le référentiel et le référentiel propre du corps = vitesse spatiale du corps dans le référentiel.
Donc : <math>\ L = L_0.\sqrt{1-\frac{v^2}{c^2}} \approx \ L_0.(1 - \frac{v^2}{2c^2}) = L_0 - \frac{L_0}{2c^2}v^2</math> par l'approximation aux petites vitesses devant <math>\ c</math>.
En comparant au lagrangien classique <math>\ L = \frac{1}{2}mv^2</math> (qui n'est pas réellement modifié par l'ajout de la constante <math> L_0 </math>) , on obtient : <math>- \frac{L_0}{2c^2} = \frac{1}{2}m </math> , d'où <math>\ L_0 = -mc^2</math>
Conclusion : dans un référentiel galiléen quelconque, le lagrangien est
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\ L = -mc^2.\sqrt{1-\frac{v^2}{c^2}} </math>
}}
où <math>\ v </math> est la vitesse spatiale du corps dans ce référentiel.
====Impulsion et énergie====
*Par définition de l'impulsion <math> \vec{p}</math>, on a : <math>\vec{p} = \frac{ \partial L}{\partial \vec{v}} = \frac{m\vec{v}}{\sqrt{1-\frac{v^2}{c^2}}}</math>
:L'énergie est définie par : <math>\ E = \vec{p}.\vec{v} - L </math>
:On obtient : <math>\ E = \frac{mc^2}{\sqrt{1-\frac{v^2}{c^2}}}</math>
:En particulier, pour <math>v=0</math>, l'''énergie au repos'' est <math>\ E = mc^2</math>
*En exprimant l'énergie en fonction de l'impulsion, on obtient : <math>\ E^2 = p^2.c^2 + m^2.c^4</math> ou encore <math>\ m^2.c^4 = E^2 - p^2.c^2 </math>
*On remarquera que bien qu'ayant la dimension d'une énergie, le ''lagrangien relativiste n'est pas l'énergie cinétique'' : cette dernière vaut
:<math>\ E_c = \frac{mc^2}{\sqrt{1-\frac{v^2}{c^2}}} - mc^2 = mc^2.\left(\frac{1}{\sqrt{1-\frac{v^2}{c^2}}} - 1\right)</math>
:On a bien <math> \ E_c \ </math> ≈ <math> \frac{1}{2} m \ v^2 \ </math> à l'approximation aux petites vitesses devant <math>\ c</math>
====Avec la quadri-écriture====
*On constate que <math>\ L.dt = -mc^2.\sqrt{1-\frac{v^2}{c^2}}.dt = -mc.\sqrt{(c.dt)^2-(dx_1)^2-(dx_2)^2-(dx_3)^2} = -mc\sqrt{dx_idx^i} = -mc.ds</math> , en utilisant l'égalité <math> v_i = \frac{dx_i}{dt}</math> et la définition adéquate de <math>\ ds</math> , appelé « temps propre » du corps.
:En utilisant le fait que <math>\ ds = \sqrt{dx_idx^i}</math> est le temps propre du corps, l'action minimisée entre deux points de l'espace-temps <math>\ S = -mc\int_A^Bds</math> montre que '''le chemin suivi par la particule pour aller du point A au point B est celui qui maximise le temps propre''', car le terme négatif <math>-mc</math> transforme la minimisation de <math>\ S = -mc\int_A^Bds</math> en maximisation de <math>\int_A^Bds</math>.
:En factorisant par <math>\ dt_0</math> , un temps quelconque paramétrant le système (et n'est donc pas obligatoirement le temps propre), on obtient : <math>\ L.dt = -mc.\sqrt{V_iV^i}.dt_0 = L_0.dt_0 </math>
:, en utilisant la quadri-vitesse pas obligatoirement propre <math>\ V = (V_0;V_1;V_2;V_3)</math> définie par : <math>\ V_i = \frac{dx_i}{dt_0}</math> ,avec <math>\ x_0 = c.t</math>.
:Le '''lagrangien relativiste''' d'une particule libre, paramétrée par le temps quelconque <math>\ t_0</math>, s'exprime donc :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\ L_0 = -mc. \sqrt{V_iV^i} = -mc.\frac{ds}{dt_0}</math>.
}}
*On se rappelle que la quadri-impulsion, comme l'impulsion, est définie par <math>P^i = \frac{\partial L_0}{\partial V_i} </math>
:D'où : <math> P^i = -mc\frac{V^i}{\sqrt{V^jV_j}} </math>
:Pour <math>i=0</math> , on obtient : <math>P^0 = -mc.\frac{(\frac{c.dt}{dt_0})}{(\frac{ds}{dt_0})} = -mc.\frac{c.dt}{ds} = -\frac{mc}{\sqrt{1-\frac{v^2}{c^2}}} = -\frac{E}{c} </math>
:Pour <math>i=1;2;3</math> , de manière similaire au cas i=0, on obtient : <math>(P^1,P^2;P^3)= \frac{m\vec{v}}{\sqrt{1-\frac{v^2}{c^2}}} = \vec{p}</math>.
:Le carré de la "norme" de la quadri-impulsion est <math>\ P^jP_j = (P_0)^2-(P_1)^2-(P_2)^2-(P_3)^2 = \frac{(-E)^2}{c^2} - (\vec{p})^2</math>, et aussi <math>\ P^iP_i = m^2.c^2.\frac{V^iV_i}{(\sqrt{V^jV_j}\ )^2} = m^2.c^2</math>
:D'où la formule déjà vue : <math>\ m^2.c^4 = E^2 - p^2.c^2</math>
*La constance de la quadri-impulsion, démontrée à partir des équations d'Euler-Lagrange, permet de montrer que l'énergie E et l'impulsion spatiale sont constantes par rapport au temps <math>\ t_0~</math>.
*La constante <math>\ P^iV_i-L</math> par rapport au temps <math>\ t_0~</math> est en fait la constante 0 ; une petite manipulation permet d'en déduire l'égalité déjà vue <math>\ m^2.c^4 = E^2 - p^2.c^2</math>.
:En calculant <math>\ P^iV_i-L</math> pour i∈{1;2;3} on retrouve l'énergie et l'égalité déjà citées.
*On montre facilement que quel que soit le temps <math>\ t_0~</math> choisi, s'il est celui d'un repère galiléen, la vitesse <math>\ V_i </math> et la "pseudo-norme" <math>\sqrt{V^jV_j}</math> sont constantes par rapport au temps <math>\ t_0~</math> : c'est une conséquence directe de la définition des repères galiléens, et du fait que le corps est libre.
*Dans le cas particulier où <math>\ dt_0</math> est le temps propre <math>\frac{ds}{c}</math>, alors les quadri-vitesse et quadri-impulsion sont propres, et on a l'égalité particulière <math>\sqrt{V^jV_j} = \frac{ds}{dt_0} = c.\frac{ds}{ds} = c</math> , qui peut être embarrassante pour l'utilisation des dérivées partielles dans le travail ci-dessus, et qui donne <math>\ P^0 = -m.V^0 = -\frac{E}{c}</math> et pour i=1;2;3 <math> \ P^i = m.V^i = p_i</math> .
===Cas d'un corps dans un [[champ électromagnétique]]===
====Sans la quadri-écriture====
Comme dans le cas classique, le lagrangien peut être défini en utilisant un potentiel électromagnétique <math> (\phi(\vec{q},t), \vec{A}(\vec{q},t))</math> :
<math> L = L(\vec{q},\vec{v},t) = \ -m.c^2.\sqrt{1-\frac{v^2}{c^2}} - e\phi(\vec{q},t) + e. \vec{v} \cdot \vec{A}(\vec{q},t) </math>
En prenant encore <math> \vec{p} = \frac{m\vec{v}}{\sqrt{1-\frac{v^2}{c^2}}} </math> , les équation d'Euler-Lagrange donnent :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> \dot{\vec{p}}= e \left(\vec{E} + \vec{v} \wedge \vec{B}\right)</math>
}}
<center>ce qui n'est pas une égalité pratique à utiliser car la dérivation de <math> \vec{p}~~</math> est laborieuse.</center>
L'impulsion est définie par : <math>\vec{P} = \frac{ \partial L}{\partial \vec{v}} = \frac{m\vec{v}}{\sqrt{1-\frac{v^2}{c^2}}} + e.\vec{A} = \vec{p} + e.\vec{A}</math>
On prendra donc soin de distinguer <math>\vec{p}</math> et <math>\vec{P}</math>
L'énergie est définie par : <math>\ E = \vec{P}.\vec{v} - L </math>
On obtient : <math>\ E = \frac{mc^2}{\sqrt{1-\frac{v^2}{c^2}}} + e.\phi </math>
Et après quelques calculs pour exprimer l'énergie en fonction de l'impulsion : <math>\ (E-e\phi)^2 = \left(\vec{P} - e.\vec{A}\right)^2.c^2 + m^2.c^4</math>
Toutes les approximations aux petites vitesses devant <math>c</math> redonnent les résultats classiques.
À partir du potentiel électromagnétique, le premier groupe des [[équations de Maxwell]] se démontre sans difficulté : l'équation de Maxwell-Faraday et l'équation de conservation du flux magnétique.
====Avec la quadri-écriture====
Le champ électromagnétique se manifeste sous forme d'un quadri-vecteur, appelé quadri-potentiel électromagnétique, <math>\ A^j </math> dont l'interaction avec la particule de charge <math>\ e</math> se manifeste sous forme lagrangienne par <math>\ e.A^j.dx_j </math>
La définition de l'action relativiste infinitésimale d'un champ électromagnétique est donc <math>\ L.dt = \ -mc.\sqrt{dx_j.dx^j} - e.A^j.dx_j </math> .
On pose <math> F^{ij} = \partial^i A^j - \partial^j A^i</math> tenseur champ électromagnétique.
En prenant <math>\ t_0</math> le temps propre de la particule, les équations d'Euler-Lagrange donnent les équations du mouvement de la particule :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>m. \frac{dV^i}{dt_0} = e.V_j.F^{ij} </math>
}}
Que l'on peut écrire aussi :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>mc. \frac{dV^i}{ds} = e.V_j.F^{ij} </math>
}}
En prenant : champ électrique = <math>\vec{E} = c.\left(F^{01},F^{02},F^{03}\right)</math> et champ magnétique = <math>\vec{B} = \left(F^{23},F^{31},F^{12}\right)</math> , on retrouve la [[force de Lorentz]] sous son écriture habituelle.
{{Boîte déroulante|titre=Démonstration des équations du mouvement d'une particule dans un champ électromagnétique|contenu=
<center>
On considère un champ de force dont l'interaction avec la particule de charge <math>\ e</math> se manifeste sous forme lagrangienne par <math>\ e.A^j.dx_j </math>
La définition de l'action relativiste infinitésimale d'un champ électromagnétique est <math>\ L.dt = \ -mc.\sqrt{dx_j.dx^j} - e.A^j.dx_j </math> .
En factorisant par <math>\ dt_0</math> un temps paramétrant le système, pas obligatoirement propre, on obtient : <math> L_0.dt_0 = ( -mc.\sqrt{V_j.V^j} - e.A^j.V_j).dt_0 </math>
Les équations d'Euler-Lagrange relativistes
:<math>~\frac{d~~ }{dt_0} \frac{\partial L_0}{\partial V_i } \ - \ \frac{\partial L_0}{\partial x_i} \ = \ 0</math>
donnent : <math>\frac{d~~ }{dt_0} (-mc\frac{2V^i}{2\sqrt{V_j.V^j}}- e.A^i) \ + e\ V_j \frac{\partial A^j}{\partial x_i} \ = \ 0</math>
En prenant dès maintenant <math>\ t_0 = </math> temps propre, on utilise l'égalité <math> \sqrt{V^jV_j} = c ~ </math> qui simplifie la dérivation
(mais on peut dériver avant si on veut plus de cohérence dans le suivi de la méthode), et on a :
<math>-m. \frac{dV^i}{dt_0} - e.\frac{dA^i}{dt_0} + e. \ V_j \frac{\partial A^j}{\partial x_i} \ = \ 0</math>
Sachant que <math>\frac{dA^i}{dt_0} = \frac{\partial A^i}{\partial x_j}. \frac{dx_j}{dt_0} = \partial^j A^i.V_j \ \ </math> et en posant <math> F^{ij} = \partial^i A^j - \partial^j A^i</math> tenseur champ électromagnétique,
on obtient l'équation du mouvement : <math>m. \frac{dV^i}{dt_0} = e.V_j.F^{ij} </math>
On retrouve le cas sans quadri-écriture en factorisant par <math>\ c.dt</math> le temps du référentiel, en prenant
<math>\ A^0 = A_0 = \frac{\phi}{c}~~</math> et <math> (A^1;A^2;A^3) = - (A_1;A_2;A_3) = - \vec{A}</math>.
</center>
}}
*On se rappelle que la quadri-impulsion, comme l'impulsion, est définie par <math>P^i = \frac{\partial L_0}{\partial V_i} </math>
:D'où : <math> P^i = -mc\frac{V^i}{\sqrt{V^jV_j}} - e.A^i</math>
:Pour <math>i=0</math> , on obtient : <math>P^0 = -mc.\frac{(\frac{c.dt}{dt_0})}{(\frac{ds}{dt_0})} - e.A^0 = -mc.\frac{c.dt}{ds} - e.\frac{\phi}{c} = -\frac{mc}{\sqrt{1-\frac{v^2}{c^2}}} - \frac{e}{c}.\phi = -\frac{E}{c} </math>
:Pour <math>i=1;2;3</math> , de manière similaire au cas i=0, on obtient : <math>(P^1,P^2;P^3)= \frac{m\vec{v}}{\sqrt{1-\frac{v^2}{c^2}}} + e.\vec{A} = \vec{p} + e.\vec{A} = \vec{P}</math>.
:De manière similaire au cas d'un corps libre, la constante <math>\ P^iV_i-L</math> par rapport au temps <math>\ t_0~</math> est en fait la constante 0, et une petite manipulation permet d'en déduire l'égalité déja vue <math>\ m^2.c^4 = (E - e.\phi)^2 - \left(\vec{P} - e\vec{A}\right)^2.c^2</math>.
====L'invariance de jauge du potentiel et du tenseur électromagnétiques====
On remarque que si à la place du quadri-potentiel électromagnétique <math>\ A^j </math>, on a le quadri-potentiel <math>\ A'^j = A^j + \partial^j \phi </math> où <math>\ \phi </math> est une fonction quelconque des coordonnées, alors le lagrangien devient <math> L' = L + \partial^j \phi.dx_j</math> et l'action <math>S' = S + \int^{t_f}_{t_i}\partial^j \phi.dx_j = S + \phi (x_i(t_f)) - \phi (x_i(t_i))</math>
En appliquant la méthode variationnelle qui fait varier le chemin en gardant les extrémités fixes, le terme <math>\ \phi (x_i(t_f)) - \phi (x_i(t_i))</math> est éliminé. Donc les deux potentiels <math>\ A^j </math> et <math>\ A'^j = A^j + \partial^j \phi </math> donnent les mêmes équations du mouvement : on appelle cela l'« invariance de jauge ».
On constate d'ailleurs que dans les équations du mouvement, le tenseur électromagnétique, terme représentant l'influence du champ électromagnétique, est bien invariant de jauge : <math> F'^{ij} = \partial^i A'^j - \partial^j A'^i = \partial^i A^j - \partial^j A^i + \partial^i \partial^j \phi- \partial^j \partial^i \phi = \partial^i A^j - \partial^j A^i = F^{ij}</math>
par le théorème de Schwarz : <math>\partial^i \partial^j = \partial^j \partial^i</math>.
===Cas d'un champ « de force »===
En physique classique, l'influence d'un corps sur un autre se transmet instantanément ; avec l'arrivée de l'électromagnétisme de Maxwell et plus encore avec celle de la relativité restreinte, l'influence se transmet au maximum à la vitesse de la lumière (dans le vide).
Ainsi, entre le corps influent et le corps infuencé, il se balade quelque chose dans l'espace, en général à la vitesse de la lumière, qui se répand dans l'espace et dont l'effet est un changement de trajectoire du corps influencé.
Suivant quelles propriétés ce ''champ'' (appelé ainsi car il a tendance à occuper l'espace) est-il créé, se déplace-t'il, est-il influencé par son environnement, etc ?
On peut répondre à ces questions à l'aide du principe de moindre action.
====Densité lagrangienne et équations d'Euler-Lagrange associées====
*Un champ est caractérisé par une étendue importante dans l'espace, on ne peut donc pas le repérer par les coordonnées <math>(t,\vec{q},\vec{v})</math>, mais on peut le repérer (ou plutôt le quantifier) par ses projections <math>\ (A_0,A_1,A_2,A_3)</math> sur les axes <math> (x_0,x_1,x_2,x_3)</math> et par les variations <math>\frac{\partial A_i}{\partial x_j} = \partial^j A_i</math> de ses projections (en supposant par avance que nous pourrons en déduire les dérivées secondes, comme dans le cas d'un corps localisé).
:Nous utiliserons donc <math>\ (A_0,A_1,A_2,A_3)</math> et <math>\frac{\partial A_i}{\partial x_j} = \partial^j A_i</math> , avec i,j ∈{0,1,2,3} pour un champ de la même manière que les quadri-coordonnées et la quadri-vitesse pour un corps localisé, les coordonnées <math>(x_0,x_1,x_2,x_3)</math> jouant le rôle de paramètres, comme seul le temps le faisait avant.
*L'action d'un champ est donc de la forme : <math>S = \int_V \Lambda(A_i;\partial^jA_i)d\Omega</math>
:Où V est le quadri-volume dans lequel on va appliquer la méthode variationnelle, <math>\ \Lambda</math> est appelé la « densité lagrangienne » et <math>\ d\Omega = dx_0.dx_1.dx_2.dx_3 = c.dt.dx_1.dx_2.dx_3</math>
*Par une démonstration semblable à celle déja vue dans le cas d'un corps localisable, et en utilisant la convention de sommation d'Einstein, on obtient les équations d'Euler-Lagrange pour la densité lagrangienne :
{{Bloc emphase|fond=#fff8ff;border:4px double #888| <math>\partial^j \frac{\partial \Lambda}{\partial (\partial^jA_i)} - \frac{\partial \Lambda}{\partial A_i} = 0</math>
}}
{{Boîte déroulante|titre=Démonstration des équations d'Euler-Lagrange pour une densité lagrangienne|contenu=
{{...}}
''bientôt...''
}}
====Tenseur impulsion-énergie d'un champ====
La densité lagrangienne <math>\Lambda = \Lambda(A_j;\partial_kA_j)</math> d'un champ étant donnée,
{{Boîte déroulante|titre=Détails|contenu=
<center>
On a : <math>\partial_i \Lambda = \frac{\partial \Lambda}{\partial A_j}. \partial_i A_j + \frac{\partial \Lambda}{\partial (\partial_k A_j)}.\partial_i(\partial_k A_j)</math>
puis en utilisant les équations d'Euler-Lagrange,
le lemme de Schwarz et enfin l'égalité <math>\partial_i \Lambda = \delta^k_{~i}\partial_k\Lambda</math>,
on obtient : <math>\partial_k (~\frac{\partial \Lambda}{\partial (\partial_kA_j)}.\partial_i A_j - \delta^k_{~i}\Lambda~) = 0</math>
</center>
}}
en posant
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>T^k_{~i} = \frac{\partial \Lambda}{\partial(\partial_kA_j)}.\partial_iA_j - \delta^k_{~i}\Lambda</math>
}}
tenseur « impulsion-énergie », on a :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\partial_kT^k_{~i} = 0</math>
}}
ce qui exprime sa conservation.
En posant : <math>\ w = T^0_{~0} = \frac{\partial \Lambda}{\partial(\partial_0A_j)}.\partial_0A_j - \Lambda</math> la densité d'énergie et <math> P^k = c.T^k_{~0} = c.\frac{\partial \Lambda}{\partial(\partial_kA_j)}.\partial_0A_j </math> pour k∈{1;2;3} les composantes du vecteur <math>\vec P</math>.
Nous avons alors les deux équations équivalentes
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\partial_kT^k_{~0} = 0</math>
}}
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\frac{\partial w}{\partial t} + div(~\vec{P}~) = 0</math>
}}
qui est l' « équation de conservation de l'énergie » : localement, la variation dans le temps de la densité d'énergie <math>\ w</math> est égale à l'opposée de la variation de densité d'impulsion par les composantes spatiales <math>\ P_i</math>.
====Densité lagrangienne d'un champ électromagnétique libre====
La densité lagrangienne du champ électromagnétique est :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\Lambda_{em} = -\frac{1}{4\mu_0}F^{ij}F_{ij}</math>
}}
{{Boîte déroulante|titre=Ébauche de démonstration de l'expression de la densité lagrangienne du champ électromagnétique|contenu=
*On cherche la densité lagrangienne du champ électromagnétique qui est composé du ou des nombres construits à partir du potentiel électromagnétique <math>\ A_i</math> i=0;1;2;3 qui sont invariants par changement de référentiel dans l'espace de Minkowski.
:La manifestation de ce potentiel est le champ électrique = <math>\vec{E} = c.(F^{01},F^{02},F^{03})</math> et le champ magnétique = <math>\vec{B} = (F^{23},F^{31},F^{12})</math>.
*On cherche donc les nombres invariants construits à partir de la matrice 4×4 associée : <math>\ F = (F^i_{~j})_{i,j=0,...,3}</math>
:Les coefficients du polynôme caractéristique sont invariants par la transformation d'endomorphismes <math> F \to U^{-1}.F.U </math>, où <math>\ U</math> est la matrice de Lorentz du changement de base dans l'espace de Minkowski, c'est-à-dire de référentiel galiléen en relativité restreinte.
:Les invariants de cette matrice par les changements de base sont les coefficients de son polynôme caractéristique <math>\ P_F(x) = det(F-x.Id)</math>
:On sait, par la définition de ses coefficients, que cette matrice est anti-symétrique : <math>\ ^tF = -F</math> , où <math>\ ^t</math> est la transposition de la matrice. On en déduit que <math>\ P_F(x)</math> est un polynôme pair. Soit : <math>\ P_F(x) = x^4 + a.x^2 + det(F)</math>.
:Par quelques calculs, on montre que <math>\ a = - F^i_{~j}F_{~i}^j = F^{ij}F_{ij} = -2.( \frac{\vec{E}^2}{c^2} - \vec{B}^2)</math> et que <math>\ det(F) = - \frac{1}{c^2} .(\vec{E}.\vec{B})^2</math>.
:Il y a donc deux nombres à examiner : <math>\ F^{ij}F_{ij} = -2( \frac{\vec{E}^2}{c^2} - \vec{B}^2)</math> et <math>\frac{1}{c^2}.(\vec{E}.\vec{B})^2</math>.
* Pour une raison de dimension , le nombre <math>\frac{1}{c^2}.(\vec{E}.\vec{B})^2</math> ne convient pas (la densité lagrangienne doit avoir la dimension d'une densité d'énergie), il faut essayer avec <math>\sqrt{\frac{1}{c^2}.(\vec{E}.\vec{B})^2} = \frac{1}{c}.|\vec{E}.\vec{B}| = \frac{+/- 1}{c}.\vec{E}.\vec{B}</math>.
:Mais on peut montrer que <math>\vec{E}.\vec{B}</math> est la différentielle d'une fonction, donc son ajout à la densité lagrangienne ne changerait en rien les équations d'Euler-Lagrange. Le nombre <math>(\vec{E}.\vec{B})^2</math> peut donc être écarté.
*De ce fait, <math>\ F^{ij}F_{ij}~</math> est l'unique candidat à être la densité lagrangienne.
*En choisissant un coefficient multiplicateur qui détermine les unités de mesure du champ électromagnétique, on prend : <math>\Lambda_{em} = -\frac{1}{4\mu_0}F^{ij}F_{ij}</math>
:le signe <math>\ -</math> se justifiant par des considérations liées à la minimisation de l'action.
}}
====Les équations du champ électromagnétique====
L'hypothèse de ce paragraphe est qu'il y a un courant de particules (voir même d'une seule particule) non influencé par le champ électromagnétique. Avec cette condition, on étudie les modifications du champ.
{{Boîte déroulante|titre=Détails pour déterminer la densité lagrangienne à utiliser|contenu=
On sait (par hypothèse) que l'interaction entre une particule et le champ se modélise, sous forme lagrangienne, par <math>\ e.A^j.dx_j </math> .
L'action est donc : <math>S = -\Sigma_{part}(~e.\int A^idx_i ) - \frac{1}{4\mu_0}\int_{t_i}^{t_f}dt \int F^{ij}F_{ij}.dx_1dx_2dx_3 = -\Sigma_{part}(~e.\int A^idx_i ) - \frac{1}{4c\mu_0}\int_V F^{ij}F_{ij}.d\Omega</math>
On a : <math> \Sigma_{part}(~e.\int A^idx_i ) = \Sigma_{part}(~e.\int A^i \frac{dx_i}{dt}.dt) = \int \rho dx_1dx_2dx_3 \int A^i \frac{dx_i}{dt}.dt = \frac{1}{c} \int_V A^i \rho \frac{dx_i}{dt}.d\Omega </math>
Ce que l'on peut écrire : <math> \Sigma_{part}(~e.\int A^idx_i ) = \frac{1}{c} \int_V A^i J_i.d\Omega </math>
Avec <math>\rho = \frac{d^3e}{dx_1dx_2dx_3}~~</math> la densité de charge et <math> J_i = \rho.\frac{dx_i}{dt} = (\rho .c~; \rho.\vec v ) = (\rho .c~; \vec j~)</math> le quadri-vecteur courant.
On obtient : <math>S = -\frac{1}{c} \int_V ( A^iJ_i + \frac{1}{4\mu_0}F^{ij}F_{ij}).d\Omega = \frac{1}{c} \int_V \Lambda d\Omega</math>
}}
La densité lagrangienne à utiliser est :
<center>
<math>\Lambda = -A^iJ_i - \frac{1}{4\mu_0}F^{ij}F_{ij} </math>
</center>
Les équations d'Euler-Lagrange donnent :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> \partial^i F_{ik} = \mu_0.J_k </math>
}}
Pour k = 0 , on obtient <math> \mathrm{div}\ \overrightarrow{E} \ = \ \frac{\rho}{\varepsilon_0}</math>
l'équation de Maxwell-Gauss ou [[équations de Maxwell|équation de conservation de la charge]].
Pour k∈{1;2;3} on obtient <math> \overrightarrow{\mathrm{rot}} \ \overrightarrow{B} \ = \ \mu_0 \overrightarrow{j} \ + \ \varepsilon_0 \mu_0 \ \frac{\partial \overrightarrow{E}}{\partial t}</math>
l'[[équations de Maxwell|équation de Maxwell-Ampère]].
De plus, à partir de <math> \partial^i F_{ik} = \mu_0.J_k </math>, et en utilisant l'anti-symétrie de <math>\ F_{ik}</math> et le théorème de Schwarz ( <math>~~\partial^{ik} = \partial ^{ki}~~</math> ), on obtient : <math> \mu_0. \partial^k J_k = - \mu_0. \partial^i J_i </math>
D'où les deux présentations de « l'équation de conservation de la charge » :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> \partial^k J_k = 0~~ </math>
}}
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>~~\mathrm{div}\ (~\vec{j}~) + \frac{\partial \rho}{\partial t} = 0</math>
}}
====Tenseur impulsion-énergie du champ électromagnétique====
{{...}}
==En [[relativité générale]]==
{{Référence nécessaire|On doit à [[David Hilbert]], en 1916, la première utilisation du principe de moindre action pour obtenir les équations de la relativité générale, notamment les équations du champ gravitationnel.}}
Pour la relativité générale aussi, les équations peuvent être obtenues sans faire appel au principe de moindre action : le [[principe d'équivalence]], exprimé sous la forme « on peut toujours trouver un référentiel annulant localement le champ de gravitation », permet de retrouver directement les équations du mouvement d'une particule ; et l'unicité de la forme du tenseur géométrique qui s'annule par la dérivée covariante, unicité prouvée par [[Élie Cartan]], permet de trouver les équations du champ de gravitation, ce qui fût la méthode originelle d'Einstein (bien que l'unicité en question n'était pas encore prouvée à l'époque).
Si les équations de la relativité générale sont données, on peut en déduire l'action permettant d'appliquer le principe. En particulier, avec les équations des géodésiques on peut retrouver la métrique <math>ds^2\,</math> associée.
===Particule===
====Particule dans un champ de [[gravitation]]====
Dans ce travail, on utilise l'hypothèse que la particule ne modifie pas son environnement : la masse de la particule ni sa position ne changent le champ de gravitation, cette masse doit donc être « petite ».
En vertu du [[principe d'équivalence]] d'Einstein, la gravitation est localement équivalente au choix d'un référentiel accéléré.
Dans le cadre de la relativité restreinte, en prenant un référentiel accéléré (coordonnées <math>\ (x'_0;x'_1;x'_2;x'_3)</math>), la perception locale est donc un champ de gravitation, et le changement de référentiel par rapport à un référentiel inertiel (coordonnées <math>\ (x_0;x_1;x_2;x_3)</math>) impose une [[métrique]] aux coefficients non triviaux : <math>\ ds^2 = (x_0)^2- (x_1)^2- (x_2)^2-(x_3)^2 = g^{ij}(x')x'_ix'_j</math>. Elle suffit pour déterminer les équations du mouvement dans ce référentiel du fait du principe de moindre action en relativité restreinte.
Le principe d'équivalence permet de dire qu'un champ gravitationnel réel (non dû au choix du référentiel) est aussi déterminé par la métrique <math>\ ds^2</math> (et la métrique est déterminée par le champ de gravitation) ; bien que l'utilisation d'une métrique <math>\ ds^2 = g^{ij}(x)x_ix_j = g^{ij}x_ix_j</math> qui ne soit pas causée, et donc pas compensable au delà d'un domaine local de l'espace-temps, par un changement de référentiel implique que l'espace-temps n'est pas euclidien (voir l'expérience par la pensée du disque en rotation, décrit dans [[relativité générale]]), et que l'on sort alors du cadre de la relativité restreinte pour construire une nouvelle théorie : la [[relativité générale]].
On peut donc rester dans la continuité de la relativité restreinte, et affirmer que l'action infinitésimale d'une particule ponctuelle, influencée par la seule gravitation, en relativité générale est :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>dS = -mc\sqrt{g^{ij}dx_idx_j}</math>
}}
où on suppose que <math>\ g^{ij} = g^{ji}</math> sans rien enlever à la généralité.
En utilisant le fait que <math>\ ds = \sqrt{g^{ij}dx_idx_j}</math> est le temps propre de la particule, l'action minimisée entre deux points de l'espace-temps <math>\ S = -mc\int_A^Bds</math> montre que, comme en relativité restreinte, c'est le temps propre pour aller du point A au point B qui est maximisé (localement) par le principe. '''Les géodésiques sont les chemins qui maximisent (localement) le temps propre de la particule'''.
Pour garder la cohérence physique, on a besoin de supposer que les <math>\ g^{ij}</math> sont continus ; pour pouvoir travailler avec des outils connus, c'est-à-dire des dérivations, mais aussi pour supposer que le champ gravitationnel est continu, on doit supposer qu'ils sont différentiables. Par la suite, pour les équations d'Einstein, il sera indispensable de supposer qu'ils sont C{{2}}.
En considérant un temps <math>t_0</math> quelconque :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\frac{dS}{dt_0} = L_0 = -mc\sqrt{g^{ij}V_iV_j}</math>
}}
On utilise toujours les équations d'Euler-Lagrange <math>\ \frac{d~~ }{dt_0} \frac{\partial L_0}{\partial V_k } \ - \ \frac{\partial L_0}{\partial x_k} \ = \ 0~~</math> après avoir divisé par le coefficient <math>\ -mc</math> ici inutile.
{{Boîte déroulante|titre=Détails de la démonstration|contenu=
<center>
On obtient : <math>\ \frac{d~~ }{dt_0}\left( \frac{2.g^{ik}V_i}{2.\sqrt{g^{ij}V_iV_j}} \right) - \ \frac{\partial^k g^{ij}.V_iV_j}{2.\sqrt{g^{ij}V_iV_j}} \ = \ 0~~</math>
En prenant dès maintenant <math>t_0 = </math> temps propre, on peut utiliser l'égalité <math>\sqrt{g^{ij}V_iV_j} = c</math> qui simplifie la dérivation <math>\ \frac{d~~ }{dt_0} </math>,
sans changer le résultat si on dérive avant, et on obtient
<math>\ \frac{d~~ }{dt_0} \left( \frac{g^{ik}V_i}{\sqrt{g^{ij}V_iV_j}} \right) = \frac{1}{c}.( \partial^n g^{ik}.V_nV_i + g^{ik}\dot{V_i}) </math>
En remarquant que <math>\partial^n g^{ik}.V_nV_i = \frac{1}{2}.(\partial^n g^{ik}.V_nV_i + \partial^i g^{nk}.V_nV_i)</math> , que nous utiliserons essentiellement par soucis d'esthétisme, et en changeant les indices pour n'utiliser que i, j et k,
Les équations d'Euler-Lagrange donnent : <math>g^{ik}\dot{V_i} + \frac{1}{2}.(- \partial^k g^{ij} + \partial^i g^{jk} + \partial^j g^{ik})V_iV_j = 0</math>
Avec l'égalité <math>\ g_{km}g^{ki} = \delta_m^i</math>
et le symbole de [[Elwin Bruno Christoffel|Christoffel]] : <math>\Gamma_m^{ij} = \frac{1}{2}.g_{km}(- \partial^k g^{ij} + \partial^i g^{jk} + \partial^j g^{ik})</math>
</center>
}}
On obtient l'équation :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\dot{V}_m + \Gamma_m^{ij}V_iV_j = 0</math>
}}
que l'on peut aussi écrire :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\frac{d^2x_k}{ds^2} + \Gamma_k^{ij}\frac{dx_i}{ds}\frac{x_j}{ds} = 0</math>
}}
ou encore :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> \frac{DV_k}{ds} = 0</math>
}}
avec la « dérivée covariante » : <math>DV_k = dV_k + \Gamma_k^{ij}V_idV_j </math> et <math>DV^k = dV^k + \Gamma^k_{ij}V^idV^j </math> , où <math>\ V_k = \frac{dx_k}{dt_0}</math> pour <math>\ t_0 = </math> temps propre.
Le symbole de [[Elwin Bruno Christoffel|Christoffel]] <math>\Gamma_k^{ij}</math> s'impose comme la manifestation de la gravitation dans les équations du mouvement.
Les équations du mouvement ne dépendent pas de la masse de la particule (nommée ainsi car nous avons négligé son étendue spatiale et son influence sur son environnement) : toutes les particules suivent les mêmes trajectoires (à conditions initiales identiques), c'est l'équation des géodésiques en relativité générale, en présence de la seule gravitation.
Toutefois, ces équations du mouvement ne sont pas valables pour une particule de masse nulle car dans ce cas, on a dès le départ <math>~dS = 0~~</math>, ce qui interdit tous les calculs menés ci-dessus ; on a aussi <math>~ds = c.dt_0 = 0~~</math> car le temps propre ne s'écoule pas pour une particule de masse nulle (voir [[Relativité restreinte]]), le terme <math>\dot{V}_m </math> ne peut en aucun cas avoir de sens. Il faut considérer l'onde associée à la particule pour avoir une équation ayant un sens, d'ailleurs la [[lumière]] était comprise comme une onde (électromagnétique) et comme une particule (le [[photon]], de masse nulle) lorsque la relativité générale a été écrite.
====Particule dans un [[champ électromagnétique]]====
De manière similaire à la relativité restreinte, la définition de l'action relativiste infinitésimale d'une particule ponctuelle de charge <math>\ e</math> dans un champ électromagnétique est <math>\ L.dt = \ -mc.\sqrt{g^{ij}dx_i.dx_j} - e.A^j.dx_j </math> .
Par des calculs parfaitement similaires, on en tire les équations du mouvement :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>m.(\dot{V}^k + \Gamma^k_{ij}V^iV^j) = e.V_j.F^{kj} </math>
}}
que l'on peut écrire :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> mc.\left( \frac{d^2x^k}{ds^2} + \Gamma^k_{ij}\frac{dx^i}{ds}\frac{dx^j}{ds} \right) = e.F^{kj}\frac{dx_j}{ds}</math>
}}
ou encore :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> mc. \frac{DV^k}{ds} = e.F^{kj}V_j</math>
}}
===Champ de [[gravitation]]===
Afin d'en déterminer la densité lagrangienne, puis les équations, il est nécessaire de développer un peu certaines considérations abordées ci-dessus, et même quelques nouvelles.
====Densité lagrangienne dans l'espace courbe====
Du fait de l'invariance de la trajectoire du champ par rapport aux référentiels d'où on l'observe, l'action qui la caractérise <math>S_g = \int L d \Omega</math> doit être invariante par changement de référentiel.
{{Boîte déroulante|titre=Détails justifiants la densité lagrangienne|contenu=
<center>
Soient <math>S_g = \int L d \Omega = \int L' d \Omega'</math> l'action dans deux référentiels différents.
On a : <math>\ d\Omega = dx_0.dx_1.dx_2.dx_3</math> et <math>\ d\Omega' = dx'_0.dx'_1.dx'_2.dx'_3 = J. dx_0.dx_1.dx_2.dx_3 </math>
où <math>\ J </math> est le [[Matrice jacobienne|jacobien]] du changement de variables.
On a : <math>J = \left | \det \left ( \frac{\part x'_i}{\part x_j} \right ) \right |~</math>
Or : <math>ds^2 = g^{ij}dx_idx_j = g'^{kl}dx'_kdx'_l \to \ g^{kl} = \frac{\part x'_i}{\part x_k} . \frac{\part x'_j}{\part x_l} g'^{ij} \to g = J^2.g'</math> , en prenant les [[Déterminant (mathématiques)|déterminant]]s.
Donc : <math>J = \frac{|g|^{\frac{1}{2}}}{|g'|^{\frac{1}{2}}}</math>
Ainsi <math>S_g = \int L d \Omega = \int L' d \Omega' = \int L' .J d \Omega \to L = L'.J \to L.|g|^{-\frac{1}{2}} = L'.|g'|^{-\frac{1}{2}} = \Lambda </math> est une constante du champ par rapport aux changements de référentiels.
L'objectif est donc de trouver les scalaires du champ, invariants par rapport aux changements de référentiels.
</center>
}}
En notant <math>\ \Lambda</math> le scalaire du champ, invariant par rapport aux changements de référentiels, la densité lagrangienne sera : <math>\ L = \Lambda .|g|^{\frac{1}{2}} </math>
====Définitions des [[tenseur]]s de [[Bernhard Riemann|Riemann]], de [[Gregorio Ricci-Curbastro|Ricci]], et de la [[courbure]]====
''À la manière d'[[Élie Cartan]]''
En termes mathématiques, l'espace quadri-dimensionnel défini par les considérations ci-dessus est une [[variété (géométrie)|variété]] C<sup>2</sup> où les quadri-vitesses sont des vecteurs appartenant à l'[[espace tangent|espace vectoriel tangent]] au point où on a dérivé, cet espace vectoriel étant muni de la métrique <math>\ g^{ij}</math>.
Rappelons que les coordonnées <math>(x_0;x_1;x_2;x_3)</math> sont les coordonnées des points de la variété, munie d'un système de coordonnées quelconque, représentant le choix arbitraire du référentiel physique de l'observateur.
La mesure de la gravitation, qui influe sur les géodésiques, peut se faire à travers la différence d'orientation entre deux vecteurs résultant du transport d'un seul vecteur d'origine par deux chemins géodésiques différents vers un même point final.
*L'équation des géodésiques <math>\dot{V}_m + \Gamma_m^{ij}V_iV_j = 0</math> est équivalente à <math>\frac{dV_k}{dt_0} = - \Gamma_k^{ij}V_i V_j</math>.
:Du fait que <math>V_j = \frac{dx_j}{dt_0}</math>, on déduit : <math>dV_k = - \Gamma_k^{ij} V_i dx_j </math> ; sachant que l'on a <math>\Gamma_k^{ij} = \Gamma_k^{ji}</math> comme on le voit à partir de sa définition, on pourrait aussi bien écrire <math>dV_k = - \Gamma_k^{ij} dx_i V_j </math>.
:De manière similaire, on obtient <math>dV^k = - \Gamma^k_{ij} V^i dx^j </math>
*Un vecteur <math> \left( A_i \right) </math> est dit ''transporté parallèlement'' le long d'une géodésique si les variations de ses coordonnées vérifient <math>dA_k = - \Gamma_k^{ij} A_i dx_j </math> quand il est déplacé de <math>\ (dx_j)_{j=0;1;2;3}</math> le long de la géodésique.
{{Boîte déroulante|titre=''Détails de la méthode d'Élie Cartan''|contenu=
*À partir d'un point M quelconque de la variété, considérons deux variations infinitésimales <math>\ d</math> et <math>\ \delta </math> le long de deux géodésiques quelconques, et considérons les deux trajets distincts qui utilisent alternativement l'une puis l'autre de ces géodésiques.
:1<sup>er</sup> trajet : <math>M(x_i) \to M_1(x_i+dx_i) \to M_2(x_i+dx_i+ \delta (x_i+dx_i))</math>
:2<sup>ème</sup> trajet : <math>M(x_i) \to M'_1(x_i + \delta x_i) \to M'_2(x_i + \delta x_i + d (x_i+ \delta x_i))</math>
:Afin que ces deux trajets aboutissent au même point, on suppose que <math>\ d \delta x_i = \delta d x_i</math>, ce qui est réalisable car les géodésiques utilisées à partir des points <math>\ M_1</math> et <math>\ M'_1</math> sont quelconques.
*Étudions les variations des coordonnées d'un vecteur <math> \left( A_i \right) </math> transporté parallèlement le long de chacun des chemins :
:1<sup>er</sup> trajet : <math>\ A_i \to A_i+dV_i \to A_i+dA_i+ \delta (A_i+dA_i) = A_i+dA_i+ \delta A_i + \delta dA_i = W_i</math>
:2<sup>ème</sup> trajet : <math>\ A_i \to A_i+ \delta A_i \to A_i+ \delta A_i+ d(A_i+ \delta A_i) = A_i+\delta A_i +dA_i+ d \delta A_i = W'_i</math>
:On a : <math> W_i - W'_i = \delta dA_i - d \delta A_i = \delta (- \Gamma_k^{ij} A_i dx_j) - d(- \Gamma_k^{ij} A_i \delta x_j) </math>
:Après quelques calculs, on obtient : <math> W_i - W'_i = (\partial^j \Gamma_i^{lk} - \partial^l \Gamma_i^{jk} + \Gamma_p^{lk} \Gamma_i^{jp} - \Gamma_p^{jk} \Gamma_i^{lp}) dx_j \delta x_k A_l</math>
*On définit le ''tenseur de [[Bernhard Riemann|Riemann]]'' par : <math>R_i^{j,kl} = \partial^j \Gamma_i^{lk} - \partial^l \Gamma_i^{jk} + \Gamma_p^{lk} \Gamma_i^{jp} - \Gamma_p^{jk} \Gamma_i^{lp} </math>
:L'égalité <math> W_i - W'_i = R_i^{j,kl} dx_j \delta x_k A_l</math> indique que ce tenseur mesure la différence entre deux vecteurs issus du même vecteur d'origine <math> \left( A_i \right) </math> par transport parallèle par deux chemins différents.
}}
*On définit le ''tenseur de [[Bernhard Riemann|Riemann]]'' par :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>R_i^{j,kl} = \partial^j \Gamma_i^{lk} - \partial^l \Gamma_i^{jk} + \Gamma_p^{lk} \Gamma_i^{jp} - \Gamma_p^{jk} \Gamma_i^{lp} </math>
}}
*Le ''tenseur de [[Gregorio Ricci-Curbastro|Ricci]]'' est une contraction du tenseur de Riemann : <math> R^{ij}=R_k^{i,kj}</math>
:Sa formule montre que c'est un tenseur symétrique : <math>\ R^{ij}=R^{ji}</math>
*La ''[[courbure]] riemannienne'' est le nombre obtenu par contraction du tenseur de Ricci : <math>\ R=g_{ij}R^{ij}</math>
*Toutes les égalités utilisées dans « ''détails de la méthode d'Élie Cartan'' » étant indépendantes du référentiel choisi, et c'est aussi le cas pour les définitions des tenseurs de Riemann et de Ricci (c'est d'ailleurs pourquoi on se permet de les nommer ''[[tenseur]]'' ). C'est aussi le cas de la courbure <math>\ R</math> qui est donc candidat pour être <math>\ \Lambda</math> le scalaire invariant du champ de gravitation.
*[[Élie Cartan]] a démontré que les scalaires invariants par changement de référentiel sont de la forme <math>\ \alpha R + \beta~</math>.
:<math>~\ \alpha </math> indique simplement qu'un changement d'unité est toujours possible, <math>\ \beta</math> permet d'introduire la [[constante cosmologique]].
====Outils analytiques====
=====Une application du [[Lois du mouvement de Newton|principe d'inertie]] dans l'espace courbe=====
Pour que notre travail soit bien une conséquence du principe de moindre action, la méthode utilisée ici consiste à déterminer les propriétés de la variété à partir de la métrique de ses espaces tangents.
*Les espaces vectoriels tangents (de dimension 4) sont munis de leur base « naturelle » {<math>\ \ { \vec{e}^{~0}; \vec{e}^{~1} ; \vec{e}^{~2}; \vec{e}^{~3}} </math>} : si <math>\ M(x_0;x_1;x_2;x_3)</math> est le point où l'on considère l'espace tangent, on pose <math>\vec{e}^{~i} = \left(~\frac{\partial x_j}{\partial x_i}~\right)_{j=0,1,2,3}</math> ; ce que l'on écrit souvent <math>\vec{e}^{~i} = \frac{\partial ~}{\partial x_i} </math>.
:Les équations des géodésiques sont des propriétés concernant les coordonnées <math>\frac{dx_i}{dt_o}</math> ou <math>\frac{dx_i}{ds}</math> de la quadri-vitesse le long de cette trajectoire, elles ne donnent pas d'indication pour la variation (la dérivation) d'un ''quadri-vecteur'' <math>\vec{e}^{~i}</math> d'un point à un autre de l'espace, ni même pour la dérivation du quadri-vecteur vitesse <math>\vec V = V_i \vec e^{~i}</math>.
:Pour cela, nous pouvons utiliser un principe physique réécrit sur mesure pour la relativité générale :
*[[Lois du mouvement de Newton|Principe d'inertie]] : le long d'une géodésique, et en l'absence d'intervention extérieure, le (quadri-)vecteur vitesse d'une particule est constant.
:C'est-à-dire : <math>d\vec V = \vec 0</math>
:On en tire : <math>d\vec V = \vec 0 = dV_i .\vec e^{~i} + V_i .d\vec e^{~i} = - \Gamma_i^{jk} dx_j V_k. \vec e^{~i} + V_i .d\vec e^{~i} </math>
:Le quadri-vecteur vitesse initial étant quelconque, on obtient :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math> d\vec e^{~i} = \Gamma_k^{ij} dx_j \vec e^{~k}</math>
}}
:En analysant les équations des géodésiques ou en tenant compte du fait que les « axes » des coordonnées ne sont pas obligatoirement des géodésiques, on ne peut pas affirmer que les coordonnées du quadri-vecteur vitesse sont constantes.
{{Boîte déroulante|titre=À propos du choix <math>d\vec V = \vec 0</math>|contenu=
*Dériver signifie « déterminer la droite qui indique la direction du mouvement ». Tout le problème est de savoir ce qu'est une droite quand le système de coordonnées est quelconque, voire dans un espace courbe ; une fois les droites déterminées, la dérivation peut être définie.
*Dans le cadre qui nous intéresse, quand l'expérimentateur est dans un espace de Minkowski et qu'il a choisi un système de coordonnées quelconque, ce qui y induit éventuellement une gravitation, les droites de la dérivation sont celles de l'espace de Minkowski, qui sont aussi celles du mouvement inertiel. À moins de définir une nouvelle dérivation, l'égalité <math>d\vec V = \vec 0</math> s'impose.
*Quand l'expérimentateur est dans un référentiel où il y a de la gravitation, et en l'absence d'information sur les causes de cette gravitation (dûe à une masse ou dûe à un référentiel accéléré, ou les deux) les seules droites auxquelles il a accès, en tant que physicien, sont celles du mouvement inertiel : la dérivation est donc définie par <math>d\vec V = \vec 0</math>.
:Mais ce choix est basé sur l'hypothèse que, dans son référentiel, le mouvement inertiel suit bien une droite. Si l'expérimentateur choisit comme droites les axes de son référentiel, il impose donc <math>d\vec e^{~i} = \vec 0</math>, le mouvement « inertiel » observé n'est pas droit (<math>d\vec V \ne \vec 0</math>) et est interprétable comme dû à une force (de gravitation).
:Ces deux choix, comme d'autres que l'on peut imaginer, ne sont valables que localement :Le premier assimile localement la gravitation à un référentiel accéléré dans un espace de Minkowski, le deuxième émet l'hypothèse d'une force dans un espace initialement droit ; deux choix qui redressent à leur manière l'espace-temps, ce qui ne peut se faire que localement.
}}
=====La dérivée covariante=====
Soit <math>\vec{A}(x) = A_i\vec e^{~i}</math> un quadri-vecteur dans l'espace tangent au point <math>\ M(x_0;x_1;x_2;x_3)</math>.
On a : <math>d \vec{A}(x) = (dA_i) \vec e^{~i} + A_i d(\vec e^{~i}) = ( \partial^j A_i + A_k \Gamma_i^{jk})\vec e^{~i}dx_j = D^j A_i .\vec e^{~i}dx_j</math>
En définissant la ''dérivée covariante'' par :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>D^j A_i = \partial^j A_i + \Gamma_i^{jk} A_k </math>
}}
Propriété :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>D^j A_{il} = \partial^j A_{il} + \Gamma_i^{jk} A_{kl} + \Gamma_l^{jk} A_{ik}</math>
}}
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>D^j A_i^l = \partial^j A_i^l + \Gamma_i^{jk} A_k^l - \Gamma_k^{jl} A_i^k</math>
}}
Et ainsi de suite avec tous les indices d'un tenseur, suivant leurs positions.
=====Où l'on retrouve les tenseurs de Riemann, etc.=====
À l'aide de la dérivée covariante, et après quelques calculs, on trouve : <math> \left( D^iD^j -D^jD^i \right) A_k = R_k^{l,ij} dx_i dx_j A_l</math>.
On obtient donc les notions déjà introduites « à la manière d'Élie Cartan ».
=====Égalités et propriétés utiles=====
*Théorème de Ricci : <math>\ D_kg^{ij} = 0~ \quad ~</math> et <math>~\quad D_kg_{ij} = 0~</math>
*En posant <math>\ g = \det (g^{ij}) \qquad </math>, on a : <math> |g| = -g \qquad g^{ij}.g_{ij} = \delta^i_i = 4 \qquad \ dg = g~g_{ij}~dg^{ij}</math>
*Théorème d'Ostrogradski : <math>\int_V \sqrt {-g}~ D_iA^i~d\Omega = \oint_{\part V} \sqrt {-g} A^i~dS_i</math> , quand <math>\ A^i</math> est un tenseur.
{{Boîte déroulante|titre=''Ébauches de démonstrations des égalités''|contenu=
{{...}}
}}
*La somme, la différence et la sommation d'Einstein de [[tenseur]]s définis dans le même [[espace tangent]] donnent un tenseur ; par contre s'il s'agit de tenseurs définis dans des espaces tangents différents, il n'est pas sûr que cela donne un tenseur.
:Par exemple : le symbole de [[Elwin Bruno Christoffel|Christoffel]] est défini à partir du tenseur métrique. L'équation des géodésiques <math>\Gamma_i^{jk}.V_k = \part^jV_i</math> nous montre qu'il peut être défini à l'aide de <math>\ \part^j V_i</math> qui, bien que tenseur, est construit par une différence entre deux tenseurs (les quadri-vecteurs <math>\ V_l(x_m) </math> et <math>\ V_l(x_m + dx_m) </math> ) définis dans deux espaces tangents différents : le symbole de Christoffel, lui, n'est pas un tenseur (sauf cas particuliers), comme on peut le montrer à l'aide de sa formule de définition.
*Une égalité tensorielle démontrée en un point quelconque, mais en utilisant un référentiel particulier, est une égalité vraie en ce point et pour tous les référentiels : c'est là le principal intérêt d'utiliser des tenseurs.
:Par exemple, en tout point il existe un référentiel en apesanteur (en chute libre dans le champ de pesanteur), c'est à dire pour lequel <math>\Gamma_i^{jk} = 0</math>. Dans un tel référentiel, on a <math>R_i^{j,kl} = \partial^j \Gamma_i^{lk} - \partial^l \Gamma_i^{jk} </math> et <math>D^j A_i = \partial^j A_i </math> quand <math>\ A_i</math> est un tenseur : ce qui est plus simple à utiliser pour justifier une égalité tensorielle qui sera vraie quel que soit le référentiel.
====Les équations d'Einstein du champ de gravitation dans le ''cas extérieur''====
Les tenseurs sont utilisés pour s'assurer que les égalités sont vraies quel que soit le point d'observation du physicien et quelque soit son référentiel. Les tenseurs ne transportent que des informations liées au point d'observation et à son espace tangent, du coup, les informations qui y sont utilisées et qui en sont produites ne sont que locales : ce sont des informations sur les tenseurs, mis à part les données universellement valables comme les constante c, G, et autres que l'on pourra y trouver.
Le premier cas des équations du champ est le cas où il n'y a pas de matière (localement) : on parle du « cas extérieur », sous entendu « à la matière ».
Dans ce cas, la seule composante de l'action est la composante du champ gravitationnel <math>\ S_g = K. \int \sqrt{-g}.R.d\Omega</math> , où <math>\ K</math> est une constante liée au choix des unitées : pour les unités MKSA, on prend <math>\ K = - \frac{c^3}{4 \pi G}</math> , le signe <math>\ -</math> étant dû au principe de minimisation de l'action.
Pour trouver les équations du champ de gravitation sous la forme de tenseurs de densité d'énergie qui soient symétriques, il est plus simple de transformer le lagrangien sous l'intégrale de l'action que d'utiliser les équations d'Euler-Lagrange. Le principe variationnel est appliqué en faisant varier les termes de la métrique <math>\ g^{ij}</math>, qui est la manifestation lagrangienne de la gravitation, d'après le principe d'équivalence tel qu'appliqué plus haut.
{{Boîte déroulante|titre=''Démonstration des équations d'Einstein dans le cas extérieur''|contenu=
<center>
En utilisant l'égalité <math>\ R = g_{ij}R^{ij}</math>, on a
<math>\delta S_g = K.\int \delta \left( \sqrt{-g}.g_{ij}.R^{ij} \right)d\Omega </math>
<math>= K [ \int \delta ( \sqrt{-g})g_{ij}R^{ij} d\Omega + \int \sqrt{-g}.\delta(g_{ij}).R^{ij} d\Omega + \int \sqrt{-g}.g_{ij}.\delta(R^{ij}) d\Omega ]</math>
On a <math>\delta (\sqrt{-g}) = \frac{-\delta g}{2\sqrt{-g}} = - \frac{1}{2\sqrt{-g}}g.g_{ik}. \delta g^{ik} = - \frac{1}{2} \sqrt{-g}.g^{ik} \delta g_{ik}</math> car <math> g^{ik}.g_{ik} = 4 \to \delta (g^{ik}).g_{ik} = - g^{ik}.\delta (g_{ik})</math>
Pour la 1{{re}} intégrale, on a <math> \int \delta (\sqrt{-g})g_{ij}R^{ij} d\Omega = \int \delta (\sqrt{-g})R d\Omega = - \frac{1}{2} \int g^{ij}.R. \sqrt{-g} .\delta g_{ij} d\Omega </math>
La 2{{e}} égalité est laissée inchangée.
Pour la 3{{e}} intégrale, pour simplifier les calculs, on se place dans un référentiel en apesanteur et on a donc <math>\ R^{ij} = \part^l \Gamma_l^{ij} - \part^i \Gamma_l^{lj}</math>. (Mais en général <math>\part^l \Gamma_l^{ij} \ne D^l \Gamma_l^{ij}</math> car le symbole de [[Elwin Bruno Christoffel|Christoffel]] n'est pas un tenseur).
D'où <math> \delta \left( R^{ij} \right) = \delta \part^l \Gamma_l^{ij} - \delta \part^i \Gamma_l^{lj} = \part^l \delta\Gamma_l^{ij} - \part^i \delta\Gamma_l^{lj} </math> en supposant que la variation des <math>\ g^{ij}</math> laisse le référentiel en apesanteur en ce point, ce qui laisse encore une infinité de variations possibles pour les <math>\ g^{ij}</math>.
Dans n'importe quel référentiel, <math> \delta \Gamma_l^{ij} = \Gamma _l^{ij} - \left( \Gamma_l^{ij} \right) '</math> où le symbole <math>\left( \Gamma_l^{ij} \right) '</math> est le symbole de Christoffel au même point que <math>\ \Gamma_l^{ij}</math> mais avec des termes <math>\ g^{ij}</math> modifiés
on a <math> \Gamma_k^{ij} .V_j = \part^jV_i \to \delta \Gamma_k^{ij}.V_j = \Gamma _l^{ij}.V_j - \left( \Gamma_l^{ij} \right) '.V_j = \part^jV_i - \left( \part^jV_i \right)'</math> ce qui est une différence entre deux tenseurs définis au même point, donc <math> \delta \Gamma_k^{ij}</math> est un tenseur (contrairement au symbole de Christoffel).
Et pour ce tenseur, dans le référentiel en apesanteur (et laissé comme tel, au point considéré, par la variation des <math>\ g^{ij}</math>), <math>\part^l \delta \Gamma_k^{ij} = D^l \delta \Gamma_k^{ij}</math>, d'où <math> \delta \left( R^{ij} \right) = \part^l \delta\Gamma_l^{ij} - \part^i \delta\Gamma_l^{lj} = D^l \delta\Gamma_l^{ij} - D^i \delta\Gamma_l^{lj}</math>
<math>\sqrt{-g}.g_{ij}.\delta R^{ij} = \sqrt{-g}.[g_{ij}.D^l \delta\Gamma_l^{ij} - g_{ij}.D^i \delta\Gamma_l^{lj}] = \sqrt{-g}.[D^l \left( g_{ij}.\delta\Gamma_l^{ij} \right) - D^i \left( g_{ij}.\delta\Gamma_l^{lj} \right) ]</math>
car <math>~\quad D_kg_{ij} = 0~</math> et aussi <math>~\quad D^kg_{ij} = 0~</math>
d'où <math>~\quad \sqrt{-g}.g_{ij}.\delta R^{ij} = \sqrt{-g}.D^l \left( g_{ij}.\delta\Gamma_l^{ij} - g_{lj}.\delta\Gamma_i^{ij} \right) = \sqrt{-g}.D^lA_l</math>.
D'où, en utilisant le théorème d'Ostrogradski, <math>\int \sqrt{-g}.g_{ij}.\delta R^{ij} d\Omega = \int \sqrt{-g}.D^lA_l d\Omega = \int \sqrt{-g}.A_l.dS^l = 0</math>
La nullité de la dernière intégrale est dûe au fait qu'elle est calculée sur l'hypersurface délimitant le volume d'intégration et au fait que les variations des <math>\ g^{ij}</math> sont nulles sur la frontière d'intégration.
On obtient : <math>\delta S_g = K. \int \left( R^{ij} - \frac{1}{2}g^{ij}R \right) \sqrt{-g}.\delta g_{ij} d\Omega</math>
Le principe de moindre action disant que <math>\ \delta S_g = 0 </math> et les variations <math>\ \delta g_{ij}</math> étant quelconques, on obtient <math>\ R^{ij} - \frac{1}{2}g^{ij}R = 0</math> , ce que l'on écrit (et démontre) souvent en baissant les indices.
</center>
}}
Les équations déduites sont :
<center>
<math>\ R_{ij} - \frac{1}{2}g_{ij}R = 0</math>
</center>
En faisant la « contraction » <math>\ g^{ij}R_{ij} - \frac{1}{2}g^{ij}.g_{ij}R = 0</math>, on obtient <math>\ R = 0</math> , ce qui ne signifie pas que l'espace est plat, mais plutôt qu'il s'agit d'une [[surface minimale]] à quatre dimensions, tendue entre les différentes masses qui y évoluent.
Les équations d'Einstein dans le cas extérieur sont donc :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\ R_{ij} = 0</math>
}}
====Les équations d'Einstein du champ de gravitation dans le ''cas intérieur''====
Le deuxième cas des équations du champ est le cas où il y a de la matière (localement) : on parle du « cas intérieur », c'est-à-dire « dans la matière ».
Dans ce cas, l'action est composée de l'action du champ gravitationnel <math>\ S_g = K. \int \sqrt{-g}.R.d\Omega</math> et de l'action de la matière, en y incluant le champ électromagnétique, que l'on écrit <math>\ S_m = \frac{1}{c} \int\sqrt{-g}. \Lambda_m d\Omega</math>.
{{Boîte déroulante|titre=''Démonstration des équations d'Einstein dans le cas intérieur''|contenu=
<center>
En utilisant la même méthode variationnelle, en sachant que <math>\ \delta \part =\part \delta</math>, en utilisant l'intégration par parties, et le théorème d'Ostrogradski qui permet d'écrire dans un référentiel en apesanteur <math> \int \part_l [ \frac{\part \left( \sqrt{-g} \Lambda_m \right)}{\part \left( \part_l g_{ik} \right)} ] d \Omega = \int \part_l [ \sqrt{-g} \frac{\part \left( \Lambda_m \right)}{\part \left( \part_l g_{ik} \right)} ] d \Omega = \int \sqrt{-g} \frac{\part \left( \Lambda_m \right)}{\part \left( \part_l g_{ik} \right)} dS_l = 0 </math>
En définissant le tenseur impulsion-énergie <math>\ T^{ij}</math> par l'égalité <math> \frac{1}{2} T^{ij} \sqrt{-g} = \frac{\part \left( \sqrt{-g} \Lambda_m \right)}{\part g_{ik}} - \part_l [ \frac{\part \left( \sqrt{-g} \Lambda_m \right)}{\part \left( \part_l g_{ik} \right)} ]</math>
On obtient : <math>\delta S_m = \frac{1}{c} \int \frac{1}{2} T^{ij} \sqrt{-g} \delta g_{ij}d\Omega</math>
D'où, en posant <math>\chi = - \frac{1}{2c.K}</math> , et on conclut de la même manière que dans le ''cas extérieur''.
</center>
}}
Les équations déduites sont :
{{Bloc emphase|fond=#fff8ff;border:4px double #888|
<math>\ R_{ij} - \frac{1}{2}g_{ij}R = \chi T_{ij}</math>
}}
Avec la contraction similaire au ''cas extérieur'', sachant que <math>\ g_{ij}g^{ij} = 4</math> et en posant <math>\ T = g^{ij}T_{ij} </math>, on a <math>\ R = - \chi T</math>. La [[courbure]] principale est donc proportionnelle à la ''densité d'énergie totale'' <math>\ T = g^{ij}T_{ij}</math> (ou ''trace du tenseur <math>\ T_{ij}</math>'').
== Notes et références ==
{{Références|colonnes=2}}
== Annexes ==
=== Bibliographie ===
* {{Landau}}
* Richard P. Feynman ; ''Le cours de physique de Feynman - Electromagnétisme (I)'', chapitre 19, InterEditions (1979), ISBN 2-7296-0028-0. Réédité par Dunod (2000), ISBN 2-1000-4861-9
* Jean-Claude Boudenot ; ''Électromagnétisme et gravitation relativistes'', ellipse (1989), ISBN 2729889361
* Jean-Louis Basdevant ; ''Principes variationnels & dynamique'', Vuibert (2005), ISBN 2711771725.
* Florence Martin-Robine ; ''Histoire du principe de moindre action'', Vuibert (2006), ISBN 2711771512.
* Edgard Elbaz ; ''Relativité générale et gravitation'', ellipse (1986).
=== Articles connexes ===
* [[Mécanique analytique]]
* [[Lagrangien]]
* [[Mécanique hamiltonienne]]
* [[Espace de Minkowski]]
* [[Relativité restreinte]]
* [[Relativité générale]]
* [[Principe de Fermat]], intéressant parallèle avec l'optique
=== Liens externes ===
* C. Cohen-Tannoudji, [http://www.phys.ens.fr/cours/notes-de-cours/cct-dea/ ''Forme lagrangienne de la mécanique quantique''], cours de 1966 à l'ENS
{{Portail|physique}}
[[Catégorie:Principe physique]]
[[Catégorie:Physique théorique]]
[[br:Pennreolenn an aktadur bihanañ]]
[[de:Prinzip der kleinsten Wirkung]]
[[en:Principle of least action]]
[[es:Principio de mínima acción]]
[[he:עקרון הפעולה המינימלית]]
[[ja:最小作用の原理]]
[[ru:Принцип наименьшего действия]]
[[sl:načelo najmanjše akcije]]