logo-CEA

Introduction

Plan du
cours

Classification Introduction à la physique des

Température
α
Saha
Γ
plasmas
Θ
Coulomb
Debye
ωp
Descriptions
S. Mazevet

Laboratoire de Structure Electronique

Département de Physique Théorique et Appliquée
Commissariat à l’Energie Atomique
Bruyères-Le-Châtel, France

Orsay, Septembre 2009

Orsay, Septembre 2009 p-1/21

 logo-CEA Table of contents

Introduction

Plan du
cours

Classification 1 Introduction
Température
α
Saha
Γ
Θ
2 Plan
Coulomb
Debye
ωp
Descriptions
3 Classification
Température
Degré d’ionisation: α
Saha
Cinétiques-corrélés: Γ
Plasmas classiques/dégénérés

4 Intéractions longues portées

Ecrantage de Debye
Oscillations dans les plasmas: fréquence plasma
Descriptions utilisées pour les plasmas

Orsay, Septembre 2009 p-2/21

 logo-CEA Définition

Introduction

Plan du
cours

Classification Un plasma est défini comme un ensemble de particules chargées

Température
α
Saha
Un plasma est électriquement neutre: nombre égale de particules
Γ
Θ chargées positivement et négativement
Coulomb
Debye Toutes les particules interagissent entre elles via l’interaction
ωp
Descriptions coulombienne
Contrairement à un gaz où les interactions sont de courte portée,
dans un plasma les interactions sont donc de longue portée
Ceci implique un comportement collectif des particules
Modélisable comme deux fluides chargés en interaction
Il existe des plasmas trés variés ±denses,± chauds, ± ionisés
Tous les corps se transforment en plasma lorsque la température
et/ou densité sont suffisement élevées

Orsay, Septembre 2009 p-3/21

 logo-CEA Applications

Introduction

Plan du
cours

Classification
Température
α
Saha
Γ
Θ
Coulomb
Debye
ωp
Descriptions

Astrophysique, 99% de la matière visible, atmosphéres stellaires,

nébuleuses,...
Fusion thermonucléaire controlée: magnétique ou par confinement
inertiel
Décharges électriques dans les gaz,...
Matière dans des conditions extrémes
Orsay, Septembre 2009 p-4/21
 logo-CEA Exemples de Plasmas

Introduction

Plan du
cours

Classification
Température
α
Saha
Γ
Θ
Coulomb
Debye
ωp
Descriptions

Les différents types de plasmas sont catalogués en utilisant quelques

paramètres clés.
Orsay, Septembre 2009 p-5/21
 logo-CEA Plan du cours

Introduction

Plan du
cours

Classification
Température
α
Saha
Γ
Θ Généralités sur les plasmas, grandeurs fondamentales
Coulomb
Debye
Collisions dans les plasmas
ωp
Descriptions
Description fluide, ondes dans les plasmas (2 cours)
Eléments de théorie cinétique (2 cours)
Instabilités paramétriques
Thèmes de recherche actuels sur les plasmas: interaction laser
plasmas, matière dans les conditions extrèmes, simulations.

Orsay, Septembre 2009 p-6/21

 logo-CEA Notion de température

Introduction

Plan du
cours Pour un plasma proche de l’équilibre, les collisions sont suffisement
Classification fréquentes pour que les lois de la mécanique statistique soient
Température
α applicables
Saha
Γ
Θ La distribution des vitesses des particules à une température T est
Coulomb
Debye
donnée par une distribution de Maxwell-Boltzmann (une
ωp
Descriptions
dimension)
1
f (u) = Aexp(− mu2 /kB T ) (1)
2
f du est le nombre de particules par m3 avec une vitesse entre u et
u + du , 12 mu2 est l’énergie cinétique et kB = 1.38 × 10−23 J/K la
constante de Boltzmann
La constante A est reliée à la densité n par
 1/2
m
A=n (2)
2πkB T

Orsay, Septembre 2009 p-7/21

 logo-CEA Notion de température II

Introduction

