Integral de Chemin

Introduction à l’intégrale
de chemin et applications
Alain Comtet
Laboratoire de Physique Théorique et Modèles Statistiques

Université Paris-Sud, Orsay
&
Université Pierre et Marie Curie, Paris
27 mars 2012
Table des matières
1 Mécanique quantique 1
1.1 Définition du propagateur de Feynman . . . . . . . . . . . . . . . 1
1.2 Expression du propagateur de la particule libre . . . . . . . . . . 3
1.3 Construction générale du propagateur . . . . . . . . . . . . . . . 3
1.4 Propriétés du propagateur . . . . . . . . . . . . . . . . . . . . . . 6
1.5 Propagateur euclidien . . . . . . . . . . . . . . . . . . . . . . . . 9
2 Fonctions de corrélation quantiques, oscillateur harmonique 11

2.1 Expression d’éléments de matrice dans la représentation de Hei-
senberg . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Etude de l’amplitude vide-vide . . . . . . . . . . . . . . . . . . . 13
2.3 Calcul de Z (j) pour l’oscillateur harmonique . . . . . . . . . . . 16
2.4 Théorème de Wick . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3 Méthodes semi-classiques 23
3.1 Fonction d’onde et trajectoires classiques . . . . . . . . . . . . . 23
3.2 Calcul du propagateur semi-classique . . . . . . . . . . . . . . . . 24
3.3 Expression du propagateur semi-classique en dimension 1 . . . . 27
3.4 Propagateur dans l’espace E . . . . . . . . . . . . . . . . . . . . . 30
3.5 Effet tunnel-Etude du double puits . . . . . . . . . . . . . . . . . 34
3.6 Approximation du gaz dilué . . . . . . . . . . . . . . . . . . . . . 36
3.7 Calcul du déterminant : généralités . . . . . . . . . . . . . . . . . 36
3.8 Evaluation du déterminant dans l’approximation du gaz dilué . . 37
3.9 Méthode de Laplace et méthode de la phase stationnaire . . . . . 39
4 Applications à la physique des polymères 41

4.1 Modélisation d’un polymère comme marche aléatoire . . . . . . . 41
4.2 Polymère dans un potentiel . . . . . . . . . . . . . . . . . . . . . 44
4.3 Elasticité des biomolécules . . . . . . . . . . . . . . . . . . . . . . 46
4.3.1 Modèle gaussien discret . . . . . . . . . . . . . . . . . . . 46
4.3.2 Modèle du ver . . . . . . . . . . . . . . . . . . . . . . . . 49
5 Introduction au mouvement Brownien 51

5.1 Histoire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
5.1.1 Physique . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
5.1.2 Mathématiques . . . . . . . . . . . . . . . . . . . . . . . . 52
5.1.3 Théorie de la spéculation . . . . . . . . . . . . . . . . . . 53
5.2 Définition d’une marche aléatoire . . . . . . . . . . . . . . . . . . 53
iii
iv TABLE DES MATIÈRES
5.3 Marche aléatoire sur réseau . . . . . . . . . . . . . . . . . . . . . 54

5.3.1 Autre dérivation . . . . . . . . . . . . . . . . . . . . . . . 56
5.3.2 Programme de simulation . . . . . . . . . . . . . . . . . . 57
5.4 Limite continue . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
5.4.1 Théorème . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5.5 Notions de processus aléatoire . . . . . . . . . . . . . . . . . . . . 59
5.5.1 Probabilités conditionnelles . . . . . . . . . . . . . . . . . 60
5.5.2 Processus de Markov . . . . . . . . . . . . . . . . . . . . . 60
5.5.3 Equation de Chapman-Kolmogorov . . . . . . . . . . . . . 61
5.5.4 Processus aléatoire stationnaire . . . . . . . . . . . . . . . 61
5.5.5 Théorème de Doob . . . . . . . . . . . . . . . . . . . . . . 61
5.6 Exemples de processus stochastiques gaussiens . . . . . . . . . . 63
5.6.1 Processus de Wiener . . . . . . . . . . . . . . . . . . . . . 63
5.6.2 Processus d’Ornstein-Uhlenbeck . . . . . . . . . . . . . . 64
5.6.3 Bruit blanc . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5.6.4 Construction de W (t) comme processus . . . . . . . . . . 65
6 Intégrale de Wiener et formule de Feynman Kac 67

6.1 Equation de la diffusion, formalisme opératoriel . . . . . . . . . . 67
6.2 Mesure de Wiener . . . . . . . . . . . . . . . . . . . . . . . . . . 68
6.3 Mouvement brownien avec absorption-formule de Feynman-Kac . 70
6.4 Formule de Feynman-Kac . . . . . . . . . . . . . . . . . . . . . . 71
7 Équation de Langevin et Fokker-Planck 73

7.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
7.2 Représentation fonctionnelle des mesures gaussiennes . . . . . . . 75
7.3 Lemme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
7.4 Equation de Langevin . . . . . . . . . . . . . . . . . . . . . . . . 76
7.5 Equation de Fokker-Planck . . . . . . . . . . . . . . . . . . . . . 77
7.6 Solution d’équilibre . . . . . . . . . . . . . . . . . . . . . . . . . . 79
7.7 Relation avec l’équation de Schrödinger en temps imaginaire . . . 80
7.8 Propriété du hamiltonien H+ et supersymétrie . . . . . . . . . . 81
7.9 Application . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
7.9.1 Processus de Wiener . . . . . . . . . . . . . . . . . . . . . 82
7.9.2 Processus d’Ornstein-Uhlenbeck . . . . . . . . . . . . . . 83
7.9.3 Processus bistables . . . . . . . . . . . . . . . . . . . . . . 84
7.10 Représentation fonctionnelle de P (x t|y 0) . . . . . . . . . . . . . 85
7.11 Limite D → O, formule d’Arrhénius . . . . . . . . . . . . . . . . 86
8 Intégrale de chemin sur un espace multiplement connexe 89

8.1 Electron dans un anneau mésoscopique . . . . . . . . . . . . . . . 89
8.2 Notions d’homotopie . . . . . . . . . . . . . . . . . . . . . . . . . 90
8.3 Intégrale de chemin sur un espace multiplement connexe . . . . . 91
8.4 Effet Bohm-Aharonov . . . . . . . . . . . . . . . . . . . . . . . . 92
8.4.1 Généralités . . . . . . . . . . . . . . . . . . . . . . . . . . 92
8.4.2 Explication simplifiée . . . . . . . . . . . . . . . . . . . . 93
8.4.3 Ligne de flux magnétique . . . . . . . . . . . . . . . . . . 94
TABLE DES MATIÈRES v
9 Statistiques quantiques 97
9.1 Rappel : problème à une particule . . . . . . . . . . . . . . . . . 97
9.2 Problème à n particules . . . . . . . . . . . . . . . . . . . . . . . 98
9.3 Calcul de la fonction de partition bosonique . . . . . . . . . . . . 100
9.4 Discussion qualitative des trajectoires . . . . . . . . . . . . . . . 101
9.5 Bornes sur les fonctions de partition . . . . . . . . . . . . . . . . 104
vi TABLE DES MATIÈRES
Chapitre 1
Formulation de la
mécanique quantique par
l’intégrale de chemin
Sommaire
1.1 Définition du propagateur de Feynman . . . . . . 1
1.2 Expression du propagateur de la particule libre . 3
1.3 Construction générale du propagateur . . . . . . . 3
1.4 Propriétés du propagateur . . . . . . . . . . . . . . 6
1.5 Propagateur euclidien . . . . . . . . . . . . . . . . . 9
1.1 Définition du propagateur de Feynman

Dans la représentation de Schrödinger, l’évolution dans le temps des états est
gouvernée par l’équation :
d
ih | ψ (t) >= H
b | ψ (t) >
dt
qui s’intègre dans le cas particulier d’un hamiltonien indépendant du temps :
i b 0
| ψ (t0 ) >= exp − H (t − t) | ψ (t) >
~
Considérons par exemple le cas d’une particule dans un potentiel. Le hamilto-
nien s’écrit
Hb = 1 pb2 + V (b q)
2m
Dans la représentation {b
q } l’équation d’ évolution devient
i b 0
< q~0 | ψ (t0 ) >=< q~0 | exp − H (t − t) | ψ (t) >
~
1
2 CHAPITRE 1. MÉCANIQUE QUANTIQUE
Z
i b 0
= dd ~q < q~0 | exp − H (t − t) | ~q >< ~q | ψ (t) >
~
soit encore
Z
ψ q~0 , t0 = K q~0 t0 |~qt ψ (~q, t) dd ~q

• Le propagateur K q~0 t0 |~qt ≡< q~0 | exp − i~H (t0 − t) |~q > permet d’évaluer
b
l’amplitude de transition entre les deux états. Prenons un état initial lo-
calisé en q~o
ψ (~q, to ) = δ (~q − q~o )

ψ q~0 , t0 = K q~0 t0 | q~o to
donne l’amplitude de probabilité de trouver la particule en q~0 à l’instant

t’ sachant qu’elle était en q~o à l’instant to .
• Le propagateur est la solution de l’équation de Schrödinger

dK ~2
i~ 0
= [− ∆q~0 + V q~0 ]K
dt 2m
qui satisfait la condition initiale

lim K ~0 t0 |~qt =< q~0 |~q >= δ d q~0 − ~q
q
0 t →t
• On observe que dans la représentation de Heisenberg cette amplitude peut

encore s’écrire < q~0 t0 |~qt > où les kets |~qt > et |q~0 t0 > sont états propres de
l’opérateur de position q̂ (τ ) aux temps τ = t et τ = t0
Preuve : Dans la représentation de Heisenberg les vecteurs d’état sont

indépendants du temps, en revanche les observables O
d H évoluent selon
i i
OH = exp Ĥ (t − to ) O exp − Ĥ (t − to )
d b
~ ~
où O
b est l’observable correspondante de la représentation de Schrödinger.
Les vecteurs d’états entre les deux représentations sont reliés par
i
|ψH >= exp Ĥ (t − to ) | ψS >
~
Donc | ~q t >= exp ~i Ĥ (t − to ) | ~q > est état propre de l’opérateur de
position q̂H (t).
Par conséquent :
i
< q~0 t0 | =< q~0 | exp − Ĥ (t0 − to )
~
et
i
< q~0 t0 |~q t >=< q~0 | exp − Ĥ (t0 − t) |~q >
~
1.2. EXPRESSION DU PROPAGATEUR DE LA PARTICULE LIBRE 3
1.2 Expression du propagateur de la particule

libre
p >< p~|dd p~ = 1 ainsi que l’identité
R
Insérons deux fois la relation de fermeture |~
d2
1 i
< ~q|~
p >= exp ~qp~
2π~ ~
Le propagateur correspondant s’écrit :

Z ip̂2 0
Ko q t |~qt = dd p~dd p~0 < q~0 |p~0 >< p~0 | exp −
~0 0
(t − t) |~
p >< p~|~q >
2m~
d Z
p2 0

1 i ~0 ~0
d ~0
i~
= d
d p~d p exp q p − ~qp~ − (t − t) δ p~0 − p~
2π~ ~ 2m~
d Z
p2 0

1 d i~
p ~0 i~
= d p~ exp q − ~q − (t − t)
2π~ ~ 2m~
L’exposant peut encore s’écrire :
 2 2 
 m ~0 − q~0
q m2 ~0
q − ~
q 
i 0
 
− (t − t) p~ −  − 2
2m~ 
 t0 − t (t0 − t) 

En intégrant sur p~ on obtient :

2

1
d
2πm~
d2 im q~0 − ~q
Ko q~0 t0 |~qt = exp
2π~ i (t0 − t) 2~ (t0 − t)
2

m
d2 im q~0 − ~q
Ko q~0 t0 |~qt = exp
2πi~ (t0 − t) 2~ (t0 − t)
iS
Nous observons que l’argument de l’exponentielle est de la forme ~ où S est
l’action calculée le long de la trajectoire classique :
τ − t ~0
~q (τ ) = ~q + 0 q − ~q
t −t
1.3 Construction générale du propagateur

Elle s’appuie sur la formule de Trotter :
h t t
in
lim e n A e n B = et(A+B)
n→∞
dont on sait donner une formulation mathématiquement satisfaisante (A et B

étant des matrices ou bien des opérateurs bornés sur un espace de Banach).
Posons A = ~1 Ĥo B = ~1 V̂ et t = −iT . On obtient

−iT H −iT V
n
< q~0 |e− ~ (Ho +V ) |~q >= lim < q~0 | e n~ e n~
iT co b
|~q >
b b
n→∞
Nous allons insérer dans cette expression n − 1 états intermédiaires. Il nous faut
donc calculer
THco T Vb
< q~i | exp −i exp −i |~qi+1 >=
n ~ n ~
Z
iT H
bo iT Vb
d~q < ~qi | exp − |~q >< ~q| exp − |~qi+1 >=
n ~ n ~
iT Hco iT
< ~qi | exp − |~qi+1 > exp − V (~qi+1 )
n ~ n~
En remplaçant le 1er terme par le propagateur libre on peut réécrire cet élément
de matrice sous la forme :
! d2
2
m im (~qi − ~qi+1 ) iT V (qi+1 )
exp exp −
2πi~ Tn 2~ Tn n ~
Nous pouvons ainsi exprimer le propagateur sous forme d’ une intégrale multiple
< q~0 |e−i ~ (Ho +V ) |~q >=

T b b
Z
iT H
bo iT Vb
lim < q~0 | exp − exp − |~q1 >< ~q1 | . . . . . .
n→∞ n ~ n ~
iT H
bo iT Vb
. . . . . . < ~qn−1 | exp − exp − |~q > d~q1 . . . d~qn−1
n ~ n ~
Il est commode de rebaptiser les variables q~0 = ~qo et ~q = ~qn . On obtient après
passage à la limite n → ∞
iT b
< q~0 | exp − H|~ q >=
~
! nd
2 Z
m i
lim exp S (~qo , ~q1 . . . ~qn ) d~q1 . . . d~qn−1
n→∞ 2πi~ Tn ~
où
n−1
X m 2 T
S (~qo , . . . ~qn ) = T
(~qi − ~qi+1 ) − V (~qi+1 )
i=0
2n n
Dans la limite n → ∞ , cette somme discrète constitue une discrétisation plau-

