BORISSENKO
BORISSENKO
BORISSENKO
Contact : ddoc-theses-contact@univ-lorraine.fr
LIENS
Ecole doctorale
Energie Mécanique MAtériaux
Thèse
présentée pour obtenir le titre de
Docteur de l’Université Henri Poincaré en Sciences des Matériaux
par
Eléna Borissenko
Membres du jury :
Président
M. François Professeur, Nancy-Université
Rapporteurs
A. Boutin Directrice de recherche au CNRS, Université Paris Sud
A. Simon-Masseron Maître de Conférences, Université de Haute Alsace
Examinateurs
E. Dooryhée Chargé de recherche au CNRS, Grenoble
M. Souhassou Professeur, Nancy-Université Directeur de Thèse
F. Porcher Maître de Conférences, Nancy-Université, Co-directrice
C. Lecomte Professeur, Nancy-Université Invité
2
à Gildas qui comprend tant d’autres choses
et
3
4
Remerciements
Je veux remercier János Ángyán (Chargé de recherche au CNRS, LCM3B) pour son
implication dans la partie modélisation de cette thèse, sa disponibilité et son aide dans les
5
contacts avec d’autres laboratoires. Je remercie William Nicolazzi (Docteur, LCM3B), pour
nos nombreux échanges et son assistance dans les analyses de calculs Monte-Carlo, ainsi que
Emmanuel Aubert (Docteur, LCM3B) et Enrique Espinosa (Professeur, LCM3B) pour nos
diverses discussions, scientifiques et autres.
Ma reconnaissance va à ma famille pour son appui durant toutes ces années, mes
parents et ma soeur, Hubert et Françoise, Pat, Fab et tous mes amis. Je remercie
chaleureusement Hubert Gaborieau (Agrégé des Lettres), pour les relectures attentives de mon
texte.
Tout particulièrement, j’aimerais remercier mon mari, Gildas, qui m’a soutenue,
supportée aux moments de doute, et sans qui je ne serais pas ce que je suis aujourd’hui.
6
Table des matières
7
III.4.1 Interprétation qualitative des spectres XANES................................................ 97
III.4.2 Spectre EXAFS des cristaux de Na,Co-X déshydratés. ................................... 99
III.4.3 Spectre EXAFS de poudre Na,Co-LSX déshydratée. .................................... 112
III.4.4 Spectre EXAFS de poudre Na,Co-LSX réhydratée. ...................................... 115
III.4.5 Bilan des analyses EXAFS............................................................................. 119
III.5 Simulations Monte-Carlo de la distribution cationique. .................................. 121
III.5.1 Objectifs, originalité des calculs. ................................................................... 121
III.5.2 Champs de forces : interactions entre les espèces.......................................... 122
III.5.3 Hypothèses de départ des simulations............................................................ 124
III.5.4 Géométrie de la charpente et stabilité des sites.............................................. 124
III.5.5 Distribution des cations dans la structure déshydratée................................... 125
III.5.6 Etude de la distribution des cations et H2O dans la structure hydratée.......... 127
III.5.7 «Dynamique» de la redistribution cationique lors de l'hydratation................ 131
III.6 CONCLUSIONS................................................................................................... 134
III.7 Composite Na,Co(bpy)3-X....................................................................................... 136
IV Etude du système (dmpNA-mordénite) par Diffraction des rayons X sur poudre......... 143
IV.1 Introduction .......................................................................................................... 143
IV.2 Structure de la mordénite : Généralités............................................................. 145
IV.3 Structure de la mordénite par DRX sur monocristal ....................................... 147
IV.3.1 Conditions expérimentales. ............................................................................ 148
IV.3.2 Structure de la mordénite déshydratée ........................................................... 149
IV.3.3 Discussion ...................................................................................................... 152
IV.4 Composite (dmpNA-mordénite) : Diffraction des rayons X sur poudre......... 155
IV.4.1 Conditions expérimentales ............................................................................. 155
IV.4.2 Spectre à 90 K ................................................................................................ 155
IV.4.3 Spectre à 300 K .............................................................................................. 162
IV.4.4 Discussion. ..................................................................................................... 165
IV.4.5 Diffractogramme de mordénites synthétisées avec et sans alcool. ................ 171
IV.5 CONCLUSION..................................................................................................... 174
V Conclusions et Perspectives .......................................................................................... 177
Perspectives............................................................................................................................ 179
ANNEXES ............................................................................................................................. 183
ANNEXE A : Rappels sur la Diffraction des Rayons X ............................................... 183
ANNEXE B : Approche théorique de l’absorption des rayons X................................ 193
ANNEXE C : Principe des simulations Monte-Carlo. .................................................. 205
Bibliographie ......................................................................................................................... 213
8
Chapitre I
« Introduction »
9
10
I Introduction
Les zéolithes, naturelles ou synthétiques, sont des aluminosilicates poreux
relativement stables thermiquement et chimiquement, et dont la géométrie des cavités et
canaux est parfaitement régulière. Leur charpente composée de tétraèdres TO4 (T=Si4+, Al3+)
porte une charge anionique qui est compensée par des cations situés dans les canaux, libres et
facilement échangeables. Grâce au champ électrostatique important dans les cavités et leur
géométrie régulière, les zéolithes peuvent adsorber des molécules hôtes avec une sélectivité
qui peut-être modifiée (contrôlée) en changeant le cation compensateur de charge. Ainsi, elles
sont actuellement largement utilisées dans les procédés d’échanges cationiques, pour la
séparation de gaz et également en catalyse.
11
ont été étudiés comme antennes artificielles pour le transport et la conversion d’énergie
lumineuse (Gfeller et al., 1998; Megelski at al., 2001; Ranjit & Kevan, 2001), pour la
réalisation de nouvelles générations de LED (Brühwiler & Calzaferri, 2004) ou comme milieu
lasant (Leike & Marlow, 1996; Reck et al., 1996, Simoncic & Armbruster, 2005). Elles ont
également été envisagées pour l'élaboration de commutateurs magnétiques pour le stockage
d’information (Co2+(bipyridine)3-zéolithe Y, Mizuno & Lunsford, 1983; Tiwary &
Vasudevan, 1998).
Dans tous ces cas, à la base des propriétés macroscopiques recherchées se trouve
l'interaction entre les espèces organiques invitées et la zéolithe (charpente, cations
extracharpente), encore mal connue, en particulier du point de vue structural.
12
Chapitre II
13
14
II Préparation des échantillons........................................................................................... 17
II.1 Etat de l’art ............................................................................................................. 17
II.1.1 Gel .................................................................................................................... 17
II.1.2 Nucléation ........................................................................................................ 18
II.1.3 Croissance ........................................................................................................ 19
II.1.4 Rôle des différents paramètres modifiant la croissance ................................... 21
II.1.5 Bilan-Conclusion.............................................................................................. 25
II.2 Zéolithe X ................................................................................................................ 27
II.2.1 Synthèse des cristaux de zéolithe X. ................................................................ 27
II.2.2 Optimisation du protocole ................................................................................ 29
II.2.3 Poudre de zéolithe X ........................................................................................ 31
II.2.4 Echange cationique .......................................................................................... 31
II.2.4.1 Capacité d’échange cationique en solution ................................................... 32
II.2.4.2 Échange cationique en colonne. .................................................................... 33
II.2.4.3 Échange cationique en bain. .......................................................................... 34
II.2.5 Caractérisations. ............................................................................................... 35
II.2.5.1 Analyses chimiques. ...................................................................................... 35
II.2.5.2 Analyse thermiques ATG/ATD..................................................................... 38
II.2.6 Déshydratation. ................................................................................................ 40
II.2.7 Préparation du composite Co(bpy)3-X. ............................................................ 41
II.2.8 Conclusion........................................................................................................ 41
II.3 Mordénite................................................................................................................ 43
II.3.1 Morphologie et taille des cristaux .................................................................... 43
II.3.2 Défauts structuraux .......................................................................................... 43
II.3.3 Propriétés d’adsorption/diffusion ..................................................................... 44
II.3.4 Protocoles de synthèse de mordénite. .............................................................. 44
II.3.4.1 Protocole 1 ..................................................................................................... 45
II.3.4.2 Protocole 2 ..................................................................................................... 46
II.3.4.3 Protocole 3 ..................................................................................................... 47
II.3.4.4 Synthèses à partir de verre (Gilbert & Mosset) ............................................. 48
II.3.4.5 Bilan des synthèses........................................................................................ 48
II.3.5 Caractérisation.................................................................................................. 49
II.3.5.1 Analyse chimique .......................................................................................... 49
II.3.5.2 Analyses thermiques : ATG et ATD. ............................................................ 49
II.3.6 Calcination, déshydratation. ............................................................................. 52
II.3.7 Préparation de composites. ........................................................................... 53
II.3.7.1 Pérylène – mordénite. .................................................................................... 54
II.3.7.2 Thio-indigo – mordénite. ............................................................................... 55
II.3.7.3 dmpNA-mordénite......................................................................................... 56
II.3.8 Conclusion........................................................................................................ 56
15
16
II Préparation des échantillons
Les premières zéolithes synthétiques ont été cristallisées par l’équipe de R. Barrer et
R. Milton dans les années 50 mais on tente toujours, depuis, de comprendre et modéliser leurs
processus de cristallisation. Cette tâche est rendue difficile, d’une part, par la complexité du
gel où une grande variété de réactions chimiques est possible, et d’autre part par le grand
nombre de variables (temps, température, pH, pression) modifiant la cinétique de croissance
cristalline. Pour cette raison, des protocoles de synthèse déjà établis ont été utilisés au cours
de ce travail et optimisés de façon à obtenir des cristaux de qualité et de taille convenable
pour des études par DRX sur monocristal au laboratoire. Le but de ce chapitre est de présenter
brièvement l’état de l’art sur la cristallisation des zéolithes afin de déterminer les paramètres
jouant un rôle important lors des synthèses visant à produire des monocristaux, puis de décrire
les protocoles mis en pratique lors de la thèse.
II.1.1 Gel
Le gel de synthèse de zéolithes est composé d’une phase solide (cristalline, amorphe)
et d’une phase liquide. La phase solide comporte des clusters, des précipités, des germes, des
cristaux et des espèces colloïdales (agglomérats amorphes). La phase liquide est une solution
aqueuse enrichie d’ions (alumino)silicatés plus ou moins polymérisés, d’ions hydroxyles OH-
ou fluorures F-, de cations (Na+, …) et du template. Les espèces amorphes sont dispersées
17
dans la solution et jouent le rôle de sources de nucléation et de nutriments pour la croissance
des cristaux (Mintova et al., 1999). L’un des grands problèmes de la compréhension de la
chimie des zéolithes réside dans la difficulté à distinguer les phases liquide (espèces
polymérisées) et amorphe (colloïde), donc de quantifier la concentration de la phase amorphe
(Cundy & Cox, 2005).
Réactions de condensation/hydrolyse
Le processus de cristallisation des zéolithes est rendu possible par les réactions de
condensation (formation de liaisons T-O-T, T = Al, Si) et d’hydrolyse (rupture des liaisons T-
O-T) des espèces silicatées et aluminosilicatées. Pour la silice, par exemple, la réaction
d’hydrolyse en milieu hydroxyde s’écrit comme suit :
II.1.2 Nucléation
Initialement, lors du mélange des réactifs, se forme un gel inhomogène, où les espèces
chimiques n’interagissent que partiellement. Sous l’effet de la température, du temps et de
l'agitation, des réactions de dissolution/condensation rendent possible l'établissement
progressif d'un équilibre entre les espèces silicatées et aluminosilicatées, résultant en une
redistribution équilibrée des produits de réaction entre les phases solide et liquide. A ce stade,
la phase liquide du gel contient de nombreuses espèces ioniques avec une large distribution de
charges/masses. La polymérisation progressive des entités silicatées et aluminosilicatées
18
modifie les propriétés du gel (concentration des espèces, viscosité). Les cations jouent un rôle
structurant: ils s’entourent d’anions (polymères silicatés et aluminosilicatés) pour former une
sphère de coordination énergiquement favorable en induisant la géométrie de clusters
particulièrs (Breck, 1974; Burkett & Davis, 1994; Brunner, 1992). Une phase pré-cristalline
localement ordonnée (nucleus) se forme. La croissance de certains nuclei au-delà de la taille
critique donne naissance à des nuclei stables (germes), précurseurs des cristaux de zéolithe.
Dans le cas des zéolithes faisant l'objet de cette thèse (zéolithe X, mordénite) on
s'attend à une nucléation primaire hétérogène à partir des espèces colloïdales issues de la
phase amorphe (Mintova et al., 1999; Cundy & Cox, 2005;). Les conditions de mûrissement
du gel influencent la taille critique des nuclei et vont être explorées pour les synthèses de gros
cristaux.
II.1.3 Croissance
Dès le début des études sur la croissance des zéolithes synthétiques, Barrer et al.
(1959) ont fait l’hypothèse que la croissance d’une structure poreuse n’est pas favorable si
elle doit se faire par agrégation successive de monomères TO4, T = Si, Al, de type « Primary
Building Unit » sur le nucleus. Ils ont donc envisagé que les unités constructives des zéolithes
intervenant lors des phases de nucléation/croissance sont en fait des anions plus complexes
(briques secondaires, «Secondary Building Unit (SBU)») comme, par exemple, des anneaux
composés de 4, 5, 6 tétraèdres TO4 (4MR, 5MR, 6MR), des anneaux doubles de 4 ou 6
tétraèdres (D4R, D6R)... Cette hypothèse est appuyée par des études de spectroscopie UV-
Raman sur des phases solides et liquides (Yu et al., 2001; Xiong et al., 2001) qui font état de
la présence de pics de vibration attribués aux liaisons Si-O-Al formant des entitées 4MR au
début de la croissance de zéolithes X. Ainsi, au cours de la synthèse des zéolithes A, X et
EMT, les anneaux 4MR s’interconnectent pour former des anneaux 6MR et des cages β. Ces
conclusions sont également supportées par des études de spectroscopie RMN (Thangaraj &
Kumar, 1990) et AFM (Anderson et al., 2001; Wakihara et al., 2004), mais sont en
contradiction avec les mesures RMN 29Si et 27Al de Shi et al. (1996) qui n’observent que les
fréquences dues aux monomères TO4 (T=Si, Al).
Dans le cas de la mordénite, Yu et al. (2001) et Twu et al. (1991) observent par
spectroscopie UV-Raman et IR-Raman respectivement les fréquences de vibrations attribuées
19
aux liaisons T-O-T des anneaux à 4 et 5 chaînons.
Mécanisme de croissance
Pour qu’un germe croisse, des unités constructives doivent pouvoir diffuser vers sa
surface, puis s’y intégrer suivant différentes directions cristallographiques. La vitesse de
diffusion et les mécanismes d’intégration des unités constructives à la surface sont donc des
paramètres limitant la croissance cristalline. Yang et Navrotsky (2002) dans leurs études
calorimétriques in situ de la cinétique de cristallisation des FAUjasites ont montré que
l’énergie d’activation pour la cristallisation (~ 66 - 72 kJ/mol) correspond typiquement aux
énergies de rupture et formation des liaisons T-O-T, et est trop grande pour être attribuée à la
diffusion des nutriments (~ 16 - 17 kJ/mol) (Barrer, 1981). La croissance des zéolithes semble
donc être principalement limitée par la vitesse d’intégration des blocs élémentaires à la
surface des nuclei.
∗
« Terrace-Ledge-Kink » (TLK) est un modèle thermodynamique de la formation d’une surface cristalline et
de l’introduction des défauts (Kossel, 1927 ; Stranski, 1928). Le principe : l’énergie de l’incorporation d’une
entité (atome, SBU) à la surface dépend du nombre de liaisons possibles avec les atomes de la surface ; sites
TLK (Figure ci-dessus) offrent la coordination différente.
20
Vitesse de croissance
Il est communément admis (Sun et al., 1995; Warzywoda et al., 1996; Yang et al.,
2006) que la croissance des gros cristaux n'est obtenue qu'en ralentissant la vitesse de
croissance et en limitant le nombre de germes. Pour ce faire, on essaye de séparer les phases
de nucléation et de croissance. On peut alors ralentir ensuite l’intégration des unités
constructives aux cristallites par :
Température, temps
21
que le rapport d’aspect longueur/largeur des cristaux formés (Feoktistova et al., 1989; Cundy
et al., 1993).
Les zéolithes sont des matériaux métastables et par conséquent, la cinétique jouent un
rôle crucial lors de la cristallisation. La formation des zéolithes obéit à la loi des
transformations successives d’Ostwald selon laquelle la première phase qui apparaît est
thermodynamiquement moins stable que celle qui la remplace (la précipitation initiale du gel
et sa transformation en zéolithe peuvent-être considérées comme un exemple des réactions
successives). Lorsque le temps de cristallisation augmente, la formation d’une phase dense
(plus stable) est la plus favorable (Lewis et al., 1995; Francis & O’Hare, 1998).
Homogénéisation, Mûrissement
Agent minéralisant, pH
La majorité des zéolithes est synthétisée en milieu fortement basique. Les ions OH-
catalysent alors les réactions en rendant les réactifs de départ mobiles. Ils favorisent les
interactions entre les différents réactifs en influençant leur réactivité chimique ainsi que les
réactions de condensation/hydrolyse lors de la formation des zéolithes (Cundy & Cox, 2005).
Depuis la fin des années 70, les anions fluorure sont utilisés comme agents
minéralisants dans les synthèses de certaines zéolithes siliceuses (Zones et al., 2005); donc en
milieux moins basiques (neutre/acide). L’utilisation des anions F- a pour avantage de diminuer
22
la vitesse de dissolution des réactifs et au final, d'augmenter considérablement la taille des
cristaux.
• le template au sens strict, ayant une géométrie autour de laquelle la charpente adapte
sa géométrie et sa configuration électronique (e.g. une molécule organique; un même
template peut conduire à la cristallisation de différents types de zéolithes),
Les synthèses des zéolithes peu siliceuses avec Si/Al<3-5 (« Low Silica Zeolite »
comme les zéolithes X) sont réalisées en présence de cations alcalins (sodium ou potassium)
comme SDA (Brunner, 1992; Cundy & Cox, 2005). Le rapport Si/Al de la zéolithe augmente
avec la diminution de l'alcalinité du gel de synthèse (Lechert et al., 1996). Actuellement, la
possibilité d'obtenir de nouvelles zéolithes peu siliceuses est limitée par le choix des agents
structurants car il n'est pas encore possible d'en synthétiser en présence de cations autre que
des cations alcalins.
Dans le cas des zéolithes siliceuses comme la mordénite, on a recours en général à des
templates organiques. La grande diversité de ces additifs a d'ailleurs rendu possible les
synthèses de nombreuses nouvelles phases zéolithiques (Zones & Hwang, 2002; Burton et al.,
2004).
Les interactions entre ces additifs organiques et la charpente ont fait l’objet d’études
théoriques et expérimentales (mais rarement cristallographiques). Lewis et al. (1995) et
Rollmann et al. (2000) concluent ainsi que les interactions template-charpente sont des
interactions de type « Van der Waals » d'où la complémentarité de la charpente et du
template. Rollmann et al. (2000) ont montré que, même si les molécules organiques ne sont
pas piégées dans le cristal en fin de synthèse, elles jouent un rôle stabilisateur pour les
structures poreuses.
23
Enfin, certaines espèces organiques (amines quaternaires) sont réputées complexer les
espèces alumineuses (Morris et al., 1991) dans le gel ralentissant ainsi leur intégration dans
des unités constructives du cristal et diminuant la vitesse de croissance.
Viscosité
L’eau, utilisée comme solvant polaire (ε=87ε0 avec ε, permittivité électrique et ε0,
permittivité électrique absolue), permet également de contrôler la viscosité, les concentrations
absolues en réactifs donc la vitesse de diffusion des nutriments et les réactions entre les
précurseurs silicatés et aluminosilicatés. Ce rôle est également partagé par les additifs
organiques. Les zéolithes sont généralement synthétisées sous pression autogène, pression qui
peut être modifiée par des vapeurs des espèces organiques.
Alors que la composition chimique des réactifs utilisés est clairement établie par le
fabricant, leur qualité (réactivité, pureté) est en général peu ou pas du tout renseignée.
Les principales sources d'aluminium sont des sels facilement solubles dans l'eau
(Al(NO3)3⋅6H2O, AlCl3⋅6H2O, Al2(SO4)3⋅6H2O, AlNaO2). L’utilisation d'oxyanions différents
permet d’obtenir des cristaux de tailles et morphologies différentes (Sun et al., 1995; Kumar
et al., 1998; Lu et al., 2005). Les sources de silicium et d'aluminium sont souvent riches en
impuretés qui peuvent être des sources de nucléation non contrôlée.
24
II.1.5 Bilan-Conclusion
Les synthèses de zéolithes étant des processus complexes, nous nous sommes basés
sur des protocoles déjà établis par Charnell (1971) et Sano et al. (2001) pour des zéolithes X
et mordénite respectivement. Ces protocoles ont été adaptés dans le but d'obtenir des cristaux
de grande taille et de bonne qualité pour nos études par Diffraction des Rayons X sur
monocristal. Dans ce but, ont été testés les sources de réactifs (silicium/aluminium), les
concentrations en additifs organiques, l'ordre de mélange des réactifs, l'influence des
paramètres température/temps de mûrissement sur le produit final et l'effet du mûrissement
séparé des solutions alumineuses et siliceuses. L’influence du temps de cristallisation sur la
taille et la morphologie des cristaux à une température fixe a également été examinée.
25
26
II.2 Zéolithe X
Protocole de synthèse
Le gel de synthèse est formé par deux solutions de silicate et d’aluminate préparées en
parallèle. Les proportions des réactifs sont calculées de manière à remplir le flacon de 60 ml.
27
soude (Carlo Erba, pureté 99 %) est dissoute dans 1,11 mol d’eau distillée, puis 0,011 mol de
silice pulvérulente (Merck, pureté 99 %) est ajoutée. On laisse mûrir la solution pendant un
temps Tm pour obtenir à la fin une solution transparente.
La solution silicatée (A) est transvasée dans un flacon en polypropylène avant qu’on y
introduise 0,085 mol de TEA. Après une brève homogénéisation manuelle de la solution, la
solution B est versée dans le flacon, puis agitée de nouveau manuellement : le gel formé
devient blanc opaque. Le flacon est placé immédiatement dans un bain thermostaté à
température de cristallisation (85°C).
Figure 2.1 : Image au microscope optique des cristaux de zéolithes X (octaèdres), P (sphères), A (cubes).
28
II.2.2 Optimisation du protocole
Temps de cristallisation
L'influence du temps de cristallisation sur la taille des cristaux a été testée en réalisant
7 synthèses identiques (même composition et protocoles de préparation). Dans tous les cas, le
gel initialement opaque décante au bout de 14 à 21 jours en deux phases : solide blanche et
liquide jaunâtre. A la fin de la décantation les cristaux de zéolithe X mesurent typiquement 60
à 100 µm. Contrairement à ce qui était attendu, la croissance des cristaux ne s'arrête pas à la
fin de la décantation. La Figure 2.2 montre que les cristaux de zéolithes X croissent en deux
étapes, une première étape quasi-linéaire et assez rapide jusqu’à la décantation, puis une étape
également linéaire mais beaucoup plus lente.
220
200
180
160
Taille, µm
140
120
100
80
60
40
10 20 30 40 50 60 70 80 90 100
Temps, jour
On peut supposer que la croissance dans la première étape se fait par consommation des
nutriments de la phase liquide du gel (croissance rapide) puis dans la seconde étape par
consommation des nutriments issus de la dissolution progressive des particules, cristallites,
agglomérats (croissance lente).
29
Tableau 2.1: Influence de la température de mûrissement du gel sur le produit de synthèse (la quantité des
cristaux augmente avec le nombre de signes +)
Ce Tableau 2.1 montre clairement que la taille des cristaux et la phase dominante de
synthèse dépendent de la température de mûrissement. Nous constatons ici que le
mûrissement à 4°C permet d’éliminer la phase "zéolithe X", tandis que l'augmentation de la
température inverse la tendance : augmentation de la quantité de cristaux de zéolithe X et
diminution des zéolithes A. La température optimale de mûrissement est de 60°C pour obtenir
des cristaux de zéolithe X de grande taille.
Nous avons entrepris une série de synthèses ne différant que par la quantité d’agent
complexant (TEA) : 100% TEA, 120% TEA, 135% TEA, en se basant sur le protocole décrit
dans le Chapitre II.2.1 et avec un mûrissement séparé des solutions siliceuse et aluminatée
pendant 24 heures à 4°C.
Tableau 2.2 : Influence de la quantité de TEA sur la taille et la quantité des cristaux (la quantité des cristaux
augmente avec le nombre de signes +)
N Cristallisation Produit
TEA (%) Temps (jour) X Taille (µm) A Taille (µm) P Taille (µm)
1 100 53 +++ 100-160 ++ 40-70 ++ 20-60
2 120 54 ++ 100-160 + 20-50 ++++ 20-60
3 135 53 + 120-175 - - ++++ 20-60
Le Tableau 2.2 montre que l’augmentation de la quantité de TEA dans les synthèses a
pour effet principal de changer le rapport entre les quantités des phases A, X et P dans le
produit final. Ainsi, dans le cas des synthèses avec 100% de TEA, la phase X est dominante,
avec de grands cristaux, et une quantité comparable de zéolithes A et P. Pour 120 % de TEA,
30
la phase P prédomine et on a moins de zéolithe X et très peu de zéolithe A. Les synthèses
avec 135% TEA donnent beaucoup de P, très peu de zéolithe X (moins que dans des
synthèses précédents) et ne donnent pratiquement pas de zéolithe A.
Pour les études structurales, nous avons utilisé les gros cristaux de zéolithe X (Figure
2.3) synthétisés selon le protocole décrit dans le Chapitre II.2.1 (mûrissement à 60°C, 100 %
TEA). La synthèse a duré 90 jours et subi des variations incontrôlées de température (2 arrêts
du four : le premier au bout de 4 jours, le deuxième, au bout de ~60 jours avec une baisse de
température de 85°C à l'ambiante); les cristaux sont ensuite restés dans le flacon avec le
surnageant jaunâtre pendant 30 jours supplémentaires à la température ambiante.
Nous avons également synthétisé la zéolithe LSX (Low Silica zeolite X) avec Si/Al~1
sous forme de poudre pure en utilisant le protocole de Kühl («Verified Syntheses of Zeolitic
Materials» Robson, ed. Lillerud, 2001, p. 153). Ces synthèses ont été utilisées pour les études
EXAFS nécessitant une phase pure de zéolithe X.
Les zéolithes X sont synthétisées avec des cations Na+ (ou Na+ et K+), et toutes les
autres formes sont obtenues par échange cationique. L’échange cationique en solution est la
technique la plus développée et la plus efficace comparée à l'échange par voie sèche (échange
à l’état solide): l'échange 2Na+ → Co2+ dans la zéolithe X est de ~85% s’il est effectué en
solution comparé à ~50% s’il est réalisé à l’état solide (Feddag & Bengueddach, 2005).
31
II.2.4.1 Capacité d’échange cationique en solution
Pendant l’échange, les cations se trouvant dans les pores des zéolithes diffusent vers la
solution et les cations de la solution diffusent pour entrer dans la zéolithe. Cette diffusion
résulte de la différence de concentration entre les cations dans la solution et le cristal. La
tendance des cations à diminuer cette différence de concentrations est limitée par le champ
électrique généré par la charpente, l’affinité et l’accessibilité des sites de la charpente aux
nouveaux cations (Helfferich, 1962). L’affinité des sites dépend de la force d’interaction des
cations entrants avec les molécules d’eau et les oxygènes de la charpente (Maes & Cremers,
1975), elle varie lors de l'échange avec la composition chimique de la zéolithe (l’affinité d’un
site est affectée par la présence d’un cation sur un site voisin).
32
ions hydratés (4,12 et 4,23 Å pour Ca2+ et Co2+ respectivement) mais pourrait résulter d’un
équilibre lié au rapport charge/taille (donc polarisabilité) des cations, aux interactions
électrostatiques entre les cations et la charpente et à la composition de la zéolithe (l'affinité de
la charpente aux cations change au cours de l'échange cationique).
33
Tableau 2.3 : Conditions expérimentales d’échange cationique en colonne
Les cristaux de zéolithes X récupérés sont de couleur rose et ont conservé leur
morphologie octaédrique. Le fritté est partiellement obstrué par un dépôt verdâtre attribué à
des précipités d'hydroxyde de cobalt. La diffraction des rayons X montre en fait
l’amorphisation complète des cristaux de zéolithes X et A bien que leur morphologie ait été
préservée.
Cette méthode a été utilisée par Bae et Seff (1999, 2000, 2000) pour des échanges des
cations sodium natifs par des métaux de transition. Ils ont obtenu des échanges partiels pour
Co2+ (Bae & Seff, 1999) et complets pour Ni2+ (Bae & Seff, 2000) et Zn2+ (Bae & Seff, 2000).
Cependant les auteurs ont également constaté la perte de la cristallinité lors de la
déshydratation, l’échange avec des H3O+, OH- et l’apparition d'aluminium extracharpente. Par
conséquence, cette technique a été abandonnée au profit de l'échange cationique en bain.
Dans cette méthode, 0,5 g de cristaux de zéolithe Na-X, Na-A et P sont plongés dans
une solution de nitrate de cobalt (Co(NO3)2, Aldrich) 0,02 molaire (pHdébut = 5,6) chauffée à
60°C sous agitation. Au bout de 2 jours d'échange, les cristaux sont filtrés, rincés
abondamment avec ~500 ml d'eau distillée à température ambiante puis séchés à 70°C. Les
cristaux roses sont accompagnés de résidus verdâtres d'hydroxydes de cobalt. En fin
d’échange, le pHfin de la solution est de 6,5. Cette procédure est répétée 5 fois avec des
solutions fraîchement préparées pour une durée totale d’échange de 13 jours car les analyses
chimiques par microsonde ont montré que seuls 72% des cations sodium sont échangés lors
du 1er échange. Après les 5 échanges en bain, la proportion d’ions Na+ résiduels tombe à
~15% et la couleur des cristaux fonce et devient bordeaux (Figure 2.5). Ces cristaux ont été
utilisés ensuite pour les études structurales.
34
Figure 2.5: Image au microscope optique des cristaux Na,Co-X échangés en bain
II.2.5 Caractérisations.
Les mesures de composition chimique ainsi que les images de contraste chimique sont
effectuées au service commun de l’Université Nancy-I sur une microsonde CAMECA SX50,
utilisée à la fois pour l'imagerie et pour la microanalyse (mode WDS). Les éléments plus
lourds que le fluor peuvent être dosés.
Dans le cas des zéolithes, la précision de l'analyse élémentaire est limitée par
l'échauffement des zéolithes sous le faisceau d'électrons qui produit un départ brusque des
éléments les moins liés (cations ou molécules d'eau) en amorphisant la structure. Pour limiter
ces dégâts, les mesures sont faites en limitant l'intensité du faisceau incident. A un moindre
degré, la précision de l'analyse est également limitée par l'utilisation des standards (oxydes,
e.g.) où les éléments à doser ont une coordination différente de celle des zéolithes. Nos
analyses ont été effectuées sous un courant de 10 nA et une tension de 15kV.
35
Cristaux Na,Co-X échangés en bain pendant 2 jours
Les dosages élémentaires sur 3 cristaux différents (1 point au centre du cristal) sont
rassemblé dans le Tableau 2.4.
Tableau 2.4 : Composition par maille des cristaux Na,Co-X échangés pendant 2 jours
Cristal Na Co Si Al O Si/Al
1 30 27,9 99 93 381 1,06
2 22 27,5 98 94 376 1,04
3 27 27,7 101 91 380 1,11
moyenne 26 27,7 99 93 379 1,07
(La teneur en oxygène est ajustée pour que les éléments dosés forment des oxydes)
Les images de contraste chimique (Figure 2.6) des cristaux issus de l'échange pendant
13 jours en bain montrent que certains cristaux paraissent plus clairs que d'autres. Le Tableau
2.5 donne les résultats des analyses élémentaires : la composition chimique est une moyenne
des mesures effectuées sur 5 cristaux « clairs » (5 points par cristal) et 4 cristaux sombres (5
points par cristal).
36
Tableau 2.5 : Composition par maille des cristaux Na,Co-X (bain, 13 jours d'échange)
Cristal Na Co Si Al O Si/Al
Claire 1 8(1) 34,0(8) 100(7) 93(2) 376 1,08
2 9(1) 32,7(7) 100(7) 93(2) 375 1,08
3 13(2) 31,3(7) 99(7) 93(2) 375 1,06
4 14(2) 31,1(7) 99(7) 93(2) 375 1,06
5 8(2) 31,2(7) 99(7) 93(2) 373 1,05
Moyenne 10(2) 32,1(7) 99(7) 93(2) 375 1,07
Sombre 1 16(3) 33,8(8) 99(7) 93(2) 381 1,06
2 17(3) 32,6(8) 100(7) 92(2) 378 1,08
3 16(3) 32,0(8) 100(7) 92(2) 380 1,09
4 19(3) 31,6(8) 100(7) 92(2) 380 1,10
Moyenne 17(3) 32,5(8) 100(7) 92(2) 380 1,08
Les incertitudes sont très importantes sur le dosage du sodium (~20%), elles sont
également élevées pour les mesures de silicium (~7%) et seulement ~2% pour les dosages du
cobalt et de l'aluminium. D'après l''analyse la seule différence entre les cristaux sombres et
clairs est la teneur en sodium (Figure 2.6) qui est quasiment deux fois plus élevée dans les
cristaux sombres que dans les cristaux clairs. La composition finale moyenne, sur tous les
cristaux, est Na14±3Co32,3±0,8Si100±7Al92±2O378×nH2O avec un rapport Si/Al de ~1,08. A cause
de la migration des cations sous le faisceau d’électrons ou l’échange cationique avec des
protons, dans la composition finale, la charge de la charpente n’apparaît pas totalement
compensée (~14 charges par maille manquantes).
Les images de contraste chimique (Figure 2.7) montrent que la distribution des
éléments chimiques dans les cristaux échangés en colonne est nettement plus inhomogène que
dans ceux échangés en bain (Figures 2.6).
37
Figure 2.7: Images de contraste chimique des cristaux Na,Co-X échangés en colonne.
Le Tableau 2.6 montre qu'aux bords des cristaux le rapport Si/Al est plus faible qu’au centre,
symptôme d'une désilicalisation des cristaux au cours de l’échange.
Tableau 2.6: Distribution des rapports Si/Al dans des cristaux Na,Co-X échangés en colonne
Cet effet de désilicalisation des cristaux a été également remarqué par Le Van Mao et al.
(1994) lors du traitement par une solution aqueuse de carbonate de sodium des zéolithes X, Y
et par Dessau et al. (1992) pour la ZSM-5. Il explique probablement l'amorphisation observée
en DRX.
38
(∅pore = 90 µm). Une rampe de chauffage suffisamment lente a été programmée afin de
résoudre des signaux enthalpiques proches, tout en donnant un signal exploitable (Figure
2.8) :
a) 5 b) 2
0
1
-10
-1
-15
-2
-20
-3
-25
-4
0 100 200 300 400 500 600 700 100 200 300 400 500 600
Température, °C Température, °C
Figure 2.8: Courbes de première déshydratation; a) Analyse thermogravimétrique (ATG), b) Analyse thermique
différentiel (ATD).
Les courbes (Figure 2.8) mettent en évidence une perte de masse totale d’environ 24%
(minitiale = 26,5 mg, mfinale = 20,2 mg). De ~50°C jusqu'à 90°C, la perte de masse (~3,5% de la
masse totale) peut s'interpréter par le départ des molécules d'eau de surface. Puis, entre
environ 90 et 160°C la perte de masse rapide (~16,5% de masse totale), en particulier dans la
zone ~115-150°C, peut être attribuée à la désorption des molécules d'eau ou de TEA. Enfin,
entre ~160 et 220°C (perte de masse de ~4,5%) et entre 220 et ~420°C (~3% de perte de
masse) on peut penser que la perte de masser correspond au départ des molécules d'eau les
plus liées.
La perte de masse dans Na,Co-X est régulière et très similaire aux courbes de perte de
masse de zéolithes X synthétisées sans TEA (Na-X; Na,K-X; Na,Rb-X (Szostak, 1992; Joshi
et al., 2002). Si on attribue l'intégralité de la perte de masse à la déshydratation, on peut
estimer la quantité d'eau adsorbée à ~237 molécules d'eau par maille, soit une composition
Na14±3Co32,3±0,8Si100±7Al92±2O378×237H2O.
Il faut noter que le nombre de molécules d’eau déterminé dans nos analyses est propre
39
à la zéolithe Na,Co-X de rapport Si/Al~1,08 étudiée ici, car le caractère hydrophile dépend du
type et de la quantité de cations compensateurs de charge. Le nombre de molécules d’eau
adsorbées sous conditions ambiantes et le comportement en désorption dépendent des
interactions électrostatiques eau-cation et des liaisons hydrogène avec les oxygènes de la
charpente et les autres molécules d’eau. D’après Joshi et al. (2002), toutes ces interactions
sont «gouvernées» par le cation compensateur de charge. Dans le cas des zéolithes échangées
avec des cations alcalins, par exemple, la vitesse de perte d’eau, la température de
déshydratation complète et le nombre de molécules d’eau dans la structure diminuent avec
l’augmentation de la taille du cation, ce qui s'explique par des interactions moins fortes entre
les cations et les molécules d’eau. Sur cette base (sans prendre en compte l'hydrolyse de la
charpente qui affecte la microporosité de la zéolithe et ses propriétés de diffusion), la zéolithe
Na,Co-X serait plus hydrophile que Na-X.
II.2.6 Déshydratation.
Le programme de déshydratation est basé sur les courbes d'ATG (Figure 2.8) et
présenté sur la Figure 2.9 :
400
350
300
Température, °C
250
200
150
100
50
0
0 20 40 60 80 100 120 140 160 180
Temps, heure
40
A la fin de la déshydratation, les capillaires sont scellés à température ambiante sous
vide secondaire. La couleur des cristaux de Na,Co-X change de rose (état hydraté) à bleu très
foncé (état déshydraté, Figure 2.10).
Figure 2.10 : Cristal de Na,Co-X (~180 µm selon [100]) déshydratée dans son capillaire
Les complexes (Co(bpy)3) étant trop grands pour être introduits dans les cavités
zéolithiques (~13 Å), les métaux de transition sont alors complexés in situ (méthode « ship in
bottle »). Nous avons tenté la formation de composite par deux méthodes: en phase vapeur et
en solution.
II.2.8 Conclusion
41
lente. Cette vitesse de croissance diminue avec la température de mûrissement et avec
l’augmentation du taux de TEA dans le gel. Les synthèses produisent en général un mélange
de trois phases, X, A et P, mais en modifiant la quantité de TEA ou le temps de mûrissement,
on peut exclure une phase ou l'autre du produit final.
Dans les FAUjasites, les échanges cationiques de Na+ par des métaux de transition de
la première série (Cr3+, Fe2+, Co2+, Ni2+, Cu2+, Zn2+) sont souvent partiels et compliqués
(Sherry, 1971; Maes & Cremers, 1975; Firor & Seff, 1978; Heo et al., 1986). Si,
dernièrement, des échanges cationiques quasiment complets avec Zn2+, Ni2+ (Bae & Seff,
2000 deux articles distinctes) ont été revendiqués, les échanges 2Na+ → Co2+ réalisés par Bae
& Seff (1999) ne sont que partiels. L'échange cationique en colonne mène à l'amorphisation
des cristaux de zéolithes X et A malgré la préservation de leurs morphologies; la
microanalyse mettant, en parallèle, en évidence la désilicalisation des cristaux, même en
volume. Dans cette thèse, les cristaux de Na,Co-X pour les études DRX ont été finalement
obtenus par échange cationique en bain. L'échange cationique est incomplet car ~15% de
sodium natifs demeurent dans la structure. L'eau représente ~24% de la masse de la zéolithe
Na,Co-X hydratée, d'après les analyses thermogravimétriques
(Na14±3Co32,3±0,8Si100±7Al92±2O378×237 H2O). Pour les études structurales précises (DRX,
XAS) les cristaux et la poudre de Na,Co-X ont été deshydratés en température sous vide
secondaire.
42
II.3 Mordénite.
43
mordénite « small pore » (Mortier, 1975 ; Shiokawa et al., 1989 ; Campbell & Cheetham,
2002 ; Simoncic & Armbruster, 2004,) mais on constate également des défauts altérant la
connectivité des canaux et la taille des pores des mordénites « large pore ». Mortier (1975),
Simoncic & Armbruster (2004) ont trouvé des défauts de croissance où environ 3 % de la
structure sont décalés de c/2 par rapport au reste. Sanders (1985), Shiokawa et al. (1989),
Campbell & Cheetham (2002) ont lié l'apparition de tels défauts à la présence d'aluminium
dans la charpente, l'importance des défauts augmentant avec la quantité d'Al (cf. études
expérimentales de Rudolf & Garcés, 1994). D'après Campbell et al. (2004) les défauts
structuraux liés au déplacement des chaînes 4MR de c/2 sont corrélés à des déplacements de
ces défauts selon [110] pour maintenir la structure. La présence de dislocations de vecteur de
Burgers ½[110] est d'ailleurs mise en évidence par Microscopie Electronique à Haute
Résolution (HRTEM) sur des mordénites commerciales (Song, 1999).
Les propriétés d’adsorption dépendent des conditions de synthèse (Sano et al., 2001 ;
Oumi et al., 2003) et des traitements post-synthèses (Lozano-Castello et al., 2006). En
particulier, les cristaux de Na-mordénite synthétisés en présence de BuOH n'adsorbent
quasiment pas l'azote, contrairement à ceux échangés avec NH4+ (Sano et al., 2001 ; Oumi et
al., 2003). Les cristaux synthétisés en présence d'autres alcools aliphatiques adsorbent plus
facilement l’azote que ceux synthétisés avec du BuOH. Les cristaux synthétisés sans alcool
ont les meilleures capacités d’adsorption (azote, butane et isobutane) (Sano et al., 2001 ;
Lozano-Castello et al., 2006).
Les cristaux de mordénite de grande taille (lmax= 42 µm) sont synthétisés selon le
protocole de Sano et al. (2001) en présence de 1-BuOH. Ce protocole a été optimisé au
LCM3B. par V. Coulibaly et A. Bouché qui ont obtenu des cristaux nettement plus grands
(140 µm × 100 µm × 30 µm) en diminuant la température de cristallisation de 160°C à 150°C,
et en introduisant la silice pulvérulente dans le gel de synthèse. Ce protocole est utilisé au
cours de ce travail (Protocole 1). Afin d'augmenter encore la taille des cristaux, nous avons
réalisé des synthèses avec un mûrissement séparé des solutions silicatée et aluminatée
(Protocole 2). Pour tester les propriétés d'adsorption, des synthèses sans alcool ont été aussi
44
préparées (Protocole 3).
II.3.4.1 Protocole 1
Dans 3,33 mol d'eau distillée préchauffée dans un bain-marie à T~60°C on ajoute sous
agitation magnétique 0,058 mol de soude (NaOH, Carlo Erba, 99 %), 0,002 mol de nitrate
d'aluminium (Al(NO3)3×9H2O, Riedel de Haën), 0,067 mol de silice pulvérulente (Merck,
99%) et 0,017 mol de silice colloïdale (Riedel de Haën, 99 %). On laisse le gel mûrir sous
agitation magnétique à Tm~60°C pendant environ une heure avant l’ajout de 0,125 mol de 1-
BuOH (Fluka). L’ensemble est transféré dans une « chemise » en téflon (volume de 130 ml),
placée elle-même dans une bombe en acier (Figure 2.11) et est déposée dans un four à
Tsynthèse=150°C. Après environ 21 jours, la totalité du gel est consommée et il ne reste plus
qu'une solution transparente et un agglomérat cristallin incolore.
Figure 2.11: Chemise en téflon et bombe en acier utilisées pour la synthèse des mordénites
Les cristaux sont rincés sur un dispositif d'ultrafiltration (∅pores ~ 0,2 µm) avec ~500
ml d'eau distillée puis laissés à sécher en étuve à 70°C. Ce protocole mène à la formation de
cristaux de la phase mordénite pure, de morphologie prismatique (Figure 2.12) et de deux
tailles (en moyenne les cristaux sont de ~140 µm × 100 µm × 40 µm et ~ 80 µm × 60 µm × 30
µm).
Diverses durées de mûrissement du gel, compris entre 0 et 14 h ont été testés. La durée
optimale est de 2-3 h (cristaux de taille maximale 200 µm × 150µm × 40µm). Un
45
mûrissement de 14 h diminue significativement la taille des cristaux (lmax ~ 80 µm).
II.3.4.2 Protocole 2
D'après Hincapie et al. (2004), les cristaux de mordénite sont plus grands lorsque la
silice est bien dissoute avant d’être mélangée avec la solution d'aluminate de sodium. Nous
avons donc modifié le protocole 1 en séparant la préparation des deux solutions de silicate et
d’aluminate.
La solution aqueuse est préparée en dissolvant 0,013 mol de soude (Carlo Erba, 99 %)
dans 0,56 mol d’eau distillée où on ajoute 0,002 mol de nitrate d'aluminium (Al(NO3)3×9H2O,
Riedel de Haën) dans un bain-marie (~60°C) et sous agitation magnétique.
46
Utilisation des différentes sources d’aluminium (Influence du contre-ion)
Le protocole précédent est utilisé pour tester l'influence du contre-ion apporté par la
source d'aluminium sur le produit final (Tableau 2.7).
Les cristaux les plus grands sont obtenus avec le contre-ion SO42- mais les cristaux
sont globalement de mauvaise qualité (cristaux cassés). Ces résultats sont en accord avec ceux
de Shaikh et al. (1993) qui ont également noté que les anions sulfates favorisent la formation
de grands cristaux prismatiques.
II.3.4.3 Protocole 3
Afin de tester les propriétés d'adsorption des cristaux, nous avons modifié le protocole
1 en supprimant l’alcool. La durée typique de synthèse est conservée (~20 jours) mais la
morphologie des cristaux est différente (prismes allongés selon b et fins selon l'axe c, Figure
2.13) et se rapproche de celle reportée par Simoncic & Armbruster (2004) et Lozano-Castello
et al. (2006) ou encore de celle des cristaux commerciaux, réputés avoir une meilleure
capacité d'adsorption que les cristaux synthétisés en présence d’alcool (Sano et al., 2001 ;
Oumi et al., 2003). La taille des cristaux (~100 µm × 70 µm × 30 µm) est homogène dans la
synthèse : elle est la même lorsque SO42- et NO3- sont utilisés comme contre-ions, et est
globalement plus faible que celle obtenue avec les autres protocoles.
Figure 2.13: Image au microscope optique d’un cristal de mordénite synthétisé sans alcool (la longeur maximale
est ~100 µm selon [010])
47
II.3.4.4 Synthèses à partir de verre (Gilbert & Mosset)
Nous avons tenté une série des synthèses selon cette technique, en utilisant des verres
commerciaux (verre sodocalcique, quartz). La taille maximale des cristaux obtenus n’étant
que de ~40 µm, cette voie de synthèse a été abandonnée.
Nous avons réalisé un total d’environ 50 synthèses de mordénite. Ces essais ont mis en
évidence quelques paramètres importants pour la croissance de monocristaux de grande taille,
à partir du protocole développé par Sano et al. (2001) et modifié auparavant par V. Coulibaly
et A. Bouché (LCM3B):
48
II.3.5 Caractérisation
Les analyses élémentaires ont été réalisées au service commun de microanalyse (UHP
Nancy-I) dans des conditions identiques à celles utilisées pour les zéolites X. La composition
chimique des différents cristaux est une moyenne des mesures effectuées sur 3 points par
cristal (Tableau 2.8).
Dans les cristaux de type 1, la quantité de sodium mesurée est inférieure à la quantité
d’aluminium, contrairement aux cristaux de type 2, 3. Les incertitudes de mesure sont de
l'ordre de quelques unités (comme dans le cas de Na,Co-X). Si on assume une compensation
intégrale de la charge par des ions Na+ (n(Na+) = n(Al3+)) dans la composition finale, on
obtient les compositions suivantes:
Les mordénites sont réputées être des zéolithes hydrophiles. Un échantillon de 30,7
mg de mordénite de type 1 a été caractérisé par analyse thermogravimétrique (ATG) / analyse
thermique différentielle (ATD) sur l’appareil 92-16.18 SETARAM par le Pr. M. François
49
(LCSM, UMR 7555, Université Nancy-I) afin de préciser la quantité d'eau absorbée. Deux
mesures ont été réalisées à 2 mois et demi d'intervalle sur le même échantillon afin d'évaluer
la reprise en eau à l'air libre.
Première analyse
Les courbes de perte de masse (TG) et de flux de chaleur (HF) (ATD) sont montrées
sur la Figure 2.14.
Figure 2.14: Première courbe de perte de masse (en bleu) et du flux de chaleur (en rouge) ; mordénites de type 1.
Seconde analyse
50
l’analyse, l’échantillon pèse 29,3 mg, ce qui correspond à une perte de masse totale de 12,9 %
(Figure 2.15).
Figure 2.15 : Seconde courbe de perte de masse (en bleu) et du flux de chaleur (en rouge) ; mordénites de type 1.
Discussion
Au cours de l'analyse n°1 (Figure 2.14) l'échantillon perd 12,2 % de sa masse totale
(perte rapide de ~7% de masse à 80-180°C , perte progressive de ~3,5% de masse entre ~180-
300°C). La courbe ATD montre que jusqu'à ~300°C les réactions sont endothermiques et les
signaux ATG s'interprètent comme une déshydratation de l'échantillon. Entre ~300 et 430°C
le pic exothermique (ATD) est attribué à la décomposition de BuOH (par comparaison avec
les analyses TG/DSC de mordénites naturelles et synthétisées sans alcool (Szostak, 1992) qui
ne mettent pas en évidence ce pic) et s'accompagne d'une perte de masse lente de ~1,3 % de
sa masse. Le chauffage jusqu'à 550°C fait perdre à l'échantillon encore ~0,4 % de sa masse
totale. Au delà de 550°C, la masse de l'échantillon demeure stable.
Dans l'analyse n°2 (Figure 2.15), l'évolution de perte de masse pour les réactions
endothermiques est sensiblement la même que celle d'analyse n°1, mais les pertes de masse se
produisent globalement à plus basse température. Dans la plage de température 340-550°C on
observe de nouveau des réactions exothermiques qui s'accompagnent d'une perte de masse de
~0,7 %. Si on compare les 2 analyses thermiques, on retrouve le même comportement à basse
température (desorption d'eau surfacique puis structurale), mais les réactions exothermiques à
plus haute température sont probablement de nature différente dans les deux cas (Figures 2.14
51
et 2.15).
Pour la calcination à l'air, les cristaux de mordénite sont simplement déposés au fond
d'une coupelle en céramique qui est placée dans un four programmable (Figure 2.16).
Après calcination, les échantillons sont déshydratés avec le dispositif utilisé pour la
préparation des cristaux de Na, Co-X. Les programmes de calcination et de déshydratation
sont établis à partir des analyses thermogravimétriques (Tableau 2.10). Après déshydratation
les tubes/capillaires contenant les cristaux sont scellés sous vide secondaire.
Tableau 2.10: Programmes de calcination sous l’air, déshydratation (~10-6 mbar) des cristaux de mordénite.
52
II.3.7 Préparation de composites.
L'adsorption des molécules organiques dans les zéolithes peut se faire par échange
cationique en phase liquide (Ramamurthy et al., 1993 ; Simoncic & Armbruster, 2005),
dopage en phase solide (Hoppe et al., 1994) ou dopage en phase vapeur (Reck et al., 1996).
Dans notre cas, l'efficacité de l'échange cationique est limitée par la faible quantité de cations
dans les mordénites; l'efficacité du dopage en phase solide est limitée par les fortes
interactions intermoléculaires et aboutit à un faible taux de dopage. Nous avons donc réalisé
les dopages en phase vapeur à l'aide d'une ampoule de Pyrex composée de deux
compartiments (Figure 2.17).
Le «dopant» organique est introduit dans un des compartiments et les cristaux de mordénite
préalablement calcinés et déshydratés sont placés dans l'autre. Afin d'éviter la réhydratation
des cristaux, l'ampoule est raccordée à une pompe à vide et le compartiment avec le produit
organique est maintenu au froid afin de limiter sa sublimation et le compartiment contenant
les cristaux est chauffé à ~200°C afin de désorber l'eau de surface éventuellement piégée lors
du transfert des cristaux dans l'ampoule. Après pompage pendant ~24 h sous ~5×10-5 mbar,
l'ampoule est scellée sous vide et placée dans un four à gradient de température.
Choix de la molécule
53
II.3.7.1 Pérylène – mordénite.
Le pérylène (C20H12, M=252,31 g/mol, point de fusion 276°C) est une molécule
fluorescente de couleur jaune, stable dans les conditions ambiantes (Figure 2.18). Ses
propriétés de fluorescence dépendent fortement de son environnement (cristal, solution) :
interactions pérylène-pérylène, pérylène-solvant et interactions avec un substrat (Wakayama
et al., 2006). Nous nous sommes intéressés à ce titre au confinement de cette molécule dans la
mordénite. Les dimensions du pérylène (longueur ~10,2 Å, largeur ~7,6-8 Å) et sa symétrie
(mmm) sont particulièrement adaptées aux dimentions et la symétrie des canaux droits de la
mordénite.
Figure 2.19: Composite (pérylene@mordénite) (type 1, la longueur des cristaux est ~140-170 µm selon [001])
Les cristaux de mordénite de type 4 (la longeur des cristaux est ~100µm selon [010])
se colorent rapidement dès le pompage. Le dopage a duré 2 mois à 180°C, sans modification
perceptible de la coloration durant les dernières 6 semaines de dopage.
54
II.3.7.2 Thio-indigo – mordénite.
Le dopage des mordénites (type 4) a duré 16 jours à une température de 180°C. Nous
avons constaté une coloration rapide des cristaux pendant les 2-3 premiers jours, puis plus
aucun changement de couleur. Au final les cristaux obtenus sont violets (Figure 2.21).
Figure 2.21: Cristaux du composite (thio-indigo@mordénite) (type 4, la longeur maximale des cristaux est ~100
µm selon [010]).
55
II.3.7.3 dmpNA-mordénite.
Lors du pompage, les cristaux de mordénite (type 1 et 2) n’ont pas changé de couleur,
mais sont devenus très rapidement jaunes-verdâtres dès que l’ampoule a été placée dans un
four à 130°C, la couleur s’intensifiant pendant 5 jours. Les 10 jours suivants, nous n’avons
pas constaté de changement de couleur.
II.3.8 Conclusion
56
Chapitre III
57
58
III Etats d’hydratation des zéolithes Na,Co-X................................................................. 61
III.1 Description de la structure de la zéolithe X (FAUjasite). ................................... 61
III.2 Problématique de l'étude. ...................................................................................... 62
III.3 Diffraction des rayons-X........................................................................................ 65
III.3.1 Echantillons étudiés.......................................................................................... 65
III.3.2 Conditions expérimentales. .............................................................................. 65
III.3.3 Structure du cristal totalement déshydraté. ...................................................... 66
III.3.4 Structures des cristaux (partiellement) hydratés. ............................................. 70
III.3.4.1 Structures des cristaux partiellement hydratés 2 et 3 ................................... 70
III.3.4.2 Structures des cristaux hydraté (4) et réhydraté (5). .................................... 75
III.3.5 Distribution des cations .................................................................................... 79
III.3.6 Effet de la déshydratation sur Na,Co-X. .......................................................... 86
III.3.6.1 Migration des cations. .................................................................................. 86
III.3.6.2 Distorsion du squelette lors de la déshydratation......................................... 90
III.3.7 Structures déshydratées des zéolithes Mn+-X................................................... 92
III.3.8 Distribution cationique dans des différentes structures échangés avec des
métaux de transitions........................................................................................................ 94
III.4 Spectroscopie d’Absorption des rayons X (EXAFS)........................................... 97
III.4.1 Interprétation qualitative des spectres XANES................................................ 97
III.4.2 Spectre EXAFS des cristaux de Na,Co-X déshydratés. ................................... 99
III.4.2.1 Spectre expérimental .................................................................................... 99
III.4.2.2 Spectre modèle ........................................................................................... 101
III.4.2.3 Stratégie d’affinement. ............................................................................... 103
III.4.3 Spectre EXAFS de poudre Na,Co-LSX déshydratée. .................................... 112
III.4.4 Spectre EXAFS de poudre Na,Co-LSX réhydratée. ...................................... 115
III.4.5 Bilan des analyses EXAFS............................................................................. 119
III.5 Simulations Monte-Carlo de la distribution cationique. .................................. 121
III.5.1 Objectifs, originalité des calculs. ................................................................... 121
III.5.2 Champs de forces : interactions entre les espèces.......................................... 122
III.5.3 Hypothèses de départ des simulations............................................................ 124
III.5.4 Géométrie de la charpente et stabilité des sites.............................................. 124
III.5.5 Distribution des cations dans la structure déshydratée................................... 125
III.5.6 Etude de la distribution des cations et H2O dans la structure hydratée.......... 127
III.5.7 «Dynamique» de la redistribution cationique lors de l'hydratation................ 131
III.6 CONCLUSIONS................................................................................................... 134
III.7 Composite Na,Co(bpy)3-X....................................................................................... 136
59
60
III Etats d’hydratation des zéolithes Na,Co-X
III.1 Description de la structure de la zéolithe X
(FAUjasite).
Idéalement, le groupe d’espace du type zéolithique "FAU", où tous les tétraèdres sont égaux
et réguliers, est Fd 3 m . Dans les structures réelles, le rapport Si/Al contrôle l’enchaînement
des tétraèdres TO4 et, en conséquence, le groupe d’espace. En particulier, lorsque Si/Al=1,
l'alternance de tétraèdres SiO4 et AlO4 dans la structure abaisse la symétrie de Fd 3 m à Fd 3
(règle de Löwenstein, 1954, interdisant l'interconnexion de deux tétraèdres AlO4). Le taux
d'aluminium dans la charpente conditionne également le nombre de cations compensateurs de
charge dans les cavités zéolithiques et les applications industrielles des zéolithes de type
"FAU" : les zéolithes ayant Si/Al<1,5 ("X"), sont utilisées pour la séparation sélective de gaz,
tandis que les zéolithes de ratio Si/Al > 1,5 ("Y") sont plutôt utilisées en catalyse.
61
Conventionnellement, dans les zéolithes de type "FAU", la distribution des cations est
décrite par différents sites identifiés selon la nomenclature suivante (Figure 3.1) :
Site I' : sur l'axe 3, à proximité d'une face 6MR du prisme hexagonal.
Sites II, II' : à proximité du centre d'une face 6MR d'une cage sodalite, sur l’axe 3
Site III : dans la supercage, près d'un anneau 4MR de la cage sodalite, sur un axe 2.
Les zéolithes X (et Y) échangées avec des métaux de transition sont largement
étudiées pour leurs propriétés en catalyse hétérogène (réactions d’oxydation, Blower & Smith,
1993; Silva et al., 2006; Li et al., 2008; Parpot et al., 2008) et d’adsorption sélective des
polluants organiques (Díaz et al., 2004). L'utilisation potentielle comme élément de base d’un
dispositif électronique (mémoire) est évoquée par Mizuno & Lunsford (Co2+(bpy)3-zéolithe
Y), 1983). C'est à cette application potentielle que s'est initialement intéressée cette thèse. Au
début des années 80, l’équipe de Lunsford (DeWilde et al., 1980; Quayle et al., 1982; Mizuno
& Lunsford, 1983) a étudié la complexation de métaux de transition Ru2+, Fe2+, Co2+ par 3
ligands bipyridine après encapsulation dans les supercages de zéolithes Y (Figure 3.2). Le
composite Co2+(bpy)3-zéolithe Y s’était avéré particulièrement intéressant car il possédait des
propriétés magnétiques différentes de celles observées à l'état libre, en particulier une
transition du spin ( 6 t 2 g 1 e g ←⎯→
T 5
t 2 g 2 e g ) induite thermiquement (voir également l'étude plus
récente de Tiwary & Vasudevan, 1998). Le complexe Co2+(bpy)3 dans (ClO4)22- (en forme de
62
sel) est dans l’état fondamental haut spin stable pour toute température, tandis que le
complexe encapsulé dans la zéolithe Y est dans l’état bas spin à 77 K et bascule vers un état
haut spin à plus haute température (ambiante dans l’étude de Mizuno & Lunsford, 1983 et
supérieure à 525 K dans celle de Tiwary & Vasudevan, 1998. A des températures
intermédiaires, on a un mélange d'états magnétiques (Mizuno & Lunsford, 1983; Tiwary &
Vasudevan, 1998).
Figure 3.2 : Schéma d’un complexe Co2+(bpy)3 dans la supercage d’une zéolithe X ou Y.
Dans tous les cas, le champ électrostatique dans les cavités modifie la distribution
électronique du complexe (Quayle et al., 1982; Sieber et al., 2000,). Quayle et al. (1982),
63
dans le cas de Fe2+(bpy)3-Y, suggèrent que la distribution électronique des ligands
polarisables (bipyridines) est influencée par la charge négative de la zéolithe qui modifie la
force du champ du ligand.
Dans ce contexte nous nous sommes intéressés au composite (Co2+(bpy)3-X) basé sur
la zéolithe X comme matrice hôte. La zéolithe X étant plus alumineuse que Y, sa charpente
porte une charge anionique plus importante compensée par un plus grand nombre de cations
compensateurs, d'où l'espoir d'une exacerbation des propriétés magnétiques du nouveau
composite. Nous nous sommes intéressés, dans une phase préliminaire, à la structure de la
zéolithe Na,Co-X, puisque malgré ses applications actuelles/potentielles, la structure
cristallographique précise de la zéolithe bicationique parente reste incertaine.
64
III.3 Diffraction des rayons-X
Cinq cristaux de Na,Co-X différents par leur état d’hydratation sont préparés selon les
protocoles décrits au Chapitre II :
Les expériences de DRX sont réalisées dans des conditions similaires (radiation
Mo(Kα), température ambiante, diffractomètre XCalibur II) pour tous les échantillons, de
façon à disposer a priori d'une même statistique de mesure (Tableau 3.1). Le groupe d’espace
Fd 3 est vérifié en se basant sur ses extinctions systématiques et la distribution statistique des
[110]
intensités des réflexions hkl / hlk (absence du miroir m de Fd 3 m ). Les intensités des
réflexions de Bragg sont intégrées avec le programme CrysAlisPro (version 1.171.29.2),
corrigées de l'absorption et moyennées avec Sortav (Blessing, 1995, 1997); les affinements
structuraux sont faits avec le programme Jana2000 (Petřίček et Dušek, 2005) (ANNEXE A).
65
Tableau 3.1: Conditions expérimentales
wRint (%) 4,54 / 22,08 6,75 / 15 4,53 / 8,37 8,04 / 16,61 7,10 / 9,63
R (%) 6,6 / 15,09 7,68 / 16,78 8,61 / 12,89 12,35 / 20,39 12,89 / 20,72/
Rw (%) 9,88 / 11,04 12,19 / 12,75 15,38 / 16,76 27,33 / 32,64 27,25/ 25,92/
Gof 4,36 / 3,52 4,82 / 3.55 6,39 / 5,28 12,50 / 8,59 18,02 / 11,64
1
w =
Avec
( ( )) (
2
σ Fobs
2
+ 0, 005 × F 2
obs
)
2 et affinement en F2.
66
conjoint des facteurs d'occupation de site et de Debye-Waller est stable (occ=16(1)%,
(Uiso=0,11 Å2), et permet d’éliminer effectivement le résidu sur la carte de Fourier différence
(Figure 3.3 b, R=6,91 % et Rw=10,94 %).
a) b)
Site III Na3
Figure 3.3 : Cartes Fourier différence autour du site III : avant (a) et après (b) la modélisation de Na3. Contour
0,1 e-/Å3
A ce stade, le seul pic sortant du bruit de fond (ρ=1,28 e-/Å3) est sur le site I' (Figure
3.4 a). La modélisation de ce pic par un cation est discutable :
d'une part ce pic se trouve sur un site de haute symétrie (position Wyckoff: 32e) où
l'accumulation d'erreurs est importante (Cruickshank & Rollett, 1953),
d'autre part, à cause de la répulsion électrostatique entre les cations, les sites I et I' ne peuvent
pas être occupés simultanément (exception est K+-X, Zhu & Seff, 2000). Le site I voisin est
occupé à 94,2 % par des cations Co, ce qui exclut la possibilité d'occupation du site I' à plus
de 5,8%.
De surcroît à ce stade la composition cristallographique est déjà en accord avec celle
obtenue par microanalyse. Par conséquent, nous avons testé deux stratégies d'affinement avec
des cations Na+ ou Co2+ sur le site I'. L’occupation du site I' par des cations sodium aboutit à
un taux d'occupation de 12(1) % (Uiso =0,06 Å2) incompatible avec la présence de Co1, et
finalement, le résidu localisé sur le site I' est modélisé comme un cation cobalt avec une
occupation de 4,5(3) % (Uiso=0,05(1) Å2). Dans les deux cas, les facteurs d'accord sont
similaires (R ~ 6,60 % et Rw ~9,88%) et la carte Fourier-différence montre que la densité
électronique est correctement modélisée (Figure 3.4 b) en écartant le doute de la présence d'un
cation sur le site I'.
67
a) b) Co1'
Co1 Co1
Site I'
Figure 3.4: Cartes de Fourier-différence avant (a) et après (b) la modélisation de la densité résiduelle sur le site I'
par un cation Co2+. Contour de 0,1 e-/ Å3
Les 2 hypothèses d'occupation de ce site par Na+/Co2+ ont également été testées par
Spectroscopie d'Absorption des rayons X au seuil du cobalt; les résultats correspondants sont
discutés au Chapitre III.4.
Après l'introduction du cation en site I', plus aucun pic de densité électronique ne sort
du bruit de fond, et la structure peut être considérée comme complète. La combinaison des
analyses chimiques par microsonde (dosage de Si et Al), de l’étude cristallographique (Co, Na
et O) et de l'électroneutralité du cristal indique une composition approximative
Na16Co38Al92Si100O384 pour le cristal étudié. La densité résiduelle s'accumule principalement
autour des oxygènes de la charpente (Figure 3.5) dans les liaisons T-O (T = Si4+, Al3+),
traduisant les interactions inter-atomiques, mais la résolution des données est trop faible
(Smax= 0,63 Å-1) pour entreprendre un affinement multipolaire.
O2
Co2
O4
Figure 3.5: Densité résiduelle dans un plan perpendiculaire à l'axe 3 et passant par le site II.
Contour de 0,1 e-/Å3.
Les détails de l'affinement et les paramètres structuraux sont rassemblés dans les
Tableaux 3.2 et 3.3.
68
Tableau 3.2 : Cristal 1, affinements et statistiques
Description des Charpente fixe, Charpente, Charpente Charpente Charpente Charpente Charpente Uij, Co1 Charpente
Réflexions
affinements Co1, +Co2 Co1, Co2 + Uij, Uij, Co1 Uiso, Uij, Co1 +Uij, Uij, Co1 Uij, Uij, Co2 Uij, Na2, + Uij, Co1 Uij, Co2
k, occ; Uiso Uiso Co1, Co2 Co2 + Uij. Co2 Uij. Co2 Uij,+ Na2 Na3 Uij, Na2, Na3, +
Uiso Na1' ou Co1'
ρmax (e-/Å3) Si : 9,95 Site II : Site II : Site II : 2,7 Site II : 2,62 Site III : 1,55 Site I': 1,28 O4 : 0,66
Al : 8,78 5,82 5,47
Observées
Atome x y z Wyckoff Occ., % U11 ou Uiso U22 U33 U12 U13 U23
Al -0,05396(5) 0,12073(5) 0,03795(5) 96 100 0,0304(7) 0,0212(7) 0,0251(7) -0,0002(6) -0,0037(5) -0,0043(6)
Si -0,05200(5) 0,03333(5) 0,12244(5) 96 100 0,0297(7) 0,0275(7) 0,0215(6) -0,0038(5) 0,0002(5) -0,0028(6)
O1 -0,1046(1) -0,0028(1) 0,1113(1) 96 100 0,045(2) 0,067(2) 0,044(2) -0,027(2) 0,012(2) -0,024(2)
O2 -0,0017(1) -0,0019(1) 0,1489(1) 96 100 0,050(2) 0,050(2) 0,025(2) 0,011(2) 0,004(2) 0,007(2)
O3 -0,0302(1) 0,0569(1) 0,0625(1) 96 100 0,038(2) 0,028(2) 0,027(2) 0,006(1) 0,003(1) -0,002(1)
O4 -0,0652(1) 0,0814(1) 0,1641(1) 96 100 0,066(2) 0,035(2) 0,037(2) -0,011(2) 0,010(2) -0,011(2)
Co1 0,0000 0,0000 0,0000 16 94,2(6) 0,0276(5) 0,0037(4)
Co2 0,2182(1) 0,2182(1) 0,2182(1) 32 65,7(9) 0,049(1) 0,011(1)
Na2 0,2380(5) 0,2380(5) 0,2380(5) 32 25(2) 0,029(5)
Na3 -0,147(1) 0,1250(0) 0,1250(0) 48 16(1) 0,11(1)
Co1′ 0,0608(7) 0,0608(7) 0,0608(7) 32 4,5(3) 0,06(1)
69
III.3.4 Structures des cristaux (partiellement) hydratés.
Dans les structures plus ou moins partiellement hydratées, les espèces extracharpente
(cations et molécules d'eau) ne sont qu'en partie localisées, et les études divergent quant à la
nature des espèces occupant les sites. Pour simplifier les comparaisons entre les structures,
des dénominations homogènes des sites occupés par des molécules d’eau ont été adoptées :
Cristal 2
Comme lors de l'étude du cristal 1, les premiers cations sont localisés successivement
sur le site II (ρ=27,51 e-/ Å3, Co2), puis sur le site I (ρ=19,02 e-/ Å3, Co1). En revanche, à ce
stade, le résidu principal est sur le site I' (ρ=11,31 e-/ Å3) et est modélisé comme
précédemment par un cation cobalt (Co1') (après l’affinement Co1 : occ.=38,4(6)%,
Uiso=0,037(3) Å2 et Co1' : occ.=31,8(6) %, Uiso=0,047(3) Å2 ; R=14,87 % et Rw=28,42 %). Le
pic restant (ρ=6,70 e-/ Å3) situé dans la cage sodalite sur le site II (ou II') à 1,24 Å de Co2,
s'affine de façon satisfaisante comme un sodium (Na2 : occ.=31(1)% Uiso=0,04(4) Å2,
R=11,70 % et Rw=20,69 %, Figure 3.6). L'unique molécule d'eau localisée dans la structure
apparaît comme un pic (ρ=1,79 e-/ Å3) dans la supercage à 1,84 Å de Co2 sur les cartes
Fourier-différence. Son introduction dans la structure et son affinement (Ow1, occ.=36(3) %,
Uiso=0,11(2) Å2, Figure 3.6), améliorent les facteurs d’accord (R=11,13 % et Rw=18,56 %).
70
Les affinements montrent que pour conserver des facteurs de Debye-Waller réalistes,
l’occupation de la molécule d’eau doit être inférieure à celle du cobalt Co2 et seulement 2/3
des cobalts apparaissent hydratés.
Le dernier atome introduit dans la structure (ρ=1 e-/ Å3) est un cation cobalt situé sur
un site général de la supercage (Co3, occ.=5,5(5) % et Uiso=0,11(2) Å2).
Figure 3.6 : Vue ORTEP (visualisation des ellipsoïdes d’agitation thermique des atomes): Coordination des
cations Co2 et Na2 le long de l’axe 3
71
Tableau 3.4 : Paramètres structuraux de Na,Co-X partiellement déshydraté (cristal 2)
Atome x y z Wyckoff Occ., % U11 ou Uiso U22 U33 U12 U13 U23
Al -0,05272(7) 0,12241(8) 0,03754(7) 96 100 0,030(1) 0,024(1) 0,026(1) -0,001(1) 0,0016(8) -0,003(1)
Si -0,05131(7) 0,03536(7) 0,12391(8) 96 100 0,032(1) 0,028(1) 0,026(1) 0,0022(8) -0,0016(9) -0,002(1)
O1 -0,1019(2) -0,0009(2) 0,1064(2) 96 100 0,068(4) 0,065(4) 0,063(4) -0,026(3) -0,003(3) -0,023(3)
O2 -0,0021(2) -0,0011(2) 0,1490(2) 96 100 0,052(3) 0,049(3) 0,033(3) 0,023(3) 0,004(3) 0,005(3)
O3 -0,0276(2) 0,0653(2) 0,0695(2) 96 100 0,069(4) 0,052(3) 0,049(4) 0,030(3) 0,026(3) 0,028(3)
O4 -0,0692(2) 0,0779(2) 0,1691(2) 96 100 0,045(3) 0,033(3) 0,033(3) -0,001(2) -0,001(2) -0,010(2)
Co1 0,0000 0,0000 0,0000 16 38,4(6) 0,037(3) 0,008(2)
Co1' 0,0608(1) 0,0608(1) 0,0608(1) 32 31,8(6) 0,047(3) 0,007(2)
Co2 0,23051(9) 0,23051(9) 0,23051(9) 32 54,6(6) 0,048(2) 0,016(1)
Na2 0,2002(3) 0,2002(3) 0,2002(3) 32 31(1) 0,04(4)
Co3 -0,184(1) 0,086(1) 0,092(1) 96 5,5(5) 0,11(2)
Ow1 0,2734(7) 0,2734(7) 0,2734(7) 32 36(3) 0,11(2)
72
Cristal 3
Les espèces extracharpente sont localisées dans le même ordre que celles du cristal 2,
également partiellement hydraté. Les cations cobalt sont localisés sur le site II (ρ=33,07 e-/Å3,
Co2) puis en site I (ρ=11,51 e-/ Å3, Co1) et en site I' (ρ=11,87 e-/ Å3, Co1'). Les cartes Fourier
différence révèlent ensuite le cation sodium (Na2) en site II (ρ=7,44 e-/ Å3) puis la molécule
d'eau Ow1 (ρ=2,34 e-/ Å3) située dans la supercage à 1,90(2) Å de Co2. L'affinement des
facteurs d'occupation de Ow1 et Co2 (Co2 : occ.=57,9(6) %, Uiso=0,028(1) Å2 ; Ow1 :
occ.=39(3)%, Uiso=0,10(2) Å2) montre qu'environ deux tiers des Co2 sont hydratés comme
pour le cristal 2.
Les deux seuls résidus de densité électronique significatifs sont sur des sites généraux
dans la supercage : le premier (ρ=1,15 e-/ Å3) se situe à 2,30 Å de l’oxygène O1 (site III'a) et
le second (ρ=1,13 e-/ Å3) se trouve à 2,78 Å de l’oxygène O2 près de la fenêtre 4MR du
prisme hexagonal (site III'b). Etant donné la faible distance (1,8 Å) entre ces deux pics et la
composition chimique cristallographique à ce stade de l'affinement (~29 cobalts et ~13
sodiums), le premier pic sur le site III'a est modélisé comme un cation cobalt (Co3, occ.=
6,1(3)%, Uiso=0,1 (fixe) Å2) et le second comme une molécule d’eau (Ow4). L'affinement
libre des facteurs d'occupation de site et d'agitation thermique de Co3 et Ow4 converge vers
une solution peu réaliste où le nombre de molécules d’eau est supérieur à celui des cations
Co3 auxquels elles sont coordonnées. Si les populations sont contraintes, l'affinement du
facteur de Debye-Waller converge vers une valeur négative pour Ow4. La contrainte à une
valeur "raisonnable" (Uiso=0,1 Å2) mène à une mauvaise modélisation de la densité autour
d’Ow4. La substitution du cation sodium Na3, à la place de Co3 sur le site III'a
(occ.(Na3)=12,9(7) % et occ.(Ow4)=24(1) %), ou la substitution de la molécule d’eau Ow4
par un sodium aboutit à des résultats similaires.
73
Tableau 3.5 : Paramètres structuraux de Na,Co-X partiellement déshydratée (cristal 3)
Atome x y z Wyckoff Occ., % U11 ou Uiso U22 U33 U12 U13 U23
Al -0,05245(7) 0,12339(7) 0,03719(6) 96 100 0,0165(9) 0,0140(9) 0,0131(9) -0,0008(8) 0,0008(7) -0,0027(8)
Si -0,05166(6) 0,03628(6) 0,12455(7) 96 100 0,0201(9) 0,0161(9) 0,0169(9) 0,0007(7) -0,0019(8) -0,0027(7)
O1 -0,1015(2) -0,0001(2) 0,1045(2) 96 100 0,056(4) 0,054(4) 0,057(4) -0,032(3) 0,002(3) -0,024(3)
O2 -0,0013(2) -0,0004(2) 0,1491(2) 96 100 0,045(3) 0,040(3) 0,023(3) 0,025(2) 0,007(2) 0,006(2)
O3 -0,0276(2) 0,0697(2) 0,0725(2) 96 100 0,055(4) 0,037(3) 0,037(3) 0,018(3) 0,016(3) 0,016(2)
O4 -0,0710(2) 0,0754(2) 0,1727(2) 96 100 0,029(3) 0,025(3) 0,029(3) -0,001(2) -0,001(2) -0,012(2)
Co1 0,0000 0,0000 0,0000 16 19,8(6) 0,025(4)
Co1' 0,0623(2) 0,0623(2) 0,0623(2) 32 24,9(6) 0,043(3) 0,014(2)
Co2 0,23001(7) 0,23001(7) 0,23001(7) 32 57,9(6) 0,028(1) 0,0087(8)
Na2 0,2008(3) 0,2008(3) 0,2008(3) 32 39(6) 0,04(4)
Co3 -0,158(1) 0,057(1) 0,151(1) 96 6,1(3) 0,1
Ow1 0,2741(8) 0,2741(8) 0,2741(8) 32 39(3) 0,10(2)
Ow4 -0,174(5) 0,083(5) 0,085(5) 96 6,1(3) 0,1
74
III.3.4.2 Structures des cristaux hydraté (4) et réhydraté (5).
Les affinements structuraux des cristaux hydratés sont instables à cause de la quantité
importante des molécules d’eau désordonnées et non localisées. Les occupations et les
facteurs d’agitation thermique sont fortement corrélés et doivent être maintenus fixes pendant
des affinements. De plus, dans le cas de l’occupation partielle d’un site, la discrimination de
la nature chimique de l'occupant (Na+/Co2+ou H2O) est difficile par DRX et, par conséquent,
les différentes hypothèses conduisent à des facteurs statistiques similaires.
Par exemple, dans la structure hydratée les cartes Fourier différence mettent en
évidence un pic résiduel sur le site I qui peut être modélisé par un cation cobalt (Borissenko et
al. (2008), R=12,57 % et Rw=31,1 6%) ou sodium (R=12,35 % Rw=27,33 %). Cette dernière
hypothèse est en accord avec la distribution des cations déterminée par les simulations Monte-
Carlo (Chapitre III.5) pour les cristaux 4 et 5.
Cristal 4
Les pics successivement identifiés sont Co1' en site I' (ρ=19,47 e-/ Å3), une molécule
d'eau Ow2 (ρ=11,35 e-/ Å3) dans la cage sodalite à 2,49 Å de Co1', le cation Na1 en site I
(ρ=8,11 e-/ Å3), le cation Co2 en site II (ρ=6,91 e-/ Å3), la molécule d'eau Ow1 dans la
supercage à 2,1 Å de Co2 (ρ=6,65 e-/ Å3) dont l’occupation est contrainte à être égale à celle
de Co2, le cation Na2 (ρ=4,29 e-/ Å3) localisé dans la supercage sur le site II (à 0,72 Å d'Ow1
et à 0,9 Å de Co2). Les cartes Fourier-différence montrent ensuite deux résidus sur les sites
III' : le premier (ρ=2,26 e-/ Å3) est introduit dans la structure comme un cation cobalt Co3, le
second (ρ=1,56 e-/ Å3) est modélisé comme un cation cobalt Co4. Les Co3 sont coordonnés à
une molécule d’eau Ow3 en site III dont l'occupation (en nombre de molécules d’eau)
contrainte à être égale à celle de Co3. Les cations Co3 et Co4 sont coordonnés à des
molécules d’eau Ow5 situées sur un site général dans la supercage et dont le facteur
d’occupation est contraint à être égal à la somme des occupations de Co3 et Co4.
75
Tableau 3.6 : Paramètres structuraux de Na,Co-X hydraté (cristal 4)
Atome x y z Wyckoff Occ., % U11 ou Uiso U22 U33 U12 U13 U23
Al -0,0534(2) 0,0367(1) 0,1243(2) 96 100 0,012(2) 0,010(2) 0,012(2) 0,001(2) -0,003(2) -0,003(2)
Si -0,0527(1) 0,1258(2) 0,0358(1) 96 100 0,015(2) 0,014(2) 0,011(2) -0,004(1) 0,002(1) -0,003(2)
O1 -0,1080(3) 0,0000(3) 0,1042(4) 96 100 0,030(6) 0,038(6) 0,033(5) -0,012(4) -0,009(4) -0,016(5)
O2 0,0000(3) -0,0022(3) 0,1454(3) 96 100 0,028(5) 0,029(5) 0,016(5) 0,017(4) -0,006(4) -0,009(4)
O3 -0,0293(4) 0,0774(4) 0,0732(4) 96 100 0,044(6) 0,029(5) 0,030(6) 0,008(5) 0,008(4) 0,013(4)
O4 -0,0736(3) 0,0747(3) 0,1779(3) 96 100 0,033(5) 0,018(5) 0,018(5) 0,000(4) 0,001(4) -0,015(4)
Na1 0,0000 0,0000 0,0000 16 25(3) 0,025
Co1' 0,0700(2) 0,0700(2) 0,0700(2) 32 44(2) 0,042(4)
Co2 0,2262(5) 0,2262(5) 0,2262(5) 32 18(2) 0,03
Na2 0,245(8) 0,245(8) 0,245(8) 32 30(2) 0,03
Co3 -0,169(1) 0,062(1) 0,189(1) 96 8,7(6) 0,05
Co4 -0,159(2) 0,073(2) 0,222(2) 96 6,3(5) 0,05
Ow1 0,267(1) 0,267(1) 0,267(1) 32 18(2) 0,1
Ow2 0,1689(4) 0,0811(4) 0,0811(4) 32 100 0,050(5)
Ow3 -0,195(6) 0,125 0,125 48 17,4 0,14(5)
Ow5 -0,157(3) 0,019(3) 0,231(3) 96 15 0,07(3)
76
Cristal 5
L'étude des structures des cristaux 1 à 4 nous informe sur les transformations
structurales accompagnant la désorption des molécules d'eau : sites préférentiels de sorption
de H2O, modification de la distribution cationique, relaxation du squelette. Ces modifications
structurales sont liées les unes aux autres : la présence de molécules d'eau entraîne une
migration des cations, qui affecte à son tour fortement le squelette de la zéolithe. L'étude de
l'irréversibilité de l'adsorption d'eau (ou d'autres molécules) présente un intérêt fondamental
pour la compréhension du comportement structural des zéolithes industrielles lors de leur
activation. Pour cette raison, nous nous sommes intéressés à la structure à l'équilibre d'une
zéolithe réhydratée après activation totale (cristal 1 réhydraté ~1 mois à l'atmosphère
ambiante).
77
Tableau 3.7 : Paramètres structuraux du cristal 5 (réhydraté)
Atome x y z Wyckoff Occ., % U11 ou Uiso U22 U33 U12 U13 U23
Al -0,0529(1) 0,1237(1) 0,0368(1) 96 100 0,021(1) 0,018(1) 0,019(1) -0,001(1) -0,000(1) -0,003(1)
Si -0,0527(1) 0,0355(1) 0,1253(1) 96 100 0,024(1) 0,021(1) 0,020(1) -0,001(1) -0,002(1) -0,003(1)
O1 -0,1054(3) -0,0000(3) 0,1087(3) 96 100 0,036(4) 0,043(4) 0,037(4) -0,016(3) 0,009(3) -0,017(3)
O2 -0,0042(3) -0,0021(3) 0,1458(3) 96 100 0,036(4) 0,037(4) 0,026(4) 0,016(3) -0,006(3) -0,004(3)
O3 -0,0310(3) 0,0699(3) 0,0748(3) 96 100 0,031(4) 0,041(5) 0,041(4) 0,013(3) 0,008(3) 0,015(3)
O4 -0,0737(3) 0,0729(2) 0,1739(3) 96 100 0,039(4) 0,026(4) 0,029(4) 0,002(3) -0,003(3) -0,012(3)
Na1 0,000 0,000 0,000 16 43(4) 0,025
Co1' 0,0700(4) 0,0700(4) 0,0700(4) 32 17(1) 0,036(7)
Co2 0,2190(3) 0,2190(3) 0,2190(3) 32 18(1) 0,030(6)
Co3 -0,1535(4) 0,0348(4) 0,1693(4) 96 15,5(7) 0,041(4)
Co4 -0,1544(8) 0,0786(9) 0,2163(8) 96 8,3(8) 0,042(9)
Ow1 0,266(2) 0,266(2) 0,266(2) 32 17 0,06
Ow2 0,0802(4) 0,0802(4) 0,1698(4) 32 100 0,063(5)
Ow4 -0,192(3) 0,166(3) 0,150(3) 96 15,5 0,09(2)
Ow5 -0,220(3) -0,024(3) 0,162(2) 96 23,8 0,11(2)
78
III.3.5 Distribution des cations
Selon les études de Takaishi (1996), l'affinité des cations pour les différents sites est
gouvernée par les interactions colombiennes entre les cations et les oxygènes de la charpente
zéolithique. Ce point peut être discuté sur la base du Tableau 3.8 qui résume la distribution
des cations sur les sites et leur coordination.
Cristal 1
Les cations Co1 occupent presque complètement le site I (94,2(6) %) où ils sont
coordonnés à six oxygène du squelette O3 avec une distance courte de 2,206(3) Å (Figure 3.7
a). Le site I' voisin situé sur l'axe 3 est occupé à 4,5(3) % par des cations cobalt Co1'; ceux-ci
sont localisés à 1,324(1) Å au dessus du plan formé par 3 oxygènes O3 (plan 6MR) (Figure
3.7b). L’occupation globale des deux sites est de ~98,7 %.
a) b)
Figure 3.7: Vue ORTEP (visualisation des ellipsoïdes d’agitation thermique des atomes): le prisme hexagonal le
long de l’axe 3 ; a) coordination des cations Co1 (site I), ; b) coordination des cations Co1' (site I').
Les cations cobalt occupent 65,7 % du site II et se situent dans le plan O2-O2-O2 du
6MR (Figures 3.8 a, 3.8 b) où ils sont coordonnés à trois oxygènes O2 (2,067(4) Å). 25(2) %
du site II est également occupé par des cations sodium Na2 localisés à 0,842 (1) Å au-dessus
de ce plan (Figure 3.8 c), l'occupation totale du site est de 91(2) %.
79
a) b)
c)
Figure 3.8: Vue ORTEP de la coordination de cations Co2 et Na2 (site II); a, c) plan contenant l’axe 3 et b) plan
perpendiculaire à l’axe 3.
Le reste des cations sodium est distribué dans la supercage sur les sites III qu'ils
occupent à hauteur de 16% seulement. Leur coordination avec deux oxygènes de la charpente
O4 est moins favorable, avec une distance longue de 2,47(3) Å, expliquant sans doute le
facteur d'agitation thermique isotrope important de Na3 (Figure 3.9).
a) b)
Figure 3.9: Vue ORTEP de coordination de Na3 (site III); a) plan perpendiculaire à l’axe 2 et b) plan contenant
l’axe 2
Cristaux 2, 3, 4 et 5.
Les informations sur la coordination des cations sont résumées dans le Tableau 3.8. Le
centre du prisme hexagonal est occupé par des cations cobalt Co1 dans les cristaux 2, 3 et par
des cations sodium Na1 dans les cristaux 4 et 5 (occupations de site respectives de 38,4(6) et
80
19,8(6) % pour Co1, 25(3) et 43(2) % pour Na1, avec des distances de coordination aux six
oxygènes O3 de 2,457(5), 2,589(5), 2,75(1) Å et 2,656(7) Å respectivement). Le site I' voisin
est occupé par des cations cobalt Co1' dont la population est également affectée par l’état
d’hydratation de la zéolithe (occupations de 31,8(6), 24,9(6), 44(2) et 17(1) % pour des
cristaux 2, 3, 4 et 5 respectivement). Le déplacement des cations cobalt (Co1') vers la cage
sodalite au-dessus du plan formé par les trois oxygènes O3 du 6MR (1,084(1), 1,034(5),
1,28(1) et 1,40 Å) est lié à sa coordination avec la molécule d’eau Ow2 localisée dans les
structures hydratées (cristaux 4 et 5) et située dans la cage β.
Approximativement 1/3 du site II est occupé par des cations Na+ dans les cristaux 2
(occ.= 31(1) %), 3 (occ.= 39 %) et 4 (occ.= 30(2) %). Dans la structure réhydratée (cristal 5)
les cations sodium ne sont pas localisés sur le site II. Les cations cobalt sont localisés dans
toutes les structures sur les sites II avec des occupations de ~55 % : 54,6(6) % et 57,9(6) %
pour les cristaux partiellement déshydratés (2 et 3 respectivement) et seulement 18 % dans les
structures hydratées (cristaux 4 et 5). Pour les cristaux 2 et 3, l’affinement montre que 2/3 des
ions Co2 interagissent fortement avec une molécule d’eau Ow1 localisée dans la supercage
(distance de coordination Co2-Ow1 ~1,84(2) et 1,90(2) Å). Dans les structures hydratées 4 et
5, l’affinement, bien que moins stable, suggère que tous les Co2 sont coordonnés à une
molécule d’eau Ow1 dans la supercage et une molécule d’eau Ow2 dans la cage sodalite avec
des distances de coordination de Co2-Ow1~2,1 Å et Co2-Ow2 ~2,48(4) Å ou Co2-
Ow2~2,13(1)Å.
Les derniers cations localisés occupent partiellement les sites III'a et III'b dans la
supercage. Dans la structure la moins hydratée (cristal 2), les cations Co3 occupent 5,5(5) %
du site III'b près de l’anneau à 4 chaînons du prisme hexagonal et sont coordonnés aux
oxygènes O1 de la charpente avec une distance de 2,58(3) Å. Dans le cristal 3, ~ 6 cations
Co3 (occ.=6,1(3) %) sont localisés sur le site III'a où ils sont coordonnés à deux oxygènes de
la charpente O1 et O4 à 2,30(3) et 2,28(3) Å respectivement, et à une molécule d’eau Ow4 à
1,8(1) Å. Dans les structures hydratée/réhydratée (cristaux 4 et 5), les cations Co2+ occupent
deux types de site III'a. Le premier est occupé par Co3 (8,7(6) % et 15,5(7) % pour les
cristaux 4 et 5) qui est coordonné, dans la structure hydratée, à l’oxygène de la charpente O4
et aux molécules d’eau Ow3 et Ow5 (avec les distances de 2,39(3), 2,35(5) et 2,59(8) Å
respectivement). Dans la structure réhydratée, Co3 est coordonné à des oxygènes O1, O4 de
la charpente et aux molécules d’eau Ow4 et Ow5 (à 2,12(1), 2,22(1), 2,22(8) et 2,39(7) Å
respectivement). Le second site III'a est occupé par Co4 (6,3(5) % et 8,3(8) % pour les
81
cristaux 4 et 5) qui est coordonné aux oxygènes O1 et O4 de la charpente et à une molécule
d’eau Ow5 (avec des distances de 2,31(4), 2,44(4) et 2,14(9) Å pour le cristal 4 et de 2,18(2),
2,28(2) et 2,39(7) Å pour le cristal 5).
Bae & Seff (1999) ont également étudié par diffraction X les structures de 3 cristaux
de Na,Co-X partiellement échangés (Tech=23, 50, 80 °C) et déshydratés (Tdesh=~293K, Pdesh
~0,1 Pa, ∆t=3 jours). Dans ce travail, les auteurs constatent un excès systématique du nombre
de cations compensateurs respectivement à la charge anionique de la charpente. Afin de
conserver l’électroneutralité, un certain nombre de groupements OH- est plus ou moins
localisé dans les structures; en plus des cations H3O+ localisés dans les cages sodalites du
cristal échangé à 50°C. Comme dans notre étude, les espèces extracharpente n’occupent que
partiellement leurs sites, et Bae & Seff (1999) ne différencient pratiquement la nature de
l’espèce que par l’analyse de leurs distances avec les oxygènes. En particulier, considérant les
distances Co2+-Ocharpente de l’ordre 2,1-2,4 Å, ils identifient les H3O+ comme occupant
préférentiellement les axes 3 où ils sont coordonnés à trois oxygènes de la charpente par des
liaisons hydrogènes (H3O+-Ocharpente ~2,9-3 Å) (Zhu et al., 1999). Les OH- sont coordonnés à
des ions Co2+ avec des distances à la charpente supérieure à la distance entre l’eau et
l'oxygène de la charpente.
Si on considère les occupations des sites trouvés dans ces deux études, on constate
que:
L’occupation du site I (prisme hexagonal) par les cations sodium est similaire pour le
cristal de Bae & Seff et le cristal hydraté 4. On retrouve qualitativement les mêmes
paramètres d'occupation, facteur d'agitation thermique et longues distances de coordination
dans les deux cas.
82
Le site I' est partiellement occupé par des Co2+ coordonnés à 3 oxygènes O3 de la
charpente et une molécule d’eau dans la cage sodalite (occ.=32(2)% pour Bae et Seff,
occ.=44(2) % dans cristal 4). Dans la publication de Bae & Seff, ces ions sont plus fortement
coordonnés avec des oxygènes de la charpente (distance de 2,25(2) Å contre 2,49(1) Å pour le
cristal 4), tandis que leurs distances de coordination avec l’eau sont similaires (2,42(5) et
2,49(1) Å pour Bae et Seff et le cristal 4 respectivement). La différence principale est
l’occupation partielle du site I' par des cations H3O+ coordonnés à 3×O3 à 2,9(1) Å et une
molécule d’eau à 2,3(2) Å dans la structure de Bae & Seff, 1999. Pour Bae et Seff, les
molécules d’eau de la cage sodalite coordonnent les Co2+(I') et H3O+(I') tandis que dans le
cristal 4, Ow2 coordonne Co2+(I') et Co2+(II).
Le site II' dans la cage sodalite est partiellement occupé (23,4(5) %) par des ions H3O+
situés à 3,0(2) Å des 3 oxygènes O2, dans les travaux de Bae et Seff, contrairement à nos
résultats.
Le site II abrite des cations cobalt et sodium dont l'occupation et la coordination sont
en relativement bon accord pour le cristal 3 et celui de Bae & Seff (1999) (Tableau 3.9).
Les sites III' de la supercage sont occupés par des cations Co2+ (occ.=13,5(9) %)
coordonnés à deux oxygènes de la charpente (Co-O1=2,26(3) et Co-O42,28(3 Å), deux
molécules d’eau à (Co-O8=1,8(1) et Co-O8=1,9(1) Å) et un ion OH- (Co-OOH-=2,1(1) Å)
selon Bae & Seff. Si on considère que cet échantillon a un état d’hydratation intermédiaire
aux cristaux 3 et 4, cela indique une tendance à une occupation progressive de la supercage
par un nombre de cations solvatés croissant lorsque l'hydratation de la zéolithe augmente
(Tableau 3.8).
En conclusion, si on considère que le cristal étudié par Bae & Seff a un état
d’hydratation intermédiaire entre les états des cristaux 3 et 4, on a globalement un bon accord
concernant les distributions des cations Na+ et Co2+ sur les différents sites et le désaccord
principal porte sur la présence d'ions H3O+ et OH- dans la structure. La présence d'ions
hydronium/hydroxyl peut s'expliquer par le procédé d'échange cationique en colonne utilisé
par Bae et Seff, qui favorise plus l’hydrolyse de la charpente que l’échange cationique en bain
utilisé dans ce travail. Néanmoins, il faut souligner une fois de plus, qu'étant donné le
désordre chimique et positionnel, il est difficile de distinguer par DRX la nature des espèces
extracharpente, en particulier entre H2O, H3O+ et OH-. Dans le cristal 4, les 32 molécules
d’eau, Ow2, de la cage sodalite sont situées à proximité des cations H3O+ (II') localisés par
83
Bae & Seff. Une partie des molécules Ow2 coordonne les cations Co2+(I') et Co2+(II), et il est
possible, que certains sites soient occupés en fait par des ions H3O+(II'). Néanmoins, les
facteurs d’agitation thermique élevés des hydroniums (Uiso=0,14(9) Å2) vont à l'encontre de
la forte stabilisation des H3O+(II') par des interaction de type "liaison hydrogène" avec les
oxygènes O2 de la charpente. L'attribution préférentielle des types chimiques H3O+, OH- et
H2O sur la base seule de la DRX est très discutable (voir aussi l’étude de la distribution des
ions d’hydronium dans Na,H3O-X par diffraction des rayons X et des neutrons par Zhu et al.,
1999).
84
Tableau 3.8 : Distribution des cations et des molécules d’eau dans les 5 structures Na,Co-X comparées à celle publiée par Bae & Seff (1999) ;les occupations des sites sont en
%, les distances interatomiques sont en (Å).
Site Cristal 1 Cristal 2 Cristal 3 Cristal (Bae & Seff, 1999) Cristal 4 Cristal 5
Ion Occ. Coordination Ion Occ. Coordination Ion Occ. Coordination Ion Occ. Coordination Ion Occ. Coordination Ion Occ. Coordination
I Co 94,2(6) 6×O3 2,206(3) Co 38,4(6) 6×O3 2,457(5) Co 19,8(6) 6×O3 2,589(5) Na 20,0(6) 6×O3 2,721(12) Na 25(3) 6×O3 2,75(1) Na 43(2) 6×O3 2,656(7)
I' Co 4,5(3) 3×O3 2,24(2) Co 31,8(6) 3×O3 2,202(6) Co 24,9(6) 3×O3 2,251(7) Co 32(2) 3×O3 2,249(15) Co 44(2) 3×O3 2,49(1) Co 17(1) 3×O3 2,52(1)
Tableau 3.9 : Distances interatomiques (Å) et angles de valence (°) pour les 5 cristaux étudiés comparés à celui publié par Bae & Seff (1999).
85
III.3.6 Effet de la déshydratation sur Na,Co-X.
Les Tableaux 3.8 et 3.9 rassemblent les caractéristiques structurales des cristaux 1-5
ainsi que celles obtenues par Bae & Seff (1999) sur un monocristal de Na,Co-X, de
composition similaire aux nôtres (déshydraté pendant 3 jours à 23°C et ~ 0,1 Pa). La
comparaison de différentes structures permet de conclure que la déshydratation sous vide et à
la température ambiante ne favorise que le départ des molécules d'eau les plus faiblement
liées, comme dans le cas du cristal 3 et la structure de Bae & Seff (1999). Le chauffage à
150°C stimule la désorption des molécules d'eau près de Co3 (cristal 2) mais n'est pas
suffisant pour obtenir une déshydratation complète qui nécessite un chauffage à 400° C
(cristal 1).
Déshydratation de Na,Co-X
• Migration de Co2+ : sites I' → I et Na : sites I →(II)→ III
En fonction de la déshydratation, une migration des cations cobalt du site I' vers le site I
(Figure 3.10a, Tableau 3.8) a été observée. Cet effet a été déjà remarqué expérimentalement
par Smolin (1989) sur la zéolithe Ca-X et par simulations Monte-Carlo (Chapitre III.5) et
s'interprète sur la base d'une modification des sphères de coordination des cations. Dans la
86
structure hydratée (cristal 4), les cobalts en site I' ont une coordination favorable avec trois
oxygènes de la charpente O3 et des molécules d’eau Ow2 situées dans la cage sodalite. Au
cours de la déshydratation, la cage sodalite se vide progressivement des molécules d’eau,
diminuant ainsi le nombre de coordination des cations. On observe un déplacement progressif
des cations Co du site I' vers le plan formé par 3×O3 du 6MR (raccourcissement des liaisons
Co1'-O3, Tableau 3.8) accompagné d'une migration vers le site I au centre du prisme qui offre
une coordination très favorable avec 6 oxygènes O3 (Tableau 3.8).
a) déshydrataion b) déshydrataion
Occupation des sites (Co ), %
100 100
2+
80 80
Site II
Site I
60 60
40 40
20 20
Site I' Site III'
0 0
4 Bae & Seff 3 2 1 1
Bae2 & Seff 33 4 5
4 2 1
Cristal Cristal
Figure 3.10 : La redistribution des cations cobalt au cours de la déshydratation des cristaux 4→1 et celui de Bae
et Seff (1999) sur les sites a) I/I' et b) II/III'.
De leur côté, les cations sodium sont repoussés du centre du prisme hexagonal (site I)
vers la supercage (sites III) probablement en passant par les sites II.
87
charpente, une molécule d’eau dans la supercage et une molécule d’eau dans la cage sodalite)
et dans la structure intermédiaire partiellement déshydratée (cristal 3), les cobalts sont
coordonnés à 3×O2 de la charpente et une molécule d’eau dans la supercage, molécule d’eau
qui est la dernière à être désorbée de la structure (Tableau 3.8).
Ainsi, pour les structures partiellement hydratées et déshydratées les plus précisément
décrites (cristaux 3, 2 et 1), la désorption des molécules d'eau n'affecte pas particulièrement
l'occupation du site II, mais s'accompagne d'un raccourcissement des distances d'interaction
entre les cations Co2 et les oxygènes O2, O4 de la charpente (Figure 3.11).
a) b) c)
Figure 3.11 : Adoption de la géométrie plane par Co2 au cours de la déshydratation: (a) cristal 3, (b) cristal 2, (c)
cristal 1.
On peut également noter que la migration des cations dépend, d’une part, de la
désorption des molécules d’eau induisant une de-solvatation des cations et leur réarrangement
et, d’autre part, des interactions de répulsion/attraction entre les cations et les atomes de la
charpente.
Réhydratation de Na,Co-X:
Dans la structure totalement déshydratée, les cations cobalt occupent principalement
les sites I/II tandis que les cations sodium sont sur les sites II/III. Au cours de la réhydratation
de la zéolithe, on observe les tendances suivantes :
88
cations sont coordonnés à des oxygènes de la charpente O3 et les molécules d’eau Ow2 à une
distance plus longue de ~2,5 Å. En parallèle, les cations sodium, dont beaucoup ne sont pas
localisés en raison de leur désordre, se déplacent du site II vers le site I laissé vacant où ils
adoptent une coordination 6×O3 (avec une distance de 2,656(7) Å). Cette tendance de
redistribution des cations Co2+/Na+ est en accord avec les résultats des calculs Monte-Carlo
(Chapitre III.5).
Figure 3.12 : Migration des cations. Déshydrataion: Co2+ suivent la flèche rouge et Na+ la flèche verte ;
Réhydratation: Co2+ suivent la flèche verte et Na+ la rouge.
89
III.3.6.2 Distorsion du squelette lors de la déshydratation.
Le rôle de la distorsion angulaire des tétraèdres composant les sodalites dans les
relaxations structurales des zéolithes a été étudié par Depmeier (1984). D'après ce dernier, les
petites fenêtres 4MR supportent un stress structural plus important que les plus grandes
fenêtres comme 6MR. Dans les structures de type "FAU", la relaxation accompagnant
l'adsorption/désorption de molécules invitées passe ainsi par l’ouverture/fermeture en
alternance des angles de rotule T-O-T de 4MR ; l’amplitude des variations angulaires dépend
des répulsions entre oxygènes de la fenêtre 6MR voisine. Dans les cristaux de Na,Co-X
étudiés, les fenêtres 4MR composées de Al, Si, O3 et O4 (Figure 3.8 a) concentrent
effectivement les distorsions angulaires maximales qui s’accentuent avec la déshydratation
(Tableau 3.9). Les distances On…Om (n, m=1…4) demeurent quasiment constantes pendant la
déshydratation et la fermeture des angles Si-O-Al vient d'une élongation des liaisons Al-O
tandis que l’ouverture des angles de rotule (en particulier Si-O4-Al) passe par le
raccourcissement des liaisons Si-O (Tableau 3.9).
90
Figure 3.13 : Tétraèdre TO4 régulier ( 4 3m ) où toutes les distances T-O sont égales.
Les interactions entre les cations les plus chargés Si4+, Al3+ et les oxygènes aboutissent
aux liaisons les plus fortes du cristal qui assurent sa cohésion et la formation de la charpente.
La charpente est solide mais elle n’est pas rigide : ainsi sa relaxation peut être possible par les
rotations des tétraèdres TO4 autour d’un axe de symétrie (Figure 3.13) mais ces rotations sont
coûteuses en énergie et se produisent après les modifications des angles de rotule. Dans les 5
structures de Na,Co-X étudiées, on ne note de rotations significatives que pour les cristaux 1
et 2 les plus déshydratés, où les cations Co2+ sont moins écrantés par des molécules d’eau.
Dans la structure déshydratée (cristal 1), les quatre angles de rotule Si-O-Al valent
128,0, 130,4, 119,9 et 158,5° à comparer avec la valeur idéale de 145° pour une structure
relaxée ou avec la structure hydratée ( <Si-O-Al> = 138,5°, Tableau 3.9). Les contraintes
engendrées lors de la déshydratation et la fragilisation de certaines liaisons T-O suite à leur
hydrolyse lors de l’échange cationique expliquent sans doute l’amorphisation partielle et
irréversible des cristaux qui se traduit par un élargissement systématique des pics de Bragg
(Figure 3.14).
91
1,0
Réflexion 1 5 7 Réflexion 6 16 18
1,0
Intensité normalisée
Intensité normalisée
0,8
0,8
0,6
0,6
0,4
0,4
0,2 0,2
0,0 0,0
6 8 10 12 14 16 18 20 22 8 10 12 14 16 18 20 22 24 26 28
1,0
Réflexion 7 5 13
Cristal 1
Intensité normalisée
0,8
Cristal 2
0,6 Cristal 3
Cristal 4
0,4
Cristal 5
0,2
0,0
6 8 10 12 14 16 18 20 22
Dans le cas de zéolithes de type FAUjasite, les synthèses ne sont faites qu'en présence
de cations sodium qui jouent le rôle de SDA. Les structures de zéolithes Mn+-X publiées
(Smolin et al., 1989; Yeom et al., 1997; Kim et al., 1997; Porcher, 1998; Bae & Seff, 1999;
Kim et al., 1999; Bae et al., 2003) sont obtenues par échange cationique nNa+→Mn+ à partir
de zéolithe Na-X. Le Tableau 3.10 compile les paramètres de maille et les angles de rotules
pour différentes structures Mn+-X (n=1, 2) complètement déshydratées (pour lesquelles les
structures sont plus précises que pour les structures hydratées). Les déviations des angles de
rotule T-O-T de la valeur idéale de 145° (Baerlocher et al., 2001) (Tableau 3.10) reflètent à la
fois la présence d'aluminium dans la structure et les forces colombiennes entre cations et
oxygènes de la charpente. Comme attendu, le stress supporté par la structure native Na-X
(Porcher et al., 1999) est moins important que par les structures échangées.
92
Par contre, l'échange avec les cations de plus petit rayon ionique que Na+ comme Li+ou Co2+
entraîne une diminution nette du paramètre de la maille (~-1,6% pour Li-X et ~-2% pour
Na,Co-X); la relaxation de la structure étant assurée à la fois par des modifications des angles
de valence et des rotations des tétraèdres TO4.
Tableau 3.10: Paramètres de maille (Å) et angles de rotule (°) pour des zéolithes Mn+-X totalement déshydratées.
Il est à noter que la diminution du volume de la maille est significative dans le cas de Na,Co-
X (cette étude) mais quasi inexistante pour des zéolithes Na-X (Porcher, 1998, Olson, 1970),
Ca-X (Smolin et al., 1989), Na,K-X, Na,Rb-X, Na,Cs-X (Shepelev et al., 1991). Cette
observation peut s’expliquer par le fait que les interactions entre les oxygènes de la charpente
et les cations compensateurs de charge sont d’autant plus fortes que le rapport charge/taille
des cations est élevé
Tableau 3.11: Distribution cationique dans les structures Mn+-X totalement déshydratées
Cette distribution des cations mono/divalents peut être reproduite par des simulations
93
Monte-Carlo. Dans les structures FAUjasite déshydratées, les sites les plus favorables sont
successivement les sites I puis II et I'. Dans les structures monocationiques divalentes (M2+-X)
tous les cations sont distribués entre les sites I/II qui suffisent à les accueillir. Lorsque le
nombre de cations augmente (structures monocationiques monovalentes M1+-X ou
bicationiques mono/divalentes M1+,M2+-X), les sites I et II sont occupés en premier lieu, puis
certains sites I', la répulsion électrostatique entre cations entraînant un déplacement des
cations du site I vers le site I' opposé, sur l'autre face du D6R (Di Lella, 2007). Au final, dans
les structures monocationiques avec des cations monovalents M1+-X, seuls les sites I'/II sont
occupés et la répulsion électrostatique entre cations gouverne l'affinité des sites (Buttefey et
al., 2001, Beauvais et al., 2004).
94
Site I: L’occupation et les distances de coordination des cations en site I sont
très similaires pour toutes les structures. Dans Na,Co-X, le site I est occupé par Na+ et par
Ni2+ et Zn2+ dans Ni-X et Zn-X, mais dans ces dernières, 4 Na+ et 2,3 Na+ respectivement sont
prévus par des analyses chimiques et n'ont pu être localisés par DRX et peuvent
éventuellement partager le site I.
Site I': Les métaux de transition occupent ~1/3 du site I' : 14,1(6) et 10,4(7)
ions Co dans les cristaux de Na,Co-X n°4 et étudié par Bae & Seff (1999); 8,4(5) Ni2+ dans
2+
Ni-X et 14,8(9) Zn2+ dans Zn-X. Ils sont coordonnées à des oxygènes O3 de la charpente et
des molécules d’eau dans la cage sodalite avec des distances ~2,45 Å.
Site III': La supercage abrite environ 14 Co2+, 30 Ni2+ et 10 Zn2+ (site III') qui
sont assez fortement coordonnés à des oxygènes O4 et/ou O1 de la charpente à ~2,2-2,4 Å et
des molécules d’eau/ions hydroxyles à 1,8-2,6 Å et 2,1 Å (Na,Co-X de Bae & Seff, 1999).
Dans Zn-X, Bae & Seff (2000) localisent, de surcroît, ~18 Zn2+ au centre de supercage sans
coordination à la charpente.
Une analyse cristallochimique plus précise des structures est difficile en raison du
désordre. On peut cependant entrevoir une tendance globale au dépeuplement des sites I et II
au profit du site I' et III' dans les zéolithes M-X quasiment hydratées ou hydratées.
95
Tableau 3.12: Distribution des espèces extracharpente dans les structures hydratées Na,Co-X (cristal 4 et Bae & Seff, 1999), Ni-X et Zn-X (Bae & Seff, 2000). Les
occupations de site sont données en nombre/maille et les distances en Å.
Site Co,Na-X (crystal 4) Co,Na-Xa (Bae & Seff, 1999) Ni-Xb (Bae & Seff, 2000) Zn-Xc (Bae & Seff, 2000)
Ion Occ. Coordination Ion Occ. Coordination Ion Occ. Coordination Ion Occ. Coordination
I Na 4,05(5) 6×O3 2,75(1) Na 3(1) 6×O3 2,72(2) Ni 2,3(3) 6×O3 2,65(1) Zn 1,7(3) 6×O3 2,70(1)
I' Co 14,1(6) 3×O3 2,49(1) Co 10,4(7) 3×O3 2,25(2) Ni 8,4(5) 3×O3 2,45(2) Zn 6,5(9) 3×O3 2,14(2)
3×Ow2 2,49(1) 1×O5 2,42(5) 1×O5 2,52(5)
H3O 9(3) 3×O3 2,9(1) Zn 8,3(8) 3×O3 2,54(2)
1×O5 2,3(2) 3×O5 2,39(2)
II Co 5,8(6) 3×O2 2,19(1) Co 15,2(9) 3×O2 2,16(1) Ni 4,7(6) 3×O2 2,38(2) Zn 7,3(4) 3×O2 2,18(2)
1×Ow1 2,1(1) 1×O6 2,09(6) Zn 3,4(5) 3×O2 2,28(2)
1×Ow2 2,48(2) H3O 23(2) 3×O2 2,97(3) H3O 22(2) 3×O2 2,94(3)
Na 9,6(6) 3×O2 2,48(4) Na 9(2) 3×O2 2,36(4) 1×O8 2,62(1) 1×O6 2,83(2)
II' H3O 8(2) 3×O2 3,0(2) H3O 18(2) 3×O2 3,07(1)
III Ow3 8,4(5) 1×Co(III') 2,35(5)
III' Co 8,4(6) 1×O4 2,39(3) Co 13,0(9) 1×O1 2,26(3) Ni 20,0(8) 1×O1 2,18(3) Zn 10,4(5) 1×O1 2,19(3)
1×Ow3 2,35(5) 1×O4 2,28(3) 1×O4 2,28(2) 1×O4 2,37(3)
1×Ow5 2,59(8) 1×O7 2,1(1) 1×O6 2,37(6)
Co 6,0(5) 1×O1 2,31(4) 1×O8 1,8(1) 1×O7 2,08(6)
1×O4 2,44(4) 1×O8 1,9(1) 1×O8 2,32(9)
1×Ow5 2,14(9) 1×O9 2,28(9)
Ni 9,6(7) 1×O1 2,09(3)
1×O4 2,22(3)
1×O10 2,31(9)
Cage Ow2 32 1×Co(II) 2,48(2) O5 5,2(4) 1×Co(I') 2,42(5) O5 31(2) 1×Ni(I') 2,52(5) O5 8,3(8) 1×Zn 2,39(2)
Sodalite 1×Co(I') 2,49(1) 1×H3O(I') 2,34(17)
Supercage Ow1 5,8 1×Co(II) 2,1(1) O6 15,2(9) 1×Co(II) 2,09(6) O6 20,0(8) 1×Ni(III') 2,37(6) O6 7(2) 1× H3O 2,83(2)
Ow5 14,4 1×Co(III') 2,59(8) O7 13,0(9) 1×Co(III') 2,10(12) O7 20,0(8) 1×Ni(III') 2,08(6) Non
1×Co(III') 2,14(9) O8 26(2) 1×Co(III') 1,83(11) O8 20,0(8) 1×Ni(III') 2,32(9) Zn 17,8(7) coordonnés
1×Co(III') 1,94(12) 1×H3O(II) 2,62(1) à la
O9 20,0(8) 1×Ni(III') 2,28(9) charpente
O10 9,6(7) 1×Ni(III') 2,31(9)
a
Cristal partiellement deshydraté (23°, ~0,1 Pa, 3jours), groupe d’espace Fd 3 .Les espèces O5, O6 et O8 sont des H2O, O7 sont des OH-.
b
Cristal partiellement deshydraté (23°, ~0,1 Pa, 3jours), groupe d’espace Fd 3 m . Les espèces O5, O7-10 sont des H2O, O6=OH-, 1,9(5) Aluminium est trouvé dans le centre
de la supercage.
c
Cristal hydraté, groupe d’espace est Fd 3 . 4 oxygènes O5 sont des OH-, le reste de O5 et les O6 sont des H2O.
96
III.4 Spectroscopie d’Absorption des rayons X (EXAFS).
Le but de l'étude était de préciser les structures obtenues par diffraction X sur les
échantillons déshydraté et hydraté, en particulier la nature (Co2+ ou autre) des espèces
extracharpente occupant les différents sites. L'EXAFS est une sonde locale, sensible à l'état
d’oxydation et la coordination (distance, nombre et type des ligands) de l'atome absorbant.
Elle permet donc d’accéder à l’arrangement structural à courte distance autour de cet
absorbeur, ici Co2+. Cette technique est donc complémentaire à la diffraction qui apporte des
informations moyennes sur l’environnement des ions Co2+ dans le cristal.
Les mesures d’absorption au seuil du cobalt ont été réalisées sur la ligne synchrotron
BL17C1 à Taïwan, en collaboration avec les Pr. Yu Wang et Dr. I-Jui Hsu. Les spectres ont
été mesurés en fluorescence dans la gamme d'énergie 7400<E<8800 eV avec un temps
d’exposition croissant avec l’énergie afin d’augmenter la statistique de mesure. Les détails sur
l’approche théorique de la technique et le dispositif de mesure sont résumés dans l’Annexe B.
Pic de pré-seuil.
97
environnement octaédrique, d'où une information sur la coordinence de l'absorbeur. La Figure
3.15 montre que les pics au pré-seuil ont approximativement tous la même valeur d’énergie
(E=7709,39 eV) quelque soit l'échantillon. L’intensité du pic, qui est plus forte pour les
échantillons déshydratés que pour l’échantillon hydraté, traduit le passage d'une coordination
tétraédrique (ou pseudotétraédrique) vers une coordination octaédrique dominante dans le cas
hydraté. Cette observation est en accord avec le changement de couleur observé
macroscopiquement pour les échantillons qui passent de bleu (Na,Co-X et Na,Co-LSX
déshydratés) ou gris-bleu (Na,Co-Xbpy), à rose (Na,Co-LSX hydratée) (Figures 2.5, 2.10).
0,12
Na,Co-Xdesh
0,10
Na,Co-LSXdesh
0,08 Na,Co-Xbpy
Na,Co-LSXhydr
µ(E)
0,06
0,04
0,02
Les spectres des standards Co, CoO, Co2O3 et LiCoO2, où Con+ est dans les états
d'oxydation n=0, 2, 3 et 3 respectivement, ont été collectés sur la ligne synchrotron BL17C1 à
Taïwan. La comparaison des spectres des standards et des échantillons zéolithiques (Figure
3.16) confirme clairement qu'en moyenne le cobalt est sous la forme Co(II) dans Na,Co-X.
98
Co
CoO
Co2O3
0,2
LiCoO2
Na,Co-LSXhydr
µ(E)
Na,Co-LSXdesh
Na,Co-X(bpy)
Na,Co-Xdesh
0,0
7705 7710 7715 7720
E, eV
L'analyse des spectres EXAFS est réalisée selon la procédure décrite par
Koningsberger (1988) en utilisant la suite des logiciels ATHENA pour l’extraction du spectre
expérimental et ARTEMIS pour des affinements (Ravel & Newville, 2005, ANNEXE B).
Le spectre mesuré est présenté sur la Figure 3.17 a. Le seuil d’absorption du cobalt
(Eseuil = 7717,39 eV) est déterminé à partir de la dérivée première du spectre (Figure 3.17 b).
7717,39 eV
a) 1,6
1,4
b) 0,20
1,2 0,15
1,0
dµ(E)/dE
0,8 0,10
µ(E), u.a.
0,6
0,4 0,05
0,2
0,0
0,00
-0,2
7400 7600 7800 8000 8200 8400 8600 8800 -0,05
E, eV 7680 7700 7720 7740 7760 7780 7800
E, eV
Figure 3.17 : a) Spectre d’absorption brut de cristaux Na,Co-X; b) dérivée première définissant Eseuil
• Bruit de fond
Nous avons trouvé un rayon de coupure Rbkg=0,98 Å en dessous duquel l’absorption
99
des rayons X est due aux transitions électroniques de l’atome central. Le bruit du fond est
ensuite modélisé par une fonction « spline » comportant 11 nœuds (k=1,25-16 Å-1).
• Normalisation
Les points définissant la fonction linéaire de pré-seuil sont choisis à Epré-seuil_1=-200 eV
et Epré-seuil_2=-50 eV, et pour l’après-seuil nous avons pris Eaprès-seuil_1=100 eV et Eaprès-
seuil_2=700 eV. Après normalisation, on obtient le spectre χ(E) (Figure 3.18).
100
a) b)
c) d)
Une estimation correcte de la gamme de résolution est importante car avec une
fréquence maximale trop faible, la précision ∆R sur les distances de coordination est médiocre
π
( ∆k ⋅ ∆R > ) tandis, qu'inversement, une fréquence maximale déraisonnablement grande
2
entraîne l'apparition de pics parasites qui rendent les affinements ultérieurs difficiles voire
impossibles.
Le spectre EXAFS modèle χcalc(R) est calculé à partir des hypothèses structurales
obtenues par diffraction X sur monocristal (paramètre de maille, coordonnées atomiques). On
génère ainsi des clusters autour de chaque atome de cobalt (Programme ATOM inclus dans
ARTEMIS, Ravel & Newville, 2005), à partir desquels on définit une succession de sphères
de coordination pour les différents cobalts (Tableau 3.13) :
101
Tableau 3.13 Sphères de coordination pour des cations cobalt (cristal 1, DRX).
L'occupation partielle de certains sites cristallographiques par des cations cobalt est
prise en compte en multipliant l’amplitude de χ(k) par le facteur d’occupation. Ce facteur est
calculé sur la base du nombre total de cations cobalt dans la structure, distribués sur les sites.
Les résultats des calculs pour les quatre hypothèses considérées dans l'étude de diffraction
sont rassemblés dans le Tableau 3.14.
Tableau 3.14 : Récapitulatif des facteurs d’occupation de sites par des cations cobalt.
102
III.4.2.3 Stratégie d’affinement.
Un des problèmes posés lors des affinements d'EXAFS est le nombre important de
variables pour un relativement faible nombre d'observations indépendantes, ce qui fait qu'un
modèle affiné et satisfaisant du point de vue structural et statistique n’est pas pour autant
unique et qu'il est donc nécessaire de prendre en considération tous les modèles possibles.
Pour limiter le nombre de variables, on utilise la technique d’isolation spectrale qui consiste à
faire l’affinement des paramètres Rj, Nj, σj2 d'une sphère de coordination en isolant un pic
d'oscillation de χ(R) correspondant. On fixe également l’occupation des sites ainsi que Nj aux
valeurs trouvées par la diffraction des rayons X (Tableau 3.13 et 3.14). Le paramètre S02 est
fixé à 0,8 (valeur tirée d'affinements successifs avec une pondération w variable en profitant
du fait que S02 ne varie pas avec w). Les variables d’affinement restantes sont les paramètres
structuraux ∆Rj, σ2j caractérisant chaque sphère de coordination et l’origine en énergie ∆E0.
Les Figures 3.20 montrent l’amplitude (gauche) et les parties imaginaire et réelle du
spectre expérimental. Les positions des pics principaux du spectre sont marquées (lignes
roses). L’observation de la partie réelle indique que les ondes responsables des pics 1, 2 et 3
ont le même type d’atome diffuseur : les pics 1 et 3 sont en phase et les maxima des pics de la
partie réelle coïncident avec les maxima des pics de l’amplitude. Le pic 2 est en anti-phase
par rapport aux pics 1 et 3, ce qui, d’après Lee & Pendry (1975) est causé par la
prédominance (« shadowing ») de la première sphère de coordination (pic 1).
a) b) c)
Figure 3.20 : Amplitude (a), partie imaginaire (b) et partie réelle (c) du spectre expérimental. Les positions des
pis sont repérées par les traits verticaux.
La Figure 3.21 montre l’allure du spectre expérimental χ(R) utilisé pour les
affinements. Les fenêtres en R* (1,63-2,28 Å, 1,63-2,57 Å, 1,63-3,05 Å et 1,63-3,62 Å)
choisies pour l'isolation spectrale traduisent l'extension progressive des sphères de
coordination incluses dans le modèle structural (Tableau 3.13).
*
Les limites Rmax, Rmin de la fonction fenêtre type « Hanning », dR=0,2 Å (ANNEXE B)
103
Figure 3.21 : Spectre χ(R) expérimental.
Figure 3.22 : Spectre expérimental (en bleu) et modèle (en rouge) après l'affinement n°1. La fenêtre de
Transformée de Fourier est en vert.
104
en site II est coordonné à 3 oxygènes à 1,98(2) Å.
Tableau 3.15 : Paramètres initiaux et variables affinées (en gras) de l'affinement n°1
Pic 2 : 2,31≤R≤2,62 Å
Les sphères de coordination n°1 (Co×6 oxygènes en site I) et n°5 (Co×3 oxygènes en
site I') peuvent contribuer au pic n°2. Etant données la faible différence entre les distances
cristallographiques (Tableau 3.13) et la résolution des données EXAFS, leurs 2 contributions
ne peuvent être résolues. La diffraction des rayons X montrant sans ambiguïté (Chapitre
III.3.3) que le site I est occupé par des cations cobalt, l'incertitude porte donc sur la présence
ou non de Co en site I'. Les deux hypothèses structurales à tester sont l'existence de cations
Co2+ en site I seulement (sphère de coordination n°1 seulement) ou la coexistence de cations
sur les deux sites (sphères n°1 et 5). Dans les deux cas, l'affinement EXAFS montre
clairement qu’il est exclu que des cations cobalt soient coordonnés à six oxygènes, ce qui
contredit le modèle cristallographique pour le site I où le cobalt occupe le centre du prisme
hexagonal (symétrie cristallographique 3 ) avec 6 oxygènes à 2,206(3) Å. D'après l'EXAFS,
les 15 cations en site I par maille sont coordonnés à 3 oxygènes seulement, ce qui traduit
probablement leur déplacement hors du centre du prisme et/ou une déformation locale de ce
dernier. Cette situation (absence de cobalt coordonné à 6 oxygènes et faible occupation du site
I' selon DRX) complique la discrimination entre les 2 hypothèses structurales (contribution
similaire des deux sphères de coordination au spectre, compensation possible du faible taux
d'occupation du site I' par d'autres variables) dont aucune ne peut être privilégiée à ce stade.
L'hypothèse 1 de répartition des cations uniquement sur les sites I et II (sphères de
coordination n°9 et 1) sera vérifiée dans les affinements suivants.
105
Pic 3 :2,62≤R≤3,05 Å
Des affinements préliminaires ayant montré que les variables ∆R, σ2 sont similaires
pour les sphères de coordination n°9 et 10, ces variables ont été contraintes. Le Tableau 3.16
résume les variables fixes et affinées, la Figure 3.23 illustre le progrès des affinements.
L’affinement n'est satisfaisant que si le nombre de coordination de la sphère n°10 est réduit de
3 (cristallographie) à 2, ce qui s'accorde bien avec les ellipsoïdes de déplacement atomique
allongées de O4 d'après la DRX Figure 3.24 suggérant un désordre de ces atomes compatible
avec la baisse apparente de coordination trouvée en EXAFS.
Figure 3.23 : Spectre expérimental (en bleu) et modèle (en rouge) après l'affinement 2.1
106
Tableau 3.16 : Paramètres fixes et variables (gras) l'affinement 2.1
Figure 3.24 : Figure ORTEP montrant la sphère de coordination de Co(II) dans un plan perpendiculaire à l’axe
3 et l'allongement des ellipsoïdes de O4.
La Figure 3.25 montre une image ORTEP du prisme hexagonal obtenu lors des études
DRX (Chapitre III.3) où les cations cobalt se situent au centre du prisme et sont coordonnés à
six oxygènes équidistants (sphère n°1), en contradiction avec les affinements EXAFS
précédents. Si on fait l’hypothèse de déplacement des cations hors du centre du prisme et/ou
de la déformation du prisme, on doit introduire une nouvelle sphère de coordination où les
cations cobalt du prisme sont désormais coordonnés avec une distance longue à trois
oxygènes (Co-Olong).
107
Figure 3.25 : Vue ORTEP du prisme hexagonal avec les cations en site I (Co1, coordination "idéale" à 6
oxygènes équidistants).
Les affinements préliminaires confirment que la variable σ2 est similaire pour les
sphères de coordination n°9 et (Co1-O)long (Tableau 3.17) correspondant respectivement aux
cations en site II coordonnés à trois oxygènes (Affinement 1) et aux cations en site I
coordonnés à trois oxygènes avec « longue distance ». D'un point de vue statistique, les
affinements 2.1 et 2.2 sont quasiment identiques (Tableaux 3.16 et 3.17, Figures 3.23 et 3.26)
et il est difficile de privilégier sur cette base seule l'une des deux hypothèses. Au final,
l’affinement 2.1 sera utilisé pour la suite.
Figure 3.26 : Spectre expérimental (en bleu) et modèle (en rouge) après l'affinement 2.2
108
Tableau 3.17 : paramètres fixes et variables (gras) de l'affinement 2.2
109
Figure 3.27 : Spectre expérimental (en bleu) et modèle (en rouge) après l'affinement 3
Bilan :
Les parties imaginaire et réelle de la transformée de Fourier du spectre EXAFS
(Figures 3.28 a et 3.28 b) permettent de juger si un ou plusieurs types d’atomes contribuent à
la sphère de coordination. Les informations sur le type atomique des diffuseurs proviennent
du fait que les déphasages induits par différents atomes se combinent linéairement dans les
parties réelle et imaginaire, contrairement à ce qui se passe dans le module |χ(E)| (Figure
3.27).
110
a) b)
c)
Figure 3.28 : Spectres χ(k) expérimental (en bleu) et modélisé (en rouge): (a) partie réelle, (b) partie imaginaire,
(c) décomposition suivant les différentes contributions.
Les Figures 3.28 a, b montrent que tous les pics significatifs (dans la fenêtre de
résolution) sont bien modélisés, l’affinement est donc complet et les valeurs des paramètres
structuraux et leurs incertitudes sont satisfaisants. La Figure 3.28 c illustre la décomposition
du spectre χ(k) expérimental en contributions provenant des 4 chemins de diffusion
considérés.
111
Tableau 3.19 : comparaison des paramètres structuraux obtenus par EXAFS et DRX
Etape 1 2 3
d’affinement EXAFS 2 2.1 2.2
Interaction Co2-O2 Co1-O3 Co2-O4 Co1-O3long Co1-O2
EXAFS
Coordination 3 3 2 3 3
R (Å) 1,988±0,007 2,17±0,03 2,913±0,007 2,920±0,007 3,39±0,01
σ2 (Å2) 0,0043±0,0004 0,018±0,004 0,0043±0,0004 0,0043±0,0004 0,0043±0,0004
DRX
Coordination 3 6 3 - 6
R (Å) 2,067(3) 2,206(3) 2,991(4) - 3,663(2)
Figure 3.29 : Spectres k3χ(k) de Na,Co-X (en rouge) et de la poudre LSX déshydratée (en bleu).
112
La Figure 3.30 montre la comparaison des parties imaginaires des transformés de
Fourier des spectres des cristaux de Na,Co-X et de la poudre de Na,Co-LSX déshydratés. Les
pics se trouvent à la même position dans les deux spectres jusqu’à ~3,6 Å, distance après
laquelle on note un désaccord.
Figure 3.30 : Comparaison des parties imaginaires des transformés de Fourier des spectres Na,Co-X et Na,Co-
LSX déshydratés.
Résultats d’affinements
L'affinement a été réalisé en suivant la même procédure que pour les cristaux de
Na,Co-X déshydratée : le bruit du fond est modélisé par une fonction « spline » composée de
11 nœuds (Rbkg=0,98 Å), la gamme de résolution s'étend de 2,5≤k≤11,5 Å-1 (Figure 3.31) et
S02 est fixé à 0,8 comme pour Na,Co-X (même état d'oxydation).
Les résultats sont résumés dans le Tableau 3.20 et illustrés par les Figures 3.31 et 3.32.
Figure 3.31 : Spectre expérimental (en bleu) et modèle (en rouge) de Na,Co-LSX déshydratée
113
a) b)
Figure 3.32 : Partie réelle (a) et imaginaire (b) des spectres de Na,Co-LSXdesh expérimental (en bleu) et
modélisé (en rouge)
La Figure 3.31 montre qu'un facteur d'occupation de 0,6 permet une modélisation
correcte du premier pic (1,63≤R≤2,25 Å) attribué à la sphère de coordination Co-O avec N=3.
Le pic autour de 2,25≤R≤2,59 Å est également modélisé par une interaction Co-O de même
manière que pour Na,Co-X. En revanche, les deux derniers pics (2,59≤R≤3,1 Å et
3,1≤R≤3,63 Å) nécessitent de réduire le nombre de coordination pour garder la même
occupation de site que pour Na,Co-X (Tableau 3.20). Comme pour Na,Co-X, tous les essais
pour introduire les contributions des sphères de coordination Co-Si/Al n’ont pas donné
d'affinements satisfaisants.
Etape 1 2 3
d’affinement 2 2.1 2.2
Interaction Co2-O2 Co1-O3 Co2-O4 Co1-O3long Co1-O2
S02 0,8
∆E0 (eV) 3,3±0,8
Occupation (%) 60 40 60 40 40
N 3 3 1,8 2,5 2
σ2 (Å2) 0,0021±0,0005 0,009±0,001 0,0021±0,0005 0,0021±0,0005 0,002±0,001
R (Å) 1,99±0,02 2,20±0,01 2,96±0,01 2,96±0,01 3,43±0,02
114
...) ont donné des résultats moins bons, voire une divergence de l’affinement.
La Figure 3.33 montre que les spectres expérimentaux des poudres Na,Co-LSX
réhydratée (en rouge) et déshydratée (en bleu) sont très différents et qu'on ne peut, dans ce
cas, utiliser le modèle structural de Na,Co-X déshydratée comme point de départ.
Figure 3.33 : Spectres χ(k) expérimentaux de Na,Co-LSX déshydratée (bleu) et réhydratée (rouge)
Dans le Chapitre III.3 nous avons montré que le désordre affectant la distribution des
cations et des molécules d’eau dans une zéolithe X hydratée limitait drastiquement la
précision de la structure. Pour les analyses EXAFS, la structure de départ est formée de la
charpente du cristal réhydraté (n°5) obtenue par DRX sur monocristal (Chapitre III.3),
combinée avec la distribution cationique obtenue par calcul Monte-Carlo (Chapitre III.5)
(Tableau 3.21). Etant donné le nombre élevé de variables à prendre en compte (occupations
de site, distance et nombre de coordinations), les paramètres d'occupation de site et le nombre
de coordination ont dû être fixés.
115
Résultats
Figure 3.34 : Spectres de poudre Na,Co-LSX réhydratée expérimental (bleu) et modélisé (rouge), contribution de
la sphère de coordination Co-Si/Al (vert).
Tableau 3.21 : Distribution des cations Co2+ sur les sites cristallographiques dans une structure réhydratée
(structure obtenue par combinaison DRX -calculs MC)
DRX EXAFS
Site Occ. (%/Ncations) Coordination Distance (Å) Occ.* (%)
I -
I' 17/5,4 3×O3 2,52(1) 16
3×H2O 2,51(1)
3×Si 3,46(2)
3×Al 3,45(2)
II 18/5,8 3×O2 2,18(1) 17
1× H2O 2,04(6)
1× H2O 2,13(1)
3×Si 3,13(2)
3×Al 3,17(3)
III' 23,8/22,8 1×O1 2,12(1) et 2,18(2) 67
1×O4 2,22(1) et 2,28(2)
1× H2O 2,24(7) et 2,39(7)
3×Si 2,73(2) et 3,49(2)
3×Al 3,56(3) et 2,77(2)
*
L’occupation en EXAFS est calculée comme un pourcentage des cations occupant un site relatif au nombre
total des cations cobalt dans la structure.
116
3×2,18+2,04+2,13 =2,14
Cobalt en site II : d(II) = 5 Å
2,18+2,28+2,39 =2,28
d(III') = 2,12+2,22+2,22+2,24 =2,20 Å et d(III') = Å
4 3
2,14+2,20+2,28 =2,21
d'où la distance moyenne d = 3( d(II) + d(III') + d(III') )=
1
3 Å
L’amplitude de ce pic est sensiblement plus élevée dans le spectre réhydraté (Figure
3.34) que dans celui déshydraté (Figure 3.31), ce qui traduit l’augmentation du nombre de
coordination ou du nombre de cations impliqués. Il peut être modélisé indifféremment avec
des paramètres (occupation, coordination) de (65%, 6) ou (84%, 4,5) ou (94%, 4). Le jeu de
paramètres (84%, 4,5) semble préférable puisque l’occupation correspond à la somme des
occupations des sites II et III′ (Tableau 3.20), qui indique également une coordination 5 des
cations en site II et une coordination 3 ou 4 (5 ou 6 selon les calculs Monte-Carlo) des cations
en site III′.
Le deuxième pic (1,73≤R≤2,8 Å) est nettement moins bien résolu dans le spectre
Na,Co-LSX réhydraté (Figure 3.34) que dans le spectre déshydraté (Figure 3.31),
probablement en raison d'un recouvrement entre pics issus de sphères de coordination avec
des distances proches mais légèrement différentes. Ce pic est affiné de façon satisfaisante
pour des couples (occupation, coordination) de (20%, 6), (30, /5), (40, 4) ou (50, 3).
Les pics suivants (a priori exploitables jusqu’à R≈4,5-5 Å) sont en fait très difficiles à
modéliser (interférences de contributions des diffusions sur des sphères de coordination
différentes pour lesquelles on ne dispose pas d'information de qualité par DRX). Les
diffusions multiples, possibles à cette distance, compliquent également l’analyse. Les
117
statistiques ne sont donc pas satisfaisantes (Tableau 3.22) comme on pouvait l'attendre.
Tableau 3.22 : Récapitulatif des affinements EXAFS de Na,Co-LSX réhydratée (Nind=10,41, Nvar=5, χ2=122,
χν2=22,6, R=0,074)
Etape 1 2 3
d’affinement
Interaction Co-O Co-O Co-Si/Al
Occupation (%) 82 60 50
S02 0,8
E0 -4±2
N 4,5 2 3
σ2 (Å2) 0,005±0,001
R (Å) 2,05±0,01 2,60±0,03 3,29±0,02
118
III.4.5 Bilan des analyses EXAFS
Globalement, on constate un très bon accord entre les structures obtenues par DRX et
EXAFS. Les structures déshydratées Na,Co-X et Na,Co-LSX aboutissent aux mêmes
conclusions:
2) la nature du cation occupant le site I' ne peut pas être déduite directement en
analysant les oscillations attribuables à la sphère de coordination Co(I')-O3 car l’occupation
de ce site est très faible d'après la DRX, et la différence entre les distances Co(I')-O3
(2,206(3)Å) et Co(I)-O3 (2,24(2)Å) est inférieure à la résolution du spectre EXAFS.
L’analyse EXAFS ne permet pas de conclure quant à la nature du cation en site I', une étude
par DAFS (« Diffracted Anomalous Fine Structure », Palancher et al., 2005) serait plus
adaptée ici.
3) il n’y a pas de cations Co2+ coordonnés à six oxygènes mais seulement à trois,
ce qui implique un déplacement des Co2+ hors du centre 3 du prisme hexagonal ou une
déformation du prisme. Cette observation est confirmée par d'autres études expérimentales et
théoriques : Verberckmoes et al. (1997), par exemple, étudiant des Na,Co-X et Y par
Spectroscopie de Reflectance Diffuse concluent que les Co2+ dans le prisme hexagonal se
trouvent en coordination pseudo-octaédrique. Dooryhee et al. (1991) observent ce type de
désaccord sur la zéolithe Na,Ni-Y déshydratée en EXAFS et DRX sur poudre. Ils concluent
que les cations Ni se trouvent en moyenne au centre du prisme hexagonal en coordination
octaédrique (DRX) alors que les analyses EXAFS indiquent localement une coordination
quatre : l’abaissement du nombre de coordination est expliqué par la distorsion du prisme
hexagonal. En se basant sur la faible valeur du facteur d’agitation thermique (étude DRX) des
cations Ni2+ dans le prisme, ils suggèrent que le déplacement des cations par rapport au centre
du prisme est moins favorable que la distorsion. Ce raisonnement peut également s'appliquer
dans le cas de Na,Co-X où le facteur de déplacement atomique des cations Co2+ en site I est
modéré (U=0,0276(5) Å2). L'étude par RMN 23Na de Na-Y par Engelhardt (1997) montre que
dans la zéolithe déshydratée, les cations du prisme font l'objet d'interactions quadripolaires et
sont déplacés par rapport au centre, allant du site I vers des sites Ia et Ib. Les simulations
Monte-Carlo de Heidler et al. (1996), visant à localiser les sites occupés par des cations Na+
et Cs+ dans la structure FAUjasite, montrent également ce déplacement des cations du prisme
hexagonal. Les auteurs concluent ainsi que ces cations sont coordonnés à trois oxygènes O3 et
119
non pas à six. Il est à noter que la majorité des simulations trouvées dans la littérature
concernent des cations monovalents. Des calculs similaires pour des cations divalents (en
particulier Co2+) afin de définir la coordination la plus favorable des cations M2+ dans la
structure Na,M-X seraient particulièrement intéressants.
120
III.5 Simulations Monte-Carlo de la distribution
cationique.
Cette étude est originale sur deux points : elle présente une structure de zéolithe X
bicationique qui possède des cations divalents (cas rarement abordé tant du point de vue de la
modélisation qu’expérimentalement) et elle met en évidence le rôle de l’hydratation et de la
flexibilité du squelette sur la redistribution des cations dans la zéolithe.
121
III.5.2 Champs de forces : interactions entre les espèces.
Dans les calculs les interactions entre particules i et j sont traitées par des potentiels
classiques à deux corps incluant des termes coulombien, de polarisation, de dispersion et de
répulsion :
L’énergie totale du système est évaluée par sommation de ces interactions (Vij) entre
les particules i, j en positions Ri et Rj et leurs images périodiques :
1 N r r r
U= ∑ ∑ Vij ( Ri − R j + n )
2 n i, j
i≠ j
Le terme n permet d’introduire des interactions entre les particules et leurs images
r
périodiques ( n = (nx Lx , n y Ly , nz Lz ) , avec nx, ny, nz entiers et Lx, Ly, Lz les limites de la boîte de
simulation ici égaux au paramètre de maille a). Les interactions électrostatiques sont calculées
par sommation d’Ewald à partir des charges partielles proposées par Mortier (Uytterhoeven et
al., 1992) pour des zéolithes de rapport Si/Al=1 (Tableau 3.24).
Selon les espèces en interaction, les potentiels prennent les formes suivantes:
Cation-Eau : cette interaction est décrite par un potentiel du type Lennard Jones et un
potentiel électrostatique. Les coefficients du potentiel de Lennard Jones sont obtenus en
utilisant les règles de croisement de Lorentz-Berthelot (ANNEXE C) des paramètres Lennard
Jones développés pour des interactions Cation-Cation, Ow-Ow, Ozéolithe-Ozéolithe (Pascual et
al., 2003).
Les constantes des potentiels de type Buckingham et Lennard Jones utilisés sont
données dans le Tableau 3.23, les charges des ions dans le Tableau 3.24.
122
Tableau 3.23 : Paramètres des potentiels
Buckingham
γ
= α ij exp(− β ij rij ) −
Buckingham
Vij
rij6
α (K) β (Å-1) γ (K Å6)
Co-Ozéolithe 6,11×107 4,6 4,06×105
Na-Ozéolithe 6,11×107 4,05 7,652×105
Lennard Jones
⎡⎛ σ ⎞ 6 ⎛ σ ⎞
12 ⎤
Vij = −4ε ij ⎢⎜ ⎥
ij ⎟
−⎜ ⎟
ij
⎢⎜ rij ⎟ ⎜ rij ⎟ ⎥
⎣⎢⎝ ⎠ ⎝ ⎠ ⎦⎥
ε (K) σ (Å)
Co-Co ou Na-Na (Dang, 1995) 50,27 2,584
Ow-Ow (Jorgensen et al., 1983) 78,03 3,1536
Ozéolithe-Ozéolithe 93,53 3,00
Cation
Co +2
Na +1
Zéolithe (Uytterhoeven et al., 1992)
O -0,85
T (T=Si ou Al) 1,2
Eau (Jorgensen et al., 1983)
O -1,04
H 0,52
123
III.5.3 Hypothèses de départ des simulations
Dans les charpentes, le rapport Si/Al est supposé égal à 1 et les atomes T (T=Si, Al) ne
sont pas différenciés. Pour respecter l’électroneutralité, le nombre de cations compensateurs
de charge dans les structures simulées (40 Co2+ et 16 Na+) excède légèrement le nombre de
cations localisés expérimentalement (Si/Al≈1,08 ; 38 Co2+ et 16 Na+). En raison des
approximations faites, en particulier sur la forme des potentiels d'interaction, et de la nature
essentiellement statistique des résultats, les simulations ne donnent que les grandes tendances
des distributions des atomes extracharpente sur les différents sites.
Pour plus de simplicité, la nomenclature des sites (Figure 3.1) utilisée en DRX a été
reprise pour les études de simulation Monte-Carlo. La symétrie n'étant pas prise en compte
dans les calculs, les sites sont définis avec une tolérance de position (distance donnée entre
parenthèses) :
I' : Centre du plan formé par 3×O3 de 6MR du prisme hexagonal (1,09 Å),
124
déjà les premières tendances pour la stabilisation des cations dans les différents sites (Figure
3.35).
Charpente déshydratée
Charpente hydratée
50000
Energie électrostatique, K
0 Supercage
Sodalite
-50000 Prisme
Figure 3.35 : Coupe d’énergie électrostatique le long de l’axe 3 des charpentes hydratée et déshydratée.
Dans la structure déshydratée, les sites (I, I') du prisme hexagonal sont nettement
stabilisés par rapport aux autres. Les autres sites les plus stables de la structure déshydratée
(II'/II) et de la structure hydratée (I'/II'/II) ont des énergies comparables.
La Figure 3.36a représente les coupes selon l’axe 3 de l'énergie totale (Buckingham et
Electrostatique) et de la contribution de Buckingham seule.
125
a) b)
50000 Cobalt 50000 Prisme Cobalt
Energie (Buckingham), K
Sodium 40000 Sodium
0 Supercage
30000
Sodalite
Energie, K
-100000
10000 Sodalite
II' / II II / II' Supercage
I' I' 0
-150000
I -10000 II' II II II'
I' I I'
0 10 20 30 40 0 10 20 30 40
Distance le long de l'axe 3 (Å) Distance le long de l'axe 3 (Å)
Figure 3.36 : Coupe d’énergie le long de l’axe 3 pour les cations Co2+ et Na+ dans une zéolithe déshydratée: a)
contributions des énergies de Buckingham et Electrostatique; b) énergie d’interaction dispersion-répulsion
(Buckingham).
Globalement, les occupations de sites et les distances de coordination des cations Co2+
et Na+ trouvées par la simulation sont en très bon accord avec les études expérimentales pour
les sites les plus stables, i.e. I, II (Tableau 3.25, 3.26).
Tableau 3.25 : Comparaison de la distribution des cations dans une structure déshydratée expérimentale (DRX)
et calculée (MC)
126
Tableau 3.26 : Comparaison des distances de coordination (Å) des structures déshydratées expérimentale (DRX)
et calculée (MC)
Cobalt Sodium
Site Coordination Exp. Sim. Exp. Sim.
I 6×O3 2,206(3) 2,21 - -
I′ 3×O3 2,24(2) - - -
II 3×O2 2,067(4) 2,2 2,24(1) 2,33
III/III' 2×O4 - - 2,47(3) 2,54
Les cations en site I' ne sont pas localisés par MC, ce qui s'explique par l’occupation
complète du site I voisin beaucoup plus favorable et l'impossibilité de déformer localement la
charpente pour adapter un nouveau cation. Le faible taux d’occupation des sites III/III' qui
offrent une coordination moins favorable aux cations et la dispersion des Na+ sur de
nombreux sites est en accord avec le désordre constaté en DRX.
Il est important tout d’abord de noter que les essais préliminaires de modélisation de
l’adsorption d’eau dans Na,Co-X assumant la géométrie figée de la charpente déshydratée ne
permettent pas de reproduire la redistribution des cations observée en DRX. De nouveaux
calculs à partir de la charpente de Na,Co-X hydratée ont donc été réalisés afin de voir si la
relaxation du squelette, observée précisément en RX, affectait suffisamment la contribution
électrostatique dominante, qu'elle modifiait les résultats. Les simulations ont été réalisées
dans l’ensemble grand canonique à température ambiante pour des pressions d’eau allant de
10-12 à 100 Pa (0 à ~280 molécules d' H2O par maille). Le système arrive à l’équilibre en ~106
pas de Monte-Carlo avec une tendance au ralentissement lorsque le nombre de molécules
d’eau adsorbée augmente. Les moyennes structurales sont réalisées sur ~103-104
configurations échantillonnées parmi les 106 à 107 configurations à l’équilibre.
La structure avec 243 molécules d’eau adsorbées (T=300 K, P=10 Pa), proche du
nombre déterminé par analyse thermogravimétrique (237) est examinée ici. Si on considère
l’énergie électrostatique uniquement (Figure 3.35), les sites I' et II/II' sont les plus stables
pour les cations. La prise en compte des énergies de dispersion-répulsion (Figure 3.37 b) entre
cations et charpente montre que l’ordre de préférence des sites est I'>I≈II/II' pour Co2+ et
I'≈I>II/II' pour Na+ (Figure 3.37 a).
127
a) Cobalt b) Cobalt
60000 Sodium Prisme
5000 Sodalite Sodium
Energie (Buckingham), K
40000
2500
20000
Supercage
Supercage
Energie, K
0 0
-20000 -2500
-40000 Sodalite -5000
-60000
Prisme -7500 II II
-80000
I' I' II' II'
II' / II II / II' -10000
-100000 I' I I' I
0 10 20 30 40 0 10 20 30 40
Distance le long de l'axe 3 (Å) Distance le long de l'axe 3 (Å)
Figure 3.37 : Coupe d’énergie le long de l’axe 3 pour les cations Co2+ et Na+ dans une zéolithe hydratée : a)
contributions des énergies de Buckingham et Electrostatique; b) énergie d’interaction dispersion-répulsion
(Buckingham).
Il est à noter, néanmoins, que dans le calcul illustré par la Figure 3.37a, les interactions
Lennard Jones entre cations et molécules d’eau ne sont pas prises en compte.
D'après les simulations, les Co2+ sont solvatés en premier et se distribuent au final
(toute énergie comptée) sur les sites I', II, III/III' (Tableau 3.27). Les cations sodium sont
solvatés en dernier, ce phénomène n’affectant que les quelques cations en site II accessibles.
La barrière d’énergie entre les sites I et I' pour Na+ est faible, d’où une occupation
préférentielle du site I où les Na+ trouvent une coordination favorable avec 6 oxygènes de la
charpente ; les cations Co2+ plus chargés sont stabilisés en site I' (Tableau 3.27). La situation
est contraire à celle de la structure déshydratée, et s'explique par la déformation du squelette
au niveau de l'angle de rotule O3 : la distance de coordination Mn+-O3 en site I est trop courte
pour Na+ dans la structure déshydratée (2,21 Å) mais acceptable dans la structure hydratée
relaxée (2,66 Å).
Cette simulation met donc en évidence l’occupation du prisme hexagonal par des
cations sodium et confirme que l'effet « ion-sieving» est en partie responsable de
l’impossibilité de réaliser des échanges cationiques 2Na+→Co2+ complets.
128
Tableau 3.27: Distribution des cations et coordination dans la structure hydratée (simulation).
Site Cation Nb. Cations Charpente Eau, Cage sodalite Eau, Supercage Eau, Prisme
Coord. Dist. (Å) Coord. Dist. (Å) Coord. Dist. (Å) Coord. Dist. (Å)
I Na 6,4 6×O3 2,66
I' Co 9,9 3×O3 2,21 1×H2O 2,48 1×H2O 2,50
II Co 4,7 3×O2 2,2 1×H2O 2,48 1×H2O 2,27
Na 5 3×O2 2,21 1×H2O 2,48
III/III' Co 22,7 1×O1 2,36 3 à 4×H2O 2,40
1×O4 2,3
D'après la simulation Monte-Carlo, dans la structure hydratée, les cobalts en site I' sont
coordonnés à trois O3 de la charpente, une molécule d’eau dans le prisme et une ou deux
molécules d’eau dans la cage sodalite. Les cobalts en site II sont coordonnés à 3×O2, une
molécule d’eau dans la cage sodalite et une molécule d’eau dans la supercage (ou deux
molécules d’eau dans la supercage), les cobalts en site III' à 2 oxygènes O1 et O4 de la
charpente et 3 à 4 molécules d’eau dans la cage sodalite (Figure3.38). Ces informations
complètent la distribution obtenue par la méthode conventionnelle de synthèse Fourier-
différence lors des affinements DRX préliminaires et ont ensuite été utilisées comme point de
départ pour de nouveaux affinements structuraux qui sont restés stables. Au final, les
structures expérimentale et simulée (Tableau 3.28 et 3.29) sont relativement en bon accord,
dans la limite de la précision de la structure expérimentale (désordre structural, cations et
molécules d’eau non localisés).
Tableau 3.28 : Comparaison de la distribution des cations dans les structures expérimentales (DRX) et calculées
(MC=Monte-Carlo)
129
Tableau 3.29 : Comparaison de la coordination des cations dans les structures hydratées expérimentales et calculées
Co2+ Na+
Site Coordination Hydrateé (DRX) (Å) Réhydratée (DRX) (Å) Hydratée (MC) (Å) Hydrateé (DRX) (Å) Réhydratée (DRX) (Å) Hydratée (MC) (Å)
I 6×O3 - - - 2,75(1) 2,656(7) 2,66
I' 3×O3 2,49(1) 2,52(1) 2,21 - - 2,6
1×H2O (sodalite) 2,49(1) 2,51(1) 2,48 - - -
130
L’analyse statistique de la distribution des cations dans la structure hydratée montre
que les cations Co2+ sont bien localisés sur les sites I', II et leur distribution spatiale est plus
large autour du site III'. Statistiquement, les cations Co2+ sont mieux localisés que les Na+, en
accord avec les résultats expérimentaux (Chapitre III.3). Le désaccord principal porte sur les
cobalts en site I'. La distance expérimentale de coordination de Co(I') à O3 (2,49(1) et 2,52(1)
Å pour les structures hydratée et réhydratée respectivement) est plus longue que celle trouvée
par simulation (2,21 Å). La Figure 3.38 montre que, d'après les simulations, les Co(I') sont
coordonnés à des molécules d’eau situées dans la cage sodalite (en accord avec l'expérience)
et dans le prisme hexagonal (molécules non localisées expérimentalement).
a) b) c)
Figure 3.38 : Coordination des Co2+ (en bleu) en site I' (a), II (b), III' (c) dans une structure hydratée calculée
(~243 H2O/maille)
Concernant les molécules d’eau, les quelques 243 molécules par maille sont
distribuées dans la cage sodalite (32 coordonnant Co(I') et Co(II) ), le prisme hexagonal (10
coordonnant Co(I')) et la supercage (203 coordonnant des Co(II), Co(III/III') et le reste en
liaison avec la charpente). Expérimentalement, seulement ~75 molécules par maille sont
localisées, coordonnant toutes les cobalt en sites I', II et III'.
131
Å). Cette distribution est reproduite seulement avec la charpente hydratée (a=24,936 Å). La
flexibilité de la charpente, en particulier au niveau des angles de rotules (T-O-T), permet de
fortes interactions entre les cations et les oxygènes du squelette dans la structure déshydratée
(Tableaux 3.25 et 3.26). Probablement également, la déformation (ouverture) locale des
fenêtres 6MR facilite la diffusion des cations. Idéalement, il faudrait pouvoir suivre la
relaxation de la charpente avec l’augmentation de pression d’eau dans les simulations. Ce
type de calcul n’est malheureusement pas possible actuellement car il nécessite de développer
de nouveaux potentiels d’interaction pour la charpente qui alourdiraient sensiblement les
calculs. La meilleure solution reste actuellement d'utiliser la charpente hydratée déterminée
expérimentalement sur le même type de zéolithe (mêmes cations compensateurs de charge
que dans le calcul) afin d'évaluer la localisation des espèces adsorbées dans les simulations.
Néanmoins, l’utilisation de la charpente hydratée biaise inévitablement la localisation des
premières molécules d’eau adsorbées. L’expérience de diffraction X montre, cependant, que
le volume de la maille du cristal 3 partiellement hydraté (~18,5 molécules d’eau par maille
localisées dans la supercage, Chapitre III.3) n'est contracté que de 1,28% contre 4 % pour le
cristal déshydraté. On peut donc supposer qu’à partir d’une centaine de molécules d’eau
(quasi saturation de la cage sodalite), la charpente relaxée est adaptée si on veut suivre la
tendance de la redistribution des molécules d’eau (Figure 3.39) et des cations.
220
Distribution de molécules d'eau
200 Supercage
180
Sodalite
160
Prisme
140
120
100
80
60
40
20
0
0 50 100 150 200 250
Nombre total de molécules d'eau
Ainsi, à ~100 H2O/maille, une partie des Co2+ en site II s'est déplacé vers les III/III'
par suite de leur solvatation (Figure 3.39, 3.40). La cage sodalite est remplie au ~¾ par des
molécules d’eau mais le prisme hexagonal est vide.
132
25
10
2+
0
0 50 100 150 200 250
Nombre de molécules d'eau par maille
Figure 3.40 : Déplacement des cations Co2+ du site II vers le site III' en fonction du nombre de molécules d’eau
adsorbées
Avec l'augmentation de la pression, le déplacement des cobalts des sites II vers les
sites III/III′ se poursuit (Figure 3.40), la cage sodalite continue de se remplir, de plus en plus
de molécules d’eau rentrent dans le prisme et coordonnent les Co2+ en site I′ qui quittent le
prisme. Finalement le prisme hexagonal est occupé par Na+ et H2O.
133
III.6 CONCLUSIONS
Na,Co-X déshydratée
Figure 3.41 : Distribution cationique (en nombre de cations par maille) sur les sites cristallographiques dans la
zéolithe Na,Co-X déshydratée.
Les études d’Absorption des Rayons X au seuil du cobalt ont révélé l'absence de
cobalts hexacoordonnés dans la structure complètement déshydratée. Le prisme hexagonal est
en fait probablement déformé autour du site I, et certains cations sur ce site sont probablement
déplacés en dehors de la position idéale 3 .
134
Na,Co-X (partiellement) hydratées
Lors de la réhydratation, une petite partie des cations cobalt repeuple le site I' tandis que la
majorité migre vers le site III'. Ce résultat est en accord avec les simulations Monte-Carlo.
Ces dernières montrent le rôle primordial de la relaxation de la charpente et l’adsorption de
molécules d’eau sur la redistribution des cations. L’utilisation de la charpente relaxée
(hydratée) induit un déplacement des cations cobalt du site I vers le site I' qui n’est pas
reproduit si une charpente déshydratée est utilisée, même en présence des molécules d’eau. Le
déplacement des cobalts du site II vers III' n’est possible qu’en présence de molécules d’eau
après prise en compte de la relaxation de la charpente.
Conclusion globale
135
transition ni donc, a fortiori, de la redistribution des cations dans ces structure lors de
l’adsorption/désorption des molécules d’eau. Ce travail pourrait être un point de départ pour
une étude plus systématique du comportement des cations lors d’adsorption/désorption des
molécules d’eau.
Le but initial de cette étude était de former un composite à transition de spin par
complexation de cations Co2+ par 3 ligands bipyridine dans la supercage de la zéolithe.
136
suivantes :
• La taille des cristaux utilisés (~ 180 µm selon [100]) est trop grande : la
complexation est efficace à proximité de la surface et diminue rapidement vers
le centre du cristal. On peut imaginer que les supercages près de la surface
contiennent des complexes ce qui rend difficile la migration des ligands vers le
centre du cristal. Il serait intéressent de tenter la formation du composite à
partir de poudre Na,Co-X.
137
138
Chapitre IV
139
140
IV Etude du système (dmpNA-mordénite) par Diffraction des rayons X sur poudre......... 143
IV.1 Introduction .......................................................................................................... 143
IV.2 Structure de la mordénite : Généralités............................................................. 145
IV.3 Structure de la mordénite par DRX sur monocristal ....................................... 147
IV.3.1 Conditions expérimentales. ............................................................................ 148
IV.3.2 Structure de la mordénite déshydratée ........................................................... 149
IV.3.3 Discussion ...................................................................................................... 152
IV.4 Composite (dmpNA-mordénite) : Diffraction des rayons X sur poudre......... 155
IV.4.1 Conditions expérimentales ............................................................................. 155
IV.4.2 Spectre à 90 K ................................................................................................ 155
IV.4.2.1 Indexation - Paramètres du profil............................................................... 155
IV.4.2.2 Affinement Rietveld................................................................................... 156
IV.4.3 Spectre à 300 K .............................................................................................. 162
IV.4.4 Discussion. ..................................................................................................... 165
IV.4.4.1 Distribution des cations (90K). .................................................................. 165
IV.4.4.2 Interactions dmpNA-Charpente ................................................................. 166
IV.3.4.3 Interactions intermoléculaires .................................................................... 169
IV.4.5 Diffractogramme de mordénites synthétisées avec et sans alcool. ................ 171
IV.5 CONCLUSION..................................................................................................... 174
141
142
IV Etude du système (dmpNA-mordénite) par
Diffraction des rayons X sur poudre.
IV.1 Introduction
143
Figure 4.1 : Alignement des molécules hyperpolarisables dans les canaux d’une zéolithe. Les canaux
zéolithiques permettent de contrôler le tilt θ et d'isoler les molécules les unes des autres.
Une des premières études sur les composites (molécule-zéolithe) a été réalisée par
Hoffmann et al. (1996) sur le système pNA-AlPO4-5 où la matrice aluminophosphatée
AlPO4-5 (AFI) était elle-même non-centrosymétrique. Malgré un remplissage faible et
inhomogène des canaux de l'AlPO4-5 par les molécules de pNA, le rendement SHG dépasse
celui de KDP. D’autres études ont suivi, impliquant des matrices centrosymétriques, cette
fois, telles que MFI (système pNA-ZSM-5, Reck et al., 1996) ou MTW (système pNA-ZSM-
12, Kinski et al., 2002). Le système pNA-ZSM-5 présente un fort rendement SHG mais ce
dernier se dégrade au cours du temps (effet attribué à une photolyse ou un désalignement des
molécules de pNA). L’efficacité de pNA-ZSM-12 est très faible, probablement à cause d’un
arrangement centrosymétrique des molécules dans les canaux (d’après les études structurales
et la RMN 1H et 13
C, 4 molécules occupent le canal de diamètre 5,6×5,9 Å2 où elles ne
peuvent pas bouger).
Dans cette thèse, le choix s'est porté sur le composite (dmpNA-mordénite) car la
mordénite possède des canaux droits larges de diamètre 6,5×7,0 Å2 supérieur à la taille des
principales molécules hyperpolarisables, dont les dérivés simples de pNA, et qu'elle peut être
synthétisée sous forme de monocristaux de grande taille et bonne qualité cristalline (Chapitre
II). La molécule de dmpNA (dimethylparaNitroAniline), comme son parent pNA, est
fortement hyperpolarisable : les moments dipolaires dans l'état fondamental et excité sont
5,7(2)/13,8(7) D et 5,5(2)/13,7(5) D pour dmpNA et pNA respectivement (Sinha et al., 1990).
La longueur de dmpNA (7,9 Å) n'est que légèrement supérieure à la périodicité de la
mordénite le long des canaux droits (~7,5 Å), et peut être compensée par de légères
désorientations des molécules dans les canaux (critère de commensurabilité).
144
Ce chapitre est dédié à l'étude par diffraction X sur monocristal de la structure de la
mordénite déshydratée et à une première détermination de la structure du système dmpNA-
mordénite par DRX sur poudre.
Charpente
La topologie de la zéolithe mordénite est décrite par la structure de type MOR dans
l'« Atlas of zeolite framework types » (Baerlocher et al., 2001). La charpente est composée de
tétraèdres TO4 (T=Si, Al) qui s’enchaînent pour former deux types de canaux continus
parallèles à [001] : les grands canaux 12MR (6,5×7,0 Å2) à section presque circulaire et les
petits canaux elliptiques 8MR (5,7×2,6 Å2). Ces canaux sont interconnectés par des canaux
elliptiques 8MR (3,4×4,8 Å2) parallèles à [010] (Figure 4.2).
145
rotule « inhabituel » T2-O8-T2 de 180°. Dans le zéotype MOR, on suppose une occupation
totale des sites T par Si et les distances T-O sont de 1,61 Å, les angles O-T-O de 109,47° et T-
O-T de 145°.
146
Cations
Selon la littérature, les cations compensateurs de charge sont distribués sur des sites
cristallographiques numérotés I à VI et représentés sur la Figure 4.4.
Figure 4.4 : Distribution des cations compensateurs dans la mordénite ; vue perpendiculaire à l’axe c.
Pour le groupe conventionnel Cmcm, le site I est au centre du 8MR [001] (Wyckoff
4b), le site II à l’intersection des canaux 8MR selon [001] et [010] (Wyckoff 4c), le site III est
dans la poche 5MR (Wyckoff 4c), le site IV se trouve à l’intersection de 12MR [001] et 8MR
selon [010] (Wyckoff 4c, 8f) et le site VI se trouve dans le canal principal 12 MR (Wyckoff
8g). Le site V a été abandonné.
La structure de la mordénite a été résolue pour la première fois par Meier (1961) dans
le groupe Cmcm à partir d'une zéolithe naturelle (Si/Al=5, cations Na+, K+, Ca2+) échangée
avec des cations sodium. Notre étude cristallographique a pour but de vérifier la distribution
cationique dans un cristal de ratio Si/Al=7, précurseur de ceux utilisés pour les études sur le
système (dmpNA-mordénite)
147
IV.3.1 Conditions expérimentales.
Les reflexions ont été collectées avec une redondance >5 jusqu'à une résolution de
~0,77 Å-1. Les pics de Bragg ont été intégrés en utilisant le programme CrysAlisPro (version
1.171.29.2) puis moyennées en mmm avec le logiciel Sortav (Blessing, 1995, 1997). La
structure a été résolue avec le programme SHELXS-97 (Sheldrick, 2008) puis affinée avec
Jana2000 (Petřίček, 2005). Les conditions expérimentales et les facteurs d’accord statistique
finaux sont assemblés dans le Tableau 4.1.
148
IV.3.2 Structure de la mordénite déshydratée
Les paramètres structuraux de la mordénite sont donnés dans le Tableau 4.2, les
distances et les angles de valence dans les Tableaux 4.3 et 4.4.
∗
L’affinement sans restriction des occupations des sites I et II mène à des paramètres occupation/Uiso trop
élevés 64(1)%/0,055(3) Å2 et 49(1) %/0,075(4) Å2 pour Na1 et Na2 respectivement avec des facteurs
d’accord statistiques finaux de R=5,76 % et Rw=10,82%. Ces valeurs élevées nous suggèrent un probable
désordre des cations Na+. .
149
Tableau 4.2 : Paramètres structuraux de la mordénite
Atome x y z Wyckoff Occ. (%) U11 ou Uiso U22 U33 U12 U13 U23
Si1 0,30326(7) 0,07289(5) 0,0428(2) 16 100 0,0210(6) 0,0159(6) 0,0138(5) 0,0021(6) 0,0045(6) 0,0016(5)
Si2 0,30313(7) 0,30957(6) 0,0453(1) 16 100 0,0219(7) 0,0179(6) 0,0131(6) 0,0056(7) 0,0004(6) -0,0001(5)
Si3 0,0868(1) 0,38175(8) 0,250000 8 100 0,0121(1) 0,0161(9) 0,0186(9) -0,0018(8) 0,00000 0,00000
Si4 0,0863(1) 0,22566(9) 0,250000 8 100 0,0149(9) 0,0146(9) 0,0205(9) 0,002(8) 0,00000 0,00000
O1 0,1222(2) 0,4149(2) 0,4277(5) 16 100 0,029(2) 0,063(2) 0,022(2) -0,013(2) -0,010(1) 0,001(2)
O2 0,1234(2) 0,1939(2) 0,4268(4) 16 100 0,030(2) 0,046(2) 0,026(2) -0,003(2) -0,008(2) 0,013(2)
O3 0,2386(3) 0,1229(2) -0,0112(5) 16 100 0,055(3) 0,033(2) 0,049(3) 0,020(2) -0,003(2) 0,008(2)
O4 0,0963(3) 0,3034(2) 0,250000 8 100 0,033(3) 0,015(2) 0,090(4) -0,001(2) 0,00000 0,00000
O5 0,1707(3) 0,1940(3) 0,750000 8 100 0,042(4) 0,087(5) 0,016(3) 0,028(4) 0,00000 0,00000
O6 0,1729(3) 0,4200(3) 0,750000 8 100 0,042(4) 0,038(3) 0,012(2) -0,007(3) 0,00000 0,00000
O7 0,2732(3) 0,00000 0,00000 8 100 0,027(3) 0,020(2) 0,049(3) 0,00000 0,00000 -0,005(2)
O8 0,25000 0,25000 0,50000 8 100 0,040(4) 0,046(4) 0,085(5) -0,026(3) -0,018(4) 0,039(4)
O9 0,00000 0,4030(4) 0,25000 4 100 0,016(4) 0,028(4) 0,049(5) 0,00000 0,00000 0,00000
O10 0,00000 0,2055(3) 0,25000 4 100 0,014(4) 0,029(5) 0,057(5) 0,00000 0,00000 0,00000
Na1 0,50000 0,00000 0,00000 4 51,2(8) 0,045(5) 0,045(5) 0,018(4) 0,00000 0,00000 0,016(4)
Na2 0,50000 0,0660(5) 0,25000 4 48,8(8) 0,043(3)
Na4a 0,00000 0,145(2) 0,562(4) 8 10,6(6) 0,06
Na4b 0,00000 0,187(3) 0,75000 4 7(1) 0,06
Na6 0,175(2) 0,038(1) -0,25000 8 8,6(8) 0,021(9)
150
Tableau 4.3 : Distances (Å) et angles (°) dans la mordénite
Si-O
Si1-O1 1,617(3) Si2-O2 1,598(3) Si3-O1 1,623(4) Si4-O2 1,618(3)
Si1-O3 1,601(4) Si2-O3 1,590(3) Si3-O4 1,605(4) Si4-O4 1,594(4)
Si1-O6 1,617(2) Si2-O5 1,606(2) Si3-O9 1,626(3) Si4-O10 1,612(2)
Si1-O7 1,614(2) Si2-O8 1,584(1)
<Si1-O> 1,612(3) <Si2-O> 1,595(2) <Si3-O> 1,619(4) <Si4-O> 1,611(3)
O-Si-O
O1-Si1-O3 111,9(2) O2-Si2-O3 110,2(2) O1-Si3-O1 110,1(2) O2-Si4-O2 109,7(2)
O1-Si1-O6 106,0(2) O2-Si2-O5 106,7(2) O1-Si3-O4 112,0(2) O2-Si4-O4 110,4(2)
O1-Si1-O7 108,3(2) O2-Si2-O8 110,2(2) O1-Si3-O9 105,6(2) O2-Si4-O10 107,6(2)
O3-Si1-O6 112,3(7) O3-Si2-O5 109,3(3) O4-Si4-O10 111,1(3)
O3-Si1-O7 106,8(2) O3-Si2-O8 110,0(2)
O6-Si1-O7 111,4(2) O5-Si2-O8 110,3(2)
<O-Si-O> = 109,5(2)
Si-O-Si
Si1-O1-Si3 145,8(3) Si4-O4-Si3 167,7(4) Si1-O7-Si1 140,7(4) Si4-O10-Si4 150,4(5)
Si4-O2-Si2 145,2(3) Si2-O5-Si2 145,4(4) Si2-O8-Si2 180
Si1-O3-Si2 159,6(3) Si1-O6-Si1 147,2(4) Si3-O9-Si3 149,5(5)
8MR 12MR
Interaction Coordination Interaction Coordination
Na1-O1 4 × 2,859(3) Na4a-O2 2 × 2,64(2)
Na1-O9 2 × 2,723(4) Na4a-O10 2 × 2,64(3)
Na4b-O5 2 × 3,088(7)
Na2-O6 2 × 3,133(5) Na4b-O2 4 × 3,296(5)
Na2-O1 4 × 3,294(3) Na6-O7 2 × 2,70(2)
Na6-O3 2 × 2,74(2)
Na6-O5 1 × 3,17(3)
151
IV.3.3 Discussion
Charpente
⎧ Al (T 1) = 1,76 ⋅ x ⎧ Al (T 1) = 10,56%
⎪ Al (T 2) = 0,98 ⋅ x ⎪ Al (T 2) = 5,88%
⎪ ⎪
⎨ soit, pour x=6 ⎨
⎪ Al (T 3) = 4,20 ⋅ x ⎪ Al (T 3) = 25,2%
⎪⎩ Al (T 4) = 2,83 ⋅ x ⎪⎩ Al (T 4) = 16,98%
La tendance globale indique que le site T3 est le plus riche en aluminium et le site T2
est le moins riche. La différence entre les deux relations réside dans l’occupation de T1 et T4.
Ces conclusions sont en parfait accord avec des calculs DFT effectués récemment par Oumi et
al. (2007) mais ces derniers ont montré que la stabilité des aluminiums en sites 1 et 4
152
dépendait fortement des conditions de calcul ce qui rendait difficile à tirer des conclusions.
Sur les ~6 cations Na+ par maille, 4 se trouvent dans les canaux elliptiques 8MR :
2,05(3) cations sur le site I (Na1) (4×Na1-O1 à 2,859(3) Å et 2×Na1–O9 à 2,723(4) Å) , et
1,95(3) cations sur le site II (Na2) (2×Na2–O6 à 3,133(5) Å, 4×Na2–O1 à 3,294(3) Å) (Figure
4.5). Ces sites sont les plus favorables et donc occupés à 100%.
a) b)
Figure 4.5 : Vues ORTEP de l’environnement des cations sodium Na1 (a) et Na2 (b) perpendiculaire aux
directions [001] et [010] respectivement.
Environ 1,8 Na+ se trouvent dans le grand canal 12MR : Na4a et Na4b sont à
l’intersection du canal 12MR et du 8MR parallèle à [010]. Na4a est coordonné aux oxygènes
de la charpente 2×O2 et 1×O10 avec des courtes distances 2,64(3) Å (Figure 4.6a), Na4b est
coordonné à la charpente avec des distances plus longues (2×Na4b–O5 de 3,088(7) Å et
4×Na4b–O2 de 3,296(5) Å) (Figure 4.6b). Na6 (dans le canal 12 MR est coordonné à 5
oxygènes de la charpente : 2×Na6–O7 (2,70(2) Å), 2×Na6–O3 (2,74(2) Å) et 1×Na6–O5
(3,17(3) Å) (Figure 4.7).
153
a) b)
Figure 4.6 : Environnement de Na4a (a) et Na4b (b) ; vue perpendiculaire à l’axe b.
D’après les calculs Monte-Carlo effectués par Maurin et al. (2001) les sites
principalement occupés sont les sites I (les plus favorables) puis IV et VI. Pour Si/Al=7, les
calculs donnent ~3,7 Na+ sur les sites I+II et ~2,3 sur les sites IV+VI comparés à 4 et 1,8
respectivement d'après notre étude expérimentale.
154
IV.4 Composite (dmpNA-mordénite) : Diffraction des
rayons X sur poudre.
Dans un premier temps, nous avons tenté de localiser des molécules dmpNA (p-N,N-
dimethylnitroaniline) par les études de diffraction des rayons X sur un monocristal
(180µm×100µm×70µm) en conditions de laboratoire. Cette étude n’a pas abouti, ce qui peut
s’expliquer par un remplissage faible et inhomogène des canaux de mordénite. Pour cette
raison nous nous sommes tourné vers l’étude des échantillons de petite taille
(~40µm×30µm×20 µm) par DRX sur poudre.
Les mesures de diffraction X ont été réalisées sur la ligne BL02B2 du synchrotron
SPring-8 (Hyogo, Japon) en collaboration avec K. Kato et M. Takata (λ=1,0017(1) Å). Les
diagrammes de poudre (90K et 300K) ont été collectés selon la géométrie Debye-Scherrer en
utilisant un détecteur bidimensionnel « Plaque Image », dans la gamme de résolution
2°<2θ<75,1°. Les affinements structuraux ont été réalisés avec la suite de programmes
« FullProf_Suite» (Rodríguez-Carvajal, 1993).
IV.4.2 Spectre à 90 K
Q20 1
Figure de mérite de « De Wolff » : M 20 = avec Qhkl = 2 , Q20 est la valeur de Q pour
2 Q ⋅ N 20 d hkl
la 20ième réflexion observée et indexée, N20 est le nombre de réflexions calculées (observées et inobservées)
jusqu’à Q20 , Q écart moyen des Q (observés par rapport aux calculés) pour les 20 premières réflexions.
155
F(20)=311,5∗∗ (Smith & Snyder, 1979). Les extinctions dues au mode de réseau C et au miroir
c[010] sont vérifiées (hkl : h+l=2n et h0l : l=2n).
Pour l'étape de « Profile Matching » (Le Bail et al., 1988) préliminaire aux
affinements Rietveld, une fonction Pseudo-Voigt de type "Thompson-Cox-Hasting"
(Thompson et al., 1987) incluant une correction d’asymétrie des pics jusqu’à 2θ=15° (Finger
et al., 1994) est utilisée. Les facteurs d’accord statistique valent Rp=1,70%, Rwp=2,53%
(Figure 4.8).
1 N
∗∗ FN = ⋅
∆ 2θ N possible , N est le nombre de lignes observées, Npossible est le nombre de lignes de diffraction
possibles jusqu’au N, ∆2θ est l’écart moyen entre les positions 2θ observée et calculée. FN permet
d’estimer la qualité de mesure (en estimant la position des lignes de diffraction) et la complétude des
données.
156
donc fait en deux étapes: la première concerne l’affinement de la charpente et des cations à
grands angles (20°<2θ<74,5° qui contient l’information sur la partie la plus ordonnée) ; la
deuxième concerne l’affinement complet (2,5°<2θ<74,5°) (une fois les molécules de dmpNA
localisées).
Figure 4.9: Affinement de la charpente (2θ>20°); Yobs en rouge, Ycalc en noir, Yobs-Ycalc en bleu et les positions
des pics de Bragg en vert
Les cations sodium en site I (Na1 (1/2 0 0)) et II (Na2 (1/2 0,0722(6) 1/4), localisés à
partir des cartes de Fourier-différence, occupent les mêmes positions que dans la structure
déshydratée, et ne migrent pas lors de l’adsorption de la molécule dmpNA. L’affinement libre
des paramètres occupations/Biso des cations converge vers les valeurs 66,12(8)%/5,05(6) Å2 et
61,6(1)%/7,25(9) Å2 pour Na1 et Na2 respectivement (Rp=16,3 %, Rwp=14,5 %, RBragg=10,9
%). Dans cette configuration, environ ¼ des sites I et II voisins est occupé simultanément ce
qui est improbable étant donnée la répulsion électrostatique entre ces cations (Maurin et al.,
2001). Nous avons donc testé un modèle où les occupations conjointes des sites I et II sont
contraintes à ne pas excéder 100 %. L’affinement converge vers des occupations/Biso de
54,48(4)%/2,84(4) Å2 et 45,52(4)%/4,26(6) Å2 pour Na1 et Na2 avec des facteurs d’accord
157
statistique légèrement meilleurs (Rp=16,0 %, Rwp=14,4 %, RBragg=10,8 %).
Avec cette hypothèse, 4 des 6 cations compensateurs de charge sont localisés dans les
petits canaux de la mordénite, et le reste doit donc se trouver dans le grand canal 12 MR. Les
cartes de Fourier différence ne révèlent qu'un pic de faible densité sur le site IV en position
spéciale 4c ne dépassant guère le bruit de fond. Toutes les tentatives de modélisation par Na+
ont échouées et 2 cations par maille restent donc non localisés à cette étape. Le désaccord
entre les spectres expérimental et calculé à bas angle (Figure 4.10) est principalement dû à la
contribution de la molécule de dmpNA non encore modélisée.
Figure 4.10 : Résultat de l’affinement (2θ≥20°) du squelette avec les cations compensateurs de charge.
Avec cette contrainte de corps rigide, seuls 6 degrés de liberté restent à ajuster en
recuit simulé : la position de l’origine (X Y Z) du système de référence de la molécule et les
158
trois angles d’Euler Θ, Φ et χ entre les axes (XmYmZm) de la molécule et les axes
(XcrisalYcristalZcristal) d'un repère orthonormé cristallographique (Figure 4.11). Avec ce choix,
l’angle Θ ([-π, π]) représente la mesure directe de l’inclinaison de la molécule par rapport à
l’axe du grand canal 12 MR de la mordénite. Φ et χ ([0, 2π]) sont des angles de rotation
autour de Zcristal et Zm respectivement.
a) b)
159
Figure 4.12 : Le profil des spectres après l’introduction du dmpNA à partir des simulations (sans affinements).
Les affinements successifs des paramètres de la molécule (Θ, Φ, χ, (X, Y, Z), Uiso et
occupation), des paramètres structuraux des cations et de la charpente, aboutissent à des
facteurs d’accord statistiques Rp=11,5/11,4 %, Rwp=10/10,2 % et RBragg=6,96/6,7 %, pour les
deux solutions A/B respectivement (Tableau 4.5). Les affinements avec la stratégie de
restriction d’occupation des sites I et II par Na+, adoptée lors de l'étape précédente aboutissant
à des valeurs déraisonnables des facteurs de Debye-Waller des cations (Uiso(Na1)<0),
l'affinement a finalement été fait sans restriction. L'occupation de dmpNA converge vers
12,066(7) % proche de la valeur maximale de 12,5 % avec un paramètre d’agitation
thermique, Biso= 3,42(7) Å2 modéré et similaire à celui de Na1 (2,91(5) Å2). Le faible
désaccord entre spectres observé et calculé à bas angles provient d'une correction imparfaite
de l’asymétrie des pics (Figure 4.13).
160
Figure 4.13 : Affinements Rietveld du spectre à 90 K; Yobs en rouge, Ycalc en noir, Yobs-Ycalc en bleu et les
positions des pics de Bragg en vert
La structure finale montre que les canaux principaux de la mordénite sont quasiment
remplis par ~2 molécules de dmpNA par maille désordonnées sur 8 sites cristallographiques
différents. La composition cristallographique du composite est
Na6(O2NC6H4N(CH3)2)1,93Si48O96. Les molécules de dmpNA sont inclinées de 28,66(3)° par
rapport à l’axe c ([001]) de la mordénite ; cette conformation est favorable pour poursuivre les
études optiques du composite et tester la Génération de Seconde Harmonique (Porcher et al.,
2008).
161
IV.4.3 Spectre à 300 K
162
Figure 4.15 Affinement Rietveld du diffratogramme à 300 K ; Yobs en rouge, Ycalc en noir, Yobs-Ycalc en bleu et
positions des pics de Bragg en vert.
Collecte 90 K
Rp (%) Rwp (%) χ2 RBragg (%) RF (%)
1) «Pattern Matching» (2,5°<2θ<74,5°) 1,70 2,54
Affinement à grands angles(20°<2θ<74,5°)
2) Charpente (Si, O) 20,0 18,6 39,9 13,8 11,4
3) Charpente (Si, O) + Na1, Na2 16,3 14,5 24,2 10,9 9,99
3 bis) + occ.(Na1)+occ.(Na2)≤100% 16,0 14,4 23,9 10,8 9,79
Introduction de dmpNA localisée par « Simulated Annealing » (2,5°<2θ<74,5°)
4A) Affinements dmpNA, Si, O, Na 11,5 10 62,2 6,96 9,47
Inversion NO2/N(CH3)2
4B) Affinements dmpNA, Si, O, Na 11,4 10,2 63 6,7 9,76
Collecte à 300K
1) Profil Matching (2,5°<2θ<74,5) 1,91 2,84
3) Charpente (Si, O) + Na1, Na2 (20°<2θ<74,5) 16,4 14,8 25 11,2 10,9
4) Charpente + dmpNA, Na (2,5°<2θ<74,5°) 11,4 10,5 48,9 7,4 11,1
163
Tableau 4.6: Paramètres structuraux finaux de dmpNA-mordénite (90 K)
164
Tableau 4.7 : Distances (Å) et angles (°) de la charpente de mordénite dans dmpNA-mordénite (90K)
Si-O
Si1-O1 1,588(6) Si2-O2 1,570(6) Si3-O1 1,649(5) Si4-O2 1,627(6)
Si1-O3 1,624(6) Si2-O3 1,557(7) Si3-O4 1,565(9) Si4-O4 1,651(9)
Si1-O6 1,614(3) Si2-O5 1,613(3) Si3-O9 1,620(5) Si4-O10 1,648(5)
Si1-O7 1,620(4) Si2-O8 1,596(3)
<Si1-O> 1,612(5) <Si2-O> 1,584(5) <Si3-O> 1,621(6) <Si4-O> 1,638(7)
O-Si-O
O1-Si1-O3 113,7(6) O2-Si2-O3 110,6(6) O1-Si3-O1 112,2(5) O2-Si4-O2 112,0(5)
O1-Si1-O6 106,8(4) O2-Si2-O5 105,6(4) O1-Si3-O4 109,7(6) O2-Si4-O4 108,6(7)
O1-Si1-O7 109,8(6) O2-Si2-O8 108,5(4) O1-Si3-O9 106,7(4) O2-Si4-O10 106,5(4)
O3-Si1-O6 109,7(5) O3-Si2-O5 114,0(6) O4-Si4-O10 114,7(8)
O3-Si1-O7 107,3(5) O3-Si2-O8 109,7(5)
O6-Si1-O7 109,5(4) O5-Si2-O8 108,1(4)
<O-Si-O> = 109,5(5)
Si-O-Si
Si1-O1-Si3 144,9(4) Si4-O4-Si3 163,3(6) Si1-O7-Si1 139,6(3) Si4-O10-Si4 146,8(4)
Si4-O2-Si2 146,2(4) Si2-O5-Si2 143,0(5) Si2-O8-Si2 180
Si1-O3-Si2 160,2(5) Si1-O6-Si1 151,6(3) Si3-O9-Si3 153,5(4)
IV.4.4 Discussion.
Environ 5,2 Na+ occupent les deux sites voisins I et II dans les petits canaux 8 MR de
la mordénite. La proximité des sites I et II (2,385(8) Å) interdit a priori une occupation
globale dépassant 4 cations par maille à cause de la répulsion électrostatique. La contrainte
d'occupation aboutissant à des paramètres structuraux déraisonnable (Biso(Na1)<0), la solution
non contrainte a due être adoptée au final (Tableau 4.6), même si elle conduit à un excès
apparent de cations sur ces sites (5,2 Na+). Cette anomalie traduit sans doute un désordre de
Na2 (facteur de Debye-Waller important, distances de coordination anormalement élevées :
2×Na2-O6 (3,184(9) Å), 4×Na2-O1 (3,266(5) Å) voire une occupation partielle de ce site par
des molécules d’eau (Tableaux 4.6, 4.8, Figure 4.16).
165
a) b) c)
Figure 4.16 : Position de Na2: a) vue perpendiculaire au canal 8MR selon [001], b) vue perpendiculaire à l’axe a,
c) visualisation des bouches de pores des canaux 8MR selon [010]).
Contrairement à Na2, Na1 est bien localisé et occupe à 49,44(3) % (Biso=2,91 Å2) le
site I où il est stabilisé par six interactions avec les oxygènes de la charpente 2×Na1-O9
(2,761(9) Å) et 4×Na1-O1 (2,864(6) Å) (Figure 4.17).
166
Figure 4.18 : Distribution désordonnée des molécules de dmpNA dans les canaux 12 MR ; vue perpendiculaire à
l’axe c.
a) b) c)
Figure 4.19 : Environnement local d'une molécules de dmpNA; a) plan (a b) ; b) plan (a c) ; c) plan (b c).
167
d(H…A)<R(H)+R(A)-0,12 Å et angle (D-H…A)>100° avec R(H)=1,10 Å, R(O)=1,35 Å et
R(C)=1,55 Å (Steiner, 2000). En utilisant ces définitions, la molécule et la zéolithe forment
des liaisons hydrogène impliquant les hydrogènes H3, H4 du cycle aromatique et les
oxygènes O10 et O2 de la charpente, et les hydrogènes H5 et H9 d’un des groupements
méthyles avec les oxygènes O5 et O2 (orientation A Figure 4.20, Tableau 4.9). Dans le cas de
l’orientation inverse de dmpNA (inversion des terminaisons NO2 et N(CH3)2), la molécule
interagit à travers les hydrogènes H2, H3 et H4 du cycle aromatique avec les oxygènes O7,
O10 et O2 de la charpente, et H9 d’un des groupements méthyles avec l’oxygène O3 (Tableau
4.9).
Tableau 4.9 : Contacts entre la charpente de la mordénite et dpmNA (d(D…A)<3,4 Å, d(H…A)<2,33 Å, angle
(D-H…A)>100°) pour les deux orientations NO2/N(CH3)2
168
pas suffisamment fortes pour éviter le désordre dans le canal, elles sont, cependant, assez
importantes pour préserver le composite.
Les molécules de dmpNA sont proches les unes des autres dans les canaux de la
mordénite et peuvent interagir par des liaisons «hydrogène». A cause du miroir m[001] (s’il
existe), ces molécules polaires peuvent adopter deux orientations correspondant à l'échange
des groupes terminaux nitro- ou diméthyleamino-. Localement, il est possible d’avoir des
interactions attractives NO2…N(CH3)2 ou répulsives NO2…NO2/ N(CH3)2…N(CH3)2. Les
interactions attractives entre le donneur N(CH3)2 et l’accepteur NO2 stabilisent l’arrangement
des molécules hôtes dans les canaux. Cet arrangement augmente leur polarisabilité ainsi que
la densité électronique sur l’azote de groupement nitro- à cause de l’effet donneur du
groupement alkyle. Cet effet peut-être responsable du déplacement du pic de réflectance du
composite dmpNA-mordénite (~305 nm) par rapport à ceux de dmpNA cristallin (~365 nm)
ou au minimum d’absorbance dans l’éthanol (~285 nm) (Figure 4.21).
100 0,7
90 mordénite
0,6 dmpNA-mordénite
80
dmpNA (solide)
70 0,5
Absorbance, u.a.
Réflectance, %
60
0,4
50
0,3
40
30 0,2
20
0,1
10
dmpNA (éthanol)
0 0,0
200 300 400 500 600 700 800
Longueur d'onde, nm
Figure 4.21: Spectres de reflectance totale de la mordénite, du composite dmpNA-mordénite et des cristaux de
dmpNA; spectre d’absorbance de dmpNA dans l’éthanol; spectrophotomètre UV-visible VARIAN Cary 4000
avec la sphère d'intégration.
169
Le large pic de réflectance du composite dmpNA-mordénite (~290-320 nm, transitions
n-π*/π-π* non résolues) (Figure 4.21) reflète une diversité des interactions et de
l’environnement des molécules désordonnées dans les canaux de la mordénite comme il a été
également remarqué pour d’autres composites (Hoffmann et al., 1996 ; Schulz-Ekloff et al.,
2002). La dernière hypothèse est confortée par le fait que le minimum de réflectance de
dmpNA-mordénite est plus proche du pic d’absorbance des molécules dmpNA dans l’éthanol
(~395 nm) qu’au minimum de réflectance dans l’état cristallin (~465 nm). Notons que les
minimum de réflectance du composite étudié dmpNA-mordénite et du dmpNA-AlPO4-5
(présentant un fort effet SHG) (Hoffmann et al., 1996) se situent à ~415 nm ce qui permet de
supposer un environnement similaire des molécules dmpNA.
Dans les canaux, la formation de la chaîne dipolaire est stabilisée par des interactions
électrostatiques de type «liaison hydrogène» entre molécules. Si on compare le composite
dmpNA-mordénite étudié avec les deux composites analogues pour lesquels un effet SHG a
été observé pNA-AlPO4-5 (AFI) et pNA-ZSM-5 (MFI), on retrouve des liaisons hydrogène
molécule-molécule pour pNA-AlPO4-5 tandis que pour pNA-ZSM-5 aucune liaison
hydrogène (d(H…O) minimale de 2,62 A) n'est mise en évidence (Reck et al., 1996). Pour
dmpNA-mordénite et pNA-AlPO4-5 les liaisons hydrogène stabilisent les chaînes dipolaires,
170
alors que pour pNA-ZSM-5 les interactions électrostatiques jouent ce rôle, ce qui pourrait
expliquer la faible résistance de ce dernier à l’irradiation laser.
L’alcool introduit dans la synthèse des cristaux de mordénite conditionne leur taille,
forme et qualité (Chapitre II). Les spectres de poudre des cristaux ainsi synthétisés sont très
différents (Figure 4.22).
Figure 4.22 : Diagramme de poudre des cristaux synthétisés avec (rouge) et sans (bleu) alcool.
171
Tout comme celui synthétisé avec alcool (Chapitres IV.3.2, IV.3.3), le spectre de la
mordénite sans alcool peut être indexé dans une maille orthorhombique (figures de mérite
M(20)=21 et F(20)=47,3) : a=18,0661(5) Å b=20,3798(5) Å, c=7,5033(2) Å. Le mode de
réseau C et miroir c sont vérifiés et les deux pics à 2θ=8,2° (Figure 4.23) et 11,54° non
indexés sont probablement dûs à des impuretés. En revanche, certains pics, (ex. 111), ont une
contribution Lorentzienne nettement plus importante que le reste du diagramme (Figure 4.23).
Figure 4.23 : Elargissement du pic 111, mal modélisé avec un modèle isotrope; Yobs en rouge, Ycalc en noire.
Dans leurs études structurales de Na+- et H+-mordénites par DRX sur poudre, Rudolf
& Garcés (1994) ont mis en évidence (suggéré) l'existence de défauts de croissance où
environ 21, 12 et 9 % de la structure (pour des rapports Si/Al=6 ; 7,5 ; 9,5 respectivement),
sont décalés de c/2 par rapport au reste. Il faut noter que le taux de défauts est plus important
que les 3 % évoqués par Mortier (1975), Simoncic & Armbruster (2004). D’après Rudolf &
Garcés (1994) les défauts de croissance n’induisent pas d’élargissement des raies mais
affectent les hauteurs relatives de certaines raies comme par exemple (111/130), (241/002),
(350/201), (511/530). Effectivement, lorsqu’on compare les diagrammes de poudre des
mordénites synthétisées sans et avec alcool (Figure 4.24), l’intensité relative des raies 110,
020, 200 varie peu tandis que pour les raies (111/130), (241/002), (511/530) la tendance est
inversée : dans la mordénite synthétisée avec alcool les hauteurs des raies 111, 241, 511 sont
plus hautes que celles de 130, 002, 530 alors que pour les mordénites synthétisées sans alcool
les raies 130, 002, 530 ont quasiment la même intensité que 111 ou sont supérieures aux raies
241et 511 (Figure 4.24).
172
Figure 4.24 : Diagrammes de poudre des cristaux synthétisés avec (en rouge) et sans (en bleu) alcool. La
variation relative des raies est liée aux défauts structuraux éventuelles (Rudolf & Garcés, 1994).
173
IV.5 CONCLUSION
L’adsorption des molécules de dmpNA ne modifie pas l’occupation des sites les plus
favorables I, II par des Na+. Les paramètres structuraux montrent que les sodiums restent bien
stabilisés dans le site I, en revanche, ils sont désordonnés autour des sites II ou il y a des
molécules d’eau. Finalement, nous avons localisé seulement 4 cations compensateurs de
charge sur 6 attendus.
Les grands canaux de la mordénite sont occupés par les molécules de dmpNA
désordonnées sur 8 sites cristallographiques différents dans le canal 12 MR, et inclinées de
28,66(3)° par rapport à l’axe c. La structure cristallographique couplée à l’analyse chimique
donne la composition chimique suivante : Na6(O2NC6H4N(CH3)2)1,93Al6Si42O96, indiquant un
remplissage quasi total des canaux par dmpNA. Les chaînes dipolaires de dmpNA sont
stabilisées par des liaisons hydrogène entre la charpente et les molécules et entre molécules.
L’empilement des molécules est favorable à la Génération de Seconde Harmonique.
Les cristaux de mordénite synthétisés sans alcool ont visiblement des propriétés
d’adsorption meilleures (pour les molécules étudiées) que ceux synthétisés avec alcool
(Chapitre II) ce qui est également conforme aux études d’adsorption de Oumi et al. (2003).
L’étude du spectre de poudre sur ces cristaux met en évidence la présence de la phase
amorphe dans l’échantillon. Les études de la présence de défauts de croissance évoquée par
plusieurs auteurs (Mortier, 1975 ; Rudolf & Garcés, 1994 ; Simoncic & Armbruster, 2004) est
en cours mais si on se base sur les études de Rudolf & Garcés, 1994 (par comparaison des
spectres) les cristaux synthétisés avec alcool contiennent moins de défauts que ceux
synthétisés sans alcool.
174
Chapitre V
« Conclusions et Perspectives »
175
176
V Conclusions et Perspectives
Conclusions
Le but de cette thèse était l'élaboration et l’étude de composites basés sur la zéolithe X
et la mordénite.
Même si le projet initial d’élaborer le composite magnétique n’a pas abouti, notre
étude cristallographique de la zéolithe Na,Co-X a apporté des informations intéressantes et
nouvelles. En particulier, nous avons montré que la géométrie de la charpente et la
distribution (migration) des cations compensateurs de charge sont très sensibles à l’état
d’hydratation de la zéolithe Na,Co-X. Ce point est peu documenté dans la littérature, en
particulier en ce qui concerne les études cristallographiques sur les zéolithes bicationiques.
Nous avons mis en évidence la complémentarité entre les calculs Monte-Carlo et les études
expérimentales par diffraction des rayons X. Lorsque l’expérience ne permet pas d'identifier
avec certitude la nature chimique du cation occupant un site (désordre chimique ou
positionnel), la simulation permet d’aiguiller le choix. Réciproquement, les potentiels utilisés
en Monte-Carlo ne permettant pas actuellement d'optimiser à la fois la distribution cationique
et de relaxer simultanément la charpente, ces calculs bénéficient des informations
expérimentales sur ce dernier point. Nous avons pu confirmer ainsi la stabilisation des cations
177
sodium dans le prisme hexagonal, et expliquer en partie, dans le cas de Co2+ (et plus
généralement dans celui des métaux de transition de la 1ère série) pourquoi les échanges
cationiques dans les zéolithes Na-X n'étaient que partiels. Les fortes contraintes structurales
imposées à la charpente lors de l'échange 2 Na+ → Co2+ expliquent également l'amorphisation
des cristaux de zéolithe X lors des échanges quasi complets «en colonne». En complément,
les études locales par EXAFS au seuil du Co ont mise en évidence la déformation du prisme
hexagonal en confirmant la supposition (sous réserve) de Verberckmoes et al. (1997, 1998).
Cette étude peut servir de point de départ pour comprendre le comportement structural des
zéolithes X échangées avec autres métaux de transition de 1ère série.
178
des similarités entre le nouveau composite dmpNA-mordénite et le pNA- AlPO4-5 qui a
montré un fort effet SHG (la stabilisation des chaînes dipolaires par des interactions type
« liaison hydrogène », des interactions électrostatiques et liaisons hydrogène entre la
charpente-molécules, est suffisamment forte pour préserver le composite).
Perspectives
Zéolithe X
A plus long terme, il serait utile de développer un outil d’affinement conjoint DRX et
MC appliqué aux structures poreuses avec des invités différents. Comme il est montré dans la
thèse, le choix des espèces, difficile à définir expérimentalement, serait aiguillé par les
simulations utilisant les informations expérimentales les plus sûres. En retour on pourrait
optimiser les charges et les paramètres de potentiels d’interaction pour un système particulier
pour ensuite mettre en évidence les phénomènes physiques mis en jeu et responsables des
propriétés des matériaux.
Mordénite
179
peut essayer d’aligner macroscopiquement les cristaux de (dmpNA- mordénite) ou regarder la
possibilité de former des films de mordénites orientées. Les premières tentatives d’alignement
sur une grille de monocristaux à l'aide d'un champ électrique se sont révélées infructueuses,
mais ce point devrait être approfondi (alignement mécanique à l'aide d’un film polymère ou
par un champ électrique externe).
Afin d’augmenter le signal optique (SHG) on peut essayer de réduire l'angle de tilt des
molécules dans les canaux (initialement à ~28°). Il est probable que la présence, même
réduite, de cations compensateurs de charge en site IV et VI dans le grand canal influence la
position des molécules invitées. Un échange Na+ → Ca2+ (même rayon ionique) pourrait
permettre de libérer les canaux droits par une redistribution des Ca2+ sur les sites I, II en
laissant le grand canal libre. On peut également imaginer réduire l'angle de tilt en ajustant
davantage la taille de la molécule hyperpolarisable à celle des canaux (ou vice versa).
180
«Annexes»
181
182
ANNEXES
ANNEXE A : Rappels sur la Diffraction des Rayons X
La correction de Lorentz permet de rectifier l’intensité des réflexions due au fait que
l’intensité mesurée dépend du temps qu’un nœud reste sur la surface de la sphère d’Ewald. La
correction de polarisation tient compte de la polarisation des faisceaux incident et diffracté.
La correction d’absorption tient compte du chemin parcouru par les faisceaux incident et
diffracté dans le cristal aboutissant à une différence entre réflexions équivalentes par
symétrie. En pratique, les corrections d’anisotropie de la forme du cristal sont réalisées soit
par indexation des faces (Coppens, 1970 ; DeTitta, 1985) pour les études des mordénites, soit
par détermination de surface de transmission d’un cristal pour les études des cristaux dans les
capillaires (programme Sortav, Blessing, 1995).
Puis les intensités corrigées, Iobs(hkl), sont moyennées dans le groupe ponctuel du
cristal. L’expérience de diffraction des rayons X sur monocristal permet de mesurer des
intensités équivalentes par symétrie avec une certaine redondance pour un groupe ponctuel
donné. Après la remise à la même échelle des intensités collectées dans des conditions
expérimentales différentes (temps d’exposition, variation de l’intensité du faisceau incident),
la moyenne des réflexions équivalentes permet d’obtenir les valeurs précises des intensités
diffractées et leurs incertitudes. Pour moyenner les données, le programme Sortav (Blessing,
*
Hypothèse de théorie cinématique de diffraction : les ondes incidentes sont progressives à l’intérieur du
cristal, le cristal est infini (étude de diffraction à l’infini), l’énergie du faisceau incident est conservée, l’onde
diffusée ne rediffuse pas à l’intérieur du cristal.
183
1997) est utilisé.
∑w I i i
I = i =1
N (éq. A.2)
∑w
i =1
i
1
ou w1 = σ 2 ( I ) = Sigma éq. A.5
i
2
⎛ ⎡ ⎛ Z ⎞
2 ⎤⎞
I − Mediane( I i )
w2 = ⎜⎜1 − min ⎢1, ⎜⎜ ⎟⎟ ⎥ ⎟ avec Zmax=4 et Z = i
⎟ (éq. A.6)
⎜ ⎢⎣ ⎝ Z max ⎠ ⎥⎦ ⎟ Sigma
⎝ ⎠
⎧ sin θ ⎫
σ ( I ) = max ⎨ ESD, Q( I , ) × ESD, RMSD⎬ éq. A.7
⎩ λ ⎭
avec ESD (« Estimated Standard Deviation ») et RMSD (« Root Mean square Standard
Deviation ») :
N N
∑ wσ i
2
(Ii )
N ∑ w (I i i − I )2
ESD = (éq. A.8), RMSD =
i =1 i =1
(éq. A.9)
N
N −1 N
∑ wi
i =1
∑ wii =1
Pour juger la qualité des données, les facteurs statistiques suivants sont déterminés :
184
2
N N N N ⎛ I i− I ⎞
∑
r ∑ Ii − I
N −1 i
∑r ∑ wi ⎜⎜
H N −1 i ⎝ σ (Ii )
⎟⎟
⎠
Rint = H
N (éq. A.10) wRint = 2 (éq. A.11)
∑∑ I i
r
N ⎛ I ⎞
∑r ∑ wi ⎜⎜ i ⎟⎟
⎝ σ (Ii ) ⎠
H i
H i
Facteur de structure
Le Facteur de structure est une fonction complexe définie pour T=0 K par
r N r rr
F ( H ) = ∑ f k ( H ) exp(2πiHrk ) (éq. A.12)
k =1
r r r r
où N est le nombre d’atomes dans la maille, rk = xk a + yk b + zk c est un vecteur du réseau
r r r r
direct représentant la position de l’atome k, H = ha * + kb * + lc * est un vecteur du réseau
r
réciproque et f k ( H ) est le facteur de diffusion atomique de l’atome k.
r
Le facteur de diffusion atomique, f k ( H ) caractérise le pouvoir de diffusion
cohérente des rayons X par un atome k, il est obtenu par Transformée de Fourier de sa densité
k r
électronique ρ (r ) :
r r rr r
f k ( H ) = ∫ ρ k (r ) exp(2πiHr )dr (éq. A.13).
sin θ
Le facteur de diffusion atomique décroît en fonction de ( ), à θ=0 il vaut le nombre
λ
d’électrons (Z) de l’atome. Les atomes de numéro atomique Z voisin offrent un faible
contraste (Si/Al, Na+/H2O, ...) aux rayons X.
Les atomes dans la maille sont localisés par leurs positions moyennes (x, y z) et les
ellipsoïdes de déplacement atomique (traduisant l’agitation thermique des atomes ou leur
désordre) donnés par des déplacements temporels moyens isotropes Uiso ou anisotropes Uij.
L’amplitude de vibration des atomes du cristal augmente avec la température en diminuant la
cohérence entre les ondes diffractées par les différents atomes et en affaiblissant l’intensité (et
185
le facteur de structure) de l’onde diffractée. L’agitation thermique de chaque atome est
r N r r r
F ( H ) = ∑ f k ( H ) exp(2πiHrk )Wk ( H ) (éq. A.14)
k =1
r
[ r
] ⎡ sin 2 θ ⎤
Wk ( H ) = exp− 2π 2 uk2 H 2 = exp− ⎢ Bk
λ2 ⎦
⎥ (éq. A.15)
⎣
avec Uij=<uiuj>, i, j=1, 2, 3 sont des coefficients de tenseur de déplacement atomique repérés
par rapport aux axes cristallographiques du cristal.
Lors des études de zéolithes il est fréquent qu’un site ne soit occupé que partiellement
par un atome. Ainsi, le facteur de diffusion de l’atome est pondéré par son coefficient
d’occupation (occ.k) :
r N r r r
F ( H ) = ∑ f k ( H ) exp(2πiHrk )occ.k Wk ( H ) (éq. A.17)
k =1
D’après les équations A.15 et A.17, il est évident que, lorsque l’agitation thermique est
affinée en isotrope, une corrélation entre les paramètres occupation/Uiso complique
l’affinement structural.
Les zéolithes étudiées lors de cette thèse (zéolithe X et mordénite) ne permettent pas
186
l’étude de la densité électronique car la qualité des données est fortement limitée par le
désordre lié d’une part à la distribution des atomes Si/Al formant la charpente, et d’autre part
à la difficulté de différencier et localiser les espèces extracharpente.
Les intensités expérimentales donnent accès uniquement aux modules des facteurs de
structure ; la résolution consiste donc à déterminer leurs phases. Au cours de cette thèse, des
modèles des squelettes zéolithiques existant déjà (Porcher, 1998 ; Meier, 1961) sont utilisés
en apportant l’estimation suffisante des phases. Le travail consistait surtout à affiner le
modèle en minimisant par moindres carrés l’écart entre les modules de facteurs de structure
r r
calculés du modèle Fcalc(H) et les facteurs de structures expérimentaux Fobs(H) . La fonction
minimisée est
r
( (r)
χ 2 = ∑r w( H ) Fobs
2
H − kFcalc
2
( ))
H
r 2
éq. A.18
H
r 1
où la fonction de pondération est w( H ) =
( )
r r 2 éq. A.19
σ 2 ( Fobs
2
( H )) + 0,005 × Fobs
2
(H )
r r
Les paramètres ajustés sont le facteur d’échelle (k) entre I calc ( H ) et I obs ( H ) , les positions
des atomes (x y z), leurs facteurs d’agitation thermique (Uiso/Uij) et l’occupation des sites.
187
La progression et la qualité des affinements sont estimées par les facteurs d’accord
statistiques R, Rw et Gof (« Goodness of fit »):
r r
∑ w( H )(F ( H ) − kF )
r r r 2
∑F
r
2
obs( H ) − kFcalc
2
(H )
r
2
obs
2
calc (H )
R= H r A.20 Rw = H r
∑r Fobs2 ( H )
éq. éq.
H
∑ w( H ) F
r
4
obs
H
∑ w( H )(F )
r r r 2
r
2
obs ( H ) − kFcalc
2
(H )
A.21 Gof = H éq. A.22
N obs − N par
La qualité de l’affinement est d’autant meilleure que R et Rw sont les plus faibles et
Gof est le plus proche de 1. Les cartes Fourier-différence TF ( Fcalc − Fobs ) permettent
d’introduire dans le modèle les nouveaux atomes et de continuer l'affinement. Pour localiser
et modéliser les espèces extracharpente, sont utilisées des données expérimentales ayant
I/σ(I)>2,5-3 et dans une gamme de résolution (sinθ/λ) limitée.
188
Méthode de Rietveld.
Introduction
Fonction de profil
Le profil du spectre expérimental {2θi, Yiobs} est modélisé par le spectre calculé {2θi,
Yicalc} où la fonction Yicalc représente la somme de toutes les contributions des intensités de
Bragg au point 2θi (plus le bruit du fond) :
( )
r
Yi calc = ∑r [ S ⋅ L ⋅ A ⋅ P ⋅ F 2 H ⋅ Ω(2θ )] + bi (éq. A23)
H
L : correction de Lorentz-polarisation
A : coefficient d’absorption,
189
Ω(2θ) : fonction de profil,
Dans notre cas, le bruit du fond {2θi, bi} est déterminé par interpolation linéaire entre
~80 points échantillonnés sur le diffractogramme.
La fonction de profil Ω(2θ) choisie est la fonction de type Pseudo-Voigt proposée par
Thompson, Cox, et Hastings (Thompson et al., 1987):
2 ln 2 ⎛ 4 ln 2 ⎞
G= ⋅ exp⎜⎜ − (2θi − 2θ k ) 2 ⎟⎟ éq. A.26
Hk π ⎝ Hk
2
⎠
2 1
L= ⋅
πH k ⎛ 4(2θi − 2θ k ) 2 ⎞ éq. A27
⎜1 + ⎟
⎜ H k2 ⎟
⎝ ⎠
Pour cette fonction, les largeurs à mi-hauteur (FWHM) Hk des fonctions Gaussienne (HG) et
Lorentzienne (HL) sont contraintes, ce qui permet de relier, les largeurs des raies
expérimentales aux paramètres microstructuraux de contraintes X, U, DST ou de tailles de
cristallites Ig, Y, F(SZ), lorsque les paramètres de Cagliotti U, V, W décrivant la résolution
instrumentale sont connus:
IG
H G2 = (U + DST
2
)tg 2θ k + Vtgθ k + W +
cos 2 θ k éq. A28
⎛ Y + F (SZ ) ⎞
H L = Xtgθ k + ⎜⎜ ⎟⎟ éq. A29
⎝ cos θ k ⎠
Le paramètre de mélange η
2 3
HL ⎛H ⎞ ⎛H ⎞
η = 1,36603 − 0,47719⎜ L ⎟ + 0,11116⎜ L ⎟ éq. A30
H ⎝ H ⎠ ⎝ H ⎠
( )
H = H G5 + 2,69269 H G4 H L+2,42843H G3 H L2 + 4,47163H G2 H L3 + 0,07842 H G H L4 + H L5 éq. A31
Selon que l'on fasse un «Fit de Le Bail» préliminaire à l'extraction des intensités I(hkl)
190
pour une résolution structurale, ou un affinement structural par la méthode de Rietveld, le
mode de calcul du facteur de structure F(H) change :
Fit de LeBail : Ne disposant pas de modèle structural (sauf les paramètres de maille et
le groupe d'espace), les F(H) sont affinés librement pour tous les pics de Bragg de diagramme,
sans autre contrainte que la forme de la fonction de profil à respecter. En raison du
recouvrement entre pics adjacents, une intensité commune peut être affectée à un cluster de
réflexions adjacentes, sans redistribution interne de l'intensité.
Affinement Rietveld : Les valeurs des facteurs de structure F(H) sont calculées à
partir du modèle structural, comme dans le cas des monocristaux, puis sont corrigées par un
facteur tenant compte de la multiplicité de la réflexion.
i =1
( )
χ 2 = ∑ wi Yi obs − Yi calc (α ) éq. A.32
2
Facteurs d’accord
La qualité des affinements est contrôlée par des facteurs d’accord statistiques
suivants :
Pattern Matching :
N
∑ Yi
obs
− Yi calc
i =1
Rp = N
∑ Yi
obs
i =1 éq. A.33
( ) ⎤⎥
1/ 2
⎡N 2
⎢ i∑ wi ⋅ Yi obs − Yi calc
=1
Rwp =⎢ ⎥
( )
N éq. A.34
obs 2
⎢ ∑ wi ⋅ Yi ⎥
⎣ i =1 ⎦
Affinement Rietveld :
191
(Yi obs −bi )
1/ 2
⎡ ⎤
⎢ N −P ⎥
=⎢N
2⎥
Rexp
( )
éq. A.35
⎢ ∑ wi ⋅ Yi obs ⎥
⎣ i =1 ⎦
Rwp
Goodness of fit (Gof) S = Rexp éq. A.36
2
⎛ Rwp ⎞
χν = ⎜2 ⎟ = S 2 éq. A.37
⎜ Rexp ⎟
⎝ ⎠
Bragg (RBragg) :
r r
∑r I Hr
"obs"
( H ) − I Hcalc
r (H )
RB = H
éq. A.38
∑r I "Hr "
obs
H
RF = H
éq. A.39
∑r F "Hr "
obs
H
r
r I "obs" ( H )
F "obs"
(H ) = r
L( H )
192
ANNEXE B : Approche théorique de l’absorption des
rayons X
0,7
XANES
0,6
0,5 EXAFS
0,4
µ(E)
0,3
0,2
0,1
0,0
7400 7600 7800 8000 8200 8400 8600 8800
E(eV)
Figure B1 : Spectre d’absorption des rayons X (XAS) au seuil K de cobalt (Eseuil = 7709 eV) réalisé sur
l’échantillon Na,Co-X déshydraté
La partie XANES est due aux transitions électroniques de l’atome central lors de
l’absorption des rayons X d’énergie inférieure à l’énergie de liaison des électrons en couche
profonde. Elle renseigne principalement sur l’état d’oxydation de l'atome absorbeur et sa
coordination (tétraédrique, octaédrique …, type de ligand).
L'EXAFS est une mesure des variations de la probabilité d’absorption des rayons X
d’un atome, dues à son état physico-chimique, en particulier son état d’oxydation et sa
coordination (distances, nombre et type des ligands) ; elle permet ainsi d’accéder à
l’arrangement structural à courte distance de l’atome absorbeur. Cette technique est
complémentaire de la diffraction, qui apporte des informations à longue portée.
Les points forts de l'EXAFS, par rapport aux autres techniques, sont les suivants :
• l’ordre à longue distance n’est pas nécessaire (étude de matériaux cristallisés ou non),
193
• les mesures sont peu exigeantes quant à la préparation de l'échantillon (cependant les
mesures se réalisent dans la plupart des cas sur synchrotron).
Mécanisme d’absorption.
La base de la théorie EXAFS est l’effet photoélectrique (Figure B2): l'absorption d'un
photon X d’énergie supérieure à l'énergie de liaison d’un électron en couche profonde (Ephoton
X=Eliaison+Ecinétique) excite l'atome, entraînant l'éjection de l’électron vers les couches externes
ou vers le continuum en laissant derrière lui un « trou ». Le photoélectron est considéré
comme une onde sphérique (Figure B3) de nombre d’ondes k=2π/λ.
photoélectron
photon X
Atome absorbeur
Rj
Figure B3 : Propagation de l’onde associée au photoélectron ; onde initiale (ligne pleine), onde rétrodiffusée
(ligne pointillée).
Cette onde est alors diffusée par les atomes voisins situés à la distance Rj de l’atome
194
i 2 kR j
e
central et prend la forme χ j (k ) ∝ f j (k ) ; le déphasage 2kRj dépend de la distance entre
R 2j
l’atome central et les diffuseurs environnants. La propagation de cette onde dans les potentiels
des atomes central et diffuseurs modifie le déphasage d'une quantité (δj(k)-π/2), et son
π
i ( 2 kR j + δ j ( k ) − )
e 2
expression devient χ j (k ) ∝ f j (k ) . Les interférences constructives ou
R 2j
sin(2kR j + δ (k ))
χ j (k ) ∝ f j (k ) .
R 2j
L'analyse des spectres EXAFS se fait en séparant les contributions des atomes voisins en j
sphères de coordination concentriques regroupant les Nj atomes se trouvant à des distances
sin(2kR j + δ (k ))
voisines de Rj : χ (k ) = ∑ χ j (k ) ∝ ∑ N j f j (k ) .
j j R 2j
Les transitions électroniques qui sont à l’origine du spectre XANES sont plus
complexes à mettre en équation et à modéliser, et l’interprétation quantitative des spectres
XANES est beaucoup moins développée aujourd’hui que celle d’EXAFS. Seule une
interprétation qualitative du spectre XANES est donnée dans ce travail.
2 2Rj
N j S 0 (k ) F j (k ) −
χ (k ) = ∑ λ ( k ) −2 k σ j
2 2
2
sin[2kR j + δ j (k )]e e éq. B.1
j kR j
195
2m( E − E0 )
avec k = où E0 est l’énergie de liaison du photoélectron (origine en énergie).
h2
L'affinement EXAFS consiste à dériver les variables de cette équation par affinement
par moindres carrés du spectre expérimental χexp(k). Certaines variables sont peu affectées par
la structure locale autour de l'absorbeur et peuvent être paramétrisées à partir de composés de
référence :
relaxent vers un état p f . Typiquement le facteur S02 a des valeurs comprises entre
0,7 et 0,8. Les variations de S02 avec l’environnement étant peu importantes, S02 est
transférable entre atomes de même type, dans le même état d’oxydation. S02(k), ne
varie pratiquement pas avec k et peut être approximé par une fonction constante (Stern
dans Koningsberger, ed. Prins, Chapitre 1).
m
• Fj(k) ( F j (k ) = f (k ) ): Amplitude de diffusion effective (« effective scattering
2πh 2 k
amplitude »). C’est une fonction calculée pour chaque type d’atome, qui modélise le
potentiel de diffusion des atomes libres. F(k) varie avec l’énergie ; sa décroissance est
régulière pour les éléments ayant Z<50, plus complexe pour les éléments lourds
(Z>50) (Stern dans Koningsberger, ed. Prins, Chapitre 1).
196
constants. δ(k) décroît quasi linéairement avec k, la non-linéarité à bas k se traduisant
par un léger élargissement des pics (dans l'espace R).
• λ(k) : Durée de vie de l’état excité. La durée de vie de l’état excité est limitée par la
durée de vie du « trou » et du photoélectron lui-même et est typiquement de quelques
femtosecondes. Cette durée de vie est prise en compte par un terme de libre parcours
−2 R j
moyen e λ ( k ) représentant la probabilité pour que le photoélectron rétrodiffusé sur un
atome voisin retourne au centre de l’atome central avant que d’autres diffusions ne se
produisent ou que le « trou » ne soit rempli.
D'un point de vue pratique, les fonctions F(k), δ(k) et λ(k) sont calculées par le programme
FEFF inclus dans ARTEMIS (Ravel & Newville, 2005).
Les autres variables de l'équation EXAFS sont directement liées à la structure locale
de l'échantillon :
197
Expérience sur synchrotron : mesure au seuil du cobalt.
Figure B.4: Dispositif expérimental de mesure de spectre EXAFS sur la ligne BL17C1 (NSRRC, Taïwan).
Mesure de fluorescence
198
Figure B5 : Schéma du processus d’absorption des rayons X (à gauche) avec la fluorescence (à droite).
Pour l'absorption au seuil K, les transitions Kα sont dominantes et sont utilisées sur la
ligne BL17C1. Cependant, le signal de fluorescence contient également le rayonnement des
autres raies caractéristiques du Co ainsi que les rayonnements caractéristiques des autres
éléments de l'échantillon, en plus des rayonnements dus aux diffusions élastiques (effet
Thomson) et inélastiques (effet Compton). Ces rayonnements parasites, contribuent au bruit
de fond. La diffusion élastique étant fortement anisotrope à cause de la polarisation du
rayonnement synchrotron, cet effet est atténué en plaçant le détecteur à 90° par rapport au
faisceau incident. Les autres raies de fluorescence, d'énergie supérieure à la raie Kα du Co,
sont éliminées par un filtre formé d'une fine feuille de l'élément de numéro (Z-1) par rapport à
l'élément étudié, i.e. le Fe (Eseuil=7112 eV).
L'analyse du spectre EXAFS est réalisée en trois étapes principales (selon la procédure
décrite par Sayers, Bunker dans Koningsberger, ed. Prins, Chapitre 6) : Premièrement, le
spectre d'absorption expérimental mesuré, µ(E), est mis sous la forme que l'on peut analyser,
χ(k), en modélisant le bruit de fond et en normalisant le spectre. Deuxièmement, on crée un
modèle structural à partir des informations issues de la Diffraction des Rayons X d'où on
extrait des informations sur chaque sphère de coordination j possible autour de chaque atome
absorbeur : type d'atome diffuseur (F(k)), nombre de coordination Nj, distance Rj. On calcule
alors un spectre χ(k) modèle.
Troisièmement, les variables de l'équation d’EXAFS Rj, σj2, Nj, S02, ∆E0 sont affinées afin
que le spectre modèle approche au mieux le spectre expérimental; les paramètres F(k), λ(k) et
δ(k) sont fixes. La qualité d'affinement est estimée par les facteurs d'accord statistique, les
corrélations entre variables et, surtout, par la vraisemblance de la structure locale obtenue.
Ces étapes sont décrites en détail dans la partie qui suit.
199
Extraction du spectre χ(k) expérimental.
2m( E − E0 ) µ(E)−µ0(E)
χ(k)=f( χ(E)) avec k = et χ(E)=
h2 ∆µ0(E0) éq. B.2
χ(E) est déterminé à partir du spectre expérimental µ(E) après la soustraction de l’absorption
de l’atome central (isolé) µ0(E) et la normalisation par ∆µ0(E0) du saut d’absorption au seuil
E0 ("constante de normalisation"). La fonction finale χ(k) ne contient alors que des
informations sur l'environnement local de l'atome central car sa propre absorption est
éliminée.
µ0(E) représente l'absorption des niveaux profonds de l'atome central. Cette fonction
est difficile à modéliser car l'absorption de l'atome central varie avec l'énergie de manière
inconnue. La solution retenue est d'approximer µ0(E) empiriquement par une fonction
« spline » de manière à lisser le spectre expérimental µ(E). La fonction « spline » est
composée de fonctions polynomiales cubiques liées entre elles par des nœuds équidistants, ce
qui assure qu'une seule oscillation complète peut exister entre deux nœuds et permet de
contrôler la fréquence des oscillations exclues du spectre. Le nombre de nœuds N est variable
:
2∆kRbkg
N = 1+ éq. B.3, avec
π
Rbkg : Distance à laquelle on considère que l'absorption est due à l'atome central
Normalisation
200
première dérivée de µ(E), ce qui permet de calculer le pré-seuil, µ-(E), et le bruit de fond
après-seuil, µ0+(E), qu'on extrapole ensuite jusqu'au seuil E0.
Dans notre cas, µ-(E) est modélisée par une droite passant par deux points (E0+Epré-
seuil_1, E0+Epré-seuil_2 avec Epré-seuil_1, Epré-seuil_2 < 0), qu'on extrapole jusqu'à E0 dont l'ordonnée
-
fixe µ (E0). Dans la région post-seuil, la fonction de bruit de fond µ0(E) est affinée par une
fonction quadratique de l'énergie (f~E2) passant par deux points (E0+Eaprès-seuil_1, E0+Eaprsè-
seuil_2 avec Eaprès-seuil_1, Eaprsè-seuil_2 >0). Cette fonction quadratique est extrapolée au seuil E0
pour obtenir la valeur µ0+(E0). La normalisation ∆µ0(E0) est alors estimée comme
+ −
: ∆µ 0 ( E 0 ) = µ0 ( E 0 ) − µ ( E 0 ) éq. B4
La procédure d'affinement (logiciel Artemis (Ravel & Newville, 2005) est réalisée
dans l'espace R (χ(R)=TF-1(χ(k)) par la méthode d’isolation spectrale (affinement pic après
pic). Une fonction fenêtre est appliquée pour isoler un ou plusieurs pics (dans l'espace des R).
La troncature du spectre se traduit par des pics parasites qui apparaissent lorsque la
transformée de Fourier est appliquée et qui s'ajoutent à ceux induits par la résolution limitée
des données. Cet effet de troncature peut être vu comme la convolution d'un spectre de
résolution infinie avec une fonction fenêtre W(R)/W(k):
201
1 ∞ w
χ ( R) = ∫ k χ (k )W (k )e dk ou
i 2 kR
2π 0
1 ∞ − i 2 kR
χ (k ) = ∫ χ ( R )W ( R)e dR
2π 0
Pour les modélisations des spectres EXAFS, nous avons choisi d'utiliser des fenêtres
de type "Hanning" (ARTEMIS, Ravel & Newville, 2005):
⎧ 2 π (k − k min + dk / 2)
⎪sin ( ), k min − dk / 2 ≤ k < k min + dk / 2
⎪⎪ 2dk
W (k ) = ⎨1, k min + dk / 2 ≤ k ≤ k max − dk / 2 ,
⎪
⎪cos 2 ⎛⎜ π (k − k max + dk / 2) ⎞⎟, k max − dk / 2 < k ≤ k max + dk / 2
⎪⎩ ⎝ 2dk ⎠
avec dk=0,2 Å-1, zone de l’espace de résolution où la fenêtre n’est pas égale à 1.
Cette fenêtre en cosinus carré permet une atténuation rapide des oscillations parasites,
mais avec, en contrepartie, un élargissement des pics.
N
∑ ( fi )
2
2
χ 2
χ 2 = ind2 ∑ ( fi ) éq B.5, χν 2 =
N N i =1
éq. B.6, R = N éq. B.7, avec
Nε i =1 ν ∑ χ exp ( Ri )
2
( )
i =1
202
ε : incertitude sur la fonction f qui contient des erreurs systématiques venant de la
mesure ou introduites lors des traitement des données (cf. manuel de FEFFIT) et des
fluctuations aléatoires de la mesure (bruit).
203
204
ANNEXE C : Principe des simulations Monte-Carlo.
1 N i2 i2 i2
H = T +U = ∑ ( p x + p y + p z ) + U ( x ,..., z ) éq. C.1
1 N
2m i =1
⎡ U (r N ) ⎤ ⎡ U (r N ) ⎤
− ⎢−
N N
∫ dr A( r ) exp ⎢ ⎥ exp ⎥
A = ⎣ k B T ⎦ = ∫ dr A(r )
N N ⎣ k BT ⎦
= ∫ dr N A(r N ) ρ (r N )
⎡ U (r ) ⎤
N Z
∫ dr exp ⎢−
N
⎥
⎣ k BT ⎦
éq. C.2
⎡ U (r N ) ⎤
où Z = ∫ dr N
exp ⎢− ⎥ (éq. C.3) est une intégrale sur l’ensemble des configurations
⎣ k BT ⎦
accessibles au système, appelé également fonction de partition. ρ est la densité de probabilité
de trouver le système dans configuration r N à l’équilibre. En pratique, il est rarement possible
de calculer l’intégrale de configuration Z . On peut toutefois accéder numériquement à la
valeur moyenne A sans passer par le calcul de la fonction de partition Z à l’aide des
simulations Monte-Carlo développées pour la première fois par Metropolis (Metropolis et al.,
1953).
Algorithme de Metropolis.
moyenne A , en effectuant une marche dirigée dans l’espace des phases. Chaque
configuration va apparaître avec une probabilité égale à son poids de Boltzmann. Cette
méthode est appelée échantillonnage d’importance. Ainsi, on ne calcule pas l’intégrale (on ne
somme pas sur toutes les configurations) mais la somme sur n configurations choisies
205
où A(ri) correspond à la valeur de A dans la configuration ri.
Pour échantillonner l’espace des phases, on choisit une marche dirigée stochastique
markovienne où les propriétés sont énoncées ci-dessous :
• chaque état peut être atteint à partir d’un autre état quelconque (chaîne irréductible) en
un nombre fini d’itérations.
Le passage d’un état m vers n est réalisé en utilisant une matrice de transition Π stochastique,
ρ mπ m → n = ρ nπ n → m éq. C.5
D’où
206
acc(m → n) ρ (n) ⎛ U ( n ) − U ( m) ⎞
= = exp⎜⎜ − ⎟⎟ éq. C.8.
acc(n → m) ρ (m) ⎝ k BT ⎠
Déplacement atomique.
⎧ x n = x m + (2ξ − 1)δr
⎪
⎨ y n = y m + (2ξ − 1)δr éq. C.9
⎪ z = z + (2ξ − 1)δr
⎩ n m
étudiée ( A − A )2 .
207
on applique les mouvements de translation des cations. De plus, les sauts de cations de site en
site cristallographique sont utilisés afin d’exploiter au mieux l’espace de configuration (ce
biais permet d’accéder aux différents minima d’énergie sans franchir les barrières d’énergie).
Le problème des simulations Monte-Carlo est que le système peut rester piégé dans
une région de l’espace de configuration à cause du faible taux d’acceptation des tentatives de
mouvement, de création ou de destruction. Ainsi, pour augmenter le taux d’acceptation, un
nombre de biais est introduit dans le programme « Gibbs » afin de corriger l’expression de la
probabilité d’accepter un pas Monte-Carlo. Lors de ces simulations les biais suivants sont
utilisés :
Les détails sur les biais et les expressions de la probabilité d’acceptation dans la thèse
d’A. Di Lella (2007).
Potentiel Buckingham.
γ
VijBuckingham = α ij exp(− β ij rij ) −
rij6
208
au LCP.
⎡⎛ σ ⎞6 ⎛ σ ⎞12 ⎤
Vij = −4ε ij ⎢⎜ ij ⎟ − ⎜ ij ⎟ ⎥
⎢⎜⎝ rij ⎟⎠ ⎜⎝ rij ⎟⎠ ⎥
⎣ ⎦
Les coefficients d’interactions εij, σij sont obtenus en utilisant la règle de Lorentz-
Berthelot :
σ ii + σ jj
σ ij =
2
ε ij = ε iiε jj
Les coefficients εii, σii pour un nouvel atome (M) sont évalués à partir de la valeur de sa
polarisabilité (αM) et rayon ionique (rM) et des coefficients (εii, σii, α, r) d’un atome connu ; en
particulier, dans la thèse A. Di Lella à partir de Na (Dang, 1995) : εNa=50,27 K, σNa=2,584 Å,
αNa=1,002 bohr3, rNa=1,02 Å.
6
αM ⎛σ ⎞ rM − rNa
ε M = ε Na ⋅ ⎜⎜ Na ⎟⎟ et σ M = 6 + σ Na
α Na ⎝ σM ⎠ 2
Pour modéliser les interactions Eau-Eau, le modèle TIP4P proposé par Jorgensen et al.
(1983) pour l’eau en volume (bulk) est utilisé. Dans ce modèle, les interactions dispersion-
répulsion sont traitées par le potentiel type Lennard Jones où on considère qu’une molécule
d’eau possède un centre de dispersion-répulsion centré sur l’oxygène. Pour calculer l’énergie
électrostatique, la molécule d’eau est décrite par deux charges identiques placées sur les
hydrogènes et une troisième charge compensatrice placée sur le centre de masse de la
molécule.
209
Limites des simulations Monte-Carlo.
Le modèle énergétique utilisé pour décrire les interactions entre les particules repose
toujours sur des approximations plus ou moins importantes. La plupart du temps, les énergies
sont estimées soit de manière empirique soit par calcul quantique. Dans les deux cas les
énergies sont évaluées en conditions particulières ; ainsi, elles sont limitées pour décrire des
propriétés particulières (car elles sont validées sur certaines propriétés uniquement). D’autre
part, l’approximation des interactions à 2 corps est souvent utilisée (dans notre cas aussi),
tandis qu’il est souvent plus raisonnable de considérer des interactions à n corps. Un autre
problème est la taille du système considérée infinie dans les simulations.
210
«Bibliographie»
211
212
Bibliographie
Alberti A., Davoli P., Vezzalini G., The crystal structure refinement of a natural mordenite, Z.
Kristallogr. 175 (1986), 249-256
Anderson M. W., Agger J. R., Hanif N., Terasaki O., Growth models in microporous
materials, Micropor. Mesopor. Mater. 48 (2001), 1-9
Arivanandhan M., Sankaranarayanan K., Ramasamy P., Melt growth of novel organic
nonlinear optical material and its characterization, Mater. Lett. 61 (2007), 4836-4838
Aubert E., Etude des interactions de molecules hôtes dans des zéolithes synthétiques par
diffraction des rayons X à haute résolution, these, Université Henri Poincaré Nancy-1, 2003
Bae D., Seff K., Structures of cobalt (II)-exchanged zeolite X, Micropor. Mesopor. Mater. 33
(1999), 265-280
Bae D., Seff K., Crystal structure of zeolite X nickel(II)-exchanged at pH 4.3 and partially
dehydrated, Ni2(NiOH)35(Ni4AlO4)2(H3O)46Si101Al91O384, Micropor. Mesopor. Mater. 40
(2000), 219-232
Bae D., Seff K., Extensive intrazeolitic hydrolysis of Zn(II): partial structures of partially and
fully hydrated Zn(II)-exchanged zeolite X, Micropor. Mesopor. Mater. 40 (2000), 233-245
Bae M. N., Song M. K., Kim Y., Seff K., Crystal structure of Mn46Si100Al92O384 89H2S, a
hydrogen sulfide sorption complex of fully dehydrated Mn2+-exchanged zeolite X, Micropor.
Mesopor. Mater. 63 (2003), 21-31
213
Baerlocher Ch., Meier W.M., Olson D.H., Atlas of Zeolite Framework Types, ELSEVIER,
Amsterdam-London-New York-Oxford-Paris-Shannon-Tokyo, 2001
Barrer R. M., Synthesis and reactions of mordenite, J. Chem. Soc. (1948), 2152
Barrer R. M., Baynham J. W., Bultitude F. W., Meier W. M., J. Chem. Soc. (1959), 195
Barros M. A. S. D., Araújo Jr. I. F., Arroyo P. A., Sousa-Aguiar E. F., Tavares C. R. G.,
Multicomponent ion exchange isotherms in NaX zeolite, Latin American Applied Research 33
(2003), 339-344
Beauvais C., Guerrault X., Coudert F.-X., Boutin A., Fuchs A., Distribution of sodium cations
in faujasite-type zeolite : A canonical parallel tempering simulation study, J. Phys. Chem. B
108 (2004), 399-404
Blessing R. H., An empirical correction for absorption anisotropy, Acta Cryst. A 51 (1995),
33-38
Blessing R. H., Outlier treatment in data merging, J. Appl. Cryst. 30 (1997), 421-426
Blower C. J., Smith T. D., The gas-phase decomposition of nitromethane over metal ion-
exchanged sodium Y zeolite and sodium X zeolite, Zeolites 13 (1993), 394-398
Borissenko E., Porcher F., Bouché A., Lecomte C., Souhassou M, Single crystal structure of
fully dehydrated partially Co2+-exchanged zeolite X; comparison with partially dehydrated
partially Co2+-exchanged zeolites X, Micropor. Mesopor. Mater. 114 (2008), 155-165
Boultif A., Louër D., Indexing of powder diffraction patterns for low-symmetry lattices by the
successive dichotomy method, J. Appl. Cryst. 24 (1991), 987-993
Breck D. W., Zeolite Molecular Sieves, Wiley, New York (1974), 163
214
Brühwiler D., Calzaferri G., Molecular sieves as host materials for supramolecular
organization, Micropor. Mesopor. Mater. 72 (2004), 1-23
Burkett S. L., Davis M. E., Mechanism of structure direction in the synthesis of Si-ZSM-5 :
An investigation by intermolecular 1H-29Si CP MAS NMR, J. Phys. Chem. 98 (1994), 4467-
4653
Burton A., Elomari S., Chen C. Y., Harris T. V., Vittoratos E. S., SSZ-53 and SSZ-59 : Two
novel extra-large pore zeolites, Stud. Surf. Sci. Catal. 154 (2004), 126-132
Buttefey S., Boutin A., Mellot-Draznieks C., Fuchs A., A simple model for predicting the Na+
distribution in anhydrous NaY and NaX zeolites, J. Phys. Chem. B 105 (2001), 9569-9575
Campbell B. J., Cheetham A. K., Linear framework defects in zeolite mordenite, J. Phys.
Chem. B 106 (2002), 57-62
Campbell B. J., Welberry T. R., Broach R. W., Hawoong Hong, Cheetham A. K., Elucidation
of zeolite microstructure by synchrotron X-ray diffuse scattering, J. Appl. Cryst. 37 (2004),
187-192
Charnell J. F., Gel growth of large crystals of sodium A and sodium X zeolites, J. Cryst.
Growth 8 (1971), 291-294
Coppens P., Crystallographic Computing, ed. Ahmed F. R., Hall S. R., Huber C. P.,
Copenhagen : Munksgaard (1970), 255-270
Corval A., Trommsdorff H. P., Optical spectra and excited-state dynamics of cis-thioindigo,
J. Phys. Chem. 93 (1989), 7081-7087
215
Cruickshank D. W. J., Rollett J. S., Electron-density errors at special positions, Acta Cryst. 6
(1953), 705
CrysAlisPro Software System, Version 1.171.29.2, Oxford Diffraction Ltd, UK, 2006
Cundy C. S., Lowe B. M., Sinclair D. M., Crystallisation of zeolitic molecular sieves : direct
measurements of the growth behaviour of single crystals as a function of synthesis conditions,
Faraday Discuss. 95 (1993), 235-252
Cundy C. S., Cox P. A., The hydrothermal synthesis of zeolites : Precursors, intermediates
and reaction mechanism, Micropor. Mesopor. Mater. 82 (2005), 1-78
Dang L. X., Mechanism and thermodynamics of ion selectivity in aqueous solutions of 18-
crown-6 ether: a molecular dynamics study, J. Am. Chem. Soc. 117 (1995), 6954-6960
Davis M. E., Lobo R. F., Zeolite and molecular sieve synthesis, Chem. Mater. 4 (1992), 756-
768
De Bont P. W., Vissenberg M. J., San de Beer V. H. J., Rob van Veen J. A., van Santen R. A.,
van der Kraan A. M., An x-ray absorption fine structure study of the sulfidation behavior of
CoNaY prepared by ion exchange : the influence of physisorbed water and cobalt loading, J.
Phys. Chem. B 102 (1998), 5876-5885
Depmeier W., Tetragonal tetrahedra distortions in cubic sodalite frameworks, Acta Cryst. B
40 (1984), 185-191
Dessau R. M., Valyocsik E. W., Goeke N. H., Aluminium zoning in ZSM-5 as revealed by
selective silica removal, Zeolites 12 (1992), 776-779
216
DeWilde W., Peeters G., Lunsford J. H., Synthesis and spectroscopic properties of tris(2,2'-
bipyridine) ruthenium (II), J. Phys. Chem. 84 (1980), 2306-2310
De Wolff P. M., A simplified criterion for the reliability of a powder pattern indexing, J.
Appl. Cryst. 1 (1968), 108-113
Di Lella A., Méthodes de simulation moléculaire pour l’étude de la distribution des cations et
d’adsorption de molécules polaires dans les zéolithes, thèse, Université Paris-sud 11, 2007
Dooryhee E., Catlow C. R. A., Couves J. W., Maddox P. J., Thomas J. M., Greaves G. N.,
Steel A. T., Townsend R. P., A study of cation environment and movement during
dehydration and reduction of nickel-exchanged zeolite Y by x-ray absorption and diffraction,
J. Phys. Chem. 95 (1991), 4514-4521
Engelhardt G., Cation location in dehydrated zeolite NaY revisited : SI position is displaced
from the center of the hexagonal prism, Microporous Mater. 12 (1997), 369-373
Feddag A., Bengueddach A., Echange ionique dans les zéolites par la méthode sèche. Mise en
évidence de l’influence de quelques paramètres, J. Phys. IV France 124 (2005), 243-248
Feoktistova N. N., Zhdanov W. L., Büllow M., On the Kinetics of crystallisation of silicalite
I, Zeolites 9 (1989), 136-139
Finger L. W., Cox D. E., Jephcoat A. P., A correction for powder diffraction peak asymmetry
due to axial divergence, J. Appl. Cryst. 27 (1994), 892-900
Firor R. L., Seff K., Hydrolysis in zeolites. Five-coordinate Aluminium (III). Investigations of
two hydrated partially ion-exchanged zeolite A. Crystal structures, Ni(II)Na-A and Fe(II)Na-
A, J. Phys. Chem. 82 (1978), 1650-1655
217
Francis R. J., O’Hare D., The kinetics and mechanisms of the crystallisation of microporous
materials, J. Chem. Soc., Dalton Trans. (1998), 3133-3148
Frising T., Leflaive P., Extraframework cation distributions in X and Y faujasite zeolite: A
review, Micropor. Mesopor. Mater. 114 (2008), 27-63
Gfeller N., Megelski S., Calzaferri G., Transfer of electronic excitation energy between dye
molecules in the channels of zeolite L, J. Phys. Chem. B 102 (1998), 2433-2436
Gilbert J. E., Mosset A., Large crystals of mordenite and MFI zeolites, Mat. Res. Bull. 33
(1998), 997-1003
Hammoudi H., Bendenia S., Marouf-Khelifa K., Marouf R., Schott J., Khelifa A., Effect of
the binary and ternary exchanged on crystallinity and textural properties of X zeolites,
Micropor. Mesopor. Mater. 113 (2008), 343-351
Heidler R., Janssens G. O. A., Mortier W. J., Schoonheydt R. A., Charge sensitivity analysis
of intrinsic basicity of faujasite-type zeolites using the electronegativity equalization method
(EEM), J. Phys. Chem. 100 (1996), 19728-19734
Helfferich F., Ion exchange, New York San Francisco, Toronto, London (1962), Mc Grow-
Hill Book Company
Heo N. H., Cruz-Patalinghug W., Seff K., Crystal structure of zeolite 4A ion exchanged to the
limit of its stability with nickel (II), J. Phys. Chem. 90 (1986), 3931-3935
Hincapie B. O., Garces L. J., Zhang Q., Sacco A., Suib S.L., Synthesis of mordenite
nanocrystals, Micropor. Mesopor. Mater. 67 (2004), 19-26
218
Hoffmann K., Marlow F., Caro J., Optical characterization of organized adsorbates in zeolite
microcrystals: Polarized adsorption spectroscopy, Zeolites 16 (1996), 281-286
Hoppe R., Schulz-Ekloff G., Wöhrle D., Kirschhock C., Fuess H., Synthesis, location and
photoinduced transformation of zeolite-encaged thioindigo, Langmuir 10 (1994), 1517-1523
Iizuka T., Lunsford J. H., Active intermediates in the reduction of nitric oxide by ammonia
over a CoY zeolite, J. Am. Chem. Soc. 100 (1978), 6106-6110
Jaramillo E., Auerbach S. M., New force field for Na cations in faujasite-type zeolites, J.
Phys. Chem. B 103 (1999), 9589-9594
Jeong G. H., Kim Y., Seff K., Crystal structures of the NO and NO2 sorption complexes of
fully dehydrated fully Mn2+-exchanged zeolite X (FAU), Micropor. Mesopor. Mater. 93
(2006), 12-22
Jones J. B., Al-O and Si-O tetrahedral distances in aluminosilicate framework structures, Acta
Cryst. B 24 (1968), 355-358
Jorgensen W. L., Chandrasekhar J., Madura J. D., Impey R. W., Klein M. L., Comparison of
simple potentiel functions for simulating liquid water, J. Chem. Phys. 79 (1983), 926-935
Joshi U. D., Joshi P. N., Tamhankar S. S., Joshi V. V., Shiralkar V. P., Adsorption behavior of
N2, water, C6 hydrocarbons, and bulkier benzene derivative (TMB) on Na-X zeolite and its
K+-, Rb+-, and Cs+-exchanged analogues, J. Colloid Interf. Sci. 235 (2001), 135-143
Joshi U. D., Joshi P. N., Tamhankar S. S., Joshi V. P., Idage B. B., Joshi V. V., Shiralkar V.
P., Influence of the size of extraframework monovalent cations in X-type zeolite on their
thermal behavior, Thermochim. Acta 387 (2002), 121-130
219
K
Kacirek H., Lechert H., Investigations on the growth of the zeolite type NaY, J. Phys. Chem.
79 (1975), 1589
Kelly S. D., Ravel B., EXAFS energy shift and structural parameters, CP882, X-ray
Absorption Fine Structure-XAFS13, ed. Hedman B, Pianetta P. (2007), 132-134
Kim Y., Han Y. W., Seff K., Crystal structure of fully dehydrated fully Tl+-exchanged zeolite
X, Zeolites 18 (1997), 325-333
Kim M. J., Jeong M. S., Kim Y., Seff K., Crystal structures of fully dehydrated fully Sr2+-
exchanged zeolite X and of its ammonia sorption complex, Micropor. Mesopor. Mater. 30
(1999), 233-241
Kinski I., Daniels P., Deroche C., Marler B., Gies H., Structure and properties of the
composite zeolite silica-ZSM-12/para-nitroaniline, Micropor. Mesopor. Mater. 56 (2002), 11-
25
Kohno Y., Shibata Y., Oyaizu N., Yoda K., Shibata M., Matsushima R., Stabilization of
flavylium dye by incorporation into the pore of protonated zeolites, Micropor. Mesopor.
Mater. 114 (2008), 373-379
Koningsberger D. C., ed. Prins R., X-Ray Absorption : principles, applications, techniques of
EXAFS, SEXAFS and XANES, a Wiley-Interscience Publication John Wiley & Sons (1988)
Kossel W., Extenoling the law of Bravais, Nach Ges Wiss Gottingen (1927), 143
Kumar R., Mukherjee P., Pandey R. K., Rajmohanan P., Bhaumik A., Role of oxyanions as
promoter for enhancing nucleation and crystallisation in the synthesis of MFI-type
microporous materials, Micropor. Mesopor. Mater. 22 (1998), 23-31
220
L
Le Bail A., Duroy H., Fourquet J. L., Ab-initio structure determination of LiSbWO6 by X-ray
powder diffraction, Mat. Res. Bull. 23 (1988), 447-459
Le Van Mao R., Xiao S., Ramsaran A., Yao J., Selective removal of silicon from zeolite
frameworks using sodium carbonate, J. Mater. Chem. 4 (1994), 605-610
Lechert H., Staelin P., Kuntz C., Quantitative relations of the batch composition and the Si/Al
ratio in the product of zeolites, Zeolites 16 (1996), 149-156
Lee P. A., Pendry J. B., Theory of the extended X-ray absorption fine structure, Physical
Review B 11 (1975), 2795-2811
Leike I., Marlow F., One-dimensional phase transition of dipolar molecules inside zeolite
pores, Zeolites 16 (1996), 65-69
Lewis D. W., Freeman C. M., Catlow C. R. A., Predicting the templating ability of organic
additives for the synthesis of microporous materials, J. Phys. Chem. 99 (1995), 11194-11202
Li G., Chen L., Bao J., Li T., Mei F., A recoverable catalyst Co(salen) in zeolite Y for the
synthesis of methyl N-phenylcarbamate by oxidative carbonylation of aniline, Appl. Catal. A
346 (2008), 134-139
Li Y., Wang Y., Wang A., Message filtering characteristics of semiconductor laser as receiver
in optical chaos communication, Opt. Commun. 281 (2008), 2656-2662
Lozano-Castello D., Zhu W., Linares-Solano A., Kapteijn F., Moulijn J. A., Mircopore
accessibility of large mordenite crystals, Micropor. Mesopor. Mater. 92 (2006), 145-153
221
Löwenstein W., The distribution of aluminium in the tetrahedra of silicates and aluminates,
Am. Mineralogist 39 (1954), 92
Lu B., Oumi Y., Itabashi K., Sano T., Effect of ammonium salts on hydrothermal synthesis of
high-silica mordenite, Micropor. Mesopor. Mater. 81 (2005), 365-374
Maes A., Cremers A., Ion exchange of synthetic zeolite X and Y with Co2+, Ni2+, Cu2+ and
Zn2+ ions, J. Chem. Soc., Faraday Trans. 1 71(1975), 265-277
Marlow F., Wübbenhorst M., Caro J., Pyroelectric effects on molecular sieve crystals loaded
with dipole molecules, J. Phys. Chem. 98 (1994), 12315-12319
Maurin G., Senet P., Devautour S., Gaveau P., Henn F., Van Doren V. E., Giuntini J. C.,
29
Combining the Monte-Carlo technique with Si NMR spectroscopy : simulations of cation
locations in zeolites with various Si/Al ratios, J. Phys. Chem. B 105 (2001), 9157-9161
Megelski S., Lieb A., Pauchard M., Drechsler A., Glaus S., Debus C., Meixner A. J.,
Calzaferri G., Orientation of fluorescent dyes in the nano channels of zeolite L, J. Phys.
Chem. B 105 (2001), 25-35
Meier W. M., The crystal structure of mordénite (ptilolite), Z. Kristallogr. 115 (1961), 439-
450
Mintova S., Olson N. H., Valtchev V., Bein T., Mechanism of zeolite a nanocrystal growth
from colloids at room temperature, Science 283 (1999), 958-960
222
Mizuno K., Lunsford J. H., Electron paramagnetic resonance study of tris(2,2’-bipyridine)
cobalt (II) complexes in zeolite Y : evidence for a spin equilibrium, Inorg. Chem. 22 (1983),
3484-3486
Morris M., Sacco A., Dixon A. G., Thompson R. W., The role of an aluminium-tertiary
alkanolamine chelate in the synthesis of large crystal zeolite NaA, Zeolites 11 (1991), 178-
183
Mortier W. J., Pluth J. J., Smith J. V., Positions of cations and molecules in zeolites with the
mordenite-type framework; II Dehydrated hydrogen-ptilolite, Mat. Res. Bull. 10 (1975),
1319-1325
Nightingale E. R., Jr., Phenomenological theory of ion solvation. Effective radii of hydrated
ions, J. Phys. Chem. 63 (1959), 1391-1387
Olson D. H., A reinvestigation of the crystal structure of the zeolite hydrated NaX, J. Phys.
Chem. 74 (1970), 2758-2764
Oumi Y., Kakinaga Y., Kodaira T., Teranishi T., Sano T., Influences of aliphatic alcohols on
crystallization of large mordenite crystals and their sorption properties, J. Mater. Chem. 13
(2003), 181-185
Oumi Y., Kanai T., Lu B., Sano T., Structural and physico-chemical properties of high-silica
mordenite, Micropor. Mesopor. Mater. 101 (2007), 127-133
223
P
Palancher H., Hodeau J.-L, Pichon C., Bérar J.-F., Lynch J., Rebours B., Rodriguez-Carvajal
J., Direct localization of atoms in mixed-occupancy powder by resonant contrast diffraction,
Angew. Chem. Int. Ed. 44 (2005), 1725-1729
Parpot P., Teixeira C., Almeida A. M., Ribeiro C., Neves I. C., Fonseca A. M., Redox
properties of (1-(2-pyridylazo)-2-naphthol)copper(II) encapsulated in Y zeolite, Micropor.
Mesopor. Mater. (2008), article in presse
Pascual P., Ungerer P., Tavitian B., Pernot P., Boutin A., Development of a transferable guest
host force field for adsorption of hydrocarbons in zeolites. I. Reinvestigation of alkane
adsorption in silicalite by grand canonical Monte-Carlo simulation, Phys. Chem. Chem. Phys.
5 (2003), 3684-3693
Petřίček V., Dušek M., Crystallographic Computing System Jana 2000, Institute of Physics,
Academy of Science of the Czech Republic, 2005
Pauling L., The structure of sodalite and helvite, Z. Kristallogr. 74 (1930), 213-225
Porcher F., Souhassou M., Dusausoy Y., Lecomte C., The crystal structure of dehydrated
NaX zeolite, Eur. J. Mineral. 11 (1999), 333-343
Porcher F., Borissenko E., Souhassou M., Takata M., Kato K., Rodriguez-Carvajal J.,
Lecomte C., Synchrotron powder diffraction characterization of the zeolite-based (p-N,N-
dimethylnitroaniline-mordenite) guest-host phase, Acta Cryst. B64 (2008), 713-724
224
Q
Quayle W. H., Peeters G., De Roy G. L., Vansant E. F., Lunsford J. H., Synthesis and
spectroscopic properties of divalent and trivalent tris(2,2'-bipyridine)iron complexes in zeolite
Y, Inorg. Chem. 21 (1982), 2226-2231
Ramamurthy V., Sanderson D. R., Eaton D. F., Control of dye assembly within zeolites: role
of water, J. Am. Chem. Soc. 115 (1993), 10438-10439
Ravel B., Newville M., ATHENA, ARTEMIS, HEPHAESTUS : data analysis for X-ray
absorption spectroscopy using IFEFFIT, J. Synchrotron Rad. 12 (2005), 537-541
Reck G., Marlow F., Kornatowski J., Hill W., Caro J., Structure of dipole chains in an MFI
type molecular sieve, J. Phys. Chem. 100 (1996), 1698-1704
Reshak A. H., Auluck S., Kityk I. V., Linear and nonlinear optical susceptibilities for a novel
borate oxide BaBiBO4: Theory and experiment, J. Solid State Chem. 181 (2008), 789-795
Reyes-Esqueda J., Darracq B., García-Macedo J., Canva M., Blanchard-Desce M., Chaput F.,
Lahlil K., Boilot J. P., Brun A., Lévy Y., Effect of chromophore-chromophore electrostatic
interactions in the NLO response of functionalized organic-inorganic sol-gel materials, Opt.
Comm. 198 (2001), 207-215
Rietveld H. M., A profile refinement method for nuclear and magnetic structures, J. Appl.
Cryst. 2 (1969), 65-71
Robson H., ed. Lillerud K. P., Verified syntheses of zeolitic materials, XRD patterns,
ELSEVIER, Amsterdam-London-New York-Oxford-Paris-Shannon-Tokyo, 2001
225
Rodríguez-Carvajal, Recent advances in magnetic structure determination by neutron powder
diffraction, Physica B 192 (1993), 55-69
Rollmann L. D., Schlenker J. L., Kennedy C. L., Kennedy G. J., Doren D. J., On the role of
small amines in zeolite synthesis. 2, J. Phys. Chem. B 104 (2000), 721-726
Rudolf P., Garcés J., Rietveld refinement of several structural models for mordenite that
account for differences in the X-ray powder pattern, Zeolites 14 (1994), 137-146
Sanders J. V., Crystallographic faulting in the mordenite group zeolites, Zeolites 5 (1985), 81-
90
Sano T., Wakabayashi S., Oumi Y., Uozumi T., Synthesis of large mordenite crystals in
presence of aliphatic alcohol, Micropor. Mesopor. Mater. 46 (2001), 67-74
Schulz-Ekloff G., Wöhrle D., van Duffel B., Schoonheydt R. A., Chromophores in porous
silicas and minerals: preparation and optical properties, Micropor. Mesopor. Mater. 51
(2002), 91-138
Scott G., Thompson R. W., Dixon A. G., Sacco A., The role of triethanolamine in zeolite
crystallization, Zeolites 10 (1990), 44-50
Sethian J. D., Raffray A. R., Latkowski J., Blanchard J. P., Snead L., Renk T. J., Sharafat S.,
An overview of the development of the first wall and other principal components of a laser
fusion power plant, J. Nucl. Mater. 347 (2005), 161-177
Shaikh A. A., Joshi P.N., Jacob N.E., Shiralkar V.P., Direct hydrothermal crystallization of
high-silica large-port mordenite, Zeolites 13 (1993), 511-517
226
Shepelev Yu. F., Butikova I. K., Smolin Yu I., Crystal structures of the partially K-, Rb-,and
Cs-exchanged forms of NaX zeolite in both the hydrated and the dehydrated (400 °C) states,
Zeolites 11 (1991), 287-292
Shi J., Anderson W., Carr S. W., Direct observation of zeolite A synthesis by in situ solid-
state NMR, Chem. Mater. 8 (1996), 369-375
Shiokawa K, Ito M., Itabashi K., Crystal structure of synthetic mordénites, Zeolites 9 (1989),
170-176
Sieber R., Decurtins S., Stoeckli-Evans H., Wilson C., Yufit D., Howard J. A. K., Capelli S.
C., Hauser A., A thermal spin transition in [Co(bpy)3] [LiCr(ox)3] (ox=C2O42-; bpy=2,2'-
bipyridine), Chem. Eur. J. 6 (2000), 361-368
Silva M., Freire C., de Castro B., Figueiredo J. L., Styrene oxidation by manganese Schiff
base complexes in zeolite structures, J. Mol. Catal. A 258 (2006), 327-333
Simoncic P., Armbruster T., Se incorporated into zeolite mordenite-Na: a single-crystal X-ray
study, Micropor. Mesopor. Mater. 71 (2004), 185-198
Simoncic P., Armbruster T., Peculiarity and defect structure of the natural and synthetic
zeolite mordenite: A single-crystal X-ray study, Am. Mineral. 89 (2004), 421-431
Simoncic P., Armbruster T., Cationic methylene blue incorporated into zeolite mordenite-Na :
a single crystal X-ray study, Micropor. Mesopor. Mater. 81 (2005), 87-95
Sinha H. K., Thomson P. C. P., Yates K., Excited state dipole moments and polarizabilities of
some aromatic alkenes and alkynes, Can. J. Chem. 68 (1990), 1507-1513
227
Smith G. S., Snyder R. L., FN: A criterion for rating powder diffraction patterns and
evaluating the reliability of powder-pattern indexing, J. Appl. Cryst. 12 (1979), 60-65
Smolin Y. I., Shepelev Y. F., Anderson A. A., Atomic scale mechanism of CaX zeolite
dehydration, Acta Cryst. B 45 (1989), 124-128
Song S. G., Crystal defects of mordenite structures, J. Mater. Res. 14 (1999), 2616-2620
Stranski I. N., Zur theorie des kristallwachstums, Z. Phys. Chem. 136 (1928), 259-278
Steiner T., The hydrogen bond in the solid state, Angew. Chem. Int. Ed. 41 (2002), 48-76
Szostak R., Handbook of molecular sieves, ed. Van Nostrand Reinhold, New York (1992)
Sun Y., Song T., Qiu S., Pang W., Shen J., Jiang D., Yue Y., Synthesis of mordenite single
crystals using two silica sources, Zeolites 15 (1995), 745-753
Takaishi T., Ordered distribution of Na ions in dehydrated NaX zeolite, Zeolites 17 (1996),
389-392
Tanaka Y., Kuroda H., Shionoya S., Generation of tunable picosecond pulses in the
ultraviolet region down to 197 nm, Opt. Commun. 41 (1982), 434-436
Thangaraj A., Kumar R., NMR study of soluble aluminosilicate precursors in zeolite Y
synthesis, Zeolites 10 (1990), 117-120
Thompson P., Cox D. E., Hastings J. B., Rietveld refinement of Debye-Scherrer synchrotron
x-ray data from Al2O3, J. Appl. Cryst. 20 (1987), 79-83
Tiwary S. K., Vasudevan S., Void geometry driven spin crossover in zeolite-encapsulated
cobalt tris(bipyridyl) complex ion, Inorg. Chem. 37 (1998), 5239-5246
228
Tung A. T., Park B.-H., Liang D. H., Niemeyer G., Laser-machined shape memory alloy
sensors for position feedback in active catheters, Sens. Actuat. A 147 (2008), 83-92
Twu J., Dutta P.K., Kresge C.T., Vibrational Spectroscopic Examination of the Formation of
Mordenite Crystals, J. Phys. Chem. 95 (1991), 5267-5271
Uytterhoeven L., Dompas D., Mortier W. J., Theoretical investigations on the interaction of
benzene with faujasite, J. Chem. Soc., Faraday Trans. 88 (1992), 2753-2760
Vitale G., Mellot C. F., Bull L. M., Cheetham A. K., Neutron diffraction and computational
study of zeolite NaX: influence of SIII' cations on its complex with benzene, J. Phys. Chem. B
101 (1997), 4559-4564
Wakayama Y, Mitsui T., Onodera T., Oikawa H., Nakanishi H., Growth and fluorescence
properties of perylene nanocrystals produced by ultra-rapid evaporation, Chem. Phys. Lett.
417 (2006), 503-508
229
Wakihara T., Sugiyama A., Okubo T., Crystal growth of faujasite observed by atomic force
microscopy, Micropor. Mesopor. Mater. 70 (2004), 7-13
Walton K. S., Abney M B, LeVan M. D., CO2 adsorption in Y and X zeolites modified by
alkali metal cation exchange, Micropor. Mesopor. Mater. 91 (2006), 78-84
Warzywoda J., Dixon A.G., Thompson R.W., Sacco A., Suib S.L, The role of the dissolution
of silicic acid powders in aluminosilicate synthesis mixtures in the crystallization of large
mordenite crystals, Zeolites 16 (1996), 125-137
Xiong G., Yu Y., Feng Z.-c., Xin Q., Xiao F.-S., Li C., UV Raman spectroscopic study on the
synthesis mechanism of zeolite X, Micropor. Mesopor. Mater. 42 (2001), 317-323
Xu D., Xue D., Ratajczak H., Morphology and structure studies of KDP and ADP crystallites
in the water and ethanol solutions, J. Mol. Struct. 740 (2005), 37-45
Yang X., Albrecht D., Caro J., Revision of Charnell’s procedure towards the synthesis of
large and uniform crystals of zeolites A and X, Micropor. Mesopor. Mater. 90 (2006), 53-61
Yang S., Navrotsky A., An in situ calorimetric study of zeolite crystallization kinetics,
Micropor. Mesopor. Mater. 52 (2002), 93-103
Yeom Y. H., Jang S. B., Kim Y., Song S. H., Seff K., Three Crystal Structures of Vacuum-
Dehydrated Zeolite X, M46Si100Al92O384, M = Mg2+, Ca2+, and Ba2+, J. Phys. Chem. B 101
(1997), 6914-6920
230
Yu Y., Xiong G., Li C., Xiao F.-S., Characterization of aluminosilicate zeolites by UV Raman
spectroscopy, Micropor. Mesopor. Mater. 46 (2001), 23-34
Zhang S.-Y, Talu O., Hayhurst D. T., High-Pressure Adsorption of Methane in NaX, MgX,
CaX, SrX, and BaX, J. Phys.Chem. 95 (1991), 1722-1726
Zeng L., Zha M., Ardoino M., Franzosi P., Zanotti L., Zuccalli G., Paorici C., Solution crystal
growth of urea and derivatives for nonlinear optical applications, J. Cryst. Growth 166
(1996), 528-532
Zhu L., Seff K., Olson D. H., Cohen B. J., Von Dreele R. B., Hydronium ions in zeolites. 1.
Structures of partially and fully dehydrated Na,H3O-X by X-ray and neutron diffraction, J.
Phys. Chem. B 103 (1999), 10365-10372
Zhu L., Seff K., Cation crowding in zeolites. Reinvestigation of the crystal structure of
dehydrated potassium-exchanged zeolite X, J. Phys. Chem. B 104 (2000), 8946-8951
Zones S.I., Hwang S.-J., Synthesis of high silica zeolites using a mixed quaternary
ammonium cation, amine approach : discovery of zeolite SSZ-47, Chem. Mater. 14 (2002),
313-320
Zones S.I., Hwang S.-J., Elomari S., Ogino I., Davis M.E., Burton A.W., The fluoride-based
route to all-silica molecular sieves ; a strategy for synthesis of new materials based upon
close-packing of guest-host products, C. R. Chimie 8 (2005), 267-282
231
Etude structurale par diffraction, absorption des rayons X
et simulations Monte-Carlo de matériaux zéolithiques
Les composites organiques-inorganiques avec la zéolithe comme matrice hôte sont de plus en plus
étudiés dans le but de concevoir des matériaux multifonctionnels conjuguant les propriétés de la molécule
organique et celle de la zéolithe (ex. : élaboration de commutateurs magnétiques/optiques, doubleurs de
fréquence). Les propriétés macroscopiques recherchées (SHG, transition de spin) dépendent des interactions
entre les espèces organiques invitées et la zéolithe (charpente, cations compensateurs de charge), qui sont encore
mal connue, en particulier en raison de l'absence de structures cristallographiques précises.
Cette thèse s'est intéressée à deux types de composites principalement : (Co2+(bpy)3-zéolithe X) réputé
présenter une transition de spin en température et dmpNA-MORdénite), à propriété potentielle de génération de
second harmonique.
Dans l’étude du composite magnétique, nous nous sommes particulièrement penché sur l'évolution de la
structure de la zéolithe précurseur Na,Co-X dans différents état d’hydratation. Des études complémentaires par
diffraction des rayons X sur monocristal, EXAFS et simulations Monte-Carlo ont mis en évidence la forte
interaction entre cations cobalt et charpente aluminosilicatée ainsi que la migration des cations Co2+
accompagnant l'hydratation progressive de la zéolithe Na,Co-X.
L’étude structurale du composite dmpNA-mordénite par diffraction des rayons X sur poudre
(synchrotron SPring-8, Japon) a montré que les canaux droits de la mordénite sont presque totalement occupés
par les molécules de dmpNA, qui adoptent un empilement quasi linéaire favorable à la Génération de Seconde
Harmonique.
Mots-clé : composite organique-zéolithe, Diffraction des Rayons X monocristal/poudre, simulations Monte-
Carlo, EXAFS, synthèse hydrothermale.
232