Diffusion et diffraction Sylvain Ravy SynchrotronSoleil Ligne CRISTAL

Diffusion et diffraction Sylvain Ravy Synchrotron-Soleil, Ligne CRISTAL ravy@synchrotron-soleil. fr A. Interaction Quanton/Matière 1. Les quantons 2. Absorption et diffusion 3. Diffusion Thomson B. Diffusion par un atome 1. Interférences 2. Facteur de diffusion 3. Diffusion résonante C. Diffusion par un corps quelconque 1. Fonction de corrélation de paire 2. Amplitude diffusion 3. Cristal périodique 4. Désordres

Caractéristiques des quantons Description Neutrons Électrons Champ électromagnétique Particule E=E 0 exp(i(k. r-wt)) y ~ exp(i k. r) E=hn=hc/l l(Å)=12398/E(ke. V) l=1 Å, E=12, 4 ke. V n=1018 Hz E=p 2/2 mn l(Å)=0, 286/E 0. 5(e. V) l=1 Å, E=81, 8 me. V l=2 Å, E=20, 45 me. V E=p 2/2 me l(Å)=12, 265/E 0. 5(e. V) l=1 Å, E=150 e. V l=0. 04 Å, E=100 ke. V p p=hk=hn/c p=hk (=mv) 300 K 3. 10 -6 << 1 Énergie E Impulsion k. BT/E Photons X ~ 10 -5 Interaction Charge sth ~ Z 2 barn Moments magnétiques sd ~ 10 -6 barn Noyaux (forte) sd ~ 5 barn Moments magnétiques sd ~ 3 barn Potentiel electrostatique sd ~ 108 barn Absorption 4700 barn (Z=28, 1, 5 Å) Typique : 0, 1 -1 barn -

Rayons X Longueur d’ondes utilisées l = 0. 1 Å to 3 Å et énergie des photons E = ħw = 120 ke. V to 4 ke. V Typiquement : 4 ke. V (3. 1 Å) to 30 ke. V (0. 41 Å) 1 Å=10 -10 m, 1 e. V = 1. 6 10 -19 J, E[ke. V]=12. 398/ l[Å]

A- Interaction Quanton-Matière 1 Absorption et diffusion

Interaction Quanton-Matière « Element of modern x-ray physics » J. Als-Nielsen et D. Mc. Morrow Deux processus d’interaction Absorption et diffusion kd ki ki 2 q Le quanton disparaît (conservation E et p…) Le quanton change de direction de propagation (interférence : diffraction)

Diffusion élastique, diffusion inélastique Diffusion élastique : L’état interne des particules ne change pas durant la collision Diffusions Thomson, Rayleigh : élastiques Diffusions Compton, Raman : inélastiques Mécanique quantique II, p. 894 C. Cohen-Tannoudji, B. Diu, Frank Laloë Cohérent Incohérent kf kf ki ki 2 q La cible ne change pas d’état. En général ki = kf 2 q La cible change d’état. ki ≠ kf

Les électrons dans un atome Vide (continuum) Energie des photons X Ex: EK(Ag)=25 ke. V EK(Kr)=14 ke. V EK(Fe)=7 ke. V L 2 p 3/2 2 p 1/2 2 s Niveaux de cœur (discret) K 1 s

Effet photoélectrique Fluorescence Absorption Diffusion Compton Diffusion résonante Electron Auger -EF Ka Kb Diffusion Thomson -EII -EI Diffusion

Sections efficaces Section efficace de diffusion : (différentielle) Flux incident : en Quanton/s/cm 2 DW kd ki 2 q Unité : le barn 10 -24 cm 2 Section efficace d'absorption :

Section efficace de diffusion • Processus de diffusion kd ki Section efficace différentielle de diffusion • Fonction d’onde du quanton diffusé b(q ) : longueur de diffusion peut-être complexe : déphasage Convention b > 0 • Section efficace différentielle 2 q

A-2 Diffusion Thomson

Diffusion Thomson Photons incidents et diffusés Onde diffusée : sphérique Onde incidente : plane Dipôle oscillant

Diffusion Thomson (cohérente) Longueur de diffusion de Thomson • Approche classique, n’explique pas l’effet Compton • Onde incidente plane a e’ • L’électron est soumis à une force : a. e’ 2 q a e • Cette charge oscillante crée le champ : 2 q p -p e e’ s-s 2 q e’ 2 q e • En introduisant la longueur de diffusion bth :

Section efficace de diffusion • Rayon classique de l’électron • Si le rayonnement est polarisé : • Si le rayonnement n’est pas polarisé • Pour un électron « libre » : bth noyau négligeable car mnoyau>> m

B-Diffusion par un atome 1 Interférences

Interférence Plane wave Fresnel Fraunhofer

