PhysIII Chap-7
PhysIII Chap-7
PhysIII Chap-7
& %
- 2-
Rappel : Diffraction par un réseau
' $
Introduction
La diffraction par les Rayons X (RX) présente certaines analogies avec la diffraction
par un réseau, phénomène que nous avons vu en première année (travaux pratiques).
C’est pour cette raison uniquement que nous commençons ce chapitre par ce rappel
d’optique ondulatoire.
Nous allons supposer que les fentes du réseau sont minces et que l’écran est placé
suffisamment loin du réseau pour que nous puissions faire l’hypothèse que les rayons
issus des fentes et aboutissant en un point P de l’écran sont parallèles.
ve
rs
le
po
in tP
sur
quelques 1000 nm
l'éc
ran
d de l'ordre de
θ
différence de chemin
d entre 2 rayons adjacents
θ
-157-
Rappel : Diffraction par un réseau
' $
θ, chaque ondelette sera en retard (ou en avance) sur sa voisine d’un angle de
2π d sin θ / λ , pour une longueur d’onde λ de la lumière. Une interférence construc-
tive se produira sur l’écran si ce déphasage est un multiple entier de 2π, c.à.d.
d sin θ = m λ m = 0, 1, 2 ...
0 5 10 15 20 25 30 35 40 45
θ (degrés)
Dispersion
-158-
Rappel : Diffraction par un réseau
' $
En spectroscopie, nous aimerions avoir une grande dispersion : nous avons intérêt à
travailler avec un réseau ayant une constante d faible et avec des ordres m élevés.
Pour pouvoir définir le moment où nous ne pouvons plus séparer deux raies, nous
nous mettons dans la même situation que pour le cas de la diffraction par une fente
ou un orifice circulaire :
– Nous ne pouvons plus distinguer deux raies si le maximum de l’un se trouve au
mimimum de l’autre ;
– Le “minimum” est atteint pour ∆θm.l. (m.l. = mi-largeur) quand l’onde émise
par la fente du bas interfère destructivement avec celle émise par la fente d’en
haut : ce cas se présente dans le cas de la diffraction par une fente de largeur a si
a sin θ = λ ; dans notre cas, pour la raie centrale :
N d sin ∆θm.l. = λ ou
∆θ m.l.
λ
Pour toute autre raie à θ : ∆θm.l. = .
N d cos θ
& %
-159-
Rappel : Diffraction par un réseau
' $
m ∆λ λ
En substituant ∆θm.l. dans ∆θ = , nous avons = m ∆λ , d’où :
d cos θ N
λ
R = = Nm
∆λ
& %
-160-
La diffraction des rayons X
' $
Cathode Anode
Les longueurs d’ondes des RX sont de l’ordre de l’Angstrom (10−10 m) (à compa-
rer aux longueurs d’ondes de la lumière visible qui sont de quelques centaines de
nanomètres). Il est évident qu’avec de telles longueurs d’onde, la diffraction par un
réseau ordinaire avec d ∼ 1000 nm tel qu’utilisé en optique ne donnerait rien.
C’est en 1912 que Max von Laue se rendit compte qu’un cristal solide pouvait
former un réseau à trois dimensions sur lequel les RX pouvaient “diffracter”. Dans
ce chapitre, nous allons donner quelques relations de base de la diffraction des RX,
puis nous donnerons quelques directions pour leur utilisation en Biologie.
& %
-161-
La diffraction des rayons X
' $
Nous savons que ce sont les électrons des atomes qui forment le réseau cristallin qui
re-diffusent les ondes électromagnétiques incidentes. Cependant, dans le formalisme
que nous allons développer maintenant, nous n’allons pas traiter de la diffusion, mais
montrer que “tout se passe comme si...”
Mais cette condition n’est pas la seule : il faut aussi que les “réflexions” sur les
θ
a)
1 2 α
ce
α θ
d' a b
& %
-162-
La diffraction des rayons X
' $
θ
b)
1 2 α
e
a θ
d α c
Les deux conditions ci-dessus ne peuvent pas être satisfaites simultanément, sauf si
θ = α. La première équation se réduit à zéro, tandis que nous obtenons la loi de
Bragg avec la seconde condition.
2d sin θ = m λ avec m = 1, 2, 3, ...
Résumé : Les deux conditions ci-dessus pour une interférence constructive re-
viennent donc à celles-ci :
– L’angle que fait le rayon incident doit être le même que celui du rayon réfléchi.
– Les réflexions sur différents “plans cristallins” doivent être en phase, c.à.d. doivent
satisfaire à
2d sin θ = m λ avec m = 1, 2, 3, ... Loi de Bragg
Na+
Cl-
a0
a0
Tout se passe comme si, pour les maxima d’intensité, les RX étaient réfléchis
par une série de plans réfléchissants ou plans cristallins qui s’étendent au travers du
cristal et qui contiennent des arrangements d’atomes.
& %
-163-
La diffraction des rayons X
' $
Remarque. Les RX ne sont pas réfléchis : l’utilisation des plans cristallins facilitent seulement
l’analyse de phénomène.
La figure ci-après montre la “réflexion” de trois rayons sur trois plans successifs. En
diffraction des RX, contrairement à ce que nous avons en Optique Géométrique, les
angles d’incidence et de réflexion sont mesurés à partir du plan de réflexion et non
pas à partir de la normale à ce plan.
3 2 1
θ θ
d θ θ
d θ θ
2 1
θ θ θ θ
θ θ
d d sin θ d sin θ
θ θ
Pour expliquer les maxima d’intensité, nous disons que les rayons 1 et 2 sont tels
que leurs différence de chemin, 2d sin θ , doit être égal à un nombre entier fois la
longueur d’onde λ du RX, d étant la distance entre “plans cristallins” :
2d sin θ = m λ avec m = 1, 2, 3, ... Loi de Bragg
Remarques :
– Les rayons 1 et 2 arrivent en phase sur le “plan cristallin”.
– Les RX ne sont pas réfractés : nous n’avons pas introduit d’indice de réfraction.
Pour d’autres orientations du faisceau de RX par rapport au cristal, d’autres “plans
cristallins” “réfléchissent” les RX et nous pouvons y appliquer la loi de Bragg
& %
-164-
La diffraction des rayons X
' $
d θ
θ
d θ
θ
θ
θ
La figure ci-après montre la relation que l’on peut trouver entre la distance entre
“plans cristallins” et caractéristique du réseau cristallin (a0 est la ”maille” du réseau
cristallin) :
a0
0.5a 0
& %
-165-
La diffraction des rayons X
' $
Obtention de l’image
Echantillon
Absorbeur
de Pb
θ
2θ
Film RX
direction incidente du RX
a) b)
Pour chaque “plan cristallin”, les RX diffractés sont sur un cône d’angle au som-
met de 2θ (a)) qui coupent les plans de films RX disposés tout autour selon des
lignes bien définis (b)). Ces séries de lignes obtenues seront ensuite comparées au
prédictions théoriques : nous devons alors identifier les raies en connaissant la struc-
ture cristalline de l’échantillon étudié.
& %
-166-
Détermination de la structure d’un cristal
' $
Le problème de la phase
Jusqu’à présent, nous avons placé un seul atome, donc un seul “plan cristallin”, à
une position de la maille du cristal. Ce qui nous intéresse, ce sont des cristaux qui
ont plusieurs (des dizaines ou des centaines) d’atomes dans une maille. Comment
cela va-t-il affecter le traitement de Bragg que nous venons de voir ?
Considérons pour commencer une maille cubique face centrée (figure de gauche) et
remplaçons chacun des atomes par une molécule di-atomique. Quand la condition
de Bragg est satisfaite, les réflexions sur les atomes “pleins” sont en phase, comme
celles sur les atomes “vides”.
A chaque emplacement
Maille cubique
de la maille, nous avons
face centrée
une molécule diatomique
Cependant, les ondes diffusées sur les plans cristallins “pleins” et “vides” sont
légèrement déphasés ; par conséquent, l’amplitude résultante, et donc l’intensité,
sera diminuée par interférence.
Le problème est d’obtenir une expression générale pour la différence de phase. Nous
allons donner ci-après les points essentiels de la démarche utilisée pour déterminer
la structure de gros cristaux comme ceux que l’on rencontre en Biologie.
Différence de phase Les faces d’un cristal et les plans dans le cristal peuvent
être caractérisés dans un système de trois axes. Sur la figure ci-dessous, les trois
axes ont des longueurs a, b, et c et coupent le plan ABC. On définit
a b c
= h , = k , = l
OA OB OC
Les indices (h, k, l) sont les indices de Miller et permettent de caractériser les
plans du cristal. Un exemple de ces indices est montré sur les figures de droite
pour une maille cubique.
& %
-167-
Détermination de la structure d’un cristal
' $
A b
B
(100) (110)
O C c
z
y
x
(111)
Retour au problème
Pour revenir à la discussion du début de ce paragraphe, considérons une série
de plans cristallins pour lesquels la loi de Bragg est satisfaite. En passant d’un
plan à un autre, nous avons un déphasage de 2π pour une distance de a/h.
Pour le plan des atomes “vides”, situé à une distance x, ce déphasage sera de
2 π x (h/a) (figure a) (gauche)) ci-après.
La différence de phase
Pour un point placé à la distance x de l’origine O (figure a) (droite)), l’onde
d’amplitude E0 émis par un atome placé à l’origine sera perçue avec une ampli-
tude E1. Pour un atome placé à r1, l’onde qu’il émet et reçue en x le sera avec
une amplitude différente. Pour deux atomes placés à des positions différentes r1
et r2 , l’amplitude de l’onde résultante en x sera la somme algébrique des deux
ondes E1 et E2 en x
En étendant le raisonnement à trois dimensions, pour un atome i placé en
(x, y, z) dans la maille, le déphasage sera de
hx ky lz
ϕi = 2π + +
a b c
Clairement, la phase ϕi est liée à la position de l’atome dans la maille.
Facteur de structure En utilisant le principe de superposition pour trouver
l’amplitude résultante des diffusions de toutes les ondes diffusées sur les atomes
d’une maille, nous sommons les amplitudes individuelles en utilisant les vecteurs
de Fresnel, c.à.d. en sommant les nombres complexes fi exp(i ϕi) :
X X hx ky lz
F (hkl) = fi exp(i ϕi) = fi exp 2πi + +
i i
a b c
& %
-168-
Détermination de la structure d’un cristal
' $
a/h
a
Origine x
Changement de phase
de 2π dans la direction a
(a) (b)
Fig. 1 – Différence de phase pour deux plans diffuseurs satisfaisant à la loi de Bragg (a). Addition
de deux ondes électromagnétiques dont les origines sont déplacées l’une par rapport à l’autre (b).
F est appelé facteur de structure du cristal. La sommation se fait sur tous les
atomes. Sa valeur est déterminée par les positions des atomes et de leur facteur
de diffusion fi. Ce dernier dépend du nombre et de la distribution des électrons
dans l’atome (n’oublions pas que ce sont les électrons qui re-diffusent les ondes
électromagnétiques !) et de l’angle de diffusion θ. Nous mesurons l’intensité
totale de la lumière diffractée et non pas l’amplitude résultante de la lumière :
cette intensité est égale au module au carré de l’ampitude F (hkl) :
I = |F (hkl)|2 = F (hkl) F ?(hkl)
& %
-169-
Détermination de la structure d’un cristal
' $
On montre que les coefficients de Fourier Apqr sont égaux aux facteurs de struc-
ture F (hkl) divisés par le volume de la maille (voir par exemple The Crystalline
State L. Bragg, p. 221 et suivantes) :
1 X hx ky lz
ρ(x, y, z) = F (hkl) exp − 2πi + +
V a b c
h,k,l
avec :
X X hx ky lz
F (hkl) = fi exp(i ϕi) = fi exp 2πi + +
i i
a b c
La sommation est faite sur toutes les valeurs de h, k, l pour avoir un terme
pour chaque plan (hkl) et donc pour chaque maxima de la figure de diffraction.
L’équation donnant ρ(x, y, z) précédente permet la détermination de la
structure du cristal, puisque la structure du cristal, c’est ρ(x, y, z). Les maxima
de la figure de diffraction correspondent aux pics de la distribution densité
d’électrons ρ(x, y, z) : ces pics sont aussi proportionnels au nombre atomique.
Par conséquent, si nous connaissons les facteurs de structure F (hkl), nous
pouvons immédiatement calculer la structure du cristal. Malheureusement, les
résultats de la mesure sont les intensités sur le film RX, c’est à dire une quantité
proportionnelle au module au carré |F (hkl)|2 et nous n’avons pas les phases !
Qu’avons-nous perdu ? L’information que nous n’avons pas est
h px qy rz i
exp − 2πi + +
a b c
c’est à dire la position de l’atome dans la maille, nous laissant seulement avec le
nombre de points diffuseurs ramenés à l’origine de la maille ! En effet, le module
au carré de l’exponentielle complexe est égal à 1 !
& %
-170-
Détermination de la structure d’un cristal
' $
La série de Patterson
N X
N
X h k l
|F (hkl)|2 = fi fj exp 2πi (xi − xj ) + (yi − yj ) + (zi − zj )
i j
a b c
avec u = |xi − xj | , v = |yi − yj | , w = |zi − zj |
+∞ +∞ +∞
1 X X X hx ky lz
ρ(x, y, z) = F (hkl) exp − 2πi + +
V a b c
h = −∞ k = −∞ l = −∞
X X hx ky lz
avec F (hkl) = fi exp(i ϕi) = fi exp 2πi + +
i i
a b c
& %
-171-
Détermination de la structure d’un cristal
' $
+∞ +∞ +∞
1 X X X
2 hu kv lw
P(u, v, w) = |F (hkl)| exp − 2πi + +
V a b c
h = −∞ k = −∞ l = −∞
N X
N
X h k l
|F (hkl)|2 = fi fj exp 2πi (xi − xj ) + (yi − yj ) + (zi − zj )
i j
a b c
2’- 2
2’- 3 3’- 2
1-2 2’- 1
2’- 3’ 3-2
2
3’- 3
1-3 3’- 1
3
1
1 - 3’ 3-1
3’
3 - 3’
2-3 3’- 2’
1 - 2’ 2-1
2’
a) 2 - 3’ 3 - 2’
2 - 2’ b)
maille maille
zi zi zj
fi fi f j
yi
yi yj
xi xi xj
pic de “Fourier” pic de “Patterson”
Par la transformée de Fourier, nous sommes sensibles à position des atomes alors
que dans la série de Patterson, nous sommes sensibles aux distances interatomiques.
& %
-172-
Détermination de la structure d’un cristal
' $
Les pics fi fj , (xi − xj ) et fj fj , (xij − xi) ont même intensité et des coordonnées
opposées
Itération à partir d’une structure proche Une des méthodes utilisées est
d’essayer une structure proche et de calculer les intensités. Faire l’hypothèse
d’une structure, c’est avoir la possibilité de calculer les déphasages
hx ky lz
ϕi = 2π + +
a b c
On calcule alors les séries de Fourier en prenant les |F (hkl)| déduits de la
mesure de diffraction de RX et les déphasages calculés. Si on est sur une bonne
piste, la distribution ρ(xyz) donnera de nouvelles positions des atomes que l’on
pourra utiliser pour obtenir de nouveaux déphasages.
C’est ce qu’on appelle le remplacement moléculaire. Cette méthode peut
être utilisée pour des acides nucléiques et pour des protéines. Pour une structure
nouvelle et totalement inconnue, la méthode ne marchera pas.
Dans l’itération, plus le nombre de termes de la série de Fourier est grand, plus
précise sera l’image de la structure.
[La situation est assez semblable au cas d’un microscope ou d’un télescope :
plus on recueille de lumière dans l’instrument d’optique grâce à une ouverture
plus grande, meilleure sera la résolution. Ici, plus on inclut de termes de la série
de Fourier, meilleure sera la précision.]
La figure ci-dessous montre une image de la densité électronique du glycylgly-
cine (H2NCH2CONHCH2CO2H) projetée sur une maille ; le nombre de termes
de la série de Fourier dont on a tenu compte est aussi inscrit. Avec 40 termes,
la résolution est de 0,25 nm, avec 100 termes, elle est de 0,15 nm et avec 160
termes, elle atteint 0,11 nm.
Les méthodes directes Dans les méthodes directes, on peut essayer toutes les
combinaisons de phase et rechercher celle qui “fitte” le mieux les données, ou
& %
-173-
Détermination de la structure d’un cristal
' $
40 termes
100 termes
NH 2 O
160 termes
CH2
C H
C O
CH 2 NH
O
Fig. 3 – Densité des électrons du glycylglycine (H2 NCH2 CONHCH2 CO2 H) projetée sur une maille.
& %
-174-