Cours 5

8.2.
THÉORIE DES PERTURBATIONS DÉPENDANT DU TEMPS 325
Pab(t)
si Ω=ω−ωab = 0
si Ω=ω−ωab = 0
~1/ Ω 2
j
t
Figure 8.2.5 Allure de la probabilité de transition Pa→b
(1)
(t), à Ω = ω − ωa,b xé, selon
que Ω 6= 0 ou Ω = 0. Dans ce dernier cas, le résultat est valable tant que Pab (t) 1.
Remarque
◦ Si on traitait le champ par la mécanique quantique [CTDRG87, CTDRG88], dans
le cadre de l'électrodynamique quantique (Q.E.D.), on montrerait dans le cas 1) que
l'atome a absorbé un photon d'énergie ~ω , et dans le cas 2) on montrerait
que l'atome a émit un photon d'énergie ~ω . C'est le phénomène d'émission
stimulée.
◦ Les résultats que nous venons d'obtenir sont très analogues aux phénomène de
résonance en mécanique classique linéaire (résonances et battements). C'est le même
phénomène ici retrouvé dans le cadre des états quantiques.
◦ On a |Wab | = |Wba | montrant que dans cette description, le processus d'absorption
et émission stimulée ont la même probabilité. En réalité ce n'est par vrai, et la
diérence se calcule dans le cadre de la Q.E.D. qui considère le champ comme un
système dynamique.
Néanmoins, il est possible de calculer cette diérence sans rentrer dans la complexité
de la Q.E.D.. C'est l'argument d'Einstein qui utilise des considérations de physique
statistique pour la dynamique du champ et de l'atome. Il montre alors la nécessité
du phénomène d'émission spontanée qui est que l'atome peut transiter de Eb vers
Ea < Eb en émettant un photon , même en l'absence de champ extérieur.
Ce phénomène a un équivalent classique : la radiation de Larmor d'une charge
accélérée.
Cf Bransden p499-501. @@ (à rédiger ici).
8.2.2.4 Expression de |Wba |2

La probabilité de transition Pab (t) eq(8.2.8), dépend de l'élément de matrice :
e ~
Wba = hψb | ~xˆ.E~k,ω |ψa i
2
326 CHAPITRE 8. MÉTHODES D'APPROXIMATION ; RÉSOLUTION APPROCHÉE
On note ~ ba = −ehψb |~xˆψa i

D les éléments de matrice de l'opérateur moment dipolaire
électrique D ~ = −e ~x, et θ l'angle entre ~ ba
D ~ . L'intensité
et E de l'onde incidente I est
obtenue par la norme du vecteur de Poynting ~ ~ ∧B
S = c2 ε 0 E ~ :
0
D
~
E ε0 c ~ 2
0
I = ux d énergie incident = S = Eω : intensité de l onde incidente
2
Alors :
1 ~ 2 2
|Wba |2 = E~k,ω D~ ba cos2 θ (8.2.9)
4
I ~ 2
= Dba cos2 θ (8.2.10)
2ε0 c
Le résultat (8.2.10) montre que l'inuence du champ E ~ sur l'atome est directive à cause
2
du terme cos θ : l'absorption d'énergie ou l'émission d'énergie se fait de façon privilégiée
si le moment dipolaire électrique ~ ba
D est parallèle à ~.
E Voir gure (8.2.6).
D
E θ
Atome
Figure 8.2.6
Nous discuterons au chapitre suivant le terme |Dab |2 qui est souvent nul pour des raisons
de symétries. (règles de sélections).
8.2.3 Eet d'une onde incohérente
Nous supposons toujours l'atome soumis à l'inuence d'une onde électromagnétique à

partir de t = 0. Mais maintenant, nous supposons l'onde incohérente, c'est à dire que
c'est une superposition d'onde planes de diérentes fréquences ω, (avec des déphasages
aléatoires).
C'est le cas de la lumière émise par une ampoule par exemple, où cette lumière est la
superposition de nombreux trains d'ondes provenant de la désexcitation indépendantes de
nombreux atomes.
Au contraire un faisceau cohérent (dans une certaine limite) est produit par un dispositif
LASER.
8.2. THÉORIE DES PERTURBATIONS DÉPENDANT DU TEMPS 327
Du fait de l'indépendance statistique des ondes planes composant un faisceau incohérent,

on somme les probabilités (et non pas les amplitudes).
8.2.3.1 Question :
Dans le cas d'un faisceau incident incohérent, contenant des fréquences proches de
ωab , quelle est la probabilité de transition de l'état atomique |ψa ivers un autre état
|ψb i du spectre discret :
inc
Pa→b (t) =?
8.2.3.2 Réponse
On a Z ∞
inc
Pa→b (t) = Pab (t) dω
0
on utilise (8.2.8), (8.2.6) et (8.2.10), donnant :
2π πI
inc
Pab (t) = 2
|Wba |2 t = 2 cos2 θ |Dab |2 t
~ ~ cε0
où I est l'intensité de l'onde incidente à la fréquence ωab .

qui n'est valable bien sûr tant que P 1 (car on utilise la théorie de perturbation qui
inc
suppose cela). La croissance Pab (t) est donc linéaire. (Alors que Pab (t) était quadratique
dans le cas cohérent).
On déni :
inc
dPab 2π
Tab = = 2 |Wba |2 = constante : taux de transition
dt ~
et
Tab 0
σab = (~ωab ) : section efficace d absorption
I
Si la radiation incidente est non polarisée et isotrope, c'est à dire si les trains d'ondes
ont des polarisations et directions indépendantes, on obtient :
I
Tab = |Dab |2
3~2 cε 0
Démonstration. Il faut faire la moyenne sur les directions. @@

Pinc
a−b (t)
linéaire
t
Figure 8.2.7 Loi linéaire de la transition induite par un faisceau incohérent.
Exercice (TD) : @@ Bransden p508, la loi de somme de Thomas Reiche Kuhn.
8.2.3.3 Loi de décroissance
@@
Bransden p509, distribution Lorentzienne, partie imaginaire de E.
8.2.4 Transition vers les continuum ; La photo-ionisation ; Eet

photoélectrique
On considère toujours le même modèle que précédemment, mais cette fois ci, on s'inté-
resse à la transition de l'état fondamental d'énergie Ea (état lié) vers un état du continuum
d'énergie Eb Ea , qui est un état non lié, où l'électron de l'atome est libéré et part, voir
gure 1.5.2, page 56.
Dans l'état nal, l'atome est donc ionisé.
On notera
dn
ρ(E) = : densité de niveaux du continuum par unité de volume
dE
la densité de niveau qui est le nombre de niveaux par intervalle d'énergie et unité de volume,
que l'on supposera constante autour de Eb .
2
On supposera aussi |Wab | constante autour de Eb .
Voir gure 8.2.8.
E b =E a + h ω Continuum
....
Transition
Ea Fondamental
Figure 8.2.8
inc
Question : En sommant sur les probabilités Pa→b (t) (somme sur les états naux
b), donner Pa→cont (t) qui est la probabilité de ionisation de l'atome ?
Solution
X
inc
Pa→cont = Pa→b (t)
b
Z
inc dn
= Pa→b (t) ρ(E) dE car dn = dE = ρ(E) dE
dE
dE
R inc
Donc, comme
dΩ
= ~, Pa→cont = ~ Pa→b (t) ρ(E)dΩ,
2π
Pa→cont (t) = ρ(Eb ) |Wab |2 t
~
Le taux de transition est
dPa→cont (t) 2π
Ta→cont = = ρ(Eb ) |Wab |2
dt ~
Ce dernier résultat s'appelle la règle d'or de Fermi.
8.2.4.1 Remarques
◦ Noter que cette probabilité est proportionnelle au temps t (tant que Pa→cont 1
bien sûr), et aussi proportionnelle à la densité de niveaux ρ(Eb ).
Pa−cont (t)
linéaire
t
Figure 8.2.9 Loi linéaire de la transition induite par un faisceau cohérent vers un conti-
nuum d'états.
◦ (TD) Si l'énergie cinétique nale est très supérieure à l'énergie potentielle alors
p2
E ' 2m comme pour une particule libre. La densité de niveaux est (voir TD, ou
chapitre 1 @@) :
dn 2π
ρ(E) = = 3 (2m) V p
dE h
où V est le volume où se trouve la particule libérée (V grand et n'apparaîtra pas à
la n).
on obtient alors @@ (Bransden p519), pour la section ecace diérentielle d'io-
nisation (section ecace par unité d'angle solide) :
5
dσ ~ Z
(ω) = 32α cos2 θ
dΩ mω k aa
Pour un faisceau non polarisé, on obtient :
5
dσ ~ Z
(ω) = 16α sin2 Θ
dΩ mω k aa
où Θ est l'angle entre ~k et l'axe de faisceau incident.

◦ La section ecace totale (intégrée sur toutes les directions) est :
√ 5 7/2
16π 2 Z ~
σ(ω) = α
3 a0 mω
◦ courbes expérimentales @@
electron éjecté
E k
θ
onde incidente Atome

Figure 8.2.10 Diusion de l'onde électromagnétique polarisée rectiligne, sur un atome,
éjectant un électron dans le vecteur impulsion ~k .
Faisceau Θ
incident
Polarisé z non
polarisé
(a) (b)
Figure 8.2.11 Distribution angulaire de probabilité de l'électron éjecté lors de l'ionisa-
tion, dans le cas (a) d'un faisceau polarisé selon z
ou (b) d'un faisceau non polarisé.
Exercices
◦ TD : Photo-désintégration du deuton, Bransden p520 @@
◦ TD : Théorie élémentaire de la constante diélectrique ε(ω).
8.3 Méthodes variationnelles

8.3.1 Méthode variationnelle pour problèmes stationnaires
ref : Ballentine [L.E90]p.270.

Le résultat que nous allons montrer est très simple, mais très utile dans de nombreux
problèmes de physique quantique, où il est impossible de connaître la solution exacte.
8.3.1.1 Propriété
Théorème 8.3.1. L'énergie fondamentale de l'opérateur Ĥ est donnée par
hφ|Ĥ|φi
E0 = min (8.3.1)
φ∈H,φ6=0 hφ|φi
c'est à dire E0 est le minimum de l'énergie moyenne obtenue sur tous les états
quantiques possibles.
P
Démonstration. si le spectre de Ĥ est : Ĥ|ψn i = En |ψn i, on a |φi = n |ψn ihψn |φi, et E0 ≤ En
donc : X X
hφ|Ĥ|φi = En |hφ|ψn i|2 ≥ E0 |hφ|ψn i|2 = E0 hφ|φi
n n
hφ|Ĥ|φi
Donc E0 ≤ hφ|φi . Dans le cas de l'état fondamental φ = ψ0 , on a égalité.

Cette propriété ne dit en fait rien d'autre que : la moyenne d'une série de nombre est
plus grande que le plus petit des nombres de cette série.
8.3.1.2 Méthode variationnelle
Dans un problème où l'on cherche à estimer l'état fondamental |ψ0 i et l'énergie E0 d'un
Hamiltonien Ĥ :
1. On choisit un famille d'états quantiques |φµ i appelée famille d'essai, dépendant

d'un paramètre µ, ou de plusieurs paramètres.
8.3. MÉTHODES VARIATIONNELLES 333
hφµ |Ĥ|φµ i
2. On calcule Eµ = hφµ |φµ i
et on cherche la valeur µ∗ qui rende Eµ minimum. Voir
gure 8.3.1.
3. On espère alors que |φµ∗ i sera une bonne approximation de l'état fondamental
|ψ0 i inconnu, et que Eµ∗ sera une bonne approximation de son énergie E0 .
Eµ
Εµ∗
E0
µ∗ µ
Figure 8.3.1 Illustre la méthode variationnelle : Pour estimer E0 , on recherche le mini-

mum de Eµ = hφµ |Ĥ|φµ i/hφµ |φµ i en fonction de µ.
Remarque 8.3.2. Dans de nombreux problèmes importants (ex.. théorie B.C.S. de la su-
praconductivité) cette méthode est irremplaçable et correspond souvent à la théorie du
champ moyen. Cependant la diculté est de choisir la bonne famille d'essai sur des
critères d'arguments physiques (pour bien approcher E0 ), et aussi de calculabilité (il faut
pouvoir calculer la moyenne). Souvent les calculs se font à l'ordinateur.
8.3.1.3 Exemple : vibration anharmonique d'un atome (TD)
On reprend l'exemple (8.1.4) page 314, avec
p̂2 1
Ĥ = + kx̂2 + λx̂4
2m 2
pour λ > 0, il est impossible de calculer analytiquement l'état fondamental.
Pour la fonction d'onde d'essai, il est naturel de choisir une fonction gaussienne nor-
malisée, de largeur variable σ, voir (1.1.6), page 22 :
x2

1
hx|ψσ i = exp − 2
(πσ 2 )1/4 2σ
Calcul de l'énergie moyenne :
+∞
~2
Z 2
2 2 4
1 2 4 1 x
Eσ = hψσ |Ĥ|ψσ i = dx − −1/σ + x /σ + kx + λx 2 exp − 2
−∞ 2m 2 (πσ 2 ) σ
C'est une intégrale Gaussienne (A.1.2), et donne :

r
k w 1 3 1
Eσ = ~ + + g 2
m 4 4w 4 w
avec
~
w= √ >0
σ 2 km
et
~2 λ
g= >0
km
Recherche du minimum de Eσ On observe que pour w → 0 ou w → +∞, alors

Eσ → +∞. Cela garantit l'existence d'un minimum de la fonction E (w), que l'on trouve
en écrivant :
dE
= 0
dw
⇔ w3 − w − 6g = 0
C'est une équation du troisième degré.
Posons
G = 35 g 2
Alors si G > 1, la solution est

1 1
w = √ α+
3 α
√ √
α = G+ G−1
et si 0≤G≤1 la solution est
2 ϕ
w = √ cos
3 3
√
cos ϕ = G
Résultats numériques et observations On pose Evariat. = E (wmin ). Comparons à la

solution exacte (numérique avec ~ = 1, m = 1, k = 1), tracée sur la gure 8.1.3.
λ 0 0.1 1 10 100
Eexact 0.5 0.55915 0.8038 1.5050 3.1314
Evariat. 0.5 0.5603 0.8125 1.53125 3.19244
δE = Evar. − Eexact 0 0.0012 0.0087 0.026 0.06
cf programme plan.cc
On remarque que les résultats sont bons, avec une faible erreur, et d'autant meilleurs
que λ → 0. Cela tient du fait que la fonction d'onde d'essai est la solution exacte pour
λ = 0.
8.3. MÉTHODES VARIATIONNELLES 335
Exercices (suggestions) TD : États électron. de l'hélium par th. perturbation 1er ordre,
et methode variationnelle. Tenir compte du spin, cf Sakurai p369.
Ex : excitons en puits quantiques, cf livre Singh p300.
Approx particules indépend pb dépendant du temps et états coh. donnant la mécanique
classique.
TD : modèle de Friedrichs et applications.
8.3.1.4 (*) Formule du minimax
La formule suivante généralise la formule (8.3.1) qui donnait l'état fondamental.
Théorème 8.3.3. Si Ĥ est un opérateur autoadjoint dont le spectre de valeurs propres est
(compté avec multiplicités)
E0 ≤ E1 ≤ E1 . . . ≤ En ≤ . . .
alors le niveau n est donné par

