Équation de Boltzmann quantique

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].

La forme générale de l'équation de Boltzmann est la suivante :

${\displaystyle {\displaystyle {\frac {df(\mathbf {r} ,\mathbf {p} ;t)}{dt}}+\{f(\mathbf {r} ,\mathbf {p} ;t),H(\mathbf {r} ,\mathbf {p} )\}=\left({\frac {\partial f}{\partial t}}\right)_{coll}}}$

${\displaystyle f(\mathbf {r} ,\mathbf {p} \;;t)}$ 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 ${\displaystyle (\mathbf {r} ,\mathbf {p} )}$ de l'espace des phases. ${\displaystyle H(\mathbf {r} ,\mathbf {p} )}$ est l'hamiltonien du système ; ${\displaystyle \{\,,\,\}}$ 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

${\displaystyle \left[{\frac {\partial }{\partial t}}+\mathbf {v} \cdot \nabla _{r}+\mathbf {F} \cdot \nabla _{p}\right]f(\mathbf {r} ,\mathbf {p} ,t)=\mathbf {I} _{coll}[f](\mathbf {r} ,\mathbf {p} )}$

avec le produit scalaire ( . ) et ${\displaystyle \nabla _{r}}$ l'opérateur usuel gradient pour les positions, qui se formule de manière analogue pour les moments. ${\displaystyle \mathbf {F} }$ 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.

${\displaystyle \mathbf {I} _{coll}[f]}$ 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 ${\displaystyle \mathbf {I} _{coll}[f]}$, 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 ${\displaystyle \Delta \mathbf {r} \,\Delta \mathbf {p} }$

${\displaystyle f(\mathbf {r} ,\mathbf {p} ,t)\,\Delta \mathbf {r} \,\Delta \mathbf {p} <{\frac {2}{\hbar ^{3}}}\,\Delta \mathbf {r} \,\Delta \mathbf {p} }$

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.

La forme la plus usuelle est le terme de Uehling-Uhlenbeck (UU)^[2] pour les fermions : $${\displaystyle \left({\frac {\partial f(p_{1})}{\partial f}}\right)_{coll}={\frac {g}{2\pi ^{2}}}{\frac {h^{3}}{m^{2}}}\int ...\int d\mathbf {p_{2}} \,d\mathbf {p_{3}} \,d\mathbf {p_{4}} \,\delta (\epsilon _{1}^{2}+\epsilon _{2}^{2}-\epsilon _{3}^{2}-\epsilon _{4}^{2})\,\delta (\mathbf {p_{1}} +\mathbf {p_{2}} -\mathbf {p_{3}} -\mathbf {p_{4}} ){\frac {d\sigma }{d\Omega }}[(1-f_{1})(1-f_{2})f_{3}f_{4}-(1-f_{3})(1-f_{4})f_{1}f_{2}]}$$

g est le facteur de dégénérescence des niveaux (${\displaystyle g=2}$ pour des fermions de spin 1/2, ${\displaystyle g=4}$ pour les nucléons avec l'isospin). Les termes indicés ${\displaystyle f_{i}\equiv f(\mathbf {r} _{i},\mathbf {p} _{i};t)\;;\;\epsilon _{i}=\mathbf {p} _{i}^{2}/2m}$ 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 ${\displaystyle {\displaystyle d\sigma /d\Omega }}$. Le terme de collision UU établit la conservation de l'énergie et des impulsions dans un choc élastique (les distributions de Dirac ${\displaystyle \delta }$) 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 :

${\displaystyle {\displaystyle \left({\frac {\partial f}{\partial t}}\right)_{coll}=W^{+}(1-f)-W^{-}\,f}}$

avec le terme de gain ${\displaystyle W^{+}}$ qui fait croitre le taux d'occupation local jusqu'à la limite ${\displaystyle f\rightarrow 1}$ bloquée par le principe de Pauli ; le terme de perte ${\displaystyle W^{-}}$fait décroître le taux d'occupation^[5].

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 :

${\displaystyle (1+\theta f_{1})(1+\theta f_{2})f_{3}f_{4}-(1+\theta f_{3})(1+\theta f_{4})f_{1}f_{2}}$

avec ${\displaystyle \theta =0}$ pour les particules classiques, ${\displaystyle \theta =1}$ pour les bosons et ${\displaystyle \theta =-1}$ pour les fermions. On obtient alors :

${\displaystyle {\displaystyle \left({\frac {\partial f}{\partial t}}\right)_{coll}=W^{+}(1+\theta f)-W^{-}f}}$

Dans la suite, on se limitera à l'équation pour les fermions.

Dans l'espace des moments, pour conserver l'énergie et l'impulsion, les collisions élastiques entre 2 particules A et B, de moments respectifs ${\displaystyle \mathbf {p} _{A}}$ et ${\displaystyle \mathbf {p} _{B}}$, se traduisent par une rotation aléatoire de leur moment relatif ${\displaystyle \mathbf {p} _{A}-\mathbf {p} _{B}}$ (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 ${\displaystyle \mathbf {p} '_{A}}$ et ${\displaystyle \mathbf {p} '_{B}}$ ne soit pas occupé (voir figure). En mécanique classique, la rotation aléatoire peut donc prendre toutes les valeurs de 0 à ${\displaystyle 2\pi }$ ; en mécanique quantique, seules les rotations conduisant à des états inoccupés sont possibles et ont donc une probabilité fortement réduite.

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…).

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] :

