V. K LASTROWA I POTENCJAŁOWA TEORIA ADSORPCJI
1. Klastrowy opis mechanizmu adsorpcji
1.2. Modele klasy LBET
1.2.1. Model powiązań geometrii porów i energii adsorpcji
Podążając za opracowaniem Promotor niniejszej rozprawy i wsp. [40], w celu wyprowadzenia klasy modeli LBET przyjmuje się następujące założenia i uproszczenia, stanowiące podstawy modelu powiązań geometrii porów i energii adsorpcji w materiałach mikroporowatych:
1. Maksymalna liczba warstw k jest czynnikiem różnicującym typy klastrów, a profil energii κ służy indeksowaniu klastrów jedynie w obrębie ktego typu. Założenie V.1.2.1-1 stanowi uproszczenie w oznaczeniu typologii klastrów. Klastry różniące się maksymalną liczbą warstw k są jednocześnie zróżnicowane pod względem profili energii κ, ze względu na energię w warstwach: k i k+1. A zatem, w celu uproszczenia oznaczeń wyróżnia się klastry ktego typu, w obrębie których indeksuje się klastry o różnych profilach energii κ [41].
2. Rozkład
i
pierwotnych miejsc adsorpcyjnych ktego typu mhAk (mpκ0), uwzględniony w równaniu (V-60), dany jest funkcją eksponencjalną [41]:
1
1 k
hA hAk m
m (V-61)
gdzie: mhA całkowita liczba pierwotnych miejsc adsorpcyjnych dostępna dla cząsteczek adsorbatu w przestrzeni porów adsorbentu, [(mmol/g)]; empiryczny parametr rozkładu funkcji eksponencjalnej, odpowiadający
Zapoczątkowanie i wzrost klastra ktego typu powiązane jest z istnieniem pierwotnych miejsc adsorpcyjnych mhAk (mpκ0) o profilu energii pozwalającym ulokować k cząsteczek i wynika z całkowitej liczby miejsc dostępnych w porach adsorbentu oraz typu ograniczenia mechanizmu wzrostu klastra.
Przyjmując, iż
i
w materiałach mikroporowatych, dominujący wpływ na przebieg
i
mechanizmu adsorpcji ma zjawisko zapełniania mikroporów, całkowita liczba pierwotnych miejsc adsorpcyjnych wyznaczana jest z zależności [41]:
Formuła (V-62), wykazuje interesującą właściwość wynikającą z przyjęcia niezależności pojemności monowarstwy VhA od objętości molowej Vs. A zatem, obszar dostępny dla ulokowania cząsteczki w pierwszej warstwie VhA jest niezależny od jej rozmiarów, pod warunkiem zastosowania adsorbatów o niewielkich rozmiarach cząsteczkowych (tj. N2, H2 itd.). Tym samym, możliwe jest zastosowanie różnorodnych adsorbatów, w celu zwiększenia wiarygodności analiz [40].
Materiały mikroporowate o nieregularnej strukturze, przeważnie wykazują typ II mechanizmu ograniczenia wzrostu klastra, na skutek dominującej ultra- i supermikroporowatości [129]. A zatem, czynnik wykładniczy k-1 równania (V-61) odpowiada ograniczeniu liczby warstw w klastrze i odwzorowuje geometryczny obraz rozmiaru porów materiału. Interpretacja, wartości k odpowiada maksymalnej długość stosu cząsteczek, które mogą zostać ulokowane w porze. W przypadku porów szczelinopodobnych, wyraża zatem rozmiar poru odpowiadający długości poru, z kolei w przypadku porów kulopodobnych odpowiada średnicy [11].
Empiryczny parametr funkcji eksponencjalnej maksymalnej liczby warstw w klastrze charakteryzuje geometryczny obraz struktury porów, a zatem jest ściśle powiązany z kształtem rozkładu energii pierwotnych miejsc adsorpcyjnych.
Przyjmując założenie V.1.2.1-2 w modelach klasy LBET wprowadza się dwa, wyznaczane empirycznie niezmiennicze parametry struktury porów: i VhA.
3. Energia
i
adsorpcji pierwszej warstwy QAκ (V-54), dana jest funkcją rozkładu
Q RT
fk A wynikającą
i
ze specyficznych właściwości geometrycznych powierzchni
i
materiału [41].
W rzeczywistych materiałach mikroporowatych przyjmuje się, iż czynnikiem różnicującym energię adsorpcji w pierwszej warstwie QAκ jest parametr ZAκ, reprezentujący efekty właściwości geometrycznych pierwotnych miejsc adsorpcyjnych (V-54), [41]. Przy czym, efekty te są dominujące, gdy lokowana cząsteczka adsorbatu jest apolarna (tj. Cp = 1 = const (V-53)). Z kolei, w przypadku cząsteczek adsorbatu polarnego (Cp ≠ 1 ≠ const) mogą występować równocześnie z oddziaływaniami specyficznymi. s hA hA V V m (V-62)
gdzie: VhA pojemność monowarstwy, [(mm3/g)];
s
Rozważając nieregularność struktury rzeczywistych materiałów mikroporowatych jako źródło heterogeniczności geometrycznej (tj. również energetycznej) przyjmuje się, iż efekty właściwości geometrycznych wynikają ze zróżnicowania maksymalnej liczbą warstw k cząsteczek adsorbatu ulokowanych w klastrze [41]. A zatem, w klastrach ktego typu można przyjąć równomierny rozkład ZAκ w granicach (Zfk, ZAk), przy czym zbiór dostępnych wartości rozkładu wzrasta wraz ze wzrostem maksymalnej liczby warstw k [41].
W klastrze
i
1ego typu spodziewany jest najlepszy kontakt ze względu na cząsteczkowe
i
rozmiary dostępnego obszaru. Jeśli, zatem przyjmie się ZA = ZA1, to wartość
i
parametru ZAk ZAdla każdego k > 1, spełnia zależności: ZA ≥ ZAk i Zf1 ≤ Zfk [41]. A zatem, dla każdego
i
klastra ktego typu, gdzie K = 1, …, k, generuje
i
się dwa parametry granicy rozkładu ZAk. W celu zredukowania liczby
i
tych parametrów przyjmuje się model powiązań właściwości geometrycznych pierwotnych
i
miejsc adsorpcyjnych
i
z energią adsorpcji względem ZA [41]:
Mając
i
na uwadze powyższe rozważania, analogicznie można przyjąć równomierność rozkładu energii adsorpcji
i
QAκ (V-54) w granicach (QAk, Qfk). A zatem, energię
i
adsorpcji w pierwszej warstwie QA, uwzględniającą heterogeniczność energetyczną
i
powierzchni porów materiału, można wyrazić [41]:
Z kolei, schodkowy
i
rozkład funkcji energii adsorpcji fk
QA RT
, o przedziałach odpowiadających (QAk, Qfk), wyraża się [41]:Na podstawie
i
(V-58) oraz (V-65), definiuje się również schodkowy rozkład bezwymiarowego parametru energetycznego BAκ [41]:
A Ak def Ak Z Z , A fk d ef fk Z Z oraz A11 (V-63)
Ak k cp A p d ef A U Z Q Q Q 1 min , gdzie: QAk UpZAkQcpQAZA
1Ak
Qcp QAf UpZfkQcpQAZA
1fk
Qcp (V-64) Ak fk cp A Ak fk A k Q Z RT Q Q RT RT Q f 1 (V-65)
A
fk Ak
A A
Ak fk
cp A k B B Z B B B B f 1 ln 1 1 ln 1 , gdzie BA
BAk,Bfk
(V-66)gdzie: BA rozkład bezwymiarowego parametru energii adsorpcji w warstwie n = 1 klastra klasy κ;
cp
B bezwymiarowy parametr molowej energii adhezji idealnego kontaktu adsorbat-adsorbent: RT Q B cp d ef cp exp ; Ak
B granica rozkładu parametru energii adsorpcji BAκ:
ZA Ak cp A Ak def Ak B B RT Q B 1 exp ;Założeniem V.1.2.1-3 w modelach klasy LBET wprowadza się parametr energetyczny: QA (BA).
4. Energia
i
adsorpcji Qκn (V-52) warstw n = 2, …, k-1, wszystkich klastrów klasy κ (t.j. ktego typu) jest stała i wyraża się [41]:
Przyjęty
i
w założeniu V.1.2.1-2, dominujący wpływ zjawiska zapełniania mikroporów
i
na przebieg mechanizmu adsorpcji, pozwala przyjąć niewielki zakres zmian
i
parametru geometrycznego kontaktu adsorbat-adsorbent Zcpκn wzdłuż wszystkich
i
warstw n = 2, …, k-1. Wartość Zcpκn różnicująca energię Qκn (V-52) jest, bowiem
i
w mniejszym stopniu zależna od samej geometrii poru niż przypadkowego geometrycznego
i otoczenia pierwotnego centrum adsorpcyjnego. Stąd, w materiałach
i
mikroporowatych o nieregularnej strukturze, dominujący udział ultra- i supermikroporów
i
w całkowitej objętości porów materiału, stwarza warunki do zapewnienia
i
dobrego kontaktu Zcpκn w całej objętości poru, wzdłuż warstw n = 2, …, k-1. Skutkiem
i
tego, przyjmuje się stałość parametru Zcpκn = ZC = const, energii
i
adsorpcji
i
Qκn = QC = const, a także parametru energetycznego Bκn = BC (V-58) dla wszystkich klastrów klasy κ w warstwach n = 2, …, k-1 [41]:
w efekcie, założeniem V.1.2.1-4 wprowadza się homogeniczność warstw n = 2, …, k-1 oraz parametr energetyczny QC (BC).
5. Energia
i
adsorpcji Qκn (V-52) w najwyższej warstwie n = k, wynikająca z mechanizmu ograniczenia wzrostu klastra, dana jest formułą [41]:
W efekcie przyjęcia homogeniczności warstw n = 2, …, k-1 w klastrze (założenie V.1.2.1-4), znacząco redukuje się liczbę {Qκn; κ = 1, …, n = 2, …, k-1} parametrów modelu. Jednakże, istnienie typów ograniczenia wzrostu klastrów wymusza zróżnicowanie energii adsorpcji Qκn w warstwie n = k oraz hipotetycznej warstwie n = k+1. Stąd, w celu uproszczenia analiz numerycznych, przyjęto powiązanie energii QC homogenicznych warstw n = 2, …, k-1 z energią najwyższej warstwy n = k.
fk
B granica rozkładu parametru energii adsorpcji Bfκ:
ZA fk cp A def fk B B B 1 ; AB parametr energii adsorpcji w pierwszej warstwie QA: RT Q B def A A exp . const. Q Qκn C , QC Up(12Zpp)ZCQcp oraz QC0 (V-67) gdzie: ZC skumulowany i
parametr geometryczny korygujący wartość energii Qcp ze względu
i
na niedoskonałość kontaktów adsorbat-adsorbent, spełniający izależność: i ZC Up(12Zpp) Qcp. RT Q B C def
C exp oraz BC1 oraz QC0 (V-68)
dla typu I: Qκk QC const. oraz Qκ k1 , Bk BC
(V-69) dla typu II: QκkQkkQC,
k C k k k kk k B B B 1 1 1 i Bkk BC 0
I tak, ograniczenie mechanizmu wzrostu klastrów ze względu na wzrost konkurencyjny (tj. typ I mechanizmu), podążając za definicją stopnia pokrycia ktej warstwy w równaniu (V-59), pozwala przyjąć wartość Bκk wynikającą ze stopnia pokrycia warstwy niższej n = k-1, stąd: θκk = θk(k-1). A zatem, przyjmuje się niższą wartość Qκk niż w przypadki mechanizmu typu II, ale taką samą w obrębie wszystkich klas klastrów κ = k i ograniczoną wartością QC [41].
Z kolei, w przypadku geometrycznie lub energetycznie ograniczonego wzrostu klastrów (tj. typ II mechanizmu), explicite przyjmuje się homogeniczność w warstwie n = k. A zatem, formuła energii adsorpcji najwyższej warstwy: Qκk = QC, a ograniczenie wzrostu klastra wynika z formalnych właściwości hipotetycznej warstwy n = k+1 [41].
6. Dla
i
każdego klastra ktego typu, przyjmuje się stałą wartość parametru kształtu
i
klastra βκn dla n > 1 (V-55), odpowiadającą wartości uśrednionego całkowitego parametru β [41]:
Przyjęcie stałej wartości parametru kształtu klastra βκn ≡ β ≡ const jest następstwem przyjęcia homogeniczności warstw n > 1 i pozwala znacząco ograniczyć liczbę parametrów modelu [41]. Przy czym, uproszczenie to jest akceptowalne w materiałach mikroporowatych o nieregularnej strukturze, ze względu na ograniczenie rozmiaru porów do ultra- i supermikroporów [131]. Niemniej, należy mieć na uwadze, iż parametr β charakteryzuje głównie kształt małych klastrów, a w przypadku klastrów większych odpowiada kształtowi otoczenia pierwotnego miejsca adsorpcyjnego [41].
Wzrost liczbowej wartości parametru β interpretuje się, jako rozgałęzianie klastrów, gdzie wartość graniczna odpowiada silnie rozgałęzionym strukturom dendrytycznym [41]. W przypadku, klasycznych stosopodobnych klastrów teorii BET, wartości parametru odpowiada β = 1, co jest równoznaczne z wytworzeniem tylko jednego miejsca dla ulokowania cząsteczki adsorbatu. Podobnie, zgodnie z założeniem V.1.1-2 przyjmuje się wartości β = 1 dla cząsteczki w warstwie n = 1.
Przyjmując założenie V.1.2.1-6 w modelach klasy LBET wprowadza się kolejny wyznaczany empirycznie niezmienniczy parametry struktury porów: β.