0% found this document useful (0 votes)
189 views68 pages

Présentation PowerPoint PDF

Uploaded by

Assis Rocha
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
189 views68 pages

Présentation PowerPoint PDF

Uploaded by

Assis Rocha
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 68

CO2 maîtrisé | Carburants diversifiés | Véhicules économes | Raffinage propre | Réserves prolongées

Physisorption and pore size


analysis

Gerhard Pirngruber
Direction Catalyse et Séparation
[email protected]
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
CV Gerhard Pirngruber

 1995 – 1999 PhD in heterogeneous catalysis


 Universiteit Twente (NL)
 2000 – 2005 Oberassistent
 ETH Zurich, Inst of Chemical and Bioengineering
 Research activity
 zeolites
 mesoporous silica
 2005 – now R&D Scientist, project leader
 IFP Energies nouvelles, Catalysis and Separation Division
 Research activity
 2005 – 2012 : CO2 capture, separation of hydrocarbons
 2012 – now: hydrotreating and hydrocracking catalysts
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Porous solids in catalysis
 Supported catalysts
 e.g. noble metals, transition metal sulfides
 Porous support
 provides surface for dispersing the catalyst
 provides mechanical stability
 has an influence on diffusion of reactants and products
 has an influence on adsorption of reactants and products
 determines reactor volume
 Porous catalysts
 Zeolites Optimising porosity is often the
 Oxydes key to improving catalytic
performance.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Types of porosity

 Porous cristalline structures


 Zeolites,
 Metal Organic Frameworks

 Ordered or disordered stacking of small particles


generates an interparticles porosity
 oxides (silica, alumina, etc.)
 activated carbon
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Characterisation of porosity

 Physisorption (N2, Ar, Kr)


 surface area
 pore volume
 pore size distribution
 Hg intrusion
 pore volume
 pore size distribution
 surface area
 Imaging methods
 Transmission Electron Microscopy
 Secondary Electron Microscopy
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Terminology
 Adsorption is a surface phenomenon.
 Adsorbent = solid that provides a surface for
adsorption.
 Adsorbate = molecule adsorbed on the surface.
 Adsorptive = molecule susceptible of being adsorbed.

 Physisorption: no formation of chemical bonds, no


transfer of electron density
 Chemisorption: formation of a chemical bond
(exchange of electrons)
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Terminology

 IUPAC distinguishes 3 catagories of pore sizes


 Micropores < 2 nm
 Mesopores 2 – 50 nm
 Macropores > 50 nm
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Forces involved in physisorption
 van der Waals forces
 amplified by multiple interactions with atoms or pore wall (i)
5000
   ij 12   ij 6 
4000 U LJ  4        
 r   r  
3000  
Energy

2000
adsorbate
1000 i j
0

-1000
0 2 4 6 8 pore wall
Distance

 maximal when size of adsobate (j) close to pore size

pore wall
adsorbate
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Forces involved in physisorption
 Electrostatic forces
 Any non-symmetric charge distribution in the adsorbent
generates an electric field
 adsorbates with an electric moment (dipole, quadrupole)
interact with the electric field
 undesirable in analysis of porosity/pores size because
surface chemistry/chemical composition enter into the game
 N2 has a small quadrupole moment → use of Ar, Kr is
preferable for precise micropore size analysis
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Adsorption in micropores
 Gradual filling of pore volume starting
from surfcae
Saturation zone
Adsorbed amount

Transition zone
additional adsorbate-
adsorbate interactions

Initial slope depends on adsorbent-


adsorbate interactions
Pressure
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Adsorption in mesopores

monolayer multilayer onset of pore filled by


adsorption adsorption capillary capillary
condensation condensation
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Adsorption in macropores

 Opposite pore wall is too far away to influence


adsorption and provoke capillary condensation
 First monolayer, then multilayer adsorption
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
How to measure an adsorption isotherm?
 Volumetric method
Réservoir étalon : Vres Réservoir étalon : Vres

Etat initial : PriHe, TriHe Etat initial : PriN2, TriN2


