Mécanique quantique dans l'espace des phases

Le principal avantage de la formulation dans l'espace des phases est qu'elle fait apparaître la mécanique quantique aussi proche que possible de la mécanique hamiltonienne, en évitant le formalisme des opérateurs dans l'espace de Hilbert^[6]. Cette formulation est de nature statistique et offre des connexions logiques entre la mécanique quantique et la mécanique statistique classique, permettant une comparaison naturelle entre les deux (voir limite classique et physique semi-classique). La mécanique quantique dans l'espace des phases est souvent choisie dans certaines applications de l'optique quantique ou dans l'étude de la décohérence et d'une gamme de problèmes techniques spécialisés ; cependant ce formalisme est moins couramment utilisé que celui de Schrödinger par exemple ^[7].

Les idées conceptuelles sous-jacentes au développement de la mécanique quantique dans l'espace des phases se sont ramifiées dans des développements mathématiques tels que la déformation-quantification de Kontsevich (voir la page en anglais formule de quantification de Kontsevich) et la géométrie non commutative.

Summarize Timeline Top Qs Fact Check

La distribution dans l'espace des phases f(x, p) d'un état quantique est l'équivalent de la fonction d'onde dans la représentation de Schrödinger : c'est une distribution de quasi-probabilité («quasi» car elle ne possède pas toutes les propriétés des distributions classiques de probabilités). La distribution dans l'espace des phases peut être traitée comme la description fondamentale et primitive du système quantique, sans aucune référence aux fonctions d'onde ou aux matrices de densité, bien qu'il existe des outils théoriques permettant de passer d'une représentation à l'autre (voir infra).

Il existe plusieurs façons différentes de représenter la distribution, toutes interdépendantes^[8],[9]. La plus remarquable est la distribution de Wigner, W(x, p), découverte en premier^[5]. D'autres représentations (dans l'ordre approximativement décroissant de prévalence dans la littérature) incluent : représentations P de Glauber–Sudarshan ^[10],[11] et Q de Husimi^[12], représentations de Kirkwood–Rihaczek, Mehta, Rivier et Born–Jordan^[13],[14]. Ces alternatives sont plus utiles lorsque l'hamiltonien prend une forme particulière, comme l'ordre normal pour la représentation-P de Glauber-Sudarshan. Étant donné que la représentation de Wigner est la plus courante, cet article s'y tiendra généralement, sauf indication contraire.

La distribution dans l'espace des phases possède des propriétés proches de la densité de probabilité dans un espace des phases à 2 dimensions par degré de liberté du système : soit six dimensions pour une particule isolée dans l'espace classique et 6N dimensions pour un système de N particules. La distribution est à valeur réelle, mais pas nécessairement positive, contrairement à la fonction d'onde qui est généralement à valeur complexe. Nous pouvons interpréter la probabilité de se situer dans un intervalle de position, par exemple, en intégrant la fonction de Wigner sur tous les moments et sur l'intervalle de position :

${\displaystyle \operatorname {P} [a\leq X\leq b]=\int _{a}^{b}\int _{-\infty }^{\infty }W(x,p)\,dp\,dx.}$

