Facteur de diffusion atomique

De testwiki
Aller à la navigation Aller à la recherche

En physique, le facteur de diffusion atomique est une mesure de la puissance de diffusion d'une onde par un atome. Il dépend de la nature des interactions entre l'onde et l'objet diffusant, donc de la radiation incidente : en général, il s'agit de la diffusion élastique de rayons X, de neutrons ou d'électrons. En cristallographie, les facteurs de diffusion atomiques sont utilisés pour calculer le facteur de structure d'une réflexion donnée lors de la diffraction par un cristal.

Définition générale

Schéma montrant les relations entre vecteur d'onde incident, vecteur d'onde diffusé, vecteur de diffusion et angle de diffusion
Vecteur de diffusion : différence entre les vecteurs d'onde des faisceaux diffusé et incident.

Le facteur de diffusion atomique f(K) apparaît lorsque l'on considère l'effet d'un objet (atome ou ion) sur une onde incidente ψi(r), de vecteur d'onde ki et de longueur d'onde λ comparable à la taille de l'objet, dans l'hypothèse de la diffusion élastique. L'onde ψf(r) diffusée au point r par l'objet situé au point origine s'écrit :

ψf(r)=ψi(r)f(K)eikfr

kf est le vecteur d'onde de l'onde diffusée, de même longueur d'onde λ mais de direction différente. Les normes des deux vecteurs d'onde sont donc égales :

||ki||=||kf||=k=2πλ.

L'angle entre les vecteurs d'onde incident et diffusé est l'angle de diffusion 2θ et K est le vecteur de diffusion[note 1], différence des vecteurs d'onde des faisceaux incident et diffusé ki et kf :

K=kfki,||K||=4πλsinθ.

Le facteur de diffusion atomique est la transformée de Fourier de la distribution de densité spatiale ρ(r) de l'objet diffusant, de l'espace réel vers l'espace réciproque. Il est défini par :

f(K)=0ρ(r)eiKrdr.

Du fait des propriétés de la transformée de Fourier, la largeur de la fonction f(K) dans l'espace réciproque est une fonction inverse de celle de la distribution de densité ρ(r) dans l'espace réel : des objets larges conduisent à une décroissance rapide du facteur de diffusion.

Cas des rayons X

Figure montrant les facteurs de diffusion atomique pour les rayons X de l'oxygène, du chlore, de l'anion ClModèle:Exp et du cation KModèle:Exp
Facteurs de diffusion atomique (cas des rayons X) de l'oxygène (bleu), du chlore (vert), de l'anion ClModèle:Exp (magenta) et du cation KModèle:Exp (rouge).

Les rayons X sont diffusés par le nuage électronique de l'atome et l'amplitude de la diffusion des rayons X augmente avec le numéro atomique Z. En conséquence, les rayons X ne permettent pas bien la détection des éléments légers tels que l'hydrogène et il y a très peu de contraste entre des éléments voisins dans le tableau périodique.

Pour la diffusion des rayons X, ρ(r) représente la densité électronique autour du noyau de l'atome. Le facteur de diffusion atomique est alors la transformée de Fourier de cette quantité. Il s'exprime en électrons (eModèle:Exp). L'approximation d'une distribution sphérique des électrons est habituellement suffisante en diffraction des rayons X[1].

Les premiers calculs du facteur de diffusion atomique furent effectués sous l'hypothèse de la diffusion Thomson et pour un atome de symétrie sphérique contenant Z électrons indépendants[2]Modèle:,[3]Modèle:,[4]Modèle:,[5]Modèle:,[6]Modèle:,[7]. Ce facteur de diffusion atomique approché f0 s'écrit alors comme la somme des facteurs de diffusion de chaque électron :

f0(K)=4πi=1Z0ρi(r)sinKrKrr2dr=4πi=1Z0ρi(r)J0(Kr)r2dr

ρi(r) est la densité électronique de l'électron i et J0 est la fonction de Bessel sphérique d'ordre 0. En négligeant les interactions entre les électrons, le facteur de diffusion atomique est donc une grandeur réelle. En particulier, dans le cas d'un vecteur de diffusion nul, le facteur de diffusion atomique est égal au numéro atomique de l'atome : f0(0)=Z. Cette formule est bien adaptée pour les éléments à faible numéro atomique et pour la diffusion de rayons X à courtes longueurs d'onde.

Interpolations