!
hφ|Ĥφi
En = min max
L⊂H,dimL=n φ∈L hφ|φi
Remarque 8.3.4. Cette formule est utilisée en chimie quantique.
(1) (2)
Exercice 8.3.5. (*) (ref : Arnold [Arn76]chap.5). Si Ĥ1 ≤ Ĥ2 montrer que En ≤ En .
Application : dans Rn , si 2 ellipsoides E 0, E centrés en 0 sont tels que E 0 ⊂ E , montrer que
les 1/2 axes vérient tous
a0j ≤ aj , ∀j = 1 . . . n
Démonstration. @@écrire.
Troisième partie
Mécanique quantique avancée
337
Chapitre 9
Statistiques quantiques et décohérence
Ce chapitre est une introduction aux notions de physique quantique statistique et aux
processus de décohérence. Ces notions sont essentielles pour décrire les systèmes physiques
complexes et en particulier pour aborder la physique statistique. Ce formalisme est aussi à la
base des théories récentes qui décrivent la décohérence, c'est à dire 'apparance classique du
monde à partir des phénomènes quantique (intéraction et couplage avec l'environnement).
9.1 Description d'un ensemble statistique d'états quan-

tiques par un opérateur densité
9.1.1 Dénition de l'opérateur densité
Probabilité et hasard en mécanique classique : A priori, en mécanique classique, un

système peut être en principe totalement connu et si il est régit par des lois déterministes,
on connait son passé et son futur parfaitement. Par exemple l'état d'une particule est
spéciée par sa position et son impulsion (~x, p~).
Cependant, à cause des phénomènes de chaos déterministe (sensibilité aux condi-
tions initiales), cette connaissance parfaite est illusoire et impossible. On montre que le
chaos est comme une source de hasard. Une description probabiliste du système devient
indispensable même pour un système à peu de degrés de libertés (voir cours de licence, mé-
canique analytique). Dans le cas le plus simple on pourra considérer des états de la particule
P
(~xi , p~i ) ayant chacun une probabilité Pi , avec i Pi = 1, ou mieux, utiliser une distribution
de probabilité P (~ x, p~) dxdp sur l'espace de phase pour décrire de façon probabiliste l'état
de la particule. On appelle cela un ensemble statistique d'états classiques.
A fortiori pour un système physique composé d'un grand nombre de particules, la des-
cription probabiliste devient indispensable, d'autant plus que des phénomènes d'ensemble
(dit collectifs) apparaissent et ne s'expliquent que dans la description probabiliste. C'est la
théorie de la physique statistique classique, expliquant les phénomènes de transition de
phase, etc.
339
340 CHAPITRE 9. STATISTIQUES QUANTIQUES ET DÉCOHÉRENCE
Probabilité et hasard en mécanique quantique : En mécanique quantique, il y a

une autre source de hasard, qui est un hasard qui apparait lorsque un système quantique
intéragit avec son environnement. Dans l'état des connaissances actuel ce hasard n'est pas
le fruit de mécanismes sous jacents connus, on le considère comme un hasard spontané.
Rappelons :
Armation 9.1.1. le postulat de la mesure : une particule quantique se décrit par
un vecteur ψ∈H dans un espace quantique. Si Â est un opérateur (autoadjoint) associé
à une mesure (observable), dont le spectre est (aj )j , alors la mesure de la grandeur A sur
l'état ψ donnera le résultat aj avec la probabilité

1
kPj ψk2 hψ|Pj ψi

pψ (aj ) = 2 =
kψk hψ|ψi
où Pj est le projecteur spectral sur l'espace propre de la valeur propre aj . De plus, si la

valeur observée est aj , alors après la mesure, l'état quantique est décrit par Pj ψ (le passage
ψ → Pj ψ s'appelle le collapse ou réduction de l'état quantique).
On peut reformuler ce postulat de la façon suivante :
Armation 9.1.2. Postulat de la mesure : après la mesure, le système quantique

ψ est décrit par la superposition statistique des états Pj ψ chacun ayant la probabilité
kPj ψk2 hψ|Pj ψi
pψ (aj ) = kψk 2 = hψ|ψi .
Insistons sur le fait queces probabilités pψ (aj ) pour toutes les observables pos-
sibles Â caractérisent le système physique au sens où ce sont les seules valeurs
accessibles à l'expérience que l'on peut confronter à la théorie. Par exemple, on peut dé-
duire la valeur moyenne de l'observable A est
2
X hψ|Âψi
hAiψ := pψ (aj ) .aj =
j
hψ|ψi
1.
Démonstration. Rappelons la preuve de la dernière égalité, utilisant que Pj est un projecteur orthogonal,
† 2
donc que Pj = Pj et Pj = Pj :
kPj ψk
2
hPj ψ|Pj ψi hψ|Pj† Pj ψi hψ|Pj2 ψi hψ|Pj ψi
pψ (aj ) = 2 = = = =
kψk hψ|ψi hψ|ψi hψ|ψi hψ|ψi
P P hψ|Pj ψi P hψ|aj Pj ψi hψ|Âψi

2. rappelons la preuve de la dernière égalité : j pψ (aj ) .aj = j hψ|ψi .aj = j hψ|ψi = hψ|ψi
P
car Â = j aj Pj .
9.1. DESCRIPTION D'UN ENSEMBLE STATISTIQUE D'ÉTATS QUANTIQUES PAR UN OPÉRAT
Ce que l'on veut rajouter dans ce chapitre est qu'il peut aussi y avoir du hasard qui se
rajoute, dû à la méconnaissance du système physique. Au lieu d'un unique état quantique
ψ ∈ H, on voudrait donc pouvoir considérer un ensemble statistique d'états
P ψi ∈ H,
i = 1, 2, . . . avec des probabilités respectives pi , vériant i pi = 1. D'après le postulat de
la mesure ci-dessus, après une mesure de l'observable Â on aura donc le résultat aj avec la
probabilité :
X
p (aj ) = pi .pψ (aj )
i
Rappelons encore (car c'est fondamental) que la connaissance des probabilités p (aj ) pour
toutes les observables Â caractérise le système physique. Par exemple on déduit la valeur
moyenne de A par :
X X hψi |Âψi i
hAi = p (aj ) aj = pi
i i
hψi |ψi i
Cette présence de deux sources de probabilité est mal commode et la dénition suivante
va permettre de simplier le formalisme (en gardant l'essentiel).
Dénition 9.1.3. Pour un ensemble statistique d'état ψi ∈ H, i = 1, 2, . . . avec des

probabilités respectives pi , l'opérateur densité associé est
X |ψi ihψi |
ρ̂ := pi Pψi , Pψi := (9.1.1)
i
hψi |ψi i
(aussi appelé matrice densité ou état quantique), où Pψi est le projecteur orthogonal
de rang 1 sur l'état ψi .
Dans le particulier d'un seul état ψ, on a ρ̂ = Pψ , on dit que c'est un état pur.
La propriété suivante montre que l'on peut exprimer les probabilités p (aj ) à partir de
l'opérateur densité ρ̂.
Proposition 9.1.4. L'opérateur densité ρ̂ caractérise uniquement le système physique :

pour toute observable Â , la probabilité d'observer la valeur propre aj est donnée par :
p (aj ) = Tr (ρ̂Pj ) (9.1.2)
où Pj est le projecteur spectral de Â (dénit plus haut). En particulier la valeur moyenne

est
hAi = Tr ρ̂Â
Inversement la connaissance des probabilités p (aj ) pour toute observable Â détermine ρ̂.
Démonstration. On a (utilisant que Tr (|ψihϕ|) = hϕ|ψi)

!
X |ψi ihψi |
Tr (ρ̂Pj ) = Tr pi Pj
i
hψ i |ψi i
X hψi |Pj ψi i X
= pi = pi .pψ (aj ) = p (aj )
i
hψi |ψ i i i
P
Et comme Â = j aj Pj ,
X X
Tr ρ̂Â = aj Tr (ρ̂Pj ) = aj p (aj ) = hAi
j j
Inversement, montrons que on peut reconstruire ρ̂ à partir de la connaissance des probabi-

lités p (aj ) pour toute observable Â : on verra ci-dessous que
ρ̂
admet une décomposition
P
spectrale ρ̂ = i ρi Pφi . Prenons Â = Pφi ; on obtient Tr ρ̂Â = ρi , ce qui permet de
reconstruire ρ̂.
Exemple 9.1.5. Illustrons le fait que l'opérateur densité ρ̂ caractérise le système quantique
sur l'exemple d'un spin 1/2. L'espace quantique est H = C2 . Les ensembles statistiques
1 1
|+z i, p+ = , |−z i, p+ = ,
2 2
et
1 1
|+x i, p+ = , |−x i, p+ = ,
2 2
donnent tous les deux la même matrice densité ρ̂ = 12 Id et sont donc indentiques. Sachant
que
1
|±x i = √ (|+z i ± |−z i)
2
cela peut paraitre paradoxal. En eet en mécanique classique un mélange de pièces noires
et blanches n'est pas équivalent à un mélange de pièces (demi noire/ demi blanche) et
(demi blanche/ demi noire).
Proposition 9.1.6. Un opérateur ρ̂ est un opérateur densité si et seulement si c'est un

opérateur autoadjoint dont les valeurs propres réelles ρj vérient
X
0 ≤ ρj ≤ 1 et ρj = 1.
j
|φj ihφj |
En notant (φj )j une base orthonormée de vecteurs propres de ρ̂, et Pφj = hφj |φj i
, on peut
écrire : X
ρ̂ = ρj Pφj
j
et interpréter ρ̂ comme l'opérateur densité associé à l'ensemble statistique des états φj

avec les probabilités respectives ρj .
Démonstration. Supposons que ρ̂ est un opérateur densité (9.1.1). Il est clair que ρ̂† = ρ̂
est autoadjoint. Ses valeurs propres réelles ρj vérient

X X |ψi ihψi | X
ρj = Tr (ρ̂) = pi Tr = pi = 1
j i
hψi |ψi i i
Il reste à montrer que ρj ≥ 0. Pour cela si v∈H est un vecteur quelconque alors
X hv|ψi ihψi |vi X |hψi |vi|2

hv|ρ̂vi = pi = pi 2 ≥ 0
i
hψi |ψi i i
kψi k
(on dit que c'est un opérateur positif) donc en particulier pour v = φj ( associé à la
valeur propre ρj ) on a
0 ≤ hφj |ρ̂φj i = ρj .
La réciproque est immédiate.
Remarque 9.1.7.
P
◦ La condition j ρj = 1 s'écrit aussi
Trρ̂ = 1.
◦ ρ̂ est un état pur si et seulement si ses valeurs propres sont
ρ1 = 1, ρ2 = ρi≥2 . . . = 0
Exercice 9.1.8. Montrer que ρ̂ est un état pur si et seulement si
Trρ̂2 = 1
9.1.2 Formulation du postulat de la mesure avec la matrice densité
Le postulat de la mesure énoncé en 9.1.1 peut être réformulé de la façon suivante :
Proposition 9.1.9. le postulat de la mesure : Si un système quantique est décrit

par l'opérateur densité ρ̂, si Â est un opérateur autoadjoint (observable) associé à une
mesure alors après la mesure de la grandeur A, l'opérateur densité devient
X
ρ̂0 = Pj ρ̂Pj (9.1.3)
j
où Pj est le projecteur spectral sur l'espace propre de la valeur propre aj de Â. C'est à
dire que ρ̂ devient diagonale en blocs par rapport à la décomposition spectrale de Â.
Démonstration. Avant la mesure, on a une distribution statistique d'états

P ψi avec probabi-
lité pi . La matrice densité correspondante est ρ= ipi Pψi . D'après le postulat de la me-
sure, Armation 9.1.2, après la mesure, chaque état ψi donne une distribution statistique
0 hψ |P ψ i
d'étatsψi,j = Pj ψi avec la probabilité pi,j = pi hψi i |ψj i ii . La matrice densité correspondante
est donc

0
X X hψi |Pj ψi i Pj |ψi ihψi |Pj X
ρ = pi,j Pψi,j
0 = pi = Pj ρ̂Pj
i,j i,j
hψi |ψi i hPj ψi |Pj ψi i j
On a utilisé que hPj ψi |Pj ψi i = hψi |Pj ψi i.
Remarque 9.1.10. L'opérateur (9.1.3) est diagonal par blocs. D'après l'interprétation sta-
tistique de l'opérateur densité, il décrit une superposition statistique d'évènements j, et la
probabilité de l'évènement j est la trace du bloc et donne :
Tr (Pj ρ̂Pj ) = Tr ρ̂Pj2 = Tr (ρ̂Pj ) = p (aj )

d'après (9.1.2). C'est exactement ce que donne le postulat de la mesure.
Dans la section suivante nous évoquerons la théorie de la décohérence qui propose

des mécanismes qui permettent d'expliquer de façon dynamique comment l'opérateur den-
sité peut évoluer d'un état initial quelconque ρ̂ vers l'expression (9.1.3) dans le cadre de
l'interaction avec un système de mesure et avec l'environement.
9.1.3 Opérateur densité pour un système à deux états
Un système à deux états est un système quantique dont l'espace quantique est H = C2 ,
de dimension 2. C'est le cas par exemple du spin 1/2.
On rappelle que ˆ = (σ̂x , σ̂y , σ̂z )

~σ sont les matrices de Pauli dénies en (4.3.5) :

0 1 0 −i 1 0
σx = , σy = , σz =
1 0 i 0 0 −1
Proposition 9.1.11. Dans l'espace quantique H = C2 , un opérateur densité ρ̂ s'écrit

toujours sous la forme
1 ˆ

ρ̂ = Id + P~ .~σ
2
caractérisé par un vecteur P~ ∈ R3 , vériant P~ ≤ 1, appelé vecteur polarisation. De
plus, ρ̂ est ~ = 1. On dit que ρ̂ est l'état non polarisé

un état pur si et seulement si P
si P~ = ~0 soit ρ̂ = 21 Id.
Démonstration.
Toute matrice hermitienne 2 × 2 de trace 1 peut s'écriresous la forme
1 ˆ
Id + P .~σ avec P ∈ R . Les 2 valeurs propres sont alors ρ± = 2 1 ± P~ . La condition
~ ~ 3 1
2
supplémentaire 0 ≤ ρ± ≤ 1 donne P~ ≤ 1.
9.1.4 Equation d'évolution
Proposition 9.1.12. Si Ĥ (t) est l'opérateur Hamiltonien qui fait évoluer chaque état
dψ
de l'ensemble statistique selon i~ dti = Ĥ (t) ψi (t) alors l'opérateur densité associé évolue
selon
dρ̂ −i h i
= Ĥ, ρ̂
dt ~
appelée équation de Liouville quantique.
∂f
Remarque 9.1.13. Cela ressemble à l'équation de transport de Liouville dt = {H, f } en
mécanique analytique [Fau10c].
Démonstration. On a
X |ψi (t)ihψi (t) |
ρ̂ (t) = pi
i
hψi (t) |ψi (t)i
où ψi (t) = Û (t) ψi (0)dÛ

dt
= −i
~
Ĥ · Û (t). On a aussi hψi (t) |ψi (t)i = hψi (0) |ψi
avec (0)i.
Remarquer que les poids de probabilité pi sont constant au cours du temps. Donc
X Û (t) |ψi (0)ihψi (0) |Û (t)+

ρ̂ (t) = pi
i
hψi (0) |ψi (0)i
= Û (t) ρ̂ (0) Û (t)+
donc
dρ̂ dÛ dÛ +

= ρ̂ (0) Û (t)+ + Û (t) ρ̂ (0)
dt dt dt
−i −i h i
= Ĥ Û ρ̂ (0) Û + − Û ρ̂ (0) Û + Ĥ = Ĥ, ρ̂
~ ~
9.1.5 Entropie de l'ensemble statistique
Lire
3 history of entropy sur wikipedia.
On commence par dénir l'entropie de Shannon (ou entropie de Kolmogorov-Sinaï)
pour un système dynamique simple que l'on appelle dynamique de Bernouilli.
9.1.5.1 Entropie de Shannon pour une dynamique de Bernouilli
Considérons un ensemble de n évènements possibles {1, 2 . . . n} ayant chacun la pro-

P
babilité pi , i = 1, . . . n vériant i pi = 1. (C'est une loi de probabilité). Soit N ≥ 1.
Considérons une suite de N réalisations aléatoires indépendantes :
~i := (i1 , i2 , . . . , iN ) (9.1.4)
où chaque valeur ij ∈ {1, 2 . . . n} apparait avec la probabilité pij . Comme les évènements de
la suite ~
i sont indépendants, la probabilité pour que cette suite ~i apparaisse est le produit :

p ~i = pi1 .pi2 . . . . .piN (9.1.5)
Remarquons qu'il y a nN = eN log n suites possibles.

Avant de donner le résultat général de Shannon, présentons deux exemples simples an
de développer l'intuition.
3. (ref : wikipedia history of entropy) L'histoire raconte que en 1949 Claude Shannon qui travaillait
en télécomunication au BELL labs a développé la théorie qui suit, en dénissant une fonction H qu'il
voulait appeler incertitude. Il rencontre le physicien-mathématicien Von Neumann. Ce dernier lui fait
remarquer qu'il devrait appeler sa fonction entropie pour deux raisons : 1) car cela correspond au concept
de l'entropie développée depuis Boltzmann et Gibbs 1872 qui mesure le désordre, et 2) car il deviendrait
célèbre en donnant enn une explication claire de ce qu'est l'entropie.
Exemple 9.1.14. Considérons par exemple le cas très simple
n = 3, p1 = p2 = 0.5, p3 = 0 (9.1.6)
Pour N donné, il y a 3N = eN log 3 suites possibles. Les suites ~

i ne contenant pas le chire
3 , de la forme
(1, 2, 2, 1, . . . 2) (9.1.7)
et ont chacune la probabilité
1 N
p ~i = = e−N log 2
2
Le nombre total de telles suites est
N = 2N = eN log 2
2 N

(ce qui est très peu par rapport au total car 2N /3N = 3
→0 pour N → ∞). Les suites
contenant le chire 3 comme
(1, 2, 3, 2, 1, . . . 2) (9.1.8)
sont chacune de probabilité nulle car p3 = 0.
Exemple 9.1.15. Considérons le cas un peu moins simple avec n=2 et
p1 = 2/3 ' 0.7 > p2 = 1/3 ' 0.3

Nous allons traiter cet exemple de façon non rigoureuse, mais cela va nous faire deviner
N
le résultat présenté dans le théorème de Shannon qui va suivre. Il y a 2 = eN log 2 suites
possibles. La suite qui est individuellement la plus probable est
~i = (1, 1, 1 . . .) (9.1.9)

avec p ~i = pN
1 = e
−N log(1/p1 )
, mais elle est unique (on dit que c'est une suite non
standard). La suite
~i = (2, 2, 2, 2 . . .) (9.1.10)
est au contraire très peu probable mais elle est aussi unique. On dit aussi que c'est une
suite non standard. Au contraire une suite standard est de la forme :
~i = (1, 2, 1, 1, 2, 1, 2, 1, 1, 1, 2, . . .) (9.1.11)
avec un nombre d'apparition de 1 qui est N1 ' p1 N et de 2 qui est N2 ' p2 N . On estime
que sa probabilité d'apparition est donc de l'ordre de

ps = p ~i = pN 1 N2
1 p2 = e
N1 log p1 +N2 log p2
= e−N (−p1 log p1 −p2 log p2 )
= e−N S
avec
S = −p1 log p1 − p2 log p2 > 0
appelée entropie. Par ailleurs, on estime que le nombre de telles suites standard est le
nombre de façon de choisir N1 éléments parmi N, soit :
N!
NN ' CNN1 '
N1 !N2 !
En utilisant la formule de Stirling log N ! ' N log N + . . . cela donne
NN ' elog N !−log N1 !−log N2 ! ' eN log N −(p1 N ) log(p1 N )−(p2 N ) log(p2 N )
' eN (−p1 log p1 −p2 log p2 ) = eN S
Il apparait que NN · ps ' 1 et donc, avec probabilité 1, toutes les suites qui apparaissent
sont standard. Remarquons nalement que S < log 2 donc
NN = eN S 2N = eN log 2
ce qui signie qu'il y a très peu de suites standard.
Le théorème qui suit de Shannon et MacMullen précise ce résultat un cas plus général.
Théorème 9.1.16. (Shannon, MacMillan 1949)(cf Khinchin p.17, Zinmeister p.22).

Pour une loi de probabilité (p1 , p2 , . . . pn ) données, la quantité suivante :
X
S=− pi log pi (9.1.12)
i
appelée l'entropie de Shannon, est telle que pour N 1, parmi toutes les nN = eN log n
suites possibles de la forme (9.1.4), seulement le nombre
NN = eN S
d'entre elles peuvent eectivement apparaitre et la majorité d'entre elles ont la probabilité
ps = e−N S
Elles sont appelées suites standard (comme (9.1.11)).

Parmi les autres suites (au total improbables) appelées suites non standard, il y a des
suites en grand nombre mais de très faible probabilité (comme (9.1.10)), ou des suites plus
probables mais en très petit nombre comme (9.1.9).
L'entropie vérie :
0 ≤ S ≤ log n (9.1.13)
Le maximum d'entropie est pour la loi uniforme où tous les évènements sont équi-
pi = 1/n donnant S = log n. Au contraire, le minimum d'entropie est
probables :
S=0 obtenu pour un évènement totalement prévisible comme p 1 = 1 , p2 = . . . = pn = 0 .
(remarquer que pi log pi → 0).
pi →0
Remarque 9.1.17.
◦ On peut aussi retenir la formule :
1
S= log NN
N
qui montre que (par dénition) l'entropie mesure le taux de croissance exponentiel
du nombre NN des suites standard de longueur N .
1 1
◦ Dans l'exemple (9.1.6), l'entropie est S = − log
2 2
− 21 log 12 = log 2 conformément
à (??).
◦ (*) Un énoncé plus précis (mathématiquement) est que pour N → ∞, on peut diviser
N
l'ensemble des n suites en deux ensembles : A (suites standard) et B (suites non
standard) tels que (rappel o (1) désigne un terme qui tend vers 0 pour N → ∞).

∀~i ∈ A, e−N S(1+o(1)) ≤ p ~i ≤ e−N S(1−o(1)) (9.1.14)
eN S(1−o(1)) ≤ ]A ≤ eN S(1+o(1)) (9.1.15)

X X
p ~i = 1 − o (1) , p ~i = o (1) (9.1.16)
~i∈A ~i∈B
le facteur o (1)
qui apparait provient de la loi des grands nombres. En utilisant
1
le théorème central limite qui est plus précis, on peut remplacer o (1) par
N 1/2−ε
pour tout ε > 0.
Démonstration. (*) (Zinmeister p.22). L'astuce est de faire apparaitre une somme de Bir-
kho de variables aléatoires indépendantes. D'après (9.1.5) on écrit
P
p ~i = pi1 .pi2 . . . . .piN = e j log pij
P
= e− j Xj
en posant
Xj := − log pij
Les variables aléatories Xj sont indépendantes, de même loi, et de moyenne qui est préci-
sément l'entropie :
n
X
hXj i = pi (− log pi ) = S
i=1
Considérons la somme moins la moyenne :
N
1 X
S := (Xj − S) (9.1.17)
N j=1
ainsi on peut exprimer :

p ~i = e−N S(1+S)
La loi (faible) des grands nombre (cf Proposition 9.1.18) dit que S converge vers
0 au sens ou pour tout α>0 (arbitrairement petit) et pour N → ∞,
proba (|S| < α) = 1 − o (1) (9.1.18)
Donc on dénit les suites standard par

n o
Aα : = ~i tq −α<S <α
n o
= ~i tq p ~i ∈ e−N S(1+α) , e−N S(1−α)
Cela donne (9.1.14). On a d'après la loi des grands nombres

X
p ~i = 1 − o (1)
~i∈Aα
donnant (9.1.16). Par conséquence
(]A) .e−N S(1+α) ≤ 1 − o (1) et (]A) .e−N S(1−α) ≥ 1 + o (1)