Hélium Azote
Etat final : PrfHe, TrfHe Etat final : PrfN2, TrfN2

Bilans Matières : He et N2

 PrHe
i f
PrHe V
nHe   i  f  res

 TrHe TrHe  R

Température nHe  
V
 Pc f  Pci cel
RTm
ambiante

Cellule mesure : Vcel


 PrNi 2 PrNf 2  Vres
nN 2   i  f 
Température  TrN 2 TrN 2  R
 
N2 liquide V
77 K nN 2  Pc f  Pci cel  nNads2
RTm

Etat initial : PciHe, TciHe Etat initial : PciN2, TciN2


Etat final : PcfHe, TcfHe Etat final : PcfN2, TcfN2

Masse d'adsorbant : Msol


Réservoir d'azote
© IFP

liquide à 77 K Quantité adsorbée : Qads=nN2ads/Msol

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Material balance

 Reservoir : 1 Sample cell : 2


 without adsorption (He), sample cell initially under
vacuum
pinitV1 p finalV1 p finalV2 V2 V1 pinit  p final
   
RT1 RT1 RT2 T2 T1 p final

 with adsorption (N2) determined with He

pV1 pV2
quantity initially ntot    nads
present in the RT1 RT2
reservoir
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Experimental difficulties

 Precise measurement of pressure required


 difficult at low pressure
 difficult to measure a micropore distribution
 Isotherm !!
 T2 in principle 77 K (temperature of liquid N 2 )
 part of the cell is not plunged into liquid N2, but is at ambient
temperature
 temperature gradient must be the same as during the initial
calibration of the volume with He
 pay attention to the level of liquid N 2
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Example zeolite NaX

220

powder
200
Adsorbed volume ( cm3/g)

180

160 after shaping with a


N2 at 77 K NaX (ADS 57701, Powder)
binder
140 N2 at 77 K NaX (ADS 48100) SPX 3003

N2 at 77 K NaX (ADS44956) SPX 3001


120

100
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
P/P0

 Micropores are entirely filled at very low pressure.


 Isotherm is totally flat once micropores are filled.
© IFP

 saturation zone
Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Examples - ZnO
180
Quantité adsorbée (cm3 STP/g) 160 capillary condensation
140 in the mesopores
120

100

80
multilayer adsorption on the ZnO

surface of the particles


60

40

20

0
0 0.2 0.4 0.6 0.8 1
p/p0

 no micropores
 mesopores generated by stacking of inidividual
particles
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Comparison of two ZnO samples
180

160
Quantité adsorbée (cm3 STP/g)

140

120

100 65811
80 70299

60

40

20

0
0 0.2 0.4 0.6 0.8 1
p/p0

 Surface : 65811 > 70299


 Pore size : 70299 > 65811
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Example – Al2O3
300

capillary condensation
Quantité adsorbée (cm3 STP/g)

250
in mesopores
200

150

100
multilayer adsorption
50
on surface of particles
0
0 0.2 0.4 0.6 0.8 1
p/p0

 no micropores
 high surface area, broad pore size distribution
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Comparison Al2O3
500
450
Quantité adsorbée (cm3 STP/g)

400
350
300 79999
250 61399
200 47148

150
100
50
0
0 0.2 0.4 0.6 0.8 1
p/p0

 Surface : 79999 > 61399 > 47148


 Pore size : 47148 > 61399 > 79999
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Quantitative treatment of N2 isotherms

 Pore volume: Vmicro, Vmeso, Vtotal


 t-plot
 Dubinin-Radushkevitch
 Surface area
 BET (Brunauer-Emerett-Teller)
 Langmuir
 Pore size distribution
 Micropores
 Horwath-Kawazoe, Saito-Foley
 Mesopores
 BJH (Brunauer-Joyner-Halenda)
 DFT (Density Functional Theory)
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Determination of the surface area

 BET (Brunauer Emmet Teller): multilayer adsorption


 first layer: adsorption on the surface
 all the other layers are considered like a condensed liquid

Kcondensation
Ksurface

nm = number of surface adsorption sites (monolayer)