Définitions… k : vecteur d’onde (k = 2 p/l) w: pulsation, fréquence angulaire (w=2 pn=2 p/T) k. r – wt +j : phase

Différence de phase en champ lointain 2 q kf 1 ki 2 q r Différence de phase entre les trajets 2 et 1 : D f = ki r - kf r Df = - q r Ne dépend que du : Vecteur de diffusion

Le vecteur de diffusion q ki kf q q r Les ondes diffusées par Tous les points d’un plan deux points distants de r tels que orthogonal à q q. r = m 2π diffusent en phase interfèrent constructivement

Expression de l’onde résultante En pratique, on utilise la notation exponentielle a est l’amplitude complexe

L’amplitude complexe Son amplitude au carré donne l’intensité Sa phase est celle de l’onde diffusée /t à l’origine et elle est perdue ! Sous cette forme, les calculs de diffraction reviennent à calculer des séries de Fourier (discret) ou des transformées de Fourier (continu)

B-2 Facteur de diffusion atomique

Facteur de diffusion atomique • Diffusion par une distribution de charge dr=re(r)d 3 v ki r Onde diffusée est la somme des ondes diffusées par les électrons du volume d 3 r 2 kf 1 • Longueur de diffusion atomique re(r) est la densité électronique de l’atome q vecteur de diffusion

Facteur de diffusion atomique f (q) : transformée de Fourier de r(r) f (q 0) = Z f (q ) = 0 f(q) sont tabulés : Tables internationales de cristallographie q = 4 p sinq/l

…et la diffusion Compton hq hki hkd Diffusion inélastique (incohérent) 2 q -hq Longueur d’onde Compton Å Compton plus intense aux grands angles

Neutrons: Longueur de diffusion Pseudo-potentiel de Fermi Courte portée (fm). b ~ 5 fm s=4 pb 2 ~ 3 barn s. X= Z 2 barn • b ne dépend pas de q • Eléments légers • b dépend de l’isotope • b peut être négatif

Limitations de la théorie (RX) f(q) (Thomson) dépend : i) de la densité électronique ii) du vecteur de diffusion ne dépend pas (directement) de l'énergie. . . Mais !. . .

DAFS : Diffraction Anomalous Fine Structure Intensité des réflexions (111) et (222) du Cu (25 mmx 200 nm) au voisinage du seuil K du Cu. H. Stragier et al. , PRL 69, 3064 (1992)

B-3 Diffusion résonante

Facteur de diffusion résonant f (q, w)= f 0 (q)+ f ’(q, w)+ i f ’’(q, w) f ’(q, w) associé à la dispersion (indice n) f ’’(q, w) associé à l’absorption (théorème optique) Les deux sont reliés par les relations de Kramers-Kronig (anomal)

• Modèle oscillateur amorti • Force de rappel • On pose mw 02 r, force de friction -Kv r=r. Oe-iwt • Le champ diffusé s’écrit • La longueur de diffusion est modifiée Diffusion résonante

Indice de réfraction a n ki kr kt a’ bk=m/2 On montre que pour les rayons X Réfraction : d ~10 -5 Absorption :

Relation de Kramers-Kronig (1926) Le principe de causalité donne une relation entre f' et f" P : valeur principale de Cauchy

Relation de Kramers-Kronig-3 Calcul théorique 2 électrons K w-3 KK à partir de l'absorption f ’’Au f’Au Z=79 ln(1 -xk 2) "Correction" peut être très importante

C-Diffusion par un corps quelconque 1 -Fonction de corrélation de paire

Fonctions de corrélation de paire va=V/N Volume atomique moyen Fonction de corrélation de paire dépendante du temps d 3 r r, t Moyenne spatiale, statistique, temporelle G(r, t) : TF dans le temps et dans l’espace par diffusion de neutron O t=0 Fonction de distribution de paire g(r) G(r, t=0) Diffusion des rayons X : TF de g(r) Fonction de corrélation densité-densité :

La fonction de distribution de paire Pics premier voisin deuxième voisin etc. Largeur du pic : fluctuation de distance Intégrale du pic : nombre de voisins

Les trois types d’ordre • Comportement à grande distance de g(r) définit trois types d’ordre : • Ordre à courte distance • g(r) ~ exp(-|r|/x) • x : longueur de corrélation • Ex : verre, liquide • Ordre max. à 1 D • Quasiordre à grande distance • g(r) ~ |r|-h • Pas d’échelle de longueur • Ex : Smectiques, cristaux 2 D • Ordre max. à 2 D • Ordre à grande distance • g(r) n’a pas de limite à l’infini ! • Ex : Cristaux • Pics de Bragg exp(-|r|/x)