Les valeurs des facteurs de diffusion atomiques pour chaque atome et ion sont tabulées en fonction de sin(θ)/λ dans les tables internationales de cristallographie[8]. En pratique, il est plus simple d'utiliser des fonctions d'interpolation numérique pour décrire les facteurs de diffusion atomiques. La fonction d'interpolation la plus souvent utilisée est valable pour sin(θ)/λ compris entre 0 et Modèle:Unité, cas qui correspond à la résolution de la plupart des expériences de diffraction :

f0(sinθλ)i=14aiebi(sinθλ)2+c.

Les valeurs de a1, b1, a2, b2, a3, b3, a4, b4 et c sont également données dans les tables internationales de cristallographie.

Pour sin(θ)/λ compris entre 2 et Modèle:Unité, le facteur de diffusion atomique est mieux interpolé par la formule logarithmique suivante[9] :

ln[f0(sinθλ)]i=13aisi.

Diffusion anomale

Modèle:Loupe En général, le facteur de diffusion atomique pour les rayons X est une grandeur complexe. La partie imaginaire du facteur de diffusion atomique ne devient importante que lorsque la longueur d'onde des rayons X incidents est proche d'une arête d'absorption de l'atome. On parle dans ce cas de diffusion anomale. Pour en tenir compte, deux termes correctifs sont ajoutés au facteur de diffusion atomique :

f=f0+f+if

f et f décrivent les variations en amplitude et en phase du facteur de diffusion par rapport à f0, dues à l'absorption des rayons X.

Atomes non isolés

Schéma d'une liaison σ entre deux atomes
Liaison σ entre deux atomes : localisation de la densité électronique.

Dans un cristal, les atomes forment des liaisons chimiques entre eux. Ces liaisons déforment le nuage électronique des atomes et la symétrie de l'objet diffusant n'est plus strictement sphérique. L'effet de la déformation du nuage électronique sur le facteur de diffusion atomique est généralement faible, mais il doit être pris en compte lors de l'étude de la densité électronique dans un cristal.

La densité électronique d'un atome non isolé est calculée à partir de celle de l'atome isolé en ajoutant des termes perturbatifs. Le facteur de diffusion de l'atome non isolé est alors la transformée de Fourier de cette densité électronique perturbée.

On distingue trois types de contributions à la densité électronique, provenant d'électrons différents :

  • les électrons de cœur, ne participant pas aux liaisons chimiques, ont une densité sphérique ρc identique à celle dans l'atome isolé ;
  • les liaisons chimiques ont deux effets sur la densité ρv des électrons de valence :
    • le changement de potentiel change la force de l'interaction des électrons de valence avec le noyau atomique, modifiant le rayon de leur distribution spatiale ;
    • la forme de la distribution électronique n'est plus sphérique.

La densité électronique totale d'un atome est la somme de ces trois contributions :

ρ(r)=ρc(r)+ρvs(r)+ρva(r)

les exposants s et a indiquant les contributions sphériques et asphériques des électrons de valence, respectivement.

La contribution sphérique des électrons de valence s'écrit simplement comme la contraction ou l'expansion d'un facteur κ de la densité électronique de valence ρv,0 dans l'atome isolé :

ρvs(r)=κ3pvρv,0(κr)

pv est la population d'électrons de valence contribuant à la partie sphérique de la densité électronique.

Pour la contribution asphérique des électrons de valence, la densité électronique, exprimée dans un système de coordonnées polaires, peut généralement être décomposée en une partie radiale et une partie angulaire :

ρva(r)=ρR(r)ρΘ(θ,φ)=κ'3ρR,0(κr)ρΘ(θ,φ).

La fonction angulaire ρΘ(θ,φ) est développée sur la base des harmoniques sphériques Y,m(θ,φ), définies avec les polynômes de Legendre associés P,m(cosθ) et avec m par :

Y,m(θ,φ)=(1)m(2+1)4π(|m|)!(+|m|)!P,m(cosθ)eimφ.

Comme la densité électronique d'un atome est une grandeur réelle, l'utilisation des harmoniques sphériques réelles normalisées (notées y,m), combinaisons linéaires des harmoniques sphériques, est préférée[10]. Les électrons dont la distribution n'est pas sphérique donnent donc lieu à un terme correctif pour la densité électronique de l'atome :

ρva(r)==0κ'3ρR,0(κr)m=p,my,m(θ,φ)

