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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

CHAPITRE 6. EVALUATION DES PROPRIETES D’ADSORPTION DES FAUJASITESNaX et NaY VIS-A-VIS DU CO : SIMULATIONS MONTE CARLO DANS L’ENSEMBLEGRAND CANONIQUEceux issus <strong>de</strong> calculs DFT et <strong>de</strong> donnés <strong>de</strong> la littérature. Les résultats correspondants à cessimulations sont représentés dans le paragraphe suivant.IV. SIMULATIONS MONTE CARLO DANS L’ENSEMBLEGRAND CANONIQUE DE L’ADSORPTION DE CO DANS LESFAUJSITES NaX et NaYIV.1. MéthodologieLes isothermes d'adsorption ont été calculées via un co<strong>de</strong> <strong>de</strong> calcul développé aulaboratoire en utilisant l'algorithme Monte Carlo dans l’ensemble Grand Canonique (GCMC)qui permet le déplacement (translation et rotation), la création et la <strong>de</strong>struction <strong>de</strong>s espècesadsorbées. Ces simulations consistent à évaluer le nombre moyen <strong>de</strong> molécules adsorbéesdont le potentiel chimique est égal à celui <strong>de</strong>s molécules <strong>de</strong> la phase gaz pour une températureet une pression donnée.Les simulations ont été réalisées à 300 K pour un nombre <strong>de</strong> pas Monte Carlo (MC) <strong>de</strong>3.5*10 6 en utilisant la maille unitaire <strong>de</strong> chaque faujasite. Chaque pas MC peut correspondre àune translation, rotation, création ou annihilation d’une molécule <strong>de</strong> CO. On montre qu’enchoisissant ce nombre <strong>de</strong> pas MC, cela conduit à une convergence <strong>de</strong> l’énergie totale et dunombre <strong>de</strong> molécules adsorbée pour une pression donnée. Le réseau <strong>de</strong> la zéolithe a étéconsidéré rigi<strong>de</strong> eu cours <strong>de</strong> l’adsorption comme nous l’avons indiqué précé<strong>de</strong>mment.Deux cas ont été ensuite considérés pour les cations <strong>de</strong> sodium selon qu’on les considère fixesou mobiles. Je présente tout d’abord les résultats correspondant au cas <strong>de</strong>s cations fixes, quisont ensuite comparés à ceux obtenus dans le cas <strong>de</strong>s cations mobiles.La sommation d'Ewald a été utilisée pour calculer les interactions électrostatiques etl'interaction <strong>de</strong> van <strong>de</strong>r Waals a été calculée en utilisant un rayon <strong>de</strong> coupure <strong>de</strong> 12 Å.De plus, comme il a été montré expérimentalement que l’adsorption <strong>de</strong> CO ne peut seproduire que dans la supercage, <strong>de</strong>s atomes fictifs, ayant <strong>de</strong>s rayons <strong>de</strong> van <strong>de</strong>r Waalsappropriés (3 Å), ont été placés dans les cages sodalites pour éviter que les molécules <strong>de</strong> COaccè<strong>de</strong>nt à cette région <strong>de</strong> la faujasite, lui imposant ainsi <strong>de</strong> rester uniquement dans lessupercages. Le calcul <strong>de</strong>s enthalpies d’adsorption à 300K, h ads 0, a été effectué à faible taux217

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

Saved successfully!

Ooh no, something went wrong!