${\displaystyle f_{equ}={\frac {W^{+}}{W^{+}+W^{-}}}={\frac {1}{1+\exp \left({\frac {\epsilon -\mu }{T}}\right)}}}$

Celle-ci satisfait la statistique de Fermi-Dirac. Elle exprime la probabilité pour les particules d'avoir une énergie ${\displaystyle \epsilon }$ en fonction de 2 paramètres, le niveau de Fermi ${\displaystyle \mu }$ et la température ${\displaystyle T}$. 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.

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).

Elles sont basées sur le fait que la distribution quantique dans l'espace des phases ${\displaystyle \textstyle f(\mathbf {r} ,\mathbf {p} ,t)\geq 0}$ d'un système de A fermions, supposée positive, peut être approchée par celle d'un système de ${\displaystyle \textstyle nA}$ (pseudo-)particules suivant des trajectoires classiques (avec ${\displaystyle \textstyle n\gg 1}$)^[8].

Cela revient à projeter la fonction sur une base dépendant du temps ${\displaystyle g(\mathbf {r} ,\mathbf {p} ,t)}$, par convolution avec une fonction de poids ${\displaystyle w(\mathbf {r} ,\mathbf {p} )}$, indépendante du temps et strictement positive ${\displaystyle f(\mathbf {r} ,\mathbf {p} ,t)=w(\mathbf {r} ,\mathbf {p} )*g(\mathbf {r} ,\mathbf {p} ,t)\,,}$ où les fonctions ${\displaystyle g(\mathbf {r} ,\mathbf {p} ,t)}$ 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 ${\displaystyle g(\mathbf {r} ,\mathbf {p} ,t)}$.

En pratique, on effectue un échantillonnage de type Monte-Carlo et ${\displaystyle \textstyle f(\mathbf {r} ,\mathbf {p} ,t)\approx {\frac {1}{n}}\sum _{i=1}^{nA}g(\mathbf {r} -\mathbf {r} _{i}(t))\;g(\mathbf {p} -\mathbf {p} _{i}(t))\,,}$ où les ${\displaystyle (\mathbf {r_{i}} ,\mathbf {p_{i}} )}$ 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.

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 :

1. La distribution des électrons est spatialement homogène à une approximation raisonnable (donc toute dépendance x peut être supprimée) ;
2. Le potentiel externe est fonction uniquement de la position et isotrope dans l'espace p, et donc ${\displaystyle \mathbf {F} }$ peut être mis à zéro sans perdre davantage de généralité ;
3. Le gaz est suffisamment dilué pour que les interactions à trois corps entre les électrons puissent être ignorées.

À 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].

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 : ${\displaystyle 3\,fm/c=10^{-23}s}$. 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.

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...).

• (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).

• Équation de Boltzmann
• Mécanique quantique dans l'espace des phases
• Espace des phases
• Mécanique hamiltonienne

• Portail de la physique