Un cristal donne des pics de diffraction… et réciproquement ! Ordre à grande distance : diffraction Rayons X Électrons Cristal de C 60 Quasi-cristal Existence de taches de Bragg Largeur limitée par la résolution Neutrons Sinon : diffusion diffuse Diffusion répartie uniformément Eau Cristal liquide smectique

2 -Amplitude de diffusion

Diffusion par un corps de structure quelconque • Approximation cinématique • Intensité du faisceau diffracté, négligée devant celle du faisceau incident • Pas de diffusion multiple, pas de perte d’intensité (ne conserve pas l’énergie, § Born) • Approximation valable pour de petits cristaux (mosaïque) • Pas d’approximation Théorie dynamique • Approximation statique • Fréquences des rayons X : 1018 Hz • Fréquences des vibrations atomiques : 1012 Hz (THz) z t • L’énergie des rayons X est grande devant l’énergie des excitations élémentaires t

Calcul de l’intensité diffusée Dans l’approximation cinématique : • L’intégrale est étendue à tout le cristal • r = r(t). . . est l’amplitude complexe de diffusion C’est la transformée de Fourier de la densité électronique totale Intensité diffusée

Fonction de diffusion S(q) La fonction de distribution de paire est mesurable par diffusion des rayons X ou des neutrons Dans le cas d’un cristal g(r) est la fonction de Patterson

Exemple Argon liquide 85 K Neutrons Argon solide c. f. c 6. 5 Å 3. 75 Å 5. 31 Å

3 -Les cristaux périodiques

Les cristaux ne sont pas tous périodiques !

Cristaux apériodiques Incommensurables Propriété locale (ex : polarisation) possède une périodicité incommensurable avec celle du réseau. Méthodes de cristallographie Doivent être changées l et a sont dans un rapport irrationnel Quasi-cristaux (Al. Pd. Mn) 10 1 2 9 72° 36° • Deux types de « tuiles » • Règles d’accord 8 7 3 6 5 4

Le cristal périodique Un cristal : motif associé à un réseau Na = Atome Na. Cl Groupe d’atomes C 60 Molécule Nucléosome Macromolécule Motif Cristal

3 -Le cristal périodique Noeuds périodiquement espacés Paramètres de maille a, b, c, a, b, g c b a b g a

Relation de BRAGG Plans réticulaires : miroirs q q d Si le différence de chemin optique est un multiple entier de la longueur d’onde il y a interférence constructive

La relation de Bragg donne une condition géométrique pour qu’il y ait diffraction mais aucune information quantitative sur l’intensité du faisceau diffracté

Diffusion par un cristal périodique Introduction : Cristal de N*N*N mailles Contenant un atome de facteur de diffusion f • L’amplitude de diffusion est : • Calcul d’une somme géométrique N=8 Fonction de diffusion max. si 1 2 3 4 5 qx

Conditions de Laue - 1 Cristal quelconque • Densité électronique totale rtot(r) • Approximation cinématique • Périodicité parfaite • Densité électronique d’une maille r(r) =

TF de rtot(r) Conditions de Laue - 2 ×

TF de s(r) Transformée de Fourier de la fonction « volume » L 6 0. 88 p/L 2 p/L Exemple d’un cube

Conditions de Laue - 3 • Chaque nœuds du RR remplacé par une fonction S(q) • Taille du cristal >> paramètre de maille : Intensité maximum q appartient au RR

Facteur de structure TF de la densité électronique de la maille On néglige les électrons de liaison : approximation sphérique • h, k, l, indices de Miller, • uj, vj, wj, coordonnées réduites de l’atome (rj = uj a + vj b + wj c) Ex : 2 atomes identiques en +ua et -ua

RX-Neutrons Pour les neutrons même expression, on remplace f par b … et on ne considère que la diffusion élastique

Intensité Diffractée Atome Facteur de diffusion Motif Facteur de structure Réseau direct Cristal Réseau réciproque Facteur de forme • Position des taches de Bragg : réseau • Intensité : motif • Profil : Forme des grains/cristaux

Pour mesurer S(q), il faut que les interférences puissent se former sur toute la taille du cristal Mesure de S(q) Particules d’Au sur substrat Si. O 2 1 mm Petit cristal (~ 1 mm) Faisceau X cohérent (synchrotron 3 e génération) Images SEM Intensité autour de la réflexion (1, 1, -1) mesurée en faisceau cohérent à l’Advanced Photon Source de l’Argonne National Laboratory. D ’après I. Robinson et al. , Phys. Rev. Lett. 87, 195505 (2001)

S. Labat, N. Vaxelaire, IM 2 NP Marseille

Construction d’Ewald Interprétation géometrique de la diffraction • • Diffusion élastique : ki=kf=2 p/l Vecteur de diffusion q appartient au RR Sphère d’Ewald kf q 2 p/l ki Cristal O Origine du RR Condition de diffraction : noeud du RR sur la SE