nads = total number of adsorbed molecules
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
BET equation
P 1 C 1  P 
   
 BET equation Vads P0  P  VMC V M C  P0 
K surface
C
K condensation Recommended
range of p/p0 :
 VM = monolayer volume 0,05 - 0,35 max.
 Vads = adsorbed volume
Plot p/Vads(p0-p) vs. p/p0

1/[V.(P/P0-1)]

 Calculation of surface area:
NAv : Avogadro number
PNT. = (C-1)/(VM.C)

 a : area of N2 molecule (16,2 Ǻ2)
O.O. = 1/(VM.C)

0,1 0,2 0,3 0,4 P/P0


p0VM
S BET  a. .N A S BET (m2 / g )  4.355 Vm (cm3 STP / g )
R  273K
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Exemple ZnO
Isothermes BET plot
4.355 slope
180
S BET  CBET  1
160 0.05
slope  int ercept int ercept
Quantité adsorbée (cm3 STP/g)

0.045
140
0.04 y = 0.1196x + 0.0005
120
0.035

prel / Vads(1-prel)
100 65811 0.03
80 70299 65811
0.025
70299
60
0.02
40 0.015
20 0.01 y = 0.0594x + 0.0003
0 0.005
0 0.2 0.4 0.6 0.8 1 0
p/p0 0 0.1 0.2 0.3 0.4 0.5
prel = p/p0

SBET (m2/g) CBET


65811 74 136
70299 36 237
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Example NaX

 The multilayer adsorption 0.0007

model does not apply well to 0.0006 y = 0.0064x - 3E-05

microporous solids when the 0.0005

prel / Vads(1-prel)
pores are filled. 0.0004

 Choose a lower pressure 0.0003

range for microporous solids. 0.0002

 p/p0 = 0.05 – 0.10 0.0001

C constant may be negative.


0
 0.0000 0.0200 0.0400 0.0600 0.0800 0.1000 0.1200

Does not make physical sense prel = p/p0

 lower the pressure range


4.355
further S BET   685m 2 / g
0.0064  0.00003
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Langmuir Surface Area

 Langmuir theory
 is a model of monolayer adsorption
 this hypothesis is never fulfilled
 Do not use Langmuir surface areas !!
 SLangmuir is always > SBET because multilayer
adsorption is treated like monolayer adsorption.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Determination of the pore volume
300

250
Volume adsorbé (cm3 STP/g)

200

Vads = 240
150
cm3 STP/g
100 → Vtotal
Vads = 178
cm3 STP/g 50
Vmeso = Vtotal - Vmicro
→ Vmicro
0
0 0.2 0.4 0.6 0.8 1
p/p0

 In some cases possible without using any model.


© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Conversion gas volume – pore volume

 Volume of adsorbed gas correponds to a certain


number of moles
Vads p
 STP = Standard Temperature Pressure nads 
 T = 273 K, p = 1 bar = 100 kPa RT
 The adsorbed phase is considered like a liquid phase
 Density of liquid N2 at 77 K : 0.807 g/ml

mN 2 nads  M N 2
Vp  
 N 2,liq  N 2,liq

Vp (ml / g )  1.5468 103 Vads (cm3 STP / g )


© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
t-plot

 Transform adsorbed amount in an average thickness


(t) of the adsorbed layer
 For multilayer adsorption on a flat surface Vads
t
S
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Universal curve t vs. p/p0
 The isotherms (Vads/SBET) of
many low surface area oxides
(SiO2, Al2O3, ZrO2, TiO2, MgO)
form a universal curve
t = f(p/p0). Vads
t
 Numerical description of that S BET
curve
 valid for t = 3.5 – 10 Å
 and p/p0 = 0.1 – 0.8
13.99
Harkins t ( )1/ 2
p
Jura  log  0.034
p0
1/ 3
 
 5 
Halsey t  3,54 
 ln  P0  
  P  
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Isotherm expressed as Vads = f(t)
180 180
160 160
140 limit of validity 140 Vads
S
Vads (cm3 STP/g)

Vads cm3 STP/g