Plan du L’énergie cinétique moyenne des particules dans cette distribution

cours
est R∞ 1 2
2 mu f (u)du
Classification
Température
α
Eav = −∞ R∞ (3)
Saha
−∞
f (u)du
Γ
Θ
Coulomb En utilisant une intégration par partie, on trouve, à une dimension:
Debye
ωp 1
Descriptions Eav = kB T (4)
2
Ce résultat peut être facilement étendu à trois dimensions en
utilisant
 
1 2 2 2
f (u, v, w) = A3 exp − m(u + v + w )/kb T (5)
2
L’énergie cinétique moyenne est alors Eav = 32 kb T soit 12 kb T par
degré de liberté
Il est courant, en physique des plasmas, d’exprimer la température
en unité d’énergie: 1eV = 11604K
Notion de température pour les électrons Te et les ions Ti
Orsay, Septembre 2009 p-8/21
 logo-CEA Degré d’ionisation

Introduction

Plan du
cours

Classification
L’état d’ionisation d’un plasma est lié à sa température T et sa
Température
α
Saha
densité n
Γ
Θ
Coulomb
Pour un atome A dans un plasma, une collision ionisante est du
Debye
type
ωp
Descriptions e + A → e + A+ + e (6)
A l’équilibre, le plasma contient donc ne électrons, ni ions et n0
neutres par unité de volume.
Pour un plasma globalement neutre: ne = ni = n
Le degré d’ionisation, α est définie par
n
α= (7)
n0 + n

Orsay, Septembre 2009 p-9/21

 logo-CEA Equilibre d’ionisation: Equation de Saha

Introduction

Plan du A cause des collisions, les atomes, molécules, ou ions présents dans
cours

Classification
le plasma peuvent être ionisés si la température est telle que
Température
α
Saha
kB T > Ui /10 (8)
Γ
Θ
Coulomb où Ui est le potentiel d’ionisation.
Debye
ωp
Descriptions
Si le plasma est à l’équilibre thermodynamique, l’ionisation est
contre-balancée par la recombinaison
Cet équilibre est décrit par l’équation de Saha
ni+1 ne gi+1 ge (2πme kB T )3/2
= exp[−(Ui+1 − Ui )/kb T ] (9)
ni gi h3
g sont des facteurs de dégénérescence énergétique et ge = 2
ni est la densité d’atomes dans leur ieme état d’ionisation
h est la constante Planck
(2πme kB T )3/2
h3 correspond à la longeur d’onde thermique d’un
électron
Orsay, Septembre 2009 p-10/21
 logo-CEA Equilibre d’ionisation II

Introduction

Plan du
cours

Classification Le terme qui participe le plus est exp[−Ui /kb T ]

Température
α
Si Ui >> kb T faible ionisation, α → 0
Saha
Γ Si Ui << kb T fortement ionisé α → 1
Θ
Coulomb Application Numérique:
Debye
ωp Azote à température ordinaire n0 = 3.1025 m−3 , Ui = 14.5eV,
Descriptions
T=300K nn0i = 10−26 .
Typiquement, α commence à être significatif lorsque kb T > Ui /10
Permet de distinguer les plasmas faiblement et fortement ionisés
Avec plusieurs espèces (atomes, ions,..), il faut traiter les équations
d’évolution pour toutes les espèces et les mécanismes associés:
Physique Atomique à l’équilibre et hors équilibre
Avec ne = ni = αn, le facteur n déplace l’équation en faveur de la
recombinaison: l’équation de Saha n’est pas valable à forte densité

Orsay, Septembre 2009 p-11/21

 logo-CEA Plasmas cinétiques-corrélés

Introduction

Plan du La distinction entre un plasma cinétique et un plasma corrélé se