donnant (9.1.15).
Il reste à justier la loi des grands nombres utilisée dans (9.1.18). (réf : cf wikipedia).
Proposition 9.1.18. La loi faible des grands nombres. Supposons que (Xj )j∈N
sont des variables aléatoires indépendantes, de même moyenne hXi et même variance
σ 2 = (X − hXi)2 . Pour N ≥ 1, on pose
N
1 X
SN := (Xj − hXi)
N j=1
qui est la moyenne temporelle centrée. Alors pour tout α>0 (arbitrairement petit) et
pour N → ∞,
proba (|SN | > α) = o (1)
Démonstration. Notons pour simplier Z = |S|2 et p (Z) dZ la densité de probabilité. On

a
proba (|S| > α) = proba |S|2 > α2

Z ∞ Z ∞
Z
= p (Z) dZ ≤ p (Z) dZ
α2 α2 α2
Z ∞ Z ∞
1 1
= 2 Zp (Z) dZ ≤ 2 Zp (Z) dZ
α α2 α 0
1 1
= 2 hZi = 2 |S|2
α α
où h.i désigne la moyenne. (Ce dernier calcul s'appelle inégalité de Tchebychev). On

a h(Xj − hXi)i = hXi − hXi = 0 et donc pour j 6= k h(Xj − hXi) (Xk − hXi)i =
h(Xj − hXi)i h(Xk − hXi)i = 0. Alors d'après (9.1.17)
1 X 1 X
|S|2 = 2
h(X j − hXi) (X k − hXi)i = 2
(Xj − hXi)2
N j,k N j
1 2
= σ →0
N
On déduit que proba (|S| > α) = o (1).
@@Tracer la loi de proba de X en faisant apparaitre les évènements extremes et la conver-

gence vers la moyenne@@
Exercice 9.1.19. (voir TD) Considérons les évènements 4 possibles i ∈ N et considérons la

loi de probabilité de Boltzmann
1 −βi
pi = e
Z
1
P β = kT , T
où est la température et Z est un facteur de normalisation de sorte que
i pi = 1.
1. Calculer Z en fonction de β (aussi appelé fonction de partition).
2. Calculer l'entropie S en fonction de β.

3. Donner l'expression de S (les premiers termes du developpement limité) pour β→0
(haute température)
4. Donner l'expression de S (les premiers termes du developpement limité) pour β→

∞ (basse température)
9.1.5.2 Entropie en mécanique quantique
Revenons à la mécanique quantique et à un ensemble statistique d'états quantiques

décrits par l'opérateur densité ρ̂ .
Pour caractériser l'incertitude sur le résultat de mesure d'une observable Â sur l'état
ρ̂, il faut considérer la loi de probabilité p (aj ) = Tr (Pj ρ̂) donnée en (9.1.2).
Dénition 9.1.20. Conformément à la dénition (9.1.12), l'entropie d'un ensemble

statistique d'états quantique décrits par l'opérateur densité ρ̂ , relatif à la mesure
d'une observable Â est :
X X
S ρ̂, Â = − p (aj ) log p (aj ) = − Tr (Pj ρ̂) log Tr (Pj ρ̂)
j j

Le minimum de S ρ̂, Â sur toutes les observables Â possibles est appelé entropie de
l'état ρ̂ :
S (ρ̂) := minÂ S ρ̂, Â
4. Ce modèle correspond par exemple aux niveaux d'énergie de l'oscillateur harmonique. En eet l'os-
Ei = ~ω i + 12 avec i ∈ N, et d'après la loi de Botltzman,

cillateur harmonique a les niveaux d'énergie la
probabilité pour que la particule soit au niveau i d'énergie Ei s'écrit comme dans cet exercice :
1 − Ei 1 ~ω(i+1/2) 1
0
e kT = 0 e− kT = e−βi
Z Z Z
~ω 1 ~ω
1 − 2kT
avec β= kT et Z = Z0 e .
Proposition 9.1.21. L'entropie de l'état ρ̂ est donnée par

X
S (ρ̂) := − ρi log ρi = −Tr (ρ̂ log ρ̂) (9.1.19)
i
où ρi sont les valeurs propres de ρ̂. Il est obtenu pour Â = ρ̂, c'est à dire que S (ρ̂, ρ̂) =
S (ρ̂). En particulier ρ̂ est un état pur si et seulement si S (ρ̂) = 0.
Remarque 9.1.22.
1. Pour expliquer l'écriture dans (9.1.19) : si Â est un opérateur auto-adjoint et si

f : R → R est une fonction analytique, i.e. qui admet un developpement de Taylor
f (x) = f (0) + xf 0 (0) + x2 f ”(0)
2
+ . . . convergent, alors on utilise ce développement
pour dénir
f ” (0)
f Â = f (0) + Âf 0 (0) + Â2 + ...
2
(qui est une série convergente si Â est un opérateur borné). De plus, Â se diagonalise
Â = Û DÛ −1
où Û est un opérateur unitaire et D = Diag (a1 , a2 , . . .) est une matrice diagonale

avec les valeurs propres aj sur la diagonale. Remarquons que pour tout n≥1

Ân = Û DÛ −1 Û DÛ −1 . . . Û DÛ −1 = Û Dn Û −1
donc on a :
f Â = Û f (D) Û −1
où f (D) = Diag (f (a1 ) , f (a2 ) , . . .) est une matrice diagonale. En particulier on a

Tr f Â = Tr Û f (D) Û −1 = Tr Û −1 Û f (D) = Tr (f (D))
X
= f (aj )
j
Le cas f (x) = −x log x et Â = ρ̂ donne eq.(9.1.19).
2. En mécanique classique (ou en théorie de l'information), une entropie nulle signie

par dénition que les évènements qui vont survenir sont totalement prévisible (il n'y
a qu'une possibilité). En mécanique quantique, S (ρ̂) = 0 signie que ρ̂ est un état
pur ψ, et non pas une superposition statistique d'états. Cela n'empèche pas que si
l'on fait des mesures d'observables Â, les résultats sont non prévisibles en général
car la distribution de probabilité dépend de ψ et de l'observable Â. Cela correspond

au fait que l'entropie S ρ̂, Â est non nulle en général.
Par exemple si un spin 1/2 est dans l'état pur √1 (|+z i + |−z i) (ρ̂ =
|+x i =
2
|+x ih+x |) on a S (ρ̂) = 0 selon (9.1.19). Si on mesure le spin selon z avec l'observable

1
Ŝz alors on a les probabilités p± = 2
, et la dénition ci-dessus donneS ρ̂, Ŝz =
log 2. Si on mesure le spin selon x avec l'observable
Ŝ
x alors on a les probabilités
p+ = 1, p− = 0, et la dénition ci-dessus donneS ρ̂, Ŝx = 0 = S (ρ̂)
9.2 Opérateur densité partielle pour un système com-

posé
9.2.1 Rappels sur les systèmes composés
Il est assez courant en physique qu'un système soit composé de sous systèmes plus
élémentaires. On a vu que si un système est composé de deux sous systèmes A et B, et si
l'espace quantique du système A est HA (respect. HB pour B) alors l'espace quantique du
système total est l'espace produit tensoriel :
Htot = HA ⊗ HB
Remarque : le terme sous système est ici au sens large. Ce peut être des dégrés de liberté
interne comme le spin.
Rappelons que si (|ei i)i forme une base o.n. de HA et |fj i une base o.n. de HB alors
|ei i ⊗ |fj i = |ei , fj i

forme une base o.n. de Htot = HA ⊗ HB et donc dans le cas d'espace de dimension nie,
dim Htot = dim HA . dim HB
Exemples de systèmes composés :

◦ l'atome d'hydrogène est formé de proton et électron. Alors
Hatome H = Hproton ⊗ Helectron

◦ Pour électron (de spin 1/2) :
Helectron = L2 R3 ⊗ C2

◦ Pour deux particules de spin 1/2 (on ne décrit que les spins)
H = C2(1) ⊗ C2(2)
L'état quantique du système total peut être :
|ψ1 i = |+z i ⊗ |+z i

9.2. OPÉRATEUR DENSITÉ PARTIELLE POUR UN SYSTÈME COMPOSÉ 355
(i.e. les deux particules ont le spin up selon z) ou encore
1
|ψ2 i = √ (|+z i ⊗ |+z i + |−z i ⊗ |−z i)
2
◦ Un exemple qui va nous intéresser particulièrement est : A : système étudié (par

exemple une molécule) et B : le reste de l'univers, appelé environnement. Natu-
rellement dans cet exemple extrème on ne connait pas le système B (l'univers...),
mais on peut connaitre comment il intéragit avec A (ce seront les particules de gaz
ou molécules voisines de A, ou photons de lumière incidents,...) :
Htot = HA ⊗ Henvironnement
Notons IdA l'opérateur identité dans l'espace HA . Soit |fj i un vecteur de base xé de HB .
Les opérateurs suivant serviront dans la dénition qui va suivre :
IdA ⊗ hfj | : Htot → HA
IdA ⊗ |fj i : HA → Htot
Dénition 9.2.1. Si T̂ : Htot → Htot est un opérateur linéaire dans Htot = HA ⊗ HB

alors on dénit sa trace partielle par rapport à B comme étant l'opérateur

TrB T̂ : HA → HA
dénit par X
TrB T̂ := (IdA ⊗ hfj |) ◦ T̂ ◦ (IdA ⊗ |fj i) (9.2.1)
j
aussi noté plus simplement X

TrB T̂ = hfj |T̂ |fj i
j
Cet opérateur aussi appelé opérateur partiel, ne dépend pas du choix de la base o.n.
(fj )j .
P
Remarque : de même on dénit TrA T̂ := i hei |T̂ |ei i : HB → HB .
Nous allons voir maintenant l'intérêt de ces opérateurs partiels.
Proposition 9.2.2. Si ρ̂ est l'opérateur densité décrivant un système dans l'espace Htot =
HA ⊗ H B , alors l'opérateur densité partiel du sous système A est dénit par
ρ̂A := TrB (ρ̂) .
ρ̂A A au sens où pour toute observable de la forme Âtot = Â⊗IdB (c'est à

décrit le système
dire n'agissant pas sur la partie B ), les probabilités des résultats de la mesure sont donnés,
d'après (9.1.2), par
p (aj ) = Tr (ρ̂ (Pj ⊗ IdB )) = TrA (ρ̂Â Pj )
où Pj : HA → HA , j = 1, 2, . . . sont les projecteurs spectraux pour l'observable Â, Tr est la
trace dans Htot et TrA la trace partielle dans HA . Autrement dit, l'opérateur densité partiel
caractérise complètement le sous système A pour ce qui concerne des mesures portant sur
lui.
Démonstration. D'après (9.1.2) on a
p (aj ) = Tr (ρ̂ (Pj ⊗ IdB )) = TrA (TrB (ρ̂ (Pj ⊗ IdB )))
On note |fl il HB . On remarque que (utilisant la relation de fermeture)
une base de
X X
TrB (ρ̂ (Pj ⊗ IdB )) = hfl |ρ̂ (Pj ⊗ IdB ) |fl i = hfl |ρ̂|fl0 ihfl0 | (Pj ⊗ IdB ) |fl i
l l,l0
X X
= hfl |ρ̂|fl0 i (Pj ⊗ hfl0 |fl i) = hfl |ρ̂|fl iPj
l,l0 l
= TrB (ρ̂) Pj = ρ̂A Pj

Donc p (aj ) = TrA (ρ̂Â Pj ).
Example 9.2.3. Considérons deux spins 1/2.
Htot = C2(1) ⊗ C2(2)

1. Pour l'état
|ψ1 i = |+z i(1) ⊗ |+z i(2) = |+, +i (9.2.2)
qui est un état pur dans Htot , on a
ρ̂ = |ψ1 ihψ1 | = |+, +ih+, +|

et (on détaille la dénition (9.2.1))
ρ̂A = TrB (ρ̂) = (I ⊗ h+|) |+, +ih+, +| (|I ⊗ |+i) + (I ⊗ h−|) |+, +ih+, +| (|I ⊗ |−i)
(9.2.3)
= |+ih+|
Cet état partiel est un état pur. Sa matrice densité (dans la base |+i, |−i) est

1 0
ρ̂A ≡
0 0
Son entropie est S (ρ̂A ) = 0.

2. Pour l'état
1 1
|ψ2 i = √ (|+z i ⊗ |+z i + |−z i ⊗ |−z i) = √ (| + +i + | − −i) (9.2.4)
2 2
qui est un état pur dans Htot , on a
1
ρ̂ = |ψ2 ihψ2 | = (| + +i + | − −i) (h+ + | + h− − |)
2
1
= (| + +ih+ + | + | + +ih− − | + | − −ih+ + | + | − −ih− − |)
2
et
1
ρ̂A = TrB (ρ̂) = (|+ih+| + |−ih−|) (9.2.5)
2
C'est état partiel n'est pas un état pur. Sa matrice densité (dans la base |+i, |−i)
est diagonale
1

0
ρ̂A ≡ 2
1
0 2
C'est un état mélangé (non pur). Son entropie est
1 1 1 1
S (ρ̂A ) = − log − log = log 2
2 2 2 2
C'est l'entropie maximale pour un système à deux états.
L'exemple précédent nous amène à la proposition et dénition suivante :
Proposition 9.2.4. Considérons un état pur d'un système composé ψ ∈ Htot = HA ⊗HB .
◦ Si l'état partiel ρ̂A = TrB ρ̂ est non pur, alors ρ̂B = TrA ρ̂ est aussi non pur, et on
dit que ψ est un état enchevetré ou corrélé.
◦ Au contraire, si ρ̂A = TrB ρ̂ = |ψA ihψA | est un état pur, alors ρ̂B = |ψB ihψB | est
aussi un état pur, et ψ = |ψA i ⊗ |ψB i est un état produit (non enchevétré).
Exemple : l'état partiel (9.2.3) est pur donc ψ1 dans (9.2.2) n'est pas enchevétré. Par
contre (9.2.5) n'est pas pur, donc (9.2.4) est un état enchevétré ou corrélé.
Remarque : en anglais, enchevétré = entangled.
9.2.2 Décomposition de Schmidt
Pour mettre en valeur la proposition suivante, remarquons que si |ei ii est une base de
HA et |fj ij est une base de HB , donc |ei , fj ii,j est une base de Htot = HA ⊗ HB et un état
ψ ∈ Htot se décompose dans cette base par une double somme :
X X
|ψi = |ei , fj i hei , fj |ψi = ψi,j |ei , fj i
| {z }
i,j i,j
ψi,j
Proposition 9.2.5. Décomposition de Schmidt d'un état pur du système total.

Siψ ∈ Htot = HA ⊗ HB alors il existe une base |ei ii=1→I de HA et une base |fj ij=1→J de
HB (qui dépendent de ψ ) telles que
min(I,J)
X
|ψi = ψi |ei , fj i (9.2.6)
i
2
P
avec des composantes positives ψ1 ≥ ψ2 ≥ . . . ≥ 0, i ψi = 1. Par dénition, le
vériant
rang de Schmidt de ψ est le nombre de composantes ψi non nulles.
∗
Démonstration. En utilisant l'identication (métrique) HA ≡ HA le vecteur ψ ∈ HA ⊗
∗
HB s'identie à ψ ∈ HA ⊗ HB et donc à un opérateur linéaire ψ ∈ L (HA → HB ). La
décomposition de Schmidt (9.2.6) est une application directe de la décomposition en valeur
singulières de cet opérateur ψ ∈ L (HA → HB ). Pour les détails, voir [Fau10b].
Conséquence pour l'opérateur densité ρ̂ = |ψihψ| de cet état pur :
Corollaire 9.2.6. Les opérateur densité partiels sont diagonaux dans la base de Schmidt :
X
ρ̂A = TrB (ρ̂) = ψi2 |ei ihei |
i
X
ρ̂B = TrA (ρ̂) = ψi2 |fi ihfi |
i
et donc pi = ψi2 sont des distributions de probabilité, et l'entropie de ces états partiels est
X X
S (ρ̂A ) = S (ρ̂B ) = − pi log pi = − ψi2 log ψi2 .
i i
En particulier, ρA (et ρB ) sont des états purs et ψ = |e1 i ⊗ |f1 i est un état produit si et
seulement si S (ρ̂A ) = S (ρ̂B ) = 0. Sinon ψ est un état enchevétré.
9.2.3 Un modèle simple de décohérence
Voir le TD Modèle simple de décohérence induit par l'environement. Ce mo-

dèle est issu d'un article de Zureck (90').
Annexe A
Formules
A.1 Analyse et intégrales

A.1.1 Intégrales Gaussiennes
Formule 1 Z +∞ √
exp −X 2 dX = π

(A.1.1)
−∞
R +∞ 2
Démonstration. Soit I= e−x dx. Alors
−∞
Z +∞ Z +∞ Z +∞ Z +∞
−x2 −y 2 −(x2 +y 2 )
I2 = e dx e dy = e dxdy
−∞ −∞ −∞ −∞
Z ∞ Z 2π Z ∞ ∞
−r2 −r2 1 −r2
= drr dθe = 2π re dr = 2π − e =π
0 0 0 2 0
√
donc I= π.
Formule 2 Soit
Q(x) = Ax2 + Bx + C
A, B, C ∈ C et < (A) > 0 pour la convergence.
Alors
+∞
r
B2
Z
π
exp (−Q(x)) dx = exp −C + (A.1.2)
−∞ A 4A
361
362 ANNEXE A. FORMULES
Preuve : @@
Formule 3. Pour n ≥ 0,
Z +∞ √ (2n − 1)!!
2
x2n e−x dx = π
−∞ 2n
où
(2n)!
(2n − 1)!! := (2n − 1) (2n − 3) . . . 5.3.1 ==
2n (n!)
R +∞ 2
appelée double factorielle. (par symétrie, pour une puissance impaire
−∞ x2n+1 e−αx dx =
0).
R +∞ 2
Démonstration. Il faut dériver la formule −∞ e−αx dx par rapport à α et faire α=1 à la n :
soit
+∞
Z r
−αx2 π
Iα := e dx =
−∞ α
On a
+∞
dn Iα √

−1 −3 2n − 1
Z
2 n −αx2

= −x e dx = π ... −
dαn −∞ 2 2 2
donc Z +∞ √ (2n − 1)!!
2
x2n e−αx dx = π
−∞ 2n
A.1.2 Transformée de Fourier
Dans R3 , on a la fonction de Green de la particule libre :
1 eik|x−y|
hx| |yi = (A.1.3)
−∆ + k 2 4π |x − y|
Démonstration. On pose k = iκ ∈ iR. Utilisant la relation de fermeture en base d'impulsion
(et ~ = 1),
Z
1 1
hx| |yi = d3 p~hx|pihp|yi 2
−∆ + k 2 (p + κ2 )
Z
1 1
= 3 d3 p~ei~p.(x−y) 2
(2π) (p + κ2 )
Pour calculer cette intégrale, on passe en coordonnées sphériques et utilise la formule des
résidus pour obtenir
1 e−κ|x−y| eik|x−y|
hx| |yi = =
−∆ + k 2 4π |x − y| 4π |x − y|
A.2. ALGÈBRE 363
A.2 Algèbre
A.2.1 Séries
Série arithmétique et géométrique @@
A.2.2 Diagonalisation d'une matrice 2×2

Formule générale : On considère la matrice

a b
A=
c d
Soit
∆ = (a − d)2 + 4bc.
Si ∆ 6= 0 alors la matrice A se diagonalise
1 :
√ √ √
a+d+ ∆
!
a−d+ ∆ a−d− ∆ 0√
A = P DP −1 , P = , D= 2
a+d− ∆
,
2c 2c 0 2
Remarque : la matrice P contient les vecteurs propres en colonne, et ils sont dénis à
la multiplication près par un nombre.
Cas particulier : Si

0 b |b| − |b| |b| 0
A= , P = , D=
b 0 b b 0 − |b|
A.2.3 Relations de commutation
On a :
[q̂, p̂] = i~
dF (p̂)
[q̂, F (p̂)] = i~
dp̂
dF (p̂)
[F (q̂), p̂] = i~
dq̂
preuve : Voir [CBF]p 172.
1. Avec le logiciel gratuit xcas de calcul formel (pour l'obtenir, taper xcas dans google). Et dans
xcas, écrire : A:=[[a,b],[c,d]]; D:=egvl(A); P:=egv(A); On vérira que simplify(P*D*inv(P));
redonne bien la matrice A .
Si les opérateurs A, B vérient :
[A, [A, B]] = [B, [A, B]] = 0
alors on a la formule de Glauber
exp (A) exp (B) = exp ([A, B] /2) exp (A + B) (A.2.1)
preuve : Voir [CBF]p 174, ou TD.

Autre preuve : poser C = i [A, B] ; Alors les trois opérateurs (A, B, C) vérient l'algèbre de
Weyl-Heisenberg (comme q̂, p̂, i~) :
C = [A, B] , [A, C] = [B, C] = 0
Considérons les matrices :

     
0 1 0 0 0 0 0 0 1
a =  0 0 0 , b =  0 0 1 , c= 0 0 0 
0 0 0 0 0 0 0 0 0
qui vérient aussi la même algèbre. On procède comme dans la section 2.2.6.4, on vérie sur les
matrices 3×3 que ea eb = ec/2 ea+b .
A.2.4 Algèbre des matrices de Pauli
On a
~ ~σ .B
~σ .A ~ = A.
~B ~ + i~σ . A~∧B
~ (A.2.2)
valable même si ~ B
A, ~ sont des opérateurs.
Si ~u vecteur unitaire :
exp (iϕ~σ .~u) = cos ϕ + i sin ϕ ~σ .~u
A.2.5 Relations sur les matrices
Pour une matrice A (ou un opérateur à Trace) on a :
eA = eTr(A)

det
A.2.6 Inverse d'une matrice 2×2

a b
Si A= et Det (A) = ad − bc 6= 0 alors
c d

−1 1 d −b
A =
Det (A) −c a
A.3. CALCUL DIFFÉRENTIEL DANS R3 365
A.2.7 Relations de commutations
Pour 3 opérateurs A, B, C on a
[AB, C] = A [B, C] + [A, C] B
Démonstration. On écrit :
[AB, C] = ABC − CAB = A ([B, C] + CB) − CAB

= A [B, C] + [A, C] B
A.3 Calcul diérentiel dans R3

références : Jackson [Jac75], début et n.
Ici f (~x) représente une fonction. V~ (~x) représente un champ de vecteurs.
A.3.1 Rappels sur le calcul diérentiel vectoriel.
Voir Feynmann Electromagnétisme, chap 2.
◦ Pour un champ scalaire V (x, y, z), le champ vectoriel W ~ =− −→

gradV a pour coordon-
∂V −−→ ~ .
nées : Wx =
∂x
, Wy = ∂V∂y
, Wz = ∂V∂z
. On écrit aussi : gradV = ∇V
~ = (Ex , Ey , Ez ) sont les coordonnées du champ vectoriel E(~
◦ Si E ~ x), alors le champ
−→
~ = rotE ~ a pour coordonnées :
W
Wx = ∂Ez
∂y
− ∂E
∂z
y
∂Ex
Wy = ∂z
− ∂E
∂x
z
∂Ey ∂Ex
Wz = ∂x
− ∂y
~ −
→ ~ ∧E
~ ~ =
on
s'en rappelle
en écrivant : rotE = ∇ avec l'opérateur vectoriel ∇
∂
, ∂, ∂
∂x ∂y ∂z
.
◦ Et
~ = ∂Ex ∂Ey ∂Ez
divE + +
∂x ∂y ∂z
~
qui se retient en écrivant divE ~ E
= ∇. ~.
◦ Le Laplacien ∆V est le champ scalaire :
∂ 2V ∂ 2V ∂ 2V
∆V = + +
∂x2 ∂y 2 ∂z 2
Propriétés importantes :
−→ −−→ ~
rot gradV =0
−→~
div rotE =0
−−→
div(gradV ) = ∆V
−→− → −−→
a) = grad (div(~a)) − ∆~a
rotrot (~
Formules de Stokes
◦ Si γ est un chemin dont a, b sont les extrémités, alors
−−→
Z
grad (f ) .d~
l = f (b) − f (a)
γ
◦ Si S est une surface dontγ est le contour, alors

Z Z
−→ ~
rot V .~n d2 s = V~ d~l
S Γ
◦ Si V est un volume dont

Z le bord est la surface
Z S = ∂V , alors
~
div(V ) d3 v = V~ .~n d2 s
V Σ=∂V
n Σ
b
S ds2
a γ
Lemme de Poincaré : Réciproquement les formules suivantes
~ f
grad = 0 ⇒ f = cste
~ f
~ ~u = 0 ⇒ ∃f, tq ~u = grad
rot
~ ~u
v = 0 ⇒ ∃~u, tq ~v = rot
div ~
ne sont pas toujours vraies. Elles sont vraies sur un domaine D de l'espace qui est une
boule (contractible en un point), ne contenant donc pas de trous.
A.3.2 En coordonnées sphériques :

−−→ 1 1
grad(f ) = (∂r f ) ~
er + ∂θ f ~eθ + ∂ϕ f ~eϕ (A.3.1)
r r sin θ
A.3. CALCUL DIFFÉRENTIEL DANS R3 367
A.3.3 Relations

~a ∧ b ∧ ~c = (~a.~c)b − ~a.~b ~c
~ ~
−→ → − → −
a ∧ ~b = ~a div ~b − ~b (div ~a) + ~b. ∇ ~a − ~a. ∇ ~b
rot ~
−→ − → −→
a ∧ ~b = rot(~a).~b − ~a.rot(~b)
div ~ (A.3.2)
−→ ~ −−→ −→
rot f.V = grad(f ) ∧ V~ + f rot(V~ ) (A.3.3)
div (f ~a) = grad (f ) .~a + f.div(~a) (A.3.4)
2
div(~
er ) =
r
→ − 1 ~a⊥
~a. ∇ ~er = (~a − (~a.~er ) ~er ) ≡
r r
Annexe B
Solutions des exercices
B.1 Chapitre 1
Exercice 1.1.2 page 24
Z Z
1 2 2 1
dX exp −X 2 = 1

< x0 , p0 , σ|x0 , p0 , σ >= dx exp − (x − x0 ) /σ =
(πσ 2 )1/2 (π)1/2
x−x0
C'est une intégrale Gaussienne, voir eq.(A.1.1), avec le changement de variable X= σ
.
Exercice 2 page 43 Il y a une innité d'autres possibilités plus ou moins explicites. En

voici un exemple simple. (Wk )k déni par
W0 = V0
1 1
Wk = √ (Vk + V−k ) , W−k = √ (Vk − V−k ) , si k>0
2 i 2
q q
2 k2πx 2 k2πx

On remarque que Wk (x) = L
cos L
, W−k (x) = L
sin L
pour k > 0.
Exercice 4 page 51 < x|x̂|φ >= x < x|φ >. On a en eet

Vérions la relation
0 0 0 0
R R 0 0 0
hx| x |x >< x |dx |φi = x δ(x − x)hx |φi = x hx|φi.
Remarque : une telle écriture n'est pas correcte mathématiquement car |xi ∈
/ H. L'outil
mathématique approprié est le théorème spectrale. L'idée est d'introduire le projecteur
Rx R +∞
P̂x = ” −∞ |xihx|dx” et d'écrire de façon équivalente x̂ = −∞ x dPx .
369
370 ANNEXE B. SOLUTIONS DES EXERCICES
Exercice 5 page 53 on a
Z
ψ̃x0 ,p0 ,σ (p) =< p|x0 , p0 , σ >=
dx < p|x >< x|x0 , p0 , σ >
!
(x − x0 )2

(p0 − p) x
Z
1 1
=√ dx exp i exp −
2π~ (πσ 2 )1/4 ~ 2σ 2
!
2
−
r
σ xp x p (p p )
= √ exp i 0 0 exp −i 0 exp − 0
(B.1.1)
~ π ~ ~ 2(~/σ)2
où l'intégrale Gaussienne se calcule par la formule (A.1.2). Remarquer la similarité de

l'expression obtenue pour ψ̃x0 ,p0 ,σ (p) avec celle de ψx0 ,p0 ,σ (x), mis à part facteur de phase
constant.

1. Partant de hAi (t) = hψ(t)|Â(t)|ψ(t)i = hψ|Û (−t)Â(t)Û (t)|ψi avec Û (t) = exp −iĤt/~ .
dhAi(t)
On dérive en t=0 :
dt
= hψ| iĤ
~
Â − Â iĤ
~
+
dÂ
dt
|ψi donnant le résultat souhaité.
dhHi
2. Pour Â = Ĥ supposé indépendant du temps, cela donne
dt
= 0 qui montre la
conservation de l'énergie moyenne au cours du temps.
Dh iE
1 1 1
3. On a d hx̂i /dt = i~
x̂, Ĥ = i~2m
h[x̂, p̂2 ]i = m
hp̂i où la dernière égalité utilise
Dh iE
une relation de la page 363. De même d hp̂i /dt = i~1 p̂, Ĥ = 1
i~
h[p̂, V (x̂)]i =
dV
− dx
(x̂) .
B.2 Chapitre 2
Exercice 2.1.5 page 94 par dénition, < x|T̂λ |ψ >=< x − λ|ψ >, ∀ |ψ >. Donc <
x|T̂λ =< x − λ|. Donc T̂λ+ |x >= |x − λ >. Ensuite, T̂λ est un opérateur unitaire, et
T̂λ+ = T̂λ−1 = T̂−λ . donnant T̂−λ |x >= |x − λ >.

−i

1. Û (t0 ) = exp ~
Ĥt0 . En eet si on dérive |ψ (t + t0 )i = Û (t0 ) |ψ (t)i, par rap-
∂|ψ(t+t0 )i −i
Ĥ Û (t0 ) |ψ (t)i = −i

port à t0 , on obtient
∂t0
= ~ ~
Ĥ|ψ (t + t0 )i qui est
l'équation de Schrödinger. Par ailleurs pour t0 = 0, on a bien Û (t0 = 0) = Iˆ.
2. On a
∂hx|ψx0 i ∂ ∂ ∂ ∂
= ψ (x − x0 ) = − ψ (x − x0 ) = − ψx0 (x) = − hx|ψx0 i
∂x0 ∂x ∂x ∂x ∂x
0
−i ∂ −i
= −i~ hx|ψx0 i = hx|p̂ψx0 i
~ ∂x ~
B.2. CHAPITRE 2 371
∂ ∂|ψx0 i
−i

car p̂ ≡ −i~ ∂x . Donc
∂x0~
p̂|ψx0 i. Notant |ψx0 i =
= T̂x0 |ψi, comme dans la
∂ T̂ (x0 )
= −i −i

question (1), on déduit que
∂x0 ~
p̂T̂ (x0 ) et donc T̂ (x0 ) = exp ~
p̂x0 .
3. Soit ψ̃ (p) = hp|ψi. En suivant le même calcul que ci-dessus, le générateur des
∂ ∂
boosts est b̂ = −i~ . Or on sait que x̂ = i~ (cela se vérie sur l'équation
∂p ∂p
hp|x̂|xi = xhp|xi = x Cste e−ipx/~ = i~ ∂p ∂
hx|pi). Donc b̂ = −x̂. On déduit que
−i

B̂ (p0 ) = exp ~ (−x̂) p0 .
4. On déduit que e e
Â B̂
= eB̂ eÂ e[Â,B̂ ] et donc que e−Â e−B̂ eÂ eB̂ = e[Â,B̂ ] . Dans le cas
= T̂ (x0 ) = exp −i p̂x0 , soit Â = −i
Â

qui nous interesse, e
~ ~
p̂x0 , et eB̂ =
−i
(−x̂) p0 , soit B̂ = −i ˆ , on ob-

B̂ (p0 )h = exp
i ~ ~
(−x̂) p0 . D'après [x̂, p̂] = i~Id
−i
x0 −i

tient Â, B̂ = − p0 [p̂, x̂] = −ix0 p0 /~ = −iS/~, donnant la relation
~ ~
recherchée.
S/ (2π~) est le nombre de cellules de Planck (2π~) contenues dans la surface S. Ce
nombre intervient dans la formule de Weyl, voir cours.
Exercice 2.2.7 page 111 (voir aussi [DGLR89] p.378).

1

1. On a En = ~ω n + 2
1 1
Pn = exp (−En /(kB T )) = e−α(n+1/2)
Z Z
avec
~ω Θ
α= =
kT T
P∞
Calcul de la constante Z : il faut P = 1 (normalisation des probabilités). On
P n=0 −αn
n
a la série géométrique Sα := n≥0 e = 1−e1−α donc
∞
X 1 −α/2
1= Pn = e Sα
n=0
Z
donc
e−α/2
Z=
1 − e−α
Alors
∞
X ~ω X −α(n+1/2) 1
hEx i = Pn En = e n+
n=0
Z n 2
!
~ω −α/2 X −αn 1 X −αn
= e ne + e
Z n
2 n
Il nous faut calculer
X d e−α
ne−αn = − (Sα ) =
n≥0
dα (1 − e−α )2
Donc
1 + e−α
α
~ω ~ω kΘ Θ
hEx i = = coth = coth
2 1 − e−α 2 2 2 2T
2. Pour l'énergie moyenne d'un atome à 3 dim :
hEi = hEx i + hEy i + hEz i

kΘ Θ
=3 coth
2 2T
R
Pour une mole d'atomes, l'énergie moyenne est hU i = N hEi avec N = k
nombre
d'Avogadro. Alors
3RΘ Θ
hU i = coth
2 2T
et donc la capacité calorique molaire est
2
d hU i 3RΘ (−1) −Θ Θ 1
C= = = 3R
2 sinh2 2T
Θ 2 Θ

dT 2T 2 2T sinh 2T
1 −Θ/T
Pour T Θ, on a C∼ T2
e → 0. On remarque que
r
~ω ~ K
Θ= =
k k m
où K est la constante de raideur de la liaison entre atomes. Donc Corps dur
(diamant)⇔K élevé ⇔Θ élevé.
Pour T Θ on a C → 3R. Cette valeur s'obtient aussi par le théorème d'équiparti-
tion de l'énergie (valable en mécanique classique pour des Hamiltoniens quadratiques
seulement).
B.2. CHAPITRE 2 373
1. On a λx = 2L/a, λy = 2L/b, λz = 2L/d, avec a, b, d ∈ N∗ kx = λ2πx =

entiers. Donc
π
1/2
L
a, etc.... La fréquence de ce mode (a, b, d) est ωa,b,d = ck = c kx2 + ky2 + kz2 =
2 1/2
a b2 d2
πc L2 + L2 + l2 . L'énergie du vide quantique dans la boite est alors (en pen-
sant aux deux états de polarisations possibles)
X 1
E (l) = 2 ~ωa,b,d
a,b,d>0
2
La divergence de E (l) est due aux hautes fréquences ω ; appelée divergence ultra-
violette.
2. On a
X
E (l) = ~ ωe−ω/ωc
a,b,d>0
1/2
πc l 2

et ωa,b,d = l L
(a2 + b2 ) + d2 . Comme (l/L) 1, on peut traiter a, b comme
des variables continues dans la somme (approximation de Riemann d'une intégrale),
et donc
XZ
E (l) ' ~ da db ωe−ω/ωc
d>0
Ensuite, on utilise des coordonnées polaires (a, b) → (u, θ), c'est à dire (a2 + b2 ) = u2
et dadb = ududθ, et θ = 0 → π/2. Alors
π X Z ∞
E (l) ' ~ du u ωe−ω/ωc
2 d>0 0
1/2
πc l 2

Finalement, le changement de variable u→ω= l L
u2 + d2 , donne ωdω =
πc 2

udu L
et
π L 2 X Z ∞
E (l) ' ~ dω ω 2 e−ω/ωc
2πc d>0 ω0
Z ∞
~L2 X d2
= dω e−αω
2πc2 d>0 dα2 ω0
~L2 X d2

1 −αω0
= e
2πc2 d>0 dα2 α
πc
−e−α l
πc d
πc −αω0
e−α
P P
avec ω0 = d, et α = 1/ωc . Ensuite d>0 e = d>0
l = πc =
l 1−e−α l
1
πc (suite géométrique). Donc
eα l −1
~L2 d2

1 1
E (l) = πc
2πc2 dα2 α eα l − 1
~cπ 2 L2 d2

1
=
2l3 dx2 x (ex − 1)
avec x = απc/l = πc/(ωc l).

3. Ensuite
1 1 1 1 1
x2 + O x3

= 2− + −
x (ex − 1) x 2x 12 30 × 24
donc
d2

1 6 1 1
= − 3− + O (x)
dx2 x
x (e − 1) x 4 x 15 × 24
et
4 3 !
~cπ 2 L2

ωc l ωc l 1
E (l) = 6 − − + O(1/ωc )
2l3 πc πc 15 × 24
~cπ 2 L2

ωc 4 ω 3
c 1 1
= 6 l− − + O(1/ωc )
2 πc πc 15 × 24 l3
4. On a
~cπ 2 L2

ωc 4 1 1 1
U (l) = E (l) + E (L − l) = 6 L− + + ...
2 πc 15 × 24 l3 (L − l)3
~cπ 2 L2

ωc 4 1 1
' 6 L− + ... ,
2 πc 15 × 24 l3
2 2
FCasimir (l) = − dU = − ~cπ2 L 3 1

pour L l. Donc
dl 15×24 l4
+ ... et pour ωc → ∞ , les
termes suivants s'annullent, donc
~cπ 2 L2 1
FCasimir (l) = −
240 l4
1. Utiliser la relation (A.2.1), qui s'applique car [q̂, p̂] = i~ Id, [a, a+ ] = Id.
2. Utiliser les expressions (2.2.29), et pour montrer la dernière ligne, le fait que a|0 >=
n
0, donc exp(−za)|0 >= n (−za)
P
n!
|0 >= |0 >.
αa
3. Calculer d(e |qp >)/dα en α = 0.
B.2. CHAPITRE 2 375

1. On vérie que cette représentation par des matrices 3 × 3 vérie bien les règles
+ + +
(2.2.15). Par exemple, on calcule : [ϕ (a) , ϕ (a )] = ϕ (a) ϕ (a ) − ϕ (a ) ϕ (a) =

... = ϕ Iˆ ; etc...
 
0 1 0
1 + 0 0 1 0
2. On a Q̂ = √ (a + a ), donc ϕ −iα Q̂ = (−iα ) √  0 0 1 , et ϕ exp −iα Q̂ =
2 2
0 0 0
 0
√ 02

1 (−iα ) / 2 −α /4√
exp ϕ −iα0 Q̂ = 0 1 (−iα0 ) / 2 . Pour cette dernière égalité,
0 0  1   
0 1 0 0 0 1
2
on utilise le fait que, posant M =  0 0 1 , on a M =  0 0 0  et
0 0 0 0 0 0
n
M = 0 pour n > 2. Donc exp (λM ) = n≥0 n! (λM ) = M + λM + 12 λ2 M 2 =
n 1 0
P
 
1 λ λ2 /2
 0 1 λ . On procède de même pour les quatre autres expressions. Finale-
0 0 1
ment, on calcule le produit de matrices,
1 :

ϕ exp −iγ 0 Iˆ ϕ exp −iα0 Q̂ ϕ exp −iβ 0 P ˆ ϕ (exp (−iθn̂))
1 A0 B 0
 
=  0 e−iθ C 0 
0 0 1
1 2 i 0

A0 = −z 0 e−iθ , B 0 = −iγ 0 + |z 0 | − = z 2 , C 0 = z 0
2 2
que l'on identie avec (2.2.37), pour obtenir (2.2.34).
Exercice 2.3.2 page 133 EF ' 3eV , VF ' 106 m/s.
Exercice 6 page 133

1. Un mode occupe le volume élémentaire
∆3~x∆3~k = (2π)3 dans l'espace de phase
~x, ~k . Considérons un intervalle de fréquence dν . D'après ω = 2πν = ck , cela
correspond à
2π
dk = dν
c
et à un volume dans l'espace ~k de
Vk = 4πk 2 dk
1. Pour ces produits de matrices, ainsi que ces exponentielles de matrices, on aura intérêt à utiliser un
logiciel de calcul formel (comme Maple, Mathematica ou Xcas qui est un logiciel gratuit et libre).
(volume de la sphère de rayon k et épaisseur dk ). Donc dans un volume V et un

intervalle de fréquence dν contiennent
dn = 2 (V Vk ) / (2π)3
modes. Le facteur 2 tient compte des deux états de polarisation possibles d'un mode
(droite/gauche). Donc

2V 4πk 2 (2π)
c
dν 8πV ν 2
dn = = dν
(2π)3 c3
P∞ 1
2. On a
PhNmode i = N =0 PN N , avec la probabilité PN = Z exp (−EN /kT ). Comme
1 = N PN , on déduit que la constante Z est donnée par
X X
Z= exp (−EN /kT ) = exp (−α (N + 1/2))
N N
X
= e−α/2 e−αN
N ≥0
| {z }
S
e−αN = 1
P
avec α = (~ω) /kT et la série géométrique S= N ≥0 1−e−α
. Alors
∞
X
hNmode i = PN N
N =0
1X
= N exp (−α (N + 1/2))
Z N

1 −α/2 dS 1 dS 1
= −α/2 e − = − = α
e S dα S dα e −1
1
= hν/(kT )
e −1
appelée distribution de Bose-Einstein. Ensuite
dN dN dn dn
= =< Nmode >
dν dn dν dν
1 3
3. u(ν) = V
(hν) dN
dν
= 8πhν
hν
c3 e kT −1
Exercice 7 page 135

B.2. CHAPITRE 2 377
1. Un mode occupe ∆3~x∆3~k = (2π)3 . Considérons un intervalle de fréquence dν . Pour

les ondes S, cela correspond à dk =
2π
vS
dν , et un volume dans l'espace ~k de Vk =
2
4πk dk (volume de la sphère de rayon k et épaisseur dk ). Donc dans un volume V
v 3
et un intervalle de fréquence dν = S dk contiennent dnS = 2 (V Vk ) / (2π) modes.
2π
Le facteur 2 tient compte des deux états de polarisation possibles d'un mode. Donc
2 (4π) V ν 2 dν
dnS =
vs3
de même pour les ondes P :
(4π) V ν 2 dν
dnP =
vP3
2. D'après l'hypothèse d'équidistribution entre les modes, le rapport d'énergie est égal
au rapport du nombre de modes :
3
EP dnP 1 vS 1 1
= = = 3 '
ES dnS 2 vP 2 (1, 73) 10
Exercice 9 page 136

1. On a T̂q B̂p = exp (−iq p̂/~) exp (ipx̂/~). Utiliser (A.2.1), avec A = −iq p̂/~, B =
ipx̂/~, [A, B] = (qp/~2 ) (−i~) donnant eA eB = eB eA e[A,B] donc T̂q B̂p = e−iS/~ B̂p T̂q ,
avec S = qp qui est la surface du carré dans l'espace de phase concerné par les
translations. En terme de mécanique analytique, c'est une action. Les opérateurs
commutent lorsque S = nh, n ∈ N c'est à dire lorsque la surface contient un nombre
entier de cellules de Planck.
1

A B
2. Utiliser (A.2.1) :e e = eA+B e 2 [A,B] . On a B = i(P2 − P1 ) Q̂ − Q1 − i(Q2 −

Q1 ) P̂ − P1 , A = i(P3 − P2 ) Q̂ − Q2 − i(Q3 − Q2 ) P̂ − P2 , et 21 [A, B] =
→
− → −
− 2i (P2 −P1 )(Q3 −Q2 )+ 2i (Q2 −Q1 )(P3 −P2 ) = − 2i 23∧ 12 = −i S/~ :produit vectoriel
1
et S est la surface hachurée (dans les unités q, p). On a S/~ =
2
(P2 Q3 − P1 Q3 + P1 Q2 − Q2 P3 + Q
On a donc

D̂2,3 D̂1,2 = e−iS/~ exp (A + B) = e−iS/~ exp i(P3 − P1 )Q̂ − i(Q3 − Q1 )P̂ exp (i (−P2 Q1 + Q2 P1 − P3 Q
. Par ailleurs,

D̂1,3 = exp i(P3 − P1 ) Q̂ − Q1 − i(Q3 − Q1 ) P̂ − P1

= exp i(P3 − P1 )Q̂ − i(Q3 − Q1 )P̂ exp (i (−P3 Q1 + Q3 P1 ))
. Finalement, D̂2,3 D̂1,2 = D̂1,3 e−iS/~ ei2S/~ = e+iS/~ D̂1,3 .

3. Par récurrence.
Exercice 10 page 137

1. H(q, p) = H(q0 , p0 ) + (q − q0 ) . ∂H
∂q
+ (p − p0 ) . ∂H
∂p
+ o(∆q, ∆p), avec ∆q = q − q0 ,
∆p = p − p0 .

2. Ĥ|q0 , p0 >' H(q0 , p0 ) Id + (q̂ − q0 ) . ∂H
∂q
+ (p̂ − p 0 ) . ∂H
∂p
|q0 , p0 >. (On a remplacé
q et p par les q̂ et p̂).
3. Ainsi

Û (∆t)|q0 , p0 > = exp −iĤ∆t/~ |q0 , p0 >

∂H ∂H
' exp −i H(q, p) Id + (q̂ − q0 ) . + (p̂ − p0 ) . ∆t/~ |q0 , p0 >
∂q ∂p
Û (∆t)|q0 , p0 >= exp(−iE∆t/~)D̂0,1 |q0 , p0 >, avec q1 − q0 = ∂H

donc
∂p
∆t et p1 −
∂H
p0 = − ∂q ∆t sont les déplacements du point (q0 , p0 ) par la dynamique classique, et
E = H(q0 , p0 ) est l'énergie du point classique, conservée au cours du temps. Sur la
durée ∆t innitésimale, le paquet d'onde est donc translaté dans l'espace de phase,
comme le point classique.
4. Alors en sommant des déplacements innitésimaux,
N −1
Û (t)|q0 , p0 >' exp (−iEt/~) Πi=0 D̂i,i+1 |q0 , p0 >= exp (−iEt/~) exp (iS/~) D̂0,M (t) |q0 , p0 >
, où M (t) est l'évolution classique.

RT 0
iE(t+t0 −t0 )/~ −iEt0 /~ T +t iEt/~
Û (t0 )|φ >= 0
R
5. On a
0
e Û (t + t )|q 0 p 0 >= e t0
e Û (t)|q0 , p0 >
0
dt = e−iEt /~ |φ >, à condition que e
iET /~
Û (T )|q0 , p0 >= |q0 , p0 >, car alors
Z T +t0 Z t0
iEt/~
e Û (t)|q0 , p0 > dt = eiEt/~ Û (t)|q0 , p0 > dt
T 0
.
6. D'après ci-dessus, la condition s'écrit : exp (iS/~) D̂0,M (T ) |q0 , p0 >= |q0 , p0 >, où S
est la surface de la trajectoire (action). Or M (T ) = M (0), donc il faut exp (iS/~) =
1, soit
S = hn, n ∈ N
appelée règle de quantication de Bohr-Sommerfeld : la surface doit contenir
un nombre entier de cellules de Planck.
Pour l'oscillateur Harmonique, cela donne (voir section précédente) En ' ~ωn. (Il
manque l'indice 1/2).
R +∞ R +∞ 1
Exercice 2.3.3 page 140 J = −∞ −∞ 2π~
(Voir [CBF] p571). On adx dp |x, p, σihx, p, σ| =
R +∞ R +∞ 1 2

−∞ −∞ 2π~
dx dp D(z)|0ih0|D+ (z). Or D(z)|0i = exp − |z|2 exp (za+ ) exp (−za) |0i =
2

exp − |z|2 eza |0i. Or eza |0i = n≥0 z n!
+ + n a+n n
|0i = n≥0 √z n! |ni.
P P
B.2. CHAPITRE 2 379
0
z n z n −|z|2
dxdp √ 0 2
|nihn0 |. Mais d'une part dx∧dp = dQ∧dP et dQdP z n z n e−|z|
P R R
Donc J= n,n0 h n!n0 !
e h
=
2 0 0
dρdθ e−ρ ρn+n eiθ(n−n ) = n! δn,n0 en coordonnées polaires (ρ, θ) et intégration sur ρ par
R
P ˆ
parties. Donc J= n≥0 |nihn| = I . De façon plus élégante, cette relation de fermeture est
une simple conséquence du lemme de Schur (théorie des groupes) voir [Per86] p.15.
Exercice 8 page 136 @@ revoir, voir correction de JJB @@
1. Le Hamiltonien classique est :
p~2 e2
H(~x, p~) = −
2m 4πε0 r
Soit une valeur E<0 xée. On cherche le volume de l'espace de phase des points :
ΣE = {(~x, p~)/H(~x, p~) ≤ E}
On utilise naturellement les coordonnées sphériques : p~ ≡ (p, ϕp , θp ) et ~x ≡ (r, ϕ, θ).

x = r2 dr sin θdθdϕ)
Alors (utilisant d~
Z
V (E) = V ol (ΣE ) = d~xd~p
ΣE
Z
2
= (4π) 2
r2 p2 drdp
p e2
E≥ 2m − 4πε
0r
∞
e4
Z
2 1
= (4π) p2 dp
3 (4πε0 )2
3
0 p2
2m
−E
e4 m3/2 π e4 m3/2 π

2
= (4π) √ = √
3 (4πε0 )2 4 2(−E)3/2 12 2 (−E)3/2 ε20
Ensuite d'après (2.3.2),
V (E) e4 m3/2 π
N (E) = = √
h3 12 2 (−E)3/2 ε20 h3
4 3/2 dN
2. On déduit la densité de niveaux semi-classique √ e m 5/2π 2 , qui
ρsc (E) = =
8 2(−E) ε0 h3 dE
est une expression valable lorsque les niveaux deviennent denses, c'est à dire pour
E → 0− .
4
3. A partir de l'expression Enexact = − nε12 avec ε1 = me
2~2
= 13, 6 eV , on a (sans oublier
2 2 2
la dégénérescence du niveau En qui est n ) dN = n dn donc dN/dE = n dn/dE =
1 m3/2 e6
.
2 23/2 ~3 (−E)5/2
@@
B.3 Chapitre 3
Exercice 3.1.2 page 149 on a ~x = R−1~x0 et la valeur de la fonction d'onde n'est pas
0 0
changée lors de la rotation : ψ (~x ) = ψ(~x), donc :
ψ 0 (~x0 ) = ψ(R−1~x0 ) (B.3.1)
(situation semblable au cas de la translation, gure (2.1.1)).

0 0
x ) =< ~x0 |R|ψ >=<
On peut obtenir le même résultat avec la notation d'opérateurs : ψ (~
R+~x0 |ψ >=< R−1~x0 |ψ >= ψ(R−1~x0 ), où on a utilisé le fait que l'opérateur rotation est
unitaire.
Exercice 3.3.1 page 169 Électrons bidimensionnels dans un champ magnétique

constant, avec un potentiel périodique : le spectre fractal de Hofstadter. 2

1. B ~ = (∂x Ay − ∂y Ax ) ~ez = B~ez .
~ = rot A Alors
1 ~
2
H(x, px , y, py ) = p~ − eA + V (x, y)
2m
2 2 !
1 eB eB
= px + y + py − x + V (x, y)
2m 2 2
2. Pour les unités, [Q] ≡ [p ]

√ x ≡ √ [px ] ≡ 1 (pensant à d~
p ~)
≡ e~v ∧ B et
[ ~eB ]
[px ]m[px ]/m

dt
[px ]
[q] ≡ [eBX] ≡ ([px[p]/m)m
x]
≡ 1. Donc Q̂, P̂ , q̂, p̂ sont sans dimension. A partir de
[x̂, p̂x ] = i~ et [ŷ, p̂y ] = i~, on calcule

h i
Q̂, P̂ = i, [q̂, p̂] = i ~ef f
et les autres commutateurs sont nuls. On a nalement :
φ0
2π~ef f = .
φ
3. On obtient
~ω p p
2 2

H (Q, P, q, p) = P + Q + V −X p + hef f P , −X q + hef f Q
2
2. Ce spectre a été étudié dans D. Hofstadter, Energy levels and wave functions of Bloch electrons in
rational and irrational magnetic elds Phys.Rev.B 14 ,2239, (1976). Il est (partiellement) observé expé-
rimentalement dans : Albrecht et al. Evidence of the Hofstadter Fractal energy spectrum in the Quantized
Hall Conductance Phys.Rev.Lett. 86 ,147 (2001). Chercher hofstadter buttery dans Google.
B.3. CHAPITRE 3 381
4. Si V = 0, alors H (Q, P, q, p) = ~ω
2
(P 2 + Q2 ) est un oscillateur harmonique. Ses
niveaux d'énergie (niveaux de Landau) sont donc :

1
En = ~ω n + , n = 0, 1, 2, . . .
2
mais les variables (q, p) n'interviennent pas. Cela montre que chaque niveau est
inniment dégénéré (comme la dimension de l'espace de Hilbert des opérateurs
q̂, p̂).
(La suite est facultative)
5. La périodicité est 1 selon q et selon p. Les générateurs des translations sont respec-
tivement (p̂, −q̂), et [q̂, p̂] = i ~ef f . Par conséquent
T̂q = exp(−ip̂/~ef f ), T̂p = exp(iq̂/~ef f ).
−i2π
6. La relation de Glauber donne
h i T̂q T̂p = e[−ip̂/~ef f ,iq̂/~ef f ] T̂p T̂q = e 2π~ef f T̂p T̂q . Donc
T̂q , T̂p = 0 si et seulement si
1
=N ∈N : entier
2π~ef f
φ 1
Cela donne
φ0
= 2π~ef f
=N . Autrement dit, il doit y avoir un nombre entier de
2
quanta de ux à travers la surface élémentaire X .
2
−7
A.N. : N =
X eB
h
, avec X = 2.10 m, B = 0, 21 T., h = 6, 6.10−34 J.s., e =
−14 −19
1, 6.10−19 C , N = 10 6,6.10
4 1,6.10
−34
0,21
' 2.
7. Avec cette condition, on peut (comme dans la théorie de Bloch), chercher les états
stationnaires parmi les vecteurs propres communs de T̂q et T̂p qui vérient, T̂q |ψi =
eiθ1 |ψi et T̂p |ψi = eiθ2 |ψi. Les états stationnaires ainsi trouvés |ψn (θ1 , θ2 )i vont
dépendre des paramètres continus (θ1 , θ2 ) et d'un indice discret n. Cela donnera un
spectre en bandes. Chaque bande est indicée par n.
8. Avec les notation ci-dessus, appelons H (θ1 , θ2 ) l'espace des fonctions de Bloch véri-
iθ iθ iθ
ant T̂q |ψi = e 1 |ψi et T̂p |ψi = e 2 |ψi. Comme T̂q |ψi = e 1 |ψi, la fonction d'onde
ψ(q),est périodique (à une phase près), de période 1. Par conséquent, sa transfor-
1 iθ
mée de Fourier ψ̃(p) est discrète, avec un pas ∆p = hef f = . Or T̂p |ψi = e 2 |ψi,
N
ce qui signie que ψ̃(p) est elle même périodique, de période 1. Dans l'intervalle
p ∈ [0, 1[ il y a 1/∆p = N composantes, et l'état|ψi n'a donc que N composantes
indépendantes. L'espace des fonctions de Bloch H (θ1 , θ2 ) est donc de dimension N .
Avec la perturbation V , le premier niveau de Landau perd sa dégénérescence et il
apparait donc N bandes de Bloch. (Cela est aussi vrai pour les autres niveaux de
Landau).
Bz
y
1/5 1/4 1/3 1/2 2/3 3/4 1

X x Flux magnetique −1
Φ0 /Φ
9. La condition 2π~ef f =
N
1
est exceptionelle. Par contre, si 2π~ef f =
h
eBX 2
= φφ0 =
a
b
∈ Q, on considère une cellule de base de côtés X et Y = aX , donc de surface
XY = aX 2 . Pour cette cellule de base, la constante eective 2π~0ef f = eB(XY
2π~
)
= 1b .
0
−1
Donc 2π~ef f = b ∈ N est entier, les résultats précédents s'appliquent et on
déduit que le spectre est formé de b bandes. Ainsi lorsque φ/φ0 varie continuement,
le dénominateur b (et donc le nombre de bandes) varie discontinuement, ce qui
explique l'aspect fractal du spectre obtenu. Voir gure.
B.4 Chapitre 4
Exercice 11 page 188 le résultat (4.5.1) montre que cette relation est vraie pour les
~ˆ.
h i
trois vecteurs de base ei de R3 , ce qui peut s'écrire Ŝei , Ŝej = i~Ŝei ∧ej = i~ei ∧ ej .S
~ˆ =
h i P h i P h i P
Alors ŜU~ , ŜV~ = i,j Ui Ŝ ei
, Uj Ŝ ej
= i,j Ui U j Ŝe i
, Ŝej
= i~ i,j Ui Uj ei ∧ ej .S
~ˆ.

~ ~
i~ U ∧ V .S
Exercice 4.4.1 page ??

1. Soit un état de spin 1/2 quelconque noté :
|ψ >= a|+z > +b|−z >
avec a, b ∈ C. Si b 6= 0, alors on écrit |ψ >= b|ϕi avec |ϕi = z|+z > +|−z >
et z = a/b. De façon très générale, les probabilités associées au résultat d'une
mesure sur un état quantique, ne change pas si l'état est multiplié par une constante
complexe. Plus précisement, si |ψ >= b|ϕi, et Â est une observable (opérateur) :
hψ|Â|ψi |b|2 hϕ|Â|ϕi hϕ|Â|ϕi

= 2 =
hψ|ψi |b| hϕ|ϕi hϕ|ϕi
Par conséquent |ψi et |ϕi décrivent les même états physiques

B.4. CHAPITRE 4 383
2. On calcule

hψ|Sbx |ψi 1 ~ 0 1 a ~ 2< (z)
sx = = 2 2 a, b =
hψ|ψi |a| + |b| 2 1 0 b 2 1 + |z|2
de même :
hψ|Sby |ψi ~ 2= (z) hψ|Sby |ψi ~ |z|2 − 1

sy = =− , sz = =
hψ|ψi 2 1 + |z|2 hψ|ψi 2 1 + |z|2
~ 2
s2 = hSx i2 + hSy i2 + hSz i2

Donc s2 = s2x + s2y + s2z = 2
. (Ne pas confondre avec
hS i = ~2 12 12 + 1 ).
2

~ ~ ~
3. En coordonnées sphériques (s, θ, ϕ), sx = 2
sin θ cos ϕ, sy = 2
sin θ sin ϕ, sz = 2
cos θ,
on obtient
1 + 2 s~z 1 + cos θ θ
zz = sz = = cotg2
1−2~ 1 − cos θ 2
et !
2
s
xsy 1 cos2 2θ
z = 1 + |z| −i = 1+ (cos ϕ − i sin ϕ) sin θ
~ ~ 2 sin2 2θ
donc z = cotg 2θ e−iϕ .

4. D'après les schémas :
Exercice 12 page 189
1. On a vu en cours, que les générateurs Ŝx , Ŝy , Ŝz s'expriment respectivement dans la

base o.n. (|+z i, |−z i) par les matrices de Pauli ~
2
σx , ~
2
σ y , ~
2
σz ; ce sont des matrices
hermitienne 2 × 2 de trace nulle. Par conséquent, l'opérateur ŜU ~ = Ux Ŝx + Uy Ŝy +

ˆ
~ .S
~ ∈ R, U ~ = (Ux , Uy , Uz ) ∈ R (générateur des rotations du spin 1/2)
3
Uz Ŝz = U
~
s'exprime par la matrice σU avec
2 ~

~ .~σ = Uz Ux − iUy
σU~ = Ux σx + Uy σy + Uz σz = U
Ux + iUy −Uz
2. Soit M ∈ SU (2). On a (a) : M + = M −1 et (b) : det(M ) = 1. Écrivons M =

exp (−iλG), avec λ ∈ R, et G matrice 2 ∗ 2 qui est un générateur ,i.e. élément de
+
l'algèbre de Lie su(2). Pour λ → 0, la relation (a) donne (1 + iλG ) = (1 + iλG),
+
soit G = G . La relation (b) donne 1 = det(M ) = exp (−iλ T race(G)), donc
T race(G) = 0. (La relation matricielle très générale det(M ) = exp (−iλ T race(G)),
se démontre en écrivant G dans une base propre, où elle est diagonale). Ainsi l'al-
gèbre de Lie su(2) du groupe SU (2) est formé par les matrices hermitienne 2 × 2 de

uz ux − iuy
trace nulle forme. Comme expliqué en (1), G est de la formeG = =
ux + iuy −uz
~u.~σ , et donc les matrices de Pauli forment une base de cette algèbre.
3. Soit ϕ : Ô → MÔ l'application qui à un opérateur de Hspin associe la matrice 2 × 2,

qui est son expression dans la base (|+z i, |−z i). On a vu que ϕ ŜU~ = ~
σ ~ et
2 U
que
ϕ :Rspin →
su(2)
est un isomorphisme d'algèbres de Lie. Pour R̂ ∈ Rspin , on

−iŜu
~ α −iŜu
~ α
= exp − 2i σU~ qui montre que

a ϕ R̂ = ϕ exp = exp ϕ
~ ~
ϕ : Rspin → SU (2) est un isomorphisme de groupes.
Exercice 13 page 191 Voir eq(B.4.1). Cette fois ci l'opérateur de rotation agit à la fois
0
sur l'espace ordinaire et sur le spin : ψ± (~x0 ) =< ~x0 , ±|R̂|ψ, ~s >=< R+~x0 , R+ ± |ψ, ~s >=<
R−1~x0 , R−1 ± |ψ, ~s >, c'est à dire :
−1
0
(~x0 ) ~x0 )

ψ+ −1 ψ+ (Respace
0 = R̂spin (B.4.1)
ψ− (~x0 ) −1
ψ− (Respace ~x0 )
B.5 Chapitre 7
Exercice 6.5.1 Diusion nucléon-pion page 277. D1/2 ⊗D
On a vu que
h 1/2 =i D0 ⊕D1
qui est la décomposition de Clebsh-Gordan. ~ = 0 (in-
Si Ĥ est un opérateur vériant Ĥ, S
variance par rotation) alors d'après le Lemme de Shur, Ĥ exprimé dans la décomposition
D0 ⊕ D1 est de la forme bloc:
λIˆ
 
0
Ĥ =  |{z}
0 µIˆ  , λ, µ ∈ R
|{z}
D0 D1
A.Iˆ + B.S
~1 .S
~2 = A + B ~2 J(J + 1) − 3

Or on a vu que l'opérateur s'exprime dans cette
2 2
même décomposition par
E0 = A − 34 B~2

0
0 E1 = A + 14 B~2
Il sut donc de prendre A, B tels que A − 43 B~2 = λ et µ = A + 14 B~2 .

Plus généralement, pour le couplage Dj1 ⊗ Dj2 = ⊕jj=|j
1 +j2
2 −j1 |
Dj , un opérateur invariant
B.5. CHAPITRE 7 385
s'exprime comme λj Iˆ sur chaque espace Dj de la décomposition. Il y a donc N = j2 + j1 −

|j2 − j1 | paramètres indépendant λj . Comme l'opérateur de Casimir S ~ 2 est ~2 j (j + 1) Iˆ
sur un espace Dj (et distingue donc les espaces Dj ) on déduit que un opérateur invariant
2 N −1
peut s'exprimer sous la forme Ĥ = µ0 Iˆ + µ1 S ~ 2 + µ2 S
~ 2 + . . . + µN −1 S~2 . Il est
k
possible de relier les (µj )j et les (λj )j . De façon équivalente on pourrait utiliser S ~ 1 .S
~2
k
à la place de S ~ 2 , les coecients seraient diérents.
Exercice 6.5.2 Diusion nucléon-pion page 283.
1. On aHj=1 ⊗ Hj=1/2 = Hj=1/2 ⊕ Hj=3/2 , de dimensions 3×2 = 2+4 = 6 (d'après

dimHj = 2j + 1).
2. On obtient :
3 3
|pπ + i = | ; i
2 2
0
p 3 1 p 1 1
|pπ i = 2/3| ; i − 1/3| ; i
2 2 2 2
−
p 3 1 p 1 1
|pπ i = 1/3| ; − i − 2/3| ; − i
2 2 2 2
+
p 3 1 p 1 1
|nπ i = 1/3| ; i + 2/3| ; i
2 2 2 2
0
p 3 1 p 1 1
|nπ i = 2/3| ; i + 1/3| ; − i
2 2 2 2
− 3 3
|nπ i = | ; − i
2 2
3. Dans un espace de représentation irréductible, Û doit agir comme l'identité à une

constante complexe près (Lemme de Shur). Donc dans la décomposition Hj=1/2 ⊕
Hj=3/2 :
A3/2 Iˆ4

0
Û ≡H3/2 ⊕H1/2
0 A1/2 Iˆ2
avec A3/2 , A1/2 ∈ C.
4. On obtient :
2 2
σ pπ + → pπ + hp, π + |Û |p, π + i

= = A3/2
√ √ 2
− 0
0 −
2 2 2
σ pπ → nπ = hn, π |Û |p, π i = A3/2 − A1/2
3 3
2
2 1 2
σ pπ − → pπ − hp, π − |Û |p, π − i

= = A3/2 + A1/2
3 3
2 2
(a) si A3/2 A1/2 , alors σ (pπ + → pπ + ) = A3/2 , σ (pπ − → nπ 0 ) ' 2
9
A3/2 ,
2
σ (pπ − → pπ − ) ' 19 A3/2 .
2
(b) si A3/2 A1/2 , alors σ (pπ + → pπ + ) ' 0, σ (pπ − → nπ 0 ) ' 2
9
A1/2 , σ (pπ − → pπ − ) '
4 2
9
A1/2 .
2
(c) siA3/2 = A1/2 , alors σ (pπ + → pπ + ) = A3/2 , σ (pπ − → nπ 0 ) ' 0, σ (pπ − → pπ − ) '
2
A3/2 .
5. L'expérience donne :
σ (pπ + → pπ + )
= 4.33 ' 4.5
σ (pπ − → nπ 0 )
σ (pπ + → pπ + )
= 8.47 ' 9
σ (pπ − → pπ − )
On est donc proche de la situation A3/2 A1/2 .
6. Durée de vie τ ' ∆E ~

= 6.510−24 s. L'espace de degré interne d'isospin de cette
particule ∆ est Hj=3/2 , de dimension 2j + 1 = 4. D'après l'écriture |j = 3/2, m =
3/2i = |p+ , π + i, etc... on déduit que |∆++ i = |j = 3/2; m = 3/2i, |∆+ i = |j =
3/2; m = 1/2i,|∆0 i = |j = 3/2; m = −1/2i,|∆− i = |j = 3/2; m = −3/2i.
Index
adjoint, 28 constante de Planck, 72
algèbre de l'oscillateur Harmonique, 127 corps noir, 133
algèbre de Lie complexiée, 103 corrélées, 148
algèbre de Lie d'opérateurs, 103 courbure de la connexion, 161
algèbre de Lie de l'oscillateur Harmonique,
102 degré de liberté spatial, 146
algèbre de Weyl-Heisenberg, 124 degrés de liberté, 112
analyse micro locale, 76 densité d'état, 56
analyse semi classique, 76 densité de niveaux, 132
approximation harmonique, 99 densité de probabilité, 61, 63, 84
approximation linéaire, 99 diusion isotrope, 305
autoadjoint, 29 diusion s, 305

dimension nie, 41
base d'énergie, 63 dimension innie, 41
base de Fourier, 43, 52 dispersion, 31
base de position, 50, 63 distribution de Bargmann, 141
base orthonormée, 41 distribution de Dirac, 50, 74
Bohr-Sommerfeld, 138 distribution de Husimi, 141
boost, 94 domaine de dénition d'opérateur, 27
bra, 74 double fentes de Young, 58
bra, 25 double puits de potentiel, 38
dualité onde-corpuscule, 79
capacité calorique des matériaux, 110 dé-corrélés, 149
catastrophe infra-rouge, 117 décomposition de Cartan, 104
catastrophe ultraviolette, 117 dégénérescence, 47
cellule de Planck, 130 dégénérée, 55
changement de Jauge, 153 dérivée covariante, 153, 158
changement de référentiel Galiléen, 94 détecteur, 83
choix d'une Jauge, 162
Chromo-dynamique, 166 E.D.P., 75
Classique, 10 Eet Aharonov-Bohm, 155
collapse, 340 eet photo-électrique, 115
collapse de la fonction d'onde, 80, 84 ensemble statistique d'états, 341
commutateur, 70 entropie de Shannon, 349
condition de Laue, 302 entropie statistique, 132
connexion sur un bré vectoriel, 157 environnement, 57
387
388 INDEX
espace bosonique, 223 group, 90

espace de Hilbert, 74 groupe, 20
espace de Hilbert, 24 groupe commutatif, 20
espace des phases, 13 groupe de l'oscillateur Harmonique, 127
espace des phases, 31 groupe de Lie de dimension 1, 90
espace dual, 25 groupe de Weyl-Heisenberg, 125
espace fermionique, 223 Groupes de Lie, 90
espace produit tensoriel, 146, 148 générateur, 88, 90
espace propre, 47, 81 géométrie diérentielle, 156
espace quantique du spin, 175
Hamilton, 10
espace somme directe, 148
Hamiltonien, 13, 26
espace vectoriel, 72
harmonique, 99
espace vectoriel, 20
hasard quantique, 84
espace vectoriel complexe, 20
hermitique, 29
espace vectoriel réel, 20
holonomie, 161
espaces de Sobolev, 25
Everett, 66
impulsion, 13, 72
incertitude, 67
facteur gyromagnétique de l'électron, 117
intégration par partie, 29
fentes de Young, 155
Invariance de Jauge classique, 153
bré vectoriel, 157
Invariance de Jauge quantique, 154
l quantique, 25, 42
inégalité de Cauchy-Schwartz, 74
fonction d'onde, 72
inégalité de Heisenberg, 71
fonction d'onde, 19
inégalité de Schwartz, 71
fonction de Dirac, 49
isolé, 57
fonction de Hermite, 105
isolée, 19
fonction de partition, 352
isomorphisme de Riesz, 25
fonctions d'Hermite, 105
force, 12 ket, 72
force de Casimir, 118 ket, 20
force de Lorentz, 151
force nucléaire faible, 166 l'écart type, 67
force nucléaire forte, 166 Laser, 115
force potentielle, 12 Le problème à deux corps, 14
forme linéaire, 25 lemme de Poincaré, 156
formule de Gutzwiller, 139 Limite classique-quantique, 64
formule de transport, 298 localisation forte de Anderson, 56
formule des traces semi-classique, 139 loi (faible) des grands nombre, 350
fréquence de Bohr, 321 loi de Boltzmann, 110

Loi de Newton, 12
gaz de photons, 133 loi de Planck, 133
gluons, 166 loi de probabilité, 346
gravitation, 168 Loi de Thomas Fermi, 130
INDEX 389
Loi de Weyl, 130 opérateur impulsion, 93

Lorentzienne, 50 opérateur linéaire, 27
opérateur partiel, 355
Many world interpretation, 66
opérateur positif, 343
matrice densité, 341
opérateur rotation, 149
matrices de Pauli, 179
opérateurs d'échelle, 106
Modèle d'Einstein, 110
Opérateurs de création et d'annihilation, 103
modèle de Debye, 112
opérateurs de création et d'annihilation, 106
modèle standard de la physique des parti-
opérateurs diérentiels, 27
cules, 101
opérateurs vectoriels, 145
moment d'ordre n, 69
orthogonales, 23
moyenne expérimentale, 84
oscillateur harmonique, 15, 32, 98
moyenne statistique, 84
multiplicité, 47, 55
paquet d'onde Gaussien, 73
mécanique classique, 12
mécanique quantique, 72
paquet d'onde Gaussien, 21
particule, 19, 72
nabla, 146 particule libre, 54
niveaux d'énergie, 77 phonons, 101
normalisation, 46 photon, 101
normalisé, 24 photons, 115
normalisée, 74 polarisation, 345
norme, 23, 73 polynome d'Hermite, 105
Notation de Dirac, 28 position d'une particule, 12
notation de Dirac, 20 postulat de la mesure, 58
noyau de Schwartz, 51
postulat de la mesure, 80
potentiel Harmonique, 98
observable, 59, 83, 90
potentiel quadratique, 54
observables, 81
potentiel scalaire, 150
onde plane, 51, 72
potentiel vecteur, 150
onde plane, 21
principe d'équivalence, 14
onde stationnaire, 54, 76
ondes sismiques, 135 principe de correspondance, 76, 150
opérateur, 26 principe de superposition, 76
opérateur adjoint, 28 principe de superposition, 30

opérateur autoadjoint, 46 probabilité, 84
Opérateur d'impulsion, 26 problème inverse, 302
opérateur d'évolution, 88 produit scalaire, 73
Opérateur de position, 26 produit scalaire, 23

opérateur de spin, 178 produit scalaire de Hilbert-Schmidt, 27
opérateur de translation, 92 produit tensoriel, 146
opérateur densité, 341 projecteur orthogonal, 44
opérateur densité partiel, 356 projecteur spectral, 59
opérateur identité, 27 propagateur, 88
390 INDEX
Quantication du champ électromagnétique, théorie des connections sur brés vectoriels,

114 156
quantication géométrique, 139 théorie des distributions, 50
Quantique, 10 théorie du champ moyen, 333
quarks, 166 théorème d'Erhenfest, 70
Théorème de Wigner-Eckardt, 286
rang de Schmidt, 358 trace partielle, 355
rayonnement du corps noir, 115 trajectoire de la particule, 12
relation d'incertitude, 71 transformation active, 149
Relation d'incertitude temps-énergie, 91 transformation innitésimale, 91
relation de fermeture, 82 transformation passive, 149
relation de Parceval, 53 transformée de Fourier, 52, 79
relation de polarisation, 89 trivialisation du bré, 162
relativiste, 10
relativité générale, 12 unitaire, 89
relativité restreinte, 12
valeur moyenne, 67, 81
renormalisation, 117
valeur propre, 45
représentation dans l'espace de phase, 139
valeur propre simple, 47
représentation de Heisenberg, 96
variance, 67
représentation de Husimi, 32, 139
vecteur colonne, 42
représentation de Schrödinger, 96
vecteur dual, 25
représentation en impulsion, 52
vecteur dual métrique, 74
représentation en position, 50
vecteur propre, 45
règle de Leibnitz, 163
règle de quantication, 138
xcas, 105
réduction de l'état quantique, 340
réduction du paquet d'onde, 63, 80, 84 Équations de Hamilton, 13
réduction symplectique, 238
élément du groupe, 90
section de l'espace bré, 159 éléments de matrice, 44
simplement connexe, 156 énergie cinétique, 13
spectre, 45 énergie noire, 117
spectre continu, 49 énergie potentielle, 13
spectre du corps noir, 133 équation de Liouville quantique, 345
spectre ponctuel, 56 équation de Schrödinger, 72, 75
sphère d'Ewald, 302 équation de Schrödinger, 10, 30

symbole de l'opérateur, 97 équation de Schrödinger stationnaire, 53, 76
systèmes intégrables, 13 équation linéaire aux dérivées partielles, 75

équations de Maxwell, 112
temps de décohérence, 63 état cohérent, 123
température d'Einstein, 110 état cohérent du champ, 115
théorie de Jauge, 156, 157 état de position, 74
théorie de Jauge électro-faible, 157 état fondamental, 77
INDEX 391
état non polarisé, 345

état physique, 62
état pur, 341
état quantique, 341
état quasi-classique, 73
état stationnaire, 54
états physiques, 74
évolution temporelle, 30
392 INDEX
Bibliographie
[ADR82] A. Aspect, Dalibard J., and Roger G. Experimental tests of Bell's inequalities
using time-varying analysers. Physical Review Letters, 49 :18041807, 1982.
[AM76] Neil W. Ashcroft and D. N. Mermin. Solid State Physics. Brooks Cole,
January 1976.
[Arn76] V.I. Arnold. Les méthodes mathématiques de la mécanique classique. Ed. Mir.
Moscou, 1976.
[Asp76] A Aspect. Proposed experiment to test the non-separability of quantum

mechanics. Physical review D, 14 :19441951, 1976.
[B.] Parisse B. Logiciel libre de calcul formel. Taper xcas dans google.
[Bac] H. Bacry. Lessons on group theory.
[Bas86] J.L. Basvedant. MÃ©canique quantique. Ellipses, 1986.
[BC89] B.H. Bransden and C.J. Joachain. Introduction to quantum mechanics. Long-
man, 1989.
[Bel64] Bell J. On the Einstein-Podolsky-Rosen paradox. Physics, 1 :195200, 1964.
[BPLB13] R. Bach, D. Pope, SH. Liou, and H. Batelaan. Controlled double-slit electron
diraction. New Journal of Physics, 15(3) :033018, 2013.
[CB73] Y. Choquet-Bruhat. Distributions : theorie et problemes. 1973.
[CBF] C. Cohen-Tannoudji, B. Diu, and F. Laloe. Mécanique quantique.
[Cla88] Claude Cohen-Tannoudji. Cours 1988-89-90 : Cohérences quantiques et

dissipation. Cours de Claude Cohen-Tannoudji au Collège de France.
http ://www.lkb.ens.fr/users/laloe/public_html/PHYS/cours/college-de-
france/, 1988.
[CTDRG87] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg. Photons et atomes.

1987.
[CTDRG88] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg. Processus d'interac-

tion entre photons et atomes. ..., 1988.
+
[DEC 96] D. Giulini, E. Joos, C. Kiefer, J. Kupsch, I.O. Stamatescu, and H.D. Zeh.
Decoherence and the Appearance of a Classical World in Quantum Theory.
Springer-Verlag, 1996.
393
394 BIBLIOGRAPHIE
[DGLR89] B. Diu, C. Guthmann, D. Lederer, and B. Roulet. Physique statistique. 1989.
[DM02] Dietrich Dehlinger and M. W. Mitchell. Entangled photon apparatus for the
undergraduate laboratory. American Journal of Physics, 70 :898902, 2002.
[EN99] K.J. Engel and R. Nagel. One-parameter semigroups for linear evolution
equations, volume 194. Springer, 1999.
[EN06] K.J. Engel and R. Nagel. A short course on operator semigroups. Springer,
2006.
[F. 00] F. Faure. ExposÃ© sur la quantication gÃ©omÃ©trique. http ://www-

fourier.ujf-grenoble.fr/~faure, pages 112, 2000.
[Faua] F. Faure. Cours de Mécanique quantique pour Master M1 de physique.

http ://www-fourier.ujf-grenoble.fr/~faure/enseignement.
[Faub] F. Faure. Films d'animations d'ondes quantiques. http ://www-fourier.ujf-

grenoble.fr/~faure/enseignement/meca_q/animations.
[Fau10a] F. Faure. Cours de géométrie et topologie pour la physique pour Master M2

de physique. http ://www-fourier.ujf-grenoble.fr/~faure/enseignement, 2010.
[Fau10b] F. Faure. Cours de MathÃ©matiques pour la physique. Niveau Master 1.

http ://www-fourier.ujf-grenoble.fr/~faure/enseignement, 2010.
[Fau10c] F. Faure. Cours de Mécanique Analytique pour Licence L3 de physique.

http ://www-fourier.ujf-grenoble.fr/~faure/enseignement, 2010.
[Fey63] R. Feynman. Le cours de physique de Feynman, Mécanique quantique. 1963.
[GS90] V. Guillemin and S. Sternberg. Geometric asymptotics. Amer Mathematical

Society, 1990.
[Gut91] M. Gutzwiller. Chaos in classical and quantum mechanics. Springer-Verlag,

1991.
[Har02] S. Haroche. Cours de mÃ©canique quantique. http ://.., 2002.
[HB] Michel Hulin and Odile Betbeder. ThÃ©orie des groupes appliquÃ©e Ã la
physique.
+
[HTS 01] R. Hennino, N. Tregoures, M Shapiro, L. Margerin, M Campillo, B. van Tig-
gelen, and R. Weaver. Observation of equipartition of seismic waves. Phys.
Rev. Letters, 86 :34473450, 2001.
[Jac75] J.D. Jackson. Classical electrodynamics. Wiley, 1975.
[J.J85] J.J. Sakurai. Modern Quantum mechanics. Benjamin Cummings Publishing

Company, 1985.
[L.E90] L.E. Ballentine. Quantum Mechanics. 1990.
[LGA] J. Leys, E. Ghys, and A. Alvarez. Dimension. Videos.

http ://www.dimensions-math.org.
[lS02] Dossier Pour la Science. L'art du Secret. La cryptographie. Dossier hors

SÃ©rie, juillet/oct 2002.
BIBLIOGRAPHIE 395
[Mar02] A. Martinez. An Introduction to Semiclassical and Microlocal Analysis. Uni-

versitext. New York, NY : Springer, 2002.
[Mel95] Richard B. Melrose. Geometric scattering theory. Stanford Lectures.

http ://www-math.mit.edu/ rbm/book.html, 1995.
[Mes64] A. Messiah. mécanique quantique. 1964.
[Mil94] P. Milonni. The Quantum Vacuum, an introduction to quantum electrodyna-

mics. Academic Press, 1994.
[MR90] Philippe Martin and Francis Rothen. ProblÃ¨mes Ã N-corps et champs quan-
tiques. Presses Polytechniques et Universitaires Romandes (PPUR), 1990.
[Nak03] M. Nakahara. Geometry, topology and physics. Institute of Physics Publi-

shing, 2003.
[NGWH01] Nicolas Gisin, Grégoire Ribordy, Wolfgang Tittel, and Hugo Zbinden. Quan-
tum Cryptography. to appear in Reviews of Modern Physics, e-print : quant-
ph/0101098, 2001.
[Omn00] R. Omnes. Comprendre la mÃ©canique quantique. EDP Sciences, 2000.
[Paz83] A. Pazy. Semigroups of linear operators and applications to partial dierential

equations, volume 198. Springer New York, 1983.
[Per86] A. Perelomov. Generalized coherent states and their applications. Springer-

Verlag., 1986.
[PRK95] G. Papanicolaou, L. Ryzhik, and J. Keller. Stability of the p to s energy ratio

in the diusive regime. 1995.
[Roo06] Borwein Choi Rooney. The Riemann hypothesis. Springer, 2006.
[RS72] M. Reed and B. Simon. Mathematical methods in physics, vol I : Functional

Analysis. Academic press, New York, 1972.
[Sak67] J.J. Sakurai. Advanced quantum mechanics. 1967.
[Sch66] L. Schwartz. Theorie des distributions. 1966.
[Seg95] C. Segal. Lectures on Lie groups and Lie Algebras. 1995.
[SI00] Gustafson S. and Sigal I. Mathematical concepts of quantum mechanics. Sprin-

ger, 2000.
[Ste94] S. Sternberg. Group theory and physics. Cambridge University Press, 1994.
[Tay96a] M. Taylor. Partial dierential equations, Vol I. Springer, 1996.
[Tay96b] M. Taylor. Partial dierential equations, Vol II. Springer, 1996.
[WDG90] W. M. Zhang, D. H. Feng, and R. Gilmore. Coherent states : theory and some
applications. Rev. Mod. Phys., 62 :867, 1990.
[Woo92] N.M.J. Woodhouse. Geometric quantization. Clarendon Press, Oxford, 1992.
[YCB82] M. Dillard-Bleick Y. Choquet-Bruhat, C. Dewitt-Morette. Analysis, mani-

folds and physics. North-Holland, 1982.
396 BIBLIOGRAPHIE
[Zur91] W.H. Zureck. Phys. Today, 44 :36, 1991.
[Zwo12] M. Zworski. Semiclassical Analysis. Graduate Studies in Mathematics Series.

Amer Mathematical Society, 2012.
BIBLIOGRAPHIE 397
#script qui selectionne les chapitres pour faire des chiers pdf individuels.
pdftk cours.pdf cat 1-15 387-end output cours_intro.pdf
pdftk cours.pdf cat 19-85 output cours_chap1.pdf
pdftk cours.pdf cat 361-367 output cours_annexe.pdf
pdftk cours.pdf cat 369-386 output cours_solutions.pdf
pdftk cours.pdf cat 1-367 387-end output cours_complet_sans_solutions.pdf

Cours 5

Transféré par

Droits d'auteur :

Formats disponibles

Cours 5

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours 5

Transféré par

Droits d'auteur :

Formats disponibles

8.2.

THÉORIE DES PERTURBATIONS DÉPENDANT DU TEMPS 325

8.2.2.4 Expression de |Wba |2

On note ~ ba = −ehψb |~xˆψa i

8.2.3 Eet d'une onde incohérente

Nous supposons toujours l'atome soumis à l'inuence d'une onde électromagnétique à

Du fait de l'indépendance statistique des ondes planes composant un faisceau incohérent,

on utilise (8.2.8), (8.2.6) et (8.2.10), donnant :

où I est l'intensité de l'onde incidente à la fréquence ωab .

Démonstration. Il faut faire la moyenne sur les directions. @@

Exercice (TD) : @@ Bransden p508, la loi de somme de Thomas Reiche Kuhn.

8.2.3.3 Loi de décroissance

8.2.4 Transition vers les continuum ; La photo-ionisation ; Eet

Le taux de transition est

où Θ est l'angle entre ~k et l'axe de faisceau incident.

onde incidente Atome

8.3 Méthodes variationnelles

ref : Ballentine [L.E90]p.270.

Théorème 8.3.1. L'énergie fondamentale de l'opérateur Ĥ est donnée par

8.3.1.2 Méthode variationnelle

1. On choisit un famille d'états quantiques |φµ i appelée famille d'essai, dépendant

Figure 8.3.1  Illustre la méthode variationnelle : Pour estimer E0 , on recherche le mini-

8.3.1.3 Exemple : vibration anharmonique d'un atome (TD)

On reprend l'exemple (8.1.4) page 314, avec

Calcul de l'énergie moyenne :

C'est une intégrale Gaussienne (A.1.2), et donne :

Recherche du minimum de Eσ On observe que pour w → 0 ou w → +∞, alors

Résultats numériques et observations On pose Evariat. = E (wmin ). Comparons à la

8.3.1.4 (*) Formule du minimax

La formule suivante généralise la formule (8.3.1) qui donnait l'état fondamental.

alors le niveau n est donné par

Remarque 8.3.4. Cette formule est utilisée en chimie quantique.

9.1 Description d'un ensemble statistique d'états quan-

Probabilité et hasard en mécanique classique : A priori, en mécanique classique, un

Probabilité et hasard en mécanique quantique : En mécanique quantique, il y a

l'état ψ donnera le résultat aj avec la probabilité

kPj ψk2 hψ|Pj ψi

où Pj est le projecteur spectral sur l'espace propre de la valeur propre aj . De plus, si la

Armation 9.1.2. Postulat de la mesure : après la mesure, le système quantique

P P hψ|Pj ψi P hψ|aj Pj ψi hψ|Âψi

Dénition 9.1.3. Pour un ensemble statistique d'état ψi ∈ H, i = 1, 2, . . . avec des

Proposition 9.1.4. L'opérateur densité ρ̂ caractérise uniquement le système physique :

p (aj ) = Tr (ρ̂Pj ) (9.1.2)

où Pj est le projecteur spectral de Â (dénit plus haut). En particulier la valeur moyenne

Démonstration. On a (utilisant que Tr (|ψihϕ|) = hϕ|ψi)

Inversement, montrons que on peut reconstruire ρ̂ à partir de la connaissance des probabi-

Proposition 9.1.6. Un opérateur ρ̂ est un opérateur densité si et seulement si c'est un

et interpréter ρ̂ comme l'opérateur densité associé à l'ensemble statistique des états φj

X hv|ψi ihψi |vi X |hψi |vi|2

◦ ρ̂ est un état pur si et seulement si ses valeurs propres sont

Exercice 9.1.8. Montrer que ρ̂ est un état pur si et seulement si

9.1.2 Formulation du postulat de la mesure avec la matrice densité

Le postulat de la mesure énoncé en 9.1.1 peut être réformulé de la façon suivante :

Proposition 9.1.9. le postulat de la mesure : Si un système quantique est décrit

Démonstration. Avant la mesure, on a une distribution statistique d'états

On a utilisé que hPj ψi |Pj ψi i = hψi |Pj ψi i.

Tr (Pj ρ̂Pj ) = Tr ρ̂Pj2 = Tr (ρ̂Pj ) = p (aj )

d'après (9.1.2). C'est exactement ce que donne le postulat de la mesure.

Dans la section suivante nous évoquerons la théorie de la décohérence qui propose

8.2.3 Eet d'une onde incohérente

Nous supposons toujours l'atome soumis à l'inuence d'une onde électromagnétique à

8.2.4 Transition vers les continuum ; La photo-ionisation ; Eet

Figure 8.3.1 Illustre la méthode variationnelle : Pour estimer E0 , on recherche le mini-

Armation 9.1.2. Postulat de la mesure : après la mesure, le système quantique

Dénition 9.1.3. Pour un ensemble statistique d'état ψi ∈ H, i = 1, 2, . . . avec des

où Pj est le projecteur spectral de Â (dénit plus haut). En particulier la valeur moyenne

Proposition 9.1.9. le postulat de la mesure : Si un système quantique est décrit

Dans la section suivante nous évoquerons la théorie de la décohérence qui propose

ce qui signie qu'il y a très peu de suites standard.

Elles sont appelées suites standard (comme (9.1.11)).

ainsi on peut exprimer :

Donc on dénit les suites standard par

Démonstration. Notons pour simplier Z = |S|2 et p (Z) dZ la densité de probabilité. On

où h.i désigne la moyenne. (Ce dernier calcul s'appelle inégalité de Tchebychev). On

Dénition 9.1.20. Conformément à la dénition (9.1.12), l'entropie d'un ensemble

2. En mécanique classique (ou en théorie de l'information), une entropie nulle signie

p+ = 1, p− = 0, et la dénition ci-dessus donneS ρ̂, Ŝx = 0 = S (ρ̂)

Dénition 9.2.1. Si T̂ : Htot → Htot est un opérateur linéaire dans Htot = HA ⊗ HB

aussi noté plus simplement X