Mécanique quantique dans l'espace des phases

La formulation de la mécanique quantique dans l'espace des phases place les variables de position et d'impulsion sur un pied d'égalité dans l'espace des phases. En revanche, la représentation de Schrödinger utilise soit la représentation dans l'espace des positions, soit la représentation dans celui des impulsions (voir la page espace des positions et des impulsions)[1]. Les deux principales caractéristiques de la formulation de la mécanique quantique dans l'espace des phases sont que l'état quantique est décrit par une distribution de quasi-probabilité (au lieu d'une fonction d'onde, d'un vecteur d'état ou d'une matrice de densité ) et que la règle de composition des opérateurs est redéfinie.

La théorie a été entièrement développée par Hilbrand Groenewold en 1946 dans sa thèse de doctorat [2], et indépendamment par Joe Moyal[3] ; chacun s'appuyant sur des idées antérieures d'Hermann Weyl [4] et d'Eugene Wigner[5].

Principes de base modifier

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.

Distribution dans l'espace des phases modifier

La distribution dans l'espace des phases f(xp) 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(xp), 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 :

 

Soit Â(xp) 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 Â(xp) 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]

 

NB : malgré son apparence et ses usages, W(xp) 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 h3, 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 à h3 ; ceci rend ainsi de telles "probabilités négatives" moins paradoxales.   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  . 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.

Produit-étoile (-★) d'opérateurs modifier

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 »   et de dérivée «droite »  . Pour une paire de fonctions f et g, les dérivées gauche et droite sont définies comme

 

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

 

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 :

 

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

 

Par exemple, le produit- de deux gaussiennes devient :

 

ou celui de la distribution de Dirac

 

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]

 

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

Évolution temporelle modifier

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 :

 

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

 

ou, pour la fonction de Wigner en particulier,

 

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.

Exemples modifier

Oscillateur harmonique simple modifier

 
La distribution de quasi-probabilité de Wigner Fn(u) pour l'oscillateur harmonique simple avec a) n = 0, b) n = 1, c) n = 5

L'hamiltonien pour l'oscillateur harmonique simple à une dimension dans la représentation de Wigner – Weyl est

 

L'équation aux valeurs propres- pour la fonction de Wigner indépendante du temps, se lit alors

 

Considérons d'abord la partie imaginaire de l'équation aux valeurs propres-★,

 

Cela implique que l'on peut écrire les valeurs propres- comme des fonctions d'un seul argument :

 

Avec ce changement de variables, il est possible d'écrire la partie réelle de l'équation aux valeurs propres- sous la forme d'une équation de Laguerre modifiée (et non pas l'équation d'Hermite !), dont la solution introduite par Groenewold[2] fait intervenir les polynômes de Laguerre [18]

 

avec les valeurs propres- associées

 

Pour l'oscillateur harmonique, l'évolution temporelle d'une distribution de Wigner arbitraire est simple. Une distribution initiale W(x, p; t = 0) = F(u) évolue selon l'équation d'évolution ci-dessus, pilotée par l'hamiltonien de l'oscillateur ; c'est une simple rotation globale dans l'espace des phases[2],

 
 
Fonction de Wigner d'un état cohérent (état fondamental déplacé par rapport à l'origine de l'espace des phases) pour l'oscillateur harmonique (voir le lien)

Typiquement, un état cohérent ("forme en cloche") d'un état quantique d'énergie Eħω peut aussi représenter une quantité macroscopique et apparaître comme un objet classique (un simple oscillateur mécanique) tournant uniformément dans l'espace des phases(voir les figures animées). L'intégration sur toutes les positions de départ à t = 0 de tels objets donne une configuration indépendante du temps similaire aux états propres- statiques F(u) ci-dessus ; elle fournit donc une visualisation intuitive de la limite classique pour les systèmes à grande action[7].

Moment angulaire des particules libres modifier

Supposons qu'une particule se trouve initialement dans un état gaussien de produit d'incertitude minimum, avec les valeurs moyennes de position et d'impulsion toutes deux centrées à l'origine dans l'espace des phases. La fonction de Wigner pour un tel état se propageant librement est

 

α est un paramètre décrivant la largeur initiale de la gaussienne, et τ = m/α2ħ .

Initialement, la position et l'impulsion ne sont pas corrélées (la valeur moyenne du produit position-impulsion est nulle).

Cependant, au fur et à mesure que l'état évolue, la position et l'impulsion deviennent de plus en plus corrélées, car les parties de la distribution plus éloignées de l'origine en position nécessitent une plus grande quantité de mouvement pour être atteintes : et asymptotiquement,

 