cours
fait en comparant l’énergie cinétique à l’énergie d’interaction
Classification
Température
Coulombienne
α
Saha
Γ
Energie cinétique Ecin = 23 kb T
Θ
Coulomb Energie d’interaction Coulombienne
Debye
ωp Z 2 e2
Descriptions
Uint = (10)
4π0 d
avec d distance entre les deux particules
Définition:
→ Ucin >> Uint Plasma cinétique, comportement de type gaz
parfait
→ Ucin << Uint Plasma correlé, les forces électrostatiques
modifient le comportement des particules chargées
Uint
On définit le paramètre de couplage Γ = U cin
→ Γ << 1 plasma cinétique
→ Γ >> 1 plasma corrélé
Orsay, Septembre 2009 p-12/21
 logo-CEA Plasmas cinétiques-corrélés II

Introduction

Plan du
cours On définit la longueur de Landau: longueur d’approche de 2
Classification
particule d’énergie kb T
Température
α
Saha
Γ Z 2 e2
Θ
Coulomb
r0 = (11)
4π0 kb T
Debye
ωp
Descriptions
Le paramètre de couplage devient

Z 2 e2 1 r0
Γ= × = (12)
4π0 d kb T d
La longueur de Landau permet également de classer un plasma
r0 << d → cinétique
r0 >> d → corrélé
Le plasma est corrélé à faible température ou haute densité
Le plasma est cinétique à haute température ou faible densité

Orsay, Septembre 2009 p-13/21

 logo-CEA Plasmas classiques/dégénérés

Introduction

Plan du Si la distance entre 2 particules est petite, les fonctions d’onde se

cours
recouvrent et l’on ne peut plus négliger les effets quantiques
Classification
Température
Ces effets vont se manifester d’abord sur les e car ils ont une
α
Saha fonction d’onde plus étendue
Γ
Θ
Coulomb
L’extention spatiale d’une particule est donnée par la longueur
Debye thermique de de Broglie
ωp
Descriptions
h
λth = √ (13)
mkb T
 1/3 √ r
d 3 mkb T T
= = (14)
λth 4πn h TF
avec TF la température de Fermi
Définition: → d >> λth ou T >> TF : plasma classique, n faible
ou T élevée → d << λth ou T << TF : plasma dégénéré
Lorsque le plasma est dégénéré, les effets quantiques sont
importants. Il faut considérer la statistique de Fermi-Dirac
Les électrons dans un métal est un exemple de plasma dégénéré
Orsay, Septembre 2009 p-14/21
 logo-CEA Intéractions longues portées

Introduction

Plan du
cours

Classification
Dans un gaz parfait, les molécules n’interagissent que si elles
Température
α
Saha
s’approchent à une distance de quelques AA
Γ
Θ
Coulomb
Dans les plasmas, les particules chargées interagissent à longue
Debye portée.
ωp
Descriptions
Le potentiel d’interaction est en 1/r

Q
E= (15)
4π0 r2
Nous allons étudier deux effets importants:
→ Ecrantage de Debye : distance à laquelle la présence d’une
charge se fait sentir
→ Oscillation plasma : fréquence fondamentale pour un plasma

Orsay, Septembre 2009 p-15/21

 logo-CEA Ecrantage de Debye

Introduction

Plan du Si on insère une charge dans un plasma, elle attire les charges
cours

Classification
opposées qui vont l’écranter
Température Si la température du plasma est nulle, l’écrantage sera parfait
α
Saha
Γ Si on insère une grille maintenue au potentiel constant φ0
Θ
Coulomb Les ions sont supposés fixes alors que les électrons s’arrangent en
Debye
ωp fonction de leur température
Descriptions
La fonction de distribution pour les électrons est
1
f (u) = Aexp[−( mu2 + qφ)/kb T ] (16)
2
Intégrer f (u) sur u et en notant que ne (φ → 0) = n∞ , on obtient
ne = n∞ exp(eφ/kb Te ) (17)
En prenant maintenant l’équation de poisson à une dimension
d2 φ
0 ∇2 φ = 0 = −e(ni − ne ) (18)
dx2

Orsay, Septembre 2009 p-16/21

 logo-CEA Ecrantage de Debye II

Introduction

Plan du Utilisant l’expression de ne , l’équation devient

cours

Classification
d2 φ
   
