Untitled
Untitled
FOUNDING EDITOR
MARTIN J. SCHICK
Consultant
New York, New York
SENIOR ADVISOR
ARTHUR T. HUBBARD
Department of Chemistry
University of Cincinnati
Cincinnati, Ohio
ADVISORY BOARD
P. SOMASUNDARAN
Henry School of Mines
Columbia University
New York, New York
7 . Anionic Surfactants (in two parts), edited by Wamer M. Linfield (see Volume
Part I I Adsorption from the Vapor/Gas Phase onto Nanoparticles and Porous
Solids
13. Characterization of Energetically Heterogeneous Surfaces from Experimental
Adsorption Isotherms 295
James P. Olivier
14. Computer Simulations of the Structural and Thermodynamic Properties of
Adsorbed Phases 319
William Stede
15. Surface Heterogeneity Effects o n Adsorption Equilibria and Kinetics:
Rationalizations of the Elovich Equation 355
WladyslawRudzinski and Tornasz Panczyk
16. Single- and Multicomponent Adsorption Equilibria of Hydrocarbons on
Activated Carbon: The Role of Micropore Size Distribution 391
K. and D. Do
17. Surface and Structural Properties of Modified Porous Silicas 443
Michal Kruk and Mietek Jaroniec
18. Nanodimensional Magnetic Assembly of Confined O2 473
Katsumi Kaneko
19. Heat of Adsorption of Pure Gas and Multicomponent Gas Mixtures on
Microporous Adsorbents 501
Shivaji Sircar and Madhukar B. Rao
art 111 Adsorption from the Liquid Phase onto Nanoparticles and Porous Solids
20. Surface Chemistry of Activated Carbon Materials: State of the Art and
Implications for Adsorption 529
Ljubisa R. Radovic
21. Charge Regulation at the Surface of a Porous Solid 567
Boris V. Zhrnud and Lennart Bergstrom
26. Hydrous Metal Oxides as Adsorbents for Aqueous Heavy Metals 675
Russell Crawford, David E. Maimvaring, and Ian H . Hording
lan H. Harding Centre for Applied Colloid and BioColloid Science, Swinburne
University of Technology, Melbourne, Australia
. Impens
Chemistry Laboratory of Adsorption and Catalysis,
University of Antwerp, Antwerp, Belgium
xia Liu Cominco Research. Cominco Ltd., Trail, British Columbia, Canada
Shivaji Sircar Corporate Science and Technology Center, Air Products and
Chemicals, Inc., Allentown, Pennsylvania
I. Introduction
II. Experimental
III. Results and Discussion
A. Liquid-phase modification of silica gel with tetrachlorosilane
and the first ammoniation
B. Second and higher reaction cycles
C. Porosity study
IV. Conclusions
References
I. INTRODUCTION
Due to the widespread use of ceramic coatings, several synthesis techniques have
been developed in recent decades. The majority of these techniques, such as che-
mical vapor deposition (CVD) and physical vapor deposition (PVD) and their
variants, are focused on the synthesis of flat coatings. Recently, the preceramic
polymeric synthesis route has offered the possibility to impregnate preceramic
materials into porous matrices prior to pyrolysis in order to create coated or
composite materials This technique, however, leads to pore filling and altera-
tions of the original substrate texture. The impregnation does not guarantee a
chemical bond between the substrate and the preceramic material. This can lead
to an unstable composite material after pyrolysis. T o avoid this. another technique,
called chemical surface coating (CSC) can be used. Based on the successive
chemisorption of gaseous reagents in a cyclic way at relatively low temperatures,
the technique offers clarity in the reaction mechanisms [4,5] and the evolution of
porosity [6] during the preceramic synthesis. It has been shown that the preceramic
coating is chemically bound to the substrate surface and that is has grown in a very
homogeneous way in spite of the very irregular substrate texture. The so-called
coating precursor is afterwards pyrolyzed, and a chemically bound ceramic coating,
e.g., or BN; is formed on silica gel. Characterization of the thermally con-
verted Si- N or B-N polymer shows that a very thin homogeneous silicon
oxynitride or boron nitride coating, respectively, is formed, whereas the original
substrate texture is practically unchanged [7,8].
In the procedure described in this chapter, a preceramic polysilazane coating is
synthesized on silica gel via CSC using S i c &and NH3, but the chlorosilylations are
performed in the liquid phase. In this way, a very small amount of the hazardous
reagent is used in comparison with gas-phase reactions, and the total precursor
synthesis can be performed at room temperature.
The solid-gas reaction of with thermally pretreated silica gel has been
optimized, and the reaction mechanisms were studied earlier [4]; a reaction tem-
perature of 663 K is needed to modify all the surface hydroxyl groups. T o reduce
the reaction temperature, and Hair [9] studied the effect of the presence of
during reaction with At room temperature, a
very reaction yield was obtained. The reaction mechanism was deduced from
the of silica proposed by Blitz et [see reaction
Modification of
Ammoniation
The creation of chlorosilyl groups on the substrate surface is very important in the
synthesis of a chemically bound preceramic polysilazane coating, because the
amounts and types of chlorosilyl species present influence the reaction behavior
of the reactions. In Table 1 , a comparison of the elemental analysis data is shown
for the gas- and liquid-phase synthesis routes as a function of the number of
reaction cycles. Infrared spectroscopy of the modified silica gel was used to study
the reaction yield. In Figure the infrared spectra of pure silica (973 spectrum
a) and silica modified with in both the gas and liquid phases (spectra are
shown. When toluene (spectrum c) or cyclohexane were used as solvents, all the
huge bands between 3000 and 2000 c m ' and the sharp bands superimposed on the
Si 0 vibrations between 1500 and 800 exactly correspond to
chloride, which was formed on the silica surface and could not be removed by
severe washing. When dichloromethane was used (lines d and e) no salt formation
was detected, and very fast removal of the solvent was possible due to its volatility.
This solvent's high Lewis accepting property is shown to promote the reaction of
silica gel with octadecyldimethylchlorosilane a t room temperature The nitro-
gen-containing base used in this study has to be a tertiary amine, avoiding any
polymerization reactions with Triethylamine was chosen, as its conjugated
base has a relatively high value (10.8). When the base was added in a catalytic
amount (5 1/20 relative to the OH groups present), the reaction amount
was very low (Fig. 1, spectrum e). The reason is that after protonation the base has
no further influence on the reaction, which implies that it does not act as a catalyst,
as was erroneously reported in the literature An excess of base (800 resulted
in almost complete modification of the surface hydroxyl groups (Fig. 1, spectrum
d). As a blank, the reaction was also carried out under similar conditions but
Liquid phase route Cl 1.53 (0.05) 2.7 (0.1) 1.6 (0.1) 1.6 (0.3)
N
In Fig. 2, the spectra of the samples after the first (spectrum a) and
the second liquid chlorosilylation are shown for different reaction conditions.
the second chlorosilylation on, dichloromethane could not be used and was
replaced by cyclohexane, because dichloromethane reacts with the amino groups.
The addition of to the reaction mixture (spectrum b) leads to the
deposition of triethylammonium chloride a t the sample surface, which could not be
removed with the solvent used or by sublimation, because thermal treatment also
leads to condensation reactions of the amino groups Therefore, the
are performed without the addition of triethylamine. Spectrum c of Fig. 2 is
that of the preceramic precursor after the second chlorosilylation. In the Si-
region, maxima are observed at 624 and This means that both
and groups are formed. The total disappearance of the asymmetric N-H
bending vibration at assigned to groups indicates a total
modification of these groups. Therefore, an estimation of the ratio of mono
and bidentate reactions could be calculated Elemental analysis and iteration
showed that of the groups are modified by bidentate binding and
by monodentate binding.
nitrogen, chlorine, and weight uptake during the four reaction cycles for both the
liquid- and gas-phase synthesis routes is shown in Fig. 3.. In general, all the curves
are linear except for the chlorine uptake in the liquid-phase route. where a deflec-
tion is observed at two cycles. However, this trend is not followed by the
uptake, as the ammoniation is not disturbed by solvent hindrances. From this
figure, it is clear that in the liquid case a smaller amount of polymer is formed
on the silica surface, whereas all the gas-phase slopes are higher than in the case of
liquid-phase modification. The total increase in weight gain after four reaction
cycles if 0.55 and 0.8 g per gram of unmodified silica in the liquid and gas synthesis
routes, respectively.
From the rough analytical data of Table 1, the amounts of
and can be calculated (Table 2). These data can give an idea of the
cross-linking in the polymeric preceramic material. The amount of reflects
roughly the ratio between the increasing weight of the two synthesis routes (weight:
= 0.69; whereas the amount of unreacted
groups is relatively lower in the liquid-phase synthesis route, even when the
increased weight ratio is taken into account = 0.47). This result was
expected from the trend difference in the and uptakes. Another important
factor is the percentage of primary amines that are not modified into In
the gas phase, only 26% of the total population has not reacted,
whereas 36% is unreacted in the liquid-phase synthesis route. This means that a
lower cross-linking efficiency is reached in the liquid-phase route. This value, in
combination with the lower uptake of nitrogen in the liquid-phase synthesis route,
shows that a gas-phase synthesis is more suitable for the synthesis of a dense cross-
linked precursor.
With each modification cycle during the synthesis of the structure, the
amount of product left at the silica pore walls increases. To study the effect on the
substrate porosity of both types of precursors, adsorption-desorption isotherms
are recorded at 77 K as a function of the number of reaction cycles. The isotherms
presented in Fig. 4 are from the liquid-phase synthesis route and are all of type IV
according to the (BDDT) classification
Therefore, both the substrate and the coated material are mesoporous materials.
The total pore volume decreases with increasing number of reaction cycles. The
type of hysteresis changes from type to according to the IUPAC classifica-
tion Type corresponds to the pore shape of pure silica gel, whereas the
type indicates ink-bottle or irregularly shaped pores. The evolution in the hysteresis
shows the occurrence of pore shape modification during the formation of the
polymer. This phenomenon earlier observed was for the gas-phase synthesis
route The behavior of the constants for the liquid-phase reaction shows
the same trend as in the gas phase; from the first cycle on, an important
decrease is observed, indicating total monolayer coverage.
Using the calculation methods pointed out by et al. a difference
between the theoretically expected and experimentally observed pore volume
decreases can be used to calculate the amounts of pore blocking and regular
reaction cycle
A similar trend was observed when the decreases in pore length were calculated
using integration of the pore volume distributions obtained by the BJH method
as explained earlier However, this method is valid exclusively for cylind-
rical pores. The results are shown in Fig. 6. A systematically pore length
value is obtained for the liquid-phase route. These data show that more pore
blocking occurs in this synthesis route.
The calculated amounts of pore-blocking and pore narrowing and their contri-
bution to the total pore volume decrease are listed in Table 3. The results are based
o n the data of Fig. 5 . Systematically higher absolute and relative amounts of pore
blocking are seen for the liquid-phase synthesis route in spite of the lower uptake of
N and per reaction cycle, resulting in a lower coating homogeneity. Moreover,
from the higher absolute and relative narrowing values in the gas-phase synthesis
route, it can be concluded that a higher coating efficiency is obtained using the
phase synthesis route.
Successive reactions of in the liquid phase and in the gas phase were
performed on silica gel. The reaction parameters for the
solvent choice and base promotion) were optimized to get as much uptake
as possible. The chemical composition of the polysilazane layer as a function of the
number of reaction cycles is obtained by using a combination of elemental and
infrared analysis.
A comparison of these results with the composition of a polysilazane coating
fully synthesized in the gas phase showed less uptake of N and and a greater
number of groups due to the presence of the solvent. This implies that the
(gas phase)
0.70 0.00 0.00
0.65 0.12 0.05
0.58 0.17 0.12
0.49 0.20 0.21
0.39 0.26 0.31
(liquid phase)
0.70 0.00 0.00
0.66 0.13 0.04
0.61
0.57 0.26 0.13
0.53 0.33 0.17
PV = pore volume.
"The experimental errors on the pore volumes range from 0.0005 to The errors in the
percentages are smaller than 1
I. Introduction
I I . Experimental Methods
A. Synthesis of aluminum-free clays
B. Synthesis of H-magadiite and H-kenyaite
C. Silica-intercalated layered silicate
D. Analytical methods
III. Results and discussion
A. Synthesis of aluminum-free clays
B. Octylamine/octylammonium layered silicate gel
C. Reactions with TEOS
D. Gelation effects
IV. Conclusions
References
of and
H-magadiite was prepared by titrating with 0.1 using the
method of Beneke and Lagaly A suspension of 40g of per
of deonized water was stirred for 1 h. The suspension was then
with 0.1 N to a final pH of 2.0 and then maintained at the same pH value for 1
week in a refrigerator. H-magadiite was recovered by filtering, washed with deio-
nized until chlorine free, and then dried in air at The preparation of H-
kenyaite was the same as that of
. Analytical Methods
Identification of samples was carried out by X-ray poweder diffraction (Rigaku
diffractometers equipped with radiation) and scanning electron
microscopes (Jeol, The chemical compositions of synthetic magadiite
and kenyaite were determined with an energy-dispersive X-ray spectrometer (Link
and analysis (TGA). EDX analysis was
carried out using an electron beam with the sample of carbon-coated pellets. For the
quantitative analysis of silica and sodium, quartz and sodium chloride were used as
references. analyses were performed in air on a Dupont 9900
analyzer. Basal spacings of samples were determined from the
001 X-ray powder diffraction using a equipped with Cu
radiation. Samples of layered silicate gel and TEOS-
intercalated layered silicate gels were prepared by smearing a thin film across a
microscope slide and then recording diffraction patterns of wet samples. The pro-
ducts gelled by EtOH suspension. base-catalyzed reaction. and acid-catalyzed reac-
tion were prepared by depositing a gelled suspension on a glass plate and allowing
the suspension to air dry at The calcined silica-pillared magadiite and kenyaite
for X-ray diffraction analysis were prepared by calcining the gelled products on glass
plates at in air for 4 Nitrogen adsorption-desorption isotherms were
determined by ASAP 2000 at All samples were outgassed at
under a vacuum for 4 Surface area was determined using the
Emmett-Teller (BET) equation. The pore size distributions of silica-pillared pro-
ducts were determined by using the (BJH) equation
time because pure kenyaite can be synthesized within only a short range of reaction
time.
Figure 3 shows XRD of synthetic samples obtained at in the
higher molar ratio of + = 20 and
= with reaction time. can be synthesized directly from
amorphous silica, and magadiite was not observed as a n intermediate.
The higher molar ratio accelerates the formation of kenyaite. Silica is first
dissolved in alkali solution a t elelvated temperature, and then the nucleation of
kenyaite occurs gradually. Therefore, the ratio in the solution is a very
important factor in the formation of kenyaite. Also, the higher content of silica
in the composition of starting materials may cause a large ratio in the
solution. These conditions could accelerate the nucleation of kenyaite that has a
Temperature
5 Thermogravimetric analyses of (a) and (b) Na-kenyaite.
Ref, 30.)
The unit cell of Na-kenyaite deduced from TGA and the analyses of
sodium and silica is which compares favorably with the
approximate composition supported by the early work of
Lagaly and Beneke
E D X analysis.
Source: Ref. 30.
H-magadiite (air-dried)
H-kenyaite (air-dried)
magadiite gel
kenyaite gel
TEOS-intercalated magadiite
TEOS-intercalated kenyaite
kenyaite gel, the gallery height is 2.15 indicating that this value is similar t o the
gallery height of H-magadiite. This result indicates that the expan-
sion of the gallery height of the layered phase is caused by the formation of
of ion and molecules between the silicate
layers.
Reactions with
The reactions of magadiite gel with T E O S afford
TEOS-intercalated magadiite with well-ordered basal spacings. This product exhi-
bits reflections correspoding to a basal spacing of 3.80 nm (gallery height, 2.68 nm).
The gallery height for the siloxane intercalated derivative is about larger
than the gallery height of observed in the
magadiite precursor gel. the case of TEOS- intercalated kenyaite, the basal
spacing of this product is (gallery height, The gallery height
for the siloxane intercalated derivative is 0.39 n m larger than the gallery height of
2.15 n m observed in magadiite precursor gel. Kwon
et reported the effects of organic solvents o n the intercalation of octylamine
into According to this result, the expansion of gallery height could be
effected not only by the solvation of intercalated octylamine but also by the steric
effects of solvated solvent molecules and solvent itself. In the
layered silicate gel reacting with excess TEOS, T E O S could be exchanged
with solvated free octylamine molecules in the interlayers a n d hydrogen-bonded
ion on the inner surface. The interlayer and structure of
TEOS could cause a n increase in gallery height.
reflections corresponding to basal spacings of 2.50, 4.02 and 3.42 for magadiite
and 3.09. 4.23, and for kenyaite. The silica-pillared products gelled by
base- and acid-catalyzed reactions exhibit a large increase in the gallery height,
compared with the gallery height of the product gelled by E t O H sus-
pension. A distinctive increase in the gallery height is related to the size and struc-
ture of pillared silica, which could be associated with the amount of intercalated
TEOS, the condition (catalyst type, solvent composition, p H , etc.), and the
rate of gelation. The calcination of the intercalated derivatives a t
yields silica-intercalated products, showing a 1-3 A reduction of basal spacings
(Table 5).
by base and acid catalysts could minimize the release of TEOS outside
the layered phase during the of the intercalated TEOS, because
time is markedly reduced. The hydrolytic polycondensation of intercalated TEOS
by acid o r base catalyst could form silica clusters of a highly branched and stiff
network structure. These conditions could bring about the effective of
intercalated TEOS in the interlayer and contribute to develop the size and structure
of pillared silica clusters that bring about the the large expansion of the gallery
height of pillared kenyaite.
Adsorption- desorption isotherms of nitrogen were obtained a t Several
pore characteristics calculated from them are listed in Table 5. The specific surface
areas were calculated by using the BET equation from the adsorption isotherm
below = 0.1. The calcined silica-pillared products have high surface areas of
533-845 depending on the catalyst type. As the gallery height of the calcined
sample decreases, the specific surface area is increased owing to the development of
micropores. The total pore volume of the base-catalyzed sample is the largest
among the three pillared samples, indicating that the average pore size of that
is the largest.
Figure 6 shows the pore size distributions of silica-pillared The sam-
ple treated with the base-catalyzed has the greatest microporosity (dia-
meter 2 and shows a sudden increase in volume near 3.6 nm with a
narrower pore size distribution. These results indicate that abrupt by base
catalyst is closely related to the formation of a more uniform pore. O n the other
hand, the sample gelled by the suspension of E t O H has a broader pore size dis-
tribution. In the case of by a n acid catalyst, decreases and
mesoporosity increases between 4.0 and nm compared with with no
catalyst. It is interesting that acid- and base-catalyzed products are very similar in
physical properties to surfactant-templated mesoporous silica Figure 7
shows the pore size distributions of silica-pillared The base-catalyzed
sample shows a sudden increase in volume near 2.2 with a
tively narrow pore size distribution. The acid-catalyzed sanple shows a narrow pore
size distribution with a sharp peak near 1.8 These values are similar to the
gallery heights of the samples. These results indicate that the pore size of base- or
acid catalyzed samples is closely related to the gallery height. The abrupt of
acid and base catalysts is closely related to the formation of uniform pores.
On the other hand, the sample gelled by the suspension of EtOH has a broader pore
size distribution and more (diameter 2.0 nm). Therefore, the
or base-catalyzed reaction not only allows for high surface area through optimiza-
tion of silica pillaring but also permits the synthesis of materials with narrow pore
size distributions.
I. Introduction
I I . Molecular Self-Assembly
III. Direct Silanization
IV. Silica Coatings
V. Applications
A. Self-assembled monolayers
B. Silanized monolayers
VI. Summary
References
Magnetic particles have shown great promise as a base material for engineering
carriers in separation science and technology. A potential area of application of
magnetic particles is as supports for catalysts. Magnetic particles used in these areas
are called magnetic carriers, although a distinction is sometimes made between
magnetic carriers and magnetic tags. The former applies to magnetic particles
significantly larger than the targets so that the targets coat on the particles (carriers)
and the latter refers to magnetic particles smaller than the targets so they (the tags)
deposit on the targets In this chapter, both are referred to collectively as
magnetic carriers.
The objective of magnetic carrier technology is to confer the magnetic property
to a naturally nonmagnetic target so that the target can be separated from the
stream using magnetic separators. This is particularly important for separations
Magnet
, ,
-287 -283 -279 -170 -166 -162 -158
Binding Energy Binding Energy
.2 Narrow-scan XPS spectra for the elements of interest on with
assembled (a) MHA and (b) acid, in comparison with (c) untreated
particles.
A
in water
in toluene
.5
Zeta potential of without silanization and with silanization by
APTES in water and in toluene.
XPS
DRIFTS Fe in
Before leaching leaching
60.3 4.5 --
two
sol-gel
DL
Coating method
60.3
Dense liquid 2.8
Sol-gel 1.1
Two-step nil
nil
(a)
-955 -950 -945 -940 -935 -930 -925 -375 -370 -365
Energy Binding Energy
Narrow-scan XPS spectra for
8 elements of interest on with various
treatments: (a) Non; (b) MHA coating; MHA coating by loading; and
(d) MHA coating by Ag
It has been found that reactive amine groups are capable of attracting and adsorb-
ing such as arsenate, dichromate, tungstate, and
at the microgram per liter concentration level These species can then
be effectively from an effluent by using magnetic separation magnetic
carriers to detoxify the effluent. magnetic carriers with reactive amine groups
are found to confer specific adsorption for heavy metal ions, including
mercury, lead, copper, zinc, manganese, and monovalent silver in contaminated
aqueous and nonaqueous solutions Copper loading the silanized magnetic
particles of reactive groups is shown Fig. 9. For comparison, copper
loading on silica-coated magnetic particles is also shown in this figure. It is evident
that at a pH copper an be removed effectively from concentra-
tion solutions 100% removal from a With increasing initial
copper ion concentrations, the copper removal efficiency decreases as copper load-
ing approaches the capacity limit of about of per gram of particles. In
contrast, the copper loading on the silica-coated magnetic particles is significantly
lower than on the silanized particles, suggesting the important role of reactive
amine groups in application. It should be noted that the metal ion loading is
pH-dependent. a recent communication we demonstrated the selective
loading of various metal ions such as copper zinc by varying the solution
at which loading is performed.
It is evident that the of magnetic carriers in separation science and technol-
ogy, whether in the area of biological cell or industrial effluent detox-
ification, offers great flexibility. Other areas of potential applications of magnetic
carriers with reactive amine groups include immobilization of enzymes, drug deliv-
ery, and biological sensors, among others, A detailed review of these applications
can be found elsewhere The application can be extended further to other areas
as long as a key-lock relationship (see Fig. 10) can be developed. The lock varies
from antigen or as in biological applications to metals or toxic species as
in environmental applications, the key could be an antibody or specific
The important consideration for these applications is to further
improve the stability of coatings, which significantly reduce the operating
90
0.15
magnetic carriers
80
70
0.10
60
a
50
0.05
silica coated
30
0.00
0.0 0.1 0.2 0.3 0.4
Equilibrium concentration of
FIG. 9 Copper loading on left-hand axis) and removal by right-hand axis)
APTES on silica-coated particles in comparison to silica-coated
particles left-hand axis) (loading at pH 5.3).
Three basic requirements for magnetic carriers in separation technology are fine
particle size, and chemical stability in addition to tailored
functionality. Two different approaches for engineering nanosize
carriers, molecular self-assembly and silica coating followed by
have been described. Using molecular self-assembly, a bolaamphiphile of two dif-
ferent functional groups (16-mercaptohexadecanoic acid) was found to form a
packed and stable o n nanosize producing a surface of reactive
thiol and/or functionality. Following a novel two-step silica coating pro-
cess, reactive amino groups were successfully on particles.
Examples have been given for magnetic carriers of tailored functionalities in a
variety of applications, including industrial effluent detoxification, recovery,
cell separation, drug delivery, and the preparation of biological sensors.
I. Introduction 52
I I . Acid-Base Properties of Solid Surfaces 53
A. General definitions 53
B. Extension to solid surfaces 56
C. The acid-base behavior of species confined in porous solids 58
111. Origin of Acid-Base Properties of Solid Surfaces 59
A. Ideally polarizable surfaces: metals 60
B. Open surfaces: oxides 61
C. Surfaces with variable composition: composite oxides 71
D. Layered materials: clays 75
E. Porous materials: zeolites and zeotypes 76
F. Microporous materials: activated carbons 78
G. Pseudo-liquid materials: heteropolyacids and salts 80
IV. Measuring Acid-Base Properties of Solid Surfaces 81
A. Nature, number, and strength of acid-base surface sites 81
B. Experimental methods 81
C . Theoretical concepts 88
V. Manifestation of Acid-Base Properties of Microporous Solids 91
A. Solid/gas interface: the role of acid-base sites in heterogeneous
catalysis 91
B. Solid/liquid interface: acid-base reactions in interlayers of
layered minerals 93
VI. Summary and Conclusions 94
References 95
Bronsted-Lowry
When formulated the classical views on acids and bases in 1887, his
theory dealt exclusively with reactions in aqueous media He defined as acids
all substances that upon reaction with water increase the concentration of
ion in solution; those that increase the hydroxyl ion ( O H ) concen-
tration were defined as bases.
In 1923 Bronsted [2] and [3] independently generalized these definitions
and emphasized the crucial role of protons. Acids were defined as substances able
to donate protons and bases as substances able to accept protons. Water was
longer required as a solvent, and bases no longer needed to donate
The following equilibrium defines the pairs of conjugate acids and bases and is
fundamental to the Bronsted-Lowry theory:
Acid A and its conjugate base B are two different but interrelated states in which a
chemical system may exist. The acid is the proton-rich (and proton donor) state,
while the base is the proton-deficient (and proton acceptor) state.
The reaction in the forward direction leads to the formation of a base in system 1
the conjugate base of and an acid in system 2 the conjugate acid of
An important consequence is that the position of equilibrium depends on the
relative acidity basicity of conjugate acid-base pairs in the two systems,
and 2.
2. Lewis
In a further generalization, Lewis advanced a definition that was no longer
restricted to proton-exchange reactions. An acid was defined as any substance
that can accept a pair of electrons from a donor substance, and a base as any
substance able to donate an electron pair and form a "dative" covalent bond.
these more general the acid-base interaction includes the formation of cova-
lent bonds and applies to any chemical reaction in which an addition compound
is formed through a coordinative bond:
According to the above definition, both hydrogen bonding and bonding are
base interactions.
Between Lewis acid-base are possible displacement reactions that for-
mally resemble the reaction from the Bronsted-Lowry theory:
Thus the acids of the Bronsted-Lowry theory may be viewed as Lewis acid-base
in which the acidic component is (for example, can be viewed as
an The Bronsted-Lowry bases are also Lewis bases. One can
therefore conclude that the Lewis definition subsumes the definition
of acids and bases.
The heats of formation with many Lewis acids have been measured
experimentally in the gas phase using the technique of ion cyclotron resonance
The absolute proton affinity (PA) of a gaseous base molecule has a precise
thermodynamic definition in relation to the negative of the enthalpy variation for a
hypothetical reaction of attachment of an isolated proton to a molecule M in the
gas phase,
in the absence of any interactions with the The corresponding free energy
variation is referred to as the gas-phase basicity. Many compilations of various
experimental results have been made I], and several useful correlations for ratio-
nalizing the acid-base chemistry have been revealed. Quantum-mechanical calcula-
tions for the enthalpy of in the gas phase generally give
overestimated values
to Solid
Bronsted and Lewis acids and bases are frequently encountered in the chemistry o f
solid surfaces, but in some instances these terms are used in a slightly different
context from their original proposals. From the point o f view o f the Bronsted
definition, a solid is acidic i f it is able to donate (or at least transfer partially) a
proton to a basic probe molecule to which the proton bonds (or with which it, at
least, becomes partially associated). According to the Lewis definition, a solid is
acidic i f it is able to accept an electron pair a basic molecule and also a
coordinative bond with it. The transfer o f an electron pair in the meaning o f the
Lewis definition should be clearly distinguished from surface processes in
which electron transfer ( o f one or two electrons) occurs without coordination.
Extension o f the above definitions to acid-base processes on solid surfaces is
further complicated by the "in between two phases" location o f the acid-base sites.
For example, little is known about acid-base interactions at inter-
faces. It has long been recognized that molecular interactions across the interface
between condensed phases may be split into physical and chemical terms. Physical
interactions (van der forces)contribute to the non-ideality o f fluids and have
been traditionally considered for interfaces. However, modern theories explain
interfacial phenomena on solid surfaces, such as adhesion or wetting, in terms o f
chemical interactions. Moreover, the Lewis definition o f acidity is so comprehen-
sive that it can easily be accepted that most chemical interactions at solid surfaces
may be effectivelydescribed as acid-base interactions A continuously growing
literature reinterprets interactions at interfaces in o f acid-base
properties. For example, their role was shown in relation to adhesion
in film-substrate or fiber-matrix systems as well as in wood and paper processing,
The void space within porous materials constitutes a region with special properties,
many of them not completely understood. At the level of molecular interactions in
nanoscale-sized spaces it is difficult to distinguish between "pure" acid-base pro-
cesses and perturbations due to the confined environment. Both the solid matrix
and the molecules confined to the internal open spaces may show altered acid-base
behavior in comprison with their "free-space" counterparts.
T o illustrate this we select two examples that correspond to two classes of
interfacial phenomena, both important from the practical point of view.
A interface is present in the microporous environment inside the chan-
nels and cages characteristic of zeolite structures. All catalytic chemistry based on
the acid properties of crystalline (zeolites) is affected by two addi-
tional factors known as the "solvent" and sieve" effects. The chemical
contribution, the solvent effect, describes the interaction at the molecular level
between micropore-confined molecules and the coordination field created by all
atoms that constitute the pore walls. The molecular sieve or geometical effect
accounts for the perturbation of the dynamics of confined molecules and is parti-
cularly important in relation to the selectivity of catalytic reactions occurring over
this class of microporous solid acids. Both these effects may act concurrently to
distort the intrinsic acid-base properties of molecules confined in the micropore
spaces. In the case of alkane cracking, numerous observations demonstrate that
various zeolites show very different specific rates for the same catalytic reaction.
Two opposing views were expressed to account for this: (1) The -(OH)- Si
positions in zeolite lattices span a very broad range of intrinsic acid strengths
that determine various reactivities, or (2) all acidic sites are virtually identical, and
differences in reactivity result from geometric constraints in the micropores
The perturbation due to geometrical factors may be so important that the existing
experimental techniques cannot separate the effects. so any connection to the
base index of the solid zeolite surface becomes obscure It is known, for
example, that the heat of adsorption of gaseous hydrocarbons on zeolites depends
on the chain length or the diameter and shape of the micropores Catalysts with
comparable numbers and strengths of acid sites but with different pore geometry
The development of a scientific basis to describe the chemistry of solid surfaces has
been impeded by their ill-defined nature; many factors can be responsible for
observed effects. An understanding of surface chemistry at the same level
Partially dehydroxylated
surface
Lewis Lewis
immersion
Fully
surface
IV. V. VI.
Transformation of surface structures and their acid-base characteristics
on oxide surfaces with various degrees of dehydroxylation.
in terms of acid-base sites associated with individual surface ions located at various
positions related to all possible terminations of the bulk crystal lattice.
1. Surface Morphology
The morphology of oxide crystallites depends largely on their preparation method
as well as on the prevailing crystalline phases. In general, the higher the tempera-
ture of calcination, the coarser and more crystalline is the final powdered oxide. For
most transition metal oxides the surface area cannot exceed certain limits that
reflect the microcrystalline organization of the solid For these materials of
low porosity the average pore size is much larger than molecular diameters. Their
surface is "open" for molecules impinging from the gas or liquid phase, and on the
molecular scale the interaction zone can be viewed as a "flat" surface.
Materials with "open" surfaces are ideal for characterization of intrinsic
base properties of individual surface sites. The interaction of impinging molecules is
essential bimolecular, and generally the surface does not show any peculiar
The molecules impinging on the surface can have either a Lewis base or Lewis
acid character. For instance, CO with a filled orbital at the carbon end displays
Lewis basicity, while shows Lewis acidity The primary interaction is
strictly localized on oxides, and surface are formed with the acidic or
basic centers of the solid. An example is the nondissociative adsorption of
on exposed centers on oxide surfaces:
(VIII)
Another example of dissociative adsorption by concerted action of a Lewis
base couple is bonding of alkenes on metal oxides, as in reaction (IX). However,
dissociative adsorption of water on ideally dehydroxylated oxide surfaces is by far
the most important example to cite, in view of the formation of surface
groups and their tremendous impact on the reactivity of most oxides under
or close to ambient conditions [reaction (X)].
three) were regularly found on transitional aluminas; the larger number is the result
of a distribution of ions between octahedral and tetrahedral positions, which
discriminates between O ‘ vibrations. Also, ceria is likely to present only densely
Ã
packed planes on its surface because of its type of lattice, and thus only two
(instead of three) types of hydroxyl groups are revealed in I R spectra
A dynamic, rather than static, picture of partially hydrated oxides is more
appropriate, as confirmed by proton NMR and dielectric spectroscopy studies
The results of revealed two types of protons o n partially
droxylated oxide surfaces: those associated with physisorbed water and those from
chemisorbed water. About 85% of the total protons are present in the first two
molecular layers. Subsequent NMR studies showed that a permanent
proton exchange occurs between chemisorbed and physisorbed water.
In conclusion, partially dehydroxylated oxide surfaces exhibit a large inventory
of surface O H groups and water molecules together with Lewis acidic and Lewis
basic sites with coordinative (structures and of Scheme 1). The
hydroxyl population is the source of protons that cause enhanced surface electrical
conductivity and catalytic activity. It is significant that the increase in the conduc-
tivity value is paralleled by increases in either the amount of weakly bound protons
or their mobility Almost all metal oxides are active in catalytic
of which is one of the least demanding reactions in terms of the require-
ments for the acid strength of active sites Studies on several oxide systems
show that the activity is lost after extensive dehydration and is partially restored by
5.4 (CC)
7.7 (CC)
9.5 (CC)
7.2 (S)
6.7-7.8 (DL)
3.1-3.3 (DL)
6.0-6.3 (DL)
1.2 (DL)
Data predicted theoretically by the MUSIC model and measured experimentally by P A D and other
techniques.
DL = double-layer model; CC = constant capacitance model; spectroscopic methods.
Pauling's bond valence rule has been used in various attempts to correlate acidic
properties of surface hydroxyls with their formal charge, as expressed by the above
relationship For example, the assignment of stretching vibrations of
surface hydroxyls already mentioned is based on calculations of the formal
charge on surface hydroxyls using Pauling's rule. Based on the possible relationship
between acidity and the stretching frequencies of hydroxyls in similar structures
the above reference to Pauling's rule illustrates, in fact, a relationship between
structure and properties (acidity) for surface hydroxyls on partially dehydrated
oxides. Pauling's valence bond rule served as a basis for the MUSIC model of
proton binding to surface hydroxyls on fully hydrated oxides. The results of
these semiempirical approaches are in remarkable agreement with exact quan-
tum-mechanical calculations based on local Hartree-Fock models and with
many experimental trends.
A comparison of Figs. 1 and 2 shows that FTIR spectra of isolated hydroxyls on
oxides have their counterpart in of the fully hydroxylated surfaces. The
parallelism is seen with respect to both the number and strength of hydroxyl groups
present on the surface. A plot of stretching frequencies groups (from Fig.
1) versus the values revealed by of several oxides (Fig. 2) is shown in Fig.
3. Despite some scatter of the data, a linear trend is evident. At this stage one can
only speculate that the slope of the straight describes the intensity of the
electrostatic at the oxide surface. It apparently depends on the sample's
chemical nature (possibly, through the cation's electronegativity or
radius ratio) and/or its surface structure (anion packing, preferential development
of crystal planes, surface imperfections). In any case, the correlation shown in Fig.
3 is similar to a relationship found for zeolites in which the increase in Bronsted
acidity correlates with decreasing bond strength and can be monitored
through the corresponding decrease in the stretching frequency in FTIR spectra
ions, the degree of exposure of surface cations, and their preference for defect sites.
For example, in the oxygen lattice is more loosely packed, the frequency of
stacking faults is greater, and the fraction of cations in tetrahedral positions is
slightly smaller than in although the local ionic configuration of
cations is the same in both modifications. Due to these structural differences, the
acidity, surface reactivity, and catalytic activity of are significantly higher
than those of
Surface imperfections of crystallographic structure may occur in small
oxide crystallites. Metal or oxygen vacancies are typical imperfections that cause
a local charge imbalance. Therefore, the acid properties may be strongly affected by
the crystalline size, in parallel with (and perhaps as a consquence of) variation of
electronic properties. For anatase particles with sizes in the range of 75-100 A , the
formation of new acid sites was observed with higher strength on smaller particles
1. Mixed Oxides
The acidity of materials can be controlled by mixing oxides with different proper-
ties and in various ratios. If phase separation does not occur, substitution of cations
of the minor component oxide in the host oxide lattice results in bulk phases of true
mixed oxides. Another method used to the acid-base properties of oxide
increase of acidity
increase of basicity
Activated carbons consist of small graphene-like layers that are spatially distribu-
ted without regular organization and are separated by most of which
are slit-shaped. The most reactive sites for the attack of molecular oxygen are the
edges of the graphene layers and defect sites. Upon oxygen attack, and depending
on reaction conditions, two types of carbon surface oxides are Because of
their practical importance, surface oxides of carbons have been the subject of many
studies
Basic carbon oxides are formed upon the exposure to oxygen of a clean carbon
surface obtained by outgassing to followed by cooling to room tem-
perature. Aqueous suspensions of these carbons have a strongly basic pH, reducing
properties toward noble metals, and anion-exchange properties. Acidic carbon
oxides are when the carbon is oxidized a t lower temperatures
the presence of water vapor during the oxidation process increases the
amount of acids formed. When dispersed in water, these carbons have acidic pH
and cation-exchange properties. Acidic surface oxides can also be formed by oxida-
tion at room temperature with aqueous solutions of nitric acid,
hypochlorite, or hydrogen peroxide.
(XVII) (XVIII)
Solid-state heteropoly compounds have interesting catalytic properties and
have found applications in and acetylene hydration and alcohol dehydration
Adsorption of nonpolar molecules (hydrocarbons) reaches saturation at
amounts lower than monolayers ("surface-type" adsorption). However, it is char-
acteristic of heteropoly that adsorption of polar molecules (alcohols,
ethers, amines) takes place in amounts that exceed by many times
where [B] is the concentration of the unreacted base. The acid strength of surface
sites is estimated by comparing the color of several indicators with different
values adsorbed on the solid sample, while the number of sites in various
ranges is estimated from the of base (in most cases,
needed to block the sites. The suffers several limitations. The
most serious is the fact that solid samples are in nonaqueous media
benzene) but the values of indicators used for comparison are aqueous phase
values, and the medium-transfer corrections for are ignored. Also, equilibrium
may be achieved with difficulty and the method does not distinguish between
Bronsted and Lewis sites unless hindered bases, which are selective for
Bronsted sites, are used as
An alternative method proposed by Contescu et al. is aqueous titration
combined with a suitable algorithm for deconvolution of the measured proton
binding curve into a continuous spectrum of apparent values. Even neglecting
corrections for charge accumulation at the interface. the content of the
proton affinity distribution (PAD) function is complete with respect to both the
number and the apparent strength of proton-binding sites at the inter-
face. Agreement with the number of surface acidic sites was demonstrated by
comparing the PAD data for activated carbons with the results obtained indepen-
dently using an established analytical method (Boehm titration The validity
of the evaluation of values was demonstrated by comparing the experimental
results for oxides with "open" structures such as alumina with the acid strength of
surface predicted by the MUSIC model
Another method for quantitative determination of acid sites on solid catalysts
is calorimetric nitration with amines This method uses the very simple
concept that the enthalpy of adsorption for simple bases is a direct measure of
acid strength. If a variation of the heat of adsorption with surface coverage is
where is the specific rate constant for catalysis by an acid whose acid dissocia-
tion constant is G and are constants for a series of similar catalysts.
It was assumed by several authors that a similar correlation should exist for
heterogeneous catalysis but the difficulty of measuring the number and
strength of Bronsted acid sites on solid catalysts at the same time hindered any
experimental proof of a complementary relationship. Recently, the occurrence of a
Bronsted-type relationship was demonstrated experimentally by Contescu
for isomerization of 1-butene on a series of mixed oxide catalysts. The number
and strength of Bronsted acid sites were simultaneously determined from the pro-
ton affinity distribution of solid acid catalysts by potentiometric titration For
a homologous series of catalysts, such as or
the areal reaction rates on catalysts of various molar compositions were directly
correlated with the number of acid sites of a particular type (Fig. identified as
the bridging hydroxyl in Si- OH- M structures. The significant fact was that
active sites identified in each binary oxide system were characterized by nearly
constant acid strength values (expressed on the scale of water), which shows
that the acid strength in Si- OH- M structures is indeed determined by the
electronegativity of M. Moreover, when the logarithms of specific reaction rate
constants were plotted versus the values for the active sites found on catalysts
from various binary systems, the linear correlations shown in Figs. and were
found This is the equivalent in heterogeneous catalysis of the Bronsted rela-
tionship, first developed in homogeneous catalysis.
The fact that a relationship to Eq. (8) was experimentally found for solid
acid catalysts has both practical and theoretical significance. a practical point
of view, it demonstrates that at least one experimental exists, namely the
PAD method, that affords a direct and quantitative determination of both the
number and strength of Bronsted acid sites on surfaces of porous oxide catalysts.
Concepts
I11 recent literature a considerable lack of agreement can be found as to which
acidity scale is more suitable as a reference for molecules probing the acidity of a
solid surface (gas-phase or aqueous-phase scale). This dichotomy presents a pro-
blem that is difficult to resolve and has interpretive implications for experimental
results.
The difficulty of establishing a unique acidity scale for solid surfaces has been
mentioned several times. Even the definition of "acid strength" of solid surface sites
is not thermodynamically consistent. It is appropriate to analyze some of the
theoretical concepts involved and to address some of the solutions advanced by
various authors who tried to solve this problem.
The definition of acidity in terms of the Lewis theory has a broad meaning that
encompasses the Bronsted definition as a particular case. In practice, especially in
the field of heterogeneous catalysis, it is important to differentiate between
Bronsted or Lewis acid sites on solid catalysts, since they may show completely
different catalytic behavior.
Limiting the discussion to acidity, there have been two opposing views
expressed regarding the question of which is the most appropriate acidity scale for
solid surfaces. Some authors advocate the use of the acidity scale of gaseous species;
others support the use of the scale of aqueous solutions.
in the absence of any interactions with the medium. The latter requirement is
difficult to satisfy for accurate measurements even in vacuum. However, if the
absolute proton affinity for a single base, could be independently
in the gas phase, then a relative scale of basicities could be established by
measuring the equilibrium constants of a series of gas-phase reactions with
T o bring to closure the plethora of possible examples on this subject we cite two
more. Each illustrates a specific domain of intrinsic practical importance in which
the acid-base properties of microporous solids play a dominant role in determining
a broad range of applications. The examples examine the and
interfaces, respectively. For the interface, the interdependence of
base properties and microporous structure is revealed by the catalytic properties of
zeolites, especially by their shape-selective properties. It is impossible to review all
the literature in this area, and only a few recent views are summarized. For the
interface, we focus on the role of acid-base chemistry in interlayers of
clay minerals.
Zeolites are widely used as molecular sieve adsorbents and solid acid catalysts in the
oil refining industry. The most important characteristic of these microporous solids
is their acidity, in particular its interplay with the chemical environ-
ment.
There are many differences between traditional acid-base chemistry in solutions
and acid functions on zeolites, but the general aspects of the chemistry between
The scope of this review was intentionally restricted to those and liquid/
solid interfaces where the acid-base properties at a molecular level can best be
explained and related to particularly important applications.
In the field of heterogeneous catalysis, in which most solid acids or bases used as
catalysts also have porous structure, the questions of primary interest are (1) how
to characterize in a consistent manner the number and strength of surface acid-base
couples, ( 2 ) how to distinguish between Lewis and Bronsted acid sites, and (3) how
to prepare solids with high surface area, selected micropore size, and stable surface
structures with desired acid properties. The progress made in this field is commen-
surate with its economic importance. Now it is possible to select and prepare
catalysts with an acid strength adequate for almost any specific acid-catalyzed
reaction. This has been made possible by the development of acid-base character-
ization methods that provided the basis for understanding the close relationship
between the atomic-scale structure and the properties of surface acid sites. The
progress in materials science and the development of molecularly engineering cat-
alysts have extended the spectrum of inorganic materials used as solid acid cata-
lysts. The new generation of zeolites and zeotypes with controlled molecular sieve
properties and an acid strength as high as 10 has induced a revolution in the
field of oil refining and petrochemistry. When catalysts with even stronger acidities
are required, the use of sulfated oxides is promising, although more work is needed
to stabilize these against premature deactivation.
Synthetic layered hydroxides with hydrotalcite-like structures are used as het-
erogeneous catalysts for base-catalyzed reactions. One of their major drawbacks for
technological application is their low stability during hydrothermal and/or thermal
treatments. Recent attempts to incorporate large, preferably multicharged anions
into their interlayer spaces introduced new routes for preparation of more stable
materials. The effect of intercalated polymetalates or transition metal complex
anions in the framework of layered hydroxides is twofold: They increase the gallery
heights and thermal stability on one hand and introduce additional Lewis or
Bronsted acidic sites in a basic framework on the other hand.
Oxides, especially those of Si, Fe, and Mg, are essential components of the
earth's minerals; their surface sites may have either Bronsted or Lewis acid-base
properties. The weathering of rocks and the formation of soils are processes at the
interface between minerals and natural waters or solutes from the
medium in which solid phases are either formed, altered, or dissolved.
Dissolution of solid minerals is controlled by slow surface reactions rather than
I. Introduction
II. Electro-Optical Effects in Suspensions of Spherical Particles
A. Electro-optical methods
B. Characterization of the samples
C. Electro-optical responses of spherical particles
D. Discussion
I I I . Electro-Optical Effects in Suspensions of Anisometric Particles
A. Experimental results
B. Discussion
References
Colloidal electro-optical studies started with the discovery of the Kerr effect [I].
The theories of colloidal electro-optical phenomena were strongly influenced by the
classical molecular electro-optics theory [2], particularly by the part concerning
dipole molecules The latter explains the display of negative Kerr constants
by the existence of permanent dipole moments perpendicular to the particle sym-
metry axis. The first formal theory pretending to explain the extremely high
values of the Kerr constant for colloids attributed them to particle volume.
Experimental studies demonstrated that the processes causing electro-optical phe-
nomena in colloids are much more diverse than those in molecular systems
Deforming the ionic atmosphere surrounding the colloidal particle, the electric
field, causes, on the one hand, particle electrokinetic motion (hence flow-induced
particle orientation), while on the other hand it creates an induced dipole moment,
polarization-induced particle orientation, the Kerr effect.
modes Our research is directed to the possibility that they could grow to be
large-scale and very significant in the hydrodynamic frequency range and could
influence the optical response of the system both directly and through acoustic
orientation of the particles. Thus for interacting systems self-orientation and col-
lective orientation modes are coupled with self-translation and collective transla-
tion modes, resulting in a frequency behavior of the electro-optical
responses. Although the precise mechanism of low-field electro-optical sign reversal
is not yet well understood and our studies are in the initial stage, we shall try to
demonstrate here the close similarity of the electro-optical behavior in the hertz
range of spherical and anisometric particles, hence the importance of
tic modes for the electro-optical anomalies in charged colloids.
An essential feature of this study is the choice of model samples. They are rigid
particles in highly stabilized suspensions (liquid and crystal phases), hence inter-
ference of field-induced deformation and aggregation of particles is avoided. The
highly monodisperse isotropic spheres are amorphous silica, thus neither form nor
intrinsic optical anisotropy can be related to the particles, and surface effects
attributed to latexes are also avoided. The anisometric particles are thin rods of
the clay mineral paltgorskite of axis ratio 20. The sample is and its
basic advantage is the great amount of electro-optical data accumulated in the
course of its application as a model for research into the origin of the electro-optical
phenomena. It displays a very large Kerr effect at high dilution and all features of
low-frequency anomalous behavior in the semidilute regime.
The electro-optical responses of colloidal systems were detected by reflection
spectroscopy and by the variations of scattered light intensity. Induced variations
in the reflection spectra of colloidal crystals turned out to be very sensitive to
electroacoustic effects, while the experimental setup of the second method permits
direct comparison with numerous data on anisometric colloids. Investigations into
reflection spectroscopy were carried out by Okubo and coworkers a t Gifu
University (Japan), and the light scattering studies by Stoimenova, Okubo, and
coworkers at the Bulgarian Academy of Sciences. The essential parts of the results
are due to our joint research.
General Aspects
Conventionally the electro-optical effects of interacting colloidal systems are
presented by the statistical average of the product of the optical effect of a single
particle and the orientation distribution function
a =< >
where the optical part is a function of the components of the particle-scattering
matrix. Thus, for field-induced birefringence and dichroism, the optical factor is
determined by the diagonal elements of at forward scattering:
An' -
while the induced variations of scattered light intensity are followed by any pre-
ferred combination of the components of light polarization) as a function of
the observation angle:
The above presentation illustrates the similar bases of the electro-optical meth-
ods y l ] . The components of the scattering matrix can differ in their sensitivity to
the various sources of particle anisotropy, but they yield similar information on the
induced anisotropy fluctuations, particularly on their relaxation. This scheme is
true in so far as noninteracting particles are concerned. Though the correlations
in particle orientation have been accounted for in recent electro-optical the-
ories, the influence of the applied field on the space distribution of the scattering
is generally neglected.
In interacting systems the optical and orientation factors in a are no longer
separable quantities. The induced optical effect is determined both by the single
particle scattering [form factor and by the pair distribution function [structure
factor the latter being direction-dependent Since in addition to orienta-
tion the electric field causes particle translation, even for spherical particles the
radial distribution function and the "static" structure factor attain
dependent induced anisotropy. The deformed surface potential also contributes to
this effect.
Formally introducing modified values for and (in the general case the
process involves quantities), we can write for the steady-state elec-
trical light scattering,
of the
Highly monodisperse spherical particles were the main object of our studies. The
colloidal silica spheres CS-81 and CS-91 used in the investigation are products of
Catalyst Chemicals Ind. Co. (Tokyo). By electron microscopy (TEM,
Type the spherical mean diameter and standard deviation were
determined to be 103 13.2 for CS-81 and 110 4.5 for The charge
density of strongly acidic groups on the spheres determined by
titration with a autobalance precision bridge were 0.38 and 0.48
c m , respectively. Before use the samples were purified with an ultrafiltration cell
and treated continuously on a bed of and anion-exchange resins
20-50 mesh]. The resins were used for the deionization of
water needed for the suspension preparation. The suspensions were obtained by
dilution of the stock solution with water and were further kept in contact with the
ion-exchange resins for a week in order to achieve the crystal state (conductivity
below ohm-' at Crystals in the dilute suspensions were very large
which is an advantage for light scattering observations of the dynamic
modes.
Turbidity of liquid like and crystal-like suspensions is very low compared to that
of systems due to the diminished intensity of multiple scattering for visible
light. Particularly low turbidity is observed for highly dilute and exhaustively deio-
nized samples This made it possible to follow the static and dynamic light
scattering behavior of the silica suspensions in gaseous, liquid, and crystalline states
Data analysis was performed with the and histogram
including Marquadt analysis. A single exponential decay was found only in the
colloidal gases, corresponding to free particle translation diffusion. The diffusion
coefficient is always lower than the value calculated with the Stokes-Einstein equa-
tion with the true diameter of spheres and correlates well with the effective diameter
including the atmosphere. The liquid samples display bimodal decay, corre-
sponding diffusion coefficients (differing by a factor of 3) being related to
and long-term self-diffusion of particles. static light-scattering variations for
the crystal state display very sharp peaks in their scattering vector dependence. The
analysis of crystal structure dynamics yields three apparent diffusion coefficients
that differ from each other by an order of magnitude: the basic one related to
particle (restricted) Brownian diffusion, a slow mode assigned to propagating
(0.1 %), below 1 Hz; for semidilute samples they cover the hertz interval, increasing
to higher frequencies with increasing volume fraction. In parallel a dc component is
observed, which, depending on volume fraction and ion content, can cover the
entire frequency range (up to 100 kHz), typically related to orientation effects. In
the semidilute regime the low-frequency sign reversal is observed, in close similarity
to the low-field responses of anisotropic systems. At higher sphere concentrations
(above 5% volume fraction), the responses are again positive.
The field intensity dependence of the ac and dc components of the responses in
the intermediate concentration range are complicated: linear at very low fields and
quadratic or rapidly saturating at higher fields On a larger scale they show
oscillatory behavior. In Fig. 4 they are shown for the transition concentration of
sphere volume fraction for moderate field intensities. The component
increases with field intensity and can be observed up to higher frequencies.
Modulated signals are obtained up to in the range of the negative dc
component.
Figure 5 shows typical responses for this particular sample (0.5% sphere volume
fraction). A basic difference in the electro-optical behavior of spherical compared
to anisometric particles is the linear response a t low field intensity. The modulated
signal follows the frequency of the applied field with a phase difference of
With an increase in frequency the ac component relaxes, changing in both ampli-
tude and phase and passing to double frequency, while the dc component changes
-2
0
0 1 2 3
log (Hz)
Frequency dependence of the ELS responses from deionized suspensions of
silica spheres (sphere volume fraction 0.5%) for different field intensities. 150
from a negative value to a small positive value. The decay process displays two
basic relaxation times: one, associated with the positive kilohertz effect, on the
order of milliseconds and the second, related to the low-frequency modulations,
on the order of seconds. Additional slower modes can be observed upon applica-
tion of longer pulses, as can be seen in the kilohertz response. This figure concerns
the liquid state of the sample.
The electro-optical responses are strongly dependent on the phase state of the
colloidal system. The transition of the same sample the liquid to the crystal-
line state is accompanied by a considerable increase in the cutoff frequency of the ac
component and a displacement of the dc curve to more negative values (Fig. 6).
Similar variations can be observed by increasing of the sphere volume fraction, as
demonstrated in Fig. 3, which shows that the effective fraction determines
the process. In parallel, additional modes are observed to appear in the crystalline
state. They introduce resonance peaks at frequencies depending on the sphere
volume fraction (two fractions are presented in Fig. 6).
Typical responses of the crystal were shown earlier in Fig. 2. Below 1 Hz
they are similar to the responses of the liquid phase (Figs. 2a and 2b) and are
evidently related to induced longitudinal density waves. At low fields they follow
the field frequency (Figs. 2d and and at higher fields a harmonic frequency
appears (Fig. In the interval the responses show oscillatory decay.
Estimating the parameters involved, we identified this mode as related to shear
waves propagating in the crystal lattice Hence the registered peak frequencies
are a relative measure of the elastic modulus of the system. The kilohertz effect is of
reversed sign and displays distinctly two relaxation rates (Fig. The first one
0 1 2
log f (Hz)
Frequency dependence of the ELS responses deionized suspensions of
CS81 (field intensity for different sphere volume fractions at liquid and
crystal states of the sample.
FIG. 9 Time-resolved reflection spectra for one cycle of the applied field (0.2 Hz,
for a deionized (crystal-state) suspension of CS91 of fraction 4.4%. Upper
figure: Corresponding integrated intensities. (Adapted from Ref. 48.)
AP) V
Several figures are presented in this section to illustrate the similarity of the
frequency behavior of spherical and anisometric particles. They were obtained by a
conventional electro-optical method (electric light scattering) on a typical and
studied colloidal sample The anisometric particles are rigid rods of the
clay mineral palygorskite (attapulgate). Their average length and diameter are
200 and respectively. Particle surface charge in the aqueous suspensions
is about A more detailed investigation of the sample is presented in
Refs. 58-60.
Figure 10 is a typical picture of field intensity and frequency variations of the
electro-optical responses of highly dilute aqueous colloids. The frequency curves of
the anisotropic particles (below the overlap concentration) are significantly differ-
ent from those of the spheres, due to the large effect in the kilohertz range. The
decay of the responses of the investigated sample is nearly exponential, the time
constant of several milliseconds (2-7 corresponding to particle rotational diffu-
sion. A tail on the order of 50 ms can still be detected. The high-frequency relaxa-
tion is proved to be related mainly to surface polarization of the particles, which
causes their orientation. The low-frequency relaxation interval coincides with the
cutoff frequency of the component (hydrodynamic range), where large variations
log f (Hz)
. 12 Frequency dependence of the steady component of the ELS response from
suspension (volume fraction 0.01 for different intensities in the low-field
range.
effects shown by the more concentrated silica samples. Figure 14 shows several
typical responses of the anisometric particles.
Generally speaking, three frequency ranges are clearly distinguished in the electro-
optical dispersion curves of anisometric particles:
Below 1 H z the range of large-scale density fluctuations, which are not
averaged in the optical response and can be detected both directly and by
the induced flow orientation
2. The region around 100 Hz o n the order of particle reorientation frequency,
where translation motions reduce particle orientation
3. The Kerr line, which dominates at higher fields but only until the particle
rotation is not restricted by interaction forces
Hence, two basic factors are to be taken into account in the analysis of electro-
optical data. First is the pulse duration. A t certain values of the pulse length the
response shows steady-like behavior, which is due to the separation of the time
scales involved and the restricted sensitivity of the optical response to induced
changes in the system. Second is the choice of frequency range for the determina-
tion of the Kerr constant. The latter is related to the low-field regime, where the
acoustic modes contribute significantly. Adequate separation of orientation and
20 Hz, 60
deionized
translation effects requires analysis of the responses over a wide range of field
intensities and frequencies.
The results presented for spheres reject a general view lying a t the basis of all
schemes proposed for explanation for the low-frequency sign reversal in colloidal
electro-optics the necessary presence of anisotropic o r deformable species in the
disperse system as a condition for the detection of electro-optical responses. This is
not a necessary condition for the observation of low-frequency negative
optical responses either. The essential role of particle electrophoretic mobility (giv-
ing rise to electroacoustic effects) for the latter is evident. The results presented
concern a series of phenomena in different systems widely studied by researchers
involved in the field of colloids, etc. Below we
turn our attention to some of these phenomena.
It is to be noted, in the first place, that the acoustic modes in semidilute polymers
and colloids are quite similar hence the optoelectroacoustic phenomena can
I. Introduction
I I . Special Considerations
A. Nuclear properties
B. N M R signal narrowing effect in colloids
C. Relaxation rates in colloidal species
D. Free solution species in the colloidal system
E. Adsorption phenomena
III. N M R Studies on Colloidal Oxides and Hydroxide Systems
A. Colloidal silicas
B. Colloidal aluminas
C. Other colloidal oxides
IV. Conclusions
References
7.04
1.12
5.54
Oxygen-17 shows all the negative aspects of both and allied with a very
low natural abundance, which usually requires isotopic enrichment. Hence,
NMR has not been used very often in the characterization of colloidal oxides,
although there are notable exceptions described in the literature
Unfortunately, the properties of the magnetically active nuclei and are
very unfavorable for N M R measurements. These nuclei show a relatively low
natural abundance, low sensitivity, and a high quadrupole moment. fact, tita-
nium has received little attention from the N M R spectroscopist. Therefore, despite
the great practical importance of and associates in the colloid state, it is hardly
possible to find any reference on the use of or N M R for the characteriza-
tion of these oxides.
Effect in Colloids
Colloidal sols are arbitrarily regarded as particles falling anywhere in the range of
in diameter. However, a smaller upper limit will probably be required if
the colloidal oxide species to be observed and interpreted in a spectrum by
means of conventional Fourier transform N M R in the liquid state. First, the par-
ticles must form a stable colloidal suspension ensuring the homogeneity of the
medium. This is a basic condition to avoid a gradient of concentration in the
volume being analyzed and guarantee equal excitation of the nuclei within the
spectral range. Furthermore, the size of the particles should be sufficiently small
to allow their incoherent isotropic tumbling motions so that both dipolar couplings
and chemical shift anisotropy can be eliminated. Therefore, Brownian motion is
expected to sharpen the signals of the colloidal particle and also of related
adsorbed species. In fact, the use of colloidal suspensions for N M R analysis in
specific cases has been claimed as an alternative tool for N M R analysis of solids
replacing the magic angle sample spinning method where rapid mechanical
rotation of the solid sample is applied.
v (MHz)
spin 1/2 nuclei of such particles will often be in the threshold region of the extreme
narrowing o r far from that. Furthermore, within the usual range of NMR fre-
quency (60-750 MHz),the magnitude of may be of crucial importance for the
relaxation behavior, especially in the range of between and
2. Quadrupole Relaxation
F o r the case of quadrupolar nuclei in molecules in the extreme narrowing region,
the quadrupole coupling constant influences NMR spectra through very efficient
It is worth noting that the linewidth of this resonance decreases with increasing
This behavior is in sharp contrast to the situation described by Eq. (2) and to the
dipolar relaxation process, where the linewidth increases with increasing From
the field dependence of the chemical shift and described in (3) and
the values of 7 and for spin 5/2 nuclei in the slow motion limit can be
obtained. As the parameter provides a measure of the coupling between the
nuclear quadrupole and the electric field gradient at the nucleus, it reflects the
of the local environment of the nucleus and can give information
on the nature of the colloidal species, is a measure of the mobility of the species
involved and, in principle, can provide information on the magnitude of the particle
size when (1) is used.
Species in the
The colloidal state inevitably brings about difficulties for the experimentalist when
separation of the disperse phase from the dispersion medium is needed. This is the
case when the speciation and concentration of only the free soluble species have to
be determined. Separation of the ionic solution from the small colloidal particles
for conventional analysis is nontrivial, although separation techniques
such as ultracentrifugation. dialysis, and field-flow fractionation have been success-
fully used. If the soluble species of interest have an active nuclear spin, the liquid
NMR technique will constitute an alternative and simpler way to characterize and
quantify those species without being affected by the disperse phase. An exception is
the case where the colloidal species gives a signal that fully overlaps the sharp
resonance of the solution entity. As is quantitative, the absolute concentra-
tion of the species can be estimated based an internal reference of known
concentration but different chemical shift relative to the sample signals.
Alternatively, a calibration curve can be established from a set of external standard
solutions (preferably the same substance found in the sample) measured under the
same experimental NMR conditions as those applied to the sample.
Care should be taken when the soluble species show a significant interaction
with the surface of the disperse phase. In this case, a modification in the NMR
spectrum profile may reflect an adsorption phenomenon, which should be properly
interpreted.
where and are the fractions of free and bound water and i 1, 2 stands for
longitudinal and transverse relaxations, respectively. The average relaxation rate of
water molecules in the corresponding colloidal dispersions will give information
about the particle structure and affinity of the surface for the solvent molecule. This
is illustrated in Fig. 4, where the specific water proton spin-spin relaxation rate is
plotted as a function of surface area (different particle size) of silica and alumina
sols [The specific relaxation rate - where is the relaxa-
tion rate of pure water and is related to the overall colloidal dispersion given by
The fact that the data fall on the same line is a clear demonstration that
N M R can be used to monitor the total available surface area in a dispersion in
which is not straightforwardly achieved by other means. The straight line depen-
dencies indicate that the samples satisfy the fast exchange limit. Besides, the higher
slope of the alumina line indicates a stronger solvent-surface interaction than that
for silica. More detailed studies of the water layer can be made by using and
relaxation measurements, which are sensitive to changes in motion
Moreover, splitting of and N M R signals in surface-bound
water molecules can elucidate structural modifications on oriented particle systems
as in sol-gel transition N M R relaxation studies of water protons using para-
magnetic probes as relaxation rate enhancers also provide a convenient way to
investigate ion adsorption onto the surface of colloidal oxides
When surface species are not in fast exchange, direct observation of the spectra
can be used to probe both the surface orientation and the dynamics of the adsorbed
species.
OR OR
-Si - Si -
I
OR
Q2 Middle Units Branching Units
-Si -
- Fully Condensed
Tetrafunctional Units
Stable particulate alumina sols are precursors for many materials such as controlled
porosity catalysts, ceramics, and coatings. are commonly formed by
of by dilute monoprotic acids. They can be
prepared either by the process developed by Yoldas [4] from hydrolyzed aluminum
alkoxides or by acidic peptization of the precipitates derived from hydrolyzed
aluminum salts As will be seen these alumina sols are usually
component fluids that also contain soluble aluminum species. colloi-
dal can also be by the of hydrolyzed
cations. Therefore, a brief review the characterization of aluminum species in
solution by N M R is presented in studies on
are then discussed
1. Studies on Hydrolyzed Aluminum Species
Aluminum exists aqueous solutions as the monomeric ion
below pH 3 or as ion above At intermediate pH values, for
an molar ratio 0 and 2.5, is hydrolyzed removal of
the coordinated water followed by condensation of the O H groups
to polynuclear ions. this process leads to the growth of large
colloidal species and the formation of sol-gels, N M R spectroscopic studies
have contributed remarkably to establish speciation in hydrolyzed aluminum solu-
tions. These studies were pioneered by et and notable
have been made by Bottero et al. Bertsch et al. and more recently
Kloprogge coworkers Many species with different degrees of polymeriza-
tion and different internal hydrolysis ratios have been reported in this complex
However, only the monomer and the
have been unequivocally assigned. These species
show sharp resonances in the spectrum at 0 ppm (octahedral region) and
temperature sols a t 363 K does not transform them into the form
directly but results the of tetrahedral and octahedral species (Fig.
probably resulting from condensation, A feature of all spectra in
Fig. 9 is a small broad peak a t 7 p p m which been attributed by the to
overlapping peaks: one a t 4.2 tentatively assigned to the "dimeric" species
and the other a t assigned to colloidal species.
130.3 M H z after 60 min (Fig. 1 la) showed a sharp peak a t assigned to the
monomer and its hydrolysates (6) and In addition, a rela-
tively broad peak (linewidth at slightly overlapped with the
monomer resonance. This latter, referred to as decreased in intensity with
the process, concomitant with the growth of another extremely broad
octahedral resonance, referred to as SB2, centered a t approximately 8 p p m (Figs.
with other rigid structures such as those found in silica systems. It is also
worth noting the presence of monomer even in the semisolid gel. probably moving
freely in the interstitial water enclosed in the continuous solid structure.
It is obvious a t this point that NMR can be a powerful tool for examining
the chemistry of alumina sols, but the most interesting result obtained by
Morgado al. was the unequivocal assignment of resonance to boehmite
colloidal particles. This was possible only because of the field dependence shown by
this species. The spectrum of the sol in Figs. 1 l a and 12a are seen to be quite similar
to that in Fig. 8a (except for the very signal at and it would be
tempting to assign species as the solution species discussed in
the previous section. However, a remarkable shift to lower frequency and broad-
ening of the resonance in the boehmite sol were observed at lower
magnetic fields, as depicted in Fig. 13. Morgado could demonstrate that this
field dependence was due to a second-order quadrupole-induced shift, which is
observed in solids but only for species of limited mobility in the liquid
state, far from the extreme narrowing condition (cf. The chemical shifts
and linewidths determined a t four magnetic fields were plotted against the inverse
of the square of the nuclear resonance frequency as shown in Fig. 14 (data denoted
by solid squares). A quadrupolar coupling constant of = 22. MHz and an iso-
tropic chemical shift of were calculated respectively the slope of the
line and from the intercept in the plot of 5 in accordance with (3)
(assuming = 0). A comparison of the value with those of other aluminum
described in the literature indicated that it was intermediate between
those of sites of very high symmetry (such as the tetrahedral site of with
7 0.65 MHz) and those of sites of low symmetry (such as the octahedral site of
with = but very close to 7 values for sites of reasonably high
symmetry such as in Similarly, a value of = 3.5 x was calcu-
lated from the slope of the line in the plot of according to by
using the previously estimated value of 7. This value of correlation time was much
higher than those observed for the soluble complexes = 4.2
800
600
400
0 1 2 3 4 200
1 2 3
nuclear field (x nuclear field (x
FIG. 14 Dependence of (a) chemical shift and (b) linewidth of the colloid
octahedral resonance as a function of the nuclear magnetic field (see also Table 4).
and line) = 0.10; = 0.08; -
= 0.06. (From Ref,
Properties of the
XRD data
Code amor p hous c ) (nm)
SBS 36
SBF 27
18
Pural 200 P-200 0
B 36
B2 39
B
supplied by Condea GmbH (synthesis from hydrolysis of
gels synthesized according to Ref. 5 by neutralization of
relative to the highly crystalline boehmite P-200, assumed as 100% crystalline (0%
amorphous).
= crystallite size.
Nuclear magnetic resonance studies have proved valuable for the characterization
of other colloidal oxides and hydroxides. Some selected examples of such studies
are shown in the following paragraphs.
Pozarnsky and used N M R to propose for
the synthesis of aqueous sols. In contrast with the colloids discussed in the
previous section, the N M R signal due to colloidal species (vanadate polymers)
has mostly broadened beyond the detection limit, probably due to its very low
symmetry, which enhances the problem of quadrupole interaction. Nevertheless,
based on the solution chemistry associated with the sol, the process could
be rationalized in terms of the chain polymerization of the dioxovanadium cation,
from This polymerization could be rate-limited by the
release of the cations from decavanadate anions (-426, -512, -531 if
these anions were initially present at high concentration. Using an external refer-
ence spectrum, Pozarnsky and found different starting concentrations
of the decavanadate anion that were dependent on the synthesis method. In the
slow synthesis method, the decavanadate ion accounted for 80% of the
vanadium present, whereas in the fast synthesis method it accounted for only 20%.
Therefore, a reaction scheme that explains the N M R data (after a brief initial
transient period) could be obtained, thus rationalizing the polymerization mechan-
ism.
Judeinstein and used and N M R to analyze molecular precur-
sors formed during the sol-gel synthesis of thin films but were not very
successful in using N M R for their unequivocal characterization. Three
broad resonance peaks appeared, but a definite conclusion was not possible
because the chemical shifts might be related either to bonding distortions in a single
oligomeric molecule or to different molecular species. In both cases the broadening
of the N M R signals could be explained by short relaxation times arising from
formation of oligomeric species and exchange reactions between different chemical
species.
King et applied and N M R to provide structural information
on HTMA (hexamethylenetetramine) used as a agent in the sol-gel
I. Introduction
I I . Experimental
A. Materials
B. Methods
I I I . Results and Discussion
A. Surface density and wettability on alkyl-inimobilized glass
substrates
B. Effect of various alkyi chain lengths on contact angle and
surface energies
IV. Conclusions
References
conditions of the bonded phases. The motion of such chains is the key
to understanding their chromatographic properties Many researchers, includ-
ing and Squires have investigated dual models of monomeric alkyl
surfaces ranging from a folded or "liquid-like" orientation to a rigid or bristly
type in an aqueous environment.
Electron spectroscopy for chemical analysis (ESCA) [9] and Fourier transform
infrared (FTIR) spectroscopy with attenuated total reflectance (ATR) can also
be used for routine surface studies. FTIR spectroscopy is known to have the
sensitivity to determine the average orientation and reorientation of interfacial
chains. But it does not directly provide information on the itself. The
mobility of a solute in the neighborhood of an alkyl chain can be measured by
fluorescence spectroscopy
Direct measurement of molecular motion and dynamics is offered by NMR
studies using the spin-lattice and spin-spin relaxation times In a
general sense, the changes in spin-lattice relaxation time were consistent with
increased segmental motion as a function of distance from the surface. Sindorf and
investigated the molecular motion of n-alkylsilane bonded to silica
particles using N M R . Alkyl chain mobility increase with chain length
up to = 8) until the interactions between chains occur at longer chain length
= 18, for example). However, the main disadvantage of methods is
that they are limited to fine particles having high surface areas and are not applic-
able to plate-type surfaces.
The plate contact angle method might provide a sequential scanning
curve or hysteresis loop that can be interpreted in terms of surface mobility, reor-
ientation, solvent penetration, and intrinsic under both water and air.
Water is known to be a difficult liquid for use in contact angle studies due to its
small molecular resulting in penetration and local swelling of the solid
surfaces
According to Van et the receding contact angles on
surfaces decrease in the range of = 6 to = 12, then increase to
= 18. Because of the increasing mobility due to the decreasing polymer sur-
faces with moderate-length alkyl side chains may reorganize, which may lead to a
decreased number of hydrophobic segments exposed to water. In addition, as a
minor factor, the hydrophilic ether groups may also reorient toward water surfaces.
This result may be closely related to another study of the brittle temperature
for acrylates and methacrylates. They showed a minimum in the brittle
temperatures in the ranges of = 8-12 due to increased side chain mobility.
Judging by these cases, alkyl side chain bonded surfaces may have different motions
with different chain mobilities and lengths.
It is not easy to measure such chain length effects directly on polymer
substrates because of the packing density problem polymer swelling effects
etc. Therefore, immobilized on glass substrates can be used as a
good model system for the study of surface dynamics. Because rigid glass
surfaces are nondeformable, unlike polymer substrates, a rigid glass substrate may
have a contact angle very different from that of a polymer substrate. Pure clean
glass substrates contain groups high surface energy), and the contact
angles are zero in a pure aqueous environment. Hence, we may be sure that the
II.
. Materials
Trimethyl chlorosilane, chlorosilane, n-octyldimethyl chlorosilane,
and chlorosilane were purchased from Petrarch Systems, Inc.
and used without further purification (all 99.5%). The chemical structures of
these silanes are shown in Table 1. present as minor constituents in
reagents for the n = 1 case], contain no reactive groups.
Hence, they were removed during the rinsing procedure.
Prepurified water was further purified by passing through a reagent
water system Co.). Toluene was dried with molecular sieve particles
for 24 h before use. Methanol and dehydrated ethyl alcohol were used as received.
Corning borosilicate cover glass (2940, No. 1 24 x 50 mm, 0.16-0.19 mm
thickness) was used as substrate. As determined by ESCA, surface elemental com-
positions of this borosilicate glass were B 2.3, 28, 0 = 63, C = 4.4 atomic
percentages, respectively. Minor components below 1 were Na, and
Ti.
Cleaning Glasses
Glass slides were immersed in chromic-sulfuric acid solution at for 40
which supplied active oxygen atoms that oxidize carbonaceous materials on glass
surfaces, rinsed extensively with purified water, and then dried overnight in a clean
oven at in air. All other glassware used was also carefully cleaned by the
same process.
We examined the cleanliness of the glass slides by measuring the surface tension
of purified water(72.6 0.5 at room temperature and by verifying
that there was no water contact angle hysteresis and that the slides were perfectly
wetted.
Trimethyl chlorosilane
Me
Me
-
Me
Me
chlorosilane Si
Me
system designation
RM-2) as a function of the immersion depth. A motor drives the shelf to advance or
withdraw the water container at a controlled speed. The chamber is insulated and
maintained at room temperature and constant humidity (35% RH).
The output from the electrobalance measurement is fed to an X-Y recorder, with
the balance output (force) feeding the axis and immersion depth feeding the
axis. To measure the contact angle, the silane-treated glass slides were immersed
into or withdrawn from the purified water or mixed solvent at a constant speed of
Dipping velocity was chosen to be fast enough for convenience of
measurement but slow enough to avoid speed effects (such as hysteresis on the
clean glasses) as shown in Fig. 2. Advancing contact angles were used for this
study because they can show the surface energy condition directly, whereas reced-
ing contact angles can show the surface energy condition of the wetted surfaces.
Water-methanol solvents were made with increasing fraction. Each
solution was kept for 15 min at room temperature to equilibrate; then its surface
tension was measured. The contact angle and surface tension of each sample were
measured three times, and their error ranges were within and
respectively.
The total free energy is the sum of contributions from the different
forces at the surface,
where the superscripts d and p refer to the dispersion and polar force components,
respectively.
The surface free energy of a interface can be described as a function
of individual force components of the solid and liquid. F o r the relationship, the
following equation, which was assumed by Owens a n d Wendt is most reliable
for polymeric solids:
Combining (1) and (3) yields the well-known geometric mean equation,
By measuring the contact angle 9 of two different liquids against a solid, simul-
taneous equations are obtained that can be solved for and if the values of
and are known. With some organic liquids, the components of free energy
were evaluated using reference solids. The result shows that the components
deviated by over 10% from the mean, and they cannot be accepted as constant.
The samples of different alkyl chain lengths (n 4, 8, and 18) and different alkyl
surface concentrations obtained by varying the silane concentration
indicates no tested
error range within 0.5.
4 -3 -2 -1 0
LOG CONC (MO
-5 4 -3 -2 -1 0 -5 4 -3 -2 -1 0
LOG CONC (MOLE'S00 LOG CONC TOLUENE)
contact angles and high polar surface energy. Also the disperse surface energy term
showed a lower value. However, the total surface energy exhibited higher values.
At high concentration the bonded alkyl chains are more prevalent and are
probably relatively densely packed at the higher silane solution concentrations.
Therefore, modified surfaces with alkyi chains can show a larger contact angle.
Therefore a lower polar term of surface energy, like paraffin surfaces. However, a
higher disperse term of surface energy was shown as expected, causing the total
surface energy to decrease.
Water can usually be used as a unique solvent for the study of biomedical
materials. In other applications such as reverse-phase liquid chromatography, how-
ever, different solvents such as water-methanol mixed solvents are used. As the
fraction increases, the surface tension of the solution decreases. In mixed
solvents, the trend of contact angles as a function of concentration was generally
similar to that of the aqueous case. As the proportion of methanol increases, the
contact angle decreases dramatically as in the case of decreased concentration.
Several factors govern the exact nature of alkyi chain surfaces. First, the function-
ality of silane materials may affect their monomeric or nature.
Second, the degree of and the surfce concentration may influence the
accessibility of water molecules to the unreacted hydrophilic groups remain-
ing on the underlying glass surfaces. Third, the length of the alkyl chain is a main
factor in this work. It is expected that for moderate chain lengths, the chains are
relatively mobile There are interactions between the chains van der
forces) at long chain lengths These different chain types, folded or extended,
may affect contact due to the different outermost groups exposed
and According to I R spectroscopic and ellipsometric studies by
Porter et the long (about = 18) form a densely packed,
"crystal-like" assembly with fully extended alkyl chains. As alkyl chain length
decreases, the surface structure becomes increasingly disordered and liquidlike with
lower packing density and coverage. Fourth, the effect of an aqueous environment
is different from that of the aqueous methanol solvent environment. Although we
usually use water for biomedical studies, for chromatographic studies the mixed
solvent is very important. In a mixed solvent environment, the bonded alkyi chain
can be 'bristly" or
Figure 4 show total surface energies as a function of concentration with four
different chain lengths. Except for the n = 18 case, there are transitions in surface
energies just below This results in two contact angles for = 4
and 8 in two different relatively high (0.1 and and low and mol)
concentration ranges, and the possible reason is described below.
Figure 5 shows the advancing contact angles as a function of alkyl chain
length at four different high and Sow surface concentrations mol alkyi
toluene). With increasing carbon number the advancing contact
angles decrease until = 8 at relatively high concentrations (0.5 and 0.1
or decrease until 4 at relatively low concentrations and
FIG. 4 Total surface energies as a function of concentration with four different chain
lengths.
.5 Advancing contact angles for four different concentrations (0.5, 0.1, and
as a function of carbon number in the alkyl silanes.
Contact angles obtained with the Wilhelmy plate method and near surface energies,
d
and obtained with nonlinear programming methods were varied with chain
length surface concentrations for monochlorosilanes 1, 4,
8, and 18) in various solvents.
As solution concentration decreases, total surface energy and the polar term
increase, whereas the disperse term decreases, due to surface density with
18
I. Introduction
I I . Hierarchical Structure of Collagen Fiber
A. Hierarchical structure of collagen fiber
B. Surface structure of collagen fiber
III. Pretreatment and Inner Structure of Collagen Fiber
A. Pore size distribution determined from adsorption isotherm
of water vapor on collagen fiber
B. Pore size distribution determined from adsorption isotherm
of nitrogen gas on collagen fiber
IV. Experimental Procedures
A. SEM, TEM. and AFM
B. Adsorption method
V. Summary
References
Collagen, the essential component of the extracellular matrix, is present in the skin,
bone, tendons, and connective tissues of the higher animals. The collagen fibril is
functionally related to the tenacity in bone, the tensile strength in tendons, and the
preservation of strength in skin or connective tissues. In the native state, collagen
exists in a form. Water is essentially related to the maintenance of collagen
molecules in the native conformation. That is, the inner spaces inside the collagen
fibers functionally maintain water in the native state.
Observations of the collagen fibril surface structure were performed using
atomic force microscopy (AFM) scanning electron microscopy (SEM)
period
crosslink
collagen nonhelical telopeptide
nonhelical telopeptide
ridge
et water molecules are fixed in the microfibrils and in the hole zones
adjacent to the end of the tropocollagen. It was reported that intermolecular
space enclosed by the five tropocollagen in a unit is
This indicates that there is microspacing inside the
3. Collagen Fibril and Collagen Fiber
Microfibrils are packed on a tetragonal lattice o r quasi-hexagonal lattice to
form the collagen fibril. The collagen fibril is the fibrous entity observed by A F M
or The diameters of the collagen fibrils are close to a multiple of
The collagen fibrils form a fiber, and the fibers are bundled together
with elastic fibers and fibrocytes to form fibrous connective tissue.
FIG. 5 FEG-SEM picture of critical point-dried fibrils ( x 160,000). (From Ref. 2.)
The D period of the rat tail tendon collagen was reduced by dehydration from 68
to 64nm The micropore volume of the collagen fiber was increased by
acidification Therefore, the inner structure of the collagen fiber was also
Control collagen
Enzyme treated 11
collagen
Tanned collagen
. 13 Pore size distribution of collagen fibers. The labels C-0, CS-0, C-10, CS-10. C-
20, CS-20, and CS-30 refer to the number of days the raw collagen fibers were
alkali-treated. C denotes the alkali-treated and collagens. denot
alkali-treated, freeze-dried, and acidified collagens. 1.2
1.55 2.45 (From Ref.
The pore size distribution of collagen fiber is determined by using adsorption iso-
therms of water vapor o r nitrogen gas. The preparation of the sample. adsorption
procedure, and calculation method are mentioned below.
1. Preparation of Sample
Collagen fibers are from a n organic matrix. They are thoroughly washed
in distilled water and prepared through a process of several steps. The sample
should be dried. The fiber is dried in a vacuum desiccator over phosphorus
oxide for 3 months to completely remove water from the collagen fiber The
water inside the triple helix [5] be removed. If the water inside the collagen
fibers is not completely removed, then the amounts of nitrogen gas adsorbed will
vary widely. This wide variation in the amount adsorbed may be caused by the
phase change of the residual water inside the collagen fibers due to an extremely low
adsorption temperature of
2. Adsorption Procedure
The water vapor and nitrogen gas adsorption isotherms determined by
using automatic adsorption apparatus have been reported. O n other hand. the
manual adsorption method for the adsorption of water vapor and nitrogen gas is
discussed in Refs. 25 and 4, respectively. The adsorbed is determined by
using a Shibata moisture and surface area apparatus, Model P-850
Scientific Technology, Ltd., Tokyo, Japan) The of water vapor
and nitrogen gas adsorbed are determined by the gravimetric method and
volumetric method respectively.
The pore size distribution is calculated based on the water vapor or nitrogen
gas adsorption isotherm. Before the calculation, it is essential to obtain a n accu-
rate adsorption isotherm. The amount of nitrogen gas adsorbed on the collagen
fiber is extremely low in comparison to the amount of water vapor adsorbed on
it. Approximately. 1 x kg of collagen fiber. therefore. is put in the adsorp-
tion glass bulb. The amount adsorbed on the collagen fiber is measured at equili-
brium pressures of 6.7, 13.3, 20.0, 26.7, 33.3, 40.0, 53.3, 66.7, 80.0, 93.3, and
101.3 in order to get an accurate quantity. T h e adsorption isotherm for each
Relative pressure
. 14 Adsorption and desorption isotherms of nitrogen gas on collagen fibers. C-0,
C-20, and C-30 are collagen fibers alkali-treated for 0. 10. 20, and 30 days,
respectively. Experimental data of adsorption and desorption. (From Ref. 4.)
downward from straight lines. The downward deviations are used for the calcula-
tion of the pore and the pore size distribution of the
u p to 1.55 The pore size distributions for radii of 1.55 n m or larger are
calculated using the D H The cumulative pore size distributions cal-
culated using the MP and D H methods are shown in Fig. 12.
The surface structure of collagen fibril is visually described using the SEM,
and TEM observations. The of collagen fiber is quantitatively
described using the pore size distribution calculated from adsorption isotherms.
The depends on the achievements of other researchers. The latter is mainly
my work. It is regrettable that the SEM o r AFM observation cannot be compared
with the pore size distribution o n the same collagen sample. The functions of
micropores inside the collagen in the native state must be elucidated.
I. Introduction
II. A Conceptual Model for Sorption of Inorganic Substances at
the Oxide/Water Interface
III. The Diffusion Model
IV. Implications from the Conceptual and mathematical Models
V. Critical Review of Relevant Literature
VI. Conclusions
References
solution
porosity in reactor
apparent of adsorbent
However, for oxides,
microposous structure
intraparticle diffusion For oxides, Axe
Anderson showed transfer resistance is to
particle diffusion
a spherical particle, the species mass balance yields the
In this expression, defines the mean square hopping distance, and the probability
that the adsorbate has the energy required to h o p to the next site is the
factor with representing the activation energy. The adsorbate vibrational fre-
quency. is a function of the molecular mass and the force constant, which is
dependent o n the potential energy on the pore surface.
In this chapter, we d o not distinguish between surface and pore diffusion but
instead talk about an overall effective diffusivity. In this way, we can include all the
different approaches in our review. In summary, the mass balance model is based
on the following assumptions: The adsorbent is composed of uniform, porous
spheres with uniformly distributed sorption capacity, sorption is reversible, and
diffusion is not a function of the particle dimension. These conceptual and math-
ematical models provide information about how sorption kinetics and the general
system behavior can be influenced by changes in environmental conditions.
If inorganic sorption onto porous oxides follows the conceptual and mathematical
models presented above, certain types of behavior should be observed. F o r exam-
ple, in many aquatic environments, amorphous Fe, and oxides often pre-
cipitate in the presence of other dissolved ions. Under these conditions, adsorbates
that have a n affinity for oxide binding sites should have access to those sites during
the formation process before they become sequestered in internal pores. Therefore,
when coprecipitation occurs, intraparticle diffusion is probably not a rate-limiting
step in the sorption process. Coprecipitation conditions, however, usually are not
equilibrium conditions, because the amorphous oxides continue to age and trans-
form into more stable mineral forms. During this aging process, adsorbates that are
associated with the changing solid could become involved in several different reac-
tions. If it does not fit well in the evolving mineral lattice, the adsorbate could be
released to solution. Under these conditions, initial removal during coprecipitation
would exceed the equilibrium conditions, and aging would cause a n increase in
adsorbate concentration in solution. In contrast, if it readily substitutes for
major metal ions in the lattice, the adsorbate could become fixed in the oxide. I n
all cases, the rates of aging and of the oxide, and possible release of
adsorbate, would be controlled by intraparticle diffusion.
Adsorbate size relative to the size of pores in the oxides is also an important
factor to consider. the IUPAC classification, it seems reasonable that inor-
ganic adsorbates should have access to binding sites in pores of the size
Cd H F O (pH 7 )
Cd H F O (pH 6)
Cd
Alumina
Sr HFO
'Effective diffusivities summarized here represent an average from each study; actual values were
typically within an order of magnitude of these averages.
= hydrous ferric oxide. also known as amorphous iron oxide or ferrihydrite.
values, which in some cases incorporate the effects of porosity, tortuosity, and
adsorption. Also, some diffusivity values are estimates o r averages from the infor-
mation presented in the original paper. It is clear that pore size plays a significant
role in diffusion. In general, values for the adsorbents with what the
I U P A C characterizes as micropores and (Table 1) are around
In contrast, values for adsorbents with (Table 2) are
typically u p to six orders of magnitude larger. T o better understand the effect of
diffusion. ranging from to were used in the integrated
analytical solution to the mass balance using a constant boundary condition at the
particle surface and a particle diameter of 40 (Fig. In systems with
pores, intraparticle diffusion approaches the rates expected for bulk diffusion. As
pore size increases from micropores to mesopores to macropores, diffusion rates
The adsorbate's charge and solution composition will determine its size, which
affects the resulting type of intraparticle diffusion. Specifically, when pores are
large relative to the diameter of the ion, bulk diffusion is possible. F o r pores
that are small relative to the size of the ion, transport is limited to solid-state
diffusion. When pore and ion are about the same size, the ion interacts with the
pore wall as a result of physical forces. such as electrostatic and van der
forces Diffusion under these conditions is known as configurational diffusion
In aqueous environments, layers of water adsorb in what is typically defined
as mesopores, and interaction of the adsorbate with the pore wall results in con-
figurational or surface diffusion. The size differences among anions compared to
those among cations would suggest that anionic diffusion may fall into the
state regime for certain cases, while hydrated cation can move by both surface
and pore diffusion.
Of the studies reviewed here, only the work by al. directly exam-
ined the effects of pore size. Under similar experimental conditions. they reported
that when adsorbents with micropores and were compared. diffusivity
values for Cd and increased by factors of about 2 and 33, respectively. As
noted by those authors, the fact that is more sensitive to changes in the pore
size is likely due to the fact that hydrated selenite anion is large relative to C d
(diameters of about 6.5 A and 4.5 A, respectively). The larger selenite anion be
significantly hindered while diffusing through smaller diameter pores.
et noted that in contrast to postprecipitation adsorption
tests, uptake of onto amorphous iron oxide was significantly greater and
not diffusion-limited. Furthermore, subsequent aging of the system, they
observed desorption of presumably because it did not fit well in the structure
of the aging oxide. In contrast, there apparently are some where the
adsorbate does fit well into the structure of the aging oxide. F o r example,
et reported increased uptake of on gels aged for 205
days. The behavior could not be due to slow diffusion into the solid, however,
because in the same systems the amount of Sr, and Ba decreased with aging.
Trade-offs between diffusion and incorporation into the solid can have interesting
implications for the transformation pathways of these unstable solids.
reported a correlation between the relative affinities of anions for
and enhanced formation of hematite rather then goethite.
Because amorphous oxides are metastable minerals, over time they will crystal-
lize into more stable forms. F o r example, amorphous oxide, also known as
ferrihydrite, will restructure to hematite. However, upon dissolution of the amor-
phous Fe oxide, it would reprecipitate to goethite. These structural changes in the
oxide, which are a function of time and solution chemistry, will affect the
particle diffusion process. However, Axe studied amorphous oxide aging
with experiments that included spectroscopic analyses of the before, during,
and a t the conclusion of sorption experiments. X-ray diffraction was used to
observe the oxide aging, with and without contaminant present. The diffraction
pattern of freshly precipitated amorphous Fe oxide (Fig. 2) showed little o r no
cyrstallization; however, as it aged, more structure could be seen. Over the
patterns were analyzed for the peak widths a t half-height. After 7 weeks, there is
evidence of both hematite and goethite minerals, but the oxide without Sr had more
goethite present than the one with This crystallization increased with time; after
14 weeks diffraction pattern peaks narrowed and increased in height (less noise). At
14 weeks, there were n o significant differences between the Fe oxide with and
without The fresh, 7-week, and 14-week patterns were from one set of samples,
but in the final time series, 22 weeks. a different sample was used. F o r the 22-week
samples, the diffraction pattern of the Fe oxide without has stronger, narrower
peaks than does the pattern of oxide with An analysis of the peak widths
indicates that for the oxide without Sr there is evidence of more crystal-
lization. As a result, the 5.5-month-old samples suggest that the presence of con-
taminant inhibited the crystallization of the amorphous oxide. Interestingly. the
week-old samples showed less crystallization than the samples a t 14 weeks. Overall,
I. Introduction
I I . Electrokinetic Equations
A . Flux of the ionic species
B. Flow velocity field
C. Electric potential
D . Boundary conditions
E. Equilibrium ion density distribution
F. Equilibrium ion density distribution in spatially periodic porous
media
III . Linearized Equations
A. Perturbations
B. Simplified equations
C. Dimensionless formulation
IV. Examples
A. Porous medium submitted to external pressure and potential
gradients
B. Semi-infinite void space
C. Plane Poiseuille configuration
D. Circular Poiseuille configuration
V. Numerical Simulation of Electro-Osmotic Phenomena
VI. Porous Media
A. Regular packings of spheres and ellipsoids
B. Random packings of spheres and ellipsoids
C. Reconstructed porous media
D . Discussion
VII. Fractal Aggregates
where is the fluid permittivity, the viscosity, the surface potential, and the
electric field parallel to the wall. Overbeek has shown that Eq. (2) applies locally
everywhere in the pore space provided that E represents the local electric field. If
the double layer volume is negligible, this implies the following relation between the
transport coefficients in (I):
double layer have been systematically studied and discussed by comparison with
the previously discussed ordered packings.
Reconstructed media based o n Fontainebleau sandstones were generated and
studied for three double layer thicknesses. These important data are discussed and
compared with existing results
Section VI ends with a general relation in the limit of large double layer thick-
ness that, to the best of our knowledge, is novel. This yields a single
reduced representation of the electro-osmotic coefficients for all the studied con-
figurations as functions of a dimensionless double layer thickness.
Section VII with suspensions made of fractal aggregates, because of the
practical importance of such structures that may be present in circumstances.
The conductivity of such suspensions is determined for charged and uncharged
clusters; the electrophoretic mobility is calculated in the dilute limit.
The chapter ends with some concluding remarks.
All the relevant microscopic equations that govern electro-osmotic phenomena are
given in this section, as well as a general description of spatially periodic porous
media.
We will be concerned with dense porous materials wherein the percolation flow
velocity is normally small and the characteristic Reynolds number is much less than
with
Potential
The electric potential is given by the Poisson equation
where is the dielectric constant of the electrolyte. It is assumed that for any
nonsteady process of ion transport, transient phenomena associated with the elec-
tric potential occur so fast that one can use (4) in a quasi-steady approximation,
assuming that = , t ) satisfies (10) with time being a parameter. This
that the electromagnetic wave propagation characteristic time is
smaller than the characteristic time associated with the diffusive species
transport , and being the characteristic interstitial scale and velocity of light,
where v is the outer normal to the particle bed surface its surface charge
density, and zeta potential. The two boundary conditions, namely the Dirichlet
or Neuman condition on the electric field, are found in the literature.
with the obvious zero-flux solution, which, with the help of (4) and be
written as
where is the zeta potential of the particle surfaces. These two conditions are
related by the condition that the sum of all the charges within the unit cell is
equal to zero:
indefinitely reproduced within an infinite space (see Fig. 1). Therefore, the solution
of (6)-(11) is to be sought in the three-dimensional space with
,
where {I,. are the three basic vectors characterizing the unit cell of the porous
medium and the trio of integers belongs to is the
position vector within the unit cell.
At equilibrium. all properties appearing in (14) and (15) are spatially per-
iodic functions. depend only on the local position vector r :
i= ..., = = (18)
and may be considered only within the unit cell o n the external faces of which the
periodicity conditions should be imposed. These will be formulated in the form
=0 =0
where . . . denotes the difference between the values of function at the opposite
points r lying at the corresponding unit cell boundaries.
u= = 0
In the above, the j subscript on jump conditions signifies that the component of the
jump condition along the unit vedtor I, is written for any value of
Introduce the ionic potentials by means of the relationship (cf.
Because of the combination of Eq. (24) with Eqs. (4) and (13) for the flux
of ionic species yields
Since the ionic flux across a particle surface is nil as well as the fluid velocity, ionic
potentials are subject to the boundary conditions
and using the solution for the flow velocity field at equilibrium, one can
show that
We will further assume that all the perturbing quantities are small, which is valid
only when is small with respect to the electric field prevailing in the double layer.
A condition for V P will be established below. In addition, only sta-
tionary situations are considered. Accordingly, the time derivatives and the
order perturbation terms are dropped in (24)-(31).
Then the conservation equation for the ionic species, is reduced to a
linear differential equation.
. = 0, = 0, u = 0 on
If there is no macroscopic gradient of ionic concentration +
F r o m the equations governing the flow velocity field (34). one can devise a char-
acteristic
These scales are used to define the reduced quantities, which are denoted by a
prime:
This system (44)-(47) is to be solved subject to the boundary conditions inside the
unit cell,
In addition, under the assumption that either the surface zeta potential o r the
surface charge is small, the approximation can be used to
linearize (45) and with respect to , This was done in the numerical
The electric current density flowing through the porous medium, resulting from
the application of is
Because of the linearity of the previous system of equations and boundary condi-
tions, all solutions should be proportional to E. The conductivity tensor of the
is defined by
A last physical situation may be of general interest. Previously, the jump con-
dition for the ionic potentials was calculated from (24) assuming that the
ionic concentrations have no macroscopic gradient =
= =
If this condition is not fulfilled, for example with
+ =
then the ionic concentrations are subject to a jump condition.
Integrating twice while accounting for the boundary conditions a t the walls
= and far from it (uniform we get finally
=
1 + , tanh
1-
Note that for small potentials can be linearized into
The tangential velocity increases from 0 at the wall to the value far from it,
which is given by formula [4]
or dimensionally
If the height of a plane channel (see Fig. 2b) is assumed to be very large
compared to the thickness then ( 5 8 ) is a solution for the
undisturbed potential distribution Under this condition, the conductivity
and electrokinetic tensors can be calculated analytically.
On the other hand. if the channel is not very thick compared to the layer,
Poisson's equation must be linearized to be solved analytically when the surface
potential is assumed to be low <<
Two species are present in the fluid with = = 1 a t the same initial con-
centration = Therefore, (45) becomes
The average velocity induced by an axial external electric field can be calculated
analytically. as well as the electric flux induced by a pressure gradient,
In the above, the subscript denots the components of the vectors along the
direction of the channel. we deduce the coefficients a' and
which are equal, in agreement with theorem:
These coupling coefficients are very small = 3' in a narrow channel and
tend towards in a wide one 1).
The electric flux under a given electric field gives the conductivity coefficient,
Assume that the porous medium is composed of a series of parallel tubes (bundle of
tubes model) each of diameter (Fig. 2c) The unit cell is now formed by the
region 0 a, 0 where L is a n arbitrary length. The external field and
the pressure gradient V P are applied parallel to the tube axis. T h e capillary axis is
and the distance to this axis is
Owing to the cylindrical symmetry, (44) for reads, for two
species with 2, 1, and iff
The solution of the linearized equation involves the modified function of the
first kind
Note that
Therefore, for wide capillaries the potential along the tube axis is
Again. is verified.
The expression for is truncated to order since our expression for is
only first-order accurate as in the previous case of a plane channel. F o r wide
capillaries >> 1).
Note that these approximations are generally second-order. They are third-order
if there are two symmetric species with = This linearization is generally
called the approximation
The geometry of the porous medium is discretized into a n array of
elementary cubes of size a filled with fluid o r solid. These cubes are by
mesh points. The equations are solved in a finite difference formulation with a
second-order discretization, which yields linear problems of the generic form
The vector is the unknown field or the matrix results from the operator
in the transport equation and from the local boundary conditions and the
right-hand term B results from the driving force, from the global boundary
conditions (49) and the current value of u' when solving for or the current
value of when solving for
The fields are are determined by a program that solves in any geometry
any linear differential equation of the general form
not
converged
calculation of U
Printing
End
Ellipsoids:
= 0.208. = 2.5 x
The above permeabilities roughly constitute the averages of the values obtained
along the various directions considered in the orthorhombic packings, while the
conductivities are closer to the smallest conduction.
or equivalently, because =
0-1, ellipsoids
ellipsoids
spheres
RP, ellipsoids
RM.
2A3
Plane channel
Circular channel
=simple cubic; 0-1, = types I a n d random
packings; reconstructed media.
F o r thin double layers = the data are still more scattered. Three values
are far below the others and correspond to the least accurate computations. When
these are excluded, one obtains
It is interesting to note that the variations in the relative conductivity for the
reconstructed media are comparable to those in the granular media. F o r all the
arrangements of spherical o r ellipsoidal grains considered above. including random
packings, we obtained for = 2.46
A new important result of these computations is the relation between the coupling
coefficients and the permeability for thick double layers [cf. Eq. The case of
relatively thick double layers (fine porous media) was largely overlooked in the
literature. Attention was generally focused on the opposite limit of vanishing
length, where the matching technique of inner and outer developments was exten-
sively used (see I), which yielded the general relation (3) of Overbeek
The analytical treatment for thick double layers is a formidable task because of
the layers' overlap, which requires full resolution of the governing equations.
However, for plane and circular channels, one can use formulas (66) and (78) to
obtain
Moreover, for thick double layers. with larger than the characteristic pore size,
is nearly uniform and approximately equal to the surface potential
Therefore, replacing by and using the definition of permeability, one gets
Note. however, that the convergence in (98) is not necessarily uniform along
different directions, as illustrated below.
The limit ( 9 5 ) has already been pointed out in several cases in the previous
section. Since in most practical situations the length is unlikely to exceed
widely the typical pore size, it is essential to study the rate of convergence of the
coupling coefficients toward the limit (95) to determine how far off this limit is from
the actual value in a n intermediate range, and possibly to introduce a correction
factor. T o this end, all the results of the previous section for all the investigated
directions in all the geometric models were gathered.
In order to plot ratios against a normalized double layer thickness, a
length scale applicable to all the directions and geometries has to be defined. As
already pointed out above and also shown by the criterion for conver-
gence toward ( 9 8 ) corresponds to the covering of the pore space by the electrical
double layer rather than its thickness relative to a typical grain size. Thus, an
adequate length scale should be a measure of the characteristic pore size. Such a
length scale A was introduced in Ref 42. It is essentially a pore volume to pore
surface ratio, with a measure weighted by the local value of the electric field in
the conduction process:
This length scale A is derived from the correction to due to the interfacial
conductivity in thin layers.
An electrokinetic method was proposed to measure A , as a n alternative way
to determine the permeability via
The constant in (103) results from a least squares fit of our data (Fig. 10) Note that
(102) can be recast into
where is the translational velocity of the aggregate and its gyration radius.
This formula is in agreement with the classical Stokes formula. since for a solid
sphere = 3 and =
FIG. 12 The coefficient for uncharged aggregates as a function of the solid fraction
Each point represents an average over 10 realizations. Data are for
Sander (modified: see text), D 2.5; hierarchical with linear trajectories. D = 1.9;
Thouy and D 1.6; (+) Thouy and Jullien, D = 1.9; ( x ) Thouy and Jullien,
D 2.2: Thouy and Jullien. D 2.5; spheres, D = 3.0; (0)porous spheres,
D = 3.0; and Acrivos D = 3.0.
In the two extreme cases of vanishing or infinite double layer thickness, the poten-
tial uniform throughout the fluid, equal to 0 or respectively. The fluid
conductivity , (1 is then also uniform. Accordingly, the coefficient in
(114) tends to 0 and 1, respectively, these two limits.
The coefficient 6 defined by (1 14) is expected to depend on the four dimension-
less quantities K O , and We are interested in the limiting behavior when
and tend to 0.
Let us restrict the presentation of our numerical results to the only case that
could be completely worked out, 1. The first simplifying feature is that the
coefficient is found to be proportional to the solid concentration, as can be seen
in Fig. 13.
. Electrophoretic Mobility
Let us apply to a suspension a n electric field it will induce electrophoretic motion
of the particles with velocity This is used to define the reduced electrophoretic
mobility as
When the particle concentration is sufficiently small, v,, is equal and opposite to
the fluid seepage velocity.
Since computations are very long, the investigations on were also restricted to
Furthermore, the extrapolation to infinite cluster sizes was also impaired by
this limitation.
The reduced electrokinetic mobility defined by (1 19) is again expected to
depend o n the set of dimensionless parameters KG , and Its dependence
o n can be analyzed with the help of the calculations of and Neale
By use of a cell model, these authors generalized Henry's [7] prediction of the
mobility of a single spherical particle with radius a to concentrated systems.
Their result is cast in the form
f
Although Eq. (120) is not expected to apply directly to suspensions of fractal
aggregates, the function is a natural candidate to describe the dependence of
o n the concentration with a replaced by some characteristic radius
Thus the following form was tried:
It turns out that is close to a in all cases. even though ranges from
about 5 to 14. F o r a constant size 256, is even a slightly decreasing function
of This indicates that the particle interaction effects on mobility result mostly
from interactions between the of the aggregates and that the cluster pair
interactions are negligible. The value of R depends on the relative arrangement of
the particle within the aggregates. F o r 256, can be represented by the
linear fit
I. Introduction 259
I I . Transport Processes 26 1
A. Self-diffusion in niicroemulsions 26 1
B. Conductance of microemulsions 263
C. Transport across niicroemulsion acting as liquid membrane 272
References 275
Cosurfactant
Layer
D
D
P h a s e domain
greater than that of the AOT system. The higher in electrical conductivity
have been considered to be a consequence of the presence of elongated aggregates
in the system, inside which the charge can for a long distance without
hindrance.
2. Mechanism of Percolation
The threshold volume fraction of percolation is guided by the shell
length and the overall volume fraction of the dispersed phase. F o r zero shell length
and no interparticle attractive interaction, according to the close-packed
hard sphere 4, = 0.65; systems with strong attractive interactions end
with being lowered from 0.65 to 0.10.
During percolation, the transfer of charge occur in two ways: (1) by the
mechanism of "hopping" of surfactant ions from one droplet to another
and (2) by the mechanism of transient fusion, ion transfer, and fission
While Maitra et are in favor of the hopping mechanism, et
support transient fusion and the mass transfer mechanism. The dynamics of
the fusion model have been amply treated by use of the time-resolved fluorescence
quenching (TRFQ) method The formation of the bicontinuous state by
droplet association has also been suggested for augmenting percolation The
evidence for the existence of the globular form of droplets in concentrated water/
microemulsions given by the dynamic light scattering (DLS) method does
not support the necessity of a bicontinuous structure for the occurrence of percola-
tion Spherical droplets with strong short-range attractive interaction can
produce percolation. The hopping, transient fusion, and transfer mechanisms
are modeled in Fig. 3.
A thorough recent study on conductance percolation of water
and dodecane) systems with regard to and tempera-
ture is due to their observations on percola-
tion are in agreement with the reports of al.
According to Feldman et the conductivity and dielectric polarization
below the percolation threshold arise the motion of through
water channels and are not due to the hopping of ions. Above the threshold, the
same mechanism still dominates.
3. Theory of Percolation
The with ionic o r added electrolyte in the case of
can be considered equivalent to mixtures of conductor and
insulator. A t a constant temperature in the percolation range, the conductance
of the system on the whole is guided by the scaling law
= -
where k and are constants and is related to the conductance of the dispersed
phase.
According to a report of et the exponent is 1.9 for both
static and dynamic percolations. The scaling law for temperature percolation has
the form
= -
Fusion Fission
Transfer of ions Net transfer of
m ions
FIG. 3 (A) Hopping Ions hop in the direction indicated by the curved
arrows. (For simplicity, one ion is shown.) Ion transfer by fusion
and fission. cations are in each droplet; cations are transferred.
Scaling form:
= 1.5 - 0.333)
where is the conductance of the dispersed phase and the other terms have
already been defined. 0.333 is the threshold volume fraction for
dispersions 0.333. From mass balance and the rationale put forward by
and Fang and et structural information (particle
size and population, aggregation number of surfactant and cosurfactant per dro-
plet, etc.) on microemulsions can be realized. A set of typical results is pre-
sented in Table 1. The influence of surfactant, and cosurfactant in determining
the droplet size and amphiphile aggregation number is apparent from table.
Pal et simplified the E M T D D thoery of Granqvist and Hunderi for
dispersed particles of different states of aggregation:
For spheres,
=1 - 0.333)
For chains,
For clusters,
nent of -
tives either remained indifferent or increased the percolation threshold. The expo-
(1) was found to be 1.0, less than the expected value of 1.9. A special
droplet fusion model was proposed to explain the percolation-assisting effects of
the bile salts. In a subsequent study. Ray extended the additive effects on
the percolation of the water AOT heptane system, a unique effect of
sodium the conveision of a o
(water to a percolating one In addition to this. the scaling laws
were tested, and activation energies for percolation were evaluated. Like toluene
and xylene, benzene and naphthalene increased the percolation threshold; the aro-
matic rings are considered to block the droplet fusion process, thus hindering the
phenomenon. A representation of typical results is given in Table 2. It is seen from
the table that at 33 the threshold percolation temperature of ((without
additive) is lowered to by the presence of 0.1 whereas
of same strength does not affect the threshold temperature.
iodide (TEAI) at 0.1 dm 3 increases 9, to which decreases to
in the presence of of and 0.02 TEAI. toluene,
benzene, and naphthalene at comparable weight percents (11.5) increase to 43.0,
46.5, 46.5, and respectively. With 0.10 9, for 11%
xylene diminishes to From the analysis of the results, the values of and
of (1) and ( 2 ) were found to be much lower than the expected value of 1.9 in
the presence of additives. The modes of action of the bile salts and aro-
matics (blocker) on the droplet fusion-fission process are presented in Fig. 5. For
establishing the proposed model on a firm basis, further investigations are needed.
The dependence of the percolation temperature of a microemuision
on the concentration of triblock copolymers
and water droplets was straightforwardly suggested by Eicke on the
basis of a thermodynamic model. Two distinct ranges of the
droplet] ratio were observed to correlate the percolation temperature.
5. Miscellaneous Conductance Studies
et reported the unusual conductance behavior of 35 n-
propanol alkane and water SDS in the
presence of The conductance was found to decrease with increasing tem-
perature At lower the conductance decline has been to the difficult)
encountered by ions in hopping from droplet to dioplet At higher the con-
ductance increases due to the clustering of droplets in the
The binary mixtures of lower
Ill
Y
Surfactant 0 Counter i o n Na C
anion
In the recent past liquid membranes were employed for the separation and extrac-
tion of materials, and they can be conveniently employed for separating biological
materials Microemulsions of I and types
are considered dispersed liquid membranes that can the transfer of
soluble and water-soluble compounds, respectively, across them by trapping them
in microdroplets for convenient uptake and subsequent release. The
sions I and are called bulk liquid membranes. They are recent additions
in the field of separation science and technology. This field has been fundamentally
explored and advanced by Tondre and coworkers who worked out the
fundamentals of the transport process by studying the transfer of alkali metal
and other compounds across the They
also studied the transport of lipophilic compounds and
cene) across liquid membranes
3 Co-Ion-Dependent t of by lene
Glycol Dodecyl Ether Decane W 0
Transport
Anion
Acetate
Support from the Department of Science and Technology, New Delhi (GOI) and
the Indian Statistical Institute, Calcutta in the preparation of the manuscript is
gratefully acknowledged.
I. Introduction
I I . Chromatographic Materials
A. Chromatographic zirconia
B. Colloidal aggregates
C. Other materials
I I I . Diffusion
A. Definitions of diffusion
B. Porous media
C. Large solutes
IV. Conclusion
References
According to Fick's first law, the rate at which a species diffuses in a given direction
is proportional to the concentration gradient in that direction.
In Eq. is the diffusive flux of species A (mass or per unit area per unit
time), a vector quantity; is the binary or mutual diffusion tensor describing the
diffusion of A in a of A and B; and is the spatial concentration
gradient of A. If diffusion is isotropic, then it may be characterized by a scalar
value law is a phenomenological description of diffusion on a macro-
scopic scale. It is useful in the design and analysis of processes like chromatography
that involve nonequilibrium mass transport, as mathematical models of chromato-
graphy concern with how fast a solute penetrates into the stationary
phase, the flux.
Another definition of diffusion is given by the intradiffusion coefficient, which
quantitatively describes how far an average molecule A will travel with time
The fractional porosity accounts for the reduction in the flux due solely to the
reduced area available for diffusion. It is customary to use the fractional volume
porosity, but it has been argued tht the fractional area porosity at the material's
surface is more appropriate The definition given by (3) uses the tortuosity
to account for all other effects that alter the binary diffusivity, and it lacks any
inherent mechanistic o r predictive detail.
Here diffusion is measured only within regions that are permeable to diffusion.
the liquid-filled pore space, rather than spatially averaged over permeable and
regions as with (3). The tortuosity in is not identical
to in Eq. (3). As will be discussed below, partitioning and steric hindrance affect
the two differently.
When diffusing molecules are of a size comparable to that of a confining liquid-
filled pore, the measured diffusivity may be significantly less than that measured in
free solution due to frictional drag, steric hindrance, and differences between the
solute concentration within the pore and in the free solution (partitioning)
There are many quantitative predictions for rigid and flexible particles in variously
shaped pores Partitioning affects both equilibrium and
measures of diffusion. Because the average concentration within the pore is differ-
ent from that outside, the state of solutes within the pore differs
from that of the bulk solution in important thermodynamic quantities like activity
and partial molar volume measures of diffusion are further
reduced by steric hindrance. In addition to the area through which there is no
diffusive flux because it is occupied by the solid, there is a portion of the area
near the pore surface that is inaccessible to the center of mass of the solute because
of the finite size of the solute.
models also exist for the contribution of hydrodynamic o r frictional
drag in straight cylindrical pores. These are well reviewed in Ref. 2. If a solute
adsorbs onto the pore surface and spends some portion of its time immobile. its net
rate of transport will be reduced. Additionally, in some cases adsorption of the
solute can significantly diminish the size of the pore and may even block the pore.
preventing passage of another solute However, less is known about how to
model the hydrodynamic drag in the convex, multiply connected pores provided by
colloidal assemblies.
There is also the effect of the structure of the porous material. F o r a
equilibrium measurement of diffusion, one can consider that there is no straight
path for solutes to travel in the direction of the flux. In a n equilibrium measurement
of intradiffusion, this represents the fact that solutes are not longer subject to a
purely random walk. When a solute is near a pore surface, the probabilities for
moving in each direction are no longer uniform; certain directions are prohibited by
the pore wall. F o r technical precision, then, one should differentiate between a
structure factor and a tortuosity. A tortuosity, quantitatively describes experi-
mental results in which multiple interactions affect the diffusion. A structure factor,
q, quantitatively describes only the effect of pore space geometry and topology on
diffusion. Note that for limited conditions- when studying diffusion of small
and a passive pore surface- this allows for
where (fa is the fraction of solids =1 Equation (5) is valid for low
volume fractions of solids placed periodically the dispersed solid is treated as a
perturbation to the continuous liquid phase. Other than the required periodicity,
the structure is not specified.
Exact predicition of the diffusivity requires solution of the continuity equations
For an arbitrary structure, this is not feasible. However, for randomly placed,
overlapping spheres, a statistical argument provides a closed-form analytical solu-
tion as a limiting case to a more general series solution
For periodic structures, the continuity equations may also be solved explicitly by
expansion to yield tractable series solutions T o a first approx-
imation, these treatments are equivalent at low volume fractions of solids.
Experimental studies of diffusion of small nonadsorbing solutes within ordered
assemblies of colloids is has been well predicted by these treatments
However, these theories d o not consider solutes of nonzero size, and another
method is necessary to examine this.
A random walk simulation thus provides the structure factor for any geometry that
one defines. This technique has been used to examine the structure factor for an
infinitely solute within random arrays of freely overlapping cylinders with
various orientation distributions to study the effect of the distribution of
reactive sites within a porous catalyst and to study the effect of solute size
on the structure factor within ordered assemblies of spherical colloids
While the random walk simulation is straightforward in concept and in imple-
mentation, it is limited in application. The simulation does not contain time or
energy; thus forces cannot be incorporated, and solutes and pore surfaces cannot
interact other than to prevent overlap. For hydrodynamically interacting Brownian
particles, the apparent that would be measured at long times is
expressible as the product of a hydrodynamically determined short-time diffusivity
and a contribution that accounts for the distortion of the equilibrium dispersion
structure caused by the diffusing particle The short-time diffusivity, is an
ensemble-averaged quantity that accounts for hydrodynamic forces between part-
icles and can be evaluated from the hydrodynamic diffusion tensor using Monte
simulation It has subsequently been postulated that diffusion in a
porous medium can be described analogously as a product of a factor accounting
for structure, a factor accounting for hydrodynamic forces between diffusing
solutes and the pore and the diffusivity of the solute in a solution
The friction coefficient f has been determined either from the hydrodynamic
diffusion tensor using Monte simulation or by using one of several
simple analytical expressions for a sphere translating through a right cylindrical
pore 68-70]. The most commonly used is (11) where is the
characteristic ratio of solute size to pore size
T o aid in the design of the aggregates we make in our laboratories, we have sought
to understand the diffusion of large molecules within these porous materials as well
as other materials. N M R and light-scattering techniques have allowed researchers
to directly measure diffusion in porous media that were previously inaccessible to
experimental measurement due to small length and time scales. We have summar-
ized here the framework providing a general relation between the measurement of a
diffusivity and its use in a design equation. Hydrodynamic theories
and computer simulations have led to significant advances in quantitative predic-
tion of large molecules within porous materials that possess a periodic or statistical
regularity, which we have applied to our own materials made of aggregated sphe-
rical colloids; however, these results have underscored the necessity of having at
least a qualitative knowledge of the pore geometry-generally related to the
r
sythesis method.
This work has been supported by grants G M 45988 the National Institutes
of Health and CHE 917029 from the National Science Foundation. Additional
support has come from the Training Program of the National
Institutes of General Medical Sciences the Institute for
Advanced Studies in Technology at the University of Minnesota,
and the Regional Instrumentation Facility for Surface Analysis at the
University of Minnesota.
M. I. and
Unger, Adv. (1997).
Kresge, Roth. Leonowicz,
S. Hellring. Beck. Schenker, Olson, and
Sheppard, Chem. (1994).
A. Kurganov, Unger. and J. A (1996).
F. and Science (1996).
A. Einstein,
Leipzig. 1922.
Diffusion: Mass Transfer in Systems, Cambridge Press,
Cambridge, UK, 1984.
Chem. (1992).
L. and N. A. Sci. 49-41 (1994).
Foster, J. Butt, and Bliss, Catal. (1967).
E. P. Russell, and N.
(1968).
and (1969).
(1981).
Nitsche and K . W. Res. (1994).
Anderson and J. H. Brannon. (1981).
J. A. I. P. Petrou, and
Bioseparation (1990).
Flow and transport in chaotic media, Thesis,
Minnesota, Minneapolis, MN, 1985.
and Chem. (1990).
Johnson and Stewart, Catal. (1965).
P. Friedinan and A. Chem. Eng. (1995).
D. Keffer, A. and Davis, (1996).
Maxwell. A Treatise on Electricity and Magnetism, 1. Clarendon Press,
Oxford, 1873.
B. and P. Colloid
(1986).
Bird, E. Stewart, and E. Lightfoot. Transport Phenomena. New
York.
(1963).
P. R. Leyte, and Colloid Interface Sci.
(1991).
and Leyte, Colloid Interface (1994).
I. Introduction
A. Real and ideal surfaces
B. Origin of the adsorptive potential
C . Thermodynamic meaning of the adsorptive potential
I I . Theoretical Background
A. A short history of the integral equation of adsorption
B. Solving and using the integral equation of adsorption
III. Application of Density Functional Theory
A. The deconvolution method
B. Parameters of the model matrix
C . Results
IV. Discussion
V. Conclusions
References
In the formative years of adsorption theory. it was usual to assume that a solid
surface presented a set of equivalent sites for the attachment of a n adsorbate
molecule. Such was the idealized picture used by Langmuir [2] in developing his
classic equation. Even when the model was extended to include
adsorbate interaction as was done by d e Boer [3] in developing his nonideal
where is the potential. the depth of the potential well. the molecular separa-
tion. and the separation at which = 0. For an adsorptive molecule at the
surface, the total dispersion potential, is the of its interactions
with each atom of the solid:
the preceding section, we presented the idea that real surfaces can be, and usually
are, "energetically heterogeneous" in that they provide regions of varying attractive
force o r potential to physically adsorbed molecules. In describing a surface energy
characteristic. a n energy function that is independent of temperature and concen-
tration is desirable. Such a quantity is conveniently defined by the potential energy
difference between the lowest energy state of the adsorptive molecule in the gas
phase and its lowest energy state in the adsorbed phase. both at infinite dilution. I n
this way we obtain a quantity that is independent of the kinetic states of the
molecule in either phase and that most directly the adsorptive potential
of the system. The adsorptive potential, is therefore defined for a homotattic
surface p. as
where the term - is the thermal vibrational energy per mole of adsor-
bate in excess of the zero-point energy, and are the kinetic energy
changes on adsorption due to loss in translational and rotational degrees o f free-
dom, and the energy of interaction with all neighboring adsorbed molecules.
This latter term is clearly dependent on surface concentration and could be esti-
mated by a n equation analogous to
Another thermodynamic quantity, the heat of adsorption, can be
calculated from the temperature dependence of the adsorption isotherm by obtain-
ing the slope of the adsorption and is defined by
-
where is the moles adsorbed and V and V are the molar of the gas and
adsorbed phases, respectively. Neglecting the molar volume of the adsorbed phase
and assuming that the gas phase is ideal gives the more familiar expression
It can be shown that the differential and isosteric heats are related by
Although the concepts are somewhat older, the widely used model for describ-
ing adsorption on an energetically heterogeneous surface was first explicitly stated
by Ross and The postulates that the surface of a real solid is
composed of small patches of different adsorptive potential that adsorb indepen-
dently of one another. The distribution of adsorptive energies, among these
patches can be represented by a continuous distribution function:
is the patch (or site) frequency per unit energy interval on a surface of area
A. The distribution function must normalize to unity, as was pointed out by Hill
since we are dealing with a surface of finite extent; that is, = 1 over the
range of energies considered significant,
At any equilibrium pressure under isothermal conditions, the quantity
adsorbed per unit area, on a given surface patch will depend only on the adsorp-
tive potential of that patch according to the function
The observed total amount adsorbed. a t pressure is then the sum of contribu-
tions from each patch of surface,
Equation (9) is therefore the general form for any adsorption isotherm and corre-
sponds to (IV-4) of Ref. The function is called the kernel function o r
the local isotherm. Equation (9) is now often referred to as "the integral equation of
adsorption".
Localized no interaction.
equation
over the whole range of both positive and negative. the integral yielded the
equation in the
,
In these early efforts, the goal was to find an analytical description of using
a specific theoretical model for adsorption on a surface and a n assumed
distribution function. An alternative way to use Eq. (9) is to assume the of
both and and thence derive the distribution by inver-
sion. That this can be done and that a unique solution exists for was demon-
strated rigorously by Sips T o describe the experimental isotherm, Sips chose
for an empirical equation of the form
= +
For the solution of this equation, a new set of experimental results are demanded,
namely the differential heat of adsorption, as a function of Q. Equation (18)
still contains the unknown function This function can, however, be elimi-
nated when we know as a function of Q at absolute zero. At that temperature,
every surface patch is filled in the serial order of its adsorptive potential, beginning
with the greatest. A surface patch is therefore either completely filled or completely
The function can be replaced with the Dirac delta function,
times called in this context the condensation isotherm. This results in
Diameter
. (a) Adsorption of nitrogen a t 77 K on a catalyst reference
(b) area by using Jura
thickness curve.
Relative Pressure
(a) Distribution of for a silica-alumina catalyst support from the nitrogen
isotherm at 77 (b) data showing the M D F T fitted result.
the whole surface can be calculated from the reported distributions of using
(24) and (25) together with (16):
where a, is the area of the patch. The results of this calculation are shown in
Table 2, together with total surface areas determined with the present method. The
BET values are included for reference.
The value of for Sterling from this work, corresponding to
2.15 compares very well with the best estimate of Ross and Olivier
p. of 2.16 We see from (25) that the surface parameter intervals
of correspond to about
Ross and On
123.
2. Langmuir. 38: 2267 (1916); Phys. Rev.
R. J. Phys. (1950).
and
A. E. Trans. (1951); L.
Drain and A.
Rudzinski and Everett, of
London,
Ross and Surface (1975).
Morrison and Ross. Surface (1973).
and Morrison. Colloid Interface (1979).
Porous Mates,
in tire
1996,
A. R. Walton, and Quirke, (1989).
and (1993).
and Quirke, (1993).
Conklin. and in Science
87, of Porous Solids et
Tarazona, A (1985);
and Mo!. (1987).
and Solving
I.Introduction
II.Algorithms
III.Modeling the Interactions
IV. Specific Examples
A. Elementary systems
B. Lennard-Jones atoms o n various surfaces
C. A d s o r p t i o n the critical region of the bulk fluid
D. Freezing of adsorbed layers
E. Energies of adsorption
F. Wetting and drying
V. Discussion
References
dynamic property of interest is the energy of the adsorbed phase o r , more specifi-
cally, the molar energy of adsorption from the gas. If this gas is ideal, the
energy change will be entirely to the potential energy acquired by the adsorbed
molecules when they are transferred from the bulk gas to a location that is a t or
near the adsorbate surface. This energy is to the interactions of the adsorbate
molecules with each other and with the solid, and these two contributions can easily
be evaluated separately in the simulation (and summed to give the experimental
energy of adsorption). The molar energy of adsorption is related to the enthalpy
change by an additive factor However, it should be noted here that
these energy and enthalpy changes are not the partial molar quantities usually
reported in experimental studies. F o r example, the heat is the partial
molar enthalpy change upon adsorption. If and the molar and partial
molar energy, respectively, then is equal to + Thus, to evaluate
the partial molar quantities from the molar quantities, one must either simulate the
energy to a sufficiently high precision to allow one to differentiate it with
respect to A' or take the alternative route of evaluating the partial molar energy
from the appropriate fluctuation in the GCMC ensemble. Theory gives
where can denote or and and are the lower and upper limits of the
position variables that define the adsorbed layer. Ordinarily, these are taken to be
zero and for where 1 is the thickness of the adsorbed layer, and and for
and the in-plane size variables. The brackets indicate an average over the
configurations generated in the simulation. Note that the in leads
to an integral equation for the number density in the fluid (Ref. 7. 3.131).
Chemical equilibrium in an fluid requires that the local chemical
potential be independent of F o r a fluid in a canonical ensemble (either
molecular o r Monte one considers a particular configuration of
molecules. An additional molecule is inserted a t this molecule interacts with
the molecules and the solid initially there but does not affect the simulated
molecular configuration in any way. If the energy of this interaction is denoted
by N), where the added particle is at and the A' particles are at positions
denoted by it can be shown that
where is the density of the fluid a t point , The average is taken over many
computer-generated configurations of the adsorbed molecules in ensemble.
Insertion points can be chosen a t any point in the fluid. chemical potential
of the fluid must be independent of location, o r there could be a spontaneous
transfer of molecules from the with higher chemical potentials to those
with lower, which would mean that such systems would not be at the lowest global
chemical potential and thus would not be in thermodynamic equilibrium. This
method of simulating has been used successfully for bulk model fluids,
including nonspherical molecules where the coordinate includes as
well as translational variables. the literature, one can see definitions of quantities
such as the electrochemical potential for a fluid in a n external electric field or the
intrinsic pressure in which the contribution of the gas-solid force is omitted from
the definition. These quantities will depend on position, but one still has the
requirement that their values plus the contributions due to the external potential
must be constant.
The expressions for the pressure tensor a t a surface can be integrated to give the
solid-fluid surface tension, which we here denote as (also called the spreading
pressure):
where is the volume of adsorbed fluid. This integral can be evaluated analytically
to give
ALGORITHMS
The canonical Monte method has been altered to give averages in any
of interest. The basic idea is to generate a chain of Markov states (mole-
cular configurations) for the fluid that will eventually produce configurations that
characterize the desired statistical ensemble. T o do this, one specify the tran-
sition matrix a that takes the system from state to state by making a trial
change in state and then testing to see if the change leads the toward one
with the desired properties for the ensemble that is to be simulated. This means that
of the trial changes will be rejected, and one of the goals of the algorithm is to
maximize their acceptance probability. It can be shown that this probability
is where is a quantity that depends on the parti-
cular ensemble being simulated. To accept a trial change, a random number
between 0 and 1 is generated. If the random is less than
the trial change is accepted (note the trial is always accepted if is negative),
and if it is greater. the change is rejected. Even for acceptance ratios on the order of
0.1 or more, a chain of changes can contain upward of trials, depending on the
accuracy desired.
To generate a canonical ensemble, all the trials consist in changes in a position
variable for molecules in the computer box, and 5 is the change in the energy of
the molecule that is moved due to its interactions with the other molecules in the
fluid and with the solid. The scale of these moves is adjusted to be small enough
that 8 is not always large and positive so that molecules do not overlap too
much in the new state) but large enough that the phase space of the system is
explored in a reasonable number of trials.
Ensembles other than canonical ones can be generated by adding to the types of
moves that define the Markov chain. In the case of particular interest here, the
where and are the well depth and size parameters for the pair of sites
under consideration. It is the choice of values for these parameters that is the
most serious problem in devising a potential to represent a specific system such
as a gas-phase interacting with carbon black or silica (two of the most
popular adsorbents in the experimental literature), whose bulk phase properties
have often been simulated using this kind of model. The problem is that one has
cross interactions for the adsorption energy functions with potential parameters
for the sites in the solid that are unknown. Even if one were to guess these
parameters (for C-C or 0-0 interaction, for instance), approximate combining
rules such as the Lorentz-Berthelot rule are needed. More often than not. the
validity of the results obtained is tested by comparison of experimental and
calculated low-coverage limits of the isotherm slopes and heats of adsorption,
and, if necessary, the parameters are adjusted to give agreement. [Low-coverage
limiting data are used because theory shows that in this limit the calculations
involve only isolated molecules on the surface (see Ref. 15, At
first glance, the summations over the semi-infinite numbers of sites in the solid
would seem to be a sufficiently formidable task to prevent one from successfully
carrying out simulations based on such models. However, there are ameliorating
factors: In the first place, the Lennard-Jones functions are sufficiently
ranged to allow one to truncate the sums after the nearest 100 (or so) sites
have been included; for higher precision, the long-range parts of these sums can
be replaced by integrals. Also, the atomic arrangements at planar surfaces of
crystalline materials that are made up of a single Miller plane are periodic and
thus will give site-site gas-solid interactions that have the same periodicity as
the exposed crystal face. This can be exploited to give useful expressions for the
interaction of a gas molecule site over a surface unit cell 2.2; The
two-dimensional vectors that define the size and shape of this unit cell can be
defined as and with reciprocal lattice vectors and such that
. = 1, . = 1 and 0, . = 0. If the vector connecting the
gas-phase atomic site and the center of the unit cell in the layer of solid
atoms is denoted by where denotes the orientation of the in-plane
vector one can show that
In this expression, is the angle between the vectors and the K,, are modified
functions of the second kind, and where and are
integers. Here we have given a simplified version of the theory for the case where
there is only one atom per surface unit cell and the solid atoms are all of the same
kind. Extensions to deal with these problems are straightforward but give rather
complicated expressions. T o obtain the total energy of interaction, these energies
should be summed over all values of a , but the sums converge very quickly and
convenient approximations to their values can be derived. A convenient approx-
imation to the value of the term with g 0 in (7) can be written
the ideas presented here, since few if any porous materials exhibit internal surfaces
u p of single crystal planes.
Besides its accessibility for experimental work involving large-area adsorbents,
one of the reasons for the great popularity of materials whose adsorbing surface is
made u p of the graphite basal plane is the great simplicity of their surface. This
popularity extends not just to the free surface but also to porous materials such as
parallel-walled slits and even to the buckytubes that are presently under intense
scrutiny. the first the atomic spacing in the graphite basal plane
is so compared to adsorbate site sizes that the periodic component of the
solid energy can be completely neglected. If one keeps only the integrated site-site
energies that lead to the inverse 10-4 power part of the gas--solid energy for the
infinite plane, interesting modifications are available; in particular, it is possible to
derive a n analytic expression for the interaction of a gas a t o m with a truncated
graphite basal plane that is useful in a number of problems F o r example,
stepped graphite surface can be modeled pores that are more realistic than the
parallel-walled case can be treated, including those rectangular o r triangular
cross sections Although the grooved graphite surface that uses truncated
planes to make the grooves in the outermost surface layers might be considered
be of minor interest, one could take two such surfaces facing each other to make
parallel-walled pores with alternating narrow and wide regions with constric-
tions).
Finally, the planar graphite surface can be rolled into a cylinder to model a
The integrations that led to a 10-4 function for the infinite plane are
difficult in case, but a n analytic expression for the site-wall interaction in
a smooth cylindrical tube can indeed be derived for use in the simulations of
adsorption in this material
Well-known and widely used adsorbents that are amorphous include many
preparations of high-area oxides such as silica and titania. Models of the
adsorbing surfaces of these materials have been proposed that are based on
the ideas that first, the cations can be omitted for the site-site gas-solid poten-
tials, at least for nonpolar adsorbates, and second, that the arrangements of the
oxides is essentially that of a random close-packed array of hard spheres. An
algorithm is available that will generate such arrays with periodic boundary
conditions in two directions. If the hard spheres are then transformed into
sites with interaction that represent the oxide-gas site potentials,
the nonperiodic surface can be taken to be the adsorbing plane. In this case, the
brute force approach of summing over these oxides can be used to evaluate the
gas-solid interaction for such amorphous materials. The results of this modeling
have been analyzed in a number of papers (References to earlier work
are given in Ref. 40.) In addition, variants o n this procedure have been
explored: One can roughen the surface by randomly deleting atoms from the
layer o r o r pores can be created by deleting the oxides from a
of appropriate size and shape in the computer-generated random solid
This allows one to model pores of arbitrary shape and even pore net-
works of varying connectivity (if a computer of sufficient power is available to
handle such large systems).
By "hard interactions", one means interaction energies that jump from zero to
infinity at a specific separation distance or where gg denotes gas-gas and
gas-solid. Obviously, such systems have n o interaction energies that could give
energies of adsorption that differ zero, regardless of the geometry of the
confining wall. Nevertheless, the presence of the wall o r walls produces profound
changes in the properties of a fluid of hard spheres in the vicinity of the walls.
Adsorption occurs with a n adsorbed amount denoted by the excess number of
(relative to the bulk) per unit area of surface when the pressure of the bulk
fluid is This is always greater than zero, which means that the hard sphere fluid
always wets the wall in these A t first this appear to be a surprising
result. because there are n o attractive interactions between the wall and the fluid,
but of course these interactions are absent in the adsorbed fluid as well. The number
of atoms adsorbed per unit area is given by the general expression
where denotes the density of atoms per unit volume and adsorption on a smooth
planar surface is assumed. are needed to evaluate the number
density in the fluid at a distance from the wall. Figure 1 shows the results of such a
simulation for a high density hard sphere fluid in a large hard-wall box with
walls apart. In this case, the density of the fluid far from the walls (and
thus the bulk fluid density) is = 0.90. W e see in the bottom panel that is
essentially equal to the bulk density for 5, which means that this system
behaves as if it were adsorption o n a free surface.
Specializing to planar walls for the one has the exact relation that
= where is the local density of the adsorbate in contact with the wall
when the pressure of the hard sphere fluid is F o r hard sphere mixtures
is the sum of the individual densities for each of the components in the fluid. Thus,
the pressure of the fluid can be obtained from estimates of the intercept of the curve
of versus for example, Fig. indicates that is between 8 and 9 for
this system. This result, taken together with the calculation of a t a given from
allows one to construct the isotherm Figure 2 shows the adsorption
of a hard sphere fluid on a hard wall as a function of the bulk-phase density
The simulation points compare well with results of two theoretical calculations
based on the scaled particle theory.
One also has a n exact expression for the hard wall-hard sphere fluid surface
tension which is
where is the separation distance between a pair of hard spheres, is the pair
correlation function for a pair of fluid particles, both of which are in contact with
the wall, is the pair correlation function for atoms 1 and 2 in the fluid).
A second general route to the surface tension (or spreading pressure) has been given
in I as the integral of difference between the components of the pressure
tensor parallel and normal to the surface. The local values of both of these pressure
tensor components for a bulk-phase fluid density that is close to freezing are shown
in the top panel of Fig. 1. The perpendicular pressure tensor element is independent
of as required for mechanical equilibrium. The transverse component (which is
not uniquely defined even though its integral with respect to z is oscillates and
approaches at large as it should. A11 integral over these functions gives the
gas-solid surface tension [see Eq. The alternative route to this surface tension
via (1 requires a simulation of the pair correlation function for spheres in
contact with the hard wall. Figure 3 gives a between this function and
that for the bulk fluid, and Fig. 4 shows surface tensions as a function of the bulk
fluid density.
Finally, a number of pair correlation functions with = but with varying z
are shown in Fig. 5 and compared with the bulk pair correlation functions for
. 3 The pair correlation function for hard spheres at a hard wall plotted as a
function of the reduced separation distance between the pair of spheres for a reduced
density of 0.90. The points show the simulations; the curve is for the bulk hard sphere
fluid at the same reduced density. Ref.
same density. (The density profile for this is shown in Fig. 6.) In general,
simulations indicate that sphere pair correlation functions near a hard
the two equal are to those for the bulk fluid. Figure 7
shows the dependence of the pressure of the hard sphere in a slit pore o n the
spacing for a fixed value of the density oscillatory
behavior of the pressure is to the layering the the
The minima occur easily a n integral
of layers. Maxima are produced integral of layers cannot be
readily Note the periodicity function is slightly larger
The phase diagram for hard spheres enclosed in a hard-walled slit excess
ranging a value one layer of (or solid) = 0) to
layers shown as a of the sphere Phase boundaries and
coexistence regions are indicated. nature of the ordered phases observed is discussed
in Ref. 72.)
effort. The width of the loop clearly depends on temperature, pore geometry, and
the gas-solid interaction potential. This phenomenon is a manifestation of meta-
stable behavior. Densities a t the tops and bottoms of the hysteresis loops may be
considered as corresponding to points o n the temperature-dependent lines
for the pore fluid, and the equilibrium densities can be obtained by calculating the
in the
If one measures o r simulates adsorption at temperatures near the critical region of
the bulk adsorbate, new features appear in the isotherms, for adsorption both on
the free surface and in pores At slightly temperatures
where bulk-phase condensation is absent, multilayer formation can be observed, as
is shown by the local densities in Fig. 16. This simulation was for a simple model of
ethylene a graphitic surface (chosen so that comparison with experiment
was facilitated) at a temperature roughly 3 % higher than the bulk critical value.
N o t surprisingly. the density for the are not sharp as they are at
lower temperatures, and these layers build u p gradually compared to the low
temperature behavior. In addition, the simulated isotherm in Fig. 17 shows the
characteristic cusp in coverage versus pressure that is observed The
position and height of this cusp in the excess amount adsorbed depend o n the
choice of dividing plane that defines the adsorption excess; its existence is believed
to be due to a slow increase in absolute amount adsorbed compared to the more
rapid increase in the density of the bulk fluid increasing pressure, thus yielding
adsorption excesses that decrease with increasing pressure. The analogous simula-
tion for a graphitic slit pore shows the same phenomena. but with quantitative
features that depend o n slit width. In a more recent simulation the critical
behavior of fluid adsorbed in a graphitic slit pore was carefully studied. Although
not in contradiction to the earlier simulations. this high precision study revealed
that the pore fluid density for a rather large wall separation = 20) did not
approach the bulk-phase density in the core region of the pore far from the
walls). The explanation suggested for this is based on the observation that the
compressibility for the confined fluid was not the same as that for the bulk for a
temperature very near probably due to the fact that the capillary critical
temperature is slightly different from that for the bulk even for the large slit
17 Simulated excess adsorption in atoms per unit area (in reduced units) as a
function of pressure = for an adsorbed gas slightly above its bulk critical
temperature. The points show the adsorption for two slightly different definitions of the
excess adsorption (From Ref.
At low temperature and high density, an adsorbed fluid with attractive interactions
such as those for the Lennard-Jones potential will freeze. We have already dis-
cussed the ordered phases that have been simulated for hard sphere fluids in slit
pores with small wall spacings. The 2 D Lennard-Jones fluid also forms an ordered
solidlike phase at high density, so there is every reason to expect that such phases
will be formed from adsorbed fluid. Indeed, this is the case. Experimental and
simulation studies of this behavior have been reviewed elsewhere and
need not be repeated here. It is worthy of note that the periodic boundary condi-
tions used in simulations for fluids can give problems when a solidlike phase with
long-range order in the directions parallel to the surface is Unless the
computer box happens to match perfectly with the lattice of the adsorbed
phase, there will be difficulties in matching up the images with the real sample. A
possible way of avoiding this problem is to simulate in the isobaric ensemble (for
the pressure components parallel to the surface), since this causes the area of the
computer box to fluctuate in order to hold the parallel pressure constant.
Nevertheless, questions remain concerning the correct expression for this pressure
in the presence of an adsorption potential that varies with position parallel to the
surface. A more reliable but more expensive method is to simulate a large isolated
patch of material on the periodic surface Simulations of the effects of periodic
variations in the gas-solid potential energy on the freezing transition d o reveal
some interesting features. Many of these simulations have been done with mole-
cular dynamics because of the power of this algorithm in revealing structures even
at low temperatures where Monte simulation becomes difficult. However, we
have chosen not to cover molecular simulations in the present review on
the grounds that this will be done in the near future
Starting with an early paper in which it was shown that an argon
layer on a graphitic surface at was ordered while the system in a
considerably weaker inverse 9-3 potential was not (see also Ref. a study
of the argon-graphite at 77.5 K with and without the periodic terms in
the potential revealed the periodicity had only a minor effect on the incom-
mensurate monolayer. This is definitely not the case for krypton-graphite,
where a commensurate layer does in the presence of the periodic potential
It is clear that the difference arises from the different sizes
of the two adatoms. Krypton with = 3.60 A , = 4.04 is the distance
for the minimum in the gg pair potential) matches up rather well with the
mensurate lattice spacing of 4.26 A , whereas for Ar is 3.82 A. In fact, the solid
K r monolayer undergoes an transition driven by
thermal expansion a t a temperature that depends on the magnitude of the periodic
part of the potential. Simulations of atoms on surfaces made up of exposed
As the holding potential for a gas o n a solid surface o r in a pore becomes weaker
and weaker, a point is reached where the competition between adsorp-
tion and the condensation of the adsorbate gas into a bulk liquid (or solid) is no
longer won by the adsorption process. I n effect, this describes the changeover from
wetting to In this picture, drying is a special case of adsorption
and can occur when the potential energy in the bulk phase is sufficiently negative
to that in the to stabilize the bulk relative to the adsorbed
film. Since these two energies change differently with temperature and density. there
is the strong possibility that marginal systems will undergo transi-
tions along some line in the temperature-density plane (evidently, other thermo-
dynamic variables can also be used to characterize the transition). The nature of
this transition (first-order o r continuous) is a question of considerable interest. A
nonwetting film comprises a thin layer (usually. a fraction of a in
equilibrium with the bulk fluid, and the transition to the wetting film corresponds
to the formation of a layer of fluid o n the surface whose thickness
infinity as approaches Other transitions are possible such as thin
film (but not infinitely thick) at a less than which is called
A study of argon o n a weakly interacting surface (an inverse 9-3 potential
with parameters characteristic of solid gives the local density curves for
= 1.1 that are shown in Fig. 19. The for this system, given in the
inset of the figure, indicates only small adsorption u p to a value of the bulk density
of about 0.048 (in reduced units) in equilibrium with the adsorbed film. When the
bulk-phase density exceeds this, a pronounced change in the nature of the local
density of the adsorbed phase occurs that is reflected in the jump in the amount
adsorbed. Thick layers appear, and the density of the monolayer is much greater
than that found for the thin films. Such behavior corresponds to a
transition. In contrast to this, Fig. 20 shows that the isotherm and the local den-
sities for this system a t = 0.95 grow continuously until saturation is
reached. However, the slope of the isotherm at saturation is finite, so that the
coverage is also finite a t saturation, indicating partial wetting.
I n this review, we have tried to show three aspects of the computer simulations of
condensed matter:
The study of extremely simple systems, as exemplified by hard spheres. In this
case. the computational power needed is relatively small, and the results are of
primary interest as tests of theories. They are a t best in qualitative agreement
with experiment.
I. Introduction 355
I I . Theory 357
A. Principles of the absolute rate theory and of adsorption
equilibria 357
B. The condensation approximation approach to adsorption
equilibria and kinetics 359
C. Generalization of A R T for energetically heterogeneous solid
surfaces 364
D. Theories relating the rate of adsorption/desorption processes
to the chemical potential of adsorbed molecules 374
References 388
describing the dependence of the amount adsorbed on time Here, a and are
constants characteristic of the adsorption under investigation. Assuming
that at = 0, = 0, one arrives at the following integral form of the
equation:
where is a constant and G(9) is another unspecified function of the surface cover-
age 0. The change of the surface coverage with time is thus given by
The absolute rate theory does not lead to explicit expressions for the functions
and G(0) and does not suggest that might depend o n the surface
coverage
First a coverage dependence of the preexponential factors was introduced o n a n
empirical basis, leading to the following expression for
for 5
elsewhere
is commonly applied to represent the dispersion of adsorption energies. Its applica-
tion in Eq. (9) yields a n expression for similar to the well-known
where
Thus, is the value of a t which the isotherm derivative reaches its max-
imum. Figure 2 shows the derivatives of the isotherms drawn in Fig. 1 .
Surprisingly, the above-described practice leads to pretty exact values of 9, for
adsorption systems (adsorption regimes) reported in the literature. The rea-
son for that is that appears to be the leading of a certain expansion for
a t finite temperatures F o r a better understanding of the forthcoming consid-
erations and the approximation made, we repeat briefly the principles of that
expansion.
The integral of (9) is evaluated by parts:
It can be shown that the first term on the right-hand side of (18) is negligible
compared to the second except for very small and very high values of 9,
Therefore the second term is evaluated by expanding into its Taylor series
around the point = at which reaches its sharp maximum,
The first term on the right-hand side of (20) is just the result obtained by
applying the condensation approximation, assuming the hypothetical limit
0. Because of the symmetry of the function the second term on the
right disappears, and the first nonvanishing correction to is the second
correction term. F o r the model of adsorption, to a good approximation
is not the only function that may appear in the theoretical considerations of
adsorption o n heterogeneous solid surfaces. However, one assumes that
as, cis,'
is usually called the condensation approximation for the actual adsorption energy
distribution and is defined as
The concept of precursor states and the corresponding expression for the rate of
desorption did not gain the attention of researchers using TPD. Nor did it solve the
problem of the large variations in the preexponential factor with surface coverage.
Seebauer et reviewed the existing theoretical interpretations of
factors; none of was able to explain the large variations found in the
analysis of experimental data.
The procedures used to generalize the A R T expressions for the case of energetically
heterogeneous solid surfaces are discussed in the excellent review by
The rate expressions were written in the following For adsorption,
[
exp --
Contrary to what is found in the analysis of TPD spectra based on ART, the
preexponential coefficients and were considered to be independent of surface
coverage and of adsorption energy Also, the assumption was made that the
activation energies for adsorption and desorption, and are linear functions
of the adsorption energy
where is the value of the equilibrium pressure at which the surface coverage is
equal to This is the same pressure that appears in Eq. (34). So we replace
in (35) by its value calculated from
Replacing
- 1
by the expression in
+
+ 1
we obtain
+ T
= exp
! -
Thus, in the limit of a strong surface heterogeneity leading to the Temkin equi-
librium isotherm, the of a strongly "activated" adsorption leads us to
the rate expression
- 6) = constant
Because at = 0, 9 = 0, we have constant = 0. so that
In the CA limit T 0.
-10 -5 0 5 10 15
The behavior of the function defined in (44) for three values of the
parameters = = 0.1 - -), and =5 -).
Equations (50) and (51) are the integral and differential form, respectively, of the
Elovich equation.
Note that when linear relationship (32) between and is assumed, Eqs. (50)
and (51) still preserve the form of the Elovich equation. This raises the question as
to why along with the assumption that = and as a rectangular
function, does not lead to the Elovich equation in the limit
- This inconsistency has not been discussed to date. It seems
that the source of this problem lies in the use of the Langmuir equation
T) with the ART expression for adsorption kinetics.
We note here that there are other theories of kinetics that
offer expressions for the rate that are different from the
ART expression but also lead t o the Langmuir isotherm when = 0. It is
rather strange that adsorption systems with different kinetics of the adsorption/
processes have the same form a t equilibrium.
Let us consider the situation when adsorption runs at quasi-equilibrium condi-
tions. This means that a t a surface coverage all the surface correlation functions
In Eq. ( 5 3 ) . is the partition function for the vibrations normal to the surface,
is the vibrational partition function for vibrations parallel to the surface, and
stands for rotations, vibrations due to hindered rotations, and electronic and
spin degrees of freedom. Equation (52) can now be rewritten in the form
where
When the surface is heterogeneous, will be related to the total surface cover-
age and (54) will give the "local" surface coverage of sites having adsorption
energy The features of the function will be the same as those of function
(1 1) shown in Fig. 1 . F o r instance, the derivative will have its value
at 9 =
Meanwhile the derivative of the function defined in (44) reaches its
maximum value. at = = 0.368.
Although a comparison of these two pairs of values does not show dramatic
differences, it points to a n inconsistency when combining A R T kinetic expressions
with the Langmuir model for adsorption equilibria. W e can always imagine the
adsorption process running so slowly that conditions can be
assumed. However, even stranger is the necessity of invoking the existence of an
activation energy for adsorption and its variation on various adsorption sites. In
processes involving chemisorption, a n activation energy for adsorption can be
expected. However, an activation energy of adsorption processes of physical
adsorption is less certain. W e will show that the apparent activation energy for
physical adsorption is the result of neglecting the readsorption term in the A R T
expression, (6). Since the is obtained from a fully
rigorous consideration of adsorption kinetics should take into account both the
adsorption and desorption
Integration of (6) now yields
and
where
W e next accept the assumption that the activation energy for adsorption does not
exist, = 0, and that is equal to the isosteric heat of adsorption within the
accuracy t o a constant. Then,
Figures 5 and 6 show the behavior of the function defined in (59) for
various values of the parameters = and =
Figure 7 shows the behavior of in Eq. (59) a t three different temperatures.
We can see in Fig. 7 that once again we could use the condensation approximation
We see that (17) cannot now be solved analytically to obtain the function
In view of the fact that in (59) does not go to unity as Eq. (14)
has to be modified by replacing unity o n the right-hand side by the maximum value
of 8, which is reached in the limit Now (14) takes the form
sorption kinetics d o not play an important role. What is important here is that we
have arrived at the "Elovich behavior"' without invoking the concept of activation
energy of adsorption.
In this way we have eliminated one in developing the Elovich equa-
tion on the basis of ART. However, there are other inconsistencies in the general-
ization of the A R T approach to suit the case of heterogeneous surfaces, using
(30) o r (31). Namely, the assumptions underlying the ART approach d o not involve
a particular mechanistic model of adsorption, for instance, the one leading to the
isotherm. To be fully consistent, one has to use expression (44) for 0 in
(30). In using the Langmuir equation for 9, one uses, in fact, the function
defined in which is the limiting form of 0 defined in (59) when t
and is obtained in addition by integrating full A R T expression, (6). Thus, it
will also not be consistent to use for in (30) any one generalization of the
Langmuir equation, taking into account interactions between the adsorbed mole-
cules. The same inconsistency, of course, must be observed in all the papers on
temperature programmed desorption based on applying (31) when is
represented by the Langmuir isotherm or its generalizations (Bragg-Williams,
quasi-chemical, etc.).
Therefore we would the results of our foregoing analysis of the A R T
approach as follows:
1. The A R T approach is not related to a particular mechanistic model of adsorp-
tion systems.
2. The fact that at equilibrium the A R T approach yields the isotherm is
not sufficient proof that this isotherm should be used to represent in the
generalizations of A R T for the case of heterogeneous solid surfaces. Other
the Rate of
Potential of Adsorbed Mole
At the beginning of the 1980s a new family of approaches to
kinetics appeared. A common fundamental feature of these approaches is that they
relate the rate of kinetics to the chemical potential of the
adsorbed molecules, Although all these approaches appeared at approximately
the same time. the one proposed by Nagai in 1985 could, for historical
reasons, be considered first. This arises because the idea behind that approach
can be traced back to de Boer in 1956. In this approach the rate of desorption
is represented by
Another member of the new family of approaches using the chemical potential
of adsorbed molecules, is the approach proposed by Kreuzer and Payne in 1988
They followed the assumption by in 1938, that
the rate of adsorption should be proportional to the difference between the actual
pressure and the equilibrium pressure corresponding to the
same surface coverage. On this basis they wrote the expressions for and in
the form
and
where is the partition function of gas molecules due to the internal degrees of
A t equilibrium, so takes the form
where
If we assume that it is proper to accept the low-temperature limit for the vibra-
tional partition function then will be an expression in which in the exponent
is increased by a constant, just as it is in (68).
We wrote that both for the case when the sticking coefficient is to be
constant and when it is assumed to be equal to Kreuzer's approach leads
to the isotherm at equilibrium, when = Ward and coworkers
present a different idea that applies the statistical rate theory to the
kinetics of adsorption.
When a molecule coming from a gas phase adsorbs on a solid surface, it must
both cross a gas-phase boundary and form an adsorptive bond when it encounters
the solid. Similarly, when a molecule desorbs. it must break the adsorptive bond
and enter the gas phase. Although the collision of a molecule with a solid surface
can be treated by using either classical o r quantum mechanics. the prediction of the
rate of bond formation o r breaking requires a quantum-mechanical description
because this process involves the probability of an event occurring.
Based on quantum mechanics and Ward and
developed a new approach to kinetics called the statistical
rate theory of interfacial transport. The objective of this approach is to obtain
an explicit dependence of the rate expressions on coverage so that the other proper-
ties appearing in the rate expressions are molecular properties. The S R T approach
leads to the rate expressions
and
where is a constant.
The physical meaning of is discussed in the paper by Ward and
as well as in the more recent publications of Ward and coworkers
is the exchange rate between the gas phase and the solid surface once an isolated
system has reached equilibrium. All these approaches relating the kinetics of
processes to the chemical potential of adsorbed molecules
share a n important feature. Because the expression for the chemical potential is
determined from a mechanistic model of the adsorbed phase, one arrives at a
Let
exp
The above equation describes the rate of adsorption on sites having adsorption
energy W e also believe that the variation of from site to site is a factor
governing the kinetics of adsorption on the energetically heterogeneous surfaces.
I n two previous publications. we showed how the statistical rate theory of inter-
facial transport can be generalized for the case of heterogeneous solid surfaces.
Below we repeat briefly the principles of that method before we develop its general-
ization to take into account desorption kinetics.
The experimentally monitored surface coverage and the mean rate of adsorp-
tion are the values of 0 and averaged over all kinds of adsorption sites
taken with a n appropriate statistical weight. In the presence of a dense spectrum of
adsorption energies, that statistical weight is well described by a continuous dis-
tribution function relating the number of adsorption sites to corresponding values
of adsorption energy,
When the adsorption takes place a t quasi-equilibrium conditions, the function
in Eq. (82) can be written in the form
because the chemical potential of the adsorbed molecules will be the same through-
out the surface phase and will be a function of the average surface coverage The
explicit of this function will, of course, depend on the form of the adsorption
energy distribution
Thus, in the case of a heterogeneous solid surface, the experimentally determined
values and are given by the integral
where
= x
Equation (90) is fully correct only in the limit 0, which represents the
accuracy corresponding to the condensation approximation. This means that it
can be used a t finite temperatures only when the solid surface is strongly hetero-
the variance of is larger than the variance of the derivative
Equation (90) is an equation from which other equations can be developed to
describe the kinetics of adsorption on surfaces characterized by various adsorption
energy distributions.
The in turn is related to the surface coverage
Thus. (94) looks something the Elovich equation and was developed here
as the kinetic equation corresponding to a Temkin isotherm for adsorption equili-
bria. This agrees with common experimental observations. W e see that the Elovich
equation should apply a t conditions where the desorption rate can be neglected,
since we developed it by generalizing the adsorption rate term in (77) and
neglecting the desorption rate term in (78).
As a matter of fact. most of the experiments o n adsorption kinetics were done a t
lower coverages. In the literature we often find statements such as. "After the first
very fast kinetics, a slower kinetics is observed, which follows the Elovich equa-
tion." The first fast kinetics, a t very low surface coverages. is probably governed by
purely kinetic (collision) factors. the effectiveness of the transfer of the energy
excess to the solid phase as well as possible changes in the structure of the adsorbed
phase. The last effect is due to the saturation of chemical bonds. which may be most
essential for the first molecules adsorbed on a strongly "unsaturated," "empty"
solid surface. After a certain first portion of molecules are adsorbed. the rate of
adsorption falls to the limits it runs as a t quasi-equilibrium conditions.
In addition to the above-mentioned common deviations from the Elovich equa-
tion at very low surface coverages, the other common negative deviation was
observed a t higher coverages. T o account for that negative deviation, Higuchi et
al. proposed the following modification of (I):
Other possible reasons for the deviations at higher surface coverages have also been
put forward. F r o m our derivation of the Elovich equation it follows that the main
reason may be the neglect of the desorption kinetics inherent the form of the
equation. Moreover, there is the role of the lateral interactions between the
adsorbed molecules.
Now let us see how inclusion of the lateral interactions affects the form of the
Temkin isotherm for adsorption equilibria and the Elovich equation for adsorption
kinetics.
Using and we arrive at the isotherm For patchwise
topography, we have
The problem is, in fact, similar to that in adsorption equilibria For both
and random topography, criterion (102) leads to the condition
== =-
7
Thus, from Eqs. (100) and (103). it follows that for patchwise topography,
= exp
where
Now we solve the above quadratic equation with respect to This has only one
solution:
and (108) is consistent with the Langmuirian model for adsorption equilibria.
Equation (88) now takes the form
where
= x and = x (121)
The first term on the right-hand side of (120) represents the adsorption rate
and is identical with the expression on the right-hand side of (90). The second
term on right-hand side of (120) represents the rate of the simultaneous
desorption. At quasi-equilibrium conditions, = and can be expressed
as a function of surface coverage. Equation (120) is a generalized form of
which allows us to calculate the rate of adsorption (desorption) if the adsorption
energy distribution is known. The latter function can be determined from the
observed equilibrium adsorption isotherms. For the rectangular adsorption energy
distribution replacing in (120) by defined in (93) leads us to
the following generalized form of the equation, taking into account both
adsorption and desorption occurring simultaneously:
where
and
.9The result of computer fitting of the experimental data for the kinetics of
oxygen adsorption on tungsten reported Lopez-Sancho and de
Segovia, by our generalized Elovich equation. (132). assuming that
a = 0.394, a = 3.85 = 4.374 The pressure and
temperature kept at the experiment were equal 5 ton" and 77 respectively.
and we get
where
2
exp [
exp [ + -
Now let us remark that both the differential equation, and its integral
form, (1 correspond to the level of condensation approximation. Therefore,
their application to the equation
obtained for equilibrium from (108) and the use of the rectangular adsorption
energy distribution (10) should yield the Temkin isotherm. One can easily prove
that this is true, by putting = 0 into (122). The total equilibrium iso-
therm T) then takes the form of the Temkin isotherm
to be with (26).
The condition 1 in (131) and (132) arises from the application of the
condensation approximation and should not appear when, instead of (1
some exact solutions are applied. It can be easily checked that this condition is
equivalent to the condition that T) defined in (134) must be smaller than
unity.
-4 -2 0 2 4
log t
I. Introduction 392
I I . Adsorbed Solution Theories 393
A. The ideal adsorbed solution theory 393
B. The nonideal behavior of the adsorbed phase 396
HI. Equilibria Models Based on Potential Theory 398
A. Direct application of D R and D A equations 400
B. Potential theory approach 403
C. Other methods 405
IV. Other Models 407
A. The extended Sips mode! 407
B. Isosteric heat as a function of loading model 408
C. The heterogeneous extended Langmuir model 409
D . The multisite Langmuir model 41 1
E. The multispace adsorption model 41 1
F. Martinez's model 412
G. Keller's model 413
V. Microporous Structure and Adsorption Equilibria on Activated
Carbon 414
A. The structure of activated carbon 414
B. Potential energy and slit pore configuration 415
C. Characterizing the PSD of activated carbon 417
D . Prediction of adsorption equilibria 420
VI. Microporous Structure and Adsorption Equilibria-the Micropore
Size Distribution Model 42 1
A. MPSD model for a pure component system 422
B. MPSD model for a multicomponent system 424
for each component and for the mixture are the same,
Thus, for a system consisting of components wih specified partial pressure in the
gas phase, there are +
1 independent variables to be solved: P and
Equations (1)-(4) provide the +
1 necessary equations. Since the equations
are nonlinear and the integrals in (2) cannot be solved analytically for most
of the pure component isotherm equations, the + 1 equations have to be solved
numerically.
The IAST has the advantage that no restriction exists for the type of pure
component isotherm equation used for the prediction of multicomponent equi-
libria. The success of the IAST in predicting the multicomponent adsorption
equilibria depends on three factors: (1) the quality of pure component equilibria
data; (2) the correct choice of isotherm equation for the pure component. which
should fit the experimental data well over the entire pressure range; and (3) the
regression algorithm in deriving the parameters for pure component isotherm equa-
tions. Factors 2 and 3 were studied by Richter et Recently, and
Farooq investigated the performance of a number of popular isotherm models
in fitting pure component isotherms with various in a very compre-
hensive way. In the same article, they also examined the performance of each model
equation and compared it with the IAST in predicting the multicomponent adsorp-
tion equilibrium on activated carbon (experimental data on the ternary system
methane-ethane-propane were used).
1. The Procedure
To avoid the numerical integration procedure and iterative scheme in the IAST,
and Myers introduced an algorithm that allows for fast
of multicomponent adsorption equilibria. The pure component isotherm equation
was chosen as
Equations (6) can be analytically integrated for the spreading pressure The
equilibrium amount adsorbed of each species can be obtained from the solution
of a set of A' nonlinear algebraic equations This algorithm (called
was improved by Moon and and later refined by and Myers
It is about one to several orders of magnitude faster than the IAST, depending on
the nature of the problem,
where is the fractional coverage of the surface. The spreading pressure of (2)
can then be rewritten as
= finite value
where is the fluid phase concentration and is the adsorbed phase density
The IAST is based on the assumption that the adsorbed phase is ideal. This is valid
at low pressure, and the IAST is to give good prediction of multicomponent
adsorption equilibria at surface coverages up to about 50% At high pressure,
the model prediction tends to diverge from the experimental data (normally
predicting the experimental data), presumably due to the negative bias from
law. This is expected, because at high pressure the adsorbed phase
becomes more like a liquid, and liquid mixtures of dissimilar molecules are gen-
erally nonideal (nonideality). Sometimes problems also arise regarding the spread-
ing pressure when IAST is applied to binary consisting of strongly and
weakly adsorbing components; when the partial pressure of the 'strongly
adsorbing component approaches its saturation value, the spreading pressure of
the mixture can well exceed the saturation pressure of the weakly adsorbing com-
ponent.
The surface energetic heterogeneity further complicates this problem because
micropores of different sizes or sites with different adsorption energies will present
different adsorption affinities and selectivities toward different adsorbates.
Interactions among adsorbed molecules of the same type and among those of
different types also make this problem more complicated demon-
strated that either of the following two factors can cause erroneous predictions of
IAST even in cases where there are no lateral interactions among adsorbed mole-
cules: (1) a large difference in molecular size or (2) adsorbent heterogeneity. Under
certain circumstances, these two factors could result in an azeotrope that cannot be
accounted for by IAST, Dunne and Myers studied molecules of unequal size
adsorbing into a cavity by using a computer simulation technique,
and they found that the micropore size exclusion effect for larger gen-
erates nonideal behavior in the system.
1. Real Adsorbed Solution Theory
There are two approaches to dealing with these nonideal behaviors in the adsorp-
tion process. One accounts for the nonideality of the adsorbed phase with the
introduction of an activity coefficient; this is called the real adsorbed
solution theory, or RAST In RAST, Eq. (1) is replaced by
where
It can be seen that IAST predictions seriously deviate from the experimental data
for this system even though the adsorbed phase is known to be very close to ideal
under the experimental conditions. This deviation therefore points to the surface
heterogeneity. T o address this heterogeneity, new models continue to emerge in the
study of multicomponent adsorption equilibria.
IAST
Exp. data
Pressure (bar)
FIG. 1 Plot of selectivity of ethane to methane as a function of pressure.
Ref. 47.)
where
A RT
is the adsorption potential, V is the volume adsorbed at the relative pressure
is the saturation pressure at temperature is the affinity coefficient
of the adsorbate, is the micropore volume, and is the characterisic energy of
the adsorbent.
The plot of amount adsorbed versus A is defined as the characteristic curve. The
relationship between and pore size for adsorption in micropores was studied by
Chen and Yang The exponential term which lies in the range of
accounts for the structure of the adsorbent. Large values are found to be asso-
ciated with adsorbents that have a narrow micropore size distribution, while small n
values relate to adsorbents with a wide range of sizes and possibly with
and macropores as well.
The DA isotherm contains parameters related separately to the solid and the
adsorbates. The advantage of using the potential theory for predicting gaseous
mixture adsorption equilibria is that the pure component characteristic curves are
generally independent of temperature. Thus having obtained isotherm infor-
mation of a pure component at one temperature, the mixture equilibria at
other temperatures and pressures can be predicted. The DA equation and its
1. Bering Model
Early efforts to extend the D R equation were largely empirical. Bering et al.
extended the D R equation to binary system in a direct manner. For binary systems,
their model is
where and are volumes adsorbed for each species, and is the saturated
vapor pressure for the adsorbed mixture over which the saturated vapor has the
composition as that of the equilibrium vapor over the adsorbent. The para-
meter is the average affinity coefficient, taken as
= A-I +
T o calculate and the following Lewis relationship for a binary system
is used:
where
Equations (25) and (26) are sufficient for predicting the equilibria of binary mix-
tures. The simplicity of this method constitutes an advantage over other
based on TVFM. Besides, this model is not iterative and can be extended to
component systems in a straightforward manner.
For practical purposes, the volumetric amount adsorbed ( V ) in (25) should
be converted into molar amount adsorbed, with the equation
where
molar volume of the saturated liquid at its normal boiling point
molar volume of the species at its critical temperature
This model is validated with the binary experimental data on several systems
(1) alkane and mixtures on activated carbon; (2) mixture on
zeolite; (3) binary adsorption on BPL activated carbon. The model was compared
with three other models (IAST, G-M, and Bering et model), and it was found
that under most this gives equivalent or even superior perfor-
mance.
where = is the hypothetical pressure for pure component i that gives the
same spreading pressure as that of the mixture. The adsorbed phase concentration
for a binary system is obtained by simply solving (assuming an ideal adsorbed
phase)
The number of moles adsorbed at the hypothetical pressure can be calculated from
the equation as
Theory
The basic assumption in the potential theory approach to predicting the adsorption
equilibria of a gas mixture is that the adsorption data for all pure components on
the same adsorbent will fall on the same line when plotted on some appropriate
coordinate (coalescing) This principle is reflected in the Lewis method
and in the Grant-Manes (G-M) method and its modification by and
Danner
where and are the fugacities of components 1 and 2, and and are the
fugacities of each component at its saturation vapor pressure and at system tem-
perature T. Moon and showed that in most cases the bulk fugacity can be
replaced by pressure. Equation (33) can be used to derive the adsorbed phase
composition by using the concept of pseudoadsorbate pressure and its related
saturated liquid molar volume in an iterative manner.
According to Grant and Manes the affinity coefficient that coalesces the
potential curves of each species satisfies the relationship
where is the saturated liquid molar volume of pure component i at its normal
boiling point.
where and are the adsorption energy and vibration frequency, respectively. If
the adsorption energies of the components are independent of each other, the
following approximation is made:
V=
[ -
where is the relative pressure and is the ratio of the molar liquid
. With Eq. the phase compositions of the jth subvolume can be
calculated directly. The overall adsorbed quantity for a component is the sum of its
related quantities in each subvolume, provided the adsorption potential of that
component is larger than the limiting potential. It should be noted that this
model can describe the real adsorption behaviors when the exponents are
different.
In the case that are different owing to the steric effects, real adsoption
behaviors are to be expected. The multiphase real adsorbed solution (MRAS)
approach is used, considering the "real behavior" to comprise two factors: (1)
adsorbate-adsorbent interactions, which are taken into account by the multiphase
concept; (2) adsorbate-adsorbate interactions, which are considered by activity
coefficients and are taken as the same for the same species in each subvolume.
The MRAS model is
If we used the same analogy as that for the generalization of the Langmuir equa-
tion, we would obtain for the extended Sips model the equation
However, this direct extension corresponds only to the ideal situation where the
adsorption energies of the various components are completely
demonstrated that if the adsorption energies of each compo-
nent are highly correlated to each other by the relation
then the generalized Sips equation for a multicomponent system will take the form
This equation can also be derived by formally applying the IAST, using the Sips
equation, Eq. to describe pure component data. Hence, for adsorption of gas
mixtures onto activated carbon, the use of Eq. (50) is generally justified.
Because the exponential term a # 1 , the Sips equation does not have a correct
low pressure limit, which is contrary to experimental observations. T o overcome
this problem and extend the range of applicability of isotherms of this type, Staudt
et introduced the pressure- and temperature-dependent exponents to ensure
that Henry's law holds when pressure is sufficiently low. The new generalized form
of (48) is
Staudt et suggested the following functional form of a,, which conforms to the
above limits:
where is the critical temperature and k is a constant. With these corrections, the
new model predictions were shown to be much for both single- and
multicomponent adsorption equilibria. This of Staudt et is also applic-
able to other isotherms. For for the Toth isotherm it is
- +
- RT
where is the surface coverage, t is the parameter specific to the
adsorbate pair, (3 is a parameter characterizing the solid property, E(0) is the
isosteric heat at zero loading, a is the extent of energy heterogeneity, is a
function describing the pattern of heterogeneity, and is the adsorbate-adsorbate
interaction energy. Equation (53) has also been extended to study multicomponent
adsorption equilibria. The general form of the extended equation is
The form for the pattern of heterogeneity for the solid is assumed to be the same for
all components, and this pattern is assumed to take the form
elsewhere
The cumulative energy matching scheme is employed to address the energy
matching between different adsorbates in the adsorbed phase. With the uniform
energy distribution, it is
Equation (59) must be solved numerically. The analytical solution is available only
under certain conditions (the variance of the energy distribution is the for all
components).
The HEL model is a noniterative model and has a reasonable prediction cap-
ability. With the pure component isotherm parameters derived the
equation, the model can be employed to predict the adsorption equilibria of gas
mixtures at any other temperature. Since the L U D equation is flexible in correlating
the pure component data, the HEL model prediction for a gas mixture is generally
good. Kapoor et demonstrated in their study that the HEL model gave
satisfactory results for a number of systems, which were comparable in cases
to the results from the IAST or even HIAST
However, arbitrary assumptions such as the energy distribution and
energy matching do exist in the HEL model. Meanwhile, the performance of the
HEL model depends on the L U D parameters derived pure component data.
When the quality of pure component data is not good or there are too few data
points, it is possible that the energy parameters (especially and derived
from the L U D equation are meaningless or inconsistent between different species.
In this case, the model prediction the H E L model must be treated with care.
MSL models present the correct Henry's law limit and obey the requirement of
thermodynamic consistency. This model is often used to account for the effect of
the size difference in the study of multicomponent adsorption equilibria
The MSL model can be extended to include effects such as lateral interactions in the
adsorbed phase as well as surface heterogeneity
It was demonstrated in the study of Gusev et al. that the MSAM gives a
good representation of both ideal and nonideal behavior of the adsorbed phase for
a number of systems that cannot be described by the IAST. For example. the
ideal adsorption behavior of the ethane-methane system on BPL carbon (Fig. 1 ) is
successfully simulated by the MSAM (dotted line in Fig. 1). It is
that the MSAM does introduce some arbitrary assumptions and unrealistic pictures
in describing the adsorbed phase; as the authors put it, "This model is intended for
engineering application rather than for the fundamental mechanisms of
ponent adsorption equilibria."
A similar method that treats the adsorbed phase as a combination of two phases
or patches was used by Li and Talu and G u o Li and Talu show
that with the combination of HIAST and the multicomponent van der Waals
model the highly unusual behavior of benzene-p-xylene vapor mix-
tures adsorbing onto silicalite can be more than qualitatively predicted. The
ideality of the system results in an S-shaped isobaric binary and a
crossover point of selectivity curves a t different pressures. Their results further
support the view that surface heterogeneity is the major contribution to the
'observed'"
. Martinez's
Martinez and proposes a general isotherm model for the adsorp-
tion of both pure gas and gas mixtures. This model takes into account a broad
range of factors such as adsorbate size loss of symmetry or chemical dissociation
clustering 111, and molecular interactions in the adsorbed phase The model is
based on the combination of the kinetics arguments of which account
for the molecular size effect, and statistical thermodynamics arguments which
account for other factors such as clustering and molecular interaction. This model
has two (1) the complete and (2) the reduced form in the Flory limit.
Only the latter form was investigated in detail by Martinez and
it
where is the fugacity and is related to the adsorption energy through the
equation
= = constant
The exponent a depends on the molecular size r and the fractal dimension of the
adsorbent:
= +
and in the case of BET multilayer adsorption, it is
= - + Fx -
where z is the distance between the carbon atom and the adsorbate and
and are the cross LJ parameters, which are calculated using a geometric
mean and arithmetic mean (Lorentz-Berthelot rule), respectively,
1. The Potential
where
is the minimum interaction energy between the adsorbate molecule and a single
lattice layer, is the distance between the adsorbate and the central plane of the
pore, and is the number density of carbon centers per unit are of graphite surface
(38.2
where
where
1. Horvath-Kawazoe Method
The Horvath-Kawazoe (H-K) equation was developed by and Kawazoe
to characterize the microporous structure of an adsorbent. They used the 10-4
potential to describe the potential energy of adsorbate molecules in a slit pore and
where
=
H = pore width
= number
density per unit area
= dispersion constants
From the amount adsorbed at Eq. (74) yields the pore width thus the pore
volume at this pore size can be obtained.
The implicit assumption behind the H-K equation is that the adsorbed phase is a
two-dimensional ideal gas, Henry's law region adsorption. Cheng and Yang
modified the H-K equation by including the nonlinearity of the isotherm. The
modified H-K equation is
where
=2 - 2.78 is the adsorbate diameter
= 0.35 the distance between layers of activated carbon
c,, and can be calculated from known thermodynamics parameters
H is the pore width
Following their assumption, is to by
For a given set of adsorption equilibrium data, (81) can be used (inverted) to
derive the MPSD if a local isotherm is specified. But such an inverse has been
known to be an ill-posed problem, and presents some difficulties in application
For example, for a given adsorbate and a given local adsorption isotherm,
one would find many different distribution functions that can yield the same overall
adsorption capacity (within experimental error). Furthermore, if the isotherm equa-
tion, is applied individually for each adsorbate, we might obtain very
different for the same activated carbon! This violates the basic premise
that the MPSD is an intrinsic property of the given carbon, independent of what
adsorbate is used. Realizing this inconsistency problem often seen in the literature,
we propose a different approach here, applying Eq. (81) to multiple temperature
data of many adsorbates simultaneously while constraining the MPSD as an intrin-
sic characteristic of the adsorbent. T o achieve this, the functional dependence of
interaction energy on size has to be defined first.
This equation is derived by assuming the collision of gaseous molecules toward the
pore mouth of and such a collision rate is known from the kinetic
theory of gases to be inversely proportional to the square root of the molecular
weight and to the square root of The affinity in (85) is
specific only to the solid The theoretical value of reported as
4.26 x
The maximum adsorption capacity of species, takes temperature-
dependent form
where
The parameters and are the lower and upper limits, respectively, for the
inicropore size range. The lower limit is for a pore that has the same adsorption
potential energy as that of a single lattice layer or for a pore that has zero potential
energy. The of pores having sizes two values of is very
and hence value of will suffice. The upper is off at
for convenience, since the adsorption potential energy in pores with
widths larger than this is effectively zero,
The equations have distinct of parameters to be optimized: (1)
the characterizing carbon structure ( 2 ) the isotherm
By optimizing data
of many adsorbates on tlie same activated carbon, it is possible to derive
the results for size distribution as as relevant adsorption para-
meters for
3 4 5 7
Pore half-width
Plots of energy of interaction of ethane and propane versus pore
a
0
0 20 40 60 0 40 60 80
Pressure Pressure
. Plots of MPSD fitting and experimental data of ethane and propane on
activated carbon. (From
Ref 121
3 4 5 6 7
Pore half-width
.4 Plots of the derived micropore size distribution for activated carbon and
Nuxit activated carbon. (From Ref. 126.)
With the structural parameter and isotherm parameters obtained, we employ the
extended MPSD model to simulate the binary adsorption equilibria. As a criterion
to judge the model performance, the related simulation results from IAST are also
presented for the related systems,
1. Ajax Activated Carbon
The binary adsorption equilibria on Ajax activated carbon are measured on a
rig for the following systems (1) methane-ethane (2)
methane-propane (3) (273 and 303 K). For each binary
system, the total pressure for the bulk phase was kept at about 500 Figure
6 shows a typical plot of the amount adsorbed versus the mole fraction of methane
for the system methane-propane-Ajax activated carbon. The continuous lines are
predictions from the MPSD model, and the dotted lines are from the IAST. As the
figure shows, the predictions from the MPSD model are better than those from the
IAST. This trend is also observed for other binary systems studied.
2. Nuxit Activated Carbon
For Nuxit activated carbon, the binary equilibria of the following systems were
simulated: (1) methane--ethylene (2) ethylene-propylene (3)
ethane-propane (293 and 333 K); (4) (293 K); (5) ethane-methane
(293 K); (6) Like the case of Ajax activated carbon. the
model predictions from the MPSD model are again excellent, suggesting that the
MPSD model is appropriate to study multicomponent equilibria of hydrocarbons
on activated carbon.
of the
The MPSDs
A comparison of the MPSDs derived for the two activated carbons in Fig. 4 reveals
that Nuxit activated carbon has slightly larger average pore size and variance than
Ajax activated carbon does, although the difference is not significant. This obser-
vation reflect the general feature of the activated carbons. On the
methane (1)
propane (2)
FIG, 6 Plots of the amounts of methane and propane adsorbed versus the mole
fraction of for Ajax activated carbon. Total pressure is 500 torr, and the
temperature is 303
other hand, the MPSD derived from a KOH chemically activated nutshell-derived
activated carbon shows considerable difference in MPSD from those of commercial
carbons the physical and the reliability of the derived MPSD
deserve further investigation. The microporous structure of activated carbon is very
complex; the assumption of any form for will impose some constraint in the
inverse process of (81). Furthermore, the assumption of the slit-shaped model
of infinite extent is an oversimplified picture of a real micropore, in
which the other sources of surface heterogeneity, as we stated before, are not
accounted for. Due to the sensitivity of the LJ potential to the pore size, the
MPSD obtained will locate mainly in the region of one to two LJ diameters. So
the MPSD derived here represents the effect of micropores of all sizes on
the overall adsorption equilibria. Considering the complexity of the carbon struc-
and the two-parameter gamma distribution function, the MPSD obtained this
way should be treated qualitatively rather than quantitatively, in spite of its sig-
nificant role in describing both pure and binary adsorption equilibria on hetero-
geneous surfaces.
In its theoretical aspect, the MPSD model also provides us with some insight
into the role of surface heterogeneity. It takes into account the size exclusion effect
and sets a physical criterion for the local matching energies between different
adsorbate molecules in the micropore. Although the local extended iso-
therm in the model is simplistic in its description of true equilibria, this problem can
be overcome by applying the concept of IAST in each and the overall
adsorption equilibrium is the contribution from the local equilibria over the com-
plete M P S D range. However. the model will become iterative and very
while its performance improves only marginally. The significance of this
sacrifice for purely theoretical purposes is questionable.
However, as we discussed before, the present model represents only a primitive,
picture of equilibria in activated carbon. All other
interaction mechanisms are also lumped into the potential minimum of the adsor-
bate-pore interaction, as the model does not account for interactions with the same
adsorbate molecule, nor does it reflect the importance of surface loading as the pore
filling progresses. These suggest that more fundamental work is needed for a full
understanding of the adsorption equilibria on heterogeneous
solids.
In recent years, the behavior of fluid molecules in small pores has been well
studied by computer simulation Computer simulation can provide us
with valuable information about the microscopic behavior of adsorbate mole-
cules confined in small pores in terms of intermolecular and surface forces,
thus enabling us to understand the fundamental behaviors of adsorbates in the
potential field of the pores. The current trend in the literature suggests
that for physical adsorption in activated carbon, the adsorbate-adsorbate inter-
action and adsorbate- pore interaction are well represented by the LJ potential
theory while the model micropore is a slit-shaped channel of infinite extent.
This forms the basis for the appliction of the statistical method in adsorption
processes.
which is determined by the ratio of Henry's law constants for each species and
hence independent of the bulk fluid D F T theory, this is expressed as
Cracknell and Nicholson first studied the selectivity of the same methane-ethane
system using the G C M C technique in which the ethane molecule is treated as a
spherical fluid They also investigated the predictability of the IAST theory
for the binary mixture in pores of different sizes. Their study was carried out at the
supercritical temperature (T = 296.2 K), and the results were compared with those
of Tan and Gubbins. It was found that in most cases the D F T overestimates the
selectivities (of ethane), but there is qualitative agreement between the results from
the two methods. The following reasons are attributed for the discrepancies: (1) The
weighting function for density expansion does not have enough terms; (2) the mean
field assumption for the attractive force is inadequate.
Cracknell and Nicholson further investigated the same binary system using a
more realistic model for ethane, in which the ethane was treated as two sites
connected by a certain length of bond The simulation results showed that
the elongated shape of ethane strongly influences the pore size dependence of
selectivity and that the confinement in may cause a significant
decrease in the selectivity. Their studies were also extended to the ethane-propane
system in which the ethane and propane were considered as two- and
site LJ fluids, respectively. Again, using the simulated single-component data as
input, the IAST was found to be able to predict the selectivity with reasonable
accuracy.
et al. also studied the selectivity of the ethane-methane system with
the G C M C method in model graphite micropores of various sizes at supercritical
temperature and under a broad range of pressure. Both adsorbates were treated as
I. Introduction
I I . Methods of Characterizing Porous Silicas on the
Basis of Nitrogen Adsorption Data
A. Introductory remarks
B. Evaluation of the specific surface area and pore volume
C. Calculation of the pore size and pore size distribution
D. Estimation of the adsorption energy distribution
III. Studies of Porous Structures and Surface Properties of
Organic-Modified Conventional and Ordered Mesoporous Silicas
A. Nitrogen adsorption properties of unmodified conventional
and ordered mesoporous silicas
B. Physically coated mesoporous silicas
C. Conventional mesoporous silicas with bonded organosilanes
D. Ordered mesoporous silicas with bonded organosilanes
IV. Conclusions and Future Research Perspectives
References
Porous silicas are widely used in areas of modern science and technology as
adsorbents, catalyst supports, additives, and chromatographic packings, to men-
tion just a few Many of the applications result from their high surface area,
tailored particle size and pore size, and high mechanical, thermal, and chemical
stability Another advantage of porous silicas is their high surface reactivity
toward coupling agents, which makes it possible to obtain a broad range of
materials with desired surface properties that are useful for the separation of
. Introductory
Nitrogen adsorption at 77 is a standard and widely used method to determine the
specific surface area and pore size distribution of adsorbents
This technique is also applied to determine microporosity
However, it is rarely appreciated that analysis of low pressure nitrogen adsorption
data affords considerable insight into the surface properties of porous materials
Adsorption instrumentation capable of performing low pressure measure-
ments (down to relative pressures of about is currently commercially avail-
able. Therefore, it is convenient to carry out nitrogen adsorption runs over the
entire pressure range in order to obtain information about surface area, porosity,
and surface properties from a single adsorption isotherm. It is well known that
porous usually exhibit strong structural and surface heterogeneities, the
former caused by the presence of pores of different sizes and shapes, and the latter
resulting from surface irregularity, roughness, and the presence of impurities and
various functional groups. Adsorption processes on solid surfaces are influenced by
both structural and surface heterogeneities, and it is often difficult to recognize
which of these factors is responsible for adsorption properties observed experimen-
tally.
Our recent studies of nitrogen adsorption on mesoporous silicas
were aimed at exploring the possibility of separating pore struc-
ture and surface functionality effects in order to develop better characterization
procedures for these materials. The low pressure adsorption properties of porous
silicas were found to be markedly different from those of many other adsorbents,
such as graphitized and carbons and porous and zirconia or
even alumina However, low pressure nitrogen adsorption on numerous
conventional and ordered mesoporous silicas proceeds in essentially the
way, indicating that nitrogen probe molecules are quite insensitive to details of
silica surface and structure The similarity of the surface
valuation of the
The Brunauer-Emmett-Teller (BET) method provides a standard procedure to
calculate the specific surface area of mesoporous materials However, if an
adsorbent has an appreciable number of mesopores with sizes close to the
pore size range (2-3 nm), adsorption data from the region of capillary condensation
in these mesopores should not be used in the BET calculations In such cases,
the accuracy of surface area determination can be improved by using the compara-
tive plot method as described below. The total pore volume can be estimated
from the amount adsorbed at a relative pressure close to the saturation vapor
pressure (usually about 0.99) by converting the adsorbed amount to the corre-
sponding volume of liquid adsorbate The conversion factor is equal to
0.0015468 for nitrogen at 77 K when the adsorbed amount and the pore volume
are expressed in cubic centimeters at STP per gram and cubic centimeters per
gram, respectively, and it is assumed that the density of condensed adsorbate is
equal to the density of liquid adsorbate at the temperature of measurement. This
procedure provides the volume of pores with widths below 100-200 depending
on the pressure value used in the calculations. In the case of adsorption isotherms,
which exhibit a steep increase in the adsorbed amount at pressures close to the
saturation vapor pressure, the estimation of the total pore volume may be quite
inaccurate.
In many cases, information about structural properties of high surface area
materials can be obtained by comparing the adsorption isotherm of a solid under
study with the adsorption isotherm of a reference macroporous material of similar
surface properties The comparison is made by plotting the amount
adsorbed for the solid under study as a function of the amount adsorbed for the
reference material at the same pressure. The adsorption for the reference is
usually expressed either as the statistical film thickness of the adsorbate on the
surface (in the t-plot method) or as the standard reduced adsorption (in the
plot method) The quantity a, is defined as the ratio of the adsorbed amount at
a given pressure to the adsorbed amount at a relative pressure of 0.4.
In the case when the solid under study and the reference material have the same
surface and structural properties and differ only in the amount of surface area, the
comparative plot is linear and starts from the origin of the graph. However, usually
the plots obtained on the basis of experimental data exhibit deviations from linear-
ity, which can be analyzed to obtain information about the structural and surface
properties of a given solid. If the solid under study has they are filled at
low relative pressures, and the amount adsorbed in them does not change signifi-
cantly with further pressure increase. The increase in the gas uptake for higher
pressures results mostly adsorption on the external surface area of the sample,
Many methods to calculate the pore size are currently available (see Refs. 74, 75,
and 84 and references therein), so it is important to choose a proper tool for the
analysis of mesoporous silicas. Our studies showed that the
Halenda (BJH) method is a very useful characterization technique In
the BJH approach, the adsorption process in mesopores was considered to consist
of the formation of the adsorbed layer on pore walls and subsequent condensation
of the adsorbate in the interior of the pore. The desorption process was assumed to
where is the molar volume of the liquid adsorbate, 7 is its surface tension, R is
the universal gas constant, and T is the absolute temperature. In the case of nitro-
gen adsorption at = 8.88 and =
In the BJH method, calculations are performed starting from the
desorption isotherm data for the highest pressures attained during the experiment,
where it is that all pores are completely filled with the condensed adsor-
bate. When the pressure is subsequently lowered, a certain number of the largest
pores are by the capillary evaporation, but the walls of these pores are still
covered by the adsorbate film. The decrease in the amount adsorbed, which was
caused by the lowering of pressure, can be recalculated to obtain the volume of
these emptied pores. but in order to do that, assumptions need to be made about
pore geometry and density of condensed adsorbate. The latter can be approximated
by the density of liquid adsorbate at the saturation vapor pressure and temperature
of the As far as the geometry is concerned, Barrett
assumed that pores are cylindrical in order to make calculations more tractable.
Hence, the volume of the inner space in the pores that were emptied during the first
evaporation step is equal to the volume of liquid adsorbate corresponding to the
measured decrease in the adsorbed amount. T o arrive at the total volume of those
pores, one simply needs to account for the thickness of adsorbate on the pore walls.
The pressure for which the considered decrease in adsorption took place can now
be related to the pore radius through the Kelvin equation [Eq.
A reasoning similar to that described above for desorption from the largest
pores during the first pressure-lowering step can be applied to subsequent decreases
in the adsorbed amount. But in these cases, the decrease in adsorption is caused not
only by the capillary evaporation but also by thinning of the adsorbed layer in
pores in which the capillary evaporation has already taken place. The calculations
are repeated until a sufficiently low pressure value is reached or the cumulative pore
volume obtained using the procedure described above equals or exceeds the total
pore volume. The BJH procedure follows the sequence of events during desorption
measurements; it actually does not make any difference whether adsorption or
desorption data are used. However, different relations between the pressure and
the pore size may need to be used for capillary condensation and evaporation
It needs to be noted that the Kelvin equation [Eq. was shown to underestimate
It is well known that the adsorption of a gas on a solid surface depends the
interactions between the surface and gas molecules These interactions
can be described in terms of the adsorption energy distribution (AED) for a given
adsorbate-adsorbent The AED can be analyzed to obtain information
about the surface, about its structure and functional groups
In order to calculate one needs two fundamental pieces of information.
First, adsorption isotherms for homogeneous surfaces interacting with the adsor-
bate with specified adsorption energies need to be known Second, one needs
to account for the topography of adsorption sites random, o r inter-
mediate between these two extremes) The topography is expected to
influence the adsorption on a heterogeneous surface unless lateral interactions
(interactions between adsorbed molecules) are negligible. T o deal with these pro-
blems, one needs to make suitable approximations. Under certain assumptions
concerning the properties of the adsorbent and the adsorbed film (geometry of
the surface, number of layers of adsorbed molecules, lateral interactions, mobility
of adsorbed molecules, etc.), one can derive analytical expressions for adsorption
isotherms on homogeneous surfaces For instance, the Fowler-Guggenheiin
(FG) equation was derived for localized adsorption with lateral
where the overall adsorption and were divided by the monolayer capa-
city of the adsorbent in order to obtain the overall degree of surface coverage
= and the fraction of adsorption sites of energy
The monolayer capacity can be calculated using the standard BET
T o calculate the population of adsorption sites as a function of adsorption
energies normalized AED), the integral equation for the overall adsorption
(5) or needs to be inverted. The required overall adsorption isotherm is
directly measured experimentally, but the form of the local adsorption isotherm
must be specified. The F G adsorption isotherm was used in our calculations, and
the reasons for its selection are discussed elsewhere One of the useful
features of the F G equation is that it can easily be modified to account for differ-
ences in lateral interactions in the cases of different topographies of adsorption
sites. For patchwise topography, the standard expression for the adsorption
isotherm can be used:
CBPB-B22
LiChrospher Si-100
CBPB-T24
CBPB-B22
10
Pore Size (run)
FIG. 4 BJH pore size distributions calculated using desorption isotherms for the
unmodified silica gel LiChrospher Si-100 and samples physically coated
and chemically bonded (CBPB-B22) with CBPB. (Data from Ref. 77.)
LiChrospher Si-10
CBPB-T24
Adsorption Energy
FIG. 6 Comparison of adsorption energy distributions for unmodified
100 silica gel and physically coated (CBPB-T24) as well as chemically bonded (CBPB-
samples with similar loadings of CBPB. (Data taken from Ref. 77.)
MCM-41
36% CBPB MCM-41
CBPB-BT33
0 4 8 12 16 20
Adsorption Energy
7 of adsorption energy distributions for an unmodified MCM-41
material, MCM-41 sample coated with CBPB, and a silica gel sample that was
chemically bonded and subsequently physically coated with CBPB. (Data from Refs.
46 and 77.)
Si-100
CBPB-BT33
Relative Pressure
Results of the fitting of low pressure adsorption data for a physically coated
silica gel sample with a relatively low loading of CBPB using low pressure
data for the unmodified silica gel and the modified sample (CBPB-BT33) with large
amounts of CBPB chemically bonded and physically coated. Data for CBPB-BT33 are
assumed to represent adsorption on pure CBPB. (Data from Refs. 46 and 77.)
Adsorption Energy
9 Adsorption energy distributions for silica gel samples with bonded
octyldimethylsilyl ligands. (Data from Ref.
to nitrogen molecules. Consequently, the lack of high energy sites on the AEDs for
the octyl phases under study does not necessarily indicate the removal of these sites
from the surface, as they could be effectively shielded by the long-chain
groups.
In another study low pressure adsorption properties of polymeric
and bonded phases were examined. It was shown that the introduction of
relatively small polymeric 3-aminopropylsilyl groups led to a decrease in the num-
ber of high adsorption energy sites in the range of A slight further
decrease in the number of these sites was observed when the amine ligands were
converted to with phenyl, cyanophenyl, and octyl groups,
especially for the latter. The similarity of AEDs for and
bonded phases confirmed previous findings that it is impossible to
detect the presence of groups attached to aromatic rings of the ligands,
which indicates that these two groups exhibit essentially the interaction ener-
gies with respect to nitrogen molecules. However, it was possible to distinguish
between adsorption properties of and bonded phases,
which is also consistent with the previous studies discussed in and
It was shown earlier that polar parts of organic modifiers gave rise to adsorption
energy sites with energies up to about 12 However, the bonded
phases exhibited adsorption sites with energies up to about 15 One can
expect that the latter sites are groups present in the structure of the bonded
phase due to hydrolysis of side alkoxy groups and/or on the silica surface. This idea
was supported by the examination of the polymeric octyl phase which exhib-
ited a much smaller population of adsorption sites between 10 and 13 but
showed a pronounced peak at about 13 Since octyl chains interact with
Although it is possible to study pore size changes resulting from the modification of
conventional mesoporous silicas, as was shown in such changes are
much easier to analyze in the case of ordered mesoporous silicas due to their
narrow pore size distribution and relatively small sizes. Moveover,
there has recently been a significant and rapidly growing interest in studies of
organic-modified ordered materials It was therefore
interesting to investigate the surface modification of a good quality large-pore
MCM-41 material via chemical bonding of trimethylsilyl (TMS), butyldimethylsilyl
(BDMS), and octyldimethylsilyl (ODMS) monomeric ligands and polymeric 3-ami-
nopropyl (APS-P). hexanoyl-3-aminopropyl (HAPS-P), and octyl phases
The obtained surface coverages were in the range of 2.5-3 It was
shown that the increase in the size of bonded ligands led to a systematic shift of the
step of capillary condensation in primary toward low pressures (Fig.
10). The overall shape of the adsorption isotherm remained unchanged, showing
that the modification procedure did not change the structure of the ordered
porous support. Pore size distributions for unmodified MCM-41 and
selected modified materials (Fig. 11) indicated that the modification procedures
3 4 5 6 7 8 9 1 0
. Pore Size
FIG. 11 Pore size distributions calculated from adsorption isotherms using the BJH
method with (3) and (4) for an unmodified large-pore MCM-41 and
bonded materials with polymeric 3-aminopropylsilyl (APS-P) and monomeric
(ODMS) ligands. (Data from Ref. 59.)
used markedly decreased the pore size of the material but did not lead to an
appreciable broadening of the PSD. It was shown that the increase in the size of
bonded ligands resulted in a gradual decrease in the pore size. For example, the
initial pore size of support decreased to 4.5, 4.3, and after
modification with BDMS, and ODMS ligands, respectively. Moreover,
there was good agreement between the pore size changes in the case of different
ligands of similar sizes. For instance, the pore diameters of the materials with
monomeric and polymeric octyl phases were 3.8 and respectively.
sizes of and aminopropyl phases of similar organic chain length but different
functionality were also very similar (4.3 and 4.4 respectively). It was concluded
that when reproducible surface coverages of bonded ligands can be achieved, it is
relatively easy to predict the pore size of the modified materials on the of the
pore diameter of the support and the size of the ligand. Such predictions can be
based on the data, which were presented and discussed in detail elsewhere
Adsorption Energy
The results presented here demonstrate the great usefulness of nitrogen adsorption
measurements for the characterization of unmodified and modified mesoporous
silicas. High pressure adsorption data allow for a detailed characterization of
mesoporous structures of materials, providing their surface areas and pore size
distributions. Low pressure measurements provide an opportunity to study the
surface properties of materials and to monitor changes in these properties that
result from surface modifications. For instance, it is possible to attribute the
appearance of adsorption sites of certain energies in the adsorption energy distri-
bution to the presence of some types of surface groups as silanols. aliphatic
chains, or or polar segments of organic ligands. Low pressure data can
also be used to estimate the surface coverage of organic species on the silica surface.
Therefore. nitrogen adsorption measurements carried out over a wide range of
pressures it possible to characterize both surface and structural properties
of silicas.
However, further studies are needed to improve our understanding of the rela-
tion between the adsorption properties of modified mesoporous silicas and their
structure and surface functionality. First, additional studies of modified materials
with different functional groups would be desirable to gain additional insight about
the types of surface groups that can be detected and distinguished from one another
by nitrogen probe Second, other nonpolar or weakly polar probe mole-
cules need to be examined in order to assess their applicability and sensitivity to
details of the surface, and it is very likely that some of will exhibit properties
very useful for surface characterization. Third, it is important to improve methods
of pore size analysis. New approaches based on studies of model adsorbents with
well-defined structures MCM-41) appear to be promising and successful, but
they require further refinement to improve their accuracy for materials with differ-
ent pore geometries and surface properties. Finally, the currently proposed analysis
of the effects of the pore structure and surface functionality on adsorption proper-
ties of materials is relatively easy due to the possibility of separating
these effects. Namely, high pressure adsorption data are influenced mostly by
properties of the porous structure, whereas low pressure data reflect mainly the
I. Introduction
II. Solid Nanospace of Enhanced Molecular Field
III. Magnetic Interaction of Molecular Spins
IV. Magnetic Properties of Bulk and Two-Dimensional 0 2
V. Pore Structures and Magnetic Properties of Activated
Carbon Fibers
VI. Adsorption of in a Graphitic Nanospace
VII. Magnetic Susceptibility of Confined in a Graphitic Nanospace
VIII. Dependence of Magnetic Susceptibility of Confined on Pore
Width
IX. Effect of Coexistent on Magnetic Susceptibility of
Confined
X. Magnetic Distinction of Weakly Interacting Molecules in
Nanospace
References
II.
interaction between a molecule and a surface atom as a function of
the distance between them can be expressed by the Lennard-Jones potential,
where and are the well depth and effective for the
graphitic carbon atom complex. These cross parameters are calculated according to
the Lorentz-Berthelot rules, = + Here, and
are the Lennard-Jones parameters for a surface atom and a molecule,
respectively. The interaction potential for a molecule and a single graphite
slab is given by the Steele 10-4-3 potential
where is the vertical distance of the molecule above the surface, A is the separa-
tion between graphite layers 0.335 nm), and is the number density of carbon
atoms in a graphite layer (= 1 As the of activated carbon can
be approximated by the slit spaces between the predominant basal planes of
graphitic units, the whole interaction potential of a molecule with the
micropore of an intergraphite surface distance can be given by the equation
Consequently, we can evaluate the potential profile of the molecule adsorbed in the
graphitic Here is not the effective pore width determined by the
adsorption experiment. The difference between and is a function of and
We can estimate the effective pressure from the potential profile; molecules
confined in a slit pore of 1 nm width are presumed to be exposed to the high
pressure of 100 Therefore, the graphitic micropore can offer a high pressure
field from the macroscopic view
Here is the magnetic susceptibility. When the magnetic field is weak enough, is
independent of molecular magnetism the cgs emu system has been used
instead of the SI The volume magnetic susceptibility is a dimensionless
quantity traditionally expressed in the dimension of the electromagnetic
unit (emu) is formally the cubic centimeter. The specific magnetic susceptibility
which is obtained by dividing the volumetric magnetic susceptibility by
the density, is often used. is the of the diamagnetic ( y o ) and paramagnetic
susceptibilities.
+
The value is negative and is positive. When dominates, the sample is said
to be diamagnetic; it is repelled by the magnetic field. When is the leading
contribution, the sample is paramagnetic: it is attracted by the applied field.
where T and are the temperature and Boltzmann constant, respectively. This
combines macroscopic and microscopic quantities. However, the application of the
general equation (9) is difficult without the knowledge on = Van
introduced the approximation
where is the energy level in zero magnetic field. and are called the
first- and second-order Zeeman coefficients respectively. When
exp
T
If the first excited states are to be too high in energy to couple with the
ground state. g is in principle isotropic and equal to As is linear in the
above Van formula with = 0 can be used.
which leads to
is rewritten as
-
This is the Curie-Weiss law. is the Curie constant, and is the Weiss
constant defined by = +
A plot of for a system obeying the Curie-Weiss law gives a straight
line with a slope o f c ' . The intercept with the T axis yields both the sign and the
value of 6. In the framework of this model, a positive indicates ferromagnetic
interactions and negative 6 indicates antiferromagnetic interactions. Hence the
vs. T plot can determine the sort of magnetic interaction. In the case of ferro-
magnetic interactions the Curie-Weiss is limited to the temperature
range T Accordingly, if we use a paramagnetic probe such as an
the low temperature magnetic susceptibility measurement can provide an
essential information on the intermolecular states of confined using this
analysis.
The ground state of an molecule is the triplet and therefore gas shows
the representative paramagnetism. The effective spin concentration and con-
stant of bulk gas are 1.5 and 0 K . gas liquefies at 90.2 and
the liquid solidifies at The solid oxygen has three
K T and T The
magnetic properties of solid oxygen have been actively studied since the days of
Kammerlingh Onnes. Solid oxygen has two sharp magnetic transitions cor-
responding to a-(3 and crystallographic changes. Active studies [35 showed
that the a phase is antiferromagnetic, with two-dimensional ordering. It is
suggested that the (3 phase has a "frustrated" nature
It was indicated that the and liquid phases have a common feature
Although is a simple elemental solid. there are still many obscure pro-
blems with respect to the properties and structures of
Oxygen gas of high purity (99.99%) was sealed in a highly pure quartz tube. Its
magnetic susceptibility was measured with a SQUID magnetometer system over a
temperature range of under a magnetic field of 1 Figure 2 shows the
temperature dependence of the magnetic susceptibility of condensed The
value in Fig. 2 was determined just after the equilibration of the measuring tem-
perature. The relations on cooling and heating between 2 and almost
coincide with each other. There are two anomalies at 23.8 and 43.8 which are
ascribed to the well-known and transitions
Temp
spins by structural destruction. The insensitive behavior of the phase may come
from the disordered nature of the structure Jezowski et reported that
doping of or Ar in had a drastic effect on the thermal conductivity, but the
phase was insensitive to the doping, which has a similarity to the effect of coex-
istent
The magnetic interactions of adsorbed on the graphite surface have been
extensively studied The mono- or of molecules adsorbed on the
two-dimensional graphite surface has phases such as 5, and that are different
from the bulk a , and phases; the phase change depends on both temperature
and coverage. The and phases appear in the multilayer of adsorbed on the
flat graphite surface. Also the to phase transition leads to a magnetic anomaly.
A of the two-dimensional magnetic structure of on the graphite basal
plane is shown in Fig. 4. This diagram was obtained at nonequilibrium conditions.
In the following discussion, the magnetic properties of confined in nanospace
obtained at nonequilibrium conditons are used for comparison with published
results.
nonporous carbon black were used for the standard carbon. The subtracting pore
effect (SPE) method for the high resolution a,, plot can remove the enhanced
adsorption, giving the correct area The high resolution a,, plot
for P-1.45 has two swings linearity in the regions of 0-0.3 and
The line between the two swings passes through the origin, and its slope provides
the total surface area, that is, the sum of the and external surface
areas. The slope of the line of = 1 gives the external surface area. and its inter-
cept leads to the micropore volume. Hence the high resolution plot provides the
total surface area, the microporous surface area, and the micropore volume. The
micropore geometry of ACF can be approximated by the slit, and thus the average
slit width (micropore width was determined from the microporous surface area
and the micropore volume. The a, plots of P-0.75 and P-0.96 had the only swing
below 0.3. Their surface area and volume were determined in a
similar way. Table 1 lists the total surface area, micropore volume, and average
micropore width for the three ACFs.
SUSCEPTIBILITY OF
P-0.75 is also observed a t 0.3 1.2, but the relaxation is much faster than that
of P-0.96. Therefore, P-0.75 offers a restricted molecular environment less suitable
for cluster formation than P-0.96.
The curve of adsorbed in the micropore of P-1.45 a t = 0.021
(d nm) satisfies the Curie-Weiss law 2.6 emu and
2.3 K) as is shown in Fig. 1 l a . The curves for 0.42, however, give different
features from those of the other two samples described above. T h e curves for
= 0.42, 0.60. and 0.78 have broad peaks a t 28, 36, and respectively.
Moreover, the bulk phase appears a t 1.1 (d 0.38 n m and = 0.63). This
is a smaller coverage than with either 0.96 o r P-0.75. The changes in and
are somewhat similar to a curve for the 0 phase in the system obtained
by Kobler and Marx These authors propose that liquid and phases coexist in
the region of 1 0, 1.8 and 32 T 40 (32 K is the melting point of the 8 phase;
40 that of the 0 phase) for the system. The phase is presumed to
have a liquid crystalline structure in which molecular axes are ordered like those of
0.02
lad 0.015
gas cluster
30
,
gas + cluster
a gas
low temperature range. It is plausible that molecules confined in the slit pore
of 1 n m width cannot form a two-dimensional lattice similar to that on the flat
graphite surface. The relationship between the disappearance of the 2 D lattice
formation and the slit width should shed light o n the unusual state of
confined in a narrow space.
Water molecules can form a network structure that accepts a guest molecule to
produce a clathrate compound under application of high pressure o r low tempera-
ture Recent in situ X-ray diffraction studies showed that molecules
confined in the graphitic of A C F form an ice-like structure and
produce a clathrate compound different from the bulk one even near ambient
conditions Such a clathrate structure can isolate an spin system.
Thus, knowledge of the magnetic susceptibility of the mixture confined
in the of A C F a t different compositions and fractional fillings over the
temperature range of is necessary for understanding assembly.
The results for P-0.96 are shown here. The diamagnetic contribution of ice
(-0.699 x emu g) was subtracted from the observed value The coexistent
affects seriously the magnetic structure of even the bulk condensed as
mentioned earlier. Earlier studies showed that the y-T curve of confined
for 0.16 0.25) has no critical change in the temperature range of 2-100
and confined shows within a low range; an almost
free spin state was confirmed by the Curie-Weiss parameters close to those of
gas. Temperature independence of was observed in temperature range of
20-50 K for 0.16 which is associated with magnetism due
to spin cluster formation having an explicit relaxation. The temperature depen-
dence of of confined in the of A C F is shown in Fig. 15 as a function
of in the presence and absence of Here. the content is expressed in
terms of the molar ration R,,. Although the coexistent brings about a
paramagnetic nature as mentioned above, it does not sensitively affect the magnetic
behavior of confined for because of the free spin state even without
the coexistent O n the other hand, the addition of excess noticeably
changes the relation, which indicates the perfect isolation of by coexistent
Figures 16 and 16b show the Curie-Weiss plots for R,, and 28.7,
respectively. The coexistent bends the plot upward markedly, indicating that
the coexistent accelerates the spin cluster formation. inducing random mag-
netism. In the case of the Curie-Weiss plot bends down at for both
cases, with and without T h e spin concentration is decreased slightly by the
presence of while the constant does not change, as shown in Table 3.
The decrease of the spin concentration by the presence of supports the for-
mation of spin clusters. Association of molecules should lead to a n island
formation of molecules, producing the spin cluster. The Curie-Weiss plot in the
absence of is linear. whereas that with a n excess of is convex. The
constants from the Curie-Weiss plot are also shown in Table 3. The spin
0 10 20 30 40 50 5 10 15 20 25
Temp K Temp K
Chemisorbed species solid surfaces have been widely studied. and established
methods are available for determination of the different chemisorbed species. F o r
example, five representative species of oxygen- 0;.
and 0 - h a v eb een evidenced with the aid of I R and ESR for oxygen
o n the surface of transition metal oxides However, interactions in physical
adsorption, even micropore filling. are too weak, and we have no established
method for assigning the different physisorbed states. The interaction potential
profiles of a n molecule with slit graphite are shown in Fig. 1. Although the
potential minimum difference of three shown in Fig. 1 is only 150 K, which is
very small, it is large enough for magnetic interaction. There is a possibility that a
slight difference in a weakly interacting system can be measured by the magnetic
method. The magnetic method is quite sensitive, and we identify different oxygen
states in micropores with the aid of magnetic susceptibility measurements at low
temperature, as described above. The magnetic behaviors of micropores can
provide a n adsorbed mechanism. The average intermolecular distance d is calcu-
lated from the fractional filling using the assumption of isotropic and homogeneous
I. Introduction 50 1
II. Heats of Adsorption by Gibbsian Thermodynamic Model 502
III. Heat Removal from and Supply to an Isothermal Adsorption System 506
IV. GSE Model for Mass and Energy Balances in Adiabatic
Adsorption Columns 507
V. Indirect Estimation of Isosteric Heat of Adsorption 508
A. Pure gas isosteric heat 508
B. Multicomponent gas isosteric heat 510
C. Shortcut method 510
VI. Direct Estimation of Isosteric Heat of Adsorption 513
VII. Characteristics of Pure Gas Isosteric Heat of Adsorption 516
A. Homogeneous adsorbent 518
B. Heterogeneous adsorbent 518
C. Limiting value of isosteric heat of adsorption 520
D . Special cases of adsorbent heterogeneity 520
E. Simple models for isosteric heat of adsorption 52 1
VIII. Characteristics of Multicomponent Gas Isosteric Heats of Adsorption 521
A. Limiting values of isosteric heat of adsorption 524
B. Simple models of multicomponent gas isosteric heats of adsorption 524
C. Example of model heterogeneity for a binary system 526
References 527
tinuously along the length of the adsorber during various steps of the and TSA
processes The local gas adsorption equilibria and kinetics within the adsorber
depend on the local gas pressure (P), the gas and adsorbent temperature and
the gas-phase composition (mole fraction of component Thus, it is critical
that accurate knowledge of the heats of adsorption and desorption of the adsor-
bates as a function of and be available so that correct temperature changes
within the adsorber can be estimated during the process steps. A theoretical sensi-
tivity analysis of the performance of a simple PSA process (Skarstrom cycle) for
separation of a binary gas mixture using a detailed mathematical process model
showed that small errors in the estimation of heats of adsorption of the adsorbates
and their variation with adsorbate loadings could significantly alter the process
performance Furthermore, the study clearly demonstrated the pitfalls of ignor-
ing adsorber
Traditionally. a variety of heats of adsorption and desorption for pure and
gas-solid systems have been defined by using thermodynamic
models Experimental techniques have also been developed to measure
these heats These generally use the actual amounts adsorbed as the
primary variables for representing the extents of adsorption of the adsorbates.
Unfortunately, the Gibbsian surface excesses (GSE), and not the actual amounts
adsorbed, are the only true experimental variables for measuring the extent of
adsorption In view of this fact, a detailed thermodynamic model for
component gas adsorption equilibria using GSE as base variables has already been
developed
This chapter (1) elaborates on the thermodynamic properties of multicomponent
heats of adsorption and desorption using the GSE framework, (2) demonstrates the
practical utility of these heats by developing mathematical models to describe mass
and energy balances in adsorption columns using GSE, (3) describes direct and
indirect experimental methods for estimating the heats, (4) provides experimental
and theoretical examples of the heat of adsorption of pure gas and binary gas
mixtures on microporous adsorbents (energetically homogeneous and hetero-
geneous), and ( 5 ) comments on the complex nature of heats of
adsorption. Only solids are considered in this chapter because of their
practical values.
Figure 1 shows the GSE model for equilibrium adsorption from a bulk
gas of i components N) characterized by and
The adsorption system contains a unit amount of an inert adsorbent whose total
helium void volume is It is assumed that all pores of the adsorbent are
accessible to the nonadsorbing helium gas. The dotted line in Fig. 1 represents the
interface separating the Gibbsian adsorbed phase and the bulk gas phase. It
is arbitrarily located inside the actual bulk gas phase. The adsorbed phase has a
.1
Schematic of Gibbsian surface excess (GSE) model for gas
adsorption.
where and can be the internal energy the molar enthalpy the
molar entropy the molar Gibbs free energy (g). and the molar volume of the
adsorbed and bulk gas phases, respectively, is the corresponding excess function.
It should be noted that is equal to zero by this definition, and thus
(= + is equal to
For pure gas adsorption, Eq. (3) can be used to calculate T) by measuring
the pure gas surface excess as a function of gas-phase pressure at constant
T (adsorption isotherm for pure gas). For a multicomponent gas mixture, Eq. (3)
can be used to calculate by measuring the surface excess of component
as a function of total gas-phase pressure (P) at constant T and The
variable of (3) in this case is the total surface excess for all components
R is the gas constant.
The total internal energy and the total enthalpy of the entire
adsorption of Fig. 1 at and are given by
It can be shown from Eq. that for a closed system (constant with an ideal
gas phase, one has
For the special case where is not a function of and Eq. (18) simplifies to
The negative sign for the quantity in Eq. (20) represents the case where heat is
removed from the system in order to maintain system isothermality. Equations
and (20) can be combined to obtain
Equation (21) was previously derived by many authors using the actual amount
adsorbed instead of the GSE as the base variable. The equation demonstrates
the practical value of the differential heats defined by Eq. (9).
Adsorbent Particles
Equation (25) represents the instantaneous column energy balance using the
framework. The key thermodynamic variables in the equation include q, and
defined by Eqs. (9) and respectively. The last term on the right-hand side of
Eq. (25) vanishes for isobaric process steps (constant P).
The quantity in (25) represents the instantaneous rate of change of
surface excess of component i at Z and t. This replaces the local rate of adsorption
term in the conventional energy balance equation derived by using the actual
amount adsorbed as the base variable. A simple model to describe the rate of
change of GSE is the analog of the conventional linear driving force model
Equation (27) shows that the isosteric heat of adsorption of pure gas i can be
estimated at a given value of P and T or T and by (1) measuring a set of pure gas
adsorption isotherms a function of P at constant T ) a t different tempera-
tures, (2) generating a set of adsorption isobars as a function of a t constant
P) from the isotherms or measuring them independently, and (3) estimating the
isotherm and isobar slopes at the desired values of P and The use of Eq. (27) for
calculating for pure gas is, however, very sensitive to errors in measuring the
isotherm and isobar slopes. On the other hand, it is very useful to develop analytical
expressions for as functions of T and the adsorption isotherms can be
described by analytical models
It can be seen that the experimental protocol to generate data for use in (29)
and (30) can be extremely time-consuming and impractical. We are not aware of
any such data in the published literature. The degree of complexity increases enor-
mously when the gas phase has three of more components.
It is also unfortunate that Eq. (1 1) cannot be easily integrated for
nent adsorption as in the case of pure gas adsorption because no multi-
component adsorption experiment can be designed to generate as functions of T
at constant = 1. 2, . . . , N). Consequently, there is no easy way to estimate
multicomponent (i 2) isosteric heats of adsorption by measuring as functions
of P. and
Equation (33) shows that (1) the difference between and q, increases as P or
increases (isotherm slope decreases) and (2) a, because R T ) is generally
positive.
For a binary system one gets
It must be very clear from the above discussions that the available indirect methods
for estimation of may not be very useful for 2) adsorption
systems. A calorimetric method is necessary where the heat evolved (or consumed)
due to specific changes in the adsorption system is directly measured. Recently, a
heat flux microcalorimeter was developed to measure the heat of
adsorption Figure 4a is a schematic diagram of the calorimeter. It consists
of a thermostated cell whose walls are surrounded by thermopiles for directly
measuring the heat flow signals through the cell walls. Figure 4b shows the typical
profile of thermopile response (in microvolts) versus time during an experiment.
The area under the curve of Fig. 4b gives the total energy removed from the
calorimeter during an experiment. It is necessary to calibrate the response curve
using a known heat generation system within the cell.
The cell is connected with instruments for measuring the gas-phase pressure and
composition. It is also connected to a thermostated gas injection system (dosage
side) through a valve. A unit amount of the adsorbent is placed in the cell (sample
side), and it is initially saturated with pure gas i at and or a gas at
and The experiment consists of very slowly introducing a very small quan-
tity of pure gas i (AN,, into the sample side from the dosage side through
the valve. The temperature of the gas entering the sample side is monitored
as a function of time in conjunction with the thermopile response profile. The
dosage side pressure and temperature are also as functions of time.
The system temperature is at both the start and end of the experiment.
The amount of pure component introduced during the experiment is given by
where is the rate of introduction of pure gas into the side. It can
be calculated by knowing the volume of the dosage side and the rate of
change of gas density in the dosage side during the experiment. t*
is the time interval of the experiment. The amount of energy (enthalpy,
introduced into the sample side of the calorimeter during the experiment
is given by
and are, respectively, the changes in the surface excess and the
phase partial pressure of component in the sample side. is the helium void
volume of the sample side.
An energy balance for the side of the calorimeter (open
no shaft yields [1
where is the total change in internal energy of the sample side of the
calorimeter during the experiment. is the total of heat removed
from the sample side of the through the during the experi-
ment. Equations (1 and (42)-(45) be combined to get
at constant T
where is the local surface excess isotherm of pure gas on a homogeneous patch
of sites characterized by the isosteric heat of adsorption is the normalized
probability density function for the distribution of isosteric heats of adsorption:
The lower and upper limits of in the energy distribution function are generally
assumed to be zero and infinity for mathematical convenience.
Equations (48) and (49) have been used (1) to calculate the pure gas isosteric
heat probability distribution function for a heterogeneous adsorbent from
of
Statistical thermodynamics of adsorption dictate that the pure gas adsorption iso-
therm of component be linear at the limit of 0:
where is the isosteric heat of adsorption of pure gas in the Henry's law region.
is a constant. Equation (52) is derived by integrating Eq. assuming that
is independent of T (generally valid). It follows from (51) that is a constant
in the Henry's law region.
for
Table 3 lists several simple, analytical, adsorption isotherm models and the corre-
sponding isosteric heats for pure gas adsorption. These models have been fairly
successful in describing pure gas adsorption on microporous adsorbents (type I
isotherms
The homogeneous Langmuir (MSL) model assumes that a,
adsorption sites are occupied by each molecule of adsorbate The saturation sur-
face excess for the pure adsorbate i is given by which is not a function of
T. is the fractional coverage for pure component is the
gas-solid interaction parameter for component b, is a function of T only (Table
3), and is a constant. The isosteric heat of adsorption of pure component i
for this model is independent of coverage. The MSL model reduces to the original
homogeneous Langmuir model [4] when = 1 .
The MSL model has been extended to incorporate lateral interactions in the
adsorbed phase In this case, the adsorbent is energetically homogeneous, but
increases linearly with increasing due to lateral interaction. is the pure gas
lateral interaction energy parameter for component This model reduces to the
original homogeneous Fowler-Guggenheim model [7] when a, 1. Furthermore, it
reduces to the Langmuir equation when a, = 1 and = 0.
The Toth model is derived for pure gas adsorption on a heterogeneous adsor-
bent without lateral interactions The parameters and b, are identical to
those for the Langmuir model. The parameter which is a function of T only,
is the heterogeneity parameter. The adsorbent is heterogeneous if 0. The
adsorbent heterogeneity increases as the value of the parameter k decreases. The
model reduces to the homogeneous Langmuir model when = 1. The Toth model
can be derived using (48) and (49) by assuming that the adsorbent is patchwise
homogeneous, consisting of a collection of Langmuir sites with a skewed
like distribution of site isosteric heats of adsorption The for the Toth model
decreases as increase. An inconsistency of the Toth model is, however, that
as 1.
ill.
There are very few published experimental data on multicomponent gas isosteric
heats of adsorption Even those are limited to the study of binary gas
adsorption systems. Figures 8 and 9, respectively, show the isosteric heats of
adsorption of and from their binary mixture on zeolite crystals
measured as functions of loading (surface excess) at approximately constant
0 1 2 3 4 5
= 0 ; constant P
The Henry's law constant in Eq. (53) is equal to the Henry's law constant of
pure gas [Eq. and it is a function of T only. The variable in Eq. (54) is a
different Henry's law constant, and it is a function of both P and It follows from
(1 and (54) that
The three pure gas adsorption models described by Table 3 have been extended to
cover adsorption equilibria. Table 4 gives the multicomponent
compared to Thus the and pure gas heats are comparable only for
component 1 at higher values of . These results show that it may not be correct to
use pure gas isosteric heats of adsorption as functions of its coverage to represent
the isosteric heat of adsorption of the components from a gas mixture at the same
coverage for process design purposes.
I. Introduction
I I . Surface Chemistry of Activated Carbons
III. Adsorption of Inorganic Solutes
A. Adsorption of cations
B. Adsorption of anions
IV. Adsorption of Organic Solutes
A. Role of electrostatic interactions
B. Role of dispersion interactions
C. Toward a comprehensive adsorption model
V, Summary and Conclusions
References
It is well known that the uptakes of gaseous, vapor, or liquid adsorbates by differ-
ent carbon adsorbents can vary over a remarkably wide range. Many of these
variations are due to differences in the physical surface properties, to surface
area effects, as well as to different degrees of molecular sieving because of their
widely varying pore size distributions. However, when such factors are allowed for,
large uptake differences per unit of total (or accessible) surface area still remain
These are especially important in liquid-phase adsorption; they have been well
documented, but they are not well understood. They may also be important in
vapor- and gas-phase applications. This undesirable state of affairs is surprising,
having in mind the fact that of the four technologically relevant adsorbents acti-
vated carbon is by far the most important one, with estimated worldwide sales of
about being an order of magnitude higher than those of its competitors
.1 Correlation (or lack thereof) between the surface area of an activated carbon
and the uptake of (a) water vapor at 298 K and (b) pyridine at 423 K. (
area (from adsorption at 77 K); Dubinin-Radushkevich (DR) surface area (from
adsorption at 298
FIG. 2 Principal features of the activated carbon surface that give rise to surface
charge development in aqueous solution.
figure deserve special attention: (1) the origin of positive charge and (2) the
relative contributions of functional groups vs. graphene layers to the development
of surface charge. In Fig. 3 one such graphene layer is depicted, albeit in two
dimensions only. It consists of saturated basal plane sites, (2) unsaturated edge
sites, (3) edge sites saturated with heteroatom-containing functional groups, and (4)
electrons delocalized over the entire graphene layer. The small size and poor
alignment of the crystallites make them highly porous and reactive. The high
osity, especially in the micropore range (with equivalent diameters of
exposes large surfaces to gaseous and liquid environments. The crystallite faces
are composed of fused benzene rings (graphene layers). The functional groups at
the crystallite edges are ubiquitous because of carbon's great affinity for oxygen;
they possess a wide range of thermal and chemical stability. Most of them are acidic
carboxyl and groups), but some may bc basic carbonyl plus ether,
or groups). Some edges are terminated by thermally stable groups,
while others are much more reactive because are terminated either with free
radicals unpaired electrons) or with somewhat stabilized electron pairs.
In contrast to more ordered carbons heat-treated carbon blacks), where the
overwhelming majority of the surface consists of basal planes, in activated carbons
the ratio of edge to basal plane carbon atoms can be as high as 10-20%.
Three possible sources of the positive charge are highlighted in Fig. 2: basic
oxygen-containing functional groups pyrones or chromenes), protonated
amino groups (for nitrogen-rich carbons), and graphene layers acting as Lewis
bases and forming electron donor-acceptor (EDA) complexes with molecules.
When nitrogen-containing functional groups are not abundant, which is the case in
+ + + OH'
and at pH = =
Taken together, these results strongly suggest that (1) there is an important
contribution of the graphene layers to the positive surface charge of carbon
bents and (2) the widely held belief (especially reflected in the literature reviewed in
that surface charge develops as a consequence of both deprotonation and
protonation of the same oxygen functional groups should be discarded.
Protonation of aromatic carboxyl groups in aqueous solution is highly unlikely.
[Thus, for example, only acids whose strength exceeds - 7, HI
= can protonate aromatic carboxyl groups, while is, of
course, a much stronger acid than and is the
dominant species in water.] Instead, protonation of the basal plane is consid-
ered to be much more important for the applications of interest here.
The use of titrations is commonly considered to be an alternative to electro-
phoresis for the estimation of surface charge. Actually, in the case of solid carbons,
and especially in the case of high surface area activated carbons, this turns out to be
a complementary method and failure to make the distinction can be mislead-
ing. Electrophoretic mobility is considered to be a response to primarily external
surface charge on the particles, while titrations probe both the external and internal
surfaces of particles. For convenience I shall refer to the titration-derived para-
meter (pHpzc) as the point of zero charge, while the electrokinetically derived
parameter will be referred to as the isoelectric point. I acknowledge that
a more fundamental distinction between these two parameters requires a dis-
cussion of specific adsorption (in an ideal situation, the two should be identical in
the absence of specific adsorption), but such a discussion is typically not included in
adsorption studies, especially of organic compounds. In most cases of practical
interest, because of preferential (diffusion-controlled)
air-induced oxidation carboxyl group generation) on the external surfaces of
carbon particles. Hence, the smaller the difference between the two values, the more
homogeneous the distribution of surface charges. This is illustrated in Table 1. A
dramatic example of this difference between and has been shown
(albeit implicitly) by Dixon et for the case of a commercial activated carbon,
10.2 vs. 3.5.
ox
16h
+ + OH-
In addition to giving rise to a positive surface charge, the graphene layers parti-
cipate in TI-TI interactions with aromatic adsorbates. These are known to
be sensitive to substituent effects on the aromatic ring The substituent
effects on the carbon adsorbents are examined in detail in and
Suffice it to note here that there are at least three ways in which the
density of a graphene layer can decrease (1) deliberate of
the edges with oxygen functional groups, by controlled oxidative treatments; (2)
accidental (and more or less inevitable) decoration of the edges with oxygen
functional groups upon adsorbent exposure to ambient air; and (3) additional
localization of electrons at the edges by forming "in-plane pairs" with the
0 20 40 60 80 100 120
Time of atmospheric exposure (days)
"The activated carbons A, P were prepared from almond shells, olive stones, and peach stones,
respectively. Their aqueous suspensions all had pH 7.0.
surface area, total pore volume, and micropore volume (in pores 7.5 were determined
by adsorption and mercury porosimetry.
'Calculated using Langmuir adsorption capacities obtained at 293K and pH 5.
Source: Ref. 67.
when pH was above (even though some of this increase was attributed by
the authors to hydroxide precipitation). The authors calculated for "a
hydrated radius of approximately 25 A which is larger than the reasonable size of
hydrated cadmium ions."
These results indicate that a much smaller fraction of the adsorbent surface is
typically covered by inorganic adsorbates than by organic (especially aromatic)
adsorbates (see IV). A logical inference is that, in contrast to the mechanism
of adsorption of organic solutes (see IV), the graphene layers are not involved
in the adsorption of inorganic cations from dilute aqueous solution. Instead,
cation-exchange processes or interactions with functional groups at the graphene
layer edges are primarily responsible for the adsorption of metals. The study by
et is interesting in this context. These authors found that "saturation
adsorption [of an activated charcoal cloth] increases slightly as the p H increases
from 3 to 5.5," which is consistent with a decrease in electrostatic repulsion or a
higher degree of ion exchange. They assumed that at the saturation adsorption
capacity at pH 5.5, mercury exists as a species (even though Baes and
Mesmer suggested that neutral is predominant) and thus calculated
a surface coverage of 3.5% for the unhydrated cation; however, when a single
sphere of hydration was assumed to surround it, an unusually high coverage of
90% was obtained. One possible reason for such a large uptake is that this
carbon (whose surface chemistry was not described) had a very high surface func-
tional group concentration; perhaps a more likely explanation is that the metal
concentration used was much higher (by several orders of magnitude) than those
typically encountered in wastewater streams.
Clearly, further fundamental studies are needed to confirm these conclusions.
For example, it will be interesting to find out whether much higher surface cov-
erages can be accomplished on a carbon whose maximum number of
exchangeable) adsorption sites, - is not only created but
also made electrostatically accessible by adjusting the solution chemistry (and in the
The adsorption of anions is a much more controversial topic than the adsorption of
species. It is quite intriguing that in many cases the uptake does not appear
to be much higher, as one might expect if the positively charged graphene layers
(see Figs. 2 and 3) were electrostatically accessible and "active" in adsorption.
Indeed, in a review of anion adsorption by Hingston the use of carbonaceous
adsorbents is not even mentioned.
For example, and Hancock cite and discuss a study of
and in which it was "calculated that the surface coverage amounted to less
than one gold atom per 5 n m 2 of available carbon surface." They thus "excluded
basal planes, carboxylic acid groups and basic oxides as so-called 'points of attach-
ment' and proposed that the adsorption sites were probably quinone-type groups"
or "special type of In sharp contrast, Klauber notes that "it is
difficult to convert the pmol loadings into a sensible figure for aurocyanide
coverage" but proceeds to calculate (based on X-ray photoelectron spectroscopic
results) that the "extrapolated loading of 2000 pmol approximately equates to a
0.65 monolayer equivalent"; he thus concludes that "the graphitic planes of the
activated carbons [are] the location onto which the anions adsorb."
Table 3 summarizes the results of Solar et which clearly illustrate the
dramatic improvement in molybdenum adsorption once the electrostatic repulsion
is eliminated. Thus, for example, the surface coverage of the as-received carbon
black increased by more than an order of magnitude at pH despite the
uncertainties regarding the surface area occupied by the adsorbed species and the
complex speciation diagram of these high values could well be inter-
preted to represent a large fraction of the adsorbent surface.
The need to purify both municipal and industrial wastewaters places increasingly
stringent demands on the performance of activated carbon adsorbents. Reliance on
the physical surface properties high surface area, adequate pore size distribu-
tion) is thought to be particularly insufficient for the cost-effective removal of
aromatic pollutants and the subsequent reactivation of the spent adsorbent. A
thorough discussion of the opportunities that judicious design of surface chemistry
offers in these applications is presented elsewhere together with a comprehen-
sive and critical review of the relevant literature. Here I summarize one of the key
issues: the relative importance of electrostatic versus dispersive interactions in the
adsorption of aromatics,
It is well known that methylene blue and other dyes are used to estimate the
adsorption capacity of carbons and even to characterize their porous structure
This implies that nonspecific dispersion adsorbate-adsorbent interactions are
dominant in these systems and that the entire surface is accessible to the adsorbate.
However, as Fig. 7 illustrates, the adsorption capacity depends very much on the
chemistry of the carbon surface and the pH; unless the uptake by this carbon is
determined at a relatively high pH, the measured surface area can be underesti-
mated by a factor of 2 from 500 to This means that electrostatic
interactions can be very important as well.
Based on extensive experimental evidence regarding the importance of and
surface chemistry there is no question that the arguments discussed in the
previous section in the context of adsorption of inorganic solutes are applicable
also-and indeed are required--for understanding the adsorption of organic
solutes, many of which are weak electrolytes. What does need careful consideration
here are the answers to the following two questions:
1. Which specific features of carbon surface chemistry need to be known, in both
qualitative and quantitative terms, in order to be able to optimize an activated
carbon for a specific pollution control task?
2. What are the relative roles of electrostatic and dispersive interactions in the
various adsorbate-adsorbent systems of interest?
As-received
Filtrasorb 400
Oxidized
Filtrasorb 400
H e a t treatment in at for 4 h; oxidation in boiling for 1 h
Source: Ref. 127.
WVB
CG6
Nonspecific dispersion interactions between the adsorbent and the adsorbate are
the dominant driving force for the adsorption of gases and vapors. Therefore, in the
absence of molecular sieving effects, the entire carbon surface is "active'" in the
removal of gaseous or vapor-phase pollutants. It was seen in that adsorp-
tion of inorganic solutes takes place predominantly on specific sites on the carbon
surface, by ion exchange, and that very high concentrations of solute or sur-
face functional groups are necessary for achieving high removal efficiencies. Here I
show how the interplay of electrostatic and dispersion interactions governs the
uptake of aromatic pollutants on activated carbons.
My interest in this topic was sparked by intriguing initial results for the adsorp-
tion of benzoate, and fumarate anions on chemically modified (oxidized vs.
nitrided) activated carbons For example, both oxidation (which reduces the
of the carbon) and nitriding (which raises the were found to sup-
press the uptake of benzoic acid; this is contrary to the expectations based on
electrostatic arguments alone. On the basis of these and other findings, my co-
workers and I proposed the following mechanisms of adsorption
1. Adsorption of benzoate (aromatic) anions occurs primarily on the basal plane
of carbon, and the electron-withdrawing effects of nitrogen and carboxyl
functional groups suppress the interaction of the basal planes with the
adsorbate's aromatic rings.
2. Adsorption of aliphatic anions occurs also predominantly on basal plane,
and the same electron-withdrawing effects enhance its interaction with the
carboxyl anions.
More recently, we studied in more detail the relative importance of dispersion and
electrostatic adsorbate-adsorbent interactions as we examined the uptake of
methylene blue and nitrophenol as well as that of benzoic, oxalic, and fumaric
acids In agreement with our earlier findings, we concluded that while electro-
static interactions are important, dispersion interactions appear to be dominant in
the adsorption of aromatic solutes. On the other hand, electrostatic repulsion
appears to be much more important for the adsorption of aliphatic anions.
Further validation of these conclusions was achieved by contrasting the
behavior of three chemically different carbons in adsorbing two vastly different
aromatic solutes; nitrobenzene is a very weak Lewis acid that possesses the
GCW Americas),
"BET surface area.
'Micropore volume obtained from the equation applied to the
adsorption isotherm.
treatment in inert atmosphere at 1223
'Treatment in at 363 K for 6
with at 873 K for 3 h.
Source: Ref. 138.
Reduced concentration
0.0
0.000 0.001 0.002 0.003 0.004 0.005
Reduced concentration
Reduced concentration
Reduced concentration
In the above isotherm equation, K is the equilibrium constant for the molecules
(due to London dispersion interactions), while K' is the equilibrium constant for
the ions and has the form
exp exp
In the above expression, which consists of two driving forces for adsorption, is
the valence of the charged species, F is the Faraday constant, is the interfacial
electrical potential due to the development of surface charge, and is the potential
arising from the dispersive adsorbent-adsorbate interactions. Implicit in this
expression is the realization that the co-ionic solute species an
cation interacting with the surfaces a t pH can adsorb only if the (attrac-
tive) adsorption potential is greater than the repelling electrostatic potential
Explicit exclusion of co-ionic species from the surface is included in the
model for a heterogeneous adsorbent surface, which conforms to the Freundlich
isotherm:
with
I. Introduction
I I . Characterization of Porous Structure
A. Nitrogen adsorption
B. Mercury porosimetry
C. Limitations of the cylindrical pore model
D . Determination of surface charge density and surface potential
III. Thermodynamic Background of Charge Regulation
IV. Surface ionization
A. Protolytic properties of partly oxidized silicon nitride: Heuristic
ideas
B. Distribution of surface potential and surface charge
over pore sizes
C. Application of charge regulation model for determination of
equilibrium constants of surface reactions
References
When two different phases containing some carriers of electricity, such as electrons
or ions, are brought into contact with each other, charge separation occurs, leading
to the formation of an electrical double layer This is a rather general contact
phenomenon that occurs in a variety of systems, not necessarily involving a solid.
Known examples are emulsions containing an ionogenic surface-active component,
or solution systems. Indeed,
the origin of the charge and the mechanism of charge transfer may be different,
When talking about a "porous solid", it is usually meant that the pore system has a
fixed geometry that is not affected by the surrounding environment. Although more
alleged than proven and not thermodynamic at all, this assumption of absolute
rigidity makes it possible to bring together results of porosimetric measurements
normally performed with a dry solid with potentiometric and electrokinetic data
that bring information about the surface charge and potential of the solid in con-
tact with an aqueous electrolytic solution.
Compared to nonporous solids, charge regulation at the surface of porous solids
reveals some new features if the mean pore size is (1) large enough to justify
consideration of the surface charge, which is discrete by its very nature, as con-
tinuously smeared over the surface, but (2) small enough to ensure sufficient
Gas adsorption is by far the most common technique used to gain information
about the porous structure from the isotherms. As adsor-
bent, nitrogen is the most commonly used. The treatment of results here is based on
a cylindrical pore model.
The adsorption isotherm is measured at 77 gradually increasing the nitrogen
pressure, from zero up to the saturation pressure, Indeed, the boundary
points are excluded because of instrumental limitations. When the relative pressure
exceeds a certain limit, capillary condensation commences. resulting in filling
pores with liquid adsorbate. As continues to increase, larger and larger pores
are filled, until all porous space is flooded at 1. Then, by reversing the course
of events, the desorption isotherm can be obtained.
Assuming that the surface is perfectly "wetted" by adsorbate, the contact
angle is zero, the maximum pore radius, filled at the point of the isotherm,
corresponding to the relative pressure can be found from the Kelvin equa-
tion,
where stands for the mole amount of nitrogen (per unit mass of solid) adsorbed
at relative pressure
Plotting vs. yields the so-called differential pore size distribution, prob-
ably the most widely used characteristic of a porous solid.
In a more rigorous treatment, a correction should be made to account for the
finite thickness of a multilayer molecular film preadsorbed on the pore walls prior
to the capillary condensation. Most modern instruments automatically execute all
these routines.
The same data can be used for an evaluation of the specific surface area. On
noting that the internal surface area of cylindrical pores with radius in the range
to + is
As, =
and substituting expressions (2) and (3) for and respectively, the well-known
Kiselev formula for the specific surface area, follows at once:
Mercury porosimetry makes use of the fact that the pressure, required to intrude
liquid mercury into a cylindrical pore with nonwettable walls is inversely pro-
portional to the radius, of the pore. This is a direct consequence of the
Young-Laplace equation,
and expressing from Eq. (7). Further processing of data follows essentially in
the same ways as in the case of nitrogen capillary condensation described in the
previous section.
Like any idealized model, the cylindrical pore model has many limitations, knowl-
edge of which is essential for proper interpretation of experimental data. In parti-
cular, the following should be kept in mind:
1. Except when cylinders are plugged a t one end, the model fails to explain the
existence of hysteresis between the adsorption and desorption branches of
isotherm. The pore size distributions calculated from these two branches
may differ considerably. A typical example is given in Fig. 1. The desorption
branch always gives a more narrow distribution shifted to the fine-pore
region, which is, in general, in better agreement with that calculated from
mercury intrusion data. At the same time, unless different pores are isolated,
the desorption isotherm may be influenced by some network effects, such as
pore blocking, which do not affect the adsorption branch.
2. The actual pore size may be underestimated at the lower end of the
ous region because of some decrease in surface tension of capillary liquid with
increasing curvature of its
3 In gas adsorption measurements, there is a critical lower value of relative
pressure, 0.4 for nitrogen, below which capillary condensate can-
not exist as a separate phase, no matter whether the pore system extends to
finer pores or not; this is a consequence of the existence of so-called tensile
strength limit As a result of this, a very narrow false maximum may
appear on the differential pore size distribution at its fine-pore end, typically
between 17 and 20 A in the case of nitrogen.
0.8
0.6
0.4
0.2
0.0
Neither the surface charge density nor the surface potential at the
interface can be measured directly. They are to be retrieved from pH potentiometric
o r electrokinetic data on the basis of certain theoretical concepts and models.
unit area of the surface is equal to the of protons (hydroxyl ions) bonded or
released at that area
where is the Avogadro number, V the volume of the solution, the mass of the
dispersed solid, the water dissociation constant, the initial concentration
of protons determined by the amount of added, and the equilibrium
concentration due for measurement. For insoluble solids in the absence of concur-
rent ion adsorption, can be identified as the number of elementary charges, per
unit surface area: that is, can be treated as the surface charge density.
2. Electrokinetic Measurements
In principle, both the potential and charge density at the surface of a porous solid
can be calculated from electrokinetic data such as the electro-osmotic transfer rate
or the conductivity of a porous plug. Interpretation of this type of experimental
data is based on the solution of the so-called electrokinetic equations
and
But again, the obtained in way cannot be directly compared to the charge
density that appears in the charge regulation model.
Notwithstanding the existing limitations of the interpretation of experimental
results, electrokinetic measurements are a valuable tool for locating, paradoxically
in the present context, the uncharged state of ionizable surfaces. As shown later,
no proper interpretation of pH potentiometric charge data is possible
without knowledge of the conditions under which the net surface charge becomes
zero, commonly specified by indicating the point of zero charge or the isoelectric
point.
and between the chemical potentials and the mole fractions, of the components,
Although preferable from a theoretical viewpoint, the use of mole fractions as unis
of concentration may be inconvenient in the case of heterogeneous systems. Since,
however, mass balance is usually implied, concentration units can be changed as
appropriate.
It is expedient to express via the increments,
The same result can be arrived at in another way. Introduce a so-called adsorp-
tion phase, which is composed of adsorbed components and the solid itself, referred
to as the adsorbent. Next, it is necessary to define excess quantities. To d o so, a
reference system is required Normally, the latter is taken to be composed of
the same components, and under the same conditions, as the system in study
but without adsorption. Then the excess quantities, Qf, are defined as
= -
where is some quantity measured in the main system and is the same quantity
measured in the reference system. Chemical potentials of all components except
adsorbent can be chosen to be the in both systems. Then, according to general
where
from which. after comparison with Eq. the same result as expressed in (19)
follows at once.
Now that Eq. (17) is available and well substantiated, it is possible to proceed
with its application to the problem of charge regulation. The underlying idea here is
that for solids, especially those containing ionizable surface groups, brought
into contact with an aqueous electrolytic solution, the main contribution to 7 is
made by the free electrostatic energy, of the interface,
=
where can be calculated according to the formula
where the apparent equilibrium constants, K,, have been introduced. Then, making
use of the evident identity,
where the prime reminds us that the sum should be taken over surface-attached
components only. The latter equation together with the mass conservation principle
constitute the basis of charge regulation.
Before proceeding further, it is of concern to demonstrate the logic of the basic
equation (30) in explaining equilibrium properties of a few real processes. Consider,
for instance, two proteolytic reactions that occur at the surface of amine-modified
silica in contact with aqueous electrolytes:
Example 1. Protonation of surface amino groups,
Here =
3 .
and therefore this time is expected to increase
together with This somewhat trivial, through, for the equilibrium
constant of the reverse reaction is just the reciprocal of that for the direct reaction.
Of course, the above conclusions are so evident that there would not probably
be any need to resort to Eq. (30). They are intended to demonstrate its general
consistency.
A',
SiOH
Such a set of reactions responsible for charging the surface constitutes a charge
regulation model. Writing down the mass action law for each of the reactions,
where the subscript is added to emphasize the pore size dependency of the corre-
sponding quantities. This result can also be obtained by using charging
process, (29).
Now let be the pore size distribution normalized as follows:
determines the average free electrostatic energy of the solid whose specific surface
area is given by
The right-hand side of Eq. (53) is a function of as can be seen after retrieving the
concentrations of protonated amino groups and deprotonated from
( 3 5 ) . By taking into account the evident relations,
and the fact that sinh = and sinh = the existence of a finite
real root of the transcendental equation sinh = f ( u ) is established. The finiteness
of the root necessitates equality of to thus simplifying the task of
finding that root to solving the quadratic equation
where
where the primes over the summations recall as usual that the summation is limited
to surface-attached components.
Figure 6 compares the pH dependencies of surface charge density and surface
potential calculated either neglecting or allowing for the above-mentioned effect.
+ +
F o r a given pore size the solution to this charge regulation problem is found from
the system
where
FIG. 7 Titration data for mesoporous silica gel. The characteristics of the sample are
specified in the caption to Fig. 1 . Ionic strength: 0.001 0.01 and (A)
0.1 Discrete points are the experimental results; solid lines are the best-fit
curves calculated according to the charge regulation model based on reaction set (59)
6
with the following parameters: = 4.6 = 1.05 x
= 1.0 = = 2.371 A. (From Ref. 28.)
where is the pore size distribution defined in (48). The above equation can
be rewritten in a finite-difference representation reflecting the discreteness of experi-
mental data,
Ill
40
Ill
30
0)
20
0)
10
= (A) -
where v is the volume of the solution, is the of the sample, and is the
product of the solubility of orthosilicic acid and its dissociation constant. Now,
and as well as some other parameters of interest, can be used as fitting para-
meters and evaluated by minimizing the difference, usually expressed as the squared
deviation, between the theoretical and experimental results
Knowledge of intrinsic equilibrium constants of surface reactions allows
lation of surface specification as a function of pH and ionic strength of the solution,
sometimes leading to quite unexpected results. As a n example, Fig. 8 demonstrates
the percentage of different species a t the surface of phosphate-doped silica gel
It is the charge regulation that allows protonated silanols to exist in appreciable
amounts u p t o p H as high as 6-7 (for pure silica gel, the amount of protonated
silanols becomes negligible a t a p H of about 4 a t the same time dramatically
reducing the dissociation ability of phosphate groups. Thus, the charge regulation
model proves to be effective in the analysis and modeling of protolytic and ionic
equilibria a t charged interfaces.
I. Introduction
I I . Surface Charging Mechanisms
A. Unequal hydration
B. Surface hydrolysis
C. Adsorption
D. Isomorphous replacement
I I I . A Hybrid Triple-Layer Model
IV. Site-Binding Model
V. Balance Equations
VI. Determination of Site-Binding Constants
VII. Prediction of Surface Charge Density
VIII. Effect of Ion Adsorption on the Zeta Potential
IX. Prediction of Surface Species
X. Current Status
XI. Summary
References
Nanosized particles and porous materials play an important role in many techno-
logical applications. They have been used widely as biological sensors, integrated
electronic devices, catalyst supports, and adsorbents for detoxification of industrial
effluents and domestic water, decoloring of processing water, and purification of
pharmaceuticals and proteins. In Chapter 3, Xu showed some typical applica-
tions of functionalized superparamagnetic particles in biological cell separation and
industrial effluent detoxification. Porous media. on the other hand, constitute an
..
tration of metal ions present in addition to the nature of the surface itself As a
result, the interaction of and porous media with targets such as metal
ions is anticipated to be affected by changing either surface properties or solution
conditions. It is therefore important to understand the charging mechanism and to
be able to predict the variations of surface charge (the sign and magnitude) in
response to changes in solution.
In this chapter, a brief summary is presented of various charging mechanisms,
each with a typical example. The progress on the development of a surface
plexation model to describe quantitatively the adsorption of or anionic
species is then reviewed, with emphasis on contributions from our research labora-
tory. A general approach to modeling the charged particles and predicting surface
species distribution is described. Finally, a few examples are given to illustrate the
applications of the approach to an understanding of molecular adsorp-
tion at a mineral surface.
Hydration
When an ionic solid crystal, such as is brought into contact with pure water, it
often becomes negatively charged although the solid particle is electrically neutral.
The negative surface charge indicates the existence of excess surface anions (I-).
The asymmetrical distribution of surface cations and anions may be attributed to
unequal hydration between cations and anions. A stronger hydration of a cation
results in a lower activity of that cation in the solution phase and hence a
surface deficiency. The surface charge can in this case be predicted by comparing
the hydration free energy of cations and anions. Provided that the concentration of
is increased to above the surface positively charged
because of the increased activity of in the solution and hence excess on
the surface. I n this example, and I are referred to as potential-determining
ions (pdi's). Most inorganic solids suspended in an aqueous solution in which they
have some solubility are charged according to this mechanism. Other examples
include and
It is well known [3] that natural silica particles are negatively charged in pure water,
which can be attributed to surface hydrolysis followed by the pH-dependent dis-
sociation of surface hydroxyl groups. When a silica matrix is broken in air, the
surface concentration of and are the same to meet the requirement of
electrical neutrality. When the fractured surface is brought into contact with water,
these charged sites, called unsaturated sites, are readily hydrolyzed to form neutral
surface species of releasing either O H or Adjusting the pH, chan-
ging and O H activity, will cause the neutral surface species to dissociate,
giving rise to a net surface charge. In this case, and O H are pdi's and the
pH at which the particle is neutral surface charge density a = 0) is called the
point of zero charge (pzc). Negative surface charge is created by an acidic
A
To use the surface complexing phenomenon in various technological applications,
it is essential to develop a predictive tool using a surface complexation model. For a
charged surface, the electrostatic field of the particle is known to affect the dis-
tribution of ions in the interfacial region and the probability of a given species
reacting with the surface. Therefore, the contribution of electrostatic free energy to
the surface binding of various species must be considered. To incorporate the sur-
face charge effect on the surface complexing process, a few models have been used
including the diffuse double-layer the constant-capacitance model
(CCM), the basic Stern model (BSM), the triple-layer model (TLM) and
recently the four-layer model (FLM) With these models, the effect of surface
charge on surface complexation has two aspects: (1) ion distribution (usually
distribution of ionic species given by exp where is the
charge of the adsorbing ions; elemental charge; electrical potential at the
location of adsorbing ion; the Boltzmann constant; and temperature), and
(2) electrostatic contribution = to the intrinsic (chemical) free energy
of adsorption. The TLM in combination with surface complexation has been used
extensively in developing a predictive tool for understanding the acid-base and
metal-binding properties of oxide colloids and is considered to be satisfactory for
most metal oxide and hydroxide systems It has been shown that only the
triple-layer model incorporated with proper surface complex reactions can repro-
duce both the surface charge and surface potential relationships
However, a recent comparison among various models showed that the
TLM cannot give a satisfactory description of all the major adsorption phenomena
such as charging curves, pH and salt dependency of cation adsorption, and
exchange simultaneously. In this study, we describe a different
approach based on, but not limited by, the existing triple-layer model shown sche-
matically in Fig. 2.
In this model, there exist two charge-free layers and a diffuse double layer. The
potential determining ions are located on a surface, designated by subscript 0, and
are responsible for developing surface charge The surface potential at this
plane drops linearly across the inner Helmholtz layer to at the inner
Helmholtz plane (denoted by where most specifically adsorbing ions along
with pair formers are located. (Note: A fraction of specifically adsorbing ions, such
as zinc ions, may locate on the surface and hence contribute to surface charge.)
Further out in the solution phase, a second charge-free layer is present due to the
adsorption of specific adsorbing ions along with ion pairs, with a further linear
drop in the potential across this layer but at a reduced gradient. The plane dividing
the compact layer and diffuse layer is referred to as the outer Helmholtz plane,
where the diffuse layer starts (denoted by d). The potential at the outer Helmholtz
plane coincides with the potential at the end of the diffuse double layer
Based on the classical electrical double layer theory, the potential drop across
the inner Helmholtz layer is related to the surface charge by
= -
where is the integral capacitance of the inner Helmholtz layer, given by
in which and are the dielectric constant in a vacuum and relative dielectric
constant, respectively, and is the thickness of the layer. The potential drop across
the outer Helmholtz layer, on the other hand, is given by
= - -
in which is the integral capacitance of the outer Helmholtz layer. Finally,
according to the Gouy-Chapman-Stern diffuse layer theory, the potential at the
outer Helmholtz plane is related to the charge density by
in which and are the bulk concentration and charge, respectively, of the symme-
trical electrolyte counterions in the diffuse layer. Based on the definition of differ-
ential capacitance of the diffuse double layer, in (2) can be expressed as
Equations (1)-(4) are the foundations of electrical double layer theory and are
often used in modeling the adsorption of metal ions at interfaces of charged solid
and electrolyte solutions. In a typical TLM, the outer layer capacitance is often
fixed at a lower value = 0.2 whereas inner layer capacitance can
be adjusted to between 1.0 and 1.4 It should be noted that the
plane model (TPM) is a variation of the classical triple-layer model, in which the
outer layer capacitance is not fixed. Although the physical presentations of the
TLM and TPM are identical as shown in Fig. 2, both involve a surface layer
an inner Helmholtz plane and an outer Helmholtz plane where the
diffuse double layer starts, a one-step protonation process 1 approach) is,
in general, assumed in the TPM, in contrast to a two-step protonation process
2 approach) in the Another distinct difference is that pair-forming ions
are assumed to be on the outer Helmholtz plane in the TPM but on the inner
Helmholtz plane in the TLM. In our study, the outer layer capacitance is allowed
to vary while the pair-forming ions are placed on the inner Helmholtz plane with a
complete set of surface complexation reactions being considered. Therefore, our
approach represents a of the TPM and
T o successfully model the surface charge and surface potential requires an accurate
description of surface complexing reactions. For illustration purposes, we use
lerite, a sulfide mineral, as an to show the surface site-binding model
The first step in the model development is to identify the surface species.
In the present case these are and with
senting surface. The second step involves considerations of all species present in
the inner Helmholtz plane for the system of interest. For sphalerite in a simple
electrolyte solution as in our example, the species includes -
and With all the species being identi-
fied, the complexing reactions associated with these species need to be specified and
the mass action law applied. These are summarized in Table 1 for the current example.
The constants in Table 1 are site-binding constants, and the subscript
represents concentration of surface species, which is related to the corresponding
bulk concentration through = distri-
bution. Replacing all surface concentrations in Eqs. (1)-(6) by this relation, and
rearranging them, the concentration of individual species can then be expressed in
terms of the concentrations of two basic surface species and the
bulk concentrations of ionic species and and the electrical poten-
tial at the locations of ions through corresponding site-binding constants. These
relations are summarized in Table 2 for ease of reference.
+ +
t
0 0
+ +
Ill
Ill t
Ill
t
3:
+
3:
+ + 0 0
a
Ill Ill
Copyright 1999 by Marcel Dekker. All Rights Reserved.
+ +
2
+ 3:
33 0
N
0
N
Ill Ill Ill
ace
an
Ill
Copyright 1999 by Marcel Dekker. All Rights Reserved.
After establishing the surface reactions (5)-(lo), Tables 1 and
the charge densities at arious planes can be derived by using speciation concentra-
tions. The surface charge, for example, is given by the summation of all ionic
species on the surface,
Assuming that surface species are distributed among the total number of avail-
able surface sites, one final condition that has to be satisfied is the
surface mass balance given by
+ + +
+ - + -
One of the major uses of setting up these fundamental relations and constraints is
to determine site-binding constants of surface complexing reactions These
constants are the foundation for various applications, such as predicting the surface
charge density and potentials, identifying adsorption characteristic of metal ions on
solids, and calculating the distribution of surface species. In our hybrid TLM
approach, the basic information required to solve this complex problem is a set
of experimental or (at least as many data points as the parameters
to be estimated minus balance constraints imposed, five for the above example).
In this chapter we use as our input variables. Although it is possible to
measure using ion-sensitive field effect transistors a common approach is to
use the zeta potential as an approximation of due to its access
experimentally. The possible error involved with this approximation was addressed
in detail earlier The effect of shifting the slip plane, where is measured, a
known distance Z into the solution on the fitted site-binding constants can be
examined using the exact solution of the Poisson equation
In the above example, a hybrid TLM was developed for a sphalerite in NaCl
solutions. This approach can be readily applied to other systems. T o illustrate its
versatility, experimental data on the boehmite water interface, generously provided
FIG. 4 Comparison between the predicted. using triple-layer and hybrid triple-layer
models, and measured surface charge density and zeta potential of in
solutions.
With the surface binding constants being determined the adsorption of ferrous
ions on sphalerite can be predicted with the use of commercial packages such as
MINEQL, developed by et o r S O I L C H M , developed by Sposito
and Coves Here we focus o n the effect of ferrous ion adsorption o n the
sphalerite zeta potential. Using the best-fit surface binding constants the
zeta potentials of sphalerite as a function of ferrous ion concentration at various
pH values are calculated and the results are given as solid lines in Fig. 5 , along with
the measured values under the same conditions. Here, the excellent agreement
between the two is evident. F o r all p H values studied, a similar effect of ferrous ions
o n of sphalerite was observed, the zeta potential increases initially with ferrous
I. SUMMARY
In this chapter, a novel hybrid triple-layer model was described. Three significant
differences between the hybrid and conventional triple-layer models were noted:
1. N o assumption was made with regard to 1 or 2 approach in contrast to
the conventional triple-layer model in which an assumption of iep
+ valid only for a system of symmetrical binding of
formation ions, is often involved,
2 The capacitance of the outer Helmholtz layer was allowed to vary (following
the differential capacitance of the diffuse double layer), in contrast to being
fixed in the conventional
3. The surface charge density was not required as input; rather it was calculated
from the fitted equilibrium constants of surface complexing reactions. An
improved fit was clearly demonstrated with the hybrid triple-layer model
compared to the conventional one.
The sphalerite-ferrous surface binding constants in solutions were deter-
mined. Applications of the hybrid triple-layer model in predicting electrokinetics,
surface charge density, ion adsorption, and surface solution speciation were
described and illustrated with examples.
I. Introduction
II. Simplified Lattice Ion Hydration Theory
I I I Nonequilibrium Electrophoretic Measurements
A. Electrokinetic phenomena-electrophoresis
B. Laser-Doppler electrophoresis
C. Results and discussion
IV. Extended Lattice Ion Hydration Theory
A. Case A (KF, RbF, CsF, CsCI, CsBr, and CsI)
B. Case B (LiF, NaF, LiCl, NaCI, LiBr, NaBr, Lil, Nal, and Kl)
C. Case C (KCl, RbCl, KBr, RbBr, and Rbl)
D . The sign of the surface charge
V. Oxygen Defect States
VI. Conclusions
References
in such systems the double layer is collapsed and electrostatic particle interaction
forces do not become significant except at very close separation distances.
Furthermore, Roman et [4] failed to show any direct experimental evidence to
validate the hypothesis. Conventional electrophoresis could not be used to study
soluble salts because of the difficulty in obtaining colloidal sized particles and also
because of the high conductivity of the solutions In the streaming
potential technique was evaluated without success. Later, Pizarro [6] had similar
problems in providing experimental evidence regarding the role of surface charge in
the flotation of alkali salt systems. Since then, the advent of the laser-
Doppler electrophoretic technique has made possible the characterization of certain
interesting systems including soluble salt systems that were hitherto not investigated
due to the lack of an appropriate experimental technique.
Based on nonequilibrium electrokinetic measurements for alkali halides by laser-
Doppler electrophoresis reported by Miller et al. the sign of the surface charge
of alkali halides in their saturated brines has been established. The surface charge
results thus determined are generally as expected from simplified lattice ion hydra-
tion theory. However, there are seven
KI, and RbI-when the results from nonequilibrium electrophoretic measurements
by laser-Doppler electrophoresis are compared with the predictions from lattice ion
hydration theory. Previously it was observed that in almost each case inaccu-
rate predictions occur when the difference between and anionic gaseous
hydration energies is small. However, more recent theoretical analysis has shown
that when the surface is considered to be partially hydrated, dependent on the
properties of the lattice ions such as the ionic radius (steric effect) and the hydration
number (hydration effect), the sign of the surface charge of all alkali halides except
for KC1 can be accurately predicted by this extended lattice ion hydration theory.
Finally, it has been shown by Yalamanchili and Miller that oxygen defects
present in the KC1 lattice can strongly influence its surface charge characteristics
and the negative deviation of KC1 from theory can be explained by the presence of
oxygen defects. Further, Veeramasuneni et showed that oxygen defects
deliberately introduced into positively charged can reverse the surface charge.
The fundamental aspects of, and some recent contributions to, the analysis of
surface charge of alkali halides in their saturated solutions include simplified lattice
ion hydration theory, nonequilibrium electrophoresis extended lattice ion
hydration theory and oxygen defect states and these issues are discussed
in the following sections.
If the hydration energy of the surface lattice cation is more negative than the
hydration energy of the surface lattice anion, as described by reactions then the
surface of M X will carry a negative charge. The converse is also true. The complete
analysis requires knowledge of the hydration free energies for gaseous ions [reac-
tions and the lattice energy for surface ions [reactions The hydration free
energy, for individual ions can be determined from the Born expres-
sion given by
where
N = Avogadro's number
z valence of the ion
= charge on the ion
= radius of the ion
dielectric constant of the water
The lattice energy for surface ions can be determined from Born-Mayer theory
as expressed by the equation
where
= lattice energy
constant for the crystal
= characteristic unit cell dimension
= a parameter related to the repulsive forces in the crystal that is evaluated
from compressibility measurements
Equations (4) and (5) can be used to obtain the hydration free energy of the
surface cation or anion [Reaction as follows:
Depending on whether the hydration free energy of the cation or the hydration free
energy of the anion is more negative, the corresponding surface ion will hydrate to
a greater extent, leaving the surface either positively or negatively charged as men-
tioned above. Of course, for symmetric uni-univalent nonreactive ionic solids hav-
ing a structure in which cation and anion are interchangeable, the terms are
There are four basic types of electrokinetic techniques that can be used to evaluate
the zeta potential under a given set of conditions. The common feature of these
electrokinetic measurements is the relative motion between a charged surface and
its surrounding fluid, with an electric field either applied or generated from the
movement. Perhaps the most familiar type and most widely used electrokinetic
technique is electrophoresis, which involves the movement of charged particles
under the influence of an applied electric field
The second technique, known as electro-osmosis, causes a fluid in contact with a
stationary charged surface to move in response to an applied electric field. The third
technique, which can be regarded as the converse of electro-osmosis. is known as
the streaming potential technique. The streaming potential involves the measure-
ment of the electric potential generated when a liquid is forced under pressure to
move in contact with a stationary charged surface. Finally, the sedimentation
potential technique is a method that is, in a way, the converse of electrophoresis.
A sedimentation potential arises across the suspension when charged particles settle
in a stationary liquid. (29) explained the basic principles of these electro-
kinetic techniques in a nutshell as shown in Table 2.
The electrophoretic technique is of particular interest because it has been used
for measuring the electrophoretic mobilities of alkali halide particles. In this case,
charged particles suspended in the electrolyte are transported toward the electrode
of opposite charge when an electric field is applied across the electrolyte. However,
viscous forces oppose this movement, and when an equilibrium is reached between
electrical attraction and viscous drag, the particle moves with a constant velocity.
This velocity depends on the voltage gradient, the dielectric constant, and the
particle's zeta potential. The particle's velocity divided by the electric field is gen-
erally known as its electrophoretic mobility.
The zeta potential can be determined from the electrophoretic mobility by
using the Henry equation in which the electrophoretic mobility is seen to be a
function of the double layer thickness,
where
= electrophoretic mobility
= dielectric constant
= viscosity of the medium
a = particle radius
= double layer thickness
For large values of for high electrolyte concentrations, Smoluchowski
modified the Henry equation as
=
This equation can be used to simply calculate the zeta potential from the measured
electrophoretic mobility. For small values of Eq. (7) becomes
For many years, electrophoretic measurements have been carried out using a
simple optical microscope to observe the individual particles moving between two
electrodes when a potential gradient is imposed. The actual electrophoretic mobi-
lities are determined by manually timing the movement of particles over a fixed
measuring distance. This microscopic method, which is still being used by some
researchers. suffers from several limitations:
. First of all, the method is very slow and tedious. For this reason, a trained
operator must devote considerable time and effort to accumulate a meaningful
set of data from which a mean mobility can be calculated. Operator fatigue and
eyestrain are real factors that could affect the outcome of these
The availability of low power lasers that provide bright and coherent light sources
together with sophisticated signal analysis instrumentation have led to the devel-
opment of an electrophoretic light scattering method to measure electrophoretic
mobilities of particles in suspension. This method combines the principles of elec-
trophoresis with laser-Doppler velocimetry, a method used for measuring the speed
of moving particles by analyzing the Doppler shifts of scattered light such
moving particles. The new generation of laser-Doppler instruments have trans-
formed electrophoresis a tedious technique to one that can be applied
routinely to obtain fast, accurate data without extensive operator training.
Doppler electrophoretic offer several advantages over classical micro-
scopic methods.
Advantages of Laser-Doppler Electrophoresis
Laser-Doppler electrophoretic measurements can be made very rapidly, in a matter
of few seconds rather than a few hours, and over a wide range of particle concen-
trations. Systems suitable for these measurements typically have particle concen-
trations of up to 10,000 ppm. These features are very important in systems that are
either chemically or physically unstable or settle rapidly.
Small values of the electrophoretic mobility [even smaller than 0.1
can be measured using the laser-Doppler electrophoretic technique. This aspect is
particularly important when the mobilities involved are small in magnitude and are
in the region close to the point of zero charge.
The traditional microscopic method is complicated by the electro-osmotic velo-
city due to the negatively charged glass walls of the cell. There is only one particular
region between the wall and the center of the cell, which is called the stationary
layer, where the particles move with a velocity that is due solely to their charge. It is
difficult to consistently focus the microscope in this layer with classical methods.
However, the laser-Doppler electrophoretic technique overcomes this difficulty by
accurately measuring the particle's velocity in the stationary layer.
Beam Selection
Shutter
\ Intensity M o n i t o r
which is between 100 and 1000 Hz, is always present in the signal, so that if the
particles have no net velocity the scattered light collected by the photomultiplier
will have exactly that frequency. In this way, particles with low or zero mobility can
be accurately measured.
The nonequilibrium measurements for alkali halide particles take only a few
seconds. Since the measurement is done on a large number of dissolving particles,
a distribution of electrophoretic mobilities is obtained, and a distinct eletrokinetic
characterization of the alkali halide particles is expected.
Figure 2 presents a schematic diagram of the AZ4 standard cell. This cell con-
sists of a 4 m m diameter quartz capillary. The platinum electrodes are in compart-
ments at each end of the cell. A semipermeable membrane separates the electrodes
from the suspension sample to prevent contamination of the electrodes by the
sample. The electrode chambers are filled with an electrolyte that is at least as
conducting as the sample itself. Polarization of the electrodes is reduced by simply
reversing the direction of the applied electric field a t about 2 Hz. The measurement
of particle velocity from many field reversals is added to obtain a good
noise ratio. A technique called "duty cycling", which allows the voltage to be
switched off for a short time between cycles is used to avoid heating of the suspen-
sion sample, especially for high conductivity solutions Switching off the
applied voltage allows the heat in the sample to dissipate. The problem of
tro-osmosis is overcome by measuring the velocity of particles at a point in the cell
where the electro-osmotic velocity is zero.
4. Experimental Procedure for Laser-Doppler Electrophoresis
Electrophoretic mobility measurements for various soluble salts were performed by
the laser-Doppler electrophoretic technique. The Zetasizer 3 with AZ4 standard cell
was used for measuring the electrophoretic mobilities.
Membrane
Electrodes
.2
Schematic diagram of the capillary electrophoresis cell for the Zetasizer 3.
(From Ref.
The soluble salts were sized by dry sieving, and the 65 x 100 fraction was
used in this study for nonequilibrium electrophoretic mobility measurements. The
electrode chambers were filled with an appropriate electrolyte of required ionic
strength. The electrophoretic mobilities of the alkali halides were measured in the
following manner:
1 . The salt particles were added to of high purity water in a beaker.
2. The of salt added (from 0.5 to 2.5 g depending on density and solu-
bility) was such that the Zetasizer 3 system would signal when an ideal level of
particles was present in the solution after injection.
3. The suspension was stirred gently for a moment 1 then quickly injected
into the cell, and the electrophoretic mobility was measured.
Each measurement required approximately 15 The starting size of the salt
particles is a critical factor. If the salt particles are too small, they dissolve before
the measurement can be completed. On the other hand, if they are too large,
sedimentation occurs in the capillary. Great care was taken to avoid the injection
of air bubbles into the cell. The experiment was performed at least three times for
each salt under the same conditions to determine the reproducibility of the mea-
surements.
In the case of reference equilibrium measurements which were possible for
and the salts (-400 mesh) were allowed to equilibrate at natural pH for at
least 60 min prior to electrophoretic measurements.
mobilities that was possible for and The results are presented in Table 4
and show the nonequilibrium electrophoretic mobilities to be equivalent in sign to
the equilibrium electrophoretic mobilities. Again, the electrophoretic mobilities
reported in Table 4 were determined from the composite electrophoretic mobility
data.
It should be mentioned that diffusiophoresis, the movement of particles in sus-
pension under the influence of an externally imposed solution concentration gra-
dient, can affect the nonequilibrium measurements of alkali halide particles. When
a particle of an alkali halide salt is dissolving in solution, it is possible that a
concentration gradient will be established in the volume adjacent to the outer
boundary of the double layer. The concentration gradient thus established can
influence the measured electrophoretic mobilities by laser-Doppler electrophoresis.
However, in the case of alkali halide particles, the slope of the concentration
gradient varies with time because of the rapid dissolution of these particles in
solution. Therefore, due to the small concentration gradient close to the particle's
surface, the influence of diffusiophoresis on measured electrophoretic mobilities is
expected to be of little significance. Furthermore, the results from equilibrium
electrophoretic mobility measurements for and (no diffusiophoresis) sug-
gest that at least the sign of the surface charge, if not the magnitude, as measured
by laser-Doppler electrophoresis is not affected by diffusiophoresis due to concen-
tration gradients in the outer layers of the double layer. Further, it is interesting
that the measured electrophoretic mobility for is negative, whereas that for the
thermodynamically stable is positive. It would appear that the hydration
of lattice ions does have a significant influence on the surface charge of alkali
halides and other soluble salts
KF
Chloride
KC1
Bromide
Iodide
KI
where is the distance between the centers of two charges, and The electro-
static pair interaction (interaction energy) between a charged atom and a polar
molecule presented by Israelachvili is
where is the dipole is the distance between the centers of charge of the
ion and a polar molecule, is the angle subtended by the dipole and the line joining
the centers of charge of the ion and the polar molecule, is valence, is electronic
charge, is permittivity of free space, and is dielectric constant As men-
tioned above, Eq. (17) gives the free energy for the interaction of a charge and a
B Significant; No; F
behavior
1.4 limited by
effect
NaBr
Lit
Source: Ref. 9
Alkali ion
Attractionat distance CD and
Repulsion at distance and
ion
Attractionat distance and
at distance AD and
Alkali ion
Attraction at distance CD and
at distance C
no
interaction energy
Alkali dipole Halide dipole -
Alkali
halides Attraction Repulsion Attraction Repulsion
NaCI
Source: Ref. 9.
2 at distance and
KC1 -62.61 -23.8 16.87 -69.54 -51.45 17.13 10.27 -24.05 -45.49
-58.52 -19.41 15.94 -6 1.99 -51.45 16.15 19.74 - 15.56 -46.43
-62.61 -20.8 15.3 -68. I -46.96 15.73 21.48 -9.75 -58.36
-58.52 - 19.95 14.54 -63.93 -46.96 14.91 20.6 -1 -52.48
-58.52 -15.59 12.58 -61.53 -46.96 13.19 16.54 - 17.23 -44.3
Att = attraction; Rep = repulsion.
Source: 9.
Table 10 compares the results obtained for the extended lattice ion hydration
theory with the results for the nonequilibrium electrophoretic mobility measure-
ments. It is evident that the extended lattice ion hydration theory provides for an
KF
Chloride
Bromide
Iodide
shown in Table 11. The higher oxygen levels in Optovac KC1 appear to be due to
the use of relatively less pure KC1 as a starting material for crystal growth and the
crystal growth procedure employed. For example, if the crystal is grown under
ambient atmosphere, oxygen in the atmosphere could react with the KC1 melt,
forming various defect states that are discussed in the following paragraphs.
Unfortunately, Optovac, Inc. was unwilling to share information regarding their
crystal growth procedures. However, it is most likely that one or all the reasons
mentioned above could account for the presence of higher levels of oxygen in
Optovac KC1 crystals. Nevertheless, it seems that the difference in oxygen levels
in the University of Utah and Optovac crystals is responsible for the difference in
the electrokinetic behavior of these KC1 crystals.
Oxygen can be present in the crystal in two different forms, namely,
and centers If the crystals are grown in an ambient atmosphere, KC1 melt
reacts with oxygen to form centers) as follows:
The symbol in (20) represents a vacancy in the crystal lattice. The introduc-
tion of excess potassium ions into the KC1 lattice creates F centers ( e at anion
vacancies) in the crystal lattice as mentioned in the previous section. Therefore
(20) can be rewritten as
It is clear from (21) that the number of defects will be reduced by the
presence of F centers in the crystal lattice, thereby forming the defect pairs.
The presence of and defects can be characterized by their typical
absorption frequencies in the region Figures 7 and 8 show the
absorption spectra of the Optovac KC1 and the University of Utah respec-
tively, in the region. The presence of defects in the Optovac KC1 can be
clearly seen from Fig. 7, as indicated by the band positions at 205 and 244 nm. On
the other hand, Fig. 8 shows that the University of Utah KC1 crystal is free of
defects, as indicated by the absence of any absorption bands in the region of
Fisher et reported that both and defects in KC1 have a
characteristic absorption band in the range of 200-250 It should be mentioned
that defects could easily be converted into defects if the crystal were
exposed to prolonged radiation even in the natural environment. However, neither
the University of Utah nor Optovac crystals show any characteristic absorption due
to defects. Typically, defects have multiple absorption bands a t 440, 286,
215, and and clearly these multiple bands are absent in Figs. 7 and 8,
which indicate that the University of U t a h and Optovac crystals had not been
exposed to long-term radiation, as they are free of defects.
Therefore it can be concluded that the absorption bands in Fig. 7 for the
Optovac KC1 are due to the presence of defects. while Fig. 8 indicates that
the University of Utah KC1 crystal is free of any oxygen defects, which was
expected, as the University of U t a h crystal was grown with extreme care to
avoid oxygen contamination.
In any case, the spectra for Optovac KC1 and University of U t a h KC1
and their nonequilibrium electrophoretic mobility results suggest that defects
strongly influence the surface properties and are responsible for the deviation of
normal from predictions of lattice ion hydration theory.
The concentration of oxygen defects was varied to examine their influence on the
electrokinetic properties of Optovac and reagent grade KC1 particles were
thermally treated in a nitrogen atmosphere containing less than 1 ppm oxygen in a
tube furnace. Similarly, the KC1 crystal grown at the University of Utah was
thermally treated in an oxygen atmosphere. Table 12 presents the results from
these thermal treatment experiments. It can be noted this table that when
the University of Utah particles were thermally treated at for about 27 h
in an oxygen atmosphere, the sign of the surface charge was reversed from positive
to negative due to oxygen incorporation into the crystal lattice. In the case of
Optovac and reagent grade KC!, the sign of the surface charge reversed from
negative to positive when these particles were thermally treated in the presence
of nitrogen having less than 1 oxygen. The sign of surface charge reversal in
the case of the University of Utah KC1 indicates that oxygen has diffused into the
crystal, as suggested by the increase in oxygen level from 61 to
147 9 ppm, thus accounting for the change in the electrokinetic behavior. These
results show that when KC1 is relatively free of oxygen defects, its electrokinetic
behavior is predicted according to lattice ion hydration theory. However, it is clear
from these studies that KC1 deviates from lattice ion hydration theory when oxygen
defects are introduced into the crystal lattice.
Similar observations can be made for both Optovac and reagent grade KC1
particles. Both Optovac and reagent grade KC1 (samples with high oxygen content)
exhibited a negative surface charge before thermal treatment. When these KC1
particles were thermally treated in the presence of nitrogen having less than
1 ppm oxygen, the sign of the surface charge of these particles reversed from nega-
tive to positive as expected from lattice ion hydration theory. Presumably, oxygen
present in the crystal has diffused out during the thermal treamtment due to low
oxygen pressure in the atmosphere above the bed of particles. In fact. the
Vis spectrum of a KC1 plate (Optovac, 5 x 5 x 1 mm), thermally treated at about
in presence of a nitrogen (oxygen 1 atmosphere, revealed that
oxygen indeed been removed from the KC1 plate. as indicated by the absence of
200
WAVELENGTH
In any case, spectra for University of Utah NaCI and Optovac NaCl
and the results from nonequilibrium electrophoretic mobility measurements suggest
that the negative surface charge of the University of Utah NaCI is due to the
presence of defects. It has been demonstrated that oxygen defects can be
created in alkali halides to them negatively charged, thus confirming the
previously reported results for oxygen defects in These results demonstrate
that oxygen defect states in NaCI and KC1 crystal lattices can change the sign of the
surface charge and alter the flotation response of these salts and that of other alkali
halide salts is expected.
Support by the DOE Basic Science Division, grant DE-FG-03-93ER14315, for this
work is gratefully acknowledged. Thanks are also extended to Dr.
Yalamanchili and D r . Yuehua H u for their suggestions and helpful discussions.
I. Introduction
I I . Historical Perspective
III. Structure Within the Adsorbed Layer
A. Experimental verification
B. Observations of hydrophobic interactions
IV. Adsorption with Attractive Electrostatic Interactions
V. Adsorption with Repulsive Electrostatic Interactions
VI. Adsorption of Mixed Surfactants
VII. Predictive Models of Surfactant Adsorption
VIII. Summary
References
Surfactants are used in a wide variety of applications such as ore flotation, cleaning,
polymerization processes, and pharmaceuticals and agriculture. The usual role for
the surfactant is to modify interfacial properties, whether they be
or interfaces. T o be effective in any of these applications,
however, the surfactant must adsorb strongly at the interface. In addition to its
concentration at the interface, the conformation of the surfactant at the interface is
also an important factor. While the influence of solution properties such as con-
centration, ionic strength, and pH on adsorption are well known, the
properties of the other phase also exert a significant influence.
Perhaps the most varied combinations occur in solid-liquid systems because of
the variety of properties possessed by the solid phase. The most obvious property of
the solid phase to affect surfactant adsorption is the or hydrophobic
STRUCTURE WITHIN
Intuitively one would expect that a surfactant would adsorb onto a hydrophobic
surface with its hydrophobic portion pointed toward the hydrophobic surface and
the hydrophilic portion pointed toward the aqueous phase. This adsorption is
driven by the "hydrophobic attraction" between these two hydrophobic moieties.
For a hydrophilic surface, adsorption should take place with the opposite orienta-
tion, as observed by Furlong and
When ionic surfactants adsorb onto an uncharged, hydrophobic surface, the
charge at the solid surface increases, giving the surface a more hydrophilic char-
acter. As the adsorption continues, the concentration of charged groups a t the
outer edge of the adsorbed layer becomes great enough to block any further
adsorption due to the electrostatic repulsions between the adsorbed surfactant
and the surfactant molecules approaching the surface. This is opposite the expecta-
tion for adsorption onto hydrophilic surfaces In this latter situation. surfac-
are expected to first adsorb with the hydrophilic portion in contact with the
solid surface, frequently as a result of electrostatic attractions. This leaves the sur-
face covered with an adsorbed layer with its hydrophobic ends exposed and thus
looking very much like a hydrophobic surface. Subsequent adsorption will then
occur onto this adsorbed layer, with the hydrophobic ends of the next layer of
oriented toward the first layer. The final result, at saturation, is an
adsorbed bilayer at the surface. Further adsorption is blocked by electrostatic
interactions between the hydrophilic ends of the
FIG. 1 Amount of sodium dodecyl sulfate adsorbed onto polystyrene latex surfaces as
a function of surface charge density. Hoeft et Hoeft and
Zollars Brouwer and Zsom Gwin
F MIX
Since hydrophobic interactions appear to play a significant role in determining
surfactant adsorption it is possible that these interactions may also alter the manner
in which a mixture of surface-active materials might adsorb. In a very early inves-
tigation, Orr and found that when mixtures of anionic and nonionic
surfactants were adsorbed onto a polymeric latex particle there was a decrease in
the adsorption area per molecule for the anionic surfactants. Hay et how-
ever, observed a decease in total adsorption for sodium dodecyl sulfate and
octyltetraethylene glycol onto graphitized carbon.
In general, however, it is observed that the presence of a nonionic surface-active
material will increase the adsorption of ionic surfactants. This was directly
observed using in situ sum-frequency spectroscopy for the adsorption of sodium
dodecyl sulfate and by et where more highly ordered films
are observed for dodecanol-sodium dodecyi sulfate mixtures than for pure sodium
dodecyl sulfate. The reason for this is usually given as a screening of the electro-
static interactions between the head groups on the sodium dodecyi sulfate by the
nonionic material, allowing for a closer approach by the head group while still
maintaining strong hydrophobic interactions within the adsorbed layer. This, in
turn, allows more compact adsorbed layers and greater adsorption.
A possible reason for the decreased adsorption of some nonionic materials has
been suggested by et In a study of the adsorption of sodium
dodecyl sulfate and cellulose onto hydrophobized silica sur-
faces, they found that the amount of adsorbed cellulosic derivative was reduced
when sodium dodecyl sulfate is also adsorbed. A possible reason for this is again
the effect of hydrophobic associations. The sodium dodecyl sulfate, being more
surface-active, will tend to adsorb onto the hydrophobic surface more strongly. It
will also associate with the cellulose in solution, rendering this
nonionic material more With both the adsorbing surface and the
10
.
5
-0
0)
From reviewing all of these sources a picture begins to emerge about the current
state of knowledge about the adsorption of ionic surfactants onto hydrophobic
surfaces. For adsorption onto an uncharged surface the adsorbed associ-
ate with the surface via a hydrophobic interaction. This results in molecules being
adsorbed in a tail-down conformation. At low concentrations the molecules appear
to lie flat on the surface but become progressively more vertical as the concentra-
tion increases. In this vertical position significant lateral interactions between the
adsorbed molecules occur, leading to increased adsorption and, ultimately,
layer formation. On surfaces that themselves possess a structure, such as graphite;
this may also induce a structure in the adsorbed layer (adsorption epitaxy?) along
the individual grains of the surface as observed by and Ducker
When ionic adsorb onto a hydrophobic surface that also contains
charged groups of the same sign as the surfactant, electrostatic interactions begin to
play a role. At low surface charge densities the initial adsorption will take place as
in the case of adsorption onto an uncharged surface, with the surfactant lying
nearly flat on the surface. As the surface charge density is increased, the adsorbed
surfactant will be forced into a more vertical orientation. In this orientation the
potential exists for greater adsorption due to the opportunity for more hydrophobic
associations via either surface-surfactant or surfactant-surfactant contacts. At still
higher surface charge densities the electrostatic repulsions so strong that no
observable adsorption occurs.
For the adsorption of ionic onto hydrophobic surfaces of opposite
charge, a couple of different scenarios are possible. At low surface charge densities,
the initial adsorption may take place primarily through hydrophobic associations
as in the case of an uncharged surface, with an occasional adsorbed molecule
occurring in a head-down orientation due to electrostatic attractions between the
surface group and the surfactant head group. Once the surface charge has been
neutralized, subsequent adsorption should occur as if on an uncharged surface.
This should result in coverage with the potential for some structuring
of the layer around the sites where a molecule has adsorbed in the
down orientation. At high surface charge densities, the initial absorption should
occur predominantly via electrostatic attractions between the surface groups and
the surfactant head groups. As the adsorption proceeds, the surface will become
more and more hydrophobic due to the tails of the adsorbed extending
toward the aqueous phase and charge neutralization. Ultimately the surface will
become so hydrophobic that another adsorbed layer will start to form due to
hydrophobic associations between the surfactant tails. This ultimately leads to
the formation of an adsorbed bilayer,
When more than a single surfactant species is present, the situation becomes
more complex again. It seems clear, however, that whether the presence of other
surface-active species increases or decreases the tendency for adsorption depends
I. Introduction
II. Metal Complexing Ability of Humic Acid
A. Acid-base properties of humic acid
B. Metal-humic acid complexation
I I I . Adsorption of Metal Ions by Immobilized Humic Acid
A. Preparation of composite biopolymer adsorbent
B. Metal adsorption ability of HA-M
IV. Adsorption of Metal Ions onto Insolubilized Humic Acid
A. Preparation of insolubilized H A
B. Effect of insolubilization on metal-HA complexation
V. Practical Application of IHA
A. Cadmium-midgut gland complexation
B. Removal of cadmium from M G by competitive adsorption
References
0 0.5 1.0
Added amount of
where represents acidic groups of type on HA. The acid dissociation con-
stant is defined as
where and are the numbers of type 1 and 2 acidic groups, respectively. Their
sum should be equal to the total number of acidic groups determined from the
conductometric titration,
Magnetite
acidic groups on HA, which serve as the metal complexing sites, is higher than
that of most conventional inorganic adsorbents.
The acid dissociation constants of type 1 and 2 acidic groups differ by about two
orders of magnitude. Therefore, the complexation of metal ions with HA
can be written as
+ (- i= or 2
where represents the metal ion. The complexation constant is
defined as
where
where represents the fraction of type i acidic groups complexed with metal ions.
The number of metal ions complexed with 1 of HA, is expressed as
where and represent the numbers of type 1 and 2 acidic groups available
for metal complexation, respectively.
. Preparation of Composite
An aqueous solution containing HA, acid (AA), and activated carbon pow-
der (ACP) was well mixed. It was extruded, dropwise, onto a plastic and dried
at Then it was in an excess of solution for 1 The membranes
of HA-entrapped calcium alginate had a diameter of about 5 mm and a thickness of
about 0.3 mm. Using this immobilizing method, we obtained a membranous adsor-
bent (HA-M) consisting of 40% HA, 20% AA, and 40% on a dry weight
basis.
The fraction of acidic groups complexed with metal ions, 0, and the number of
metal ions complexed with 1 g of biopolymer, can be expressed as
where
and
where represents the number of acidic groups available for metal complexation
on 1 g of biopolymer. The complexation experiments were carried out
using the membranous calcium alginate gel (containing no H A or ACP) as the
adsorbent. The equilibrium parameters for the metal-AA complexation,
and were determined in the same manner as those for HA, and they listed
in Table 1 . The number of carboxylic groups available for metal complexation on
AA, is about 30% smaller than the total number of carboxylic groups on AA
determined from the conductometric titration, et al. reported that
the primary mechanism of of alginate was the dimerization of
sequences, with specific site binding of calcium ions between alginate chains.
This process reached completion at a calcium level of about 30% of the total cation
requirement of the alginate. Therefore, it can be considered that about 30% of the
3.04 -3.17
Lead 3.02 -5.13"
Data from Ref. 8.
The heating of raw humic acid above a certain temperature causes dehydration
(condensation) between the acidic groups and lowers its solubility in water (Fig. 8).
Therefore we tried to prepare an HA by heating (IHA).
reparation of
acid was pretreated with excess solution and converted to a calcium
salt it was made insolubilized by heating at for 1 Heating at
higher temperature or for a longer period reduced the acidic groups on HA to a
considerable extent The HA thus obtained (IHA) was ground in an agate
mortar with a pestle and washed repeatedly with 1 solution and hot
distilled water. Finally, IHA was converted from the hydrogen form to the sodium
by soaking in 1 solution, rinsed with distilled water, and dried at
IHA prepared in this way was virtually insoluble in water up to pH 10 at
and there was only 7% dissolution even in boiling water
The internal organs of most marine organisms contain valuable materials such as
DHA, proteins, and essential acids However, one of the "big
three" toxic heavy metals, cadmium. is concentrated in the internal organs of
organisms throughout the food chain. Many studies have been carried
out to study the accumulation of cadmium in the internal organs of edible marine
organisms such as mussels snails crabs whelk scallops
and squid For example, cadmium concentration in the mid gut gland
(MG) of the scallop which is one of the most popular
edible shellfish, is as high as 200 on a dry weight basis For the effective
use of these valuable materials, the cadmium removal process should be estab-
lished. Cadmium in the internal organs can be removed by treatment in acidic
solutions. Cadmium ions are released as a result of the protonation of
ing proteins in the organs. Therefore, the proton concentration of the acidic solu-
tion must be kept at a high level, and the subsequent treatment of
containing, low pH wastewater is required. Moreover, treatment in a strongly
acidic solution may denature these valuable biomaterials.
In a recent study, we proposed a new method (competitive adsorption method)
for the removal of toxic heavy metal ions from acid-sensitive and
applied it to the removal of cadmium from the M G of the scallop
In this method, a competitive adsorbent (IHA) is added to a weakly
acidic solution containing M G to lower the cadmium concentration in the liquid
phase. The lowering of the cadmium concentration in the liquid phase promotes the
release of cadmium ions from M G , and the released ions are immediately
adsorbed by IHA. In this manner the competitive adsorption method enables the
simultaneous removal of cadmium from the scallop and from the liquid
phase under mild acidic conditions.
where and represent the amounts of M G and IHA added to the system.
respectively, and are the amounts of cadmium ions compiexed with M G
and IHA, respectively, and is the total cadmium concentration in the system.
Sliced M G was used in these competitive adsorption experiments. Since 46% of the
cadmium binding protein in M G was constantly dissolved in the liquid phase, we
modified Eq. (13) to
where represents the cadmium concentration (the sum of the free cadmium ions
and the cadmium ions complexed with the dissolved protein) in the liquid phase.
The pH dependence of in the competitive adsorption systems containing 30
3
dm of M G and 1.0, 2.0, or of IHA is presented in Fig. 9. For compar-
ison, the results of acid treatment (in the absence of IHA) are also plotted in Fig. 9.
A successive approximation method was applied to numerically estimate from
(1 and (14). The experimental data (symbols) agreed well with
the predicted values (solid lines). The pH dependence curves of have minimum
values at pH 5. The competitive adsorption method enabled the simultaneous
removal toxic cadmium from both the liquid and M G phases at pH 5. The
increase in in the pH 5 region was attributed to the increase of free cadmium
ion, while the increase in the p H 5 region was attributed to the increase in the
amount of cadmium- protein complex dissolved in the liquid phase.
I. Introduction
A. The use of hydrous metal oxides as adsorbents
B. Surface characterization
C. Removal of heavy metals from solution
II. Experimental Techniques
A. Reagents
B. Background electrolyte
C. Metal ion solutions
D. Colloidal substrates
E. Precipitation, adsorption, and coprecipitation studies
F. Determination of metal ion removal
111. Precipitation, Adsorption, and Coprecipitation
IV. Comparison of H C O and H F O as Substrates
A. Adsorption
B. Coprecipitation
V. The Influence of Competing Metal Ions
A. Cr(III) with H F O
B. Zn(II) with H F O
C. Ni(II) with H F O
D. Zn(II) with HCO
E. Ni(II) with HCO
VI. The Influence of a Complexing Ligand
A. The influence of ammonia on the solution chemistry of aqueous
heavy metals
B. The influence of ammonia on the adsorption of Zn(II)
C. The influence of ammonia on the coprecipitation of Zn(II)
VII. Modeling of Adsorption and Coprecipitation
A. Introduction
The removal of heavy metal ions from both natural water supplies and industrial
wastewater streams is becoming increasingly important as awareness of the envir-
onmental impact of such pollutants is fully realized. In particular, the likelihood of
such metal ions precipitating out of solution and/or coating other materials can
have a profound effect on both aqueous and nonaqueous environments. There is
considerable evidence in the literature that the primary mechanism for transporta-
tion of metal contaminants in aquatic systems is the movement of suspended par-
ticulate material containing the adsorbed pollutant metals It is also known
that a strong correlation exists between the concentration of trace metals in the
(aquatic) environment and the extent to which those metal ions adsorb onto col-
loidal substrates present in the environment A similar correlation between the
concentration of trace metals in the (aquatic) environment and their precipitation
behavior is not so clear. There is, then, a well-founded need to study
related phenomena in order to understand and predict the behavior of toxic metals
in the environment.
Industrial companies are increasingly being required to account for the fate of
all chemical species, whether deliberately added or present as by-products, at all
stages of industrial, mining, and manufacturing processes. The processes of pre-
cipitation, adsorption, and coprecipitation. apart from directly controlling the eco-
nomics of many chemical processes, are also often involved in the cleanup of
industrial wastewater Thus there is also a well-founded need to study adsorp-
tion-related phenomena in order to understand and predict the behavior of toxic
metals in industry.
Metal ions that are prevalent in industrial effluent and have a considerable
environmental impact include and Waste solutions of
these metals are routinely found at levels that are in excess of acceptable disposal
limits. Such waste solutions are produced, for example, by electroplating and metal
pretreatment printed board manufacture
leather tanning and wool scouring and anodizing among
others These three metal ions are studied in detail in this chapter.
3. Electrokinetic Properties
The points of the two substrates used in this study were
by using a DELSA 440 and found to be at pH 8.2 for H F O and at pH 9.2 for
HCO. Both and are completely adsorbed at pH values considerably
less than either of these values (see Figs. 3-51, indicating that the surface is
positively charged during adsorption. Electrostatic forces, in this case, are unfavor-
able to adsorption, and there must be a strong specific interaction for adsorption to
occur.
All reagents used to prepare solutions were of analytical reagent (AR) quality
unless otherwise specified.
. Metal lo
Metal ion solutions were prepared from the nitrates to provide stock solutions that,
when diluted, produced concentrations of 50 ppm with respect to the adsorbing or
coprecipitating metal under investigation and with respect to the metal ion
Two types of solid substrate surfaces were used in this study: amorphous hydrous
oxide (HFO) and amorphous hydrous oxide (HCO). Both
substrates were prepared by slowly increasing the acidic pH of either an
nitrate or nitrate solution of concentration with respect to
the metal ion in question. Metal ion adsorption or coprecipitation experiments
using amorphous hydrous metal oxide substrates are generally described in terms
of the concentration (ppm) of metal ion that was used to form the colloidal adsor-
bent rather than its corresponding specific surface area
The extent to which metal ions were removed from solution was assessed by per-
forming experiments carried out in three modes: direct precipitation, adsorption,
and coprecipitation. In all cases, the initial p H of the starting metal nitrate solu-
t i o n ( ~w
) as low pH 3.5 or less). Atomic absorption analysis confirmed that no
appreciable quantities of insoluble hydroxide precipitates were present in the metal
nitrate solutions at these low pH values.
Direct precipitation profiles of single metal ions were obtained by increasing the
pH of nitrate solutions of the metal ion in question. The pH and extent of pre-
cipitation were measured after equilibration at the final pH for 40 The experi-
ment was repeated with a fresh sample for each p H measured.
In the adsorption experiments, the colloidal substrate (either H F O or HCO) was
prepared by increasing the pH of an or nitrate solution of
concentration with respect to metal ion, and allowed to equilibrate for a
period of 5 When studying the removal behavior of single metal ions, the
metal ion of interest was then added to the pH-adjusted colloidal suspension and
equilibration was allowed to take place, after which time the equilibrium pH of the
suspension and extent of adsorption were measured.
When studying the removal behavior of the various metals in multiple
systems, the metal ions under investigation were premixed in a separate container
before being added to the reaction vessel containing the pH-adjusted colloidal
substrate. This procedure was followed in order to avoid the possibility of prefer-
ential adsorption of any of the metal types that may have occurred if the metal ions
were added to the suspensions separately. Again, the experiment was repeated with
a fresh sample for each pH measured; equilibration time was 40
In the coprecipitation experiments, the p H of the solution containing both the
or nitrate solution and the metal under investigation
was increased. Again, the pH and extent of coprecipitation were measured after
The metal ion removal achieved during the various precipitation, adsorption, or
coprecipitation experiments was evaluated by determining the concentration of the
metal remaining in solution after equilibration had taken place. Each sample for
analysis was removed from the reaction vessel, filtered through a 0.22 Millipore
nitrocellulose filter into a glass vial containing two drops of concentrated nitric
acid, and analyzed using a Varian atomic absorption spectrometer
with a nitrous flame. The quantity of metal ion precipitated,
adsorbed, and/or coprecipitated was assumed to be given by a simple mass balance,
the mass added less the mass found in solution.
The influence of a second or third metal ion on the adsorption and coprecipitation
removal of metal ions is given in this section, initially using HFO, and then HCO,
as the adsorbing or coprecipitating colloid.
90 HFO
alone
with Ni
with Zn
with Ni, Zn
alone
with
with Ni
with Cr, Ni
is, however, a slight deviation at pH values less than 7. Again. the presence of
enhances metal ion removal; however, the enhancement is small and occurs
only at low degrees of metal ion removal.
alone
with
The easiest way to use the constants presented in Table 1 is to graphically plot
concentrations of the various species as a function of pH. Given the number of
species involved, this is best performed by using three separate graphs, as shown in
Fig. 18. This figure will be referred to when discussing the species likely to be
present at any given pH value. For example, at 6, the species in highest con-
centration are and (Fig and these species are relatively inde-
pendent of pH near that The species (Fig. 18b) and
(Fig. are not present in high concentration, but the shapes
of their curves indicate that they become more and more important as the is
increased.
A large number of models exist for modeling the adsorption of metal ions onto
colloidal substrates, and most of them can be adapted to include coprecipitation.
They all involve "fitting" parameters the equilibrium constant for adsorption
or the chemical specific free energy for adsorption) in order to match experimental
and theoretical results. Therefore, there remains some skepticism over the ability of
current theories to fully account for the adsorption process. Nevertheless, these
"fitting" parameters are generally held constant as a function of pH. Thus the trend
in adsorption (or coprecipitation) with pH is a guide to the success of the model. It
is not the purpose of this chapter to detail all the adsorption models or to compare
them. Rather, one model is used to test its robustness in modeling precipitation,
The James-Healy model has been used extensively to model the adsorption
behavior of a series of aqueous heavy metals onto a variety of substrates. However,
there is still considerable uncertainty in many of its underlying concepts. The model
is based on determining the free energy of adsorption for all potential adsorbing
species i For example, in the study of adsorption in the presence
of a background electrolyte, such species are the free metal ion, its hydro-
lysis species, and its nitrate complexes. For each species, is divided into
three components: and The first component,
is associated with the replacement of part of the solvation sheath that surrounds a
charged ion in solution, with a solid surface. The second component, is
associated with bringing a charged ion into close proximity to a charged surface.
The third component, is associated with specific chemical interactions
not accounted for in the first two terms. Such specific interactions may include
hydrogen bonding, crystal field stabilization, hydration forces, and chemical bond-
ing. It was also originally envisaged that it could account for errors in the calcula-
tion of and but aspects of all three terms remain open to further
definition.
The James-Healy has been used to model the adsorption of
and onto of onto of onto
amphoteric polystyrene latexes of onto and goethite of
and onto synthetic goethite of and
complexes onto and of
onto and of mercury complexes onto In each of these cases, the
model was used in the original form proposed by James and Healy although
recent advances in physical chemistry can improve the algorithms used to calculate
the various free energy components.
There has been little improvement to the calculation of , by James
and Healy despite the simplicity of the approach, which involves a discrete,
rather than a continuous, change in solvation energy. Although this may be con-
sidered a weakness, it is appropriate because there is increasingly more evidence
that interfacial properties are of a discrete nature
By contrast, there is considerable evidence to suggest that the Nernst equation,
used to calculate the surface potential and thus is inappropriate in the case
of hydrous metal oxides et al. proposed the use of zeta
potentials rather than surface potential to calculate (and also
and this forms the basis of the modified model used throughout this chapter.
Of the greatest concern is the uncertainty associated with the origin of
This term is used largely as a fitting parameter, although there are some indications
of its phenomenological basis. and James for example, showed
that complexation of a metal ion with a unidentate such as chloride,
The adsorption and coprecipitation profiles for removal using are very
similar. There is a small but significant difference between adsorption and
cipitation for and a large difference for using the same substrate. In
all cases, coprecipitation was as efficient as or efficient than adsorption, which
was in turn more efficient than precipitation alone.
In all cases, was removed at a lower pH than which was removed
at a lower pH than consistent with their hydrolytic behavior.
Hart, (1982).
A. Jenne, Trace element sorption by sediments in soils: sites and processes, in
Symposium on Molybdenum on the Environment and Peterson,
Marcel Dekker, New York, 1977, Vol. 2, 425.
Krauskopf, Geochim. 9:1 (1956).
P. Hartz, and Environ. (1987).
Waste Treatment with and Other
Data Corporation, Park Ridge. 1977.
Introduction to Waste Water Treatment Processes, Academic, New
York, 1977.
P. I. H. Harding, and D. Sep. Technol.
(1992).
C. A. Money and Cranston, Reductions in Tannery Effluent by
Implementation of the Process with Liquor Recycling and a Rapid
I.Introduction
II.Theoretical Considerations and Modeling
III.Experimental Considerations
IV. Investigations of Adsorption of Various Ions
A. Adsorption of H+ and OH-
B. Adsorption of cations
C . Elemental anions and o x o anions
D. Complex ions
E. Organic ions
V. Concluding Discussion
References
THEORETIC
In quantitative descriptions of adsorption it is usually expected that equilibrium is
reached, and thus the laws of thermodynamics may be used to describe the system,
so the equation
0
where is the free enthalpy of reaction, must be valid for any reaction, leading
to
where
standard free enthalpy of reaction
standard enthalpy of reaction
standard entropy of reaction
R gas constant
reaction temperature
equilibrium constant
product of activities of substances i
The activities of the substances are related to the concentrations by
where/, is the activation coefficient (in the case of ions, for ion pairs).
The following relation between the activity coefficient and ionic strength J is
derived from the extended theory for dissolved ions.
will cause the greatest differences between different adsorbed ions, leading to
physisorption in the case of small values and to chemisorption, with the formation
of chemical bonds, for large values. Partial models may be constructed for weak or
strong interactions.
+ OH- + (8)
0 2 4 6 8 1 0 1 2 1 4
value
-1
1 2 3 4 5 6 7 8 9 1 0 1 1 1 2
value
obtained directly by the titration method (Fig. 3a) and from model calculations
(Fig. 3b)
reports that surface charge pH, and its derivative with
respect to pH vs, pH can often be described on basis of pseudo.
homogeneity. In the case of transition (and of gibbsite) this may, at least
partially, be caused by the fact that these oxides are remarkably soluble, which will
lead to more homogeneous surfaces, as the most active centers, having the highest
energies, will dissolve first.
The adsorption of protons and hydroxyls necessarily leads to charges on the
surface, which are balanced by counterions of opposite charge in the solution.
Assumptions about the structure of the charged layers, which also contain other
adsorbed species, lead to different models double- or triple-layer models).
The point of zero charge (solution pH at which the surface is uncharged) is
related to the constants H I ) and by
pzc = - HI) +
where
provided that no surface sites are occupied by other ions. In the presence of back-
ground ions, an extrapolation to zero concentration is possible. H2) and
can be disregarded at this pH; in most publications they are not used
at all, but they are introduced here as they are to be expected theoretically. The pzc
can be determined independently from electrokinetic measurements provided the
isoelectric point (iep) is equal to the pzc. This reduces the number of constants to be
derived from the titration experiments, or it gives a tool to test the reliability of the
titration provided the pzc results are trustworthy themselves.
A t this point one may ask whether a second constant can be defined in analogy
to the ion product in the liquid water phase:
0.3 M
--
r
Triple-layer model None
a Triple-layer None
Nat. prod.' Electrophoresis
prod.' potential
prod, = product
adsorption of other ions. Table 2 contains some adsorption constants for protons
and hydroxyl ions. Other cannot be compared directly, as they d o not
contain electrostatic terms, they are not intrinsic values.
The concentration of surface centers is taken to be equal to the number of
hydroxyl groups in a monolayer o n the surface other values obtained
by acid-base titration tritium exchange and equilibration experiments
combined with calculations ("long time Table 3 gives
a n overview of numbers of centers determined by various researchers. Most values
are found between 1.0 x and 3.2 only that of [7]
(1.7 x clearly lower, possibly because of the back-titration tech-
nique. So the assumption that the number of centers is equal to the number of
surface O H groups seems to be acceptable, though in principle different of
Number of Centers
Ref. Oxide .Method
Titration
exchange
Exp + mod
Exp + mod
Exp + mod
OH groups
OH groups
OH groups
The numbers of publications about the adsorption of different cations onto alu-
mina are very different depending on the practical importance of the adsorption
process. A lot of is done with hazardous ions, or ions of known catalytic
activity, with the ions themselves or after their reduction to finely dispersed metal
particles. the other the literature about other ions, being of scientific
interest only. is limited.
Any specific adsorption of ions takes place in competition with that of or
provided it takes place at the same centers:
where denotes a free surface center. Here both surface and specia-
tion of the metal ions in solution must be known, with being equivalent to
covalent forces are absent, electrostatic forces will be dominant. Another effect
may be of greater importance here, too: As there is no evidence of strong specific
interaction, one has to consider the corresponding anion also; if the anion is
adsorbed specifically, in the inner layer of the double- o r triple-layer system,
the corresponding cation will be enriched in the outer compensating for
the charge.
0 NaCl
0.001 NaCl
NaCl
0.05 M NaCl
of ions between the surface and the liquid phase in 0.05 M Assuming a BET
surface of 100 this value on the order of looks reasonable
at 8 and a liquid concentration of 0.1 In contrast, shows a 100-fold
increased distribution coefficient, and a 10"-fold one.
The adsorption of alkaline earths onto was described by Huang and
Stumm Examples of adsorption depending on pH are shown in
Adsorption onto a slightly negative surface may lead to a reversal of the surface
charge, the surface speciation curves are shifted to higher p H values by
alkaline earth ions compared to those of pure alumina. From the corresponding
experiments it was concluded that ions are chemisorbed, or at least bound in
an inner surface layer. The logarithms of equilibrium constants for the
reaction
are: Mg, -5.43;Ca, -6.06; Sr, -6.27; Ba, -6.59. Adsorption on fresh alumina gel
was described by et As expected, the adsorption tendency is the
same as for alumina, though a direct comparison is not possible because no specific
surface can be defined. The results are discussed in terms of values.
Measurements with Ba and on show that is constant for
concentrations below 1 x while it increases linearly with log above
this concentration. At low concentration the value for both is 10, increasing
about one unit per decade of concentration for Ba above
SOM-OH
jump experiments with sorption data and found log of deprotonated surface
groups with of 9.71 and 9.59 L
Hachiya et al. found a linear relationship between the logarithms of
intrinsic adsorption constants and the hydrolysis constant of the metal ions for
Cu, Pb, Co. and They described pressure jump experiments with
showing fast and a slow process attributed to two surfaces sites - with
and The relaxation times decreased with metal ion
and proton concentrations. The kinetic constants for adsorption and desorption are
(7.24 2.0) and (3.1 0.9) The order of the
adsorption rate constants corresponds to that of the water release reaction from the
aquo complex in the liquid phase.
The and values (negative logarithms of the adsorption constant and the
solubility product) are given in Table 4. In the case of Co and it is difficult to
separate the effects of precipitation and adsorption.
Further results were obtained with and et al.
discussed adsorption isotherms of and C o ions on aluminum
hydroxide at different temperatures and derived the thermodynamic values shown
in Table 5.
2(= OH) + +
The constants are given in Table 6. Copper adsorption at "low pH" in contrast to
that of Ca, Mg, and is mentioned by Different values for
the adsorption of on hydroxide gibbsite and boehmite (5.7) may
be caused by the different numbers of centers.
The second transition ion group contains Zn, Cd, and which are known to
be strong poisons. Besides the works cited already. and give
curves for the removal of of different initial concentrations by different
amounts of Brady shows the temperature dependence of Cd and Pb
adsorption. The speciation of Zn and Cd in the liquid phase necessary for mode!
calculations is reported in Refs. 4 and 67 with remarkable differences, especially
for the species. Liu et [4] observed a linear relation of the logarithms of
adsorbed and dissolved amounts of and et al. report
extremely slow equilibration of with alumina in the presence of
and Evans studied the adsorption of Zn and Cd ions from very dilute solutions
on
hydroxide was used as an adsorbent in a of publications.
Coprecipitation of and was by Simon et al. in the
absence and presence of Shiao et demonstrate that cadmium
distribution coefficients between liquid and surface are practically uninfluenced by
the concentration, while has a remarkable influence. According to
et and Packter and Derry coprecipitation leads to the forma-
tion of Zn Besides investigating various ions and co-
workers describe sorption on especially that of Zn, at
different temperatures. They observed isotherms and derived thermody-
namic data. Kinetic studies were carried out by Simon et using
graphic techniques. In their experiments equilibrium was reached after 3 with
two consecutive first-order reactions, in contrast to the findings of Lo and
Anions have attracted much less interest than cations. Some of them, such as
nitrate, are not specifically adsorbed, they are not able to replace surface OH
groups. is adsorbed to a maximum of 2.7 at pH 2.5 from a
0.1 M solution Figure 14 shows the pH dependence of adsorption for
some anions. It must be stressed that, in contrast to the figures for cations, the
adsorption of anions cannot be compared directly, as the measuring conditions are
In principle, all ions form complexes in solution, normally aquo complexes. Here
we discuss complexes with ligands other than They are treated together,
whether cations or anions, as their charge depends on the number and charge of
the ligands for a given central ion.
The inorganic ligands that were investigated are and Brunelle
gives an overview of the tendency of elements of Groups 8, and to form
or anionic complexes and uses the adsorption of or anionic com-
plexes for catalyst preparation
Amino ligands, not changing the charge, were used in the following investiga-
tions. Spielbauer al. described aquo complexes, using Pd speciation
in the liquid phase described by Rasmussen and Jorgensen containing
for comparison. They interpret their results on the basis of the triple-layer
model (TLM) of Davis unfortunately without experimental values
above p H 11.5, where the adsorbed amount should decrease according to the
model calculation.
Contescu et al. gave adsorption isotherms for for pH 8 and
and one for pH 10.6 was given by Schwarz et al. Fuerstenau and
showed the adsorption of Cu, and C o (1 x M) at
,
c( moles
Organic ions have attracted the most interest in the form of and anionic
surfactants, though other organic ions may be a problem too.
and Motekaitis reviewed the adsorption of aliphatic acids (maleic,
citric, and acid), salicylic acid, and a number of sulfone acids
acid)-
and nitrogen-containing complex ligands (NTA, EDTA, DTPA). Equilibrium con-
stants are given for the relevant equilibria. An example of the application of organic
acids (oxalic acid, salicylic acid, is aimed for better dispersion of Pd on a
catalyst carrier by coadsorption The adsorption of a surfactant
(TTAC = chloride) is shown in the form of an
adsorption isotherm for pH 10 by Huang and Somasundaran that of anionic
surfactants and sodium octylbenzenesulfonate) by
et al. and et (Fig. 17). The load of
2
2x molim shows that the whole surface is covered, assuming that
an area of 0.28-0.83 nm2 is covered by one molecule. The anionic
show the same normal pH dependence as an inorganic anion, with greater adsorp-
tion at lower pH values 10).
The adsorption of nitrophenols onto alumina pretreated with various acids is
shown in Refs. 124 and 125, unfortunately without giving the final pH values.
First, one point may be mentioned, which was not discussed before: If the intention
is to prepare a surface with an ionic adsorption layer, as for catalyst preparation,
then adsorption from the gas phase may be an alternative. In this way organic ions
may be adsorbed onto an alumina surface dehydrated at a higher temperature.
Infrared spectra reveal that a surface aluminum salt was formed under these con-
ditions Halogenated surfaces can be formed in the same way
Complexes can be adsorbed in this way too, though normally neutral complexes
are used here because of the better volatility.
Although it is not the intent of this chapter to thoroughly discuss model calcula-
tions, they also have to be regarded to a certain extent. In recent years a number of
models have been developed that are broadly based on the surface
model of et which, in agreement with the spectroscopic results,
shows a chemical interaction at the surface. Some recent extensions of this model
try to close the gap between adsorption, precipitation, and coprecipitation with
different assumptions - surface precipitation, surface coprecipitation, the forma-
tion of a continuously variable surface layer or competitive precipitation
of pure or mixed oxides A revised physical
model based on very different assumptions was published recently by
and All models give fairly good fits to experimental values.
Robertson and Lecky demonstrate that equilibrium constants are "condi-
tional," they vary the experimental conditions, which means either that
not all experimental data were sufficiently considered or they were not correctly
implanted into the model.
The fact that different models lead to a sufficient fit of quantitative experimental
data shows that additional information be used for discrimination.
Spectroscopic data reveal the formation of chemical bonds between adsorbed
ions and the surface, which led to certain convergence to surface complex models,
and extended to considering adsorption in the case of high ion concen-
tration the solution.
A serious problem is the fact that for most models the number of model para-
meters is great and the information obtained the experiments comparatively
poor. Many parameters have to be taken from other sources or have to be esti-
mated (or even disregarded, because it is supposed that they have small influence
an estimation). Activity coefficients in the liquid phase are usually omitted by using
a "sufficient" concentration of background electrolyte. The surface activity coeffi-
cients are usually not mentioned at all. most only one type of adsorption
A. and A. (1986).
and
K. Chattoraj, Surf. Tech.
Schulthess and Sparks, Soil (1988).
S. Howard, and Han, (1993).
Colloids Surf. A
C. Ye, Solid State
Shiao, and E. Meyer, (1981).
D. Jackson, and J. K. Soil Am, J.
(1976).
Chang, Lo, Yasunaga, Colloid
Interface (1997).
R. and Evans,
Parvulescu, and E. Rev.
35: (1990).
Chang, Lee, and Lo., (1993).
and Czech. Phys. B
Cozar, and Rev. Phys. (1982).
Y. Zeng, and React. Kinet. Lett, (1991).
Rudin and Motschi, Colloid Interface
Hachiya. Sasaki, Y. Saruta, Mikami, and Phys.
(1984).
Hachiya, Sasaki. Ikeda, Mikami, and
and
L. Villa, React, Catal. Lett. (1975).
de Bokx, A. Wassenberg, and W. Catal. (1987).
and A. Lycourghiotis, (1994).
and A. Lycourghiotis, Faraday Trans. (1993).
St. R. G. Brown, and A. Parks, Colloid Interface
(1997).
A. Scheidegger, and Sparks, Colloid
(1997).
E. Brown, A. and
(1989).
A. and A. Schwarz, Catal.
Stumm and Rev. Geophys. (1990).
He, , and Martens, Soil
Am.
A. C. Colombo, and Buondonno, Soil Sci. Ain. (1991).
Petkovic, and Technol.
I.Introduction
II .
Protein Structure
III .
Polystyrene Latex Particles
IV. Interaction of Proteins with Polystyrene Latex
A. Adsorption
B. Conformation of adsorbed protein
C. Stability of adsorbed protein layers
D. Controlling protein adsorption
V. Conclusions
References
Proteins are linear chains of amino acids in which individual residues are linked by
bonds between the carboxyl group of one acid and the amino group
of the next. This is the primary sequence of the protein. The side chains of some of
the amino acids undergo enzyme-catalyzed modifications such as glycosylation and
phosphorylation, which serve to alter the functionality of the protein. Folding of
the polypeptide chain as a result of hydrogen bonding between C 0 and N-
groups of the amino acids gives rise to a-helices and and sheets, with
structured regions being linked by flexible loops. This is the secondary structure.
Other types of secondary structure have also been identified, but in globular pro-
teins these forms predominate. Compared to bond energies ( > 200
these hydrogen bonds are relatively weak (around
The extent of ordered structure in different protein species is very variable.
Some proteins, and in the context of this chapter the caseins are a prime example,
adopt a relatively open, expanded-coil conformation. In contrast, globular pro-
teins are almost spherical, with a packing density at the interior of the molecule
being equal to that found in crystals and close to the theoretical value for
packed spheres. This high packing density means that short-range interactions
such as hydrogen bonding and van der forces are important in stabilizing
the secondary structure. In general, hydrophobic amino acid residues tend to be
buried in the interior of the protein, whereas hydrophilic charged groups are
located on the outside of the protein in contact with the aqueous phase.
However, polar and amino acids are fairly equally distributed at the sur-
face of globular proteins.
An additional type of covalent bond, the disulfide bridge, can form between
individual cysteine residues and is also important in determining the stability of
many proteins. Proteins with disulfide bridges usually contain at least two, proteins
containing single disulfide bridges being relatively rare. Both inter- and intramole-
cular disulfide bridges are found; the addition of thiol-containing reagents such as
Since polystyrene latex particles are generally hard spheres, protein adsorption to
them is solely a surface phenomenon. Methods for the production of polystyrene
latex particles have been reviewed by et
The diameter of polystyrene particles can range from less than 0.1 to more
than 1000 By varying the polymerization conditions, it is possible to produce
polystyrene latex particles in a broad range of sizes and degrees of monodispersity.
Larger beads are usually produced by suspension polymerization. In this process,
agitation breaks the monomer phase into individual droplets suspended within the
dispersion phase, which also contains a surfactant to stabilize the droplets. An
initiator is also incorporated into the phase, and by varying both the
rate of agitation and the temperature, polymerization is induced. The polystyrene
particles produced by this technique typically have mean diameters between 50 and
1000 with a relatively broad distribution. Emulsion polymerization, where
surfactant-stabilized monomer micelles are polymerized by initiator molecules dis-
solved in the dispersion phase, is used to produce extremely uniform, very
spherical particles with diameters up to 0.1 In phase inversion polymerization,
the monomer, initiator, and stabilizer are all dissolved in an inert solvent. When the
growing polymer chains reach a certain critical size, phase inversion occurs, and the
chains aggregate to form stable monodisperse particles. The size range of particles
produced by this method is between those of particles produced by the other tech-
niques.
The polarity and density of charged groups on the surface of latex particles
can be varied by altering the chemical composition of the monomer mix.
Alternatively, it is also possible to change the surface by allowing acrylate mono-
mers to adsorb to preformed particles and subsequently polymerizing these in
The surface charge density of beads can be determined by measuring their
average zeta potentials and electrophoretic mobilities or by conductometric
Size, and smoothness are usually determined by microscopy or
photon correlation spectroscopy.
The various stages in the interaction of a protein molecule with a surface are shown
schematically in Fig. 1. Initially, protein molecules in the bulk phase approach the
latex surface. Adsorption may then occur, and in some instances this is followed by
rearrangement of the protein structure; both of these steps are controlled through a
combination of ionic and hydrophobic interactions and hydrogen bonding. Finally,
in some systems the structurally modified protein molecules be released back
into the bulk phase. These stages are now considered in more detail.
7.0
Urea
dephosphorylated
-Casein Urea
K-Casein,
nonglycosylated
K-Casein,
glycosylated
K-Casein
7.0 2-ME
Phosvitin
Urea
thicknesses determined by PCS except for those marked*. where SAXS was
than that measured by solution depletion techniques, the difference being explained
by the authors as being due to loosely associated protein that is not detected by the
method. Adsorption of to latex particles at p H 7.0 determined using
showed that adsorption times were less than 10 min and the maximum
value of was 4.75
2. Influence of pH and Charge on Protein Adsorption
Proteins are electrolytes with the side chains of their acid,
acid, tyrosine, histidine, and arginine residues being
involved in acid-base equilibria. The net charge on a protein is therefore dependent
on pH. At the point (pi), the net charge on the molecule is zero, and at
pH values below and above their proteins are positively and negatively charged,
respectively. The interaction of a protein with a surface will therefore vary with
and, in general, the maximum value of is attained at or near the of the protein.
Obviously, electrostatic interactions between the surface and the protein and
between adjacent adsorbed protein molecules are therefore important.
0-Casein,
dephosphorylated
nonglycosylated
K-Casein,
glycosylated
Ase
a-Lactalbumin
Lysozyme
Fibrinonectin
Fibrinogen
Immunoglobulin
Lipolase
Savinase
Bovine serum
glycosylation of this protein increases its net negative charge by one unit for each
sialic residue attached to the polypeptide chain and, as expected, also increases
the thickness of adsorbed protein layer on a negatively charged surface
(Table 1).
Another way to investigate electrostatic interactions is to change the surface
characteristics. al. prepared a series of latices of the same
particle size but different surface structures and measured the adsorption of human
serum albumin (HSA), cycochrome, myoglobin, and horseradish
oxidase. In general, hydrophilic surfaces had little affinity for proteins. Generally,
increasing the hydrophobicity or the surface charge increased the affinity of pro-
teins for the surface. Electrostatic interactions were found to be more important in
determining the adsorption of polar and relatively small protein molecules. If
adsorption of a protein onto a surface was largely determined by electrostatic
interactions, it would be expected that these should be maximal in the region
where the protein and the surface are oppositely charged. However, measurements
of the adsorption of HSA and y-globulin onto carboxylated surfaces showed that
this was not the case This may be due either to conformational changes in the
protein structure induced by adsorption to the surface, which weakens electrostatic
interactions between the protein and the surface, or to deionization of groups on
the surface of the protein, or to both.
In addition to the polarity of the surface, the density of charged groups on the
surface also affects protein adsorption At carboxyl group densities less than
2
1 n m , adsorption of bovine serum albumin to latex was both dominated by
hydrophobic interactions and relatively sensitive to pH. At densities above two
carboxyl groups per square nanometer, hydrogen bonding predominated irrespec-
tive of pH. The surface coverage of lysozyme and (a-lactalbumin
from which the calcium ion bound at the high affinity binding site has been
removed and which is therefore less stable than the native protein) on negatively
charged latices has also been shown to be dependent on the density of charged
groups on the surface of the support. At pH 7, the adsorption of lysozyme, which is
positively charged at this pH, is higher on a more negatively charged latex whereas
that of apolactalbumin, which is negatively charged, is lower, demonstrating the
importance of electrostatic interactions to protein adsorption Gallisteo and
Norde also showed that the maximum value of was pH-dependent, being
highest at a pH value just below the p i of the individual protein and decreasing at
pH values to either side of the As expected, the influence of the density of
negative charge on the surface loading of these two proteins was also pH-depen-
dent, with loading on the more highly negatively charged surface being greater at
pH values below the of each protein.
As part of a study aimed at the rational development of the
influence of surface characteristics on the interaction of fragments
obtained by peptic digest of rabbit polyclonal has been investigated
The latex particles differed in their hydrophilic character and in the nature and
density of their surface ionic groups. The most important force governing the
adsorption of the protein onto the surface was shown to be hydrophobic inter-
actions, although in low ionic strength buffers electrostatic interactions were also
found to play a part. On anionic latices, maximum values were obtained at p H
2. Globular Proteins
Due to their more highly folded structure, globular proteins can potentially
undergo more extreme conformational changes when adsorbed to surfaces, and
interpretation of these structural changes is consequently more At the
simplest level, if the molecular dimensions of a protein are known, then, purely on
the basis of the measurements of the surface coverage, information on
the orientation of the molecules on the surface can be obtained (Fig. 4). and
their fragments apparently adsorb onto hydrophobic surfaces in an end-on
orientation but provided the proteins have a high charge density opposite that of
the surface, the orientation on negatively charged hydrophilic latex switches to
on in order to maximize electrostatic interactions
More detailed information on the conformation of adsorbed proteins has been
obtained by using a variety of physical techniques. Although fluorescently labeled
HSA and have been shown to undergo conformational
changes consistent with the loss of internal fluidity when adsorbed onto
latexes in general the application of spectroscopic techniques to
the analysis of protein conformation at solid interfaces is complicated by the scat-
tering and adsorption of light caused by the particles themselves. This can be
reduced either by using very small particles that adsorb and/or scatter light to a
much smaller degree or by using particles that display negligible adsorption and
scattering. Using the latter approach, circular dichroism (CD) and fluorescence
anisotropy were used to detect changes in the conformation of the proteolytic
enzyme 309 resulting from its adsorption to a novel hydrophobic
fluoroalkoxy fluorocarbon latex At = the fluorescence of the
adsorbed protein was quenched and an increase in the a-helix content was detected
by This finding was in contrast to the decrease in a-helix content that is usually
detected when proteins adsorb onto hydrophilic silica particles. At =
which represents full surface coverage, the fluorescence and a-helix content of the
adsorbed protein were essentially the as in solution, indicating that in this
instance the extent of conformational changes in the structure of adsorbed proteins
is dependent on the density of packing of the molecules on the surface.
In addition to spectroscopic techniques, microcalorimetry, which measures tran-
sition enthalpies (AH) and transition temperatures (To.) as proteins undergo coop-
erative unfolding, has also been applied to the study of conformational changes
resulting from protein adsorption. Since PS latex particles show no thermal transi-
tions over the temperature range of interest for protein structural analysis, their
presence does not interfere with this technique. HSA, a three-domain protein whose
HIGH
Temperature
Proteins adsorbed onto polystyrene latex particles have proved to be good models
for the interfacial behavior of proteins in processed foods. The behavior of caseins
on latex particles appears to be very similar to that at the interfaces of
emulsions, at the interface of foams, and at the surface of casein micelles.
Core funding for the Hannah Research Institute is provided by The Scottish Office
Agriculture, Environment and Fisheries Department.
A. A. and Ed.
(1994).
Dalgleish and J. Leaver, in Food Polymers,Gels and Colloids Dickinson,
Royal Society of Chemistry, Cambridge, UK, 1991. pp. 1
G. Dalgleish, Colloids Surf. B: Biointerfaces (1993).
Dalgleish, Colloids Surf. (1990).
V. Brooksbank, Davidson, S. Horne, and Leaver. Chem.
Faraday Trans. (1993).
A. Mackie, and A. North, Faraday Trans.
(1991).
J. Leaver, Brooksbank, and Colloid Interface
(1994).
A. Mackie. and Dann, in Food Polymers,Gels and Colloids
Dickinson, Royal Society of Chemistry, Cambridge, UK, 1991, pp.
A, Kondo and K. Colloid Interface Sci. (1992).
Norde, Cells Mater. (1995).
Bale. and Sutton, Colloid Interface Sci.
(1988).
Ortega Galvez and
(1996).
Duinhoven, Poort, van der Voert, Agterof, Norde, and
Lyklema, Colloid Interface (1995).
Duinhoven, Poort, van der Voert, Agterof, W Norde, and
J. Colloid Interface Sci. (1995).
I. Introduction
II. Adsorption onto Porous Materials
I I I . Improving Dissolution
IV. Effect of Adsorption on Chemical Stability
V. Fluorescence Spectroscopic Studies of Interaction
VI. Therapeutic Use of Porous Materials
VII. Conclusions
References
When poorly water soluble drugs are administered in solid dosage forms, the dis-
solution rate is often considered a rate-limiting factor in the adsorption process.
The rate of dissolution of these drugs depends on the effective surface area and the
energy state within the drug crystals The amorphous form of the drug can
exhibit a much higher dissolution rate and better bioavailability than the crystalline
form. Various techniques have been used to reduce the crystallinity or to obtain the
amorphous form of a crystalline drug The formation of an amorphous solid is
usually achieved by freeze drying, spray drying, solid dispersion, condensation from
the vapor, and milling. The rapid cooling of indomethacin liquid causes the glass
formation of the drug et [4] reported that amorphous
prepared by pulverization and lyophilization showed different dissolution behavior.
According to Yamaguchi the amorphous 16-membered macrolide com-
pound prepared by different conditions of spray drying showed different degrees of
solubility and stability.
The changes in the X-ray diffraction patterns of mixtures of 30% flufenamic acid
(FFA) and 70% of either florite R or florite S after storage at 0% RH and are
shown in Fig. 1 In the mixture of F F A and florite the intensities of X-ray
diffraction peaks due to FFA crystals gradually decreased with period of storage,
and a halo pattern was observed after 5 weeks. In contrast, these X-ray diffraction
peaks were still observed in the FFA-florite mixture even after storage for 8
40 80 40 120
Temperature Temperature
is of the form
Although this calculation was rather too simple, a considerably large value com-
parable to the total amount added was obtained, which could explain the rapid
disappearance of crystals of benzoic acid.
Drug molecules that penetrate into the capillaries have higher sticking probabil-
ities (adsorption rates) than those colliding with flat surfaces because they undergo
a great number of collisions with the capillary walls Thus, these
easily lose their kinetic energy during the collisions and rapidly become fixed on the
walls, which provide the sink condition at the surface of the porous powders.
T o verify that these phenomena take place via the gaseous phase, Konno et
studied the rate of amorphization using crystalline substances with various
vapor pressures. Benzoic acid and its derivatives were used as the crystalline sub-
stances, and activated carbon as the adsorbent. After each crystalline substance was
mixed with activated carbon in a porcelain mortar, the changes in the X-ray dif-
fraction patterns during storage were examined, and the time required for the
crystal peaks to disappear was determined. The relationship between the vapor
pressure of the crystalline substances and the time required to become amorphous,
plotted on a logarithmic scale, is shown in Fig. 4. At atmospheric pressure and
benzoic acid, which has the highest vapor pressure of the substances used,
became amorphous within only 1 whereas p-hydroxybenzoic acid, whose vapor
, , ,
10
Time
Plots of vapor pressure of the crystalline substances versus the time required to
become amorphous in mixtures with activated carbon. At and atmospheric
pressure; at and reduced pressure. 1, Benzoic acid; 2, o-hydroxybenzoic
acid; 3, acid; 4. p-hydroxybenzoic acid; 5, o-chlorobenzoic acid; 6,
acid; 7, p-chlorobenzoic acid. (From Ref. 13.)
dency was recognized at reduced pressure at The results suggest that the
phenomenon of the gradual loss at room temperature of the crystalline
properties of some organic pharmaceuticals mixed with adsorbents takes place
via the gaseous phase and is induced by the vapor pressure of the pharmaceuticals
themselves
Activated carbon, the residue from the destructive distillation of various organic
materials, was found to be capable of forming an amorphous solid. After storage of
a mixture of aspirin and activated carbon at for 7 days, the aspirin changed
from crystalline to amorphous. In the condition of reduced pressure, the
zation, the transfer of organic molecules to adsorbent surface took place more
rapidly as the mean free path of the organic molecules was much longer than at
pressure
In the presence of porous and nonporous silicon dioxides, organic crystallites
were converted into the amorphous state. Kim et al. reported the kinetic of the
transformation of the crystalline drug in the presence of various types of silicas
using DSC.
The relative enthalpy change obtained with the 1 : 1 drug-silica mixtures after
storage at is shown in Fig. 5. These thermal data clearly indicate that the rate
of crystalline transformation into the amorphous state was markedly dependent on
the types of silica employed. The mixtures with a mesoporous silica, such as silica
Time
.7 Dissolution profiles of FFA from mixtures with MAS into various dissolution
media at 20% FFA, fresh; 20% stored for 5 h at a reduced
pressure at A, pH 1.2 solution (JP 1st fluid); 3.0 citrate buffer solution;
pH 4.0 acetate buffer solution; pH 5.0 acetate buffer solution. (From Ref. 16.)
T i m e (min)
The need to ensure product quality and hence stability of the drug concerned in the
product is obvious. Decomposition leads to a decrease in drug potency, and the
decomposition products may be unpleasant or toxic. Stability is therefore an
important issue with pharmaceutical products.
The hydrolytic decomposition of aspirin has been widely studied in aqueous
media. It has been shown that the rate of hydrolysis of solid aspirin is proportional
to the relative humidity of the atmosphere The effect of additives on this
hydrolysis has been investigated, and the presence of additives has been found to
affect the decomposition rate of aspirin.
The kinetics of aspirin hydrolysis were studied in mixtures containing 5%
aspirin. The pore diameters of varied from 75 to Plots of the
order rate constants at 50° under the two different relative humidities versus the
pore diameter of CPG are shown in Fig. 9 In the CPG mixtures where the pore
diameter was less than the rate constants of aspirin decomposition were
independent of CPG pore diameters in the cases of both RH and 79% RH.
However, the rate constant in the mixtures with CPG of large pore size decreased
with increasing CPG pore diameter. The pore diameter profiles of the rate constant
indicate that the stability of aspirin was closely related to changes in the crystal-
P o r e diameter
Effect of pore diameter (log scale) on the first-order rate constant for the
decomposition of aspirin at Aspirin concentration is 5%. 0% RH; 79%
(From Ref. 20.)
Relative
Source: Ref. 24
in good agreement, and they can be considered the critical amounts of water when
aspirin hydrolysis takes place in the heated mixtures.
Edwards investigated the hydrolysis of aspirin in aqueous solution at 17°
and demonstrated that the aspirin is hydrolyzed by general acid-base catalysis and
water molecules for ionic and nonionic aspirin, comprising six simultaneous reac-
tions involving O H , and for ionic and nonionic aspirin. The intrinsic
hydrolysis rate constant in the heated mixture was comparable with the hydrolysis
rate constants of the two-element reactions of nonionic aspirin and or ionic
aspirin and in aqueous solution. As aspirin molecules would be adsorbed onto
the pore surface of PCC in the molecular state, as a possible mechanism of aspirin
hydrolysis in the mixture with PCC it was suggested that the aspirin is dispersed
monomolecularly in the heated mixture and reacts with water molecules rather than
by acid-base catalysis.
. 12
Change in fluorescence spectrum of 1% pyrene-PCC by storage at
= Left: Changes within 1 day. Right: Changes within 2 weeks.
(From Ref.
control and bleomycin solution groups. In contrast, the size of the tumor signifi-
cantly decreased in the bleomycin-small silica group 0.01, p 0.05). Lymph
node metastasis occurred in ail the animals of the control group, but the incidence
of such metastasis was 64% 1) in the bleomycin solution group and only 9%
11) in the bleomycin-small silica group. A new drug delivery system using silica
particles might be useful for the control of tumor and lymph node metastasis.
When mixed with porous materials, organic medicinal molecules changed their
arrangement in the crystal lattice from the crystalline state to a disordered phase,
accompanying changes in pharmaceutical properties. This transformation pro-
ceeded through the sublimation of organic from the crystals and adsorp-
tion onto the porous materials via the gaseous phase. The adsorption equilibrium
depended on the temperature, water content, and surface properties of the porous
materials. The resulting solid phase consisted of continuous disordered states invol-
ving a liquid-like molecular arrangement near the surface, an amorphous portion
including a vitreous state, and imperfect crystals where the crystal growth was
highly inhibited in the pore. Since solid dosage forms such as tablets, capsules,
and pills are muiticomponent systems containing the active ingredient, diluents,
binders, lubricants. and other pharmaceutical additives, it is important to recognize
that adsorption phenomena have a great influence on the effectiveness, stability,
and safety of dosage forms.