> = av. N
P(R)
=
1/R
where R = Si/Al. For compositions in the range R > 2 i t is a good approximation to take P$i as equal to this average value but at lower values of R there is considerable dispersion in the individual probabilities within the set { P^} . This is illustrated in Figure 7 which depicts the elements of {P^i^ for the formation of sodalite cages from an ensemble of 6R's for which Kflp = 0.5, as a function of composition. Figure 7 was obtained by an explicit enumeration of a l l possible sodalite
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
Silicon and Aluminum Ordering of Zeolites
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
MELCHIOR
1.5
2.0
2.5
-R-
Figure 7. Probabilities for formation of Si-O-Al linkages for various S i types upon formation of a sodalite cage from 6-ring sub-units, as a function of composition. The dashed l i n e i s the average for a l l sites < .> ι = 1/R. p
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
256
INTRAZEOLITE
CHEMISTRY
cage arrangements which can be formed from Μ , M , and P 6R's at S i / A l = 13/11, J.4/10, and 15/9, together with the assumption that Pjj£ = P(R) = 1/R f o r R > 2. The d i s p e r s i o n w i t h i n { Pfli } i s a r e f l e c t i o n o f the strong c o n s t r a i n t represented by Lowenstein's r u l e when the c o n n e c t i v i t y o f the s o d a l i t e cage i s taken i n t o account. In s t a t i s t i c a l terms, the p r o b a b i l i t y P ^ i o f adding an S i - O - A l linkage to the i t h s i atom i s s t r o n g l y dependent on the linkages formed i n the 3
2
2
(N-l)UL s t e p .
The s t a t i s t i c s of s o d a l i t e cage formation depicted i n Figure 7 apply whether the 6R u n i t s are considered as i s o l a t e d 6R s or as the hexagonal faces of a set of D6R u n i t s . Thus, the p r o b a b i l i t i e s i n F i g u r e 7 represent Ρβ£ f o r pathway I but P^£ f o r pathway I I . The remaining step which must be considered f o r each o f these pathways i s the formation o f D6R's from the 6 - r i n g faces o f a s o d a l i t e cage ( P 4 1 f o r Pathway I) or from i s o l a t e d 6R s u b - u n i t s (?3£ f o r pathway I I ) . Although the s t a t i s t i c s f o r these two i s o l a t e d steps are i d e n t i c a l , g i v e n the same mixture of M and P 6 - r i n g s , the steps d i f f e r i n an important r e s p e c t . The d i f f e r e n c e a r i s e s from the f a c t that the 6 - r i n g faces i n v o l v e d i n the step SC •> f a u j a s i t e f o r p a t h way I are not simply the members of the 6R ensemble involved i n the previous step 6R SC. Formation of a s o d a l i t e cage from four 6R s r e s u l t s i n the formation o f an a d d i t i o n a l set of four 6 - r i n g faces which at some compositions may have a d i f f e r e n t mix of M2 and P arrangements. T h i s can be seen most c l e a r l y by c o n s i d e r i n g the composition S i / A l = 1.67. As has already been noted, f o r K^p = 0, t o p o l o g i c a l c o n s t r a i n t s force the arrange ment designated as YX i n F i g u r e 1. Note that w h i l e YX i s formed from one M and three P 6R u n i t s , two of the four a d d i t i o n a l 6 - r i n g faces have a M arrangement. These M 6 - r i n g faces must be included i n the s t a t i s t i c s of the subsequent step i n which f a u j a s i t e i s formed. T h i s expansion can be described by an e f f e c t i v e r a t i o K ^ p . and must be r e f l e c t e d i n the p r o b a b i l i t i e s Ρ4£ f o r pathway I . T h i s p o i n t i s c r u c i a l i n that i t provides the opportunity o f d i s t i n g u i s h i n g pathway I and pathway I I .
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
f
2
2
1
2
3
2
2
2
F i g u r e 8 shows the D6R arrangements which have been considered f o r the c a l c u l a t i o n o f P4£ f o r pathway I and of Ρβ£ f o r pathway I I over the composition range 1 S i / A l 2, there i s s i g n i f i c a n t d e v i a t i o n f o r compositions i n the lower range R < 2. T h i s d e v i a t i o n i s most apparent f o r the composition R 1.67. Table I I g i v e s the t h e o r e t i c a l r e s u l t s at t h i s composition f o r K$fp = 0 and K^p = 0 . 5 as w e l l as the d i s t r i b u t i o n obtained by a l e a s t - s q u a r e s i n t e r p o l a t i o n of the experimental d a t a . n
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
CHEMISTRY
r
β
Ordered S o d a l i t e B u i l d i n g U n i t s Formulation of l o c a l S i , A l o r d e r i n g i n f a u j a s i t e i n terms o f pathway I , 6R -. SC -> f a u j a s i t e , represents an extension and refinement o f our previous treatment i n terms o f ordered SC s u b - u n i t s ( 4 ) . As we have noted, pathway I f o r K^p = 0 corresponds to formulation of f a u j a s i t e from the SC arrangements of F i g u r e 1. A c t u a l l y the set o f SC u n i t s s p e c i f i e d by K^p = 0 i s more e x c l u s i v e than that s p e c i f i e d p r e v i o u s l y ( 4 ) since there are minimum A l , A l next nearest neighbor SC u n i t s which cannot be formed using only M and P 6 R s . Moreover, c o n s i d e r a t i o n of 6R ensembles with K^p φ 0 provides a s t r a i g h t f o r w a r d approach to the study o f d e v i a t i o n s from complete l o c a l o r d e r . C a l c u l a t i o n of t h e o r e t i c a l 4-neighbor d i s t r i b u t i o n s [ S i ( n A l ) ] for pathway I i s a more i n v o l v e d task than f o r p a t h way I I , e s p e c i a l l y f o r a 6R ensemble with K^p ψ 0 . The f i r s t step i s to convert the 6R ensemble to a weighted average of s o d a l i t e cage s u b - u n i t s by a p p l i c a t i o n of the p r o b a b i l i t y set Ρβ£ which f o r t h i s pathway i s e x e m p l i f i e d by F i g u r e 7 f o r KMP = 0 . 5 . A d d i t i o n a l complexity a r i s e s because, i n a d d i t i o n to the 3-neighbor d i s t r i b u t i o n s f o r each S i s i t e , we r e q u i r e a c a l c u l a t i o n o f the r e s u l t a n t M / P r a t i o Kjjp., needed i n the next s t e p . T h i s value o f Kj.jp. i s then used to c a l c u l a t e 3
2
1
2
2
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
15.
MELCHIOR
Silicon and Aluminum Ordering of Zeolites
259
Table II Theoretical and Experimental A l Neighbor Distributions for Faujasite at Si/Al = 1.67
Theoretical (R = 1.67) [ S i ( n A l ) ] Pathway I
6R .. SC
KMP = 0 KMP = 0.5 Pathway II
n s
4^ 3^ 2 , 1, 0
Faujasite
[.125, .318, .390, .167, 0] [.161, .304, .345, .168, .023]
6R ... D6R + Faujasite
KMP = 0 KMP - 0.5
[.142, .288, .400, .170, 0] [.171, .285, .351, .172, .020]
Experimental (R = 1.67). [,130, .324, .369, .163, .015]
.0btained by computer f i t interpolation of data in Table I (11).
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
260
INTRAZEOLITE
the s t a t i s t i c a l distribution of D6R linkages and the associated probabilities used to convert the weighted average 3-neighbor distributions to the 4-neighbor distributions for faujasite. Theoretical A l neighbor distributions have been calculated for pathway I with K^p in the range 0-0.5 for the compositions N • 14 and N =15 (R = 1.40 and 1.67). The results for R - 1.40 are closely similar to those calculated for pathway I I . For R = 1.67, however, significantly better agreement with experiment i s obtained for pathway I. This can be seen in Table II which compares results for pathways I and II with the (interpolated) experimental distribution obtained by 2 9 i NMR. The crucial comparison is the relative intensities of Si(4Al) and Si (3A1) for the two pathways. Table II shows that good agreement i s obtained for pathway I for K^ps? 0. From these results we conclude that the observed 29$i NMR data are well accounted for by assuming pathway I and postulat ing essentially complete ordering in the 6R ensemble, equivalent to assuming distributions calculated for the ordered sodalite cages in Figure 1. Table III summarizes the results of this calculation for the compositions R = 13/11, 14/10, 15/9, 16/8, 17/7, and 18/8. The f i r s t line (a) for each composition in Table III i s the distribution calculated for the "pure" compositions, with an i n f i n i t e l y narrow distribution in N ]. For R = 13/11, sodalite cage X missing one A l atom is used, corresponding to a faujasite framework formed from and M 6R sub-units in the ratio 3:1. For compositions R 14/10, 15/9, and 16/8, respectively, only sodalite cages XY, YX, and Y are used, corresponding to frameworks formed from an ensemble of and P. 6R units. For these compositions the effective M /P ratio Kj.jp. must be used to calculate P ^ for the s t a t i s t i c a l distribution of D6R linkages. For XY at R = 14/10, MP. MP " > noted above, for YX at R - 15/9 Kflp. has the value 0.5. It i s precisely this observation that distinguishes pathways I and II. For R > 2 i t has been assumed that the constraint of the framework connectively i s sufficiently weak that P £ 2s 1/R. It should be noted that even for R >> 2 local order persists despite this assumption because we have assumed ordered p2 6R sub-units. The theoretical distributions (a) for these pure compositions may be compared with the results in (4). Even for integral values of N |, the associated 6R ensembles are composed of a mixture of 6R's with three, two, and one A l atoms and i t i s reasonable to assume a microscopic disper sion in Ν χ about the average Ν χ at each composition. The theoretical distributions in line (b) for each composition are obtained by applying a binomial distribution (1:2:1) to the pure distributions, in order to simulate this microscopic dispersion. The expected local dispersion in Ν χ would best be included in A1
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
CHEMISTRY
A1
S
A
2
β
2
K
2
=
K
0
b u t
a s
N
A
Α
Α
Α
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
15.
MELCHIOR
261
Silicon and Aluminum Ordering of Zeolites
Table I I I
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
Theoretical A l Neighbor Distributions for Ordered Sodalite Cage Building Units Si/Al
4 Al
3Al
2 Al
1 Al
0 Al
1.18
(a) .615 (b) .615 (d) (.605)
.308 .250 (.253)
.000 .077 (.087)
.000 .019 (.023)
.077 .038 (.030)
1.40
(a) (b) (c) (d)
.286 .317 .306 (.317)
.357 .339 .344 (.345)
.286 .245 .263 (.234)
.071 .080 .076 (.081)
0 .018 .013 (.024)
1.67
(a) (b) (c) (d)
.117 .125 .131 (.130)
.333 .317 .328 (.324)
.383 .392 .363 (.369)
.167 .167 .165 (.163)
0 0 .013 (.015)
2.00
(a) (b) (c) (d)
.000 .027 .014 (.026)
.250 .240 .237 (.243)
.500 .454 .479 (.442)
.250 .263 .257 (.263)
.000 .016 .013 (.032)
2.43
(a) (b) (c) (d)
.000 .000 .005 (.013)
.139 .148 .245 (.134)
.429 .420 .416 (.412)
.371 .362 .364 (.370)
.061 .070 .071 (.073)
3.00
(a) (b) (c) (d)
.000 .000 .002 (.000)
.074 .083 .079 (.058)
.333 .330 .326 (.356)
.444 .434 .435 (.468)
.148 .156 .158 (.123)
(a) Theory for sodalite cage sub-units, pure compositions. (b) Theory for binomial d i s t r i b u t i o n (1:2:1) of pure compositions. (c) Theory for pure compositions with K^jp 0.25. (d) Interpolated experimental data (unnormalized). β
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
262
INTRAZEOLITE
CHEMISTRY
the intermediate step 6 R -.• SC, i . e . , i n the c a l c u l a t i o n o f 3 neighbor d i s t r i b u t i o n s f o r the s o d a l i t e cage s u b - u n i t s . Prelimi nary r e s u l t s suggest that the approximation used to o b t a i n (b) i n Table I I I may overestimate somewhat the e f f e c t o f d i s p e r s i o n in T h e r e f o r e , we have included f o r comparison ( l i n e c, Table I I I ) the t h e o r e t i c a l r e s u l t s f o r 1 . 4 < R < 3 . 0 f o r an i n f i n i t e l y narrow d i s t r i b u t i o n i n Ν^χ but with a small value KMP = 0 . 2 5 . to i l l u s t r a t e the e f f e c t o f s l i g h t d e v i a t i o n s from complete P order i n the 6 R ensemble. The i d e a l theory would i n c l u d e both sources o f l o c a l d i s o r d e r , which refinement i s l e f t to a future p u b l i c a t i o n . Table I I I a l s o i n c l u d e s on l i n e (d) f o r each composition the i n t e r p o l a t e d experimental r e s u l t obtained from a computer f i t to the ^ S i NMR data (11). The p r e d i c t e d d i s t r i b u t i o n s (b) i n Table I I I are compared with the experimental data i n F i g u r e 4 . The agreement o f theory and experiment i s very good across the e n t i r e composition range, lending strong support to the c o n c l u s i o n that f a u j a s i t e i s formed from an ordered 6 R ensemble v i a the pathway 6R •+ SC -». f a u j a s i t e .
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
2
2
S i , A l Ordering i n ZK4 I t has a l r e a d y been noted that the ordered s o d a l i t e cage model which d e s c r i b e s the S i , A l d i s t r i b u t i o n i n f a u j a s i t e q u i t e w e l l does not account f o r the Z K 4 % i NMR data over the range 1 < R < 3 . I t has been shown that the consequences o f c o n s i d e r ing ordered 6 R s u b - u n i t s are c l o s e l y r e l a t e d to a d e s c r i p t i o n i n terms o f ordered s o d a l i t e cages. I t seems reasonable, t h e r e f o r e , to consider the formation o f Z K 4 by one or another o f the pos s i b l e pathways i n v o l v i n g 4 R ' s ( c f . F i g u r e 2 ) . In the f o l l o w i n g we d i s c u s s the o r d e r i n g i n ZK4 i n terms o f the pathway Τ - . 4 R .. D4R Z K 4 . For the case o f ZK4 no d i s t i n c t i o n can be made between t h i s pathway and the a l t e r n a t i v e 4 R - . S C - . Z K 4 . 2
C a l c u l a t i o n o f the p r e d i c t e d [ S i ( n A l ) ] d i s t r i b u t i o n f o r ZK4 constructed from an ensemble o f 4 R s u b - u n i t s i s f o r m a l l y s i m i l a r to the foregoing treatment of f a u j a s i t e based on 6 R s u b - u n i t s . The f i r s t step i s to c a l c u l a t e the s t a t i s t i c a l mix of D 4 R u n i t s which can be formed over the composition range 1 < R < 3 . In order that the & 4 next nearest neighbor i n t e r a c t i o n s be m i n i mized we exclude S14 u n i t s from our 4 R ensemble. The a p p r o p r i ate mixture o f D 4 R ' s i s then assembled to form a ZK4 c r y s t a l , subject to the c o n s t r a i n t o f Loewenstein's r u l e as i t a p p l i e s to the topology o f s o d a l i t e cages formed from s i x D 4 R ' s . Figure 8 shows the c a l c u l a t e d p r o b a b i l i t i e s P 4 Î o f adding an aluminum neighbor f o r the step D4R•+ Z K 4 , a r e s u l t obtained i n the same manner as was F i g u r e 7 f o r f a u j a s i t e . For the present case the t o p o l o g i c a l c o n s t r a i n t i s much l e s s severe because the 4 R s u b - u n i t imposes no l o c a l order beyond Loewenstein s r u l e . The number of p o s s i b l e s o d a l i t e cage arrangements i s much greater and the r e s u l t i n g d i s t r i b u t i o n corresponds to what we might term 1
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
15.
MELCHIOR
Silicon and Aluminum Ordering of Zeolites
263
"Loewenstein s t a t i s t i c a l , " subject to no additional constraint. We find that the distributions [(Si(nAl)] calculated i n this manner deviate s i g n i f i c a n t l y from the observed distributions for ZK4 at low S i / A l . That the detailed s t a t i s t i c a l application of Loewenstein's rule does not account for the observed ZK4 data implies that a further constraint must be imposed i n the step D4R ZK4, i . e . , there are next nearest neighbor interactions which have not been considered. On the other hand, reasonably good agreement i s obtained over the composition range 1.4 £ R £ 3.0 by taking P4i to be equal to i t s average = 1/R for a l l S i sites (Figure 9, dashed l i n e ) . I t can be shown that the averaging of
1.0
1.5
2.0 Si/Al
2.5
3.0
Figure 9. Probabilities for formation of Si-O-Al linkages for various S i sites upon formation of a sodalite cage from 4-ring sub-units i n ZK4.
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
264
INTRAZEOLITE
P^£ corresponds to a nonstatistical mixture of sodalite cages formed by D4R ZK4. Not unexpectedly, this s t a t i s t i c a l bias can be expressed in terms of a preference for para arrangements of A l atoms in the eight 6-rings formed by the combination of six D4R units. Consider the following example for one composi tion Si/Al = 1.4. The s t a t i s t i c a l mixture of sodalite cages, that used to calculate in Figure 9 at this composition,can be characterized as having a M /P ratio of 5:1. As can be seen by Figure 9 the probabilities for different S i environ ments are significantly different from the average = 1/R. However, for a non-random combination for which the resulting M2/p2 £ i 1 i t can be shown that a l l S i sites have P4£ within 5 percent of the average value of 1/R. This result implies a significant next nearest neighbor constraint on the arrangements favored in ZK4 crystallization, comparable to that i n faujasite. In principle at least we could carry out a detailed calculation of Έ^ι for D4R ZK4 based on some formulation to express the preference for P 6-rings leading to a ZK4 structure composed of p a r t i a l l y ordered sodalite cages. As yet such calculation has not been achieved. 2
s
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
CHEMISTRY
2
:
2
Conclusions The observed [Si(nAl)] distributions for faujasite have been satisfactorily accounted for by a model which assumes an ensemble of P2 ordered 6-rings assembled into a faujasite l a t t i c e v i a the pathway 6R -> SC •. faujasite, subject only to the topological constraints imposed by s t r i c t adherence to Loewenstein s rule. Rationalization of the [Si(nAl)] distributions for ZK4 i s less straightforward because the A l , A l next nearest neighbor inter action (responsible for the preference for P 6-rings in faujasite) must be taken into account for the formation of sodalite cages from 4R or D4R sub-units. The data for ZK4 and, more importantly, comparison of the ZK4 data with the faujasite data leads to the conclusion that ZK4 i s formed by a pathway involving 4-rings as the fundamental unit. Both ZK4 and faujasite exhibit short range s t a t i s t i c a l order which can be rationalized in terms of a competition between the s t a t i s t i c s of the step in which the sodalite cage i s formed and the tendency of the system to avoid A l , A l next nearest neighbor interactions. The basic difference between the two systems can be understood in terms of the difference in the severity of the constraint imposed by Loewenstein.s rule for the formation of a sodalite cage from four ordered 6R's in the case of faujasite or from six 4R's or D4R's in the case of ZK4. The local order in faujasite can be considered as "frozen in" by the topological constraint imposed by the ordered 6R sub-units. As such this order i s more a property of the sub-unit than a manifestation of next nearest neighbor interactions in the faujasite crystal. By contrast the next nearest neighbor interaction in ZK4 cannot be 1
2
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.
15.
MELCHIOR
Silicon and Aluminum Ordering of Zeolites
265
accommodated i n a s i n g l e ordered s u b - u n i t and the r e s u l t i n g l o c a l order may be more a p p r o p r i a t e l y considered a property o f the z e o l i t e framework.
Literature Cited 1. 2.
Downloaded by UNIV OF PITTSBURGH on March 11, 2016 | http://pubs.acs.org Publication Date: May 17, 1983 | doi: 10.1021/bk-1983-0218.ch015
3.
4. 5. 6. 7. 8. 9. 10. 11.
Lippmaa, E. T.; A l l a , Μ. Α.; Pehk, T.J.;Engelhardt, G. M.; J . Am. Chem. Soc., 1978, 100, 1929. Englehardt, G.; Lohse, V.; Lippmaa, E.; Tarmak, M.; Magi. M.; Z. Anorg. A l l g . Chem., 1982, 482, 49. Klinowski, J.; Ramdas, S.; Thomas, J . M.; Fyfe, C. Α.; Hartman, J . S.; J . Chem.Soc.,Faraday Transactions, 1982, 78, 1025. Melchior, M. T.; Vaughan, D. E. W.; Jacobson, A. J.; J . Am. Chem. Soc., 1982, 104, 4859. Lippmaa E.; Magi M., Samoson; Α., Tarmak; M., Englehardt, G.; J . Am. Chem. Soc., 1981, 103, 4992. Bursill, L. Α.; Lodge, Ε. Α.; Thomas, J . Μ.; Cheetham, Α. Κ.; J . Phy. Chem., 1981, 85, 2409. Melchior, M. T.; Vaughan, D. E. W.; Jarman, R. H.; Jacobson, A. J.; Nature, 1982, 298, 455. Jarman, R. H.; Melchior, M. T.; Vaughan, D. E. W.; Chapter __, in this book. Loewenstein, W.; Am. Miner., 1954, 39, 92. Hexem, J . G.; Frey, M. H.; Opella, S. J.; J . Chem. Phys. 1982, 77, 3847. Jarman, R. Η., personal communication.
RECEIVED
February 25, 1983
Stucky and Dwyer; Intrazeolite Chemistry ACS Symposium Series; American Chemical Society: Washington, DC, 1983.