eφ
Température
α 0 2 = en∞ exp −1 (19)
Saha dx kb Te
Γ
Θ
Coulomb
Par linéarisation, on obtient
Debye
ωp
Descriptions
d2 φ n∞ e2
0 2
= φ (20)
dx kb Te
la solutions est du type

φ = φ0 exp(−|x|/λD ) (21)

λD est la longueur de Debye

 1/2
0 KTe
λD = (22)
ne2

Orsay, Septembre 2009 p-17/21

 logo-CEA Ecrantage de Debye: interprétation

Introduction

Plan du Si une charge est plus loin que λd , elle ne voit pas le potentiel φ0
cours

Classification
Lorsque la densité augmente λD diminue
Température
α
Lorsque la température augmente, λD augmente: sans agitation
Saha
Γ
thermique, toutes les charges viennent sur le potentiel
Θ
Coulomb La longueur d’onde de Debye permet de parler de quasi-neutralité
Debye
ωp
pour un plasma.
Descriptions
Si la taille du plasma, L, est telle que L >> λD , on voit que le
plasma sera globalement neutre
Conditions de validité: Il faut que le nombre de particules dans un
volume de Debye soit >> 1:
 3/2
4 3 4 0 KTe
ND = πn0 λD >> 1 → π n0 >> 1 (23)
3 3 ne2
 3/2
d
>> 1 → d >> r0 (24)
r0
Dans un plasma cinétique

Orsay, Septembre 2009 p-18/21

 logo-CEA Oscillations dans les plasmas: fréquence plasma

Introduction

Plan du
cours

Classification
Si les e sont déplacés par rapport aux ions considérés comme fixes,
Température il y a apparition d’un champ électrique qui tend à les ramener à
α
Saha
Γ
l’équilibre et maintenir la neutralité
Θ
Coulomb Les e étant légers, ils se mettent à osciller avec une fréquence
Debye
ωp charactéristique dite fréquence plasma ωp
Descriptions
Ces oscillations sont suffisement rapides pour que les ions n’aient
pas le temps de réagir
En considérant les ions fixes, la température électronique nulle et
en l’absence de champs magnétique, le système à résoudre est:

dve
mne dt = −ene E
 eq. du mouvement
0 ∇.E = 0 ∂E/∂x = e(ni − ne ) eq. de poisson (25)
 ∂ne

∂t + ∇.(ne ve ) = 0 eq. de continuité

Orsay, Septembre 2009 p-19/21

 logo-CEA Oscillations dans les plasmas: fréquence plasma II

Introduction

Plan du On suppose de petites perturbations ne = n0 + n1 , ve = v0 + v1

cours

Classification
et E = E0 + E1
Température Comme le plasma est quasi-neutre et à l’équilibre avant la
α
Saha
Γ
perturbation
Θ
Coulomb
∂n0 ∂v0 ∂E0
Debye ∇n0 = v0 = E0 = 0; = = =0 (26)
ωp
Descriptions
dt ∂t ∂t
Le système devient alors

∂v1
 me ∂t = −eE1

0 ∇E1 = en1 (27)
 ∂n1

∂t + n0 ∇.v1 = 0

Les solutions du type n1 = n1 ei(kx−wt) satisfont

1/2
∂ 2 n1 n0 e2

− ωp2 n1 = 0 avec ωp = (28)
∂t 0 me

Orsay, Septembre 2009 p-20/21

 logo-CEA Descriptions utilisées pour les plasmas

Introduction

Plan du
cours

Classification
Température
α
Saha
Γ
Θ
Une description particulaire (incomplète): équation du
Coulomb mouvement...
Debye
ωp
Descriptions
Théorie cinétique: équation de Boltzman. Particulièrement utilisée
pour déterminer les coefficients de transport
Théorie fluide: le plasma est considéré comme un ou deux fluides
(e+ions). Permet de déterminer un ensemble de données
macroscopiques comme la pression, la température,...
Ces descriptions sont équivalentes et complémentaires

Orsay, Septembre 2009 p-21/21