120 120
100 100
t
65811
80 80

60 60

40 40

20 20

0 0
0 0.2 0.4 0.6 0.8 1 0 0.5 1 1.5 2
p/p0 t (nm)

t = thickness that the


13.99
t ( )1/ 2 adsorbed phase
p
 log  0.034 would have on a non-
p0 porous material
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Interpretation of t-plots

 Adsorption on a non-porous solid

Vads

Vads
S
p/p0 t

 Adsorption on a microporous solid t


 micropores rapidly filled
 then adsorption on external surface Vads
Vads
S
t
Vmicro
p/p0
© IFP

t
Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Example zeolite NaX
Isotherm t-plot
300 300
limit of validity
250 250
Volume adsorbé (cm3 STP/g)

Vads (cm3 STP/g)


200 200

150 150

100 100

50
Vmicro
50

0 0
0 0.2 0.4 0.6 0.8 1
0.00 0.50 1.00 1.50 2.00
p/p0
t (nm)

 Vmicro = 175 cm3 STP/g = 0.271 ml/g


© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Surface BET vs. Surface t-plot
70 St  plot (m2 g 1 )  1.5468  PNT (cm3 g 1nm1 )
60

50
Vads cm3 STP/g

40 y = 45.11x + 0.965

30

20
y = 21.83x + 0.685
10

0
0 0.2 0.4 0.6 0.8 1 1.2
t (nm)

SBET (m2/g) St-plot


65811 74 70
70299 36 34
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Micro- and mesoporous solids
Solide non poreux : Adsorption sur Surface Externe - Courbe t = t(P/P0)

------------------ Pressions relatives P/P0 croissantes ------------------

µP mP

Solide poreux : Adsorption sur Surfaces Interne & Externe


Vol. Ads.

Solide Non Poreux Solide µ et m-Poreux


Vol. Ads.

Adsorption sur
Surface Externe
Adsorption sur Surface Externe VµP+VmP

VµP Adsorption sur


Surface Interne
(mésopores) => SB.E.T.
© IFP

Epaisseur t Epaisseur t

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
t-plot may be ambiguous

 Zeolite beta with supermicropores


 no clear linear region in the t-plot
 micropore volume depends on the interval chosen for
extrapolation
200
220 H-MCBeta
H-MCBeta steamed and leached

Vads cm /g STP
N2 ads cm /g STP

180
200

160

3
3

180
140

160
120

140 100
0.0 0.2 0.4 0.6 0.8 1.0 0 2 4 6 8 10 12

p/p0 Thickness - Halsey [A]


© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Méthode de Dubinin-Radushkevitch

 Théorie de remplissage de micropores


Vads
Vmicro

 exp   A / E 
2
 A   RT ln
p
pvap, sat

 A = potentiel d'adsorption
 E = energie caracteristique de l'adsorbant
 pvap,sat = pression vapeur saturante
 pvap,sat = p0 (1 atm) pour N2 à 77K
 Equation linéarisée
2
 p
log Vads  log Vmicro  D log 
 p0 
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Example NaX
-0.51
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
-0.52
Pression élevée
-0.53 Potentiel faible
Adsorption dans mesopores
log10(Vol. liquide)

-0.54

-0.55 Pression faible


Potentiel élevée
-0.56
Début remplissage
micropores
-0.57

-0.58
log10(P/Po)^2

Vmicro = 10^(-0.548) = 0.283 ml/g


© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Average pore size

 From geometrical rules


 Cylindrical Pore:
V = pore volume V  r h
2
 S 2 4V
  d
 S = surface S  2rh V r S

 Spherical pore :  4 3
V  r S 3 6V
 3  d
S  4r 2 V r S

 Slit pore :
V  lhr S 1 2V
  d
S  lh V r S
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Micropore size distribution
 Based on the relation between adsorption strength
and the ratio between adsorbate size and pore size
 condition: adsorption controled by van der Waals interactions