p,m est la population d'électrons contribuant au terme ,m de la somme.

Finalement, la densité électronique d'un atome non isolé s'écrit[11] :

ρ(r)=ρc(r)+κ3pvρv,0(κr)+=0κ'3ρR,0(κr)m=p,my,m(θ,φ).

Dans un cristal, la densité électronique d'un atome est soumise aux contraintes de symétrie imposées par le groupe ponctuel de la position de Wyckoff sur laquelle il se trouve. Ces contraintes indiquent quels termes ,m peuvent être présents dans l'expression de la densité électronique. Les paramètres κ, κ, pv et p,m peuvent être obtenus par affinement de structure si les données ont été collectées avec une très haute résolution.

Cas des neutrons

Il existe deux différentes interactions des neutrons avec la matière : avec le noyau des atomes et avec leurs moments magnétiques.

Interactions nucléaires

Graphique montrant les parties réelle et imaginaire de la longueur de diffusion b de l'isotope de l'europium 151Eu en fonction de l'énergie des neutrons incidents.
Parties réelle (vert) et imaginaire (bleu) de la longueur de diffusion b de l'isotope Modèle:ExpEu en fonction de l'énergie des neutrons incidents. La plage d'énergie des neutrons thermiques est indiquée en orange. Graphique effectué à partir des données publiées en référence[12].

La diffusion nucléaire d'un neutron libre par un noyau atomique s'effectue par le biais de l'interaction forte. La longueur d'onde des neutrons thermiques (quelques Å) et froids (jusqu'à quelques dizaines d'Å) généralement utilisée est de 4 à 5 ordres de grandeur supérieure à la taille du noyau (quelques fm). Les neutrons libres du faisceau incident se comportent comme une onde plane. Pour ceux qui sont soumis à une diffusion nucléaire, le noyau atomique agit comme une source ponctuelle secondaire et diffuse les neutrons sous forme d'onde sphérique (bien qu'il s'agisse d'un phénomène quantique, ceci peut être visualisé en termes classiques par le principe de Huygens-Fresnel). Dans ce cas ρ(r) représente la distribution de densité spatiale du noyau, qui est un point infinitésimalement petit comparé à la longueur d'onde des neutrons (distribution de Dirac, δ). La distribution de Dirac est contenue dans le pseudo-potentiel de Fermi Vf(r) qui décrit la diffusion d'un neutron par un noyau :

Vf(r)=h22πmnbδ(r)

h est la constante de Planck, mn la masse du neutron et b la longueur (ou amplitude) de diffusion[note 2]. La transformée de Fourier d'une distribution de Dirac est la fonction constante 1 : l'amplitude de diffusion est indépendante de l'angle de diffusion 2θ. Par abus de langage, il est souvent dit que les neutrons « ne possèdent pas de facteur de diffusion atomique ».

Comme l'interaction est nucléaire, chaque isotope possède une amplitude de diffusion différente. La transformée de Fourier est proportionnelle à l'amplitude de l'onde sphérique, qui a la dimension d'une longueur. Ainsi, l'amplitude de diffusion caractérisant l'interaction d'un neutron avec un isotope donné est appelée « longueur de diffusion » b. Elle est exprimée en femtomètres (fm). Les longueurs de diffusion neutroniques varient de façon erratique entre éléments voisins dans le tableau périodique et entre isotopes d'un même élément. Elles ne peuvent être déterminées qu'expérimentalement, car la théorie des forces nucléaires n'est pas adaptée pour calculer ou prédire la valeur de b à partir des autres propriétés du noyau[13]. Les valeurs des longueurs de diffusion de neutrons thermiques (de longueur d'onde entre Modèle:Unité et Modèle:Unité) pour chaque isotope sont données par le National Institute of Standards and Technology[14] ainsi que dans les tables internationales de cristallographie[15].

La longueur de diffusion b est fortement dépendante de l'énergie des neutrons incidents, en particulier dans le cas des neutrons chauds (d'énergie supérieure à Modèle:Unité, voir la figure). L'énergie des neutrons incidents nécessaire pour exciter le noyau est en général supérieure à celle utilisée dans les expériences de diffraction (entre 0,3 et Modèle:Unité pour l'isotope Modèle:ExpEu, par exemple), c'est la raison pour laquelle on observe peu d'effets de diffusion anomale avec les neutrons. Pour les quelques isotopes concernés, comme le Modèle:ExpB, b n'est pas une grandeur réelle mais est complexe :