Sphère d’Ewald http: //marie. epfl. ch/x-ray/

Laue Bragg 2 p/l=k q q=Qhkl O dhkl Si Qmh, mk, ml appartient à la sphère d’Ewald :

Résoudre une structure Mesurer les positions des taches de Bragg Determiner les paramètres de réseau Mesurer l’intensité des taches de Bragg Résoudre le problème des phases pour déterminer la densité électonique

Problème des phases On ne mesure que l’intensité |Fhkl|2 d’une réflexion de Bragg Les phases ne peuvent pas être obtenues expérimentalement mais par calcul. Résolution Les intensités mesurées sont telles que : ||Qhkl|| < Qmax. < 4 p/l kd q 4 p/l ki rtot(r) Sphère de résolution 1. 15 p/Qmax : l/4 ) est convoluée par une fonction de largeur 1. 15 p/Qmax est la résolution ( mini =

Intensité intégrée a’ : vitesse de rotation du cristal • Facteur de Lorentz • Facteur de polarisation Sphère d’Ewald da S(q) 2 p/l d. W q d 3 q s q=Qhkl q os c a qdd d 3 q q q Rayons x

Théorie dynamique-1 Diffraction sur des cristaux parfait • Théorie dynamique (M. Von Laue, P. Ewald, G. Darwin) • Dépend de la géométrie de diffraction • Même conditions de diffraction (Laue, Bragg) à la réfraction près… q q Pouvoir réflecteur (géométrie de Bragg) Th. Cinématique Th. dynamique

Théorie dynamique-2 Pdyn. < Pcin. Réflectivité Extinction secondaire : Grain B moins illuminé que A q q A Cristal mosaïque Idéalement imparfait (Petits cristaux, Poudres) B q « Rocking curves » Extinction primaire : Interférences négatives entre faisceaux diffusés n fois Réflectivité Courbe de Darwin 100 % L L : longueur d’extinction q

3 -Désordre

Les deux types de désordres • Désordre de deuxième espèce (de type liquide) : La fonction de distribution de paire g(r) tend vers 0 Pas de pics de diffraction • Désordre de première espèce (cristal désordonné) : Ordre à grande distance Pics de diffraction (modifiés) ET Diffusion diffuse

• Agitation thermique • À un instant donné, pas de périodicité parfaite • Périodicité rétablie en moyenne Cristal désordonné : 1 -La température • Structure moyenne périodique • Moyenne statistique Moyenne temporelle (Hypothèse ergodique) • Désordre d’orientation • Ex : C 60, cristaux plastiques T=300 K c. f. c. Kroto et al. 1985

réel 0: 2 -Les défauts • Cristal Dimension • Lacunes, intersticiels • Défauts topologiques • Induisent des déformations qui concernent l’environnement atomique local, comme le nombre de voisins acune ntes 00 K) lorés Intersticiel • Plasticité slocation coin (Impureté) • Dopage des semi-cond. • Couleur des joyaux • Plasticité Désinclinaison Faute d’empilement Joint de grain

Expérience Simulation Diffusion diffuse thermique Thermal Diffuse Scattering Si 300 K RX // <111> Fausses couleurs, Échelle log. RX // <100> M. Holt, Phys. Rev. Lett 83, 3317 (1999)

Cristal de C 60 Kroto et al. 1985 T=300 K c. f. c. R. Moret, P. Launois, S. Ravy

Les cristaux désordonnées Réflexions de Bragg Diffusion diffuse ? Expression de l’amplitude complexe Ordre à grande distance mais contenu des mailles dépend de Ruvw

Calcul de l’intensité instantanée N(m) nombre de terme de la somme : N(m)= N : moyenne sur l’espace et le temps Expression générale de l’intensité diffusée

Diffusion diffuse Fn écart à la valeur moyenne de Fn ID(q) : Diffraction Facteur de structure moyen IDD(q) : Diffusion diffuse

Conservation de l’intensité diffusée totale ID diminue IDD augmente Diffraction + Diffusion diffuse

Deux types de désordre Désordre de déplacement Désordre de substitution

Facteur Debye-Waller Un atome à l’origine Cristal harmonique Intensité diminuée du facteur e-2 W Im N q grand, T grand Re Re Facteur Debye-Waller

Facteur Debye-Waller-2 Maille contenant n atomes en rj IDe -2 W Vibrations isotrope Diffraction permet de mesurer :

Exemple : Critère de Lindemann Solide fond quand : Aluminium c. f. c. a=4. 04 Å Fusion

À suivre…

Théorème optique Mécanique quantique II, p. 940 C. Cohen-Tannoudji, B. Diu, Frank Laloë Ombre : Interférence entre onde incidente et onde diffusée