only
 Mathematical models
 Horvath-Kawazoe
 Saito-Foley Saito, Foley, AICHE Journal 1991, 37, 429.

 Input parameters
 pore geometry: slit-shaped, cylinder, sphere
 parameters of the Lennard-Jones potential well: ε and σij
 not always well known for atypical solids
 Input data
high precision isotherm at very low pressure !
© IFP


Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Mesopore size distribution

 via the theory of capillary condensation


A B C->D

D
Vads
E

D->E E->F F->A


F

C
A
B

P/P0
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Kelvin equation
 Capillary forces in the pore lower the vapor pressure of the
condensed liquid
 Pressure at the concave side of an interface is higher than at the
convex side.
overpressure counteracts the surface tension, which tries to collapse the

interface area
 Application to interface between gas and pressure
adsorbed liquid film in a pore difference Δp
 Pressure in liquid is lower than gas pressure.
 Means that chemical potential in liquid surface
is lower, in other words tension
that the equilibrium vapor pressure
is lower.
 Capillary condensation occurs at a
lower pressure than condensation.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Kelvin equation
 Mathematic formalism
 Work against interfacial tension = change in free enthalpy

dA    dncapliq  capliq  dnliq  liq


dVcapliq
dncapliq  dnliq 
Vm
dVcapliq dVcapliq
dA      capliq   liq
Vm Vm
dA Vm    dVcapliq (  capliq  liq ) pcap
capliq  0  RT ln
dA p0
( capliq  liq )   Vm   psat
dVcap, g liq  0  RT ln
p0
pcap dA
RT ln  Vm
© IFP

psat dVcap, g
Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Generalized Kelvin equation

pcap V  dA dV/dA = change in volume per change in


ln  m  interface area
psat RT dV Vm = molar volume

 dV/dA depends on the pore geometry


 Sphere: dV/dA = r/2
 Cylinder dV/dA = r
 Slit dV/dA = d distance between slits
 Relation to curvature of the pore
pcap 2 Vm
ln 
psat r RT
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Capillary condensation

p p

2rp
vl
pc  p0  exp( )
RT  (rp  t )
2(rp – t)
 Multilayer adsorption on the surface: layer thickness t 
 At a certain effective pore radius r p-t, capillary condensation
occurs
© IFP

 Pore is filled. There is a step in the adsorption isotherm.


Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Example isotherm with hysteresis

Type H1
1200
capillary
1000
evaporation
Volads (cm3/g)

800

600

400
capillary condensation
200

0
0.0 0.2 0.4 0.6 0.8 1.0
p/p0
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Application of the Kelvin equation
1200

1000
Harkins-Jura t-plot equation
Volads (cm3/g)

800

p 13.99
600 log  0.034  2
p0 t
p  4.14
400 log 
p0 rm

200 Kelvin equation

0
0.0 0.2 0.4 0.6 0.8 1.0 Cylinder !!
p/p0

pc/p0 t(pc) rm rc rp= rc+t


Ads 0.671 8.2 23.9 11.9 20.2
Des 0.565 7.0 16.7 16.7 23.7
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Origin of the hysteresis

Adsorption: Desorption:
cylindrical meniscus hemisperical meniscus
rm = 2reff rm = reff

 The lower rm, the lower is the pressure of capillary


condensation/evaporation
 rm,ads = 2 rm,des
 Capillary evaporation at a lower pressure than capillary
condensation
 Hysteresis loop in isotherm
© IFP

L.H. Cohan, JACS 60 (1938)


Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Origin of the hysteresis
 Truth seems to be complicated than in Cohan’s
theory.
 Cohan’s theory, based on the shape of mensiscus,
theory suggests that different vapor liquid equilibria
exist in a spherical or cylindrical meniscus, which
leads to hysteresis => is an equilibrium picture.
 Molecular simulations (DFT) suggest that adsorption
branch is not in thermodynamic equilibrium, but is a
metastable state.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Hysteresis and metastability

 Adsorption branch comprises metastable states.


 grand free energy of filled pore is lower than
that of the empty pore in the hysteresis
 Desorption branch is in thermodynamic
equilibrium.
 Reason for metastability in adsorption branch
 barrier of nucleation by formation of a liquid
