Équation de Klein-Gordon
# Équation de Klein-Gordon

L'équation de Klein-Gordon, parfois également appelée équation de Klein-Gordon-Fock, est une version relativiste de l'équation de Schrödinger décrivant des particules massives de spin nul, sans ou avec charge électrique, établie indépendamment en 1926 par les physiciens Oskar Klein et Walter Gordon. C'est un exemple d'équation aux dérivées partielles dispersive.

## L'équation de Klein-Gordon

### Dérivation

L'équation de Klein-Gordon standard (sans champs électromagnétique) peut être obtenue de plusieurs façons. Une méthode consiste à écrire une formulation covariante de l'équation d'Euler-Lagrange, et une autre consiste à partir de l'invariant relativiste donnant l'énergie d'une particule isolée telle que :

$E^{2}\ =\ p^{2}\ c^{2}\ +\ m^{2}\ c^{4}=\ {\vec {p}}.{\vec {p}}\ c^{2}\ +\ m^{2}\ c^{4}$

où E est l'énergie totale de la particule, p sa quantité de mouvement, m sa masse propre, et c la célérité de la lumière dans le vide. On applique alors à cette équation énergétique le principe de correspondance de la mécanique quantique : $E\leftrightarrow i\hbar {\frac {\partial }{\partial t}}$ et ${\vec {p}}\leftrightarrow -i\hbar {\frac {\partial }{\partial {\vec {x}}}}=-i\hbar {\vec {\nabla }}$. On obtient alors l'équation dite de Klein-Gordon en faisant agir la relation obtenue sur une fonction d'onde:

$-\ \hbar ^{2}\ {\frac {{\partial }^{2}\Psi ({\vec {r}},t)}{{\partial }t^{2}}}\ =\ -\ \hbar ^{2}\ c^{2}\ \nabla ^{2}\ \Psi ({\vec {r}},t)\ +\ m^{2}\ c^{4}\ \Psi ({\vec {r}},t)$

Cette équation se réécrit sous la forme suivante :

$\nabla ^{2}\ \Psi ({\vec {r}},t)\ -\ {\frac {1}{c^{2}}}\ {\frac {{\partial }^{2}\Psi ({\vec {r}},t)}{{\partial }t^{2}}}\ =\ {\frac {m^{2}\ c^{2}}{\hbar ^{2}}}\ \Psi ({\vec {r}},t)$

On peut également utiliser le formalisme relativiste (en unités naturelles soit $\hbar$ = 1 et c = 1) :

$(\partial _{\mu }\partial ^{\mu }+m^{2})\Psi =0$

avec la convention :

$\partial _{\mu }=\left({\frac {\partial }{\partial t}},{\vec {\nabla }}\right)$

$\partial ^{\mu }=\left({\frac {\partial }{\partial t}},-{\vec {\nabla }}\right)$

### Difficultés d'interprétation

Les solutions de l'équation de Klein-Gordon présentent de sérieuses difficultés d'interprétation dans le cadre de la mécanique quantique originelle, théorie censée décrire une seule particule. Si l'on cherche par exemple à construire une densité de probabilité de présence qui vérifie l'équation relativiste de continuité :

${\frac {\partial \rho ({\vec {r}},t)}{\partial t}}\ +\ {\vec {\nabla }}\cdot {\vec {j}}({\vec {r}},t)\ =\ 0$

on obtient inévitablement les grandeurs suivantes :

$\rho ({\vec {r}},t)\ =\ N\ \Im {\mathfrak {m}}\left(\ {\overline {\Psi }}({\vec {r}},t)\ {\frac {\partial \Psi ({\vec {r}},t)}{\partial t}}\ \right)$

${\vec {j}}({\vec {r}},t)\ =\ Nc^{2}\ \Im {\mathfrak {m}}\left(\ {\overline {\Psi }}({\vec {r}},t)\ {\vec {\nabla }}\ \Psi ({\vec {r}},t)\ \right)$

où ${\overline {\Psi }}$ est le complexe conjugué de $\Psi$, et $N$ est une constante arbitraire. Or, cette densité $\rho$ n'est pas positive partout, donc ne peut représenter une densité de probabilité de présence.

Le cadre pertinent pour interpréter cette équation quantique relativiste sans difficultés est la théorie quantique des champs

## De l'équation de Klein-Gordon à l'équation de Dirac

Le fait que la densité $\rho$ ne soit pas positive partout provient du fait que cette densité contient une dérivée première par rapport au temps, comme l'a remarqué Dirac en 1928. Ceci est lié au fait que l'équation de Klein-Gordon contient une dérivée temporelle seconde.

### Approche naïve

Pour obtenir une équation relativiste du premier ordre en temps, on peut songer à quantifier directement l'expression :

$E\ =\ {\sqrt {p^{2}c^{2}+m^{2}c^{4}}}$

La procédure de quantification canonique conduit alors à l'équation :

$i\,\hbar \ {\frac {\partial }{\partial t}}\Phi ({\vec {r}},t)\ =\ {\sqrt {\ -\,c^{2}\,\hbar ^{2}\,\nabla ^{2}\,+\,m^{2}\,c^{4}\ }}\ \Phi ({\vec {r}},t)$

