11.07.2015 Views

Zalfa NOUR Modélisation de l'adsorption des molécules à fort ...

Zalfa NOUR Modélisation de l'adsorption des molécules à fort ...

Zalfa NOUR Modélisation de l'adsorption des molécules à fort ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

CHAPITRE 2. NOTIONS FONDAMENTALES SUR LA MODELISATION MOLECULAIRE– THEORIE DE LA FONCTIONELLE DE LA DENSITE ET SIMULATIONS MONTE CARLOaugmentent considérablement les temps <strong>de</strong> calculs. Notons qu’une polarisation moyenne estincluse implicitement lors du processus <strong>de</strong> paramétrisation puisque les charges atomiques sontgénéralement choisies pour conduire à un moment dipolaire plus grand que celui observé pourune molécule isolée (par exemple, le moment dipolaire effective pour H 2 O dans l’état soli<strong>de</strong>est <strong>de</strong> 2.5 D comparé à 1.8D dans la phase gaz). Calcul <strong>de</strong>s interactions électrostatique – Sommation d’EwaldUn simple calcul direct du potentiel électrostatique nécessite <strong>de</strong>s temps <strong>de</strong> calcul qui<strong>de</strong>viennent rapi<strong>de</strong>ment insurmontables. C’est pourquoi, nous avons utilisé dans notre travailune métho<strong>de</strong> spécifique qui est la métho<strong>de</strong> <strong>de</strong>s sommations d’Ewald 52qui convertit lasommation classique en trois séries dont chacune converge plus rapi<strong>de</strong>ment que la sommationinitiale comme nous le verrons au fur et à mesure.La métho<strong>de</strong> d’Ewald consiste à ajouter à chaque charge ponctuelle q i du système unedistribution <strong>de</strong> charges , <strong>de</strong> même valeur mais <strong>de</strong> signe opposé à la charge q i , qui formeécraniun écran pour diminuer l’interaction entre les charges voisines. Cette distribution <strong>de</strong> chargeest exprimée à l’ai<strong>de</strong> d’une gaussienne <strong>de</strong> largeur centrée sur chaque charge ponctuelle q i :Le terme coulombienEécrani3qi2 2( r) exp( r )3Équation 802q qi, j 0 ijij4 rse décompose alors en une somme <strong>de</strong> trois termes :1- Le premier terme correspond à une sommation dans l’espace direct <strong>de</strong> l’énergieélectrostatique due à la <strong>de</strong>nsité <strong>de</strong> charge i (r) :1 qiqjE erfc(r)Équation 814 r1 ij0 i j i ijoù erfc(x) représente la fonction d’erreur complémentaire définie par :x22erfc ( x) exp( t) dt12Équation 82Pour <strong>de</strong> gran<strong>de</strong>s distances, erfc(rij) tend vers 0 et on constate que la sommation intervenantdans E 1 converge plus vite qu’une somme en r -1 . Le fait d’introduire <strong>de</strong>s distributions <strong>de</strong>charges gaussiennes a donc bien permis d’écranter les interactions à longue distance.76

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

Saved successfully!

Ooh no, something went wrong!