Soit Â(x, p) un opérateur représentant une observable, il peut être projeté dans l'espace des phases sous la forme A(x, p), via la transformée de Wigner ; l'opérateur Â(x, p) peut être retrouvé à partir de A(x, p) par la transformée de Weyl (voir l'article en anglais Weyl-Wigner transform).

La valeur moyenne de l'observable par rapport à la distribution de l'espace des phases est ^[3],[15]

${\displaystyle \langle {\hat {A}}\rangle =\int _{-\infty }^{\infty }dp\int _{-\infty }^{\infty }dx\;A(x,p)W(x,p).}$

NB : malgré son apparence et ses usages, W(x, p) n'est pas une véritable distribution de probabilité conjointe, car ses différentes régions ne représentent pas des états mutuellement exclusifs, comme l'exige le troisième axiome de la théorie des probabilités. De plus, elle peut, en général, prendre des valeurs négatives même pour des états purs, en violation du premier axiome ; l'exception notable est la distribution des états cohérents qui est strictement non-négative.

On peut prouver que les régions d'une telle valeur négative sont "petites" : elles ne peuvent pas s'étendre à des régions compactes de volume supérieur à quelques h³, et disparaissent donc dans la limite classique. Elles sont masquées par le principe d'incertitude, qui ne permet pas une localisation précise dans les régions d'espace de phase de volume inférieur à h³ ; ceci rend ainsi de telles "probabilités négatives" moins paradoxales. ${\displaystyle {\hat {A}}}$ dans l'équation ci-dessus, peut être interprété comme une espérance mathématique dans l'espace de Hilbert liée à un opérateur, aussi dans le contexte de l'optique quantique, cette équation est-elle connue sous le nom de théorème d'équivalence optique.

Une approche alternative à la mécanique quantique basée sur l'espace des phases cherche à définir une fonction d'onde (pas seulement une densité de quasi-probabilité) sur l'espace des phases, généralement au moyen de la transformée de Segal-Bargmann. Pour être compatible avec le principe d'incertitude, la fonction d'onde dans l'espace des phases ne peut pas être une fonction arbitraire, sinon elle pourrait être localisée dans une région arbitrairement petite de l'espace des phases ; aussi la transformée de Segal-Bargmann est une fonction holomorphe de ${\displaystyle x+ip}$. Il existe par ailleurs une densité de quasi-probabilité associée à la fonction d'onde dans l'espace des phases ; c'est la représentation Q d'Husimi de la fonction d'onde de position.

Dans la formulation de l'espace des phases, L'opérateur binaire non commutatif fondamental qui remplace l'opérateur standard de produit est le « produit-étoile », représenté par le symbole ★ ^[2]. À chaque représentation de la distribution dans l'espace des phases est associé un produit-étoile caractéristique différent. Dans la suite, sont présentées les propriétés du produit-★ pour la représentation de Wigner-Weyl.

Pour des raisons de commodité de notation, nous introduisons la notion de dérivée « gauche » ${\displaystyle {\overset {_{\gets }}{\partial }}_{x}}$ et de dérivée «droite » ${\displaystyle {\vec {\partial }}_{x}}$. Pour une paire de fonctions f et g, les dérivées gauche et droite sont définies comme

${\displaystyle {\begin{aligned}f{\overset {_{\gets }}{\partial }}_{x}g&={\frac {\partial f}{\partial x}}\cdot g,\\f{\vec {\partial }}_{x}g&=f\cdot {\frac {\partial g}{\partial x}}.\end{aligned}}}$

La définition différentielle du produit-étoile est

${\displaystyle f\star g=f\,\exp {\left({\frac {i\hbar }{2}}\left({\overset {_{\gets }}{\partial }}_{x}{\vec {\partial }}_{p}-{\overset {_{\gets }}{\partial }}_{p}{\vec {\partial }}_{x}\right)\right)}g,}$

où l'argument de la fonction exponentielle peut être interprété comme une série de puissance. Des relations différentielles supplémentaires permettent d'écrire cela en termes de variation des arguments de f et g :

${\displaystyle {\begin{aligned}(f\star g)(x,p)&=f\left(x+{\tfrac {i\hbar }{2}}{\vec {\partial }}_{p},p-{\tfrac {i\hbar }{2}}{\vec {\partial }}_{x}\right)\cdot g(x,p)\\&=f(x,p)\cdot g\left(x-{\tfrac {i\hbar }{2}}{\overset {_{\gets }}{\partial }}_{p},p+{\tfrac {i\hbar }{2}}{\overset {_{\gets }}{\partial }}_{x}\right)\\&=f\left(x+{\tfrac {i\hbar }{2}}{\vec {\partial }}_{p},p\right)\cdot g\left(x-{\tfrac {i\hbar }{2}}{\overset {_{\gets }}{\partial }}_{p},p\right)\\&=f\left(x,p-{\tfrac {i\hbar }{2}}{\vec {\partial }}_{x}\right)\cdot g\left(x,p+{\tfrac {i\hbar }{2}}{\overset {_{\gets }}{\partial }}_{x}\right).\end{aligned}}}$

Il est également possible de définir sous une forme intégrale de convolution^[16], essentiellement par la transformée de Fourier :

${\displaystyle (f\star g)(x,p)={\frac {1}{\pi ^{2}\hbar ^{2}}}\,\int f(x+x',p+p')\,g(x+x'',p+p'')\,\exp {\left({\tfrac {2i}{\hbar }}(x'p''-x''p')\right)}\,dx'dp'dx''dp''.}$

Par exemple, le produit-★ de deux gaussiennes devient :

${\displaystyle \exp {\big (}{-a}(x^{2}+p^{2}){\big )}\star \exp {\big (}{-b}(x^{2}+p^{2}){\big )}={\frac {1}{1+\hbar ^{2}ab}}\exp \left(-{\frac {a+b}{1+\hbar ^{2}ab}}(x^{2}+p^{2})\right),}$

ou celui de la distribution de Dirac

${\displaystyle \delta (x)\star \delta (p)={\frac {2}{h}}\exp \left(2i{\frac {xp}{\hbar }}\right).}$

Les distributions pour des états propres d'énergie sont appelées états propres-★ ou fonctions propres-★. Celles-ci sont obtenues, de manière analogue à l'équation de Schrödinger indépendante du temps, par l'équation aux valeurs propres-★ ^[17],[18]

${\displaystyle H\star W=E\cdot W,}$

où H est l'hamiltonien, une simple fonction de l'espace des phases, le plus souvent identique à l'hamiltonien classique.

L'évolution temporelle de la distribution spatiale des phases est donnée par une modification quantique du théorème de Liouville^[3],[9],[19] (flot hamiltonien). Cette formulation résulte de l'application de la transformation de Wigner à la version quantique de l'équation de Liouville appliquée à la matrice densité, l'équation de von Neumann qui fait intervenir les commutateurs quantiques :

${\displaystyle i\hbar {d \over dt}{\hat {\rho }}=[{\hat {H}},{\hat {\rho }}]\,,}$

Donc, pour une représentation donnée avec son produit-★ associé, l'équation d'évolution temporelle d'une distribution spatiale dans l'espace phases est donnée par

${\displaystyle {\frac {\partial f}{\partial t}}=-{\frac {1}{i\hbar }}\left(f\star H-H\star f\right),}$

ou, pour la fonction de Wigner en particulier,

${\displaystyle {\frac {\partial W}{\partial t}}=-\{\{W,H\}\}=-{\frac {2}{\hbar }}W\sin \left({\frac {\hbar }{2}}({\overset {_{\gets }}{\partial }}_{x}{\vec {\partial }}_{p}-{\overset {_{\gets }}{\partial }}_{p}{\vec {\partial }}_{x})\right)H=-\{W,H\}+O(\hbar ^{2}),}$

où {{ , }} est le crochet de Moyal : transformée de Wigner du commutateur quantique, tandis que { , } est le crochet de Poisson classique^[3].

Ceci donne une illustration concise du principe de correspondance : cette équation se réduit manifestement à l'équation classique de Liouville à la limite ħ → 0. Dans l'extension quantique du flot hamiltonien, la densité de points dans l'espace des phases n'est cependant pas conservée ; le flux de probabilité apparaît "diffusif" et compressible^[3]. La notion de trajectoire quantique est donc une question ouverte ^[20] (voir le film pour le potentiel Morse, ci-dessous, pour apprécier la non-localité du flux de phase quantique).

NB Compte tenu des restrictions imposées par le principe d'incertitude sur la localisation, Niels Bohr a vigoureusement nié l'existence physique de telles trajectoires à l'échelle microscopique. Au moyen de trajectoires formelles dans l'espace des phases, le problème d'évolution temporelle de la fonction de Wigner peut être rigoureusement résolu en utilisant la méthode de l'intégrale de chemin ^[21] et la méthode des caractéristiques quantiques (voir page en anglais ^[22], bien qu'il existe de sérieux obstacles pratiques dans les deux cas.