En raison de la présence d'une racine carrée sur l'opérateur aux dérivées partielles spatiales, cette équation semble a priori bien peu commode à résoudre. On sait aujourd'hui donner un sens mathématiquement précis à l'opérateur ${\sqrt {-\,c^{2}\,\hbar ^{2}\,\nabla ^{2}\,+\,m^{2}\,c^{4}\ }}$ : c'est un opérateur pseudo-différentiel, qui a notamment la particularité d'être non-local — c'est-à-dire que ${\sqrt {-\,c^{2}\,\hbar ^{2}\,\nabla ^{2}\,+\,m^{2}\,c^{4}\ }}\Phi$ dépend des valeurs de $\Phi$ ailleurs que sur un voisinage de ${\vec {r}}$.

### L'équation de Dirac

Dirac recherchera alors une autre équation relativiste du premier ordre en temps et en espace. Il commencera par essayer d'établir une relation de dispersion du type :

$E\ =\ {\vec {\alpha }}\cdot {\vec {p}}\,c\ +\ \beta \ m\,c^{2}$

entre l'énergie, la masse et l'impulsion. Il réussira et, après quantification canonique, obtiendra finalement une équation qui porte aujourd'hui son nom, l'équation de Dirac, et qui décrit très bien les fermions de spin un-demi comme l'électron. Le cadre pertinent pour interpréter cette équation quantique relativiste sans difficultés est encore celui de la théorie quantique des champs.

## Utilisation de l'équation de Klein-Gordon en théorie quantique des champs

Si l'équation de Klein-Gordon ne permet pas de décrire une particule comme l'électron autour du noyau, elle permet de décrire un ensemble de particules de spin 0 dans le cadre de la théorie quantique des champs. L'inconnue Ψ n'est alors plus une fonction d'onde mais un opérateur agissant sur un vecteur d'état appartenant à un espace de Fock. Ce type d'espace décrit l'état d'un système quantique comprenant plusieurs particules dont le nombre peut varier.

La seule particule fondamentale de spin 0 est (en date de septembre 2016) le boson de Higgs. Le champ de Higgs satisfait l'équation de Klein-Gordon.

Il existe des particules composites que l'on peut représenter comme des particules de spin 0 comme les mésons π ou les noyaux de certains atomes comme le carbone 12.

## Aspects historiques

### La querelle des origines

Il est amusant de remarquer que, selon Dirac, Schrödinger aurait d'abord écrit l'équation relativiste dite aujourd'hui de Klein-Gordon, ceci pour tenter de décrire l'électron au sein de l'atome d'hydrogène. En effet, la lecture du premier mémoire de Schrödinger publié en février 1926 montre que celui-ci a déjà essayé une équation d'onde relativiste, mais ce premier mémoire ne contient pas l'équation écrite explicitement. Les prédictions obtenues n'étant pas conformes aux résultats expérimentaux assez précis obtenus par Paschen dès 1916, Schrödinger se serait alors aperçu que c'était l'équation non relativiste — dite aujourd'hui de Schrödinger — qui donnait le bon spectre pour l'hydrogène (après inclusion des effets de spin de façon ad hoc). Schrödinger n'a publié son équation relativiste que dans le quatrième mémoire de 1926.

Entre-temps, plus précisément entre les mois d'avril et septembre 1926, pas moins de cinq autres articles, publiés indépendamment, contenaient l'équation dite aujourd'hui de Klein-Gordon. Les auteurs de ces cinq articles, dont les références figurent dans la bibliographie, sont : Klein, Gordon, Fock, de Donder et van den Dungen, et enfin Kudar.

Enfin, dans son second article de 1926, Fock introduit également la procédure de couplage minimal, décrivant le couplage de la particule massive de charge électrique $e$ à un champ électromagnétique externe donné, représenté par un quadri-potentiel $A_{\mu }$.

### Détails techniques

**Couplage minimal**

Pour une charge $e$ en présence d'un champ électromagnétique externe donné, représenté par le quadri-potentiel $A_{\mu }$, la prescription de couplage minimal de Fock conduit à substituer à la quadri-impulsion $p_{\mu }$ la quantité suivante :

$p_{\mu }\quad \longrightarrow \quad p_{\mu }\,-\,e\,A_{\mu }$

Introduisons explicitement les composantes temporelle et spatiales de la quadri-impulsion :

$p^{\mu }\ =\ {\begin{pmatrix}{\frac {E}{c}}\\{\vec {p}}\end{pmatrix}}$

$p_{\mu }\ =\ {\begin{pmatrix}{\frac {E}{c}}\\-\ {\vec {p}}\end{pmatrix}}$

et du quadri-potentiel :

$A^{\mu }\ =\ {\begin{pmatrix}{\frac {V}{c}}\\{\vec {A}}\end{pmatrix}}$

$A_{\mu }\ =\ {\begin{pmatrix}{\frac {V}{c}}\\-\ {\vec {A}}\end{pmatrix}}$

