15.07.2013 Views

Ce document est le fruit d'un long travail approuvé par le jury de ...

Ce document est le fruit d'un long travail approuvé par le jury de ...

Ce document est le fruit d'un long travail approuvé par le jury de ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

3.3. Induction 65<br />

Dans l’approche <strong>de</strong> Jensen, la distance intervenant dans <strong>le</strong> tenseur d’interaction T ab<br />

tu <strong>est</strong><br />

modifiée suivant l’équation 2.21. Pour chaque type d’interaction, il faut donc déterminer <strong>le</strong><br />

<strong>par</strong>amètre a permettant <strong>de</strong> reproduire au mieux l’énergie <strong>de</strong> polarisation à courte distance.<br />

Dans <strong>le</strong> cas <strong>de</strong> la fonction d’atténuation <strong>de</strong> Tang et Toennies, l’énergie d’interaction<br />

<strong>est</strong> multipliée <strong>par</strong> la fonction fn(r) (voir équation 2.22). Le <strong>par</strong>amètre a du développement<br />

<strong>est</strong>, comme dans <strong>le</strong> cas <strong>de</strong> l’approche <strong>de</strong> Jensen, déterminé pour chaque interaction. Stone<br />

et collaborateurs [91, 92] ont proposé <strong>de</strong> déterminer ce <strong>par</strong>amètre à <strong>par</strong>tir <strong>de</strong>s énergies<br />

d’ionisation <strong>de</strong>s fragments en interaction. L’application <strong>de</strong> la fonction d’atténuation à<br />

l’énergie d’interaction d’induction (équations 3.11) donne :<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

U MM<br />

ind<br />

1 <br />

= 2 A B=A ∆Qat T ab<br />

tu fn(rab) 1/2Qb u<br />

∆Qa′ <br />

t ′ = −<br />

B=A αa′ a<br />

t ′ t<br />

fn ′(r) = 1 − exp(−ar) n ′ (ar)<br />

k=0<br />

k<br />

k!<br />

T ab<br />

tu fn ′(ra ′ b) 1/2 (Q b u + ∆Q b u)<br />

(3.12)<br />

Le facteur d’atténuation considéré dans cette approche revient à multiplier <strong>le</strong> tenseur<br />

d’interaction T ab<br />

tu <strong>par</strong><br />

<br />

fn(r ab). Stone et Millot [228] ont montré que la fonction d’atténu-<br />

ation <strong>de</strong> Tang and Toennies pouvait être raisonnab<strong>le</strong>ment tronquée à l’ordre n ′ = 6 2 . Un<br />

moyen pour éviter <strong>le</strong> calcul <strong>de</strong> la racine carrée <strong>de</strong> la fonction d’atténuation a été proposé<br />

dans un artic<strong>le</strong> <strong>de</strong> Meredith et Stone [229]. La version simplifiée <strong>de</strong> la fonction d’atténua-<br />

tion tronquée à l’ordre 6, f6(r), <strong>est</strong> :<br />

En conclusion, l’énergie d’induction modè<strong>le</strong> <strong>de</strong>vient :<br />

U MM<br />

ind = 1 <br />

2<br />

[f3(0.725b; r)] 2 ≈ f6(b; r) (3.13)<br />

A<br />

<br />

B=A<br />

∆Q a t T ab<br />

tu f3(β, rab)Q b u<br />

∆Q a′<br />

<br />

t ′ = − α<br />

B=A<br />

a′ a<br />

t ′ t T ab<br />

tu f3(β, ra ′ b)(Q b u + ∆Q b u)<br />

(3.14)<br />

où β = b/0.725 <strong>est</strong> <strong>le</strong> coefficient d’atténuation modifié <strong>de</strong> Tang et Toennies. Le terme<br />

f3(β, rab) s’écrit explicitement :<br />

f3(β, rab) = 1 − exp(−βrab)<br />

<br />

1 + βrab +<br />

(βrab) 2<br />

2<br />

+ (βrab) 3<br />

<br />

6<br />

(3.15)<br />

2. De ce fait, seuls <strong>le</strong>s calculs <strong>de</strong> puissances inférieures ou éga<strong>le</strong>s à 6 doivent être effectués. Il <strong>est</strong><br />

remarquab<strong>le</strong> que cel<strong>le</strong>s-ci soient déjà évaluées dans <strong>de</strong>s fonctions tel<strong>le</strong>s que <strong>le</strong> Lennard-Jones 12–6 ou <strong>le</strong><br />

Buckingham exp–6.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!