qui traduit l'étalement du paquet d'onde libre dans l'espace des coordonnées).

En effet, il est possible de montrer que l'énergie cinétique K de la particule devient asymptotiquement uniquement radiale, en accord avec la notion standard de mécanique quantique de l'indépendance d'orientation du moment angulaire (non nul) d'un état fondamental, soit [23]:

 
 

Potentiel de Morse modifier

Le potentiel de Morse est utilisé pour décrire approximativement la structure vibrationnelle d'une molécule diatomique.

L'évolution au cours du temps de la fonction de Wigner pour le potentiel de Morse U(x) = 20(1 − e−0.16x)2 en unités atomiques (a.u.). Les lignes solides représentent les courbes de niveau du Hamiltonien H(x, p) = p2/2 + U(x).

Effet tunnel quantique modifier

L'effet tunnel est un effet quantique caractéristique où une particule quantique, n'ayant pas suffisamment d'énergie pour passer au-dessus d'une barrière, a la possibilité de la traverser. Cet effet n'existe pas en mécanique classique.

La fonction de Wigner pour le passage par effet tunnel à travers une barrière de potentiel U(x) = 8e−0.25x2 en unités atomiques (a.u.). Les lignes solides représentent les courbes de niveau du Hamiltonien H(x, p) = p2/2 + U(x).

Potentiel quartique modifier

Lorsque le potentiel n'est pas quadratique (ici il est quartique), la fonction de Wigner d'un état initial localisé va se disperser dans l'espace des phases.

L'évolution au cours du temps de la fonction de Wigner pour le potentiel quartique U(x) = 0.1x4 en unités atomiques (a.u.). Les lignes solides représentent les courbes de niveau du Hamiltonien H(x, p) = p2/2 + U(x).

Interférence de 2 états dans un puits harmonique modifier

Évolution de deux états cohérents dans un même puits de potentiel harmonique.

 
Fonction de Wigner de deux états cohérents interférents évoluant selon l'hamiltonien de l'oscillateur harmonique simple. Les projections d'impulsion et de coordonnées correspondantes sont tracées à droite et sous le tracé de l'espace des phases.