bridge across the pore
 condensation occurs when limit of
metastability is reached
limit of metastability
gas phase

liquid phase
© IFP

Monson, MMM, 2012.


Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Peterson, JCS Farad 2, 1986.
Pore size distribution – BJH model
p0 – all pores filled

p1 – capillary evaporation in
largest pore
reduction of layer thickness

p2 – capillary evaporation in
2nd largest pore
reduction of layer thickness in
both pores
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Mathematic description - BJH

rp rp
(rk  t ) 2
V  2
V p
rk rk + t rp

 In each desorption step pn-1  pn, capillary


evaporation occurs from a pore of size r pn
 The volume desorbed in that step (Vn) can be related
to the pore volume by the geometrical relation given
above.
 Complication: Reduction in layer thickness in the
pores, which were already emptied, also contributes to
© IFP

Vn. of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Characterization
The BJH-equation E.P. Barret, L.G. Joyner, P.P. Halenda, JACS 73 (1951) 373.

rpj  t j
2 2 n 1
rpn rpn 2V p
V pn   Vn   tn   Apj Ap 
(rkn  tn ) 2
(rkn  tn ) 2
j 1 rpj rp

capillary evaporation correction term


 For each desorption step the average diameter of the pore, which
undergoes capillary evaporation is calculated from the Kelvin
equation and the t-plot equation: rp = rk + t
p 13.99
log  0.034  2
p0 t
p  4.14
log 
p0 rk
 t is the change in layer thickness in each desorption step
 Vn is the volume desorbed in each step
 A plot of pore volume vs. pore radius is obtained.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Limits of the BJH-model

 BJH underestimates the pore size below 7.5 nm –


WHY?
 Model separates the adsorbed film and the capillary
condensate – not a realistic picture
 Fluid-wall interactions are neglected
 Kelvin equation may not be valid in very narrow pores
 Surface tension might increase with curvature

 One should speak about a BJH-value rather than pore


diameter.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Adsorption vs. desorption isotherm

 Adsorption isotherm – Pros and cons


 Condensation can be delayed – not in thermodynamic
equilibrium
 Cylindrical meniscus not stable – not advisable to use
adsorption branch for cylindrical pores
 Desorption branch – Pros and cons
 Evaporation not delayed, in thermodynamic equilibrium –
generally preferred
 Ink-bottle type pores: smallest openening determines the
desorption
 In a network of interconnected pores: percolation (transport)
effects determine the desorption
Conclusion: look at both and compare them
© IFP


Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Ink-bottle pores

rdes
rdes

 In desorption smallest pore openening is determining.


 When pressure of capillary evaporation is reached for
the smallest pore opening, whole pore is suddenly
emptied.
 Sudden drop in the desorption isotherm – type H2
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Pore network effects

N2 can only desorb from


pore A and C if pore B has
already been emptied.
Seaton, CES, 1991.

 Desorption branch of isotherm not in thermodynamic


equilibrium any more, due to pore blocking.
 Desorption is controlled by a percolation process.
 probability that the pore is connected to the outer surface.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Ink-bottle type hysteresis loop

Type H2

Adsorbent:
xerogel and
alcogel

• If pore size distribution is calculated from desorption branch, an


artificially narrow pore size distribution is obtained
• The adsorption branch has to be used to calculate the pore size
distribution.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Ink-bottle type hysteresis loop
Isotherm Pore size distribution
1000

900 6

800 Desorption
Vol ads (cm /g)

700
3

dV/dD
600

500 Adsorption

400 0

300
0 2 4 6 8
200
Pore size (nm)
100
0.0 0.2 0.4 0.6 0.8 1.0
p/p0

 Narrow peak in the pore size distribution of the adsorption branch


is an artefact caused by the forced closure of the isotherm at p/p 0
= 0.43
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Isotherm without hysteresis loop

700

600
Vol ads (cm /g)

500
3

400

300

200

100
0.0 0.2 0.4 0.6 0.8 1.0
p/p0
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Closure point of the isotherm

 Critical temperature inside a pore lower than in bulk:


