title: Principe de moindre action et relativité restreinte
text: # Principe de moindre action et relativité restreinte
  
  En relativité restreinte, le principe de moindre action pour un point matériel donne des équations d'Euler-Lagrange très semblables à celles de la mécanique classique, mais le lagrangien n'est plus égal à la différence entre l'énergie cinétique et l'énergie potentielle. En fait, le principe de moindre action se base dans ce cas seulement 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é (on peut généraliser à un système de points), déterminées à partir de principes généraux, tels que par exemple :
  
  * Comme la trajectoire dans l'espace-temps ne dépend pas du référentiel par rapport auquel on l'observe, l'action qui la détermine, ainsi que les fonctions qui composent l'action, sont invariantes par changement de référentiel.
  * L'indépendance de plusieurs 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 $~-mc\tau ~~$ où $~\tau ~$ est « temps propre » du trajet (temps s'écoulant dans le référentiel du corps au cours du trajet). Minimiser l'action revient à maximiser le « temps propre », du fait du signe $-$, que la masse $~m~$ est positive et constante, et de la constance de la vitesse de la lumière $~c~$. Cela est aussi lié à la « longueur généralisée » de la trajectoire, mesurée par la métrique de l'espace-temps de Minkowski.
  
  Un champ électromagnétique amène des différences de parcours entre les corps, suivant leurs charges et leurs répartitions.
  
  Et comme en physique classique, toutes les équations peuvent être obtenues sans le principe de moindre action.
  
  ## Étude d'un point matériel
  
  ### Formulation analogue à celle de la mécanique non relativiste
  
  **Point matériel isolé (corps libre)**
  
  * 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 $~x;y;z~$ (repère généralement supposé cartésien orthonormé) et sa coordonnée de temps $~t~$, subissant toutes une modification en cas de changement de référentiel galiléen.
  
      Les coordonnées d'espace d'un point matériel $M$ sont ici souvent notées $x_{i}$ (où $i=1;2;3$), mais la comparaison avec la quadri-écriture nécessite une attention particulière car cette dernière utilise des indices écrits en haut (à ne pas confondre avec des exposants) puisque ce sont des quantités contravariantes.
      On étudie ici le principe de moindre action sous la forme « classique » : $~S=\int L(x_{i}(t);v_{i}(t),t)\ dt$ , avec $~v_{i}={\frac {dx_{i}}{dt}}~$ (où $~i=1;2;3$).
  * De même qu'en mécanique non relativiste, l'homogénéité et l'isotropie de l'espace imposent que le lagrangien ne dépende que de la norme de la vitesse (ou son carré) : $~L=L(\left\Vert {\vec {v}}\right\Vert ^{2})$.
  
      L'invariance par changement de référentiel galiléen impose que l'action soit un scalaire vis-à-vis de la transformation de Lorentz ; ceci suggère une expression de la forme : $~S\propto \int ds$ avec l'élément d'espace-temps $~ds={\sqrt {(c.dt)^{2}-(dx)^{2}-(dy)^{2}-(dz)^{2}}}$, c'est-à-dire : $~S\propto \int {{\sqrt {1-\beta ^{2}}}\ dt}~$ avec $~\beta ={\frac {v}{c}}$.
      ♦ remarque : le minimum de l'action correspond ainsi à un mouvement dont la durée propre ($~d\tau ~$ telle que $~ds=c\ d\tau$) pour le point étudié est extrémale.
  * Pour les faibles vitesses : ${\sqrt {1-\beta ^{2}}}\approx 1-{\frac {v^{2}}{2c^{2}}}$ devrait logiquement redonner la limite non relativiste : $L={\frac {1}{2}}mv^{2}$. Or, les lois qui s'en déduisent (équations d'Euler-Lagrange) sont inchangées si on ajoute au lagrangien la dérivée totale d'une fonction ne dépendant que de la position et du temps, donc en particulier pour une constante.
  
      On peut donc considérer :
  
  $S=-mc\int ds$
  
  correspondant au lagrangien :
  
  $L=-mc^{2}{\sqrt {1-\beta ^{2}}}$
  
      ♦ remarque : par construction, ceci peut décrire une particule matérielle, mais la description d'un photon nécessite une étude plus approfondie.
  
  **Impulsion et énergie**
  
  * Par définition l'impulsion $~{\vec {p}}~$ relativiste correspond à : $~p_{i}={\frac {\partial L}{\partial v_{i}}}~$ (où $~i=1;2;3$) ; ainsi :
  
  ${\vec {p}}=\gamma m{\vec {v}}$
  
  $\gamma ={\frac {1}{\sqrt {1-\beta ^{2}}}}$
  
  * De même qu'en mécanique non relativiste, le hamiltonien peut être défini par : $~H=\sum _{i}\ p_{i}v_{i}-L={\vec {p}}.{\vec {v}}-L=\gamma mc^{2}={\sqrt {m^{2}c^{4}+c^{2}p^{2}}}~$ (il doit être exprimé en fonction de l'impulsion).
  
      Il correspond à une énergie, constante pour tout système dont le hamiltonien (le lagrangien) ne dépend pas explicitement du temps :
  
  $\ E=\gamma mc^{2}={\sqrt {m^{2}c^{4}+c^{2}p^{2}}}$
  
      En particulier, pour $~v=0$, l'énergie au repos est $~E=mc^{2}$.
      ♦ remarque : par différence, l'énergie cinétique correspond à : $~E_{c}=(\gamma -1)mc^{2}~$ ; $~E_{c}\approx {\frac {1}{2}}mv^{2}~$ à l'approximation des faibles vitesses ; mais le lagrangien relativiste d'un point matériel isolé en diffère d'une constante.
  
  **Interaction avec un champ électromagnétique**
  
  * L'expérience montre que les effets d'un champ électromagnétique peuvent être exprimés à l'aide d'un potentiel scalaire $\phi$ et d'un potentiel vecteur ${\vec {A}}$.
  
      Ceci correspond à un quadri-potentiel $~{\tilde {A}}=({\frac {\phi }{c}},{\vec {A}})~$ utilisé ici uniquement pour expliquer la forme du terme d'interaction dans le lagrangien. Ce terme doit dépendre de la vitesse (les effets magnétiques en dépendent) et être invariant de Lorentz.
      Pour généraliser un déplacement $d{\vec {M}}$ de coordonnées $dx_{i}$ (avec $i=1;2;3$) on peut utiliser le quadri-déplacement $d{\tilde {M}}=(c\ dt;d{\vec {M}})$.
      On peut alors proposer de décrire l'interaction d'une particule de charge $e$ par un terme d'interaction proportionnel à : $e\ {\tilde {A}}\cdot d{\tilde {M}}=e\ \phi (x_{i},t)dt-e\ {\vec {A}}(x_{i},t)\cdot {\vec {v}}\ dt$.
      Afin de retrouver dans la limite des faibles vitesses et sans effet magnétique : $L=E_{c}-U$ avec $U=e\ \phi$, on considère ici le lagrangien :
  
  $L=L(x_{i},v_{i},t)=-m.c^{2}.{\sqrt {1-\beta ^{2}}}-e\phi (x_{i},t)+e\ {\vec {A}}(x_{i},t)\cdot {\vec {v}}$
  
  * On en déduit une impulsion généralisée : $P_{i}={\frac {\partial L}{\partial v_{i}}}=p_{i}+e\ A_{i}$, soit ${\vec {P}}={\vec {p}}+e\ {\vec {A}}$ en conservant la notation ${\vec {p}}=\gamma m{\vec {v}}$.
  
      ♦ remarque : on prendra donc soin de distinguer ${\vec {p}}$ et ${\vec {P}}$.
  * De même que pour un corps libre, le hamiltonien peut être défini par : $H=\sum _{i}\ P_{i}v_{i}-L={\vec {P}}.{\vec {v}}-L=\gamma mc^{2}+e\ \phi ={\sqrt {m^{2}c^{4}+c^{2}({\vec {P}}-e{\vec {A}})^{2}}}+e\ \phi$ (il doit être exprimé en fonction de l'impulsion généralisée).
  
      Il correspond à une énergie, constante pour tout système dont le hamiltonien (le lagrangien) ne dépend pas explicitement du temps :
  
  $E_{tot}=E+e~\phi =\gamma mc^{2}+e~\phi ={\sqrt {m^{2}c^{4}+c^{2}p^{2}}}+e\ \phi$
  
      Toutes les approximations aux petites vitesses devant $c$ redonnent les résultats classiques.
  
  **Équations du mouvement**
  
  * Le mouvement est décrit par les équations d'Euler-Lagrange : ${\frac {dP_{i}}{dt}}={\frac {\partial L}{\partial x_{i}}}=e\ {\frac {\partial ({\vec {A}}\cdot {\vec {v}})}{\partial x_{i}}}-e\ {\frac {\partial \phi }{\partial x_{i}}}$.
  
      Avec la notation gradient, ceci peut s'écrire sous la forme : ${\frac {d{\vec {P}}}{dt}}=e\ {\vec {\nabla }}({\vec {A}}\cdot {\vec {v}})-e\ {\vec {\nabla }}\phi$, avec par ailleurs (puisque $x_{i}$ et $v_{i}$ sont des variables indépendantes) : ${\vec {\nabla }}({\vec {A}}\cdot {\vec {v}})=({\vec {v}}\cdot {\vec {\nabla }}){\vec {A}}+{\vec {v}}\times ({\vec {\nabla }}\times {\vec {A}})$ avec la notation rotationnel.
      En outre, puisque la particule se déplace, la variation de ${\vec {A}}$ en un point fixe donné est : ${\frac {\partial {\vec {A}}}{\partial t}}={\frac {d{\vec {A}}}{dt}}-({\vec {v}}\cdot {\vec {\nabla }}){\vec {A}}$ ; ceci donne finalement :
  
  ${\frac {d{\vec {p}}}{dt}}=e\ {\vec {E}}+e\ {\vec {v}}\times {\vec {B}}$ avec ${\vec {E}}=-{\vec {\nabla }}\phi -{\frac {\partial {\vec {A}}}{\partial t}}$ et ${\vec {B}}={\vec {\nabla }}\times {\vec {A}}$
  
      ♦ remarque : certains calculs sont plus pratiques en fonction de ${\vec {v}}$ mais d'autres sont plus simples en fonction de ${\vec {p}}$.
      ♦ remarque : attention à ne pas confondre l'énergie $E$ et le champ électrique ${\vec {E}}$.
  
  **Relation entre l'action et l'énergie impulsion**
  
  * De façon générale : $~dS=L~dt=\sum _{i}\ P_{i}~dx_{i}-H~dt~$ (où $~i=1;2;3$).
  
  ### Formulations quadrivectorielles
  
  **Notations**
  
  * Avec les notations quadri-vectorielles, un point $M$ a des coordonnées d'espace-temps contravariantes $\ (x^{\alpha })=(ct;x;y;z)$ avec $\ \alpha =(0;1;2;3)$. On note les indices quadri-vectoriels par des symboles grecs ; on utilise des indices latins pour désigner la partie spatiale des quadrivecteurs.
  
      Nous adopterons la convention de sommation d'Einstein : sommation sous-entendue pour tout indice répété en haut et en bas : $U^{\alpha }V_{\alpha }=\Sigma _{\alpha =0}^{3}U^{\alpha }V_{\alpha }$.
      On note $\eta _{\alpha \beta }$ la matrice de la métrique permettant d'abaisser les indices (les termes non diagonaux sont nuls) : $(\eta _{00}=1;\eta _{ii}=-1)$. Ainsi pour un quadri-vecteur $\ (U^{0};U^{1};U^{2};U^{3})\$ on obtient : $U_{\alpha }=\eta _{\alpha \beta }U^{\beta }\$ c'est-à-dire : $\ U_{0}=U^{0}$ et $\ U_{i}=-U^{i}$. On note $\eta ^{\alpha \beta }$ la matrice inverse permettant d'élever les indices (les termes non diagonaux sont nuls) : $(\eta ^{00}=1;\eta ^{ii}=-1)$.
      En particulier pour deux quadri-vecteurs $\ (U^{0};U^{1};U^{2};U^{3})\$ et $\ (V^{0};V^{1};V^{2};V^{3})$, on définit le produit scalaire $\ {\tilde {U}}\cdot {\tilde {V}}=\eta _{\alpha \beta }U^{\alpha }V^{\beta }=U^{\alpha }V_{\alpha }=V^{0}.U^{0}-V^{1}.U^{1}-V^{2}.U^{2}-V^{3}.U^{3}$ et la (pseudo)-norme : $\ \Vert {\tilde {U}}\Vert ^{2}={\tilde {U}}\cdot {\tilde {U}}=(U^{0})^{2}-(U^{1})^{2}-(U^{2})^{2}-(U^{3})^{2}$.
      De manière similaire, on écrira : $\partial ^{\alpha }={\frac {\partial }{\partial x_{\alpha }}}$ et $\partial _{\alpha }={\frac {\partial }{\partial x^{\alpha }}}$.
  * Pour décrire le mouvement d'un point matériel, on peut continuer à paramétrer la quadri-position en fonction du temps : $\ x^{\alpha }(t)$ (dans ce cas $\ x^{0}(t)$ est « inutile ») ; mais on peut aussi utiliser un autre paramètre, en particulier si on souhaite une formulation invariante de Lorentz.
  
      On peut choisir de paramétrer par le temps propre $\tau$ ou selon l'élément d'espace temps tel que $ds=c\ d\tau$. Certains calculs formels sont toutefois moins évidents dans ce cas : il est alors possible de choisir un paramètre $\sigma$ « quelconque », quitte à imposer ensuite $\sigma =\tau$ ou $\sigma =s$.
  * Il existe toutefois, entre les notations non relativistes et les notations quadri-vectorielles, des différences de conventions de notations pouvant causer des confusions de signes. Il semble ne pas exister de consensus dans la communauté scientifique pour privilégier un choix particulier de notations (principalement parce que cela n'a pas d'influence sur les lois physiques qui s'en déduisent). Les différentes approches sont précisées dans ce qui suit.
  
  **Méthode se cantonnant aux propriétés de base**
  
  * De nombreux physiciens définissent l'énergie-impulsion indépendamment du principe de moindre action et se limitent à utiliser les propriétés de base de ce dernier pour montrer comment les deux sont liés (par exemple S. Weinberg ; Gravitation and cosmology : « Il apparaît ... que les $x^{\alpha }$ obéissent à ces équations du mouvement et nous concluons donc que cela qualifie l'expression utilisée pour l'action comme appropriée pour ce système » ).
  
      ♦ remarque : le principe de moindre action n'est en rien indispensable à l'établissement des lois physiques (il n'est qu'une façon, parfois pratique, de les décrire).
  * Compte tenu des détails indiqués dans les autres parties, on traite ici directement le cas avec champ électromagnétique (considéré comme imposé par l'extérieur).
  * Avec $V^{\alpha }={\frac {dx^{\alpha }}{d\tau }}$, l'action peut s'écrire : $S=-mc\int ds-e\int A_{\alpha }V^{\alpha }d\tau =-mc\int {\sqrt {\eta _{\alpha \beta }dx^{\alpha }dx^{\beta }}}-e\int A_{\alpha }dx^{\alpha }$.
  * On peut utiliser : $\delta S=-mc\int {\frac {\eta _{\alpha \beta }dx^{\alpha }\delta (dx^{\beta })}{\sqrt {\eta _{\alpha \beta }dx^{\alpha }dx^{\beta }}}}~-~e\int A_{\alpha }~\delta (dx^{\alpha })=-m\int V_{\alpha }~d(\delta x^{\alpha })~-~e\int A_{\alpha }~d(\delta x^{\alpha })$.
  
      L'intégration par parties donne : $\delta S=-\left[(mV_{\alpha }+eA_{\alpha })~\delta x^{\alpha }\right]+\int {\frac {d(mV_{\alpha }+eA_{\alpha })}{d\sigma }}\delta x^{\alpha }d\sigma$.
      Pour des $\delta x^{\alpha }$ nuls aux extrémités du mouvement : $\delta S=\int \eta _{\alpha \beta }{\frac {d(mV^{\alpha }+eA^{\alpha })}{d\sigma }}\delta x^{\beta }d\sigma$.
  * Pour des $\delta x^{\alpha }$ quelconques durant le mouvement, la condition d'extremum $~\delta S=0~$ impose les équations du mouvement : ${\frac {d(mV^{\alpha }+eA^{\alpha })}{d\sigma }}=0$.
  
      On obtient en particulier une quadri-vitesse constante pour un point matériel isolé : ${\frac {dV^{\alpha }}{d\sigma }}=0$.
      ♦ remarque : on paramètre généralement par $\tau$ mais on peut ici en principe utiliser un paramètre $\sigma$ « quelconque ».
      ♦ remarque : ici encore, la description d'un photon nécessite d'utiliser un autre paramètre car $~d\tau =0~$ donc $V^{\alpha }$ n'est pas défini ; en outre, la méthode doit être adaptée pour éviter qu'un terme (quel qu'il soit) proportionnel à $d\tau$ apparaisse au dénominateur.
  
  **Formulation la plus simplement ressemblante à celle sans quadri-écriture pour un corps libre**
  
  * Ce qui précède suggère qu'on peut raisonner en paramétrant par $\tau$ (ou par $s$ qui lui est proportionnel) et que le principe de Hamilton peut s'écrire à l'aide d'un lagrangien $L(x^{\alpha }(\tau );V^{\alpha }(\tau );\tau )$ avec $V^{\alpha }={\frac {dx^{\alpha }}{d\tau }}$.
  
      De nombreux physiciens utilisent alors : $\ L=-mc\ {\sqrt {\eta _{\alpha \beta }V^{\alpha }V^{\beta }}}=-mc^{2}$. Avec $\ ds={\sqrt {\eta _{\alpha \beta }dx^{\alpha }dx^{\beta }}}$, l'action minimisée entre deux points de l'espace-temps $\ S=\int _{A}^{B}L\ d\tau =-mc\int _{A}^{B}ds$ montre que le chemin suivi par la particule pour aller du point A au point B est celui qui maximise le temps propre (le terme négatif $-mc$ transforme la minimisation de $\ S=-mc\int _{A}^{B}ds$ en maximisation de $\int _{A}^{B}ds$).
      Cette écriture peut sembler formelle dans la mesure où ${\sqrt {\eta _{\alpha \beta }V^{\alpha }V^{\beta }}}=c$ semble indiquer qu'on cherche l'extrémum d'une « constante » (fonction indépendante de la position et de la vitesse). En fait il faut distinguer la quadri-distance $ds$ et le paramètre $\tau$ (ou $s$) utilisé.
  * Les quantités candidates pour définir la quadri-impulsion « conjuguée » $P_{\alpha }$ semblent alors être (mais cela reste à préciser) : ${\frac {\partial L}{\partial V^{\alpha }}}=-mc{\frac {V_{\alpha }}{\sqrt {\eta _{\alpha \beta }V^{\alpha }V^{\beta }}}}=-mV_{\alpha }$.
  
      L'utilisation d'un autre paramètre aboutit au même résultat.
  * Pour contourner l'ambiguïté de l'écriture formelle avec ${\sqrt {\eta _{\alpha \beta }V^{\alpha }V^{\beta }}}=c$ (qui n'est qu'apparente), on peut choisir un paramètre $\sigma$ « arbitraire », en notant $L(x^{\alpha }(\sigma );{\underline {V}}^{\alpha }(\sigma );\sigma )$ avec ${\underline {V}}^{\alpha }={\frac {dx^{\alpha }}{d\sigma }}$, quitte à considérer ensuite la limite $\sigma \rightarrow \tau$ par continuité.
  * Les quantités candidates pour définir la quadri-impulsion « conjuguée » $P_{\alpha }$ semblent alors être (mais cela reste à préciser) : ${\frac {\partial L}{\partial {\underline {V}}^{\alpha }}}=-mc{\frac {{\underline {V}}_{\alpha }}{\sqrt {\eta _{\alpha \beta }{\underline {V}}^{\alpha }{\underline {V}}^{\beta }}}}=-m{\frac {dx_{\alpha }}{d\sigma }}{\frac {d\sigma }{d\tau }}=-mV_{\alpha }$. Le fait de passer par l'intermédiaire de $\sigma$ ne change rien (on peut paramétrer directement par $\tau$).
  * Or la comparaison avec la mécanique non relativiste (étayée par la description sans quadri-écriture) indique une énergie-impulsion $~p^{\alpha }=mV^{\alpha }~$ ; le signe ne correspond pas. Raisonner avec une impulsion généralisée $~P^{\alpha }=-p^{\alpha }~$ (ou l'analogue pour les composantes covariantes ; la métrique $\eta _{\alpha \beta }$ n'influence pas les équations) semble au moins contre-intuitif.
  * Cette étrangeté de signe provient de la façon un peu trop simple avec laquelle l'impulsion est définie en mécanique classique, puis généralisée à la mécanique relativiste quand on n'utilise pas les notations quadrivectorielles. Quand on écrit dans ce cadre (sans indices co/contravariants) : $~P_{i}$ “=” ${\frac {\partial L}{\partial V_{i}}}~$ (avec $~\alpha =i=1;2;3~$), on l'interprète (en précisant les indices) comme : $~P^{i}$ “=” ${\frac {\partial L}{\partial V^{i}}}~$.
  
      Or ceci devient ici incohérent, puisque dans un changement de coordonnées (ici l'exemple typique est la transformation de Lorentz), la grandeur contrevariante $~P^{i}$ subit la transformation inverse de celle à appliquer à la quantité covariante ${\frac {\partial L}{\partial V^{i}}}~$.
      On devrait donc maintenant définir ces "impulsions conjuguées" par $~{\mathcal {P}}_{i}=-{\frac {\partial L}{\partial V^{i}}}~$, afin que cela corresponde à $~{\mathcal {P}}^{i}=\eta ^{ii}{\mathcal {P}}_{i}=p^{i}=mV^{i}~$ (sans sommation sur les indices, avec $~\eta ^{ii}=-1~$).
  * Puisque $~\eta ^{00}=+1~$, pour justifier que ce signe s'applique aussi à $~{\mathcal {P}}_{0}=-{\frac {\partial L}{\partial V^{0}}}~$ (quantité qui était décrite à l'aide du hamiltonien $H$ dans la formulation non quadrivectorielle), on peut considérer la relation entre l'action et l'énergie-impulsion, avec l'écriture quadrivectorielle : $~dS=L~d\tau =-{\mathcal {P}}_{\alpha }~dx^{\alpha }~$.
  
      Ceci montre que de façon symétrique $~{\mathcal {P}}_{\alpha }=-\partial _{\alpha }S~$ (alors que la notation non quadrivectorielle distinguait les signes : $~P_{i}=\partial _{i}S~$ et $~{\frac {H}{c}}=-\partial _{0}S~$).
      ♦ remarque : ceci montre aussi que le changement de signe de $~\eta ^{\alpha \beta }~$ (signature {-+++} au lieu de {+---}) ne ferait que déplacer l'étrangeté de signe sur la coordonnée temporelle au lieu des coordonnées spatiales.
  * L'effet de la métrique peut aussi apparaitre en relativité restreinte avec des coordonnées non cartésiennes. Ainsi une transformation de Lorentz associée à un mouvement d'axe (Ox) peut par exemple s'exprimer avec des coordonnées cylindriques $~(ct;x;r;\theta )~$. On utilise dans ce cas $~ds^{2}=c^{2}\ dt^{2}-dx^{2}-dr^{2}-r^{2}\ d\theta ^{2}~$ et la métrique $\eta _{\alpha \beta }$ est remplacée par $g_{\alpha \beta }$ avec $~g_{33}=-r^{2}~$. La composante angulaire donne $~V^{3}={\frac {d\theta }{d\tau }}~$ et l'impulsion conjuguée est un moment cinétique $~P_{3}=-r^{2}{\frac {d\theta }{d\tau }}~$ (en géométrie riemannienne la méthode de Lagrange donne l'impulsion généralisée associée à la coordonnée choisie).
  
      Ceci montre qu'à partir d'un certain niveau de généralisation, les expressions des "impulsions" généralisées ne correspondent pas forcément à ce qui pourrait être simplement attendu. Les signes contre-intuitifs constatés ici n'ont donc rien de rédhibitoire.
      ♦ remarque : la relativité générale utilise d'ailleurs la métrique $~g_{\alpha \beta }~$ pour décrire la gravitation. Il est donc judicieux de commencer par préciser l'écriture quadrivectorielle correspondante en relativité restreinte avant de généraliser davantage.
      ♦ remarque : pour diminuer l'impact des ambiguïtés de signe liées à la métrique, une autre méthode (L. Landau et E. Lifchitz) peut consister à utiliser une première coordonnée imaginaire au lieu des notations co/contra-variantes $~x_{0}=ict~$ (comme en notations non relativistes, il n'y a alors pas d'indices en haut) et de même pour tous les quadrivecteurs, mais ceci ne s'adapte pas à la relativité générale (ni même aux coordonnées non cartésiennes).
  * Rien toutefois n'interdit de raisonner avec les quantités $~P_{\alpha }={\frac {\partial L}{\partial V^{\alpha }}}=-mV_{\alpha }~$, même s'il peut devenir ambigu de continuer à les nommer "impulsions" (mais cette habitude est ancrée depuis longtemps).
  
      C'est d'ailleurs le choix de certains physiciens comme Cornelius Lanczos (The Variational Principles of Mechanics).
  * Puisque le lagrangien d'un point matériel isolé est indépendant de la position, les équations d'Euler-Lagrange donnent alors pour lois du mouvement : ${\frac {dP_{\alpha }}{d\tau }}=0$ (constantes du mouvement).
  * Le carré de la "norme" de cette quadri-impulsion est en outre : $P^{\alpha }P_{\alpha }=p^{\alpha }p_{\alpha }=m^{2}~V^{\alpha }V_{\alpha }=m^{2}c^{2}$, donc conforme à la relation usuelle : $E^{2}-c^{2}p^{2}=m^{2}c^{4}$.
  * Le hamiltonien correspondant (à ne pas confondre avec celui de la formulation non quadrivectorielle) peut être défini par : $~H=P_{\alpha }V^{\alpha }-L=-mc^{2}+mc^{2}=0~$. Écrit ainsi il ne décrit pas une énergie, mais il est vrai que ce rôle est ici dédié à $P^{0}$.
  * D'un premier point de vue, on peut utiliser le fait que $L=-mc^{2}$ pour écrire formellement : $H={\frac {-1}{m}}\ P_{\alpha }P^{\alpha }+mc^{2}$ (il doit être exprimé en fonction de l'impulsion).
  
      On obtient alors : ${\frac {\partial H}{\partial P_{\alpha }}}={\frac {-2\ P^{\alpha }}{m}}=2V^{\alpha }$ ; ceci n'est pas conforme à la description de Hamilton.
      ♦ remarque : on obtient ici de façon un peu inattendue une expression sans radical.
  * Mais pour contourner l'ambiguïté de l'écriture formelle, on peut choisir un paramètre $\sigma$ « arbitraire », en notant $L(x^{\alpha }(\sigma );{\underline {V}}^{\alpha }(\sigma );\sigma )$ avec ${\underline {V}}^{\alpha }={\frac {dx^{\alpha }}{d\sigma }}$, puis considérer ensuite la limite $\sigma \rightarrow \tau$ par continuité.
  
      Le hamiltonien peut alors s'écrire : $H=P_{\alpha }{\underline {V}}^{\alpha }-L=-mc{\frac {\eta _{\alpha \beta }{\underline {V}}^{\alpha }{\underline {V}}^{\beta }}{\sqrt {\eta _{\alpha \beta }{\underline {V}}^{\alpha }{\underline {V}}^{\beta }}}}-L=-mc{\sqrt {\eta _{\alpha \beta }{\underline {V}}^{\alpha }{\underline {V}}^{\beta }}}-L=0$.
      En choisissant maintenant une paramétrisation par $\tau$ et en utilisant le fait qu'alors $L=-mc^{2}$ on peut écrire formellement : $H=-c\ {\sqrt {P_{\alpha }P^{\alpha }}}+mc^{2}$ (il doit être exprimé en fonction de l'impulsion).
      On obtient alors : ${\frac {\partial H}{\partial P_{\alpha }}}={\frac {P^{\alpha }}{m}}=V^{\alpha }$ ; ceci est conforme à la description de Hamilton.
      ♦ remarque : on obtient ici de façon un peu plus attendue une expression avec radical.
  * Tout ceci ressemble toutefois surtout à du bricolage formel, pouvant difficilement servir de base de raisonnement convaincante.
  
  **Formulation la plus simplement ressemblante à celle sans quadri-écriture dans un champ électromagnétique**
  
  * Le champ électromagnétique peut être décrit à l'aide d'un quadri-vecteur, appelé quadri-potentiel électromagnétique, $A^{\alpha }$ dont l'interaction avec la particule de charge $e$ se manifeste sous forme lagrangienne par un terme $-e~A_{\alpha }dx^{\alpha }$.
  
      En paramétrant par $\tau$ le lagrangien relativiste peut donc s'écrire :
  
  $L=-mc~{\sqrt {V_{\alpha }V^{\alpha }}}-e~A_{\alpha }V^{\alpha }$
  
  * En posant $F_{\alpha \beta }=\partial _{\alpha }A_{\beta }-\partial _{\beta }A_{\alpha }$ tenseur champ électromagnétique (et $F_{~\beta }^{\alpha }=\eta ^{\alpha \gamma }~F_{\gamma \beta }$), les équations d'Euler-Lagrange donnent les équations du mouvement de la particule :
  
  $m~{\frac {dV^{\alpha }}{d\tau }}=e~F_{~\beta }^{\alpha }V^{\beta }$
  
  * Avec le lagrangien relativiste : $~L=-mc~{\sqrt {V_{\alpha }V^{\alpha }}}-e~A_{\alpha }V^{\alpha }$, les équations d'Euler-Lagrange : ${\frac {d~}{d\tau }}({\frac {\partial L}{\partial V^{\alpha }}})\ -\ {\frac {\partial L}{\partial x^{\alpha }}}=0$ donnent :
  
         ${\frac {d~}{d\tau }}(-mV_{\alpha }-e~A_{\alpha })+e~\partial _{\alpha }A_{\beta }~V^{\beta }=0$
  * Sachant que $A_{\alpha }$ ne dépend pas explicitement du paramètre $\tau$ (le temps fait partie des $x^{\alpha }$ et n'est pas le paramètre), on peut écrire :
  
         ${\frac {dA_{\alpha }}{d\tau }}=\partial _{\beta }A_{\alpha }~{\frac {dx^{\beta }}{d\tau }}=\partial _{\beta }A_{\alpha }~V^{\beta }$.
  * En posant $F_{\alpha \beta }=\partial _{\alpha }A_{\beta }-\partial _{\beta }A_{\alpha }$ tenseur champ électromagnétique, on obtient ainsi l'équation du mouvement :
  
         $m~{\frac {dV_{\alpha }}{d\tau }}=e~F_{\alpha \beta }~V^{\beta }$.
  * On retrouve la quadri-écriture plus usuelle en élevant l'indice $\alpha$ avec $\eta ^{\alpha \beta }$.
  * Le tenseur champ électromagnétique $F_{\alpha \beta }$ comporte 16 composantes ; puisqu'il est antisymétrique par construction, les composantes diagonales sont nulles et seules six des autres sont indépendantes.
  
      On les note généralement par deux « pseudo-notations vectorielles tridimensionnelles » :
         • Champ électrique : ${\vec {E}}=c~.\left(F^{10},F^{20},F^{30}\right)$ ; par exemple : $E^{1}=-\partial _{1}\phi -{\frac {\partial A^{1}}{\partial t}}=c~\partial ^{1}A^{0}-c~\partial ^{0}A^{1}$.
         • Champ magnétique : ${\vec {B}}=\left(F^{32},F^{13},F^{21}\right)$ ; par exemple : $B^{1}=\partial _{2}A^{3}-\partial _{3}A^{2}=-\partial ^{2}A^{3}+\partial ^{3}A^{2}$.
      En décomposant la relation précédente pour l'exprimer avec ces notations, on retrouve la force de Lorentz sous son écriture habituelle : ${\frac {d{\vec {p}}}{dt}}=e\ {\vec {E}}+e\ {\vec {v}}\times {\vec {B}}$.
      ♦ remarque : vis-à-vis des transformations spatiales dans $\mathbb {R} ^{3}$, ${\vec {E}}$ se comporte comme un vecteur et ${\vec {B}}$ se comporte comme un pseudovecteur ; vis-à-vis des transformations spatio-temporelles galiléennes dans $\mathbb {R} ^{4}$ il faut appliquer la transformation de Lorentz sur l'objet global $F_{\alpha \beta }$.
      ♦ remarque : attention à ne pas confondre l'énergie $E$ et le champ électrique ${\vec {E}}$.
  * On obtient donc les équations du mouvement correctes, mais avec la quadri-impulsion « contre intuitive » : $P_{\alpha }={\frac {\partial L}{\partial V^{\alpha }}}=-mV_{\alpha }-e~A_{\alpha }$, c'est-à-dire : $P^{\alpha }=-p^{\alpha }-e~A^{\alpha }$.
  * Le carré de la « norme » de la quadri-impulsion permet d'écrire : $\ p^{\alpha }p_{\alpha }=(-P^{\alpha }-eA^{\alpha }).(-P_{\alpha }-eA_{\alpha })=m^{2}c^{2}$, donc « vaguement » conforme à la relation usuelle, mais avec la convention de signe inversée : $(-E_{tot}-e\phi )^{2}-c^{2}.(-{\vec {P}}-e{\vec {A}})^{2}=m^{2}c^{4}$.
  * Le hamiltonien correspondant peut être défini par : $H=P_{\alpha }V^{\alpha }-L=-mc^{2}+mc^{2}=0$. Écrit ainsi il ne décrit pas une énergie, mais il est vrai que ce rôle est ici dédié à $P^{0}$.
  
  **Formulation en changeant la convention de signe de l'impulsion**
  
  * Pour contourner ce problème de signe, certains physiciens (par exemple J.M. Raimond ; Électromagnétisme et Relativité ; cours à l'E.N.S.) suggèrent d'adapter la méthode de Hamilton en changeant la convention de signe de l'impulsion :
  
  ${\mathcal {P}}_{\alpha }=-\partial _{\alpha }S=-{\frac {\partial L}{\partial V^{\alpha }}}=mV_{\alpha }+e~A_{\alpha }$
  
      Bien sûr cela nécessite d'adapter les signes dans d'autres équations, entre autres celles d'Euler-Lagrange : $~{\frac {d{\mathcal {P}}_{\alpha }}{d\tau }}=-{\frac {\partial L}{\partial x_{\alpha }}}~$.
  * Par la même méthode que celle utilisée dans la partie précédente, on obtient dans ce cas : $\ p^{\alpha }p_{\alpha }=({\mathcal {P}}^{\alpha }-eA^{\alpha }).({\mathcal {P}}_{\alpha }-eA_{\alpha })=m^{2}c^{2}$, donc conforme à la relation usuelle : $(E_{tot}-e\phi )^{2}-c^{2}.({\vec {\mathcal {P}}}-e{\vec {A}})^{2}=m^{2}c^{4}$.
  * En modifiant de façon concordante l'expression du hamiltonien, on obtient de même : $H=-{\mathcal {P}}_{\alpha }V^{\alpha }-L=-mc^{2}+mc^{2}=0$.
  
  **Formulation en changeant la convention de signe de l'action**
  
  * Pour contourner le problème de signe, d'autres physiciens (par exemple B. Houchmandzadeh ; Théorie de la Relativité ; cours à l'université de Grenoble) suggèrent de changer le signe de l'action, dont en fait seule la condition d'extremum est utilisée (contrairement à ce que suggère le « principe de moindre action », qu'il faudrait peut être plutôt nommer « principe d'action stationnaire »).
  
  $L=mc~{\sqrt {V_{\alpha }V^{\alpha }}}+e~A_{\alpha }V^{\alpha }$
  
      ♦ remarque : cet effet ne peut pas être obtenu en utilisant la métrique en notation anglo-saxonne ($ds^{2}=dx^{2}+dy^{2}+dz^{2}-c^{2}\ dt^{2}$) car elle est associée à plusieurs modifications de signes qui se compensent (ce qui semble assez logique pour des notations « covariantes »).
      ♦ remarque : il faut dans ce cas changer le signe de tous les termes de l'action (y compris le terme du champ électromagnétique étudié dans la suite pour obtenir les équations de Maxwell).
  * On obtient dans ce cas naturellement :
  
  $P_{\alpha }={\frac {\partial L}{\partial V^{\alpha }}}=mV_{\alpha }+e~A_{\alpha }$
  
  * Il peut sembler « incorrect » d'utiliser pour l'action relativiste une expression dont la limite aux faibles vitesses ne redonne pas (ici uniquement pour le signe) l'expression non relativiste. Toutefois, l'expression du lagrangien n'est a priori pas unique : l'important est que ce dernier prédise des lois physiques dont la limite redonne les lois non relativistes, ce qui est le cas ici.
  
  **Formulation quadratique (lagrangien géodésique)**
  
  * Si on paramètre par $\tau$, on peut multiplier le lagrangien par ${\frac {ds}{c~d\tau }}=1$. On obtient ainsi pour un point matériel isolé le lagrangien quadratique $L=m~{\frac {ds^{2}}{d\tau ^{2}}}=m~\eta _{\alpha \beta }V^{\alpha }V^{\beta }$.
  
      Par construction, pour les faibles vitesses, l'action ainsi obtenue tend (au signe près) vers l'expression non relativiste. Par contre cela correspondrait à une impulsion inadaptée : $P_{\alpha }={\frac {\partial L}{\partial V^{\alpha }}}=2mV_{\alpha }$.
  * On convient alors généralement (par exemple : R. K. Nesbet ; Variational Principles and Methods ; P. Tourrenc ; Relativité et Gravitation) de multiplier le lagrangien par ${\frac {1}{2}}$, ce qui donne une expression de forme analogue à celle du cas non relativiste (généralement nommée lagrangien géodésique) :
  
  $L={\frac {1}{2}}m~\eta _{\alpha \beta }V^{\alpha }V^{\beta }$
  
      ♦ remarque : un avantage important de cette expression, quadratique mais sans radical, est de simplifier un certain nombre de calculs de mouvements.
      ♦ remarque : On peut reprocher à l'action ainsi définie de ne pas tendre vers l'expression non relativiste dans la limite des faibles vitesses (non seulement pour le signe, mais aussi pour le facteur ${\frac {1}{2}}$). Toutefois, le résultat tend vers la limite non relativiste d'une autre façon : sa forme est analogue, vis-à-vis du paramètre $\tau$, à celle du cas non relativiste vis-à-vis du paramètre $t$. Or $\tau$ tend vers $t$ pour les faibles vitesses. De ce fait, les équations du mouvement qui s'en déduisent tendent alors vers les lois non relativistes.
  * Pour contourner l'ambiguïté de l'écriture formelle avec $\eta _{\alpha \beta }V^{\alpha }V^{\beta }=c^{2}$ (qui n'est qu'apparente), on peut choisir un paramètre $\sigma$ « arbitraire », en notant ${\underline {L}}(x^{\alpha }(\sigma );{\underline {V}}^{\alpha }(\sigma );\sigma )$ avec ${\underline {V}}^{\alpha }={\frac {dx^{\alpha }}{d\sigma }}$, quitte à considérer ensuite la limite $\sigma \rightarrow \tau$ par continuité.
  * Il faut toutefois prendre garde que dans ce cas le paramètre $\sigma$ n'est pas « muet ». Dès lors qu'on choisit un paramètre « quelconque » (non affine en fonction de $\tau$), il faut respecter la condition : $S=\int L~d\tau =\int {\underline {L}}~d\sigma$. Ceci impose : ${\underline {L}}={\frac {1}{2}}m~\eta _{\alpha \beta }V^{\alpha }{\underline {V}}^{\beta }={\frac {1}{2}}m~\eta _{\alpha \beta }{\underline {V}}^{\alpha }{\underline {V}}^{\beta }{\frac {d\sigma }{d\tau }}$.
  
      On peut alors utiliser ${\frac {c^{2}~d\tau ^{2}}{d\sigma ^{2}}}=\eta _{\alpha \beta }{\underline {V}}^{\alpha }{\underline {V}}^{\beta }$ pour obtenir : ${\underline {L}}={\frac {mc}{2}}{\sqrt {\eta _{\alpha \beta }{\underline {V}}^{\alpha }{\underline {V}}^{\beta }}}$.
      Ceci correspond à une quadri-impulsion : $P_{\alpha }={\frac {m}{2}}{\frac {dx_{\alpha }}{d\sigma }}{\frac {d\sigma }{d\tau }}={\frac {m}{2}}V_{\alpha }$. Le résultat pour $\sigma =\tau$ ne correspond pas à la limite pour $\sigma \rightarrow \tau$.
      ♦ remarque : cela était prévisible puisque le passage au lagrangien quadratique utilise une multiplication opportuniste par ${\frac {ds}{c~d\tau }}=1$, inadéquate pour $\sigma \neq \tau$.
      ♦ remarque : la paramétrisation directement par $\sigma \neq \tau$ est possible sans le coefficient ${\frac {1}{2}}$ mais elle ne donne pas un lagrangien quadratique.
  * On peut traiter ainsi le cas plus général d'une particule en interaction avec un champ électromagnétique (sans oublier qu'il faut changer le signe de tous les termes) :
  
  $L={\frac {1}{2}}m~V_{\alpha }V^{\alpha }+e~A_{\alpha }V^{\alpha }$
  
  * Ceci correspond à l'impulsion (donnant correctement les équations du mouvement vues précédemment, ce qui en suivant sur ce point le raisonnement de S. Weinberg « qualifie l'expression utilisée pour l'action comme appropriée pour ce système ») :
  
  $P_{\alpha }={\frac {\partial L}{\partial V^{\alpha }}}=mV_{\alpha }+e~A_{\alpha }$
  
  * On peut en outre dans ce cas considérer le hamiltonien (exprimé en fonction le l'impulsion généralisée) :
  
  $H=\eta _{\alpha \beta }P^{\alpha }V^{\beta }-L={\frac {1}{2m}}\eta _{\alpha \beta }(P^{\alpha }-e~A^{\alpha })(P^{\beta }-e~A^{\beta })$
  
  * Conformément à la méthode de Hamilton, ceci redonne : ${\frac {\partial H}{\partial P_{\alpha }}}=V^{\alpha }$. En outre, dans la mesure où le hamiltonien ne dépend pas explicitement de $\tau$, il correspond à une constante du mouvement : $H={\frac {mc^{2}}{2}}$.
  * Les équations du mouvement peuvent s'écrire : ${\frac {dP^{\alpha }}{d\tau }}=-{\frac {\partial H}{\partial x^{\alpha }}}={\frac {e}{m}}\eta _{\mu \nu }(P^{\mu }-e~A^{\mu })(\partial _{\alpha }A^{\nu })=e~V^{\beta }(\partial _{\alpha }A_{\beta })$.
  
      Par ailleurs : ${\frac {dP_{\alpha }}{d\tau }}={\frac {dp_{\alpha }}{d\tau }}+e~{\frac {dA_{\alpha }}{d\tau }}$ ; or $A^{\alpha }(x^{\beta })$ ne dépend pas explicitement $\tau$ donc : ${\frac {dA_{\alpha }}{d\tau }}=\partial _{\beta }A_{\alpha }{\frac {dx_{\beta }}{d\tau }}$.
      Ainsi : ${\frac {dp_{\alpha }}{d\tau }}=e~F_{\alpha \beta }V^{\beta }$.
  
  **Description d'un photon**
  
  * La description corpusculaire d'un photon peut être étudiée à partir de celle d'une particule massive, en passant à la limite pour $m\rightarrow 0$. Ceci n'est possible qu'avec un lagrangien quadratique, sinon il apparait forcément $d\tau =0$ au dénominateur.
  
      Avec une notation quadratique, l'action d'une particule massive peut s'écrire : $S={\frac {m}{2}}\int {\frac {ds^{2}}{d\tau ^{2}}}d\tau ={\frac {1}{2}}\int \eta _{\mu \nu }p^{\mu }dx^{\nu }$ avec l'impulsion $p^{\mu }=mV^{\mu }={\frac {m}{\sqrt {1-\beta ^{2}}}}{\frac {dx^{\mu }}{dt}}$.
      Le point délicat pour un photon est l'indétermination de $V^{\mu }$ due à la limite $\beta \rightarrow 1$ quand $~m\rightarrow 0$. Puisque la relation $~E^{2}=m^{2}c^{4}+c^{2}p^{2}~$ donne la limite $~E=h\nu =cp~$ pour l'énergie-impulsion d'un photon, la limite cherchée pour $~v\rightarrow c~$ peut être obtenue en substituant : ${\frac {m}{\sqrt {1-\beta ^{2}}}}\rightarrow {\frac {h\nu }{c^{2}}}$.
      Ceci donne : $~S={\frac {h\nu }{2c^{2}}}\int \eta _{\mu \nu }v^{\mu }v^{\nu }dt~~$ avec $~v^{\mu }={\frac {dx^{\mu }}{dt}}~$ donc un lagrangien quadratique paramétré par $t$ :
  
  $L={\frac {h\nu }{2c^{2}}}\eta _{\alpha \beta }v^{\alpha }v^{\beta }$
  
      ♦ remarque : la méthode peut s'adapter à la relativité générale, mais avec un paramètre (là encore imposé) différent de $t$.
  * On peut se demander si la démarche utilisée ici peut être applicable au lagrangien avec radical, quitte à utiliser un autre paramètre $\sigma$.
  * On considère alors l'action : $~S=mc\int {\frac {ds}{d\sigma }}d\sigma =mc\int {\sqrt {{\underline {V}}_{\mu }{\underline {V}}^{\mu }}}d\sigma ~$ avec $~{\underline {V}}^{\mu }={\frac {dx^{\mu }}{d\sigma }}$.
  
      La limite peut être obtenue en substituant : $mc\rightarrow {\frac {h\nu }{c^{2}}}{\sqrt {\frac {ds^{2}}{dt^{2}}}}={\frac {h\nu }{c^{2}}}{\sqrt {{\underline {V}}_{\mu }{\underline {V}}^{\mu }}}~{\frac {d\sigma }{dt}}$.
      Finalement, ceci correspond forcément à un lagrangien quadratique : $~S={\frac {h\nu }{c^{2}}}\int {\underline {V}}_{\mu }{\underline {V}}^{\mu }{\frac {d\sigma }{dt}}d\sigma ={\frac {h\nu }{c^{2}}}\int \eta _{\alpha \beta }v^{\alpha }v^{\beta }dt$.
      Dans la mesure ou ce lagrangien quadratique donne les bonnes équations du mouvement, mais des impulsions incorrectes d'un facteur 2, on le multiplie par un coefficient ${\frac {1}{2}}$.
  * Puisque $~{\frac {\partial L}{\partial x^{\alpha }}}=0~$ ceci donne pour un photon l'impulsion constante : $p_{\alpha }={\frac {\partial L}{\partial v^{\alpha }}}={\frac {h\nu }{c^{2}}}v_{\alpha }$.
  * On obtient par ailleurs le hamiltonien (exprimé en fonction de l'impulsion) : $H=\eta _{\alpha \beta }p^{\alpha }v^{\beta }-L={\frac {c^{2}}{2h\nu }}\eta _{\alpha \beta }p^{\alpha }p^{\beta }=L$.
  
      Ceci donne, conformément aux relations de Hamilton : $~{\frac {\partial H}{\partial p_{\alpha }}}=v^{\alpha }$.
  
  ## Étude 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 influencé, il se propage quelque chose dans l'espace, en général à la vitesse de la lumière, dont l'effet est un changement de trajectoire du corps influencé.
      Suivant quelles propriétés ce champ (appelé ainsi car il décrit une grandeur physique définie en chaque point de l'espace) est-il créé, se propage-t-il, est-il influencé par son environnement, etc ?
      On peut répondre à ces questions à l'aide du principe de moindre action.
  * À 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.
  
      ♦ remarque : dans cette partie, puisqu'on ne se préoccupe pas de l'énergie-impulsion des particules chargées, on utilise la convention de signe négatif du lagrangien (il se trouve que tous les termes sont négatifs).
  
  ### Invariance de jauge du potentiel et du tenseur électromagnétiques
  
  * Si on remplace le quadri-potentiel électromagnétique $~A_{\alpha }~$ par le quadri-potentiel $~A'_{\alpha }=A_{\alpha }+\partial _{\alpha }\varphi ~$ où $~\varphi ~$ est une fonction quelconque des coordonnées, alors le terme d'interaction du lagrangien devient $~L'=L-e~\partial _{\alpha }\varphi ~.V^{\alpha }~$ et l'action $~S'=S-e~\int _{\tau _{i}}^{\tau _{f}}\partial _{\alpha }\varphi ~.dx^{\alpha }=S-e~.(\varphi (x^{\alpha }(\tau _{f}))-\varphi (x^{\alpha }(\tau _{i})))$.
  
      ♦ remarque : ne pas confondre $~\varphi ~$ avec le potentiel « scalaire » $~\phi =c~A^{0}$.
  * En appliquant la méthode variationnelle qui compare les chemins en gardant les extrémités fixes, le terme $~\varphi (x^{\alpha }(\tau _{f}))-\varphi (x^{\alpha }(\tau _{i}))~$ est éliminé.
  
      Ainsi les deux potentiels $~A_{\alpha }~$ et $~A'_{\alpha }=A_{\alpha }+\partial _{\alpha }\varphi ~$ donnent les mêmes équations du mouvement : cette propriété caractérise ce qu'on nomme « invariance de jauge ».
      On constate d'ailleurs que dans les équations du mouvement, le tenseur champ électromagnétique $~F_{\alpha \beta }~$ est invariant de jauge : $~F'_{\alpha \beta }=\partial _{\alpha }A'_{\beta }-\partial _{\beta }A'_{\alpha }=\partial _{\alpha }A_{\beta }-\partial _{\beta }A_{\alpha }+\partial _{\alpha }\partial _{\beta }\varphi -\partial _{\beta }\partial _{\alpha }\varphi =\partial _{\alpha }A_{\beta }-\partial _{\beta }A_{\alpha }=F_{\alpha \beta }~$ puisque selon le théorème de Schwarz : $~\partial _{\alpha }\partial _{\beta }=\partial _{\beta }\partial _{\alpha }$.
      ♦ remarque : la grandeur généralement nommée « champ électromagnétique » est le tenseur champ électromagnétique ; il peut se déduire du potentiel et la description des phénomènes électromagnétiques peut se faire à l'aide de l'un ou l'autre (le principe de moindre action utilise le potentiel, mais il n'est pas indispensable à la description des lois physiques) ; au sens général, rien n'interdit de considérer le potentiel comme un « champ vectoriel » mais ce n'est pas l'usage.
      ♦ remarque : cette équivalence des descriptions par le champ tensoriel ou le potentiel vectoriel a été discutée à propos de l'effet Aharonov-Bohm.
  
  ### Densité lagrangienne et équations d'Euler-Lagrange associées
  
  * Un champ correspond à une grandeur physique ayant une valeur en chaque point de l'espace, on ne peut donc pas le repérer comme un point matériel et le décrire par des coordonnées $~x^{\alpha }(\tau )~$ et leurs dérivées $~V^{\alpha }$.
  
      Le potentiel quadri-vectoriel peut être décrit par ses coordonnées $~A_{\alpha }~$ pour chaque point de l'espace : $~A_{\alpha }(x^{\beta })~$.
      Les quantités jouant un rôle analogue à celui de la vitesse sont les dérivées des coordonnées : $~{\frac {\partial A_{\alpha }}{\partial x^{\beta }}}=\partial _{\beta }A_{\alpha }$.
      Ainsi les coordonnées $~x^{\alpha }~$ du point où on décrit le potentiel jouent le rôle de paramètres (analogues du paramètre $~\tau ~$ de la position du point matériel sur sa trajectoire).
      ♦ remarque : on peut en principe aussi utiliser les coordonnées contra-variantes $~A^{\alpha }~$ mais cela complique inutilement certains calculs.
  * L'action d'un champ peut alors s'écrire sous la forme : $~S={\frac {1}{c}}\int _{\mathcal {V}}\Lambda (A_{\alpha };\partial _{\beta }A_{\alpha })~d\Omega ~$, où $~{\mathcal {V}}~$ est le quadri-volume dans lequel on applique la méthode variationnelle et où $~d\Omega =dx^{0}.dx^{1}.dx^{2}.dx^{3}=c.dt.dx^{1}.dx^{2}.dx^{3}~$ est l'élément infinitésimal de quadri-volume.
  
      La fonction $~\Lambda ~$ est appelée « densité lagrangienne » (si on applique totalement l'analogie formelle, elle pourrait en principe aussi dépendre explicitement des paramètres $~x^{\alpha }~$).
      ♦ remarque : le coefficient $~{\frac {1}{c}}~$ dans l'expression de l'action sert à définir $~\Lambda ~$ comme une densité volumique d'énergie (on intègre sur $~ct$).
  * Par une démonstration semblable à celle déjà 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 :
  
  $\partial _{\beta }{\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}-{\frac {\partial \Lambda }{\partial A_{\alpha }}}=0$
  
  * On considère ici la variation de l'action correspondant à des variations du potentiel $~A_{\alpha }$, supposées nulles aux limites du quadri-volume d'intégration $~{\mathcal {V}}~$.
  
      Cette variation de l'action peut s'écrire : $~\delta S={\frac {1}{c}}\int _{\mathcal {V}}\left({\frac {\partial \Lambda }{\partial A_{\alpha }}}~\delta A_{\alpha }+{\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}~\delta (\partial _{\beta }A_{\alpha })\right)~d\Omega ~$.
  * Le second terme peut s'exprimer : $~\int _{\mathcal {V}}\left({\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}~\delta (\partial _{\beta }A_{\alpha })\right)~d\Omega =\int _{\mathcal {V}}\left({\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}~\partial _{\beta }(\delta A_{\alpha })\right)~d\Omega =\int _{\mathcal {V}}\partial _{\beta }\left({\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}~\delta A_{\alpha }\right)~d\Omega -\int _{\mathcal {V}}\left(\partial _{\beta }({\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}})~\delta A_{\alpha }\right)~d\Omega ~$.
  
      Or le quadrivecteur : $~{\mathfrak {V}}^{\beta }={\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}~\delta A_{\alpha }~$ a un flux nul à travers l'hypersurface « bordant » le quadri-volume d'intégration $~{\mathcal {V}}~$, puisque $~\delta A_{\alpha }~$ s'y annule.
      D'après le théorème d'Ostrogradski, l'intégrale de sa divergence est nulle : $~\int _{\mathcal {V}}\partial _{\beta }{\mathfrak {V}}^{\beta }~d\Omega =0~$.
      Ainsi en simplifiant : $~\delta S={\frac {1}{c}}\int _{\mathcal {V}}\left({\frac {\partial \Lambda }{\partial A_{\alpha }}}-\partial _{\beta }({\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}})\right)~\delta A_{\alpha }~d\Omega ~$.
  * La variation de l'action est donc nulle, quelles que soient les variations $~\delta A_{\alpha }~$ au voisinage de l'extremum, si et seulement si est vérifiée la condition : $~{\frac {\partial \Lambda }{\partial A_{\alpha }}}-\partial _{\beta }({\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}})=0~$.
  
  ### Tenseur impulsion-énergie d'un champ
  
  * La densité lagrangienne $~\Lambda =\Lambda (A_{\alpha };\partial _{\beta }A_{\alpha })~$ associée au champ électromagnétique étant donnée, on peut (de manière analogue à la construction du hamiltonien associé à l'énergie d'une particule) définir un « tenseur énergie-impulsion » associé au champ :
  
  $T_{~\alpha }^{\gamma }={\frac {\partial \Lambda }{\partial (\partial _{\gamma }A_{\beta })}}.\partial _{\alpha }A_{\beta }-\delta _{~\alpha }^{\gamma }\Lambda$
  
  * De même que le hamiltonien décrit une énergie constante s'il ne dépend pas explicitement du paramètre $\tau$, le tenseur énergie-impulsion décrit une grandeur « conservée » s'il ne dépend pas explicitement des paramètres $x^{\alpha }$ :
  
  $\partial _{\gamma }T_{~\alpha }^{\gamma }=0$
  
  * On peut considérer : $~\partial _{\alpha }\Lambda ={\frac {\partial \Lambda }{\partial A_{\beta }}}.\partial _{\alpha }A_{\beta }+{\frac {\partial \Lambda }{\partial (\partial _{\gamma }A_{\beta })}}.\partial _{\alpha }(\partial _{\gamma }A_{\beta })$.
  
      Puis en utilisant les équations d'Euler-Lagrange : $~{\frac {\partial \Lambda }{\partial A_{\beta }}}=\partial _{\gamma }{\frac {\partial \Lambda }{\partial (\partial _{\gamma }A_{\beta })}}~$ et le lemme de Schwarz : $~\partial _{\alpha }(\partial _{\gamma }A_{\beta })=\partial _{\gamma }(\partial _{\alpha }A_{\beta })~$ on obtient :
         $~\partial _{\alpha }\Lambda =\partial _{\gamma }{\frac {\partial \Lambda }{\partial (\partial _{\gamma }A_{\beta })}}.\partial _{\alpha }A_{\beta }+{\frac {\partial \Lambda }{\partial (\partial _{\gamma }A_{\beta })}}.\partial _{\gamma }(\partial _{\alpha }A_{\beta })=\partial _{\gamma }\left({\frac {\partial \Lambda }{\partial (\partial _{\gamma }A_{\beta })}}.\partial _{\alpha }A_{\beta }\right)$.
      L'égalité $~\partial _{\alpha }\Lambda =\delta _{~\alpha }^{\gamma }\partial _{\gamma }\Lambda ~$ permet finalement d'écrire :
         $\partial _{\gamma }\left({\frac {\partial \Lambda }{\partial (\partial _{\gamma }A_{\beta })}}.\partial _{\alpha }A_{\beta }-\delta _{~\alpha }^{\gamma }\Lambda \right)=0$.
  * La quantité : $~w=T_{~0}^{0}={\frac {\partial \Lambda }{\partial (\partial _{0}A_{\beta })}}.\partial _{0}A_{\beta }-\Lambda ~$ décrit une « densité d'énergie » et les quantités $~\Pi ^{k}=c~.T_{~0}^{k}=c~{\frac {\partial \Lambda }{\partial (\partial _{k}A_{\beta })}}.\partial _{0}A_{\beta }~(k=1;2;3)~$ décrivent les composantes d'un vecteur « densité de courant d'énergie » $~{\vec {\Pi }}$.
  
      L'équation $~\partial _{\gamma }T_{~0}^{\gamma }=0~$ est une « équation de conservation de l'énergie » : localement, la variation dans le temps de la densité d'énergie $~w~$ dans un volume infinitésimal est égale au flux entrant de la densité de courant d'énergie $~\Pi ^{k}~$ à travers la surface entourant ce volume (d'après le théorème d'Ostrogradski).
  
  ${\frac {\partial w}{\partial t}}+div(~{\vec {\Pi }}~)=0$
  
      ♦ remarque : les composantes $T_{~0}^{\gamma }$ ne définissent par contre pas un quadri-vecteur (c'est pour cela qu'on utilise la grandeur tensorielle $T_{~\alpha }^{\gamma }$).
  
  ### Densité lagrangienne d'un champ électromagnétique libre
  
  * La densité lagrangienne décrivant le champ électromagnétique est :
  
  $\Lambda _{em}=-{\frac {1}{4\mu _{0}}}F^{\alpha \beta }F_{\alpha \beta }$
  
  * On cherche la densité lagrangienne du champ électromagnétique qui est composé du ou des nombres construits à partir du potentiel électromagnétique $~A^{\alpha }~$ qui sont invariants par changement de référentiel dans l'espace de Minkowski.
  
      La manifestation de ce potentiel correspond au champ électrique $~{\vec {E}}=c~.(F^{10},F^{20},F^{30})~$ et le champ magnétique $~{\vec {B}}=(F^{32},F^{13},F^{21})$.
  * On cherche donc les nombres invariants construits à partir de la matrice 4×4 associée : $~{\boldsymbol {F}}=(F_{~\beta }^{\alpha })$.
  
      Les coefficients du polynôme caractéristique sont invariants par la transformation d'endomorphismes $~{\boldsymbol {F}}\to {\boldsymbol {L}}^{-1}.{\boldsymbol {F}}.{\boldsymbol {L}}~$, où $~{\boldsymbol {L}}~$ 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 $~P_{F}(x)=det({\boldsymbol {F}}-x.{\boldsymbol {Id}})$.
      On sait, par la définition de ses coefficients, que cette matrice est anti-symétrique : $~^{t}{\boldsymbol {F}}=-{\boldsymbol {F}}~$ (en notant $~^{t}{\boldsymbol {F}}~$ la matrice transposée).
      On en déduit que $\ P_{F}(x)$ est un polynôme pair. Soit : $~P_{F}(x)=x^{4}+a.x^{2}+det({\boldsymbol {F}})$.
      Par quelques calculs, on montre que $~a=-F_{~\beta }^{\alpha }F_{~\alpha }^{\beta }=F^{\alpha \beta }F_{\alpha \beta }=-2.({\frac {{\vec {E}}^{2}}{c^{2}}}-{\vec {B}}^{2})~$ et que $~det({\boldsymbol {F}})=-({\frac {\vec {E}}{c}}.{\vec {B}})^{2}$.
  * Pour une raison de dimension , le nombre $~({\frac {\vec {E}}{c}}.{\vec {B}})^{2}~$ ne convient pas (la densité lagrangienne doit avoir la dimension d'une densité d'énergie), il faudrait essayer avec $~{\sqrt {({\frac {\vec {E}}{c}}.{\vec {B}})^{2}}}=\left|{\frac {\vec {E}}{c}}.{\vec {B}}\right|=\pm {\frac {\vec {E}}{c}}.{\vec {B}}$.
  
      Mais par ailleurs on peut montrer que $~{\vec {E}}.{\vec {B}}~$ est la 4-divergence d'une 4-fonction, donc son ajout à la densité lagrangienne ne changerait en rien les équations d'Euler-Lagrange. Le nombre $~\pm {\frac {\vec {E}}{c}}.{\vec {B}}~$ peut donc aussi être écarté.
      ♦ remarque : pour une particule ponctuelle, l'ajout au lagrangien de la dérivée $~{\frac {df}{d\tau }}~$ d'une fonction ne change pas les équations car cela ajoute à l'action la variation $~f(\tau _{fin})-f(\tau _{ini})~$ qui est constante lorsqu'on fait varier l'action ; de même pour une densité lagrangienne, l'ajout de la 4-divergence $~\partial _{\gamma }f^{\gamma }~$ d'un 4-vecteur ne fait qu'ajouter à l'action le flux du 4-vecteur à travers l'hypersurface limitant le domaine d'intégration, or cette quantité reste constante lorsqu'on fait varier l'action (les variations aux limites sont nulles).
  * De ce fait, $~F^{\alpha \beta }F_{\alpha \beta }~$ 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 : $~\Lambda _{em}=-{\frac {1}{4\mu _{0}}}F^{\alpha \beta }F_{\alpha \beta }$
      ♦ remarque : le signe $-$ correspond au choix lié à la minimisation de l'action, ainsi qu'aux choix des signes des autres termes de l'action (commentés dans les parties précédentes).
  
  ### Les équations du champ électromagnétique
  
  * L'hypothèse de ce paragraphe est qu'il y a un courant $~J^{\alpha }$ de particules chargées (même s'il n'y en a éventuellement qu'une seule) ; on n'étudie donc pas le champ « libre » mais directement le champ en interaction avec des charges.
  
      On suppose le courant « non influencé » par le champ électromagnétique. Avec cette condition, on étudie les variations de l'action associées aux variations du champ (décrit par le potentiel).
      ♦ remarque : ceci signifie que, pour nettement simplifier les calculs, on étudie les variations de l'action associées aux variations du 4-potentiel (sans faire varier les mouvements des particules chargées) ; de façon analogue, l'étude de particules chargées subissant un champ électromagnétique se fait en supposant le champ fixé (sans faire varier le 4-potentiel, qui pourtant peut être modifié par les courants associés aux mouvements des charges) ; cette séparation est possible car les termes associés aux deux types de variations se factorisent séparément.
  * Dans la densité lagrangienne, le terme d'interaction peut ici s'écrire : $-A_{\alpha }J^{\alpha }$. La densité lagrangienne à utiliser est ainsi :
  
  $\Lambda =-A_{\alpha }J^{\alpha }-{\frac {1}{4\mu _{0}}}F^{\alpha \beta }F_{\alpha \beta }$
  
  * Pour décrire les mouvements des particules chargées, l'interaction entre une particule et le champ se modélise dans l'action par : $-e~A_{\alpha }V^{\alpha }~d\tau =-e~A_{\alpha }dx^{\alpha }$.
  
      L'utilisation d'une densité lagrangienne nécessite par contre ici de raisonner avec une densité de charge (quitte à utiliser pour chaque particule ponctuelle une distribution de Dirac) : $~\rho ={\frac {d^{3}e}{d^{3}x}}~$ avec $~d^{3}x=dx~dy~dz=dx^{1}dx^{2}dx^{3}$.
      Le terme d'interaction de l'action peut donc s'écrire :
         $-\sum _{part}(~e\int A_{\alpha }dx^{\alpha })=-\sum _{part}(~e\int A_{\alpha }{\frac {dx^{\alpha }}{dt}}~dt)=-\int \int \rho ~A_{\alpha }{\frac {dx^{\alpha }}{dt}}~d^{3}x~dt=-{\frac {1}{c}}\int _{\mathcal {V}}A_{\alpha }~\rho ~{\frac {dx^{\alpha }}{dt}}~d\Omega$.
      ♦ remarque : les intégrations sur $~d^{3}x~$ et $~dt~$ ne peuvent pas être simplement séparées car $~A_{\alpha }~$ et $~{\frac {dx^{\alpha }}{dt}}~$ doivent être calculés au même point que $~\rho ~$ ; si on intègre d'abord sur $~d^{3}x~$ on peut utiliser (pour chaque particule) $~\rho =e~\delta ^{3}(x^{k}-X^{k}(t))~$ où $~X^{k}(t)~$ correspond à la position à l'instant considéré ; on peut obtenir une écriture plus symétrique sur $~d\Omega ~$ avec $~\rho =e~\delta ^{4}(x^{\alpha }-X^{\alpha }(t))~$ où $~X^{0}(t)=ct~$ (puisqu'ici on paramètre par $~t$).
  * On peut alors définir un quadrivecteur densité de courant : $J^{\alpha }=\rho ~{\frac {dx^{\alpha }}{dt}}=(\rho ~c~;\rho ~{\vec {v}})=(\rho ~c~;{\vec {j}})$.
  
      Ainsi le terme d'interaction peut s'écrire : $~-\sum _{part}(~e.\int A_{\alpha }dx^{\alpha })=-{\frac {1}{c}}\int _{\mathcal {V}}A_{\alpha }J^{\alpha }~d\Omega$.
      ♦ remarque : il est important de noter que, par rapport à la transformation de Lorentz, la densité de charge $~\rho ~$ n'est pas un scalaire (la charge est un scalaire, mais la contraction des volumes modifie $~\rho ~$) ; par ailleurs, le 4-objet $~{\frac {dx^{\alpha }}{dt}}=(c~;{\vec {v}})~$ n'est pas un 4-vecteur à cause de la dilatation des durées ; il se trouve que les deux effets se compensent et que le produit des deux est un quadrivecteur.
  * On obtient ainsi au total : $~S=-{\frac {1}{c}}\int _{\mathcal {V}}(A_{\alpha }J^{\alpha }+{\frac {1}{4\mu _{0}}}F^{\alpha \beta }F_{\alpha \beta })~d\Omega ={\frac {1}{c}}\int _{\mathcal {V}}\Lambda ~d\Omega$.
  * Les équations d'Euler-Lagrange donnent :
  
  $\partial _{\alpha }F^{\alpha \beta }=\mu _{0}~J^{\beta }$
  
  * On obtient :
  
         $\Lambda =-A_{\alpha }J^{\alpha }-{\frac {1}{4\mu _{0}}}\eta ^{\mu \rho }\eta ^{\nu \sigma }F_{\mu \nu }F_{\rho \sigma }$
         ${\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}=-{\frac {1}{4\mu _{0}}}\left(F^{\mu \nu }{\frac {\partial }{\partial (\partial _{\beta }A_{\alpha })}}(\partial _{\mu }A_{\nu }-\partial _{\nu }A_{\mu })~+~F^{\rho \sigma }{\frac {\partial }{\partial (\partial _{\beta }A_{\alpha })}}(\partial _{\rho }A_{\sigma }-\partial _{\sigma }A_{\rho })\right)~$ ;
         ${\frac {\partial \Lambda }{\partial (\partial _{\beta }A_{\alpha })}}=-{\frac {1}{4\mu _{0}}}\left(F^{\mu \nu }(\delta _{\mu }^{~\beta }~\delta _{\nu }^{~\alpha }-\delta _{\nu }^{~\beta }~\delta _{\mu }^{~\alpha })+F^{\rho \sigma }(\delta _{\rho }^{~\beta }~\delta _{\sigma }^{~\alpha }-\delta _{\sigma }^{~\beta }~\delta _{\rho }^{~\alpha })\right)~=-{\frac {1}{\mu _{0}}}F^{\beta \alpha }$.
  * Par ailleurs : ${\frac {\partial \Lambda }{\partial A_{\alpha }}}=-J^{\alpha }~$ donc $~\partial _{\beta }F^{\beta \alpha }=\mu _{0}~J^{\alpha }$.
  * Pour $~\beta =0~$ on obtient : $~\partial _{\alpha }F^{\alpha 0}=\mu _{0}~J^{0}$.
  
      Mais puisque $~F^{00}=0~$ il reste : $~\partial _{k}F^{k0}={\frac {\partial _{k}E^{k}}{c}}=\mu _{0}~\rho ~c$.
      Avec $~\mu _{0}c^{2}={\frac {1}{\varepsilon _{0}}}~$ ceci donne l'équation de Maxwell-Gauss : $~\mathrm {div} ~{\vec {E}}~=~{\frac {\rho }{\varepsilon _{0}}}$.
  * Pour $~\beta =k=1;2;3~$ on obtient : $~\partial _{0}F^{0k}+\partial _{i}F^{ik}=-{\frac {1}{c^{2}}}{\frac {\partial E^{k}}{\partial t}}+\partial _{i}F^{ik}=\mu _{0}~j^{k}$.
  
      Par ailleurs pour $~k=1~$ (par exemple ; de même pour les autres indices) : $~\partial _{2}F^{21}+\partial _{3}F^{31}=\partial _{2}B^{3}-\partial _{3}B^{2}=({\vec {\mathrm {rot} }}~{\vec {B}})^{1}$.
      Au total ceci donne l'équation de Maxwell-Ampère : $~{\overrightarrow {\mathrm {rot} }}~{\vec {B}}~=~\mu _{0}{\vec {j}}~+~\varepsilon _{0}\mu _{0}~{\frac {\partial {\vec {E}}}{\partial t}}$.
  * En outre $~\partial _{\alpha \beta }F^{\alpha \beta }=-\partial _{\alpha \beta }F^{\alpha \beta }=0~$ puisque $~F^{\alpha \beta }~$ est antisymétrique et que le théorème de Schwarz $~(\partial _{\beta \alpha }=\partial _{\alpha \beta })~$ implique la symétrie.
  
      On en déduit en corollaire l'équation de conservation de la charge :
  
  $\partial _{\alpha }J^{\alpha }=0$
  
  ${\frac {\partial \rho }{\partial t}}+\mathrm {div} ~{\vec {j}}~=0$
  
  ### Tenseur énergie-impulsion du champ électromagnétique
  
  * De façon générale, on peut considérer : $~T_{~\beta }^{\alpha }={\frac {\partial \Lambda }{\partial (\partial _{\alpha }A_{\mu })}}.\partial _{\beta }A_{\mu }-\delta _{~\beta }^{\alpha }\Lambda ~$ avec dans ce cas : $~\Lambda _{em}=-{\frac {1}{4\mu _{0}}}F^{\mu \nu }F_{\mu \nu }$.
  
      On obtient : $~{\frac {\partial \Lambda }{\partial (\partial _{\alpha }A_{\mu })}}=-{\frac {1}{\mu _{0}}}F^{\alpha \mu }~$ ; ainsi : $~T_{\alpha \beta }=-{\frac {1}{\mu _{0}}}\eta ^{\mu \nu }F_{\alpha \mu }.\partial _{\beta }A_{\nu }+{\frac {1}{4\mu _{0}}}\eta _{\alpha \beta }F^{\mu \nu }F_{\mu \nu }$.
  * Ainsi déterminée, l'expression du tenseur énergie-impulsion est généralement dissymétrique ; on peut toutefois montrer qu'elle peut être symétrisée en lui ajoutant la dérivée $~\partial _{\gamma }\Psi ^{\alpha \beta \gamma }~$ d'un tenseur $~\Psi ^{\alpha \beta \gamma }~$ antisymétrique par rapport à $~\alpha \gamma$.
  
      Cela ne change rien à la propriété $~\partial _{\alpha }T^{\alpha \beta }=0~$ puisque $~\partial _{\alpha \gamma }\Psi ^{\alpha \beta \gamma }=0~$ (à la fois symétrique et antisymétrique). On peut en outre montrer que cela ne modifie pas les énergies et impulsions qui s'en déduisent par intégration (les contributions de $~\Psi ^{\alpha \beta \gamma }~$ sont de même nulles).
  * Compte tenu de $~\partial ^{\mu }F_{\alpha \mu }=0~$ (on raisonne pour le champ électromagnétique seul) on peut ajouter $~{\frac {1}{\mu _{0}}}F_{\alpha \mu }.\partial ^{\mu }A_{\beta }={\frac {1}{\mu _{0}}}\partial ^{\mu }(F_{\alpha \mu }A_{\beta })$.
  
      Ceci done une formulation symétrique : $~T_{\alpha \beta }=-{\frac {1}{\mu _{0}}}\eta ^{\mu \nu }F_{\alpha \mu }F_{\beta \nu }+{\frac {1}{4\mu _{0}}}\eta _{\alpha \beta }F^{\mu \nu }F_{\mu \nu }$.
      ♦ remarque : ici encore on rencontre divers choix de conventions de signe ; par exemple les notations anglo-saxonnes changent le signe de la définition de $~\eta _{\alpha \beta }$ ; il faut adapter les signes dans l'expression précédente (et de même si on change le signe de l'action).
  * Ceci peut s'exprimer en considérant : $~F^{\mu \nu }F_{\mu \nu }=F^{0i}F_{0i}+F^{i0}F_{i0}+F^{ij}F_{ij}~$ avec :
  
         $F^{0i}F_{0i}=F^{i0}F_{i0}=-{\frac {{\vec {E}}^{2}}{c^{2}}}$
         $F^{ij}F_{ij}=2{\vec {B}}^{2}~$ (avec les symboles de Levi-Civita et de Kronecker : $~F^{ij}=\varepsilon ^{ijk}B_{k}~$ ; $~F_{ij}=-\varepsilon _{ijl}B^{l}~$ ; $~\varepsilon ^{ijk}\varepsilon _{ijl}=2~\delta _{~l}^{k}~$ ; $~-B_{k}B^{k}={\vec {B}}^{2}$).
      AInsi $~F^{\mu \nu }F_{\mu \nu }=-2~({\frac {{\vec {E}}^{2}}{c^{2}}}-{\vec {B}}^{2})$.
      ♦ remarque : on peut utiliser des indices covariants et contravariants, mais les symboles de Levi-Civita ne se comportent pas tout à fait comme des tenseurs : $~\eta ^{il}\eta ^{jm}\eta ^{kn}\varepsilon _{lmn}=-\varepsilon ^{ijk}$.
  * On en déduit la densité d'énergie : $~w=T^{00}={\frac {1}{\mu _{0}}}{\frac {{\vec {E}}^{2}}{c^{2}}}-{\frac {1}{2~\mu _{0}}}({\frac {{\vec {E}}^{2}}{c^{2}}}-{\vec {B}}^{2})={\frac {\varepsilon _{0}{\vec {E}}^{2}}{2}}+{\frac {{\vec {B}}^{2}}{2\mu _{0}}}$.
  
      De même la densité de « courant d'énergie » (associée au vecteur de poynting) : $~\Pi ^{i}=c~T^{0i}=-{\frac {c}{\mu _{0}}}\eta _{kj}F^{0k}F^{ij}={\frac {1}{\mu _{0}}}E_{j}\varepsilon ^{ijl}B_{l}={\frac {[{\vec {E}}\times {\vec {B}}]^{i}}{\mu _{0}}}$.
  