Notes et références modifier

  1. Dans la littérature, on trouvera souvent que les concepts d'impulsion, de moment linéaire, de quantité de mouvement sont employés de manière équivalente ; c'est justifié lorsque les forces dérivent d'une énergie potentielle (voir l'article quantité de mouvement). Le concept de moment est plus général ; en mécanique hamiltonienne c'est la variable conjuguée de la variable position.
  2. a b c et d Groenewold, « On the principles of elementary quantum mechanics », Physica, vol. 12, no 7,‎ , p. 405–460 (DOI 10.1016/S0031-8914(46)80059-4, Bibcode 1946Phy....12..405G)
  3. a b c d et e Moyal et Bartlett, « Quantum mechanics as a statistical theory », Mathematical Proceedings of the Cambridge Philosophical Society, vol. 45, no 1,‎ , p. 99–124 (DOI 10.1017/S0305004100000487, Bibcode 1949PCPS...45...99M)
  4. (de) Weyl, « Quantenmechanik und Gruppentheorie », Zeitschrift für Physik, vol. 46, nos 1–2,‎ , p. 1–46 (DOI 10.1007/BF02055756, Bibcode 1927ZPhy...46....1W)
  5. a et b Wigner, « On the Quantum Correction for Thermodynamic Equilibrium », Physical Review, vol. 40, no 5,‎ , p. 749–759 (DOI 10.1103/PhysRev.40.749, Bibcode 1932PhRv...40..749W, hdl 10338.dmlcz/141466)
  6. Ali et Engliš, « Quantization Methods: A Guide for Physicists and Analysts », Reviews in Mathematical Physics, vol. 17, no 4,‎ , p. 391–490 (DOI 10.1142/S0129055X05002376, arXiv math-ph/0405065)
  7. a et b Curtright et Zachos, « Quantum Mechanics in Phase Space », Asia Pacific Physics Newsletter, vol. 01,‎ , p. 37–46 (DOI 10.1142/S2251158X12000069, arXiv 1104.5269)
  8. Cohen, « Generalized Phase-Space Distribution Functions », Journal of Mathematical Physics, vol. 7, no 5,‎ , p. 781–786 (DOI 10.1063/1.1931206, Bibcode 1966JMP.....7..781C)
  9. a et b Agarwal et Wolf, « Calculus for Functions of Noncommuting Operators and General Phase-Space Methods in Quantum Mechanics. II. Quantum Mechanics in Phase Space », Physical Review D, vol. 2, no 10,‎ , p. 2187–2205 (DOI 10.1103/PhysRevD.2.2187, Bibcode 1970PhRvD...2.2187A)
  10. Sudarshan, « Equivalence of Semiclassical and Quantum Mechanical Descriptions of Statistical Light Beams », Physical Review Letters, vol. 10, no 7,‎ , p. 277–279 (DOI 10.1103/PhysRevLett.10.277, Bibcode 1963PhRvL..10..277S)
  11. Glauber, « Coherent and Incoherent States of the Radiation Field », Physical Review, vol. 131, no 6,‎ , p. 2766–2788 (DOI 10.1103/PhysRev.131.2766, Bibcode 1963PhRv..131.2766G)
  12. Kôdi Husimi (1940). "Some Formal Properties of the Density Matrix", Proc. Phys. Math. Soc. Jpn. 22: 264–314.
  13. Agarwal et Wolf, « Calculus for Functions of Noncommuting Operators and General Phase-Space Methods in Quantum Mechanics. I. Mapping Theorems and Ordering of Functions of Noncommuting Operators », Physical Review D, vol. 2, no 10,‎ , p. 2161–2186 (DOI 10.1103/PhysRevD.2.2161, Bibcode 1970PhRvD...2.2161A)
  14. Cahill et Glauber, « Ordered Expansions in Boson Amplitude Operators », Physical Review, vol. 177, no 5,‎ , p. 1857–1881 (DOI 10.1103/PhysRev.177.1857, Bibcode 1969PhRv..177.1857C, lire en ligne); Cahill et Glauber, « Density Operators and Quasiprobability Distributions », Physical Review, vol. 177, no 5,‎ , p. 1882–1902 (DOI 10.1103/PhysRev.177.1882, Bibcode 1969PhRv..177.1882C, lire en ligne).
  15. Lax, « Quantum Noise. XI. Multitime Correspondence between Quantum and Classical Stochastic Processes », Physical Review, vol. 172, no 2,‎ , p. 350–361 (DOI 10.1103/PhysRev.172.350, Bibcode 1968PhRv..172..350L)
  16. Baker, « Formulation of Quantum Mechanics Based on the Quasi-Probability Distribution Induced on Phase Space », Physical Review, vol. 109, no 6,‎ , p. 2198–2206 (DOI 10.1103/PhysRev.109.2198, Bibcode 1958PhRv..109.2198B)
  17. Fairlie, « The formulation of quantum mechanics in terms of phase space functions », Mathematical Proceedings of the Cambridge Philosophical Society, vol. 60, no 3,‎ , p. 581–586 (DOI 10.1017/S0305004100038068, Bibcode 1964PCPS...60..581F)
  18. a et b Curtright, Fairlie et Zachos, « Features of time-independent Wigner functions », Physical Review D, vol. 58, no 2,‎ , p. 025002 (DOI 10.1103/PhysRevD.58.025002, Bibcode 1998PhRvD..58b5002C, arXiv hep-th/9711183)
  19. Mehta, « Phase‐Space Formulation of the Dynamics of Canonical Variables », Journal of Mathematical Physics, vol. 5, no 5,‎ , p. 677–686 (DOI 10.1063/1.1704163, Bibcode 1964JMP.....5..677M)
  20. M. Oliva, D. Kakofengitis, O. Steuernagel, « Anharmonic quantum mechanical systems do not feature phase space trajectories », Physica A, vol. 502,‎ , p. 201–210 (DOI 10.1016/j.physa.2017.10.047, Bibcode 2018PhyA..502..201O, arXiv 1611.03303)
  21. Marinov, « A new type of phase-space path integral », Physics Letters A, vol. 153, no 1,‎ , p. 5–11 (DOI 10.1016/0375-9601(91)90352-9, Bibcode 1991PhLA..153....5M)
  22. Krivoruchenko et Faessler, « Weyl's symbols of Heisenberg operators of canonical coordinates and momenta as quantum characteristics », Journal of Mathematical Physics, vol. 48, no 5,‎ , p. 052107 (DOI 10.1063/1.2735816, Bibcode 2007JMP....48e2107K, arXiv quant-ph/0604075)
  23. Dahl et Schleich, « Concepts of radial and angular kinetic energies », Physical Review A, American Physical Society (APS), vol. 65, no 2,‎ , p. 022109 (ISSN 1050-2947, DOI 10.1103/physreva.65.022109, arXiv quant-ph/0110134)

Articles connexes modifier