b=b+ib

b décrit la variation en phase de la longueur de diffusion due à l'absorption de neutrons. Cette absorption donne ainsi lieu à de la diffusion incohérente. Dans la pratique, afin de simplifier l'analyse de données de diffraction de neutrons, il est possible de n'utiliser que des isotopes n'absorbant pas les neutrons lors de la synthèse du cristal, par exemple le Modèle:ExpB dans le cas du bore.

Interactions magnétiques

Bien que neutre, le neutron est un fermion et possède un spin. Lors de la diffusion de neutrons par de la matière condensée, la diffusion magnétique est due à l'interaction du spin du neutron avec les moments magnétiques (de spin et orbitaux) des électrons non appariés dans les orbitales externes de certains atomes : ce sont les électrons des sous-couches d pour les métaux de transition, 4f pour les terres rares et 5f pour les actinides. En général, une seule orbitale contribue à la diffusion magnétique. ρ(r) représente la distribution spatiale du moment magnétique de ces électrons non appariés autour du noyau. Comme les orbitales électroniques ont une taille comparable à la longueur d'onde des neutrons, le facteur de diffusion magnétique a une forme similaire à celui pour les rayons X. Cependant, cette diffusion magnétique n'a lieu qu'avec les électrons de valence, sans la contribution importante des électrons de cœur comme c'est le cas pour les rayons X.

Ainsi, au contraire de la diffusion nucléaire, l'objet diffusant lors de la diffusion magnétique est loin d'être une source ponctuelle et a même une taille effective plus grande que celle des objets diffusants pour les rayons X, ce qui conduit à une décroissance plus rapide du facteur de diffusion en fonction de sin(θ)/λ[16]. Une autre différence avec la diffusion nucléaire est que le facteur de diffusion magnétique ne dépend pas des isotopes mais du degré d'oxydation des atomes.

En définissant l'opérateur densité de magnétisation M^(r) qui donne le moment magnétique au point r par unité de volume, le facteur de diffusion magnétique fm d'un atome s'écrit :

fm(K)=q|0M^(r)eiKrdr3|qq|0M^(r)dr3|q

q| représente un état quantique d'un atome. Pour un vecteur de diffusion nul, le facteur de diffusion magnétique fm(0) est proportionnel au moment magnétique total de l'atome ou ion diffusant.

Graphique montrant la forme des intégrales radiales d'ordre 0, 2 et 4 pour le vanadium(IV) et le cuivre(II), à partir desquelles les facteurs de diffusion magnétique pour les neutrons peuvent être calculés
Intégrales radiales d'ordre 0 (rouge), 2 (noir) et 4 (bleu) pour le vanadium(IV) (lignes continues) et le cuivre(II) (lignes en pointillés), à partir desquelles les facteurs de diffusion magnétique pour les neutrons peuvent être calculés. Graphique effectué à partir des données publiées en référence[17].

Dans le cas où une seule orbitale, de symétrie sphérique et de moment orbital bloqué, est responsable de la diffusion magnétique, le facteur de diffusion magnétique peut s'écrire[18] :

fm(K)=J0(K)+(12g)J2(K)

g est le facteur de Landé et les fonctions J sont les intégrales radiales de la fonction d'onde radiale ψr(r) des électrons non appariés, définies par :

J(K)=4π0J(Kr)|ψr(r)|2r2dr

J étant la fonction de Bessel sphérique d'ordre .

Dans le cas où le moment orbital des électrons n'est pas bloqué, le facteur de diffusion atomique est :

fm(K)=J0(K)gSg+(J0(K)+J2(K))gLg

avec

gS=J(J+1)L(L+1)+S(S+1)J(J+1)gL=J(J+1)+L(L+1)S(S+1)2J(J+1)g=1+gS2

L étant le moment cinétique orbital, S le spin électronique et J le moment cinétique électronique total.

Si la symétrie de l'orbitale n'est pas sphérique, les intégrales radiales J d'ordre supérieur à 2 doivent être prises en compte.

Cas des électrons

Les électrons interagissent avec les potentiels électrique et magnétique des atomes dans la maille du cristal et le facteur de diffusion atomique électronique en est la transformée de Fourier[19]. Au contraire de la diffusion de rayons X ou de neutrons, la distribution de potentiel φ(r) n'est pas localisée autour d'un atome mais est une fonction continue de l'espace et dépend des autres atomes dans le cristal.

