Équation de Boltzmann quantique — Wikipédia
L'équation de Boltzmann quantique, également connue sous le nom d'équation d'Uehling-Uhlenbeck, est l'extension de l'équation de Boltzmann classique qui donne l'évolution temporelle hors d'équilibre d'un système de particules quantiques en interaction. L'équation de Boltzmann quantique se distingue formellement de celle classique par le terme de collision, qui donne l'évolution de la distribution de quantité de mouvement d'un gaz localement homogène causée par les collisions binaires entre les molécules ; ce terme de collision, qui prend en compte des contraintes de la mécanique quantique, a été initialement formulé par Lothar Nordheim (1928)[1], et par Edwin Uehling (en) et George Uhlenbeck (1933)[2].
L'équation de Boltzmann-Uehling-Uhlenbeck (BUU)
[modifier | modifier le code]Forme générale de l'équation de Boltzmann
[modifier | modifier le code]La forme générale de l'équation de Boltzmann est la suivante :
est la distribution à un corps dans la hiérarchie BBGKY[3] qui donne la densité de probabilité pour une particule quelconque du système de se situer en un point de l'espace des phases. est l'hamiltonien du système ; est le crochet de Poisson.
Lorsque les forces dérivent d'un potentiel, et en explicitant le second terme, on obtient la forme la plus usuelle de l'équation complète
avec le produit scalaire ( . ) et l'opérateur usuel gradient pour les positions, qui se formule de manière analogue pour les moments. représente les forces appliquées sur les particules de la distribution : soit forces extérieures, soit forces entre les particules dans la limite du champ moyen.
est l'opérateur de collision, qui rend compte des interactions directes entre les particules, en général les collisions binaires. La mécanique quantique doit être prise en compte dans la forme exacte de , qui dépend de la physique du système à modéliser [4].
Lorsque l'on traite des systèmes de fermions, le principe d'exclusion de Pauli impose des contraintes fortes sur la densité dans l'espace des phases (voir Physique semi-classique) qui doit satisfaire pour tout volume élémentaire
Dans un régime purement hamiltonien (sans collisions binaires), si la contrainte est vérifiée initialement, elle le sera durant toute l'évolution du système (théorème de Liouville) ; il est donc nécessaire que le terme de collisions prenne en compte cette limitation.
Terme de collision de Uehling et Uhlenbeck (UU)
[modifier | modifier le code]La forme la plus usuelle est le terme de Uehling-Uhlenbeck (UU)[2] pour les fermions :
g est le facteur de dégénérescence des niveaux ( pour des fermions de spin 1/2, pour les nucléons avec l'isospin). Les termes indicés caractérisent les particules en voie d'entrée (indices 1 et 2) et en sortie (indices 3,4) ; la section efficace (probabilité) de collision est . Le terme de collision UU établit la conservation de l'énergie et des impulsions dans un choc élastique (les distributions de Dirac ) et le respect du principe d'exclusion de Pauli.
La fonction de distribution est ici normalisée à un maximum égal à l'unité, c'est donc le nombre d'occupation local dans l'espace des phases.
En factorisant le terme de collisions, on obtient la forme synthétique :
avec le terme de gain qui fait croitre le taux d'occupation local jusqu'à la limite bloquée par le principe de Pauli ; le terme de perte fait décroître le taux d'occupation[5].
Extension au cas des bosons
[modifier | modifier le code]Le terme de collisions ci-dessus, défini pour des fermions, est un cas particulier d'une forme générale qui inclut les particules classiques et les bosons[6]. On peut généraliser les termes de blocage spécifique aux fermions dans l'équation ci-dessus, en prenant en compte respectivement les statistiques classiques, de Bose-Einstein et de Fermi-Dirac : ils s'écrivent alors :
avec pour les particules classiques, pour les bosons et pour les fermions. On obtient alors :
Dans la suite, on se limitera à l'équation pour les fermions.
Blocage de Pauli pour les fermions
[modifier | modifier le code]Dans l'espace des moments, pour conserver l'énergie et l'impulsion, les collisions élastiques entre 2 particules A et B, de moments respectifs et , se traduisent par une rotation aléatoire de leur moment relatif (voir l'article Équation de Boltzmann). Le principe de Pauli en mécanique quantique, implique que pour que la collision soit possible, l'espace des phases pour les impulsions finales et ne soit pas occupé (voir figure). En mécanique classique, la rotation aléatoire peut donc prendre toutes les valeurs de 0 à ; en mécanique quantique, seules les rotations conduisant à des états inoccupés sont possibles et ont donc une probabilité fortement réduite.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/3/34/Boltzmann_quantique_-_Blocage_de_Pauli.png/330px-Boltzmann_quantique_-_Blocage_de_Pauli.png)
Dans le cas de la figure, on considère un système hors équilibre, constitué de deux distributions de fermions identiques (cercles grisés) en interaction (la figure pourrait correspondre à une réaction nucléaire entre 2 noyaux lourds à des énergies supérieures à 70 MeV par nucléon).
Pour la paire de fermions A et B, à gauche, la totalité du cercle (en vert) est accessible avec le terme de collision classique ; seule une très petite partie (en vert à droite) est possible avec le terme de collisions UU ; la probabilité de collisions est réduite d'un ordre de grandeur.
L'équation de Boltzmann quantique s'applique donc à des systèmes fortement hors équilibre (modes collectifs de grande amplitude, collisions…).
Solutions asymptotiques
[modifier | modifier le code]Le terme de collision dans l'équation quantique de Boltzmann donne un comportement irréversible à l'évolution du système. L'équation détermine une solution d'équilibre[5] :
Celle-ci satisfait la statistique de Fermi-Dirac. Elle exprime la probabilité pour les particules d'avoir une énergie en fonction de 2 paramètres, le niveau de Fermi et la température . Elle correspond au maximum de l'entropie du système.
Comme son équivalent classique, l'équation quantique détermine une flèche du temps ; c'est-à-dire qu'après un temps suffisamment long, elle conduit à une distribution d'équilibre qui n'évolue plus macroscopiquement.
Bien que la mécanique quantique soit microscopiquement réversible dans le temps, l'équation quantique de Boltzmann donne un comportement irréversible car l'information de phase est ignorée[7], seul le nombre moyen d'occupation des états quantiques est conservé. La solution de l'équation quantique de Boltzmann est donc une bonne approximation du comportement exact du système sur des échelles de temps courtes par rapport au temps de récurrence de Poincaré, ce qui n'est généralement pas une limitation sévère, car le temps de récurrence de Poincaré peut être plusieurs fois l'âge de l'univers même dans de petits systèmes.
Résolution numérique
[modifier | modifier le code]Comme son équivalent classique, l'équation BUU n'a de solutions analytiques que dans des cas très spécifiques. Dans la pratique, on la résout par des méthodes numériques, soit par des méthodes eulériennes sur réseau, soit par des méthodes lagrangiennes (voir la section « Les méthodes de résolution numérique directe » de l'article Équation de Boltzmann).
Comme elle traite des problèmes hors équilibre avec une grande extension dans l'espace des phases; les méthodes lagrangiennes sont particulièrement adaptées. De plus, ces méthodes permettent une simulation directe du terme de collision UU (voir infra).
Méthodes pseudo-particulaires
[modifier | modifier le code]Elles sont basées sur le fait que la distribution quantique dans l'espace des phases d'un système de A fermions, supposée positive, peut être approchée par celle d'un système de (pseudo-)particules suivant des trajectoires classiques (avec )[8].
Cela revient à projeter la fonction sur une base dépendant du temps , par convolution avec une fonction de poids , indépendante du temps et strictement positive où les fonctions doivent être solutions de l'équation de Boltzmann sans collisions et représenter une base la plus complète possible de l'espace des phases[9]. La solution de l'équation est donnée par les trajectoires classiques des fonctions .
En pratique, on effectue un échantillonnage de type Monte-Carlo et où les sont les coordonnées des pseudo-particules dans l'espace des phases. Selon le contexte, les fonctions g utilisées sont soit des distributions de Dirac, soit des gaussiennes[9]. Le calcul explicite des collisions entre les pseudo-particules nécessite une renormalisation de leur section efficace par rapport à celle des particules physiques.
Applications
[modifier | modifier le code]Physique des semi-conducteurs
[modifier | modifier le code]L'équation de Boltzmann est utilisée en physique des semi-conducteurs[10]. Un modèle typique d'un semi-conducteur peut être construit sur les hypothèses suivantes à la base de l'équation de Boltzmann :
- La distribution des électrons est spatialement homogène à une approximation raisonnable (donc toute dépendance x peut être supprimée) ;
- Le potentiel externe est fonction uniquement de la position et isotrope dans l'espace p, et donc peut être mis à zéro sans perdre davantage de généralité ;
- Le gaz est suffisamment dilué pour que les interactions à trois corps entre les électrons puissent être ignorées.
Physique nucléaire
[modifier | modifier le code]À partir des années 1980, avec le développement des accélérateurs permettant d'étudier les réactions entre noyaux lourds de haute énergie, une part importante de la physique nucléaire a été consacrée à l'étude de la matière nucléaire très loin de l'équilibre ; l'étude théorique de ces réactions est très largement basée sur les équations du type Boltzmann quantique[11].
![Exemple d'utilisation de l'équation de Boltzmann quantique en physique nucléaire (réaction très inélastique Sn+Xe à 35MeV/u). Les courbes de niveau représentent la densité nucléaire projetée sur le plan de réaction ; en couleurs pour le projectile, pour montrer les échanges de nucléons lors de la réaction.](http://upload.wikimedia.org/wikipedia/commons/8/8e/Simulation_%28Boltzmann_quantique%29_de_r%C3%A9actions_nucl%C3%A9aires_entre_noyaux_lourds.gif)
La figure ci-contre donne un exemple d'utilisation de l'équation de Boltzmann quantique en physique nucléaire (réaction très inélastique Sn+Xe à 35 MeV/u).
Les courbes de niveau représentent la densité de probabilité des nucléons projetée sur le plan de réaction (cadre de 80 x 60 fermis) ; les courbes sont tracées en couleurs pour le projectile, pour montrer les échanges de nucléons lors de la réaction. L'unité de temps est le fm/c : . La solution de l'équation de Boltzmann donne accès aux échanges de nucléons et d'énergie dans la réaction et aux trajectoires des fragments produits.
Physique relativiste
[modifier | modifier le code]Aux énergies relativistes pour les systèmes baryoniques, il est nécessaire : d'une part de rendre l'équation invariante par la transformation de Lorentz et d'autre part, de tenir compte de la production de nouvelles particules dont les seuils de production sont dépassés[12] (collisions nucléaires à très haute énergie[13], milieux stellaire...).
Références
[modifier | modifier le code]- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Quantum Boltzmann equation » (voir la liste des auteurs).
- L. Nordhiem W. et R. W. Fowler, « On the kinetic method in the new statistics and application in the electron theory of conductivity », Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character, vol. 119, no 783, , p. 689–698 (DOI 10.1098/rspa.1928.0126, Bibcode 1928RSPSA.119..689N)
- (en) Uehling et Uhlenbeck, « Transport Phenomena in Einstein-Bose and Fermi-Dirac Gases. I », Physical Review, vol. 43, no 7, , p. 552–561 (ISSN 0031-899X, DOI 10.1103/PhysRev.43.552, Bibcode 1933PhRv...43..552U, lire en ligne)
- Christian Grégoire, « Effets des collisions sur la dynamique nucléaire » [PDF], sur HAL - Open science, (consulté le )
- Filbert, Hu et Jin, « A Numerical Scheme for the Quantum Boltzmann Equation Efficient in the Fluid Regime », ESAIM: M2AN, vol. 46, no 2, , p. 443–463 (DOI 10.1051/m2an/2011051, arXiv 1009.3352)
- (en) Jørgen Randrup et Bernard Remaud, « Fluctuations in one-body dynamics », Nuclear Physics A, vol. 514, no 2, , p. 339–366 (ISSN 0375-9474, DOI 10.1016/0375-9474(90)90075-W, lire en ligne, consulté le )
- (en) Bao, Markowich et Pareschi, « Quantum kinetic theory: modelling and numerics for Bose-Einstein condensation », Modeling and Computational Methods for Kinetic Equations, Series: Modeling and Simulation in Science, Engineering and Technology, modeling and Simulation in Science, Engineering and Technology, , p. 287–320 (ISBN 978-1-4612-6487-3, DOI 10.1007/978-0-8176-8200-2_10)
- Snoke, Liu et Girvin, « The basis of the Second Law of thermodynamics in quantum field theory », Annals of Physics, vol. 327, no 7, , p. 1825–1851 (DOI 10.1016/j.aop.2011.12.016, Bibcode 2012AnPhy.327.1825S, arXiv 1112.3009, S2CID 118666925)
- (en) Bernard Remaud, « Pseudo-Particle Models for Nuclear Reactions with Heavy Ions », Annales de Physique, vol. 21, no 5, , p. 503–535 (ISSN 0003-4169 et 1286-4838, DOI 10.1051/anphys:199605003, lire en ligne, consulté le )
- (en) P. Schuck, R. W. Hasse, J. Jaenicke et C. Grégoire, « Semiclassical and phase space approaches to dynamic and collisional problems of nuclei », Progress in Particle and Nuclear Physics, vol. 22, , p. 181–278 (ISSN 0146-6410, DOI 10.1016/0146-6410(89)90004-5, lire en ligne, consulté le )
- Snoke, « The quantum Boltzmann equation in semiconductor physics », Annalen der Physik, vol. 523, nos 1–2, , p. 87–100 (DOI 10.1002/andp.201000102, Bibcode 2011AnP...523...87S, arXiv 1011.3849, S2CID 119250989)
- (en) G. F. Bertsch et S. Das Gupta, « A guide to microscopic models for intermediate energy heavy ion collisions », Physics Reports, vol. 160, no 4, , p. 189–233 (ISSN 0370-1573, DOI 10.1016/0370-1573(88)90170-6, lire en ligne, consulté le )
- H.-Th. Elze, M. Gyulassy, D. Vasak et Hannelore Heinz, « Towards a relativistic selfconsistent quantum transport theory of hadronic matter », Modern Physics Letters A, vol. 02, no 07, , p. 451–460 (ISSN 0217-7323, DOI 10.1142/S0217732387000562, lire en ligne, consulté le )
- Q. Li, J. Q. Wu et C. M. Ko, « Relativistic Vlasov-Uehling-Uhlenbeck equation for nucleus-nucleus collisions », Physical Review C, vol. 39, no 3, , p. 849–852 (DOI 10.1103/PhysRevC.39.849, lire en ligne, consulté le )