On obtient alors explicitement :

$p^{\mu }\,-\,e\,A^{\mu }\ =\ {\begin{pmatrix}{\frac {1}{c}}\,\left(\,E\,-\,e\,V\,\right)\\{\vec {p}}\,-\,e\,{\vec {A}}\end{pmatrix}}$

**Spectre de l'atome d'hydrogène d'après le couplage minimal**

On repart de l'équation de dispersion relativiste d'une particule massive isolée :

$E^{2}\ =\ {\vec {p}}^{~2}\ c^{2}\ +\ m^{2}\ c^{4}\quad \Longleftrightarrow \quad p_{\mu }\,p^{\mu }\ =\ m^{2}\,c^{2}$

On y introduit le couplage minimal de Fock :

$\left(p_{\mu }\,-\,e\,A_{\mu }\right)\,\left(p^{\mu }\,-\,e\,A^{\mu }\right)\ =\ m^{2}\,c^{2}$

ce qui donne explicitement en substituant les composantes :

$\left(E\,-\,e\,V\right)^{2}\ -\ \left({\vec {p}}\,-\,e\,{\vec {A}}\right)^{2}c^{2}\ =\ m^{2}\,c^{4}$

En présence d'un potentiel Coulombien statique décrivant l'interaction de l'électron avec un proton (supposé infiniment lourd), le potentiel vecteur est nul : ${\vec {A}}=0$, et on a :

$\left(\,E\,+\,{\frac {e^{2}}{4\pi \epsilon _{0}r}}\,\right)^{2}\ -\ {\vec {p}}^{~2}c^{2}\ =\ m^{2}\,c^{4}$

On applique la quantification canonique à cette équation classique, qui devient un opérateur aux dérivées partielles :

$\left(\,i\,\hbar \,{\frac {\partial ~~}{\partial t}}\,+\,{\frac {e^{2}}{4\pi \epsilon _{0}r}}\,\right)^{2}\ -\ \left(\,-\,i\,\hbar \,{\vec {\nabla }}\,\right)^{2}c^{2}\ =\ m^{2}\,c^{4}$

d'où l'équation dépendant du temps :

$\left(\,i\,\hbar \,{\frac {\partial ~~}{\partial t}}\,+\,{\frac {e^{2}}{4\pi \epsilon _{0}r}}\,\right)^{2}\ \Psi ({\vec {r}},t)\ +\ \hbar ^{2}c^{2}\,\nabla ^{2}\Psi ({\vec {r}},t)\ =\ m^{2}\,c^{4}\ \Psi ({\vec {r}},t)$

On recherche enfin les états stationnaires d'énergie $E$ constante sous la forme d'une fonction purement spatiale multipliée par une exponentielle oscillante en temps :

$\Psi ({\vec {r}},t)\ =\ \psi ({\vec {r}})\ e^{-iEt/\hbar }$

On obtient alors l'équation aux valeurs propres :

$\left(E\,+\,{\frac {e^{2}}{4\pi \epsilon _{0}r}}\,\right)^{2}\ \psi ({\vec {r}})\ +\ \hbar ^{2}c^{2}\,\nabla ^{2}\psi ({\vec {r}})\ =\ m^{2}\,c^{4}\ \psi ({\vec {r}})$

Schrödinger fut découragé par le fait que cette équation ne donne pas le spectre correct de l'atome d'hydrogène. On obtient en effet les niveaux d'énergie suivants :

$E_{n,l}\ =\ mc^{2}\ \left[\,1\,-\,{\frac {\alpha ^{2}}{2\,n^{2}}}\,-\,{\frac {\alpha ^{4}}{2\,n^{4}}}\ \left({\frac {n}{l+{\frac {1}{2}}}}\,-\,{\frac {3}{4}}\right)\,\right]\ +\ O(\alpha ^{6})$

où le nombre quantique principal $n$ est un nombre entier strictement positif, le nombre quantique orbital $l$ est un nombre entier positif compris entre 0 et $n-1$, et $\alpha$ est la constante de structure fine :

$\alpha \ =\ {\frac {e^{2}}{4\pi \epsilon _{0}\hbar c}}\ \simeq \ {\frac {1}{137,04}}$

Le terme d'ordre $\alpha ^{2}$ est correct, mais le terme suivant d'ordre $\alpha ^{4}$, qui décrit la structure fine, n'est pas conforme aux résultats expérimentaux obtenus par Paschen dès 1916. L'expression correcte à cet ordre est en effet :

$E_{n,l}\ =\ mc^{2}\ \left[\,1\,-\,{\frac {\alpha ^{2}}{2\,n^{2}}}\,-\,{\frac {\alpha ^{4}}{2\,n^{4}}}\ \left({\frac {n}{j+{\frac {1}{2}}}}\,-\,{\frac {3}{4}}\right)\,\right]\ +\ O(\alpha ^{6})$

avec le nombre quantique $j=l+1/2$, le demi supplémentaire étant lié au spin de l'électron, qui n'est pas inclus dans l'équation de Klein-Gordon.