Il est cependant possible d'écrire la distribution de potentiel φ(r) comme étant la somme de distributions individuelles φi(r) des atomes aux positions ri :

φ(r)=iφi(rri).

Dans le cas des interactions électriques, qui sont les interactions prédominantes, la distribution de potentiel totale, périodique, est notée φe(r) et l'onde diffusée ψf(r) est donnée par l'équation de Schrödinger indépendante du temps :

2ψf(r)+8π2meh2[E+φe(r)]ψf(r)=0

m est la masse relativiste des électrons incidents, e la charge d'un électron et eE l'énergie cinétique de l'électron incident[20]. La correction relativiste de la masse des électrons est nécessaire, car le potentiel d'accélération des électrons utilisés en diffraction varie entre 10 et Modèle:Unité, donnant lieu à des vitesses v entre 20 % et 80 % de celle de la lumière, notée c. En tenant compte de cette correction, la masse de l'électron et sa longueur d'onde sont :

m=m01β2,λ=hm0c1β2β

m0 est la masse de l'électron au repos et β est le rapport v/c.

La résolution de l'équation de Schrödinger conduit à :

ψf(r)=ψi(r)+2πmeh20eik||rr||||rr||φe(r)ψf(r)dr

ψi représente l'onde incidente. Il s'agit d'une formule récurrente, puisque ψf apparaît dans l'expression de l'intégrale. L'expression du facteur de diffusion atomique pour la diffusion d'électrons à partir de l'expression donnée dans la définition générale est donc beaucoup plus compliquée que dans le cas de la diffusion de rayons X.

Approximation de Mott

Graphique montrant les facteurs de diffusion atomique de l'oxygène pour les électrons et pour les rayons X
Facteurs de diffusion atomique de l'oxygène pour les électrons (en noir) et pour les rayons X (en bleu). Graphique effectué à partir des données publiées en référence[21]Modèle:,[8].

Les facteurs de diffusion atomique électroniques, fe, peuvent être calculés à partir de ceux pour les rayons X, fX, en utilisant la formule de Mott[9], sous l'approximation d'atomes isolés et sphériques :

fe(K)=8π2me2h2ZfX(K)||K||2=me22h2(λsinθ)2(ZfX(K))

Z est le numéro atomique de l'atome ; fe est exprimé en Å. Cette formule prend en compte la diffusion élastique du nuage électronique ainsi que celle du noyau. Elle n'est plus valable aux très petits angles de diffusion : pour K=0, comme fX(0)=Z, le numérateur et le dénominateur du second terme de fe s'annulent et la fonction diverge. Dans ce cas, la formule d'Ibers[22] est utilisée :

fe(0)=4π2me23h2Zr2

r2 est la moyenne du carré du rayon de l'atome.

Une expression alternative des formules de Moth et d'Ibers laissant apparaître la permittivité du vide ϵ0 est aussi trouvée dans la littérature[23] :

fe(K)=2πme2h2ϵ0ZfX(K)||K||2fe(0)=πme23h2ϵ0r2

Interpolation

Les valeurs des facteurs de diffusion pour chaque atome et ion sont données pour sin(θ)/λ entre 0 et Modèle:Unité dans les tables internationales de cristallographie[21]. Ces valeurs peuvent être interpolées par une somme de cinq gaussiennes[24] :

fe(sinθλ)i=15aiebi(sinθλ)2

Les paramètres a1, b1, a2, b2, a3, b3, a4, b4, a5 et b5 sont également donnés dans les tables internationales de cristallographie.

Correction asphérique

Comme les rayons X, les électrons sont sensibles à la forme du nuage électronique des atomes via la distribution de potentiel φ(r). Les liaisons chimiques dans un cristal déforment la partie externe du nuage électronique, donnant lieu à une déviation de la sphéricité idéale. Dans ce cas, pour des angles de diffusion non proches de 0, la formule de Mott est toujours valable et est utilisée avec les valeurs de fX corrigées[25].

Notes et références

Notes

Modèle:Références

Références

Modèle:References Modèle:Traduction/Référence

Voir aussi

Articles connexes

Liens externes

Modèle:Portail


Erreur de référence : Des balises <ref> existent pour un groupe nommé « note », mais aucune balise <references group="note"/> correspondante n’a été trouvée