sible de l’intégrale
Z t0 m 2
~ 0
S q , ~q = ˙
~q (τ ) − V (~q (τ )) dτ (1.1)
t 2
1.3. CONSTRUCTION GÉNÉRALE DU PROPAGATEUR 5
En effet, posons
~qi − ~qi+1
~q˙i+1 = T
(1.2)
n
Il vient :
n−1
Xh m ˙2 iT
S (~qo , ~q1 . . . ~qn ) = ~qi+1 − V (~qi+1 )
i=0
2 n
Cette écriture est purement formelle car rien ne prouve que cette somme ait
effectivement une limite au sens de Riemann. Il faudrait pour celà que les che-
mins sur lesquels on intègre soient dérivables. On s’attend en fait à ce que les
chemins qui donnent une contribution dominante soient ceux pour lesquels :
2
im (~qi − ~qi+1 )
exp T
∼1
2~ n
soit
2 ~ T
(~qi − ~qi+1 ) ∼ 2
mn
2
d’où (∆~q) ∼ ∆t . L’intégrale précédente semble donc concentrée sur des tra-
jectoires non différentiables.
Oubliant pour l’instant ces subtilités, nous définissons formellement la mesure
de Feynman en posant
!n d2
m
D~q (τ ) = lim d~q1 . . . d~qn−1
n→∞ 2πi~ Tn
Nous obtenons ainsi l’expression fonctionnelle du propagateur :
Z 0
Z i t m ˙2
K q~0 t0 |~qt = D~q (τ ) exp ( ~q (τ ) − V (~q (τ )) dτ
~ t 2

Le propagateur K q~0 t0 |~qt fait donc intervenir une somme infinie sur tous
les chemins reliant les points ~q et q~0 . Chaque terme de cette somme est donné
par une phase proportionnelle à l’action classique.
t’
~q q~0
1.4 Propriétés du propagateur

1. Nous avons adopté la notation K q~0 t0 |~qt pour faire le lien avec la représentation
des éléments de matrice < q~0 t0 |~qt > dans le point de vue de Heisenberg.
L’invariance par translation dans le temps implique que K ne dépend que
de la différence t0 − t

K q~0 , t0 + τ |~q, t + τ = K q~0 t0 |~qt
R
2. L’insertion de l’identité |~q1 >< q~1 |dq~1 = 1 dans l’élément de matrice
conduit aux relations suivantes
i b 0
K q~0 t0 |~qt =< q~0 | exp − H (t − t) |~q >
~
Z
i b 0 iH
b
= < q~0 | exp − H (t − τ ) |q~1 >< q~1 | exp − (τ − t) |~q > dq~1
~ ~
Z
= K q~0 t0 |q~1 τ K (q~1 τ |~qt) dq~1
Z
K q~0 t0 |~qt = dq~1 K q~0 t0 |q~1 τ K (q~1 τ |~qt) (1.3)
1.4. PROPRIÉTÉS DU PROPAGATEUR 7
t’ t
t τ
En termes physiques cette équation exprime le principe de Huygens de l’
optique ondulatoire.
Dans la théorie de processus stochastiques l’ équation (1.3) est l’ équation
de Chapmann Kolmogorov.
(Pour donner un sens plus précis à cette correspondance, il faudrait pou-
voir interpréter K comme une probabilité de transition et non pas comme
une amplitude de probabilité)
3. Dans ce qui précède, nous avons supposé que l’état initial ainsi que l’état
final sont états propres de l’opérateur de position.
Pour des états initiaux et finaux quelconques décrits respectivement par
ψi (q) et ψf (q) la probabilité de transition s’écrit :
Z
ψf∗ q~0 , t0 K q~0 t0 |~qt ψi (~q, t) dn ~qdn q~0
4. La connaissance du propagateur nous permet d’accéder au spectre du ha-

miltonien. Pour un spectre discret il vient :
0
K q~0 t0 |~qt =< q~0 |e− ~ (t −t) |~q >
iH
c
−iH
(t0 −t) |n >< n|~q >
X
< q~0 |n >< n|e
c
= ~
n
X iEn (t0 − t)
= ψn q~0 ψn∗ (~q) exp −
n
~
Posons q~0 = ~q et intégrons sur ~q on obtient
iEn (t0 − t)
Z X
φ (t0 − t) = d~q K (~qt0 |~qt) = exp −
n
~
La transformée de Fourier de φ (t)
i ∞
Z
it
X 1
φ (t) e ~ (E+i) dt =
~ 0 n
E n − E − i
est la trace de la résolvante

1
G (z) =
H −z
On a X 1
T race G (E + i) =
n
En − E − i
Les pôles de G (z) donnent les états liés du système. On en déduit la
densité d’état
1 X
ρ (E) = Im T race G (E + i) = δ (En − E)
π n
5. Considérations topologiques
Reprenons l’expression du propagateur

Z 0
~

0 0
Z i t
K q t |~qt = D~q (τ ) exp Ldτ
~ t
Si l’espace de configuration {~q (t)} noté (C) est multiplement connexe,

nous verrons qu’il sera utile de distinguer les chemins appartenant à des
classes d’homotopie distinctes.
Par exemple si (C) = R2 / {O} les chemins fermés sont étiquetés par un
entier n ∈ Z qui représente le nombre de tours effectués par la trajectoire
autour du point O.
n=0
~q = q~0
r = r~0
0
n=1
0
n=2
Le propagateur sera donné par l’expression suivante :
Z 0
i t
X Z
K q~0 t0 |~qt = χ (g) D~q (τ ) exp Ldτ
~ t
gπ1 (c)
Pour satisfaire la loi de composition des propagateurs

Z
K q t |~qt = K q~0 t0 |q~00 t00 K q~00 t00 |~qt dq~00
~0 0
1.5. PROPAGATEUR EUCLIDIEN 9
les facteurs de phase doivent satisfaire la relation χ (g) χ (h) = χ (g ◦ h)

Cette relation exprime que les χ forment une représentation unitaire uni-
dimensionnelle du groupe d’homotopie π1 (C)
exemple :
(C) = S1 , R2 / {O}, χ = einθ
1.5 Propagateur euclidien

La physique statistique des systèmes quantiques à l’équilibre fait intervenir
la fonction de partition
Z(β) = T race e−βH

Z
= d~q < ~q | e−βH | ~q >
−βH
ainsi que l’opérateur densité ρ = e Z dont les éléments de matrice
s’écrivent
1 X
ρ(~q, q~0 ) = ψn (~q0 ) ψn∗ (~q) exp −βEn
Z n
Montrons comment exprimer ces quantités sous forme d’intégrale fonction-

nelle. Un calcul en tout point similaire au calcul précédent permet d’écrire
le propagateur euclidien en dimension d’espace d = 1.
1 t
Z Z
tH
< q 0 | exp − | q >= DE q (τ ) exp − H(q, q̇)dτ
~ q(t)=q 0 ~ 0
q(0)=q
où
mq̇ 2
H(q, q̇) = + V (q)
2
est le hamiltonien classique et DE q(τ ) une certaine mesure ”euclidienne”
sur les chemins. En posant t = β~ il vient
1 β~
Z Z
< q 0 | exp −βH | q >= DE q (τ ) exp − H(q, q̇)dτ
q(t)=q 0 ~ 0
q(0)=q
Par conséquent la fonction de partition quantique

Z
Z(β) = dq < q | e−βH | q >
s’obtient en intégrant sur toutes les trajectoires fermées parcourues dans

~
le temps β~ = kT . Dans ce cadre on peut maintenant donner un sens
mathématique précis aux différents objets que nous avons rencontrés jusque
là. Posons pour simplifier ~ = m = 1. Un chemin quelconque q(τ ) issu de
l’origine peut être représenté sous forme d’une série de Fourier
∞
!
p τ √ X sin(nπτ /β)
q(τ ) = β ξ0 + 2 ξn
β n=1
nπ
En utilisant les propriétés d’orthogonalité des fonctions fn (τ ) = cos nπτ /β

on peut réecrire le terme cinétique sous la forme
Z β ∞
1 1X 2
q̇ 2 (τ )dτ = ξ
2 0 2 n=0 n
C’est l’identité de Parseval. Convenons de définir la mesure euclidienne

DE q(τ ) par la formule
β ∞
ξ2
Z
1 dξ
√ n exp − n
Y
DE q(τ ) exp − q̇ 2 (τ )dτ =
2 0 n=0
2π 2
Le propagateur d’une particule libre (H0 = 21 p2 ) peut ainsi s’écrire

∞
ξ2
Z Y
dξ
<q|e −βH0
| 0 >= √ n exp − n δ(q − q(β))
n=0
2π 2
En séparant les intégrales sur ξ0 et ξn on obtient
1 q2
< q | e−βH0 | 0 >= √ exp −
2πβ 2β
qui est le résultat attendu pour le propagateur libre euclidien. Cette obser-
vation nous invite à construire une mesure de probabilité sur les chemins.
Pour celà nous introduisons un ensemble infini de variables aléatoires gaus-
siennes indépendantes ξn telles que E(ξn ) = 0, E(ξn2 ) = 1. La trajectoire
représentée par la série de Fourier aléatoire q(τ ) est la trajectoire d’un
processus stochastique tel que
∞
ττ0 X 1
E(q(τ )q(τ 0 )) = + 2β 2 π2
sin(nπτ /β) sin(nπτ 0 /β)
β n=1
n
où E désigne l’espérance mathématique par rapport à la mesure produit

définie plus haut. Le calcul donne
E(q(τ )q(τ 0 )) = min(τ, τ 0 )
Cette fonction de corrélation définit le processus de Wiener, processus

à trajectoires continues mais non dérivables.La formule exprimant q(τ )
comme une série de Fourier aléatoire est due à Paley et Wiener.
Chapitre 2
Fonctions de corrélation
quantiques, oscillateur
harmonique
Sommaire
2.1 Expression d’éléments de matrice dans la représentation
de Heisenberg . . . . . . . . . . . . . . . . . . . . . 11
2.2 Etude de l’amplitude vide-vide . . . . . . . . . . . 13
2.3 Calcul de Z (j) pour l’oscillateur harmonique . . . 16
2.4 Théorème de Wick . . . . . . . . . . . . . . . . . . 18
2.1 Expression d’éléments de matrice dans la

représentation de Heisenberg
On se propose d’évaluer l’élément de matrice d’un opérateur diagonal, soit par
exemple < q 0 t0 |q̂ (t1 ) |qt >. prenons pour simplifier une description dans laquelle
t = 0 est l’instant de référence, donc :
i i
q̂ (t1 ) = e ~ Ht1 q̂e− ~ Ht1
b b
i
|qt >= e ~ Ht |q >
b
d’où
< q 0 t0 |q̂ (t1 ) >=< q 0 |e−

iH
c
~ (t0 −t1 ) q̂e i~H
c
(t0 −t1 ) |q >
Insérant des états intermédiaires et faisant apparaı̂tre les propagateurs, on peut

réécrire cette expression
Z
K (q 0 t0 |q1 t1 ) q1 K (q1 t1 |qt) dq1
En terme d’intégrale fonctionnelle K (q 0 t0 |q1 t1 ) est donnée par une intégrale sur
tous les chemins tels que q (t0 ) = q 0 ; q (t1 ) = q1 de même K (q1 t1 |qt) est donné
11
12CHAPITRE 2. FONCTIONS DE CORRÉLATION QUANTIQUES, OSCILLATEUR HARMONIQU
par une intégrale sur les chemins tels que q (t1 ) = q1 ; q (t) = q. Puisque le point
q1 est arbitraire, réunissant ces 2 chemins, on intégrera sur tous les chemins
q (τ ) tels que q (t) = q ; q (t0 ) = q 0
q1 définit un point intermédiaire d’une nouvelle discrétion
t0
t1
q
q q1 q 0
Si nous réexprimons les deux propagateurs en terme d’une intégrale fonc-

tionnelle, l’action totale étant additive il vient :
Z Z t0
0 0 i
< q t |q̂ (t1 ) |qt >= Dq (τ ) q (t1 ) exp L (q, q, τ ) dτ
~ t
Essayons maintenant d’évaluer :
< q 0 t0 |T q̂ (t2 ) q̂ (t1 ) |qt >
Il convient cette fois ci de découper le chemin en trois chemins élémentaires

ordonnés dans le temps.
2.2. ETUDE DE L’AMPLITUDE VIDE-VIDE 13
t0
t2
t1
q q1 q2 q0 q
Pour t2 > t1 il vient
Z
iS
< q 0 t0 |q̂ (t2 ) q̂ (t1 ) |qt >= Dq (τ ) q (t2 ) q (t1 ) exp
~
Si t2 < t1 on rangera les opérateurs dans l’ordre chronologique
Z
0 0 iS
< q t |q̂ (t1 ) q̂ (t2 ) |qt >= Dq (τ ) q (t1 ) q (t2 ) exp
~
D’où la formule générale :
Z Z t0
0 0 i
< q t |T q̂ (t2 ) q̂ (t1 ) |qt >= Dq (τ ) q (t2 ) q (t1 ) exp L (q, q, τ ) dτ
~ t
2.2 Etude de l’amplitude vide-vide

On s’intéresse plus particulièrement aux valeurs moyennes dans le vide (c’est à
dire dans l’état fondamental) de T produits d’opérateurs. Ce sont les fonctions de Green.
Montrons comment accéder à < Ω|T q̂ (t1 ) q̂ (t2 ) |Ω > en partant de
< q 0 t0 |T q̂ (t1 ) q̂ (t2 ) |qt >
Supposons que t1 > t2 . Insérant des états intermédiaires il vient :
< q 0 t0 |T q̂ (t1 ) q̂ (t2 ) |qt >=

X
< q 0 t0 |n >< n|q̂ (t1 ) q̂ (t2 ) |m >< m|qt >
n,m
les états intermédiaires sont états propres du hamiltonien Ĥ tels que Ĥ|n >=
En |n >.
L’expression ci dessous peut encore s’écrire :
0
(q) e ~ (Em t−En t ) < n|q̂ (t1 ) q̂ (t2 ) |m >
X i
ψn (q 0 ) ψm
∗
n,m
τ →−∞
Posons t = −iτ et t0 = −iτ 0 et prenons la limite τ 0 →+∞ Il vient :
0
(q) e ~ (Em τ −En τ ) < n|q̂ (t1 ) q̂ (t2 ) |m >
X 1
ψn (q 0 ) ψm
∗
n,m
E0
→ ψ0∗ (q 0 ) ψ0 (q) e− ~ (τ 0 −τ ) < Ω|q̂ (t ) q̂ (t ) |Ω >
1 2
en effet seul le fondamental contribue dans cette limite.

De même
< q 0 t0 |qt >→ ψ0∗ (q 0 ) ψ0 (q) e−

Eo
~ (τ 0 −τ )
D’où l’expression des fonctions de Green
< q 0 t0 |T q̂ (t1 ) q̂ (t2 ) |qt >

< Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >= lim
0
t →−i∞
t→+i∞
< q 0 t0 |qt >
Fonctionnelle génératrice
Couplons le système à une force extérieure j (τ ) dépendant du temps
L → Lf = L + ~q (τ ) j (τ )
Par conséquent l’amplitude prend la forme

Z Z t0
0 0 (j) 1
< q t |qt > = Dq (τ ) exp i L + q (τ ) j (τ ) dτ
t ~
En dérivant par rapport à j (τ )
1 δ
< q 0 t0 |qt >(j) |j=0 =< q 0 t0 |q̂ (t1 ) |qt >
i δj (t1 )
pour t < t1 < t0
De même en dérivant deux fois :

2
1 δ δ
< q 0 t0 |qt >(j) |j=0 =< q 0 t0 |T q̂ (t1 ) q̂ (t2 ) |qt >
i δj (t1 ) δj (t2 )
2.2. ETUDE DE L’AMPLITUDE VIDE-VIDE 15
La formule précédente donnant l’amplitude < Ω|T q̂ (t1 ) q̂ (t2 ) |Ω > prend la
forme :
1 δ
< Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >= Z (j)
Z (0) δj (t1 ) δj (t2 )
où Z (j) = lim t0 →−i∞ < q 0 t0 |qt >(j) .

t→+i∞
Cette écriture, en grande partie symbolique nous invite à effectuer une conti-
nuation analytique au niveau de l’intégrale fonctionnelle elle même (passage en
temps imaginaire).
Considérons
iĤ
− (t2 −t1 )
< q 0 t2 |qt1 >=< q 0 |e ~ |q >
Z Z t2
i
= Dq (τ ) exp L (q, q̇) dτ
~ t1
posons t2 = iτ2 t1 = −iτ1 τ = −iτ 0
i t2 e
Z Z
H
q 0 |e− ~ (τ2 −τ1 ) |q >= Dq (τ 0 ) exp L (q, q̇) dτ 0
c
~ t1
avec Le (q, q̇) = − m 2

2 q̇ − V (q) = −H (q, q̇) le passage en temps imaginaire trans-
forme le lagrangien en le hamiltonien .
1b
1 τ =τ2
Z Z
−
H(τ2 −τ1 )
0
< q |e ~ |q >= Dq (τ ) exp − H (q, q̇) dτ
~ τ =τ1
la limite τ2 → ∞ τ1 → −∞ fait apparaı̂tre.
1 +∞
Z Z
Dq (τ ) exp − H (q, q̇) dτ
~ −∞
La formule précédente exprime que la fonction à deux points < Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >
est la continuation analytique (t1 = −tτ1 , t2 = −iτ2 ) de
R +∞
Dq (τ ) q (τ1 ) q (τ2 ) exp − ~1 −∞ H (q, q̇) dτ
R
R +∞
Dq (τ ) exp − ~1 −∞ H (q, q̇) dτ
R
On peut donner une forme plus compacte à cette dernière expression en intro-
duisant la fonctionnelle génératrice euclidienne :
Z Z +∞ Z +∞
1
Z (j) = Dq (τ ) exp − H (q, q̇) − j (τ ) q (τ ) dτ
−∞ ~ −∞
qui décrit le couplage du système à une source classique dépendant du temps.

La dérivation fonctionnelle par rapport à la source fait apparaı̂tre les valeurs
moyennes recherchées.
δ2 Z (j)
< Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >= continuation analytique de |j=0
δj (τ1 ) δj (τ2 ) Z (0)
Nous allons illustrer cette approche en traitant en détail le problème de l’os-

cillateur harmonique.
R
Remarque : le signe du terme de source est parfois + j (τ ) q (τ ) dτ
2.3 Calcul de Z (j) pour l’oscillateur harmonique

Considérons le hamiltonien :
mq̇ 2 mω 2 q 2
H (q, q̇) = +
2 2
On considère la fonctionnelle génératrice (~ = 1)
Z Z
Z (j) = Dq (τ ) exp − dτ [H (q, q̇) + j (τ ) q (τ̇ )]
Après intégration par partie, l’exponentielle s’écrit :
m d2 mω 2
Z
exp − dτ q (τ ) − + q (τ ) + j (τ ) q (τ )
2 dτ 2 2
Une discrétion de cette intégrale nous conduirait à une expression de la forme.
1
E=
q̃1 Aij qj + j̃i qi
2
où A est une matrice symétrique. Pour diagonaliser cette forme quadratique on
pose
q = Q − A−1 j
e−1 = Q
e − j̃ A
q̃ = Q e − j̃A−1
D’où
1e 1
E= QAQ − j̃A−1 j
2 2
Par extension le changement de variable fonctionnel
q = Q − A−1 j
donnera
Z
1
Z (j) = exp dτ dτ 0 j̃ (τ ) A−1 (τ 0 , τ ) j (τ )
2
Z Z hm m 2 2i
DQ (τ ) exp − dτ Q̇2 + ω Q
2 2
2.3. CALCUL DE Z (J) POUR L’OSCILLATEUR HARMONIQUE 17
Soit encore
Z
1
Z (j) = Z (0) exp dτ dτ 0 j̃ (τ 0 ) A−1 (τ 0 , τ ) j (τ )
2
Reste à calculer A−1 (τ 0 , τ )

d2
Par définition AA−1 = 1 soit −m dτ 2 + mω
2
A−1 (τ 0 , τ ) = δ (τ 0 − τ )
En passant aux transformées de Fourier il vient

0
+∞
eik(τ −τ )
Z
−1 0 1
A (τ , τ ) = dk 2
2πm −∞ k + ω2
1 −ω|τ −τ 0 |
= e
2mω
Le calcul précédent nous enseigne que

< Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >= continuation analytique
δ2 Z (j)
|j=0 = A−1 (τ1 − τ2 )
δj (τ1 ) δj (τ2 ) Z (0)
Vérifions que A−1 (τ1 , τ2 ) = 2mω

1
e−ω|τ1 −τ2 | est bien la continuation en temps
imaginaire du propagateur de l’oscillateur harmonique.
Il nous donc évaluer explicitement
< Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >
Pour l’oscillateur harmonique le hamiltonien s’écrit :

2 2
e = p̂ + mω q̂ 2
H
2m 2
Pour t1 > t2 il vient
< Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >=
< Ω|eiHt1 q̂e−iHt1 eiHt2 q̂e−iHt2 |Ω >=

b b b b
X
eiEo (t1 −t2 ) | < Ω|q̂|n > |2 e−iEn (t1 −t2 )
n
Pour l’oscillateur harmonique le seul élément de matrice non nul est celui qui
connecte l’état fondamental |Ω > au premier état excité.
Les fonctions d’onde correspondantes s’écrivent
mω 14 mω 2
< q|Ω >= ψo (q) = exp −q
π 2
√
ψ1 (q) = 2mωqψo
R −∞
D’où < Ω|q̂|n >= ψo (q) q ψ1 (q) dq = √ 1
−∞ | {z } 2mω
√ 1 ψ
2mω 1
Par conséquent
1
< Ω|T q̂ (t1 ) q̂ (t2 ) |Ω >= 2mw e−iw(t1 −t2 ) pour t1 > t2
De façon générale on obtient
1 n −iω(t1 −t2 ) o
< Ω | T q̂(t1 )q̂(t2 ) | Ω > = e θ(t1 − t2 ) + e−iω(t2 −t1 ) θ(t2 − t1 )
2mω
1 −iω|t1 −t2 |
= e
2mω
Si dans cette formule on pose t = −iτ on obtient 2mω 1

e−ω|τ1 −τ2 | qui est bien
−1
l’expression obtenue pour A (τ1 , τ2 ) par la méthode fonctionnelle.
Remarques
1. le calcul en temps imaginaire fournit des expressions mathématiquement
d2 2
saines : opérateurs elliptiques du type − dt2 + ω pour l’oscillateur har-
monique. La continuation analytique fournit ensuite automatiquement la
prescription de Feynmann.
2. les fonctions de Green < Ω | T q̂(t1 )q̂(t2 ) | Ω > contiennent de l’informa-
tion sur leP spectre du hamiltonien. En effet < Ω | T q̂(t1 )q̂(t2 ) | Ω >=
eiE0 (t1 −t2 ) n < Ω | q̂ | n >| e−iEn (t1 −t2 ) pour t1 > t2 .
Le comportement pour t → ∞ des fonctions de Green euclidiennes contrôle
l’énergie du premier état excité (relativement au fondamental).
1 −ω|τ1 −τ2 |
< Ω | T q̂(t1 )q̂(t2 ) | Ω >−→ e
2mω
pour | τ1 − τ2 |→ ∞ avec une décroissance exponentielle donnée par
3 1
E1 − E0 = ω− ω=ω
2 2
2.4 Théorème de Wick

Supposons pour simplifier que A est une matrice n×n d’élément Aij . Le courant
j est un vecteur à n composantes j1 . . . jn . La fonctionnelle génératrice s’écrit
Z Y
1
Z(j1 , . . . jn ) = dqi exp − q̃ A q − j̃q
i
2
1 1
'√ exp j̃A−1 j
det A 2
Considérons une fonction de Green à p points

Z Y
1 1
< qi1 . . . qip >= dqi qi1 . . . qip exp − q̃Aq
Z(0) i
2
2.4. THÉORÈME DE WICK 19
Par parité les fonctions de Green d’ordre impair sont nulles. Montrons que celles
d’ordre pair s’expriment en terme de la fonction à deux points
Z Y
1 1
< qi1 qi2 > = dqi qi1 qi2 exp − q̃Aq
Z(0) i
2
1 ∂ 2 Z(j)
= |j=0 = (A−1 )i1 i2
Z(0) ∂ji1 ∂ji2
La fonction à quatre points s’écrit
∂2

1 1
< qi1 qi2 qi3 qi4 > = exp j̃A−1 j
Z(0) ∂ji1 ∂ji2 ∂ji3 ∂ji4 2 j=o
= (A−1 )i1i2 (A−1 )i3i4 + (A−1 )i1i3 (A−1 )i2i4 + (A−1 )i1i4 (A−1 )i2i3
=< qi1 qi2 >< qi3 qi4 > + < qi1 qi3 >< qi2 qi4 > + < qi1 qi4 >< qi2 qi3 >
On retrouve ainsi le théorème de Wick comme une propriété élémentaire des

intégrales gaussiennes. Elle se généralise aux intégrales fonctionnelles.
Exemple : oscillateur
Z
1
Z(j) = Z(0) exp dτ1 dτ2 j(τ1 )A−1 (τ1 τ2 )j(τ2 )
2
1 δ 2 Z(j)
< Ω | T q(τ1 )q(τ2 ) | Ω > =
Z(0) δj(τ1 )δj(τ2 )
= A−1 (τ1 , τ2 )
1 −ω|τ1 −τ2 |
= e
2mω
Les fonctions de Green à 4 points sont données par le théorème de Wick
< Ω | T q(τ1 )q(τ2 )q(τ3 )q(τ4 ) | Ω >=
< Ω | T q(τ1 )q(τ2 ) | Ω >< Ω | T q(τ3 )q(τ4 ) | Ω >

+ < Ω | T q(τ1 )q(τ3 ) | Ω >< Ω | T q(τ2 )q(τ4 ) | Ω >
+ < Ω | T q(τ1 )q(τ4 ) | Ω >< Ω | q(τ2 )q(τ3 ) | Ω >
Application : calcul perturbatif de fonctionnelles non gaussiennnes
n ∞ Z n
(−1)k k Y
Z Y X
1 1
dqi exp − q̃Aq + λV (q) = λ dqi V (q)k exp − q̃Aq
i=1
2 k! i=1
2
k=1
chaque terme se calcule grâce au théorème de Wick.
Remarque
Autre dérivation de Z(j)
Introduisons les fonctions propres (réelles) de l’opérateur

d2
A = −m + mω 2
dt2
Aqn (t) = n qn (t)
On a une base de fonctions propres orthogonales
Z +∞
qn (t)qm (t)dt = δnm
−∞
Développons sur cette base une fonction q(t) quelconque ainsi que j(t)
X
q(t) = cn qn (t)
X
j(t) = jn qn (t)
Nous voulons évaluer
Z
mω 2
R +∞ m 2
q 2 − j(t)q(t)
R
Z(j) = Dq(t)e −∞ 2 q̇ + 2
On peut montrer que

+∞
Y
Dq(t) = N dcn
n=−∞
Z Y
Z(j) = N dcn e−ξ
Z
1 1X 2
ξ= q(t)Aq(t) + j(t)q(t) = cn n + jn n
2 2
X n 2jn cn

= c2n +
2 n
" 2 #
X n jn jn2
= cn + − 2
2 n n
X j2 X n 2
n jn
=− + cn +
2n 2 n
X j 2 n=+∞Y Z X n jn
2
n
Z(j) = N exp + exp − cn + dcn
n
2n n=−∞ 2 n
X j 2 +∞ Y Z n
n
= N exp dcn exp − c2n
n
2 n n=−∞ 2
X j2
n
= Z(0) exp
n
2 n
2.4. THÉORÈME DE WICK 21
Montrons que
X j2 1
Z
n
= dτ dτ 0 j(τ )A−1 (τ, τ 0 )j(τ 0 )
n
2 n 2
X qn (τ )qn (τ 0 )
où A−1 (τ, τ 0 ) = satisfait
n
AA−1 = δ(τ − τ 0 )
1
Z X j2
j(τ )A−1 (τ, τ 0 )j(τ 0 ) = n
2 2n
Chapitre 3
Méthodes semi-classiques
Sommaire
3.1 Fonction d’onde et trajectoires classiques . . . . . 23
3.2 Calcul du propagateur semi-classique . . . . . . . 24
3.3 Expression du propagateur semi-classique en di-
mension 1 . . . . . . . . . . . . . . . . . . . . . . . . 27
3.4 Propagateur dans l’espace E . . . . . . . . . . . . . 30
3.5 Effet tunnel-Etude du double puits . . . . . . . . . 34
3.6 Approximation du gaz dilué . . . . . . . . . . . . . 36
3.7 Calcul du déterminant : généralités . . . . . . . . 36
3.8 Evaluation du déterminant dans l’approximation
du gaz dilué . . . . . . . . . . . . . . . . . . . . . . . 37
3.9 Méthode de Laplace et méthode de la phase sta-
tionnaire . . . . . . . . . . . . . . . . . . . . . . . . . 39
Le qualificatif semi-classique désigne tous les phénomènes de transition entre

la mécanique quantique et la mécanique classique. L’échelle caractéristique des
phénomènes quantiques est donnée par la constante de Planck ~. Lorsque les
actions mises en jeu sont grandes devant ~ on s’attend à ce que les équations
quantiques se réduisent aux équations classiques. A première vue le problème
semble analogue à celui du passage de la mécanique relativiste à la mécanique
non relativiste, passage qui se fait simplement en prenant la limite c → ∞. La
difficulté du problème vient de ce qu’il ne suffit pas de poser ~ = 0 dans les
équations classiques. La limite ~ = 0 est hautement singulière ; pour l’analyser
il faut mettre en oeuvre des méthodes asymptotiques telles que la méthode de la
phase stationnaire ou la méthode de Laplace que nous rappelons en appendice.
3.1 Fonction d’onde et trajectoires classiques

Le passage de l’optique ondulatoire, théorie décrite par les équations de Maxwell,
à l’optique géométrique dans laquelle les entités de base sont les rayons lumineux
peut être réalisé par la méthode eikonale.
En écrivant une composante quelconque du champ électromagnétique sous la
forme
E = aeiψ
23
24 CHAPITRE 3. MÉTHODES SEMI-CLASSIQUES
On montre que les rayons lumineux sont les trajectoires orthogonales aux sur-
faces d’onde ψ = cste. Ils sont définis par la relation ~k = ∇ψ, ~ de même la
pulsation est donnée par ω = − ∂ψ∂t . La limite de l’optique géométrique, λ petit
devant les dimensions caractéristiques d des instruments d’optique, implique que
la phase varie de façon appréciable sur une distance d’ordre d. Par conséquent
l’eikonale ψ est une quantité grande.
Pour transposer cette méthode à la mécanique quantique on pose :
iS
ψ = ae ~
où S >> ~
En reportant cette expression dans l’équation de Schrödinger
∂ψ ~2
i~ =− ∆ψ + V (~x) ψ
∂t 2m
On obtient le système d’équations aux dérivées partielles
∂S 1 ~ 2 ~2 4a
+ ∇S + V (~x) = (3.1)
∂t 2m 2m a
!
∂ 2 ~
∇S
a + div a2 =0 (3.2)
∂t m
Si on néglige le membre de droite de l’équation (1), on reconnait l’équation de

Hamilton-Jacobi satisfaite par l’action classique S.
L’équation (2) s’interprète comme une équation de continuité pour la densité
~
gradS p
~
ρ = |ψ|2 = a2 et le courant classique ~j = ρ~v où ~v = =m .
m
Par conséquent quand ~ est petit on peut décrire le système en terme d’un flot de
particules transportées le long des trajectoires classiques perpendiculairement
aux surfaces d’ondes S = cste.
3.2 Calcul du propagateur semi-classique

On se propose de calculer le propagateur en dimension d’espace d=1
Z 0
c(t0 −t)
i t
iH
Z hm i
0 −
< q |e ~ |q >= Dq (τ ) exp dτ q̇ 2 (τ ) − V (q (τ ))
~ t 2
3.2. CALCUL DU PROPAGATEUR SEMI-CLASSIQUE 25
dans la limite ~ → 0.
La méthode de la phase stationnaire nous invite à rechercher les solutions qc (τ )
qui rendent stationnaire l’action classique
Z t0 hm i
Sc {q (τ )} = dτ q̇ 2 (τ ) − V (q (τ ))
t 2
avec les conditions aux limites {q (t) = q, q (t0 ) = q 0 }.

Ces solutions définissent les trajectoires classiques qc (τ ), solutions de l’équation
du mouvement
d2 ∂V
m qc (τ ) = − (qc (τ ))
dt2 ∂q
L’approximation semi-classique consiste à développer l’action autour de qc (τ )

en posant √
q (τ ) = qc (τ ) + ~η (τ )
avec η (t) = η (t0 ) = 0.
En se limitant aux termes quadratiques il vient
t0
m 2 1 ∂2V
Z
(qc ) η + O η 3
2

S {q (τ )} = S {qc (τ )} + ~ dτ η̇ − 2
t 2 2 ∂q
Nous poserons désormais Sc = S {qc (τ )}. A cet ordre le propagateur s’écrit

Z Z t0
0 b (t0 −t)
− ~i H i S~c i
< q |e |q >= e Dη (τ ) exp η (τ ) M η (τ ) dτ
2 t
2 2
d ∂ V
où M = −m dτ 2 − ∂q 2 (qc (τ )) est un opérateur différentiel de type Sturm-
Liouville.
On est donc ramené à calculer une intégrale fonctionnelle gaussienne associée à
l’opérateur de Sturm-Liouville M . Soit {φn (τ )} un ensemble complet de fonc-
tions propres de M vérifiant les conditions aux limites φn (t) = φn (t0 ) = 0. Elles
vérifient
M φn (τ ) = λn φn (τ )
ainsi que les relations d’orthogonalité
Z t0
φn (τ ) φm (τ ) dτ = δnm
t
On peut développer la fluctuation autour de la trajectoire classique sur cette

base de fonctions
∞
X
η (τ ) = cn φn (τ )
n=1
Nous pouvons définir la mesure Dη (τ ) par l’expression

N
Y
Dη (τ ) = lim dcn
N →∞
ṅ=1
Le terme exponentiel s’écrit

Z t0 ∞
X
η (τ ) M η (τ ) = λn c2n
t n=1
Par conséquent
Z Z t0
i
Dη (τ ) exp η (τ ) M η (τ ) dτ =
2 t
N Z ∞ N r
Y i Y 2π iπ
lim dcn exp λn c2n = lim e4
N →∞
n=1 −∞ 2 N →∞
n=1
λn
où nous avons supposé toutes les valeurs propres positives.
Pour se débarrasser du facteur de normalisation, on pourra évaluer le rapport
de 2 intégrales fonctionnelles
R t0 s
∞ (o)
Dη (τ ) exp 12 t η (τ ) M η (τ ) dτ
R
Y λn
t 0 =
Dη (τ ) exp 1 λn
R R
2 η (τ ) Mo η (τ ) dτ
t n=1
Lorsque le produit infini converge, on peut ainsi définir le déterminant de l’opérateur

M
∞
Y λn
detM ≡ (o)
n=1 λn
Application : propagateur de l’oscillateur harmonique
mω 2 2
V (q) = q
2
Par conséquent
d2 ∂2V d2
M = −m − = −m − mω 2
dτ 2 ∂q 2 dτ 2
Fonctions propres
nπ (τ − t)
φn (τ ) = sin ,n > 0
t0 − t
Valeurs propres " #
n2 π 2 2
λn = +m 2 −ω
(t0 − t)
Pour t0 − t < ωπ = T2 où T est la période du mouvement on remarque que toutes
les valeurs propres sont positives.
d2
En choisissant comme opérateur de référence Mo = −m dτ 2 avec les mêmes
conditions aux limites η (t) = η (t0 ) = 0 il vient

∞
" #
Y λn Y n2 π 2 − ω 2 (t0 − t)2
detM = =
(0)
n=0 λn n>0
n2 π 2
∞
!
2
Y ω 2 (t0 − t)
= 1−
n=1
n2 π 2
3.3. EXPRESSION DU PROPAGATEUR SEMI-CLASSIQUE EN DIMENSION 127
En utilisant la formule d’Euler

∞
x2

sinπx Y
= 1− 2
πx n=1
n
Il vient
sin ω (t0 − t)
det M =
ω (t0 − t)
Pour fixer les normalisations, remarquons que l’opérateur Mo intervient dans le
calcul du propagateur d’une particule libre.
Au chapitre 1 nous avons montré que
r
iH
bo m So
< q 0 | exp − (t0 − t) |q >= exp i
~ 2πi~ (t0 − t) ~
Dans l’approximation semi-classique on a
< q 0 | exp − i~H (t0 − t) |q > 1 i

b
=√ exp (Sc − So )
< q 0 | exp − iH~ 0 (t0 − t) |q > det M ~
b
Par conséquent
s
ω (t0 − t)
r
0 iH
b m Sc
< q | exp − (t0 − t) |q > = exp i
~ 2πi~ (t0 − t) sin ω (t0 − t) ~
r
mω Sc
= exp i
2πi~ sin ω (t0 − t) ~
où
mω 02
q + q 2 cos ω (t0 − t) − 2qq 0

Sc =
2 sin ω (t0 − t)
Exercices
1. Utiliser l’expression du propagateur pour calculer l’énergie et la fonction

d’onde de l’état fondamental
2. Montrer que la fonction de partition de l’oscillateur harmonique s’écrit
1
T race exp −β H = ωβ~
2 sinh 2
3.3 Expression du propagateur semi-classique en

dimension 1
Une extension du formalisme précédent permet de calculer le propagateur semi-
classique pour un potentiel arbitraire. La formule finale, établie par Van-Vleck
en 1928 s’écrit
X1 iM π ∂2S 1 i
< q 0 t0 |qt >= √exp − |− | 2 exp S (q, q 0 , t0 − t)
2πi~ 2 ∂q∂q 0 ~
où M ∈ N est le nombre de points conjugués, c’est le nombre de points de la tra-
jectoire où l’approximation semi-classique est mise en défaut ( M=2 lors d’une
réflexion sur un mur avec condition de Dirichlet ). Le signe somme signifie qu’il
faut sommer sur toutes les trajectoires classiques.
Exercice :
1. vérifier cette formule dans le cas de l’oscillateur harmonique
2. méthode des images
Remarques :
La formule de Van-Vleck cesse d’être valable lorsque le préfacteur de la for-

∂2S
mule de Van-Vleck tend vers l’infini.
∂q∂q 0
Pour dégager la signification physique de cette condition, il est intéressant de
réécrire ce préfacteur en utilisant le formalisme de Hamilton-Jacobi.
Considérons le mouvement d’une particule classique avec les conditions initiales
q (t0 ) = q0 , p (t0 ) = p0 .
Sa position au temps t est
q1 (t) ≡ q (q0 , p0 , t)
Exprimons l’impulsion initiale en fonction de l’action S, on a
∂S
p0 = −
∂q0
Par conséquent
∂2S ∂p0
=−
∂q0 ∂q1 ∂q1
∂2S ∂q1
La divergence de implique que = 0.
∂q0 ∂q1 ∂p0
Cette condition exprime le fait qu’une petite variation de l’impulsion initiale
n’affecte pas le point final. Il existe dans ce cas une infinité de trajectoires clas-
siques reliant q0 et q1 dans le temps t.
q(t) = q1
ρo
ρ(t) = q1
q (to ) = qo
3.3. EXPRESSION DU PROPAGATEUR SEMI-CLASSIQUE EN DIMENSION 129
Illustration : oscillateur harmonique
La trajectoire classique satisfaisant q (t0 ) = q0 et q (t1 ) = q1 s’écrit q (τ ) =

a cos ω (τ − t0 ) + b sin ω (τ − t0 ) 0 ≤ τ < t
où
a = q0
q1 − q0 cosω (t1 − t0 )
b=
sinω (t1 − t0 )
n
Pour ω (t1 − t0 ) = nπ, b n’est pas défini sauf si q1 = (−1) q0 . Il existe alors une
infinité de trajectoires connectant q0 et q1 .
q(τ )
qo
π
ωτ
−qo
Pour comprendre le rôle joué par ces points singuliers dans l’approche fonc-
tionnelle, revenons à l’expression des valeurs propres de l’opérateur de fluctua-
tion
" #
n2 π 2 2
λn = m 2 −ω
(t0 − t)
Si |t − t0 |ω = nπ , la nième valeur propre de cet opérateur devient nulle. Par

1
conséquent √det = ∞, l’approximation semi-classique cesse d’être valable car
l’action n’est plus une fonctionnelle quadratique dans la perturbation δq. Dans
le cas de l’oscillateur harmonique on franchit pour la première fois un point
conjugué pour n = 1 soit (t0 − t) = ωπ .
Dès que t0 − t > ωπ il apparait une valeur propre négative dans le spectre des
fluctuations.
λn
π
ω < t0 − t < 2π
ω
par conséquent le det change de signe donc

1 1 π 1
√ →√ = e−i 2 √
det −det det
Il va donc apparaı̂tre une phase supplémentaire dans le propagateur. Si on
franchit M points conjugués on obtient la phase exp −i M2π .
Remarque : en dimension supérieure l’ensemble des points conjugués définit
une surface appelé caustique. En optique, la caustique est la surface qui sépare
la région sombre de la région éclairée. On peut la définir géométriquement
comme l’enveloppe des rayons lumineux réfléchis par un miroir. Sur la caus-
tique l’éclairement est donc maximum puisqu’il y a une accumulation de rayons.
caustique
Miroir cylindrique à section
circulaire
Notons que l’optique géométrique prédit un éclairement discontinu à la tra-

versée d’une caustique. Dans le cadre de l’optique ondulatoire, on peut calculer
l’intensité lumineuse de part et d’autre de la caustique. Elle s’exprime en terme
de fonctions d’Airy qui jouent un rôle important dans la description de l’arc en
ciel.
3.4 Propagateur dans l’espace E

1
On considère la résolvante G (E) = H−E . Ses éléments de matrice sont des
fonctions méromorphes de E admettant pour pôles les valeurs propres de H
dans le cas d’un spectre discret.
Nous nous proposons de calculer le propagateur ”dans l’espace E”
1
G (q, q 0 , E) =< q 0 | |q >
H −E
3.4. PROPAGATEUR DANS L’ESPACE E 31
Il s’exprime à partir du propagateur de Feynman

1 ∞
Z
0 iEτ
G (q, q , E) = − dτ K (q 0 τ |q0) exp
i~ 0 ~
intégrale convergente pour E = E0 + i avec > 0.
Partons de la formule de Van Vleck en traitant le terme exponentiel par la
méthode de la phase stationnaire
s
1 1 iM π ∂ 2 Sc
G (q, q 0 , E) = − √ exp − −
i~ 2πi~ 2 ∂q∂q 0
Z ∞
i
dτ exp + [Sc (q, q 0 , τ ) + Eτ ]
0 ~
La phase est stationnaire pour
∂
E=− Sc (q, q 0 , τ ) |τ =t
∂τ
Par conséquent E peut s’interprèter comme l’énergie d’une particule décrivant
la trajectoire classique connectant les points q et q 0 dans le temps t. Posons
A (t) = Sc (q, q 0 , t) + Et où t est une fonction implicite de E déterminée par la
condition de stationnarité précédente. Il vient
s
1 1 iM π ∂2S
G=− √ exp − −
i~ 2πi~ 2 ∂q∂q 0
Z ∞
A(t) i ∂2S 2
ei ~ dτ exp |
2 τ =t
(τ − t)
0 2~ ∂τ
En utilisant H = pq̇ − L nous pouvons réécrire A (t) sous la forme
Z t Z q0
A (t) = (E + L) dτ = p dq
0 q
s
2π~
L’intégrale sur τ donne, à une phase près, ∂2S
. Par conséquent
∂τ 2
v
u ∂2S Z
1u ∂q∂q 0
G ∼ t ∂ 2 S exp i p dq
~ ∂τ 2
Posons W (q, q 0 , E) = p dq = S + Et et réexprimons le facteur préexponentiel

R
en fonction de E = − ∂t∂
S (q, q 0 , t).
On a
∂W ∂W ∂S
(q, q 0 , E) = t (q, q 0 , E) = (q, q 0 , t)
∂E ∂q ∂q
Donc
∂2S ∂2W 0 ∂ 2 W ∂E
= (q, q , E) +
∂q∂q 0 ∂q∂q 0 ∂q∂E ∂q 0
∂2W ∂2W ∂2S

= 0
(q, q 0 , E) −
∂q∂q ∂q∂E ∂t∂q 0
Or
∂S ∂W
=
∂q 0 ∂q 0
∂2S ∂ 2 W ∂E ∂2W ∂2S

= = −
∂t∂q 0 ∂q 0 ∂E ∂t ∂q 0 ∂E ∂t2
Par conséquent
∂2S ∂2W ∂2W ∂2W ∂2S

0
= 0
+
∂q∂q ∂q∂q ∂q∂E ∂q 0 ∂E ∂t2
∂2S ∂2W
∂q∂q 0 ∂q∂q 0 ∂2W ∂2W
∂2S
= ∂2S
+
∂t2 ∂t2
∂q∂E ∂q 0 ∂E
∂S
En utilisant ∂t = −E il vient
∂2S ∂E 1
=− = − ∂2W
∂t2 ∂t ∂E 2
Par conséquent
∂2S
∂q∂q 0 ∂2W ∂2W ∂2W ∂2W
∂2S
=− 0
+
∂t2
2
∂E ∂q∂q ∂q∂E ∂q 0 ∂E
Nous obtenons ainsi l’expression finale du propagateur à une phase près

s
X 1 ∂2W ∂2W ∂2W ∂2W W
G (q, q 0 , E) ∼ 0
− 2 0
exp i
~ ∂q∂E ∂q ∂E ∂E ∂q∂q ~
Application : formule BKW

Considérons le mouvement d’une particule dans le puit de potentiel suivant
V(q)
qo q1
q q’
q
qo q1
3.4. PROPAGATEUR DANS L’ESPACE E 33
L’ensemble des trajectoires allant de q à q 0 est constitué de la trajectoire

”directe” qq 0 , de la trajectoire réfléchie sur les points q0 , q1 , ainsi que toutes ses
répétitions.
La fonction de Green G (q, q 0 , E) s’obtient en sommant sur toutes ces trajectoires
classiques.
Pour une trajectoire se réfléchissant n fois, le terme exponentiel s’écrit
Z q0
Wn = p dq + nW
q
Z q0 I
= p dq + n p dq
q
π
En admettant qu’il apparait une phase supplémentaire de 2 à chaque réflexion,
il vient
∞ I
X in 1
exp p dq + inπ =
1 + exp ~i p dq
H
n=0
~
Les pôles de G sont donc donnés par la condition de quantification BKW .

I
1
p dq = 2 k + π~
2
Par conséquent l’énergie du système ne peut prendre que des valeurs discrètes,
celles pour lesquelles l’aire de l’espace de phase enclose par la trajectoire clas-
sique est un multiple demi entier de 2π~.
ρ
ρ2
2m
+ V (q) = E
q
3.5 Effet tunnel-Etude du double puits

1 2
V (q) = (q − µ2 )2
4
~2 d 2
H=− + V (q)
2m dq 2
−µ µ
On se propose de calculer
T T TH
< −µ, | µ, − >=< −µ | e− ~ | µ > dans la limite T → ∞
2 2
Ecrivons les états propres de H comme des combinaisons linéaires d’ états lo-
calisés dans chacun des puits
ψo (q − µ) + ψo (q + µ)
ψ1 (q) = √ énergie Eo − ∆E
2
ψo (q − µ) − ψo (q + µ)
ψ2 (q) = √ énergie Eo + ∆E
2
Si à l’instant initial la particule est localisée dans le puits de droite
ψ1 (q) + ψ2 (q)
ψo (q − µ) = √
2
A l’instant T le système sera dans l’état
1 E0 −∆E E0 +∆E

√ ψ1 (q)e−T ( ~ ) + ψ2 (q)e−T ( ~ )
2
Par conséquent
Z +∞
TH 1 E0 −∆E E0 +∆E

< −µ | e− ~ | µ >= ψo (q+µ) √ ψ1 (q)e−T ( ~ ) + ψ2 (q)e−T ( ~ )
−∞ 2
Eo T T ∆E
= e− ~ sh
~
TH Eo T T ∆E
lim < −µ | e− ~ | µ >= e− ~ sh
T →∞ ~
3.5. EFFET TUNNEL-ETUDE DU DOUBLE PUITS 35
Par l’intégrale fonctionnelle on voit que

Z Z T2 2
TH 1 mq̇
< −µ | e− ~ | µ >= Dq exp − + V (q) dτ
~ − T2 2
∂V
L’équation du mouvement mq̈ = ∂q , décrit le mouvement d’une particule dans
le potentiel −V (q).
Pour simplifier posons m=1, les équations du mouvement deviennent
q̈ = q(q 2 − µ2 )
On a l’intégrale première 12 q̇ 2 = V (q)+E. Pour construire des solutions d’action

finie il faut que pour q = µ la vitesse q̇ = 0 , par conséquent E = 0. Les solutions
s’écrivent
µ(τ − τo )
q(τ ) = +µth √
− 2
On dispose donc d’une famille à un paramètre de solutions d’action finie qui
interpolent entre ±µ
−µ µ
L’action euclidienne s’écrit
+∞ Z +∞
q̇ 2
Z
SE = + V (q) dτ = 2V (q)dτ
−∞ 2 −∞
En utilisant dτ = √ dq nous obtenons

2V (q)
Z +µ p
SE = 2V (q)dq
−µ
Problème : calcul du déterminant des petites fluctuations

d2 ∂2V
Y
det − 2 + = n
dτ ∂q 2 n
d2 ∂2V

− 2+ Xn (τ ) = n Xn (τ )
dτ ∂q 2
Difficulté : il existe une valeur propre nulle puisque la solution classique satisfait
∂V ... ∂2V
q̈c = |qcl ⇒ q c = q̇c
∂q ∂q 2
Par conséquent Xo (τ ) = q̇c (τ ) est le fondamental de H ( solution sans noeud)

donc det = 0.
Ce mode nul reflète l’invariance par translation. Si q(τ ) est solution , q(τ −τo )
est encore solution ∀τo . Par conséquent il ne faut pas calculer le déterminant
dans tout l’espace fonctionnel, mais il faut calculer séparemment la contribution
venant de la coordonnée collective (centre de masse de l’instanton).
3.6 Approximation du gaz dilué

Pour pouvoir utiliser la formule
TH Eo T T ∆E
lim < −µ | e− ~ | µ >= e− ~ sh
T →∞ ~
Il faut garder T fini dans l’intégrale de chemin et ne prendre la limite T → ∞
qu’à la fin du calcul. Or la solution classique obtenue n’est strictement valable
que pour T → ∞.
Pour contourner cette difficulté on imagine que la particule ne part pas de

+µ à −∞ mais de µ − avec une énergie V (µ − ) ∼ 0 pour t = − T2 . Elle décrit
la trajectoire approchée suivante
µ− ∈
q(τ )
τ3 τ2 τ1
− T2 T
2
∈ −µ
constituée d’un antiinstanton, d’un instanton et d’un antiinstanton. Les
points τ3 < τ2 < τ1 sont les centres des instantons sur lesquels on doit intégrer
(coordonnées collectives).
3.7 Calcul du déterminant : généralités

Z Rt
< qf | e−
TH
~ | qi >= Dq(τ )e− ~
1
o [ m2 q̇2 +V (q)]dτ
La méthode de Laplace donne

TH Sc 1
< qf | e− ~ | qi >= e− ~ r
d 2 ∂2V
det −m dt2 + ∂q 2
∂V
La solution stationnaire est donnée par mq̈ = ∂q équation du mouvement d’une
particule de masse m dans le potentiel −V
m 2
E=+ q̇ − V (q) est une intégrale première
2
Exemple d’un puit isolé
3.8. EVALUATION DU DÉTERMINANT DANS L’APPROXIMATION DU GAZ DILUÉ37
V(q)
q
TH
Calculons < q = 0 | e− ~ | q = 0 > dans la limite ~ → 0. Les équations du
mouvement donnent une trajectoire classique unique connectant q = 0 à q = 0
dans le temps T . Elle s’écrit
q(t) = 0 ⇒ q̇(t) = 0 ⇒ S = 0
TH 1 −S
< q = 0 | e− ~ | q = 0 >= q e
~
d2
det −m dt2 + mω 2
∂2V
où mω 2 = ∂q 2
Nous avons déjà rencontré ce déterminant dans l’étude de l’oscillateur har-

monique
iH(t0 − t)
r
iπ mω
< q | exp − | q >= exp −
~ 4 2π~ sin ω(t0 − t)
posant t0 − t = −iT ; sin ω(t0 − t) = shωT

i
r r
TH iπ imω mω
< q | exp − | q >= exp − =
~ 4 2π~shωT 2π~shωT
Dans la limite T → ∞
r
− T~H mω − ωT
lim < q = 0 | e |q=0>= e 2
T →∞ π~
T Eo
= ψo2 (0)e− ~
Ce résultat est bien en accord avec les formules classiques

mω 14 ~ω
ψo (q) |q=0 = ; Eo =
π~ 2
mω 14 1 mωq 2
en effet ψo (q) = exp −
π~ 2 ~
3.8 Evaluation du déterminant dans l’approxi-

mation du gaz dilué
L’opérateur de fluctuation s’écrit
d2 ∂2V
M =− + |q=qc
dτ 2 ∂q 2
avec
1 2
V (q) = (q − µ2 )2
4
∂2V
= 3q 2 − µ2
∂q 2
Les solutions instanton (ou antiinstanton) sont
µ
qc = +µ th √ (τ − τo )
− 2
Donc
∂2V 2 3µ2
= 2µ −
∂q 2 ch2 √µ2 (τ − τo )
Le terme constant
∂2V
= 2µ2
∂q 2
décrit les fluctuations quadratiques au voisinage du minima |q| = µ.
On a dans ce cas
r r r
1 mω mω − Eo T
= → e ~
det M 2π~shωT π~
~ω
où Eo = 2 est l’énergie du fondamental d’un puit isolé.
En présence d’une configuration diluée constituée de n instantons/antiinstantons

on aura
1 Eo T
√ = K n e− ~
det M
où K est une constante calculable en terme de la solution classique.
Intégrant sur les coordonnées collectives τ1 , τ2 , ... on aura

Z Z T2 2
1 q̇ Eo T
Dq(τ ) exp − + V (q) dτ = e− ~
~ − T2 2
∞ Z T Z τ1 Z τn−1 n
X 2 So
dτ1 dτ2 dτn Ke− ~
impairs − T2 − T2 − T2
∞
Eo T X T n − So n Eo T
So

= e− ~ Ke ~ = e− ~ sh KT e− ~
impairs
n!
En comparant avec la formule initiale
TH Eo T T ∆E
< −µ | e− ~ | µ >= e− ~ sh
~
on obtient
Z µ
1 p
∆E = K~ exp − dq 2V (q)
~ −µ
3.9. MÉTHODE DE LAPLACE ET MÉTHODE DE LA PHASE STATIONNAIRE39
3.9 Méthode de Laplace et méthode de la phase

stationnaire
• Méthode de Laplace (voir Copson p 39)
Soit g (x) et h (x) deux fonctions réelles continues définies sur l’intervalle
fini ou semi infini a < x 6 b telles que
1. g (x) eth(x) est absolument intégrable sur l’intervalle ∀t > 0.
2. h (x) a un maximum dans l’intervalle [a, b] pour x = a tel que h0 (a) =
0 h00 (a) < 0 et le maximum de h (x) dans tout sous intervalle
fermé de [a, b] est strictement inférieur à h (∞).
3. h00 (x) est continu.
Rb q
Alors pour t → ∞ on a a g (x) eth(x) ∼ g (a) eth(a) −2thπ2 (a) .
• Généralisation : elle permet de traiter des intégrales de la forme

Rb
a
g (x) eth(x) dx dans la limite t → ∞. On suppose que h (x) a un maxi-
mum unique en SΣ [a, b]. En divisant l’intervalle [a, b] en deux intervalles
[a, b] et [a, b] on se ramène au cas précédent.
exemple : la fonction Γ (t) et la formule de Stialing

Z ∞
Γ (t + 1) = e−x xt dx
0
posons x = tu
Z ∞
Γ (t + 1) = e−tu tt+1 ut du
0
Z ∞
t+1
=t et[logu−u] du
0
la fonction h (u) = logu − u a un maximum pour u =p1. Sur chaque

intervalle [1, ∞] [01] on obtient la même contribution e−t 2t
π
.
D’où √ 1
Γ (t → ∞) ∼ t → ∞ 2π tt+ 2 e−t
• Méthode de la phase stationnaire

Si h (x) a un point stationnaire pour c [a, b]
1. h00 (c) > 0
s
Z b
ith(x) 2π ith(c)+ iπ 1
I= g (x) e dx = e 4 + 0
a th00 (c) t
2. h00 (c) < 0

s
Z b
ith(x) 2π ith(c)+ iπ 1
I= g (x) e dx = − 00 e 4 + 0
a th (c) t
Chapitre 4
Applications à la physique
des polymères
Sommaire
4.1 Modélisation d’un polymère comme marche aléatoire 41
4.2 Polymère dans un potentiel . . . . . . . . . . . . . 44
4.3 Elasticité des biomolécules . . . . . . . . . . . . . . 46
4.3.1 Modèle gaussien discret . . . . . . . . . . . . . . . . 46
4.3.2 Modèle du ver . . . . . . . . . . . . . . . . . . . . . 49
4.1 Modélisation d’un polymère comme marche

aléatoire
Un polymère idéal est constitué d’une chaı̂ne de monomères identiques sans in-
teractions. On peut le modéliser comme un ensemble de N batonnets librement
articulés dont l’orientation relative est distribuée de facon uniforme.
Supposons qu’une des extrémités de la chaı̂ne reste fixée en r~0 . Selon les orien-
tations relatives des différents monomères, il existe un grand nombre de confi-
gurations de longueur N b dont l’autre extrémité occupe la position r~N = ~r. On
se propose de calculer la distribution de probabilité de l’extrémité du polymère.
~
r
~
ro
41
42 CHAPITRE 4. APPLICATIONS À LA PHYSIQUE DES POLYMÈRES
Hypothèse :
• La position ~rj du monomère j ne dépend que de celle du monomère

précédent.
On se propose de calculer la distribution de probabilité PN (~r | ~r0 ) de l’extrémité
~r d’un polymère fixé en ~r0 .
Probabilité de transition
1
P (~rj − ~rj−1 ) = δ (| ~rj − ~rj−1 | −b)
4πb2
Partons de la récurrence
Z
PN (~r | ~r0 ) = P ~r − r~0 PN −1 r~0 | ~r0 d3 r~0
Posons
Z
~
ΓN ~k = PN (~r | ~r0 ) eik~r d3~r

Soit λ ~k la fonction de structure
Z

~ ~ 1 sin kb
λ k := d3~reik~r δ (| ~r | −b) =
4πb2 kb
La solution de la récurrence donne

ΓN ~k = ΓN −1 ~k λ ~k
Conditions initiales
Z
~
Γ0 ~k = P0 (~r | ~r0 ) eik~r d3~r
Z
~ ~
= δ 3 (~r − ~r0 ) eik~r d3~r = eik~r0
Il vient
h iN
~
ΓN ~k = λ ~k eik~r0
Par conséquent
Z N
1 ~ sin kb
PN (~r | ~ro ) = 3 e−ik(~r−~ro ) d3~k
(2π) kb
Analysons le comportement de longues chaı̂nes dans la limite d’échelle

N → ∞, b → 0, N b2 fixé.
A k fixé on obtient
N !
2 2
sin kb (kb) N (kb)
= exp N log 1 − ∼ exp −
kb 6 6
4.1. MODÉLISATION D’UN POLYMÈRE COMME MARCHE ALÉATOIRE43
N b2
Posons 3 = s, il vient
Z
1 ~ sk2
P [~r(s) = ~r | ~r(0) = ~r0 ] = 3 d3~k
e−ik(~r−~r0 )− 2
(2π)
23
1 1 2
= exp − (~r − ~r0 )
2πs 2s
La convergence vers une distribution gaussienne résulte du théorème de la limite

~ j+1 = R
centrale. Le polymère idéal définit une marche aléatoire R ~ j + ~b. Cette
marche aléatoire converge, dans la limite d’échelle considérée, vers un mouve-
ment brownien dont nous venons de calculer la probabilité de transition.
Remarques :
1. La distance quadratique moyenne bout à bout entre les extrémités du

polymère est donnée par
2
< (~r(s) − ~r(0)) >= N b2 = 3s
La longueur totale du polymère est infinie dans la limite continue en effet
N b2 fini

⇒ Nb → ∞
b→0
~
r √
b N ~
ro
2. Dans la limite continue P [~r(s) = ~r | ~r(0) = ~r0 ] est solution de
∂P 1
= ∆ P avec lim P [~r(s) = ~r | ~r(0) = ~r0 ] = δ (~r − ~r0 )
∂s 2 s→0
Par conséquent P coı̈ncide avec le propagateur euclidien et peut donc être

représentée par l’intégrale de chemin
Z ~
r (s)=~
r Z s 2
1 d~r
P [~r(s) = ~r | ~r(0) = ~r0 ] = D~r(τ ) exp − dτ
~
r (0)=~
r0 2 0 dτ
Cette expression permet d’interpréter la variable s comme l’abscisse cur-

viligne le long du polymère (bien que dimensionnellement s ∼ b2 ).
3. Dans un polymère physique, les monomères interagissent entre eux et il

existe des contraintes d’autoévitement de sorte que
2
< [~r(s) − ~r(0)] >∼ N 2ν b2
où l’exposant ν est donné avec une bonne approximation par la formule
de Flory (d 6 4)
3
ν=
d+2
• Pour d = 1 ν = 1 r ∼ N en effet un polymère autoévitant
unidimensionnel est complètement étiré.
• Pour d = 4 ν = 36 = 21 r2 ∼ N
On retrouve un système diffusif car l’interaction est non pertinente
pour d > 4.
4.2 Polymère dans un potentiel

On considère un polymère en présence d’une surface solide. Par exemple un po-
lymère en solution au voisinage d’une paroi. On suppose que chaque monomère
interagit avec la paroi par un potentiel U (~r).
Configuration typique
~r1
~ro
Le poids statistique attaché à une chaı̂ne d’extrémités fixées ~r0 , ~r1 peut se
1
calculer en supposant le système à l’équilibre à température T = kβ .
Z Z t " 2 #
1 1 d~r
D~r(s)δ [~r(0) − ~r0 ] δ [~r(t) − ~r1 ] exp − + βU (~r(s)) ds
Z 0 2 ds
Le numérateur peut s’écrire

1
< ~r1 | e−tH | ~r0 > où H = − ∆ + βU
2
Z est un facteur de normalisation
Z
Z= d~r0 d~r1 < ~r1 | e−tH | ~r0 >
4.2. POLYMÈRE DANS UN POTENTIEL 45
Afin de caractériser l’extension spatiale du polymère, calculons la distance qua-

dratique moyenne entre les extrémités
2
R2 (t) =< [~r(t) − ~r(0)] >
2
d~r0 d~r1 (~r0 − ~r1 ) < ~r0 | e−tH | ~r1 >
R
2
R (t) = R
d~r0 d~r1 < ~r0 | e−tH | ~r1 >
Développons en fonctions propres
X
< ~r0 | e−tH | ~r1 >= ψn (~r0 ) ψn (~r1 ) e−tEn
• A basse température β → ∞, si le potentiel d’interaction avec la paroi est

attractif on aura un état lié (sous certaines conditions). En ne retenant
que cette contribution on obtient
2
d~r0 d~r1 (~r0 − ~r1 ) ψ0 (~r0 ) ψ0 (~r1 ) e−tE0
R
R2 (t) w R
d~r0 d~r1 ψ0 (~r0 ) ψ0 (~r1 ) e−tE0
d~r0 ~r02 ψ0 (~r0 ) d~r1 ψ0 (~r1 ) + d~r1 ~r12 ψ0 (~r1 ) d~r0 ψ0 (~r0 )
R R R R
= R R
d~r0 ψ0 (~r0 ) d~r1 ψ0 (~r1 )
R
le terme croisé disparait car d~r0 ~r0 ψ0 (~r0 ) = 0 si ψ0 (~r0 ) est à symétrie
radiale donc R2 (t) tend vers une constante indépendante de t.
Le polymère est dans une phase localisée
• A haute température β → 0, on peut remplacer < ~r1 | e−tH | ~r0 > par le
propagateur libre, par conséquent
2
d~r0 d~r1 exp − (~r1 −~
R r0 ) 2
2 2t (~r1 − ~r0 )
R (t) = r1 −~ 2
d~r0 d~r1 exp − 2tr0 )
(~
R
En introduisant le centre de masse
~ = ~r0 + ~r1 et la variable relative ~r = ~r1 − ~r0

R
2
il vient
2
d~r ~r2 exp − ~r2t √
R
2
R (t) = R r2
~
= 3t
d~r exp − 2t
On a donc
√ √ √
< R2 > = 3t = b N → ∞
Ceci suggère l’existence d’ une transition de délocalisation pour une température

critique T = Tc . Question : commment se comporte R2 (t) au voisinage de
T = Tc ?
4.3 Elasticité des biomolécules

4.3.1 Modèle gaussien discret
Revenons au modèle discret de N monomères de longueur b pour lequel nous
avons écrit la probabilité d’une configuration sous la forme
Z
PN (~r | ~r0 ) = P (~r − ~rN −1 )P (~rN −1 − ~rN −2 ) . . . P (~r1 − r~0 )d~r1 d~r2 . . . d~rN −1
Cette expression traduit l’indépendance statistique des différents éléments de la

chaı̂ne
~ro
~r F~
~r1
~rN −1
On suppose la chaı̂ne attachée au point ~r0 à l’une de ses extrémités et soumise
à une force F~ constante s’exercant sur l’ extrémité ~r = ~rN .
Dans l’ensemble canonique le poids statistique d’une configuration peut donc
s’écrire
e−βU (~r) P (~r − ~rN −1 )P (~rN −1 − ~rN −2 ) . . . P (~r1 − ~r0 )
où U (~r) = −F~ (~r − ~r0 ) est l’énergie potentielle (travail de la force F~ ). Ecrivons
U (~r) sous la forme
XN
U (~r) = −F~ (~ri − ~ri−1 )
i=1
4.3. ELASTICITÉ DES BIOMOLÉCULES 47
Par conséquent la fonction de partition canonique pourra se mettre sous la forme
N
Z Y
~
Z= eβ F (~ri −~ri−i ) P (~ri − ~ri−1 )d~r1 . . . d~rN
i=1
On intègre sur toutes les variables, y compris l’extrémité ~r du polymère.
Dans le modèle de N maillons rigides on a

1
P (~ri − ~ri−1 ) = δ(| ~ri − ~ri−1 | −b)
4πb2
L’intégration va se factoriser en une partie radiale triviale et une partie angu-
laire. Le ième maillon sera repéré par ses angles (θi , ϕi ) par rapport à la direction
de la force.
~
F
~ri − ~ri−1
θi
~ri−1
ϕi
Par conséquent
Z XN N
Y
Z ∼ exp β F b cos θi sin θi dθi dϕi
i=1 i=1
La fonction de partition se factorise sous la forme

Z +1 N N
4πshβF b
Z = 2π eβF bu du =
−1 βF b
Nous en déduisons l’énergie libre
F = −N kT log Z = −N kT [log shβF b − log βF b]
et l’élongation de la chaı̂ne

∂F 1
l=− = N kT βb coth βF b −
∂F F

1
= N b coth βF b −
βF b
= N bL(x)
1
où L(x) = coth x − x est la fonction de Langevin
Discussion : L’élongation est fonction du paramètre sans dimension x =

Fb
βF b = kT
• pour x → ∞ (force infinie ou température nulle) on obtient l = N b. Tous

les maillons s’orientent dans le sens de la force donnant ainsi à la chaı̂ne
sa longueur maximale
• pour x → 0 L(x) ' x3 . A petite force ou haute température on a
N b2 F
l=
3kT
On retrouve ainsi une loi linéaire avec un module d’élasticité proportionnel

à la température. C’est un comportement entropique qui traduit le fait que
le modèle repose essentiellement sur un comptage de configurations.
Expériences avec l’ADN

(voir Marko JF, Siggia ED, Science, 265,505-508,1995)
bF
L’expression précédente Nl b = 3kT fait apparaı̂tre une force typique F0 = kT
b
qui est du bon ordre de grandeur à condition de prendre pour b non pas la
distance entre bases (b = 3, 37Å) mais une distance effective 103 Å. Il vient
1 1
F0 = ×1, 6 × 10−19 × 3 = 4 × 10−14 N
40
|{z} 10 × 10−10
1
40 eV
L’accord avec la courbe expérimentale n’est pas satisfaisant. Il convient donc

de recourir à un autre modèle prenant en compte les propriétés d’élasticité de
l’ADN.
F
exp
Analogie avec le paramagnétisme

H = −~ ~ décrit le couplage d’un moment magnétique à un champ extérieur.
µB
Les moments magnétiques tendent à s’aligner dans la direction du champ, de
même que le polymère s’aligne dans la direction de la force.
~
B µ
~
4.3. ELASTICITÉ DES BIOMOLÉCULES 49
4.3.2 Modèle du ver

On modélise l’ADN comme un polymère semi-flexible. L’énergie élastique d’une
configuration de longueur l est
Z 2
B l ∂~t
Z l
E= ds − F~ ~tds
2 0 ∂s 0
~t(s)
où ~t est le vecteur tangent

le long de la chaı̂ne. Le premier terme représente
d~
t ~
n
une énergie de courbure ds = R et le second une énergie élastique de déformation
~t = d~r ⇒ F~ ~tds = F~ d~r

ds
La fonction de partition est donnée par l’intégrale de chemin
Z
Z = D~t exp −β E
D~t décrit la mesure des chemins tracés sur la sphère unité issus du point ~t0 =
~t (s = 0) et arrivant au point ~t = ~t(s = l)
~t
~to
Paramétrisant le vecteur ~t sous la forme (sinθ cosϕ, sinθ sinϕ, cosθ)

on obtient
βB l 2
Z Z l
βE = θ̇ + sin2 θϕ̇2 ds − β F cos θ ds
2 0 0
Le premier terme représente le lagrangien d’une particule de masse m = βB se
déplacant sur la sphère S2 . Le hamiltonien associé est
1 ∂2

1 1 1 ∂ ∂
H0 = − ∆2 = − sin θ +
2m 2βB sin θ ∂θ ∂θ sin2 θ ∂ϕ2
Le second terme peut s’interpréter comme un couplage dipôlaire. Le hamiltonien

complet s’écrit donc
1 ~2
L
H=− ∆ − βF cos θ = − βF cos θ
2Bβ 2Bβ
~ est l’opérateur de moment cinétique orbital. Par conséquent le poids sta-

où L
tistique d’un chemin de longueur l s’écrit
P l, ~t | 0, ~t0 ∼< ~t | exp −H l | ~t0 >

l 1 2
Hl = − ∆ − β BF cos θ
βB 2
βB = ξ a les dimensions d’une longueur appelée longueur de persistence

l 1 Fξ
Hl = − ∆− cos θ
ξ 2 kT
En l’absence de force extérieure

l
P l, ~t | 0, ~to ∼< ~t | exp − ∆ | ~t0 >

2ξ
• état fondamental de fonction d’onde Yoo et d’énergie nulle
• premier état excité Y1m d’énergie L(L + 1) = 2
Pour une chaı̂ne de longueur l telle que ~t0 //~k
1 3 l
P (l, ~t | 0, ~to ) ∼ + cos θ exp −
4π 4π ξ
Pour l’ADN la longueur de persistence ξ ∼ 500Å est beaucoup plus grande

que la distance entre bases. C’est cette propriété qui permet de décrire l’ADN
comme un modèle élastique continu.
Chapitre 5
Introduction au mouvement
Brownien
Sommaire
5.1 Histoire . . . . . . . . . . . . . . . . . . . . . . . . . 51
5.1.1 Physique . . . . . . . . . . . . . . . . . . . . . . . . 51
5.1.2 Mathématiques . . . . . . . . . . . . . . . . . . . . . 52
5.1.3 Théorie de la spéculation . . . . . . . . . . . . . . . 53
5.2 Définition d’une marche aléatoire . . . . . . . . . . 53
5.3 Marche aléatoire sur réseau . . . . . . . . . . . . . 54
5.3.1 Autre dérivation . . . . . . . . . . . . . . . . . . . . 56
5.3.2 Programme de simulation . . . . . . . . . . . . . . . 57
5.4 Limite continue . . . . . . . . . . . . . . . . . . . . . 57
5.4.1 Théorème . . . . . . . . . . . . . . . . . . . . . . . . 59
5.5 Notions de processus aléatoire . . . . . . . . . . . 59
5.5.1 Probabilités conditionnelles . . . . . . . . . . . . . . 60
5.5.2 Processus de Markov . . . . . . . . . . . . . . . . . . 60
5.5.3 Equation de Chapman-Kolmogorov . . . . . . . . . . 61
5.5.4 Processus aléatoire stationnaire . . . . . . . . . . . . 61
5.5.5 Théorème de Doob . . . . . . . . . . . . . . . . . . . 61
5.6 Exemples de processus stochastiques gaussiens . 63
5.6.1 Processus de Wiener . . . . . . . . . . . . . . . . . . 63
5.6.2 Processus d’Ornstein-Uhlenbeck . . . . . . . . . . . 64
5.6.3 Bruit blanc . . . . . . . . . . . . . . . . . . . . . . . 64
5.6.4 Construction de W (t) comme processus . . . . . . . 65
5.1 Histoire
5.1.1 Physique
1827 : L’histoire du mouvement brownien commence avec les observations du
botaniste anglais Georges Brown. En observant de fines particules en suspen-
sion dans l’eau il constate qu’elles sont animées d’un mouvement irrégulier et
51
52 CHAPITRE 5. INTRODUCTION AU MOUVEMENT BROWNIEN
incessant.
Tout au long du 19ème siècle plusieurs physiciens vont chercher à comprendre
l’origine du phénomène. Une explication plausible repose sur ”l’hypothèse” moléculaire,
les nombreux chocs entre la particule et les molécules environnantes seraient la
cause du mouvement observé.
1905 : Einstein donne une description théorique du phénomène. En utili-

sant la mécanique statistique, il calcule la densité de probabilité d’observer la
particule au point x. En dimension 1
1 x2
f (x, t) = √ exp −
4πDt 4Dt
Il en déduit que l’écart quadratique moyen < x2 (t) > croit linéairement avec le
temps t.
< x2 (t) >= 2Dt
Pour des particules sphériques de rayon a, il utilise la relation de Stokes donnant

le coefficient de friction η en fonction de la viscosité du fluide et en déduit le
kT
coefficient de diffusion D = 6πηa
1909 : Jean Perrin réalise des expériences qui permettent de mesurer D

R
et d’en déduire k = N , d’où la détermination du nombre d’Avogadro N =
23
6, 02.10 et une confirmation de ”l’hypothèse moléculaire”. En France cette
”hypothèse moléculaire” avait un certain nombre de détracteurs influents parmi
lesquels Marcellin Berthelot qui était considéré comme le pape de la chimie.
−→ Prix Nobel Jean Perrin (1926).
5.1.2 Mathématiques
En 1827 quand Brown publie ses travaux, les mathématiciens manipulent les
fonctions, les dérivent, les développent en série. La notion de fonction analytique
commence à voir le jour (eq de Cauchy Riemann 1814). Cependant les seules
fonctions qui font l’objet d’étude sont des fonctions régulières.
En 1872 Weierstrass exhibe une fonction continue non dérivable
∞
X
f (t) = an sin(bn t)
n=1
0 < a < 1 < ab

Ces objets sont alors considérés comme des ”monstres” par les mathématiciens
(Hermite). Le mouvement brownien a joué le rôle de monstre auprès des mathématiciens
selon la parabole de Benoı̂t Mandelbrodt.
”L’histoire des sciences regorge de sorciers et de contes de fées. Un
sorcier crée un monstre non par besoin ni par malice mais simplement
pour se prouver que la bête n’était point inconcevable. Le monstre
lâché, les paysans lui refusent l’entrée de leur village, car ses traits
les effraient autant qu’ils forcent leur incrédulité. Et puis un jour
une fée leur ouvre les yeux, le monstre est honnête homme et tout
prêt à les servir. On s’y habitue, on finit même par le trouver beau”.
5.2. DÉFINITION D’UNE MARCHE ALÉATOIRE 53
Le mouvement brownien comme nous le verrons plus loin est effectivement un

monstre puisqu’il s’agit d’un processus continu mais non dérivable. Sur l’inter-
valle [01] il peut être construit comme la série de Fourier aléatoire :
∞ √
X cn 2
x(t) = sin(nπt)
n=1
nπ
dans laquelle les cn sont des variables aléatoires gaussiennes i i d.
−→ Les travaux de Wiener (1923) marquent le démarrage de la théorie des

processus stochastiques. La théorie des probabilités acquiert un véritable statut
mathématique, le mouvement brownien devient un honnête homme.
−→ Travaux de Paul Lévy
5.1.3 Théorie de la spéculation

Bachelier soutient en 1900 une thèse de doctorat intitulée ” Théorie de la
spéculation” dans laquelle il présente une application du calcul des probabi-
lités aux opérations boursières. Il montre que, dans un marché idéal, la valeur
d’une action à l’instant t peut être décrite comme un processus à accroisse-
ment indépendants et homogènes dont les trajectoires sont continues. C’est la
première théorie mathématique du mouvement brownien
5.2 Définition d’une marche aléatoire

Considérons une particule dans l’espace euclidien Rd . Soit ~r0 sa position initiale
à l’instant t0 . A des temps ultérieurs t1 , t2 , . . . tn elle subit des sauts ~s1 , ~s2 , . . . ~sn
de sorte qu’ au temps tn sa position est
n
X
~rn = ~r0 + ~sk
k=1
~r5
~ro ~s1
Ce processus définit une marche aléatoire si les sauts successifs {~sk } sont des
variables aléatoires indépendantes distribuées selon la même loi.
p(~s)d~s = P rob {~s < ~sk < ~s + d~s}
On définit l’état du système par la position ~r de la particule au bout de n sauts.

Problème : calculer la probabilité de la position finale de la particule au bout
de n sauts.
Pn (~r)d~r = P rob {~r < ~rn < ~r + d~r}
Les déplacements ~sk étant indépendants, la position ~rn ne dépend pas de

toute l’histoire du processus mais seulement de la position de la particule lors
du dernier saut. On a ~rn = ~rn−1 + ~sn
Ceci constitue un exemple de processus de Markov.
L’état du système peut être soit continu (comme dans l’exemple ci dessus), soit
discret. C’est notamment le cas des marches aléatoires sur réseau que nous al-
lons maintenant discuter.
5.3 Marche aléatoire sur réseau

On considère un réseau hypercubique de dimension d consistant de sites repérés
par un vecteur ~r = (r1 , r2 , . . . rd ) ∈ Zd . Chaque site est relié à ses plus proches
voisins par un lien
On considère une particule partant du point ~r0 . A chaque itération, la par-

ticule saute d’un site ~l au site ~r = ~l + ~s avec une probabilité
p(~s) = p(~r − ~l)
Désignant par Pn (~r|~r0 ) la probabilité que la particule soit au site ~r au bout de n

sauts sachant qu’elle est partie du point ~r0 , la propriété de Markov nous donne
la récurrence
X
Pn+1 (~r|~r0 ) = p(~r − ~l) Pn (~l|~r0 )
~l
avec
P0 (~r|~r0 ) = δ~r,~r0
normalisation
X
Pn (~r|~r0 ) = 1
~
r
5.3. MARCHE ALÉATOIRE SUR RÉSEAU 55
Considérons pour commencer le cas d’un réseau unidimensionnel

+∞
X
Pn+1 (r|r0 ) = p(r − l) Pn (l|r0 )
l=−∞
P0 (r| r0 ) = δr,r0
Il est commode de faire une analyse de Fourier de cette récurrence. Introduisons
la série de Fourier
+∞
X
λ (θ) = p(r)eirθ
r=−∞
appelée fonction de structure
P
p(r) = 1
⇒ λ(0) = 1
La fonction génératrice
+∞
X
Γn (θ, r0 ) = Pn (r|r0 )eirθ
r=−∞
satisfait la récurrence
+∞
X
Γn+1 (θ, r0 ) = Pn+1 (r|r0 )eirθ
r=−∞
X
= p(r − l) Pn (l|r0 ) eirθ
r,l
X
= p(r − l)Pn (l|r0 ) ei(r−l)θ eilθ
r,l
= Γn (θ, r0 ) λ (θ)
La condition initiale P0 (r|r0 ) = δr,r0 donne Γ0 (θ, r0 ) = eir0 θ donc
Γn (θ, r0 ) = λn (θ)eir0 θ
d’où
Z π
1
Pn (r|r0 ) = exp −irθ Γn (θ, r0 )dθ
2π −π
Z π
1
= exp −i(r − r0 )θ λn (θ)dθ
2π −π
Puisque Pn (r|r0 ) ne dépend que de r − r0 on a
Pn (r|r0 ) = Pn (r − r0 |0) = Pn (r − r0 )
Z π
1
Pn (r) = e−irθ λn (θ)dθ
2π −π
Cette formule se généralise trivialement au cas d’un réseau Zd

Z in
1 ~
h
~
Pn (~r) = e −i~
rθ
λ(θ) dθ~
(2π)d
Exemples
1. Marche aléatoire symétrique sur Z
1
p(1) = p(−1) =
2
1 iθ
e + e−iθ = cos θ

λ(θ) =
2
Z π
1
Pn (r) = (cos θ)n e−irθ dθ
2π −π
Z π
1
= (cos θ)n cos rθ dθ
2π −π
On retrouve la distribution binomiale

n!
Pn (r) = 2−n n+r
n−r
2 ! 2 !
2. Marche aléatoire symétrique sur Z2

1
p(1, 0) = p(0, 1) = p(0, −1) = p(−1, 0) =
4
1 iθ1
e + e−iθ1 + eiθ2 + e−iθ2

λ(θ1 , θ2 ) =
4
1
= [cosθ1 + cosθ2 ]
2
5.3.1 Autre dérivation

La relation
~rn = ~r0 + ~s1 + ~s2 + . . . + ~sn
permet de calculer la fonction caractéristique de la variable aléatoire ~r = ~rn −~r0
en l’exprimant comme une somme de variables aléatoires indépendantes
~ n
h i
~ ~
E eiθ~r = E eiθ(s~1 +...+~sn ) = eiθ~s
h in
= λ(θ)~
Donc
Z h in
~ ~
eiθ~r Pn (~r)d~r = λ(θ)
Z in
1 ~
h
~
Pn (~r) = d
e−iθ~r λ(θ) dθ~
(2π)
5.4. LIMITE CONTINUE 57
5.3.2 Programme de simulation

Voir le site suivant : www.math.utah.edu
Questions typiques que l’on aimerait résoudre :
• quelle est la probabilité que le marcheur revienne à son point de départ.
• combien de sites distincts visite t-il en un temps N ?
• est-ce qu’un chemin typique s’intersecte ?
• distribution de probabilité de l’éloignement maximum ?
• temps moyen passé en un point donné ?
• influence d’une perturbation ?
5.4 Limite continue

L’idée de faire un zoom temporel et spatial d’une marche aléatoire discrète
permet de définir une limite continue non triviale. Cette idée est développée en
détail dans les cours de physique statistique (liens entre mécanique statistique et
théorie des champs). L’intérêt principal de cette limite continue est de montrer
le caractère universel du processus obtenu - indépendance par rapport au réseau
hexagonal, cubique... Ce caractère universel est bien entendu une conséquence
du théorème limite central.
Dans le cas particulier d’une marche aléatoire sur Z, nous venons de montrer

n
E eiθ(rn −r0 ) = [λ(θ)]
Supposons pour simplifier que la marche aléatoire soit issue de 0.
E eiθrn = λ(θ)n

Considérons une longue marche aléatoire de longueur n issue de 0, elle est définie
par les positions successives de la particule r1 , r2 , . . . rn
n
1 2 3 4 5 6 7
Effectuons un changement d’échelle, à la fois sur n pour se ramener à l’inter-

valle [01], et en espace pour que les points restent à distance finie quand n → ∞.
On pose
r[nt]
xn (t) ≡ √
n
où t prend des valeurs discrètes n1 , n2 , . . . 1. Il est sous entendu que l’on peut
étendre la fonction xn (t) pour tout t ∈ [0, 1] par interpolation linéaire.
Exemple : n = 7 devient
xn (t)
t
1
n = 100 aura l’allure suivante
xn (t)
t
1
Etude du processus aléatoire xn (t)
[nt]

iθxn (t)
θ
iθ
√ r

E e =E e = λ √ n [nt]
n
Exemple d’une marche aléatoire symétrique sur Z
1
p(1) = p(−1) =
2
λ(θ) = cos θ
nt
θ2 1 θ2 t
lim E eiθxn (t) = 1 − +O 2
= e− 2
n→∞ 2n n
Par conséquent xn (t) converge vers une variable aléatoire X(t) de loi
1 − x2
P (x, t) = √ e 2t
2πt
5.5. NOTIONS DE PROCESSUS ALÉATOIRE 59
Preuve :
Z +∞
1 x2 θ2 t
√ eiθx− 2t dx = e− 2
2πt −∞
5.4.1 Théorème
Le processus discret xn (t) converge en loi lorsque n → ∞ vers un processus
aléatoire continu X(t) caractérisé par les trois propriétés
• X(0) = 0
• ∀ 0 < t1 < t2 < t3 < . . . < tn l’accroissement X(tn ) − X(tn−1 ) est
indépendant de X(t1 ), X(t2 ) . . . X(tn−1 ).
• ∀ t2 > t1 l’accroissement Y = X(t2 )−X(t1 ) suit une loi gaussienne centrée
de variance t = t2 − t1 et de densité
1 y2
P (y, t) = √ exp −
2πt 2t
On en déduit que la distribution de probabilité de la loi jointe

{X(t1 ), . . . X(tn )} où 0 < t1 < t2 . . . < tn
est donnée par
1 x2 1 (x2 − x1 )2
√ exp − 1 p exp − ...
2πt1 2t1 2π(t2 − t1 2(t2 − t1 )
Dans la suite nous réintroduisons la constante de diffusion en posant

1 x2
P (x, t) = √ exp −
4πDt 4Dt
on vérifie que P est solution de l’équation de diffusion
∂P ∂2P
=D 2
∂t ∂x
avec condition initiale
P (x, t) = δ(x)
t→0
5.5 Notions de processus aléatoire

Un processus est défini par une famille infinie de variables aléatoires {X(t)}
indexées par une variable temporelle t ∈ R. On le caractérise par une hiérarchie
de distributions jointes
{X(t1 ), X(t2 ), . . . X(tn )} t1 < t2 < . . . < tn
P (x1 < X(t1 ) < x1 + dx1 , . . . xn < X(tn ) < xn + dxn ) ≡
Pn (x1 t1 , x2 t2 , . . . xn tn )dx1 dx2 . . . dxn
Elles satisfont manifestement les conditions
1. positivité Pn > 0
R
2. normalisation Pn dx1 . . . dxn = 1
3. compatibilité
Z
Pn dxj = Pn−1 (x1 t1 , . . . , xj−1 tj−1 , xj+1 tj+1 , . . . xn tn )
5.5.1 Probabilités conditionnelles

On définit la probabilité conditionnelle que X(t2 ) = x2 à dx2 près sachant que
X(t1 ) = x1 .
P2 (x1 t1 , x2 t2 )
P (x2 t2 | x1 t1 ) =
P1 (x1 t1 )
en intégrant sur x2
Z Z
P2 (x1 t1 , x2 t2 )dx2 = P1 (x1 t1 ) P (x2 t2 | x1 t1 )dx2
= P1 (x1 t1 ) donc
Z
P (x2 t2 | x1 t1 )dx2 = 1
Généralisation
Pn (x1 t1 , . . . xn tn )
P (xk+1 tk+1 , . . . , xn tn | x1 t1 . . . xk tk ) =
Pk (x1 t1 , . . . xk tk )
5.5.2 Processus de Markov

Il est caractérisé par la condition ∀ t1 < t2 .. < tn
P (xn tn | x1 t1 , . . . xn−1 tn−1 ) = P (xn tn | xn−1 tn−1 )
P est appelée probabilité de transition.

Un processus de Markov est donc entièrement caractérisé par les deux fonctions
P1 (x1 t1 ) et P (x2 t2 | x1 t1 )
Il suffit donc de connaı̂tre la loi de probabilité initiale et la probabilité de tran-

sition.
Exemple
P3 (x1 t1 , x2 t2 , x3 t3 ) = P2 (x1 t1 , x2 t2 )P (x3 t3 | x1 t1 , x2 t2 )

= P (x2 t2 | x1 t1 )P1 (x1 t1 )P (x3 t3 | x1 t1 , x2 t2 )
= P (x2 t2 | x1 t1 )P1 (x1 t1 )P (x3 t3 | x2 t2 )
En poursuivant cet algorithme on peut construire toutes les Pn .

5.5. NOTIONS DE PROCESSUS ALÉATOIRE 61
5.5.3 Equation de Chapman-Kolmogorov

Intégrons l’équation précédente sur x2 , on obtient
Z
P2 (x1 t1 , x3 t3 ) = P1 (x1 t1 ) P (x2 t2 | x1 t1 )P (x3 t3 | x2 t2 )dx2
En divisant par P1 (x1 t1 ) on obtient l’équation de Chapman-Kolmogorov
Z
P (x3 t3 | x1 t1 ) = P (x3 t3 | x2 t2 )P (x2 t2 | x1 t1 )dx2
Par conséquent pour se donner un processus de Markov, il faut se donner

deux fonctions
P1 (x1 t1 ) et P (x2 t2 | x1 t1 )
Celles-ci ne sont pas indépendantes

R
1. P1 (x2 t2 ) = P (x2 t2 | x1 t1 )P1 (x1 t1 )dx1
2. P doit satisfaire Chapman-Kolmogorov
5.5.4 Processus aléatoire stationnaire

Un processus aléatoire est stationnaire si sa loi temporelle est invariante par
translation :
Pn (x1 t1 , . . . , xn tn ) = Pn (x1 t1 + τ, . . . , xn tn + τ )
En particulier
P1 (x1 t1 ) = P1 (x1 ) est indépendant du temps
5.5.5 Théorème de Doob

”Un processus aléatoire gaussien stationnaire est markovien si et seulement si
sa fonction d’autocorrélation est exponentielle”
• Exprimons que X(t) est un processus stationnaire. La probabilité de tran-
sition P (x2 t2 | x1 t1 ) ne dépend que de t2 −t1 . On pose P (x2 t1 +t | x1 t1 ) =
Pt (x2 | x1 ). Elle satisfait l’équation de CK
Z
Pt+t0 (x3 | x1 ) = Pt0 (x3 | x2 )Pt (x2 | x1 )dx2
• Exprimons que X(t) est un processus gaussien. Les probabilités Pn sont

donc de la forme
1
Pn ∼ N exp − xi xj Aij (t1 , . . . tn )
2
Par conséquent
1
Pt (x2 | x1 ) = D exp − (A x21 + 2B x1 x2 + C x22 )
2
où A, B, C, D sont des fonctions de t. En particulier, en utilisant la sta-

tionnarité
1 x2
P1 (x) = √ exp −
2π 2
La normalisation
r
B2
Z
A
Pt (x2 | x1 )dx2 = 1 donne D = , C=
2π A
L’équation de consistance
Z
P1 (x2 ) = Pt (x2 | x1 )P1 (x1 )dx1 donne B 2 = A(A − 1)
• Exprimons maintenant la fonction de corrélation

Z
E(x(t1 )x(t2 )) = x1 x2 P2 (x1 t1 , x2 t2 )dx1 dx2
Z
= x1 x2 P1 (x1 t1 )P (x2 t2 | x1 t1 )dx1 dx2
Z
= x1 x2 P1 (x1 )Pt2 −t1 (x2 | x1 )dx1 dx2
Posons C(t) = E(x(0) x(t))

Il vient
Z
C(t) = x1 x2 P1 (x1 ) Pt (x2 | x1 )dx1 dx2
avec A = (1 − C 2 )−1
Par conséquent
1 (x2 − Cx1 )2
Pt (x2 | x1 ) = p exp −
2π(1 − C 2) 2(1 − C 2 )
En utilisant l’équation CK
Z
Pt+t0 (x3 | x1 ) = Pt0 (x3 | x2 )Pt (x2 | x1 )dx2
on obtient l’équation fonctionnelle
C(t + t0 ) = C(t)C(t0 ) t, t0 > 0
dont la solution est C(t) = exp −γt. On obtient ainsi le processus d’Ornstein-
Uhlenbeck
1 (x2 − e−γt x1 )2
Pt (x2 | x1 ) = p exp −
2π(1 − exp −2γt) 2(1 − e−2γt )
5.6. EXEMPLES DE PROCESSUS STOCHASTIQUES GAUSSIENS 63
5.6 Exemples de processus stochastiques gaus-

siens
5.6.1 Processus de Wiener
Le mouvement brownien issu de l’origine X(0) = 0, appelé processus de Wiener
est caractérisé par
1 1 (x2 − x1 )2
P (x2 t2 | x1 t1 ) = p exp −
4πD(t2 − t1 ) 4D t2 − t1
1 x2
P1 (x1 t1 ) = √ exp − 1
4πDt1 4Dt1
P1 (x1 t1 ) dépend de t1 ⇒ il s’agit donc d’un processus non stationnaire.
Calculons la fonction de corrélation notée indifféremment < x(t2 )x(t1 ) > ou
E(x(t2 )x(t1 ))
Z
E(x(t2 )x(t1 )) = dx2 dx1 x2 x1 P2 (x2 t2 , x1 t1 )
Z
= dx2 dx1 x2 x1 P (x2 t2 | x1 t1 )P1 (x1 t1 )
Z
1 1
= dx2 dx1 x2 x1 p √
4πD(t2 − t1 ) 4πDt1
1 (x2 − x1 )2 x2
exp − − 1
4D t2 − t1 4Dt1
Il s’agit d’un processus gaussien décrit par la matrice
t2 1
− t2 −t

1 t1 (t2 −t1 )
A= 1 1
1
2D − t2 −t 1 t2 −t1
dont l’inverse

−1 t1 t1
A = 2D
t1 t2
nous permet de calculer la covariance.

Pour t2 > t1 on obtient E(x(t2 )x(t1 )) = 2Dt1 .
La formule générale s’écrit
E(x(t2 )x(t1 )) = 2D min (t1 , t2 )
1. On vérifie que le processus n’est pas stationnaire puisque sa covariance ne

dépend pas seulement de la différence des temps.
2. Montrons que les trajectoires ne sont pas différentiables.

Pour cela considérons
Z ∞
x2

| X(t + h) − X(t) | 2
P rob >k = √ exp − dx
h 4πDh kh 4Dh
Z ∞
1 y2
=√ √ exp − dy
πD k h 4D
∀k > 0, pour tout h → 0 cette quantité tend vers 1. Par conséquent ∀k, la
différence |X(t+h)−X(t)|
h > k avec probabilité 1.
5.6.2 Processus d’Ornstein-Uhlenbeck

1 (x2 − x1 e−γ(t2 −t1 ) )2
P (x2 t2 | x1 t1 ) = p exp −
2π(1 − e−2γ(t2 −t1 ) ) 2(1 − e−2γ(t2 −t1 ) )
1 x2
P1 (x1 , t1 ) = P1 (x1 ) = √ exp − 1
2π 2
Fonction de corrélation
E(x(t2 )x(t1 )) = e−γ(t2 −t1 )

1
γ est le temps de corrélation.
5.6.3 Bruit blanc

C’est une idéalisation d’un processus dont le temps de corrélation tend vers 0.
E(η(t2 )η(t1 )) = δ(t2 − t1 )
Montrons comment relier formellement le bruit blanc au processus de Wiener

W (t).
√ Z t
W (t) = 2D η(τ )dτ
0
satisfait
W (0) = 0
Z t0 Z t
E(W (t)W (t0 )) = 2D dτ 0 dτ δ (τ − τ 0 )
o o
Z t0
= 2D dτ 0 θ(t − τ 0 )
o
Z min(t,t0 )
= 2D dτ 0 = 2D min (t, t0 )
o
5.6. EXEMPLES DE PROCESSUS STOCHASTIQUES GAUSSIENS 65
W (t) est donc un processus de Wiener. Montrons de même comment construire

le processus d’Ornstein-Uhlenbeck Y (t) à partir d’un bruit blanc.
Z t
Y (t) = dτ e−γ(t−τ ) η(τ )
−∞
satisfait
Z t Z t0
0 0
E(Y (t) Y (t0 )) = dτ e−γ(t−τ ) dτ 0 e−γ(t −τ ) δ(τ − τ 0 )
−∞ −∞
Z t
0
= dτ e−γ(t−τ ) dτ e−γ(t −τ ) θ(t0 − τ )
−∞
Z min(t,t0 )
0
= e−γ(t+t ) dτ e2γτ
−∞
−γ(t+t0 ) 1 2γ min(t,t0 )
=e e
2γ
or t + t0 − 2 min(t, t0 ) = |t − t0 |
0
donc E(Y (t) Y (t0 )) = 2γ1
e−γ|t−t | .
Par conséquent Y (t) est le processus d’Ornstein-Uhlenbeck discuté précédemment.
Remarque : la condition initiale est renvoyée à l’infini Y (−∞) = 0.
5.6.4 Construction de W (t) comme processus

Considérons l’ensemble des fonctions f ∈ L2 [0, 1] défini par
Z 1
f 2 (t) dt < ∞
0
Toute fonction f ∈ L2 peut être décomposée sur une base {fn } telle que
Z 1
fn (t) fm (t)dt = δnm
0
X
f (t) = cn fn (t)
Z 1
où cn = fn (t)f (t)
0
Supposons que les cn soient des variables aléatoires gaussiennes iid
E(cn ) = 0 E(cn cm ) = δnm
Considérons
X Z t
W (t) = cn fn (s)ds
n 0
Montrons que W (t) est un brownien issu de l’origine.

Preuve :
X Z t Z t0
0
E [W (t) W (t )] = E(cn cm ) fn (s)ds fm (u)du
n,m 0 0
XZ t Z t0
= fn (s)ds fn (u)du
n 0 0
En utilisant la relation de fermeture

X
fn (s)fn (u) = δ(s − u)
Il vient
Z t Z t0
E [W (t)W (t0 )] = ds duδ(s − u)
0 0
Z t Z t0
= ds δ(s − u)du
0 0
Z t
= dsθ(t0 − s) = min(t, t0 )
0
E [W (t) W (t0 )] = min(t, t0 )

Chapitre 6
Intégrale de Wiener et
formule de Feynman Kac
Sommaire
6.1 Equation de la diffusion, formalisme opératoriel . 67
6.2 Mesure de Wiener . . . . . . . . . . . . . . . . . . . 68
6.3 Mouvement brownien avec absorption-formule de
Feynman-Kac . . . . . . . . . . . . . . . . . . . . . . 70
6.4 Formule de Feynman-Kac . . . . . . . . . . . . . . 71
6.1 Equation de la diffusion, formalisme opératoriel

Nous avons montré dans le cours précédent que le processus de Wiener issu de
l’origine est caractérisé par sa probabilité de transition.
1 (x − x0 )2
P (x t|x0 t0 ) = p exp −
4πD(t − t0 ) 4D(t − t0 )
1 x2
et P1 (x t) = √ exp −
4πDt 4Dt
On vérifie facilement que P (xt|x0 t0 ) est la solution de l’équation aux dérivées
partielles
∂P ∂2P
=D 2
∂t ∂x
qui satisfait la condition initiale
lim P (xt|x0 t0 ) = δ(x − x0 ).
t→t0
Ceci nous suggère d’introduire un formalisme opératoriel. Posons

∂2
Ho = −D
∂x2
67
68CHAPITRE 6. INTÉGRALE DE WIENER ET FORMULE DE FEYNMAN KAC
Considérons l’équation d’évolution
∂
Uo (t, to ) = −Ho Uo (t, to )
∂t
dont la solution s’écrit Uo (t, t0 ) = e−Ho (t−t0 t > t0 .

La probabilité de transition P (xt|x0 t0 ) peut s’écrire comme le noyau de
l’opérateur U
P (xt|x0 t0 ) =< x|e−Ho (t−t0 ) |x0 >
où nous avons adopté la notation de Dirac.

L’équation de Chapman Kolmogorov
Z
dx1 P (xt|x1 t1 )P (x1 t1 |x0 t0 ) = P (xt|x0 t0 )
s’écrit
Z
dx1 < x|e−Ho (t−t1 ) |x1 >< x1 |e−Ho (t1 −t0 ) |x0 >=< x|e−Ho (t−t0 ) |x0 >
En utilisant la relation de ”fermeture” nous voyons que l’équation de Chapman-

Kolmogorov est équivalente à la relation opératorielle
e−Ho (t−t1 ) e−Ho (t1 −t0 ) = e−Ho (t−t0 )
Les opérateurs e−tH définissent une loi de semi-groupe
6.2 Mesure de Wiener

(voir I.M. Gelfand et A.M. Yaglom, Journal of Mathematical Physics, volume
1, 1960 p 48-69 et Ph. Martin, initiation à l’intégrale fonctionnelle)
Puisque la trajectoire x(τ ) d’une particule brownienne est aléatoire, il est
naturel de se poser les questions suivantes
• Comment définir la probabilité d’un chemin ?
• Comment calculer la valeur moyenne d’une fonctionnelle F [x(τ )] sur les

chemins browniens ?
Considérons l’ensemble des trajectoires issues de x0 en τ0 et astreintes à traverser

les portes [x1 , x1 + dx1 ], . . . [xn , xn + dxn ] aux temps respectifs τ1 , . . . , τn
x1 < x(τ1 ) < x1 + dx1
xn < x(τn ) < xn + dxn

6.2. MESURE DE WIENER 69
x(t)
Xo
τ0 τ1 τ2 τ3 t
cf confirmation d’un polymère
Pour calculer la probabilité de cet événement nous pouvons utiliser le ca-

ractère Markovien du processus
P (x0 τ0 , x1 τ1 , . . . xn τn )
P (x1 τ1 , . . . xn τn |x0 τ0 ) =
P1 (x0 τ0 )
= P (xn τn |xn−1 τn−1 ) . . . . . . P (x1 τ1 |x0 τ0 )
On peut ainsi définir une mesure attachée à cette classe particulière de che-
mins issus de x0 .
dµ [x(τ )] ≡ P [x1 τ1 , . . . xn τn |x0 τ0 ] dx1 . . . dxn
Nous aimerions donner un sens à la valeur moyenne de la fonctionnelle F [x(τ )]

Z
F [x(τ )]dµ[x(τ )]
Pour celà considérons une fonctionnelle qui ne dépend du chemin que par l’in-
termédiaire d’un nombre fini de points
F [x(τ )] ≡ F [x(τ1 ), x(τn−1 ), . . . x(τn )]
La valeur moyenne de cette fonctionnelle s’écrit

Z Z
d µ[x(τ )]F [x(τ )] = dx1 . . . dxn P [x1 τ1 , . . . xn τn |x0 τ0 ]F [x1 , x2 , . . . xn ]
x0 ,τ0
La mesure de probabilité attachée aux chemins browniens issus de x0 s’appelle

mesure de Wiener.
Dans le cas où le point final x(τ ) = x est prescrit on introduit la mesure de
Wiener conditionnelle en intégrant.sur tous les points sauf le dernier.
Z xτ Z
dµ[x(σ)]F [x(σ)] = dx1 . . . dxn P [x1 τ1 , x2 τ2 . . . xn τn , xτ |x0 τ0 ]F [x1 , . . . xn ]
x0 ,τ0
Z
= dx1 . . . dxn P [xτ |xn τn ] . . . . . . P [x1 τ1 |x0 τ0 ]F [x1 , . . . xn ]
On a la relation
Z Z Z xτ
dµ[x(τ )]F [x(τ )] = dx dµ[x(τ )]F [x(τ )]
x 0 τ0 x0 τ0
Exemple : considérons les valeurs moyennes de browniens issus de l’origine

Z Z
E[x(τ1 )] = dµ[x(τ )]x(τ1 ) = dx1 x1 P (x1 τ1 |00) = 0
0,0
Z
E[x(τ1 )x(τ2 )] = dµ[x(τ )]x(τ1 )x(τ2 )
0,0
Z
= dx1 dx2 P (x2 τ2 | x1 τ1 )P (x1 τ1 | 00)x1 x2
d’après le calcul du chapitre précédent on a
E[x(τ1 )x(τ2 )] = 2D min(τ1 , τ2 )
6.3 Mouvement brownien avec absorption-formule

de Feynman-Kac
Nous nous proposons de calculer
Rτ l’espérance par rapport à la mesure de Wiener
de la fonctionnelle exp − τ0 dσΩ [x(σ)]. Pour cela il est commode d’interpréter
Ω(x) comme une probabilité d’absorption au point x. Supposons que la particule
ait une probabilité par unité de temps Ω(x) > 0 d’être absorbée au point x.
−τ0
Divisons l’intervalle [τ0 , τ ] en n + 1 sous intervalles de longueur = τn+1 .
Posons
τk = τ0 + k k = 0, . . . n, τn+1 = τ
La probabilité que la particule ne soit pas absorbée le long d’une trajectoire
donnée x(σ) τ0 6 σ 6 τ est
n
Y n
X Z τ
lim (1−Ω(x(τk ))) = lim exp − Ω(x(τk )) = exp − dσΩ [x(σ)]
n→∞ n→∞ τ0
k=0 k=0
Par conséquent la probabilité PΩ (xτ | x0 τ0 ) de trouver la particule en x au

temps τ est
Z xτ Z τ
PΩ (xτ | x0 τ0 ) = dµ [x(σ)] exp − dσΩ [x(σ)]
x0 τ0 τ0
Z xτ n
X
= lim dµ exp − Ω(x(τk ))
n→∞ x 0 τ0 k=0
6.4. FORMULE DE FEYNMAN-KAC 71
Pour calculer PΩ revenons à l’approximation discrète. Le poids de probabilité

attaché à la subdivision précédente est
P (xτ | xn τn ) . . . . . . P (x2 τ2 | x1 τ1 )P (x1 τ1 | x0 τ0 )dx1 . . . dxn
n
n+1 Y n
1 1 X
= √ dxk exp − (xk+1 − xk )2
4πD 4D
k=1 k=0
On peut réécrire cette expression sous une forme suggestive ”à la Feynman”
n
n+1 Y n
1 X xk+1 − xk 2
= √ dxk exp − ( )
4πD 4D 0
k=1
Z τ 2
1 dx
= Dx(σ) exp − dσ
4D τ0 dσ
Par conséquent nous pouvons écrire PΩ sous deux formes équivalentes
Z xτ Z τ ( 2 )
1 dx
PΩ (xτ | xo τo ) = Dx(σ) exp − dσ + Ω [x(σ)]
xo τo τo 4D dσ
Z xτ Z τ
= dµ [x(σ)] exp − dσΩ [x(σ)]
x0 τ0 τ0
6.4 Formule de Feynman-Kac

Montrons que PΩ est la solution fondamentale de l’équation aux dérivées par-
tielles
∂ ∂2
PΩ (xτ | xo τo ) = D 2 PΩ (xτ | x0 τ0 ) − Ω(x)PΩ (xτ | xo τo )
∂τ ∂x
lim PΩ (xτ | x0 τ0 ) = δ(x − x0 )
τ →τ0
Démonstration : revenons à l’expression discrète
Z xτ n
X
dµ[x(τ )] exp − Ω(x(τk )) =
x 0 τ0 k=0
Z
dx1 . . . dxn < x | e−H0 | xn > e−Ω(xn )
< xn | e−H0 | xn−1 > e−Ω(xn−1 )

< x1 | e−H0 | x0 > e−Ω(x0 )
Z
= dx1 . . . dxn < x | e−Ho e−Ω | xn >
< xn | e−Ho e−Ω | xn−1 >

< x1 | e−Ho e−Ω | x0 >
où Ω représente maintenant un opérateur qui agit multiplicativement. En utili-

sant la relation de fermeture on obtient
< x | (e−Ho e−Ω )n+1 | x0 >

n+1
(τ − τo )Ho (τ − τ0 )Ω
= lim < x | exp − exp − |x0 >
n→∞ n+1 n+1
La formule de Trotter
A B
lim (exp exp )n = exp(A + B)
n→∞ n n
donne
PΩ (xτ |x0 τ0 ) =< x| exp −(τ − τ0 )(Ho + Ω)|x0 >

=< x|U (τ, τ0 )|x0 >
On fait ainsi apparaitre le propagateur de Feynman
< x| exp −(τ − τ0 )(Ho + Ω)|x0 >
qui satisfait l’équation aux dérivées partielles.

∂2

∂
PΩ (xτ |x0 τ0 ) = D 2 − Ω(x) PΩ (xτ | x0 τ0 )
∂τ ∂x
ou encore
∂
− PΩ (xτ | x0 τ0 ) = HPΩ (xτ | xo τo )
∂τ
2
∂
avec H = −D ∂x 2 + Ω(x).
Ω(x) qui était défini comme une probabilité par unité de temps pour la
particule d’être absorbée s’interprète donc comme un potentiel en mécanique
quantique.
Résumé
Z xτ Z τ
P (xτ | x0 τ0 ) = dµ [x(σ)] exp − dσΩ [x(σ)]
x 0 τ0 τ0
est solution de l’équation

∂P ∂2
= D 2 P − Ω(x)P
∂t ∂x
Chapitre 7
Équation de Langevin et
Fokker-Planck
Sommaire
7.1 Introduction . . . . . . . . . . . . . . . . . . . . . . 73
7.2 Représentation fonctionnelle des mesures gaussiennes 75
7.3 Lemme . . . . . . . . . . . . . . . . . . . . . . . . . . 76
7.4 Equation de Langevin . . . . . . . . . . . . . . . . . 76
7.5 Equation de Fokker-Planck . . . . . . . . . . . . . . 77
7.6 Solution d’équilibre . . . . . . . . . . . . . . . . . . 79
7.7 Relation avec l’équation de Schrödinger en temps
imaginaire . . . . . . . . . . . . . . . . . . . . . . . . 80
7.8 Propriété du hamiltonien H+ et supersymétrie . . 81
7.9 Application . . . . . . . . . . . . . . . . . . . . . . . 82
7.9.1 Processus de Wiener . . . . . . . . . . . . . . . . . . 82
7.9.2 Processus d’Ornstein-Uhlenbeck . . . . . . . . . . . 83
7.9.3 Processus bistables . . . . . . . . . . . . . . . . . . . 84
7.10 Représentation fonctionnelle de P (x t|y 0) . . . . . 85
7.11 Limite D → O, formule d’Arrhénius . . . . . . . . . 86
7.1 Introduction
De nombreux phénomènes peuvent être décrits de façon phénoménologique par
l’ équation de Langevin
dxi
= φi (x) + ηi (t) i = 1...N
dt
où les xi sont des variables macroscopiques représentant par exemple une densité
de charge, de matière, les φi sont des forces de dérive telles que des gradients
de température, de concentration, des champs électriques ou gravitationnels.
ηi (t) est un bruit gaussien
< ηi (t)ηj (t0 ) >= 2 δij δ(t − t0 )
73
74 CHAPITRE 7. ÉQUATION DE LANGEVIN ET FOKKER-PLANCK
1. P. Langevin, C.R. Acad.Sci 146, 530 (1908)

Mouvement d’une particule dans un fluide
~
dV
m ~ + ~η (t)
= −γ V
dt
2. Circuit électrique
L dI
dt
= −RI + V (t)
< V (t)V (t0 ) >= Γδ(t − t0 )
3. Grande diversité de phénomènes décrits dans ce cadre (systèmes déterministes

perturbés par un bruit)
Nous avons déjà rencontré l’équation de Langevin sous forme intégrée dans
deux cas particuliers
Rt√
1. x(t) = 0 2Dη(τ )dτ est un processus dont la fonction de corrélation
< x(t) x(t0 ) >= 2D min(t − t0 )
est celle du processus de Wiener qui satisfait l’équation différentielle sto-

chastique : √
dx(t) = 2Dη(t)dt
Probabilité de transition :
1 (x − x0 )2
P (x t|x0 0) = √ exp −
4πDt 4Dt
√ Rt
2. x(t) = 2γD −∞
d τ e−γ(t−τ ) η(τ )d τ
a pour fonction de corrélation

0
< x(t)x(t0 ) >= D e−γ(t−t )
Ceci définit le processus d’Ornstein Uhlenbeck de probabilité de transition.
1 (x − x0 e−γt )2
P (x t|x0 t0 ) = p exp −
2πD(1 − e−2γt ) 2D(1 − e−2γt )
Le processus est solution de l’équation différentielle stochastique

p
dx(t) = 2γD η(t)dt − γx(t)dt
ou encore
p
ẋ(t) = −γx(t) + 2γD η(t)
7.2. REPRÉSENTATION FONCTIONNELLE DES MESURES GAUSSIENNES75
7.2 Représentation fonctionnelle des mesures gaus-

siennes
La loi de distribution multinormale
√
detA 1
P (x) = n exp − x̃A x
(2π) 2 2
est caractérisée par une matrice de covariance
< xi xj >= A−1

ij
On a donc
Aik < xk xl >= δil
A un processus Gaussien de densité

Z
1
P (x(τ )) ∼ exp − x(τ )A(τ, τ 0 )x(τ 0 )dτ dτ 0
2
correspond une fonction de corrélation < x(τ )x(τ 0 ) >= C(τ, τ 0 ) telle que
Z
A(τ, τ 00 )C(τ 00 , τ 0 )dτ 00 = δ(τ − τ 0 )
Exemples
1) Le processus de Wiener est caractérisé par
1 d2
A(τ, τ 0 ) = − δ(τ − τ 0 )
2D dτ 02
Vérification : la fonction de corrélation
C(τ, τ 0 ) = 2Dmin(τ, τ 0 ) = D [τ + τ 0 − |τ − τ 0 |]
doit satisfaire
Z
A(τ, τ 00 )C(τ 00 , τ 0 )dτ 00 = δ(τ − τ 0 )
soit
1 d2
Z
− δ(τ − τ 00 ) C(τ 00 , τ 0 ) dτ 00 = δ(τ − τ 0 )
2D dτ 2
1 d2
− C(τ, τ 0 ) = δ(τ − τ 0 )
2D dτ 2
2) Bruit blanc
C(τ, τ 0 ) =< η(τ ) η(τ 0 ) >= 2Dδ(τ − τ 0 )
Donc
Z
A(τ, τ 00 )2Dδ(τ 00 − τ 0 )dτ 00 = δ(τ − τ 0 )
d’où
2DA(τ, τ 0 ) = δ(τ − τ 0 )
Par conséquent
Z
1
P (η) = exp − η(τ ) A(τ − τ 0 )η(τ 0 )dτ dτ 0
2
Z +∞
1
P (η) = exp − η 2 (τ )dτ
4D −∞
7.3 Lemme
X ∂f
< Xi f (X) >= < Xi Xj >< >
j
∂Xj
Le membre de droite s’ écrit d’après le résultat précédent
√ Z
−1 ∂f det A ∂f 1
(A )ij < >= n d~x (A−1 )ij exp − x̃ A x
∂Xj (2π) 2 ∂xj 2
Intégrons par parties
√ Z
det A 1
=+ n d~x(A−1 )ij Ajk xk exp − x̃ A x f (x)
(2π) 2 2
√ Z
det A 1
= n d~x(A−1 A)ik xk f (x) exp − x̃ A x
(2π) 2 2
√ Z
det A 1
= n d~xxi f (x) exp − x̃ A x
(2π) 2 2
=< Xi f (X) >
L’utilisation répétée de cette relation permet de montrer que les moments d’ordre
pair se factorisent en un produit de moments d’ordre deux. Les moments d’ordre
impair s’annulent trivialement.
7.4 Equation de Langevin

ẋ = φ(x) + η(t)
< η(t)η(t0 ) >= 2Dδ(t − t0 )
< η(t) >= 0
Z +∞
1
P (η)Dη = exp − η 2 (τ )dτ Dη
4D −∞
Dans la littérature mathématique l’équation s’écrit
dx = φ(x)dt + dW
en accord avec la définition du processus de Wiener
Z t
W = η(τ )dτ
o
7.5. EQUATION DE FOKKER-PLANCK 77
7.5 Equation de Fokker-Planck

Pour des conditions initiales fixées x(0) = x0 , on se propose de calculer la valeur
moyenne au temps t > 0 de g(x(t)) en effectuant une moyenne sur toutes les
réalisations du processus
x(t)
xo
t
t0
Z
< g(x(t) >= g(x) P (x t|x0 0)dx
Dans la suite on omettra la condition initiale en écrivant la probabilité de tran-

sition
P (x, t) = P (x t|x0 0)
Par définition
Z
g(x(t)) = g(x) δ (x − x(t)) dx
Par conséquent
Z
< g(x(t) >= g(x) < δ(x − x(t) > dx
On peut donc écrire la probabilité de transition sous la forme
P (x, t) =< δ(x − x(t) >
Soit
Z Z +∞
1
P (x, t) = Dη δ [x(t) − x] exp − η 2 (τ )dτ
4D −∞
Dérivons par rapport au temps

Z Z +∞
∂ 1
Ṗ (x, t) = Dη δ[x(t) − x] exp − η 2 (τ )dτ
∂t 4D −∞
Z Z
∂ 1
= Dη δ [x(t) − x] [φ(x) + η(t)] exp − η 2 (τ )dτ
∂x(t) 4D
Z Z
∂ 1
= Dη − δ [x(t) − x] [φ(x(t) + η(t)] exp − η 2 (τ )dτ
∂x 4D
Z Z
∂ 1
=− Dηδ [x(t) − x] [φ(x(t) + η(t)] exp − η 2 (τ )dτ
∂x 4D
Z Z
∂ ∂ 1
=− [φ(x)P (x, t)] − Dηδ [x(t) − x] η(t) exp − η 2 (τ )dτ
∂x ∂x 4D
∂ ∂
=− [φ(x)P (x, t)] − < η(t)δ [x(t) − x] >
∂x ∂x
Dans le cas d’une distribution multinormale on a la relation
X ∂f
< xi f (x) >= < xi xj >< >
j
∂xj
Cette relation s’étend à des processus gaussiens continus

Z
δf
< η(t)f [η] >= < η(t)η(t0 ) >< > dt0
δη(t0 )
En particulier
δ δ[x(t) − x]
Z
< η(t)δ[x(t) − x] >= < η(t)η(t0 ) >< > dt0
δη(t0 )
Z
δ
= 2D δ(t − t0 ) < δ[x(t) − x] > dt0
δη(t0 )
δ
= 2D < δ[x(t) − x] >
δη(t)
δx(t) δ
= 2D < δ[x(t) − x] >
δη(t) δx(t)
δx(t) ∂
= −2D < δ[x(t) − x] >
δη(t) ∂x
∂ δx(t)
= −2D < δ[x(t) − x] >
∂x δη(t)
Intégrons formellement l’équation de Langevin en posant
Z t Z t
x(t) = x(0) + φ[x(τ )]dτ + η(τ )dτ
0 0
La dérivée fonctionnelle s’écrit donc

Z t Z t
δx(t) δφ[x(τ )]
= dτ + δ(τ − t0 ) dτ
δη(t0 ) 0 δη(t0 ) 0
| {z }
θ(t−t0 )
La première contribution se simplifie en observant que l’équation de Langevin

est causale : x(τ ) ne dépend de η(t0 ) que si t0 < τ .
Cette contribution s’écrit donc :
Z t
δ[φ(x(τ )]
dτ
t 0 δη(t0 )
7.6. SOLUTION D’ÉQUILIBRE 79
Elle tend vers 0 pour t → t0 .

Par conséquent δx(t)
δη(t) = 0 + θ(0).
Il faut poser θ(0) = 21 . Pour nous en convaincre, prenons un bruit avec une
corrélation finie
< η(t)η(t0 ) >= 2D f (t − t0 )

R +∞
où f (t) est une fonction paire telle que −∞
f (t)dt = 1.
Les équations deviennent :
Z
δf
< η(t)f (η) >= 2D f (t − t00 ) < > dt00
δη(t00 )
Si f (η) = δ[x(t) − x] on obtient :
δ ∂ δx(t)
< δ[x(t) − x] >= − < δ[x(t) − x]
δη(t00 ) ∂x δη(t00 )
δx(t)
= θ(t − t00 ) + contribution qui s0 annule pour t → t00
δη(t00 )
Z
∂
< η(t)f (η) >= −2D f (t − t00 )θ(t − t00 )dt00 P (x, t)
∂x
Or
Z +∞ Z 0
1
f (t − t00 )θ(t − t00 )dt00 = f (t)dt =
−∞ −∞ 2
D’où le résultat final

∂
< η(t)δ[x(t) − x] >= −D P (x, t)
∂x
Par conséquent P satisfait l’équation de Fokker-Planck
∂P ∂2 ∂
= D 2 P (x, t) − [φ(x)P (x, t)]
∂t ∂x ∂x
7.6 Solution d’équilibre

Pour discuter l’approche à l’équilibre, il est commode d’introduire ”le potentiel”
dont dérive la ”force” φ(x).
∂
φ(x) = − U (x)
∂x
En dimension 1, une force dérive toujours d’un potentiel. En dimension supérieure
φi (xj ) doit satisfaire les conditions d’intégrabilité
∂φi ∂φj
=
∂xj ∂xi
Dans le cas d’une force linéaire φi = Lij xj on obtient Lij = Lji . Ce sont les
relations de réciprocité d’Onsager.
Revenons à la dimension 1, l’équation de Fokker-Planck s’écrit
∂P ∂ ∂P ∂U
= [D +P ]
∂t ∂x ∂x ∂x
∂P
Une solution d’équilibre, indépendante du temps satisfait à ∂t = 0 soit
∂P ∂U
D +P = −J
∂x ∂x
où J est une constante qui s’interprète comme le courant de diffusion.
Si J = 0, il vient Peq = exp − UD(x)
R
Cette solution d’équilibre n’existe que si Peq (x)dx < ∞
exemple : dans un champ de pesanteur U (x) = kx on a une solution normalisable
sur [a, ∞[ mais pas sur ] − ∞, +∞[ pour des raisons physiques évidentes
Des solutions avec courant peuvent intervenir dans des situations où il y a
un gradient de concentration
$P_{o} > P_{1}$
$P_{o}$ J $P_{1}$
7.7 Relation avec l’équation de Schrödinger en

temps imaginaire
Montrons comment faire apparaı̂tre un opérateur autoadjoint en posant
U (x)
P (x, t) = exp − ψ(x, t)
2D
Il vient
h i
U0
DP 0 + P U 0 = D − 2D ψ + ψ 0 exp − U2D
(x)
+ PU0
h 0
i
= + U2 ψ + Dψ 0 exp − U2D (x)
∂ h 00
U0 0
i
U0
h 0
i
(DP 0 + P U 0 ) = + U2 ψ + 2 ψ + Dψ 00 exp − U2D
(x)
− 2D + U2 ψ + Dψ 0 exp − U2D
(x)
∂x h i
00 02
= + U2 − U4D ψ + Dψ 00 exp − U2D (x)
7.8. PROPRIÉTÉ DU HAMILTONIEN H+ ET SUPERSYMÉTRIE 81
Soit
U 02 U 00

∂ψ
= −Dψ 00 + − ψ = H+ ψ
∂t 4D 2
∂P
− = LP
∂t
∂ψ
− = H+ ψ
∂t
∂2 U 02 U 00
H+ = −D + −
∂x2 4D 2
U U
L = exp − H+ exp
2D 2D
7.8 Propriété du hamiltonien H+ et supersymétrie

Exprimons H+ en terme des opérateurs différentiels du premier ordre
√ ∂ 1 ∂U
Q= D + √
∂x 2 D ∂x
√ ∂ 1 ∂U
Q+ = − D + √
∂x 2 D ∂x
On vérifie que H+ = Q+ Q
Conséquences :
1. H+ est défini positif
< n|Q+ Q|n >=k Q|n >k2 > 0
2. fondamental Q|0 >= 0
U (x)
ψ0 (x) = exp − est un état propre d’énergie nulle associé à Peq (x).
2D
Exprimons la probabilité de transition solution de

∂P ∂ ∂P ∂U
∂t = ∂x [D ∂x + P ∂x ]
limt→o P (x t|y 0) = δ(x − y)
en terme des fonctions propres de H+
ψ0 (x) X
P (x t|y 0) = ψn (x)ψn (y)e−En t
ψ0 (y) n
ψ0 (x)
Comportements limites t → ∞ P (x t|y 0) → ψ0 (y) ψ0 (x)ψ0 (y) = ψ02 (x) = exp − UD
(x)
ψ0 (x)
t→0 P (x t|y 0) → ψ0 (y) δ(x − y) = δ(x − y)
Approche à l’équilibre : elle est décrite par le spectre de relaxation constitué

de l’ensemble des En (n > 1). En pratique, on s’intéresse essentiellement à la
première valeur propre non nulle E1 .
On pose :
" 2 #
∂2 1 ∂2U

+ 1 ∂U
H+ = Q Q = −D 2 + −
∂x 4D ∂x 2 ∂x2
" 2 #
∂2 1 ∂2U

+ 1 ∂U
H− = QQ = −D 2 + +
∂x 4D ∂x 2 ∂x2
Soit ψn un état propre de H+
H+ ψn = En ψn
Q+ Qψn = En ψn
Agissant avec Q il vient :
QQ+ Qψn = En Qψn
Si Qψn 6= 0 on en déduit que Qψn est un état propre de H− .

Si Qψn = 0 ⇒ ψn = ψ0 est l’état fondamental de H+
H+ $H_{−}$
$E_{2}$ E2
$E_{1}$ Q E1
$E_{0}$
On en déduit que la partie positive du spectre de H+ coı̈ncide avec celui de H− .

Seul le fondamental de H+ n’est pas apparié.
Remarque : on passe de H+ à H− en changeant le signe de la force ∂U ∂x →

∂U
− ∂x , ce qui suggère que cette symétrie est reliée au renversement du sens du
temps.
7.9 Application
7.9.1 Processus de Wiener
√
ẋ = 2D η (t)
∂2
U (x) = 0 H+ = −D
∂x2
7.9. APPLICATION 83
1
Spectre continu donné par ψk (x) = eikx √
2π
2
H+ ψk (x) = D k ψk (x)
Z +∞
2
P (x t|y 0) = dk ψk (x)ψk∗ (y)e−Dk t
−∞
Z
1 2
= eik(x−y)−Dk t dk
2π
1 (x − y)2
=√ exp −
4πDt 4Dt
Il n’y a pas de solution d’équilibre
7.9.2 Processus d’Ornstein-Uhlenbeck

√
ẋ = −kx + 2D η(t)
∂
φ(x) = −kx = − U (x)
∂x
kx2 kx2
U (x) = =⇒ Peq (x) = exp −
2 2
La relaxation vers l’équilibre est décrite par le hamiltonien
∂2 U 02 U 00 ∂2 k 2 x2 k
H+ = −D 2
+ − = −D 2 + −
∂x 4D 2 ∂x 4D 2
V (x)
gap?En = (n + 12 )ω
x
solution d’équilibre
Pour retrouver la correspondance habituelle avec l’opérateur de Schrödinger

il faut poser
p2 1 1
= D p2 ⇒ D = ⇒m=
2m 2m 2D
mω 2 x2 k 2 x2 1 ω2 2 k 2 x2
= ⇒ x = ⇒ω=k
2 4D 2D 2 4D
On en déduit les valeurs propres
1 k
En = (n + )k − = nk
2 2
et le temps de relaxation
1
τ=
k
7.9.3 Processus bistables

Il s’agit d’étudier la relaxation lente vers l’équilibre en présence d’un potentiel
bistable du type
U (x)
Pour D → 0 soit T → 0 on s’attend à observer une relaxation très lente. Cette

relaxation est gouvernée par une valeur propre λ1 très petite. Pour la calculer
il est commode d’ utiliser la correspondance H+ → H− .
λ2
λ1 λ1
λo H−
Le fondamental de H− a une énergie λ1 qui peut être évaluée par une méthode
variationelle. De plus on est très souvent dans une situation où le partenaire
supersymétrique a un seul minimum.
V+ V−
Considérons l’exemple suivant (Bernstein et al. PRL 52,1933,(1984)) corres-

4
pondant à un cas critique U (x) = x4
Exemple : posons D = 1 U 0 (x) = x3
∂2 1 3
H+ = − + x6 − x2
∂x2 4 2
7.10. REPRÉSENTATION FONCTIONNELLE DE P (X T |Y 0) 85
∂2 1 3
H− = − + x6 + x2
∂x2 4 2
V+ V−
Méthode variationnelle pour H− en prenant les fonctions d’essai

1
ψ−1 = exp − α1 x2
2
1
ψ−2 = x exp − α2 x2
2
donne
λ1 6 1, 38
λ2 6 4, 51
à comparer avec les valeurs ”exactes” obtenues par diagonalisation d’une matrice
100x100
λ1 = 1, 37 λ2 = 4, 45
7.10 Représentation fonctionnelle de P (x t|y 0)

L’expression de la probabilité de transition en terme du hamiltonien H+ nous
suggère de lécrire sous forme d’intégrale fonctionelle
ψo (x) X
P (x t|y 0) = ψn (x)ψn (y)e−En t
ψo (y) n
Le membre de droite peut s’écrire
ψo (x) < x|e−tH+ |y > ψo−1 (y)
∂ 2
1 ∂2U 1 ∂U 2

avec H+ = −D ∂x2 − 2 ∂x2 + 4D ∂x
∂2
= −D ∂x 2 + Ω(x)
En utilisant la formule de Feynman Kac, il vient

Z t " 2 #
ψ0 (x) xt 1 ∂2U
Z
1 ∂U
P (x t|y 0) = dW exp − dτ −
ψ0 (y) y0 0 4D ∂x 2 ∂x2
" 2 #
xt t
1 ∂2U
Z Z
1 1 ∂U
= exp − [U (x) − U (y)] dW exp − dτ −
2D y0 0 4D ∂x 2 ∂x2
Revenons à la mesure de Feynman en posant
Z t
1 1 ∂U
− [U (x) − U (y)] = − ẋ dτ
2D 2D o ∂x
Z 2
1 dx
dW = Dx exp − dτ
4D dτ
Il vient
Z Z t
P (x t|y 0) = D x(τ ) exp − Ldτ
0
2
1 ∂2U

1 2 1 ∂U 1 ∂U
L= x + ẋ + −
4D 2D ∂x 4D ∂x 2 ∂x2
" 2 #
1 ∂2U

1 ∂U ∂U
= ẋ2 + 2ẋ + −
4D ∂x ∂x 2 ∂x2
2
1 ∂2U

1 ∂U
= ẋ + −
4D ∂x 2 ∂x2
" Z t 2 #
t
∂2U
Z Z
1 ∂U 1
P (x t|y 0) = D x(τ ) exp − ẋ + −
4D 0 ∂x 2 o ∂x2
Cette expression suggère de traiter D comme un petit paramètre au même titre

que ~. La méthode du col nous invite à rechercher les solutions qui rendent
stationnaires le premier terme. Les solutions d’action nulle ẋ = − ∂U∂x ne sont
autres que les solutions des équations du mouvement déterministe ! Mais ce
ne sont pas les bonnes solutions. En traitant le deuxième terme comme une
perturbation on peut retrouver aisément certains résultats classiques
voir Onsager-Machlup, Phys Rev 37, 405 (31) et Phys Rev 38, 2265 (31).
7.11 Limite D → O, formule d’Arrhénius

2
On néglige le terme en ∂∂xU2 qui est d’ordre 1.
Il reste
Z Z t 2
1 ∂U
P (x t|y 0) w Dx(τ ) exp − ẋ + dτ
4D 0 ∂x
dans la limite D → O l’intégrale est dominée par les trajectoires qui minimisent
Z t 2 Z t " 2 #
∂U 2 ∂U
ẋ + = ẋ + dτ + 2(U [x] − U (y])
0 ∂x 0 ∂x
7.11. LIMITE D → O, FORMULE D’ARRHÉNIUS 87
Supposons que U (x) = λ(x2 − a2 )2

Ux
Uo = λa2
-a 0 a
Calculons P (0, t|−a, 0) dans la limite t → ∞. Cette quantité donne la probabilité

de remonter la barrière de potentiel sous l’effet d’une fluctuation.
Rt 2
Cherchons les solutions classiques d’action S = 0 [ẋ2 + ( ∂∂xU )2 ]dτ
L’équation du mouvement
2
∂U ∂ 2 U

1 ∂U
2ẍ = 2 2
donne le mouvement dans le potentiel −
∂x ∂x 2 ∂x
∂U ∂ 2 U
2ẍ ẋ = 2ẋ
∂x ∂x2
2
∂ ẋ2

∂ ∂U
=
∂t ∂t ∂x

∂U 2
− ∂x
-a 0
2
∂U ∂U
ẋ2 − = E = O ⇒ ẋ = + > 0 à comparer à l’équation déterministe de signe opposé !
∂x ∂x
Donc
Z ∞ Z ∞
∂U 2 2
S=2 [ẋ + ( ) ]dτ = 2 ẋ2 dτ
0 ∂x 0
Z ∞ Z ∞
∂U ∂U
=2 ẋ dτ = 2 dτ = 2[U (0) − U (−a)]
0 ∂x 0 ∂τ
En regroupant les termes on obtient
"Z 2 #
t
1 2 ∂U
S= [ẋ + ]dτ + 2 (U (x) − U (y))
4D 0 ∂x
(Z 2 )
∞
1 2 ∂U
= [ẋ + ]dτ + 2 [U (0) − U (−a)]
4D 0 ∂x
1 U0
= [4[U (0) − U (−a)]] = 4
4D D
où U0 est la hauteur de la barrière. D’où la formule d’Arrhénius
Uo Uo
P ∼ exp − = exp −
D kT
Chapitre 8
Intégrale de chemin sur un

espace multiplement
connexe
Sommaire
8.1 Electron dans un anneau mésoscopique . . . . . . 89
8.2 Notions d’homotopie . . . . . . . . . . . . . . . . . 90
8.3 Intégrale de chemin sur un espace multiplement
connexe . . . . . . . . . . . . . . . . . . . . . . . . . 91
8.4 Effet Bohm-Aharonov . . . . . . . . . . . . . . . . . 92
8.4.1 Généralités . . . . . . . . . . . . . . . . . . . . . . . 92
8.4.2 Explication simplifiée . . . . . . . . . . . . . . . . . 93
8.4.3 Ligne de flux magnétique . . . . . . . . . . . . . . . 94
8.1 Electron dans un anneau mésoscopique

Considérons un système physique unidimensionnel dont l’espace de configura-
tion C est un cercle.
Exemples : rotateur plan, électron dans un anneau mésoscopique
Il est décrit par une coordonnée généralisée ϕ (t) qui satisfait 0 ≤ ϕ ≤ 2π

avec identification des points 0 et 2π.
Partant du lagrangien du rotateur plan
2
I dϕ
L=
2 dt
nous nous proposons d’évaluer le propagateur de Feynman K (ϕ0 t0 |ϕt). Nous de-
vons pour cela prendre en compte tous les chemins reliant ϕ et ϕ0 . L’ensemble
de ces chemins {ϕn (τ )} est indexé par un entier n qui est le nombre d’enroule-
ments du chemin, nombre de fois où la particule a entouré l’origine.
89
90CHAPITRE 8. INTÉGRALE DE CHEMIN SUR UN ESPACE MULTIPLEMENT CONNEXE
Essayons de relier le problème classique au problème quantique.
Les fonctions d’onde du rotateur s’écrivent
in2 t

1
ψn (ϕ, t) = √ exp inϕ − , n∈ Z
2π 2I
Par conséquent le propagateur peut s’écrire

X
K (ϕ0 t0 | ϕt) = ψn (ϕ0 , t0 ) ψn∗ (ϕ, t)
n
in2 (t0 − t)

1 X 0
= exp in (ϕ − ϕ) −
2π n 2I
En utilisant la formule de Poisson

X X
f (n) = fˆ (2πn)
n n
où fˆ est la transformée de Fourier

Z +∞
fˆ (k) = eikx f (x) dx
−∞
On obtient une expression équivalente du propagateur

2
s i I (ϕ0 − ϕ + 2πn)
I X
2 (t0 − t)
K (ϕ0 t0 | ϕt) = e
2πi (t0 − t) n
L’argument de l’exponentielle
2
I (ϕ0 − ϕ + 2πn)
Sn =
2 (t0 − t)
s’interprète comme l’action classique associée à un chemin ϕn (τ ) de nombre

d’enroulement n. Chaque propagateur partiel peut être interprété comme un
propagateur libre sur l’espace de recouvrement universel de C qui est R. Le
point ϕ0 a une infinité de préimages p−1 (ϕ0 ) = ϕ0 + 2πn.
8.2 Notions d’homotopie

Soit C un espace topologique. On munit l’ensemble des chemins fermés (ou
dont les extrémités sont fixées) d’une loi de composition qui consiste à recol-
ler les chemins. Deux chemins déformables l’ un dans l’autre seront considérés
comme équivalents. On montre que les classes d’équivalence de chemins forment
un groupe appelé groupe fondamental.
Définitions :
8.3. INTÉGRALE DE CHEMIN SUR UN ESPACE MULTIPLEMENT CONNEXE91
1. un chemin d’extrémités x et y est une application continue f : I → C de

l’intervalle I = [0 6 t 6 1] dans C telle que f (0) = x , f (1) = y
2. un chemin est dit fermé si x = y
3. loi de composition interne des chemins
f I → C d’extrémités x et y
g I → C d’extrémités y et z
f ◦ g est défini par l’application
f ◦ g (t) = f (2t) pour 0 6 t 6 12

f ◦ g (t) = g (2t − 1) pour 12 6 t 6 1
4. deux chemins f : I → C et g I → C de mêmes extrémités x et y sont homo-
topes s’il existe une application continue ϕ : J → C où J = {0 6 t, u 6 1}
telle que
ϕ(t, 0) = f (t) f
ϕ(t, 1) = g(t) x
y
ϕ(0, u) = x
ϕ(1, u) = y
g
L’ensemble des chemins fermés issus de x considérés à une homotopie près
forme un groupe pour la loi de composition des chemins.
Ce groupe, appellé groupe fondamental, est noté π1 (C, x). On montre que le
point de base x peut être choisi arbitrairement.
Exemples :
• C est le cercle S1 :
π1 (S1 ) = π1 (SO (2)) = Z
• C est une sphère avec p anses :

le groupe fondamental est un groupe à 2 p générateurs. Exemple : le tore
p=1
8.3 Intégrale de chemin sur un espace multiple-

ment connexe
Soit C un espace de configuration multiplement
connexe et π1 (C) son groupe
~0 0
fondamental. Pour définir le propagateur K q t | ~qt sur (C) on commence
par regrouper les chemins déformables les uns dans les autres, c’est à dire les
chemins appartenant à une même classe d’homotopie g. On définit ainsi un pro-
pagateur partiel
Z q~(t0 )=q~0
S
D~q (τ ) exp i
q
~(t)=~
q ~
On somme ensuite ces différentes contributions avec des phases χ (g) ∈ C où g
étiquette les différentes classes d’équivalence de chemins.
Z q~(t0 )=q~0

~0 0
X S
K q t | ~qt = χ (g) D~q (τ ) exp i
q
~(t)=~ q ~
gπ1 (C)
En utilisant le principe de superposition on montre que les phases doivent sa-

tisfaire
χ (g) χ (h) = χ (g ◦ h)
Les nombres complexes χ (g) de module | χ (g) |= 1 forment une représentation
unitaire de dimension 1 de π1 (C). Ce sont les caractères de ce groupe.
Seule la physique permet de fixer ces phases. Dans le cas d’une particule chargée
couplée à une ligne magnétique traversée par un flux φ nous montrons que
φ
χ (g) = eie ~ n
8.4 Effet Bohm-Aharonov

8.4.1 Généralités
En physique classique, les champs E~ et B
~ agissent sur les particules chargées en
créant une force
F~ = q E~ + ~v ∧ B
~
En l’absence de champ, aucune force ne s’exerce sur la particule. Il ne peut

donc y avoir d’effet physique pour une particule se déplaçant dans une région
~ =B
où E ~ = 0.
Aharonov et Bohm ont imaginé une expérience de pensée montrant que ces
considérations cessaient d’être vraies en physique quantique.
L’expérience est une variante de celle des trous d’Young dans laquelle on intro-
duit un champ magnétique localisé dans une région située entre les fentes.
Dans le dispositif expérimental (H. Batelaan et A. Tonomura, Physics Today,

Septembre 2009), on s’arrange au moyen d’un blindage pour que les particules
ne pénètrent pas dans la région hachurée.
L’expérience montre que les franges d’interférence se déplacent lorsqu’on ap-

plique le champ B,~ résultat paradoxal en physique classique. Les particules ne
pouvant pénétrer dans la région hachurée, aucun effet physique ne devrait ap-
paraı̂tre.
8.4. EFFET BOHM-AHARONOV 93
8.4.2 Explication simplifiée

Considérons deux chemins ~q1 (t) et ~q2 (t) issus du même point, passant respec-
tivement par chacune des fentes et arrivant au même point M
q1
M
q2
iS(q1 ) iS(q2 )
L’amplitude totale est e ~ +e ~
iS(q1 )
h i
i
=e ~ 1 + e ~ [S(q2 )−S(q1 )]
Ces deux chemins interfèrent de façon constructive si S (q2 ) − S (q1 ) = 2nπ~.

Etudions comment ce résultat est affecté par la présence du champ B ~ = rot ~
~ A.
~ ~q˙
L (q, q̇) → L (q, q̇) + eA
par conséquent l’action S1 devient

Z
S1 → S1 + e ~ (~q1 ) v~1 dt
A
Z
= S1 + e ~ (~q1 ) d~q1
A
De même Z
S2 → S2 + e ~ (~q2 ) d~q2
A
~ (q) d~q, intégrale prise le long du chemin fermé

H
Donc S2 − S1 → S2 − S1 + e A
q1−1 ◦ q2 . Il apparait par conséquent une phase supplémentaire ~e Ad~ ~ q = eφ et
R
~
eφ
par conséquent une modification des franges d’interférences ( sauf si = 2nπ).
~
Problème : ne faut-t-il pas inclure d’autres chemins ?
O
M
Le paragraphe suivant analyse de façon détaillée les différentes classes de

trajectoires qui interviennent dans ce problème.
8.4.3 Ligne de flux magnétique
0 x
p
On décrit le problème en terme des coordonnées polaires r = x2 + y 2 et θ.
Le potentiel vecteur associé à la ligne de flux peut s’écrire

A~ = − Ky , Kx = − K sinθ , K cosθ = K ∇θ ~
r2 r2 r r
En dehors de l’origine Hon a B ~ = rot

~ A~ = 0. Pour un lacet quelconque entourant
~
le point O on a Φ = Ad~r = 2πK. On en déduit que la ligne porte le flux
magnétique φ = 2πK. L’ensemble des chemins (C) d’extrémités ~q et q~0 peut
être indexé par un entier n Z tel que
Z
Aµ dxµ = K (θ0 − θ + 2 π n)
c
En présence d’un champ magnétique, la formule de Feynman-Kac peut s’écrire

sous la forme
!
Z Z β
m~˙2
r
< r~0 |e −βH
|~r >= D~r (τ ) exp − ~ ~r˙ dτ
− ieA
0 2
L’origine du facteur i se retrouve facilement en partant de l’expression du La-

grangien d’une particule de charge e et de masse m plongée dans un champ
magnétique.
m~r˙ 2 ~ ~r˙
L= + eA
2
Le hamiltonien correspondant s’écrit
1 2
H= ~
p~ − eA
2m
~ est le potentiel vec-
Etudions les fonctions d’onde et le spectre de H lorsque A
teur précédent.
~ = 0 le hamiltonien s’écrit
Dans la jauge div A
8.4. EFFET BOHM-AHARONOV 95
1 ~ p + e2 A
~2

H= −4 − 2eA~
2m
e2 K 2

1 2eK ∂
= −4 − 2 + 2
2m ir ∂θ r
" 2 #
1 1 ∂ ∂ 1 ∂
= − r − 2 − ieK
2m r ∂r ∂r r ∂θ
Pour simplifier les calculs posons m = e = 1. Nous allons calculer le propa-

gateur < r~0 |e−βH |~r > et montrer explicitement qu’il s’écrit comme une somme
de propagateurs partiels associés à chaque classe d’homotopie.
√
Les solutions de Hψ = Eψ régulières à l’origine sont ψE,m = √12π Jν r 2E eimθ
où E ∈ R+ , m ∈ Z et ν = |m − K|.
Elles satisfont la relation de fermeture

XZ ∞
dEψE,m r~0 ψE,m ∗ (~r) = δ ~r − r~0
m 0
Par conséquent
1
Z ∞ √ √ 0
Jν r 2E Jν r0 2E eim(θ −θ)
X
< r~0 |e−βH |~r >= dEe−βE
2π 0 m
+∞
~r2 + ~r02
0
1 X im(θ0 −θ) rr
= e Iν exp −
2πβ −∞ β 2β
Reprenons l’expression de l’intégrale de chemin

!
β
m~r˙ 2
Z Z
< r~0 |e−βH |~r >= D~r (τ ) exp − ~ ˙
− ieA~r dτ
0 2
Pour un chemin de nombre d’enroulement n nous pouvons écrire

Z β
exp i ~ ~r˙ dτ = exp ieK (θ0 − θ + 2πn)
A
0
Par conséquent le propagateur doit s’écrire sous la forme

X
< r~0 |e−βH |~r >= exp ieK (θ0 − θ + 2πn) µn
n
où µn est un propagateur partiel associé à la classe d’homotopie n. En utilisant

la formule de Poisson on montre que

~2 + r~02 Z +∞
r 0
1 rr
µn = exp − dt exp it (θ0 − θ + 2nπ) I|t|
2πβ 2β −∞ β
On vérifie que le propagateur est, comme il se doit, une fonction monovaluée

(sous la transformation θ0 → (θ0 + 2π)). En l’absence de flux magnétique K = 0,
on vérifie que
r2 + r02 X im(θ0 −θ)

1 rr0
< r~0 |e−βH |~r >= exp − e I|m| ( )
2πβ 2β m
β
On retrouve comme il se doit le propagateur libre.

Plus généralement pour eK entier, par conséquent pour un flux magnétique quantifié
< r~0 |e−βH |~r >= exp ieK (θ0 − θ) < r~0 |e−βH0 |~r >
Le seul effet du champ est donc d’introduire un facteur de phase global alors
que dans le cas non quantifié chaque propagateur partiel est pondéré par une
phase non triviale.
Remarques : la condition eK = n se réécrit en terme du flux magnétique
2πn
φ= = nφ0
e
où φo = 2π e est le quantum de flux élémentaire. Dans les unités physiques
habituelles
2π~
φ0 =
e
La section efficace de diffusion

dσ sin2 (πα)
=
dθ 2π cos2 (θ/2)
s’annule lorsque α = φ/φ0 est entier c’est à dire lorsque la ligne magnétique
porte un flux quantifié.
Chapitre 9
Statistiques quantiques
Sommaire
9.1 Rappel : problème à une particule . . . . . . . . . 97
9.2 Problème à n particules . . . . . . . . . . . . . . . 98
9.3 Calcul de la fonction de partition bosonique . . . 100
9.4 Discussion qualitative des trajectoires . . . . . . . 101
9.5 Bornes sur les fonctions de partition . . . . . . . . 104
9.1 Rappel : problème à une particule
Z(β) = T race e−βH

Z
= dq < q | e−βH | q >
Z Z t
tH 1
or < q 0 | exp − | q >= Dq (τ ) exp − H(q, q̇)dτ
~ q(t)=q 0 ~ o
q(0)=q
Posons t = β~
Z Z β~
1
< q 0 | exp −βH | q >= Dq (τ ) exp − H(q, q̇)dτ
q(β~)=q 0 ~ 0
q(0)=q
Ainsi
Z Z Z β~
1
Z(β) = dq Dq (τ ) exp − H(q, q̇)dτ
q(0)=q ~ 0
q(β~)=q
~
Il faut donc intégrer sur toutes les trajectoires fermées de période β~ = kT
97
98 CHAPITRE 9. STATISTIQUES QUANTIQUES
9.2 Problème à n particules

Pour un système de n bosons identiques, on écrira
Z
1 X
Z(n, β) = dq1 . . . dqn < qσ1 . . . qσn | e−βH | q1 . . . qn >
n!
σ∈Sn
Soit sous une forme plus compacte

Z
1 X
Z(n, β) = dq < σq | e−βH | q >
n!
σ∈Sn
Montrons que l’intégrale de configuration est la même pour les deux permuta-
tions appartenant à la même classe de conjugaison, ce sont des permutations π
et σ telles que π = rσr−1 .
Z Z
dq < σq | e−βH | q >= dq < πq | e−βH | q >
Preuve :
L’élément de matrice densité de e−βH est invariant par permutation par

conséquent
< πq | e−βH | q >=< r−1 πq | e−βH | r−1 q >
or r−1 π = σr−1 donc
< πq | e−βH | q >=< σr−1 q | e−βH | r−1 q >
Prenons comme élément d’intégration la variable
q 0 = r−1 q
Il vient
Z Z
dq < πq | e−βH | q >= dq 0 < σq 0 | e−βH | q 0 >
Par conséquent dans la somme

X 1 Z
dq < σq | e−βH | q >
n!
σ∈Sn
On pourra regrouper toutes les permutations appartenant à la même classe

de conjugaison de Sn . On est ramené à faire une somme sur les classes de
conjugaison P de permutations.
X 1 Z
C(σ) dq < σq | e−βH | q >
n!
σ∈P
Pour poursuivre le calcul nous utilisons les propriétés suivantes

9.2. PROBLÈME À N PARTICULES 99
1. Toute permutation σn peut se décomposer en produit de permutations

cycliques
exemple :

1 2 3 4 5 6 7 1 4 2 3 7 5 6
σ7 = =
4 2 7 1 3 6 5 4 1 2 7 5 3 6
Cette décomposition permet de réécrire σ7 comme le produit de quatre

permutations cycliques avec deux cycles de longueur 1, un cycle de lon-
gueur 2 et un cycle de longueur 3.
De façon générale on pourra décomposer (de façon non unique) σn en un

produit de n1 cycles Pde longueur 1, n2 cycles de longueur 2 et ns cycles
de longueur s avec sns = n
2. Deux permutations quelconques τn et πn = rτn r−1 appartenant à la

même classe admettent la même décomposition en cycles. Par conséquent
les classes de conjugaison peuvent être indexées par une suite d’entiers
(n1 , n2 , . . . ns )
exemple (énumération des classes de σ3 ) :
1 2 3
n1 = 3, n2 = 0, n3 = 0
(a) 1 2 3
1 2 3 1 2 3 1 2 3
n1 = 1, n2 = 1, n3 = 0
(b) 1 2 3 1 2 3 1 2 3
1 2 3 1 2 3
n1 = 0, n2 = 0, n3 = 1
(c) 1 2 3 1 2 3
3. Calcul du nombre de permutationsP appartenant à la classe de conjugaison

(n1 , n2 , n3 , . . . ns ) avec sns = n.
Il nous faut donc calculer le nombre de partitions de n en n1 , n2 , . . . ns

parties de longueurs respectives (1, 2, 3, . . . s).
On obtient
n!
C(n1 , n2 , . . . ns ) = Q Q ns
s ns ! ss
Q
Le premier terme du dénominateur s ns ! prend en compte des échanges
de cycles de même longueur, le second terme prend en compte les permu-
tations cycliques à l’intérieur d’un cycle de longueur s.
4. Munis de ces propriétés, la somme sur les configurations se ramène au
calcul de l’intégrale sur des échanges cycliques de s éléments.
exemple :
~q1 , . . . , ~qs → (~q2 , . . . , ~q1 )
et supposons pour simplifier que le hamiltonien est séparable
< ~q2 , . . . ~q1 | exp −βH | ~q1 , . . . ~qs >
Z q~2 ,β~ Z q~1 ,β~ s 2 Z β~ Xs
m β~ X dqi
Z
= D~q1 . . . . . . D~qs exp − + W (qi (τ ))dτ
q
~1 ,0 q
~s,0 2~ o i=1 dτ o i=1
= K(~q2 , β~ | ~q1 , 0)K(~q3 , β~ | ~q2 , 0) . . . K(~q1 , β~ | ~qs , 0)
Où
X
K(~q1 , β~ | ~qs 0) = ψn (~q1 )ψn (~qs )∗ exp −βEn
n
Après intégration nous obtenons
Z X n X
d3 ~qi < ~q2 , . . . ~q1 | exp −βHn | ~q1 , . . . ~qs >= exp −βsEn ≡ fs
i=1 n
Il est clair que n’importe quelle autre permutation cyclique de s éléments

conduirait au même résultat.
9.3 Calcul de la fonction de partition bosonique

(voir R.P. Feynman, Statistical Mechanics, a set of lectures, p 63)
X 1 Z
Z(β) = C(σ) dq < σq | e−βH | q >
n!
Z
X 1
= Q Q ns dq < σq | e−βH | q >
s n s ! s s
{ns }
Dans le cas d’un hamiltonien séparable on obtient

X 1 Y
Z(β) = Q Q ns (fs )ns
s n s ! s s s
{ns }
XY 1 ns
fs
=
s
ns ! s
{ns }
9.4. DISCUSSION QUALITATIVE DES TRAJECTOIRES 101
P
La somme est prise sur toutes les configurations telles que sns = n.
Pour lever cette contrainte et disposer d’une formule explicite il est commode
de travailler dans l’ensemble grand canonique.
∞
X
Ξ(z, β) = z n Z(β)
n=o
P
où n = sns d’où
∞ n
XY 1 fs s s
Ξ= z
ns ! s
{ns } s=1
∞ X ∞ n
Y 1 fs s s
= z
s=1 n =o s
n ! s
s
∞ ∞
Y fs s X fs s
= exp z = exp z
s=1
s s=1
s
Calculons l’argument de l’exponentielle

∞ ∞ X s
X fs X z exp −βsEn
zs =
s=1
s s=1 n
s
∞
XX (z exp −βEn )s
=
n s=1
s
X
=− log(1 − z exp −βEn )
n
d’où
Y
ΞB (z, β) = (1 − ze−βEn )−1
n
Equation d’état
pV X
= log ΞB = − log(1 − ze−βEn )
kT n
Pour un ensemble de fermions indépendants on obtiendrait

Y
ΞF = (1 + ze−βEn )
n
9.4 Discussion qualitative des trajectoires

• Extension typique des trajectoires
Réécrivons l’intégrale de chemin en prenant comme variable d’intégration

~
t= τ
m
Le terme d’action devient

2
m β~ X d~qi 1 β~ X
Z Z
dτ + W (~qi (τ ))dτ
2~ o i
dτ ~ o i
β~2 β~2
1
Z m X d~qi 2 m
Z m X
= dt + 2 W (~qi (τ ))dτ
2 o i
dt ~ o i
Le terme cinétique décrit n mouvements browniens indépendants sur l’in-

2
tervalle temporel t = β~ m . L’extension typique L des trajectoires est telle
q
2 β~2 β
que < L >∼ m donc L ∼ ~ m .
Comparons l’extension typique L à la distance moyenne d entre particules.
1
Pour une densité fixée on a d ∼ ρ− 3 .
• Estimation de la température critique Tc .
Si L d, on s’attend à ce que les effets quantiques dominent car les
échanges de particules deviennent significatifs
β
0
q
− 31
L > d définit une échelle de température T < Tc telle que ~ βc
m =ρ
2
~2 ρ 3
soit Tc =
km
qui est de l’ordre de la température de dégénérescence d’un gaz de Bose.
Dans la limite de haute température, seuls les chemins de petite longueur

vont contribuer. On s’attend donc à ce qu’il n’y ait pas d’échanges de particules.
Seule la permutation identique va contribuer
< q1 . . . qn | e−βH | q1 . . . qn >
X Z
1
Z β~ X d~qi 2
∼ exp −β W (~qi ) Dq(τ ) exp − dτ
~ 0 i
dτ
X
= exp −β W (qi )K0 (q, β~ | q, 0)n
3n
X m 2
= exp −β W (qi )
2π~β~
Donc
Z n 3n
1 X m 2
Z= d~q1 . . . d~qn exp −β W (qi )
n! i=1
2πβ~2
On retrouve ainsi l’expression de la fonction de partition classique dans l’ap-

proximation de Maxwell-Boltzmann. L’approximation s’appuie sur le fait que
9.4. DISCUSSION QUALITATIVE DES TRAJECTOIRES 103
les trajets se font en un temps très court ∼ β~. Les particules ne peuvent pas
faire d’excursions lointaines.
Vérification
3n Z n
p2
Z
1 1 2mπ 2
h i
n n −β 2m +V (q) −βV (q)
d pd qe = e dq
(2π~)3n (2π~)3n β
3n Z n
2mπ 2
−βV (q) 3
= e d q
4π 2 ~2 β
3n Z n
m 2
−βV (q) 3
= e d q
2πβ~2
Application : étude de la transition λ dans He4

C
2, 18C T
Les effets quantiques interviennent s’il y a des échanges de particules. On
peut visualiser en dessinant les cycles correspondants
A
D G
1 cycle à 3 particules
1 cycle à 4 particules
B 1 cycle à 1 particule
F
D
E H
C
Chaque cycle est associé à un polygone particulier dont le poids statistique
est
32
m m X
exp − (qi − qσi )2
2πβ~2 2~2 β
Feynman suppose les atomes distribués uniformément ainsi
X mn(ρ)d2
Z∼ exp −
2~2 β
P
où la somme porte sur tous les polygones ayant n(ρ) côtés.
Voir Feynman, Phys Rev D 91 (1953) et Statistical Mechanics chap 11 p 343

D.M. Ceperley, Path integrals in the theory of condensed Helium, Rev. Mod.
Phys. 67,279 (1995).
9.5 Bornes sur les fonctions de partition

On peut rendre rigoureuses les considérations précédentes en reprenant l’expres-
sion de Z(β)
Z
Z(β) = d~q < ~q | e−βH | ~q >
où
Z q
~(β~)=~
q Z β~
1 m ˙2
< ~q | e−βH | ~q >= D~q (τ ) exp − ( ~η (τ ) + V (~q(τ )) dτ
q
~(0)=~
q ~ 0 2
Posons ~q(τ ) = ~q + ~η (τ ) avec ~η (0) = ~η (β~) = ~0 (pont brownien)

Z ~ (β~)=~
η 0 Z β~
−βH 1 m ˙2
< ~q | e | ~q >= D~η (τ ) exp − ( ~η (τ ) + V (~q + ~η (τ )) dτ
~ (0)=~
η 0 ~ 0 2
Utilisons l’inégalité de Jensen. Pour toute fonction convexe

X X
f( αi ui ) 6 αi f (ui )
Application
Z s Z s
1 1
exp g(u)du 6 exp g(u)du
s 0 s 0
Par conséquent
Z β~ Z β~ Z β~
1 1 1
exp − V (~q + ~η (τ )) dτ = exp − βV (~q + ~η (τ )) dτ 6 exp −βV (~q + ~η (τ )) dτ
~ 0 β~ 0 β~ 0
Insérons cette inégalité dans l’expression de Z(β) et intervertissons les intégrations.

Intégrant d’abord sur ~q on voit apparaitre
Z Z
d~q exp −βV (~q + ~η (τ )) = d~q exp −βV (~q)
d’où
Z β~ Z Z β~ Z
1 1 m ˙2
Z(β) 6 dτ D~η (τ ) exp − ~η dτ d~q exp −βV (~q)
β~ 0 η
~ (0)=~
η (β~) ~ 0 2
La première intégrale fonctionnelle s’obtient en partant du propagateur d’une

particule libre. On obtient
d2 Z
m
Z(β) 6 d~q exp −βV (~q)
2π~2 β
On reconnait dans le membre de droite la fonction de partition classique

Z
1 d p2
d −β( 2m +V (~
q)
Zc (β) = d pd qe )
(2π~)d
9.5. BORNES SUR LES FONCTIONS DE PARTITION 105
D’où l’inégalité
Z(β) 6 Zc (β)
Voir K. Symanzik J.M.P. 6 (1965) 1155 et C. Thompson J.M.P. 6 (1965) 11812

Equation fonctionnelle satisfaite par la résolvante
Application à la construction des états liés

calculons pour E ∈ C − [0, ∞[
(E − Ho ) [G (E) − Go ] = (E − Ho ) G (E) − 1
= (E − Ho − V + V ) G (E) − 1
= (E − H) G (E) + V G (E) − 1
= V G (E)
D’où l’équation fonctionnelle
G (E) = G0 + G0 V G (E)
Par itération il vient
G (E) = G0 + G0 V G0 + G0 V G0 V G0 + . . . . . .
Application : montrer que à une dimension tout potentiel attractif donne un

état lié.
Il vient pour le calcul de G0
Z +∞
ψk (x0 )∗ ψk (x)
< x | Go (E) | x0 >= dk
−∞ E − Ek
1 ~2 k 2
on a ψk (x) = √ eikx Ek =
2π 2m
Pour E = − | E |< 0 il vient

+∞ 0
eik(x−x )
Z
1
0
< x | Go (E) | x >= − 2 k 2 dk
2π −∞ | E | + ~2m
+∞ 0
eik(x−x )
Z
m
=− dk
π~2 −∞ k2 + 2m|E|
~2
Cauchy donne
s r
0 m ~2 2m | E |
< x | G0 (E) | x >= − 2 exp − + | x − x0 |
~ 2m | E | ~2
Donc pour E < 0 la résolvante admet le développement

2 Z
0 m 0 m
< x | G (E) | x >= − 2 exp −k0 | x−x | + − 2 exp −k0 | x0 −x00 | V (x00 ) exp −k0 | x00 −x | d
~ k0 ~ k0
Nous allons considérer un puit très peu attractif

V(x)
par conséquent on s’attend à ce que l’état lié s’il existe soit proche de E = 0
donc k0 ∼ 0. Ne gardant que les termes les plus singuliers il vient
2 Z +∞
m m
< x0 | G (E) | x > = − + − V (x00 )dx00 + . . . . . .
~2 k 0 ~2 k 0 −∞
m 1
=− 2 R +∞
~ k0 1 + 2 m
V (x00 )dx00
~ k0 −∞
qui présente un pôle pour

Z
m
k0 = − 2 V (x00 )dx00
~
D’où l’energie de l’état lié
+∞ 2
~2 2
Z
m
E=− k =− 2 V (x)dx
2m 0 2~ −∞
• En dimension 2 , le même argument permet de montrer que tout potentiel

attractif donne un état lié d’énergie exponentiellement petite
2π ~2 2π~2
E0 = − 2
exp − pour V0 → 0
e γ mΩ0 mV0 Ω0
dans le cas d’un potentiel circulaire de profondeur − | V0 |.

γ = 0, 577 est la constante d’Euler.
• En dimension 3 , la profondeur du puit doit excéder une certaine valeur
critique pour qu’il puisse appraı̂tre un état lié.

Integral de Chemin

Transféré par

Droits d'auteur :

Formats disponibles

Integral de Chemin

Transféré par

Informations du document

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

Integral de Chemin

Transféré par

Droits d'auteur :

Formats disponibles

Introduction à l’intégrale

Laboratoire de Physique Théorique et Modèles Statistiques

2 Fonctions de corrélation quantiques, oscillateur harmonique 11

4 Applications à la physique des polymères 41

5 Introduction au mouvement Brownien 51

5.3 Marche aléatoire sur réseau . . . . . . . . . . . . . . . . . . . . . 54

6 Intégrale de Wiener et formule de Feynman Kac 67

7 Équation de Langevin et Fokker-Planck 73

8 Intégrale de chemin sur un espace multiplement connexe 89

1.1 Définition du propagateur de Feynman

donne l’amplitude de probabilité de trouver la particule en q~0 à l’instant

• Le propagateur est la solution de l’équation de Schrödinger

• On observe que dans la représentation de Heisenberg cette amplitude peut

Preuve : Dans la représentation de Heisenberg les vecteurs d’état sont

1.2 Expression du propagateur de la particule

Le propagateur correspondant s’écrit :

En intégrant sur p~ on obtient :

1.3 Construction générale du propagateur

dont on sait donner une formulation mathématiquement satisfaisante (A et B

Posons A = ~1 Ĥo B = ~1 V̂ et t = −iT . On obtient

< q~0 |e−i ~ (Ho +V ) |~q >=

Dans la limite n → ∞ , cette somme discrète constitue une discrétisation plau-

Nous obtenons ainsi l’expression fonctionnelle du propagateur :

1.4 Propriétés du propagateur

4. La connaissance du propagateur nous permet d’accéder au spectre du ha-

Posons q~0 = ~q et intégrons sur ~q on obtient

est la trace de la résolvante

Reprenons l’expression du propagateur

Si l’espace de configuration {~q (t)} noté (C) est multiplement connexe,

Pour satisfaire la loi de composition des propagateurs

les facteurs de phase doivent satisfaire la relation χ (g) χ (h) = χ (g ◦ h)

(C) = S1 , R2 / {O}, χ = einθ

1.5 Propagateur euclidien

Z(β) = T race e−βH

Montrons comment exprimer ces quantités sous forme d’intégrale fonction-

Par conséquent la fonction de partition quantique

s’obtient en intégrant sur toutes les trajectoires fermées parcourues dans

En utilisant les propriétés d’orthogonalité des fonctions fn (τ ) = cos nπτ /β

C’est l’identité de Parseval. Convenons de définir la mesure euclidienne

Le propagateur d’une particule libre (H0 = 21 p2 ) peut ainsi s’écrire

En séparant les intégrales sur ξ0 et ξn on obtient

où E désigne l’espérance mathématique par rapport à la mesure produit

E(q(τ )q(τ 0 )) = min(τ, τ 0 )

Cette fonction de corrélation définit le processus de Wiener, processus

2.1 Expression d’éléments de matrice dans la

< q 0 t0 |q̂ (t1 ) >=< q 0 |e−

Insérant des états intermédiaires et faisant apparaı̂tre les propagateurs, on peut

Si nous réexprimons les deux propagateurs en terme d’une intégrale fonc-

Essayons maintenant d’évaluer :

< q 0 t0 |T q̂ (t2 ) q̂ (t1 ) |qt >

Il convient cette fois ci de découper le chemin en trois chemins élémentaires

2.2 Etude de l’amplitude vide-vide

< q 0 t0 |T q̂ (t1 ) q̂ (t2 ) |qt >

Supposons que t1 > t2 . Insérant des états intermédiaires il vient :

< q 0 t0 |T q̂ (t1 ) q̂ (t2 ) |qt >=

en effet seul le fondamental contribue dans cette limite.

< q 0 t0 |qt >→ ψ0∗ (q 0 ) ψ0 (q) e−

D’où l’expression des fonctions de Green

< q 0 t0 |T q̂ (t1 ) q̂ (t2 ) |qt >

Couplons le système à une force extérieure j (τ ) dépendant du temps