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

CHAPITRE 6. EVALUATION DES PROPRIETES D’ADSORPTION DES FAUJASITES NaX etNaY VIS-A-VIS DU CO : SIMULATIONS MONTE CARLO DANS L’ENSEMBLE GRANDCANONIQUEdans la faujasite NaY, nous avons considéré le potentiel décrivant l’interaction du CO avecles cations des sites II (Na + SII/CO). Alors que dans le deuxième, considérée pour décrirel’adsorption de CO dans la faujasite NaX, l’interaction Na + /CO a été modélisée à l’aide desparamètres correspondant aux deux potentiels Na + SII/CO et Na + SIII’/CO. Ces deux champs deforces ont été nommés respectivement CF1 et CF2.Le paragraphe suivant va être consacré à décrire les démarches que nous avons suiviespour la paramétrisation des deux potentiels Na + /CO. Comme ceci a été fait à partir des calculsDFT sur des clusters représentants la faujasite, la première étape était donc le choix de cesclusters, ainsi que de la méthode de calcul (fonctionnelle et base). Ces choix ont été validésensuite par une comparaison des paramètres géométriques et énergétiques des clustersoptimisés en présence de CO avec des données expérimentales et/ou théoriques antérieures. Apartir des configurations optimisées ainsi obtenus, l’étape suivante a consisté à construire lescourbes d’énergie potentielle pour les deux types d’interactions: Na + SII/CO et Na + SIII’/CO, etde les paramétrer à l’aide de fonctions analytiques qui peuvent être ensuite incorporées dansl’expression du champ de forces.III.1 Paramétrisation des potentiels décrivant les deux interactionsNa + SII/CO et Na + SIII’/CO à l’aide des calculs DFTIII.1.1 Modèles et méthodes utilisés dans les calculs DFTIII.1.1.A. ModèlesDeux clusters A et B constitués de 20 tétraèdres chacun ont été découpés à partir de lamême structure périodique de NaY utilisée dans le chapitre précédent 28, 36, 37 . Les liaisonspendantes ainsi générées dans chaque cluster ont été saturées à l’aide d’atomes d’hydrogène.Dans les deux modèles, tous les atomes d’aluminium ont été substitués aux atomes desilicium à part un seul choisi pour qu’il soit situé dans la fenêtre 6T à coté d’un cation Na + ensite II pour le modèle A et dans la fenêtre 4T à coté d’un cation Na + en site III’ pour lemodèle B (Figure 1).208