Tc,pore < Tc, bulk
 Tc,pore decreases with decreasing pore diameter1
 Above Tc no condensate-vapour meniscus
 For narrow pores Tc < 77 K
 No hysteresis for filling and emptying of these pores
 Tc,pore < 77 K for pores, which show capillary
condensation at p/p0 = 0.4  closure point of isotherm
 Tc(Ar) > Tc(N2)  Ar isotherms show hysteresis when
N2 isotherms don‘t2
© IFP

[1] R. Evans, J. Phys. Condensed Matter 2 (1990) 8989.


[2] M. Thommes, R. Köhn, M. Fröba, J.Phys.Chem. B 104 (2000) 7932.
Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Summary – Pore size analysis
 BJH relies on Kelvin equation and the universal t-
curve for determining pore size distribution
 both concepts have weaknesses
 Alternative methods exist
 Density Functional Theory
 Broekhoff de Boer – improvement of BJH
 Derjaguin – concept of disjoining pressure
 surface tension concept that takes interaction with solid into
account
 Every model assumes a certain pore geometry
(cylindrical, spherical, slit-shaped) – influences the
results!!
 The adsorption and the desorption branch contain
© IFP

different
Characterization information
of porous solids – look
- Characterization of Catalysts at both.
and Surfaces - G Pirngruber 25 october 2016
Literature
 Textbooks on adsorption
 D.A. Ruthven, Principles of Adsorption and Adsorption
Processes, Wiley
 D.D. Do, Adsorption Analysis: Equilibria and Kinetics,
Imperial College Press
 Review articles discussing a more moderne view of
hysteresis effects and pore size analysis
 P.A. Monson, Understanding adsorption/desorption
hysteresis for fluids in mesoporous materials using simple
molecular models and classical density functional theory,
Microporous Mesoporous Materials 160 (2012) 47.
 B. Coasne et al., Adsorption, intrusion and freezing in porous
silica: the view from the nanoscale, Chem. Soc. Rev. 42
(2013) 4141.
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
CO2 maîtrisé | Carburants diversifiés | Véhicules économes | Raffinage propre | Réserves prolongées

Alternative explanations of the


Kelvin equation
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Physical principles of capillary
condensation

 Vapour pressure of a liquid under external pressure


External pressure  Equlibrium gas-liquid g = l
H2O vapour pressure pg,0
H2O H2O
g  Equilibrium disturbed by external
H2O pressure
dg = dl  vl dpl = vg dpg
H2O l vl dpl = RT/pg * dpg
vlpl = RT ln (pg/pg0)
pg = new vapour pressure of H2O
pg0 = orignial vapour pressure of H2O vl
pl = external pressure on liquid (H2O), pg  pg ,0  exp( pl )
which induces the change in vapour RT
pressure
© IFP

vl = molar volume of liquid (H2O)


Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Vapour pressure of a drop of liquid

2
pin  pout 
r pout pin - pout

 = surface tension overpressure


r = radius of curvature inside the drop
holds against
2 pin
pl  pout 
r
surface tension
2 tries to contract
pl  p  the drop
r
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Vapour pressure of a liquid void

2 pout pin - pout


pin  pout 
r overpressure
in gas bubble
holds against
2
pin  pl  pin
r

2 surface tension
pl  p  tries to collapse
r the void

2vl
pg  pg ,0  exp( ) Kelvin equation: vapour pressure
RT  r inside a void is lowered
© IFP

Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016
Kelvin equation
 Capillary forces in the pore lower the vapor pressure
of the condensed liquid
pcap Evaporate n moles from capillary
p
cap  0  RT ln liq  0  RT ln sat and condense on a flat surface.
p0 p0
You have to work against the
surface tension γ.
n

The work done to overcome the surface


tension is equal to change in chemical
Δl θ potential.
V psat
n    RT ln
Vm pcap
Explains why vapor pressure of the
unconfined liquid is higher than in the
© IFP

capillary unconfined liquid


capillary.
Characterization of porous solids - Characterization of Catalysts and Surfaces - G Pirngruber 25 october 2016

You might also like