ADSORPTION
AND
SURFACE
ADSORPTION N.
A.
ENERGY ON
THE
(I):
THE
EFFECT
OF
y-PLOT*
GJOSTEINt
Based on the idea that an adsorbing gas atom can interact with any of three different types of surface sites-surface, single ledge, and double ledge sites-each being characterized by discretely different adsorption free energies-a phenomenological theory is developed, governing the shape alterations of the y-plot due to adsorption from the gas phase. When single and double ledges are preferred sites for adsorption, the theory predicts that an appreciable increase in the curvature of y-plot can occur as the pressure of the system approaches the characteristic pressures of these sites. At some higher pressure, depending on the relative adsorption energies of single and double ledge sites, a marked decrease in the slope of the y-plot at the origin of the cusp will occur. When the surface sites are preferred for adsorption the curvature of the y-plot decreases, and the slope at the origin of the cusp increases as the system pressure approaches the characteristic pressure of the surface site. In addition to, and simultaneously with these changes, each point of the y-plot will shrink towards the origin of the y-plot. ADSORPTION
ET ENERGIE L’ADSORPTION
DE SURFACE (1): 1’INFLUENCE SUR LE DIAGRAMME y
DE
Elle est basee sur l’idee qu’un atome gazeux peut s’adsorber en trois genres d’emplacements, caract&i&s par des energies libres d’adsorption differentes, placement en surface proprement dit, contre un gradin peu marque (“single ledge”) au creux dun gradin accentue(“double ledge”). Les auteurs developpent une theorie qui regit les modifications de la forme du diagramme de l’energie de surface y, dues 21 I’adsorption dune phase gazeuse. Quand l’adsorption a lieu preferentiellement contre des gradins, cette theorie prevoit une augmentation sensible de la courbure du diagramme y quand la pression du systeme s’approche des pressions Pour une certaine pression plus elevee, en relation avec les caracteristiques de ces emplacements; energies relatives dadsorption, se produit une diminution marquee de la pente du diagramme y vers son origine .
Quand l’adsorption a plut6t lieu en surface proprement dite la courbure du diagramme y diminute, et la pente a l’origine augmente quand la pression du systeme d’approche de la pression caracteristique du site d’adsorption. En m&me temps, les points du diagramme y se rassemblent aux environs de l’origine. ADSORPTION
UND OBERFLACHENENERGIE (I): EINFLUSS ADSORPTION AUF DIE 1/-FLACHE
DER
Ein adsorbiertes Gasatom kann mit drei verschiedenen Oberlliichenpliitzen wechselwirken: mit Pliitzen auf der F&he, an einer einfachen und an einer doppelten Stufe; die Freie Energie der Adsorption ist fur die drei Pliitze verschieden. Diese Vorstellung wird einer phanomenologischen Theorie zugrunde gelegt, die die Anderungen der y-F&he infolge Adsorption aus der Gasphase behandelt. Wenn einfache und doppelte Kanten bei der Adsorption bevorzugt werden, kann nach der Theorie eine betrachtliche Zunahme der Krtimmung der y-Fliiche auftreten, sobald der Druck des Systems sich den charakteristisohen Drlicken dieser Oberfliichenpliitze ntihert. Bei einem hijheren Druck, der von den relativen Adsorptionsenergien der einfachen und doppelten Stufen abhangt, tritt eine ausgepriigte Abnahme der Steigung der y-Flache an der Spitze auf. Wenn die bevorzugten Pliitze der Adsorption die auf der Fl5che sind, nimmt die Krtimmung der y-Fliiche ab, und Steigung an der Spitze nimmt jedoch zu, wenn der Druck des Systems sich dem oharakteristischen Druck des Flachenplatzes nahert. Gleichzeitig mit diesen Veranderungen tritt eine zusiitzliche Schrumpfung der gesamten y-Fliiche auf.
INTRODUCTION
Recently, determining
cusps at the (100) and (111) orientations.
methodso-4)
have
been
developed
the variation
of the surface energy
In nickel
for
some minor
(y)
notably (llO), (210) and (410), but the results here are not as well established.(i)
of metals with crystallographic orientation of the surface. When interpreted in terms of the polar y-plot,(5*6) these results indicate that the shape of
In carrying a strong
cusps may occur
out experiments
likelihood
that
at other orientations,
of this type there is
impurity
atoms,
adsorbed
the y-plot for nickel and copper is similar to that predicted by Herring(b) several years ago. For these
from the gas phase or from the crystal interior, will be
metals, curved
understand the alterations in shape of the y-plot which will occur as a result of gaseous adsorption. Lacmann and Stranski”) have given a treatment for the change in surface energy of low index planes with
the y-plot appears surface except for
to be a continuously major inward-pointing
* Received July 19, 1962. t Metallurgy Department, Scientific Motor Company, Dearborn, Michigan. ACTA
METALLURGICA,
VOL.
Laboratory,
11, AUGUST
1963
Ford
present
on the surface,
gaseous adsorption, 957
and thus it is important
to
but there has been no treatment
ACTA
958
METALLURGICA,
of the effect of adsorption on the shape of the y-plot in the vicinity of a cusp. The present work is directed towards solving this latter problem. Utilizing a model of surface struoture, which has evolved from the surface energy work described above, a phenomenologi~al theory is developed, which predicts the change in shape of a cusp in the y-plot with the extent of adsorption. From the results of this analysis, an estimate will be made of the extent to which the present determinations of y-plot have been influenced by adsorption. In an accompanying paper,(*) it will be shown that dependence of surface-energy-induced thermal faceting on gaseous impurity adsorption oan be explained on the basis of this analysis. SURFACE
STRUCTURE
Low index planes Burton, Cabrera and Frank(s) and &Iullins(lo) have treated the problem of atomic roughening of an atomi~a~y smooth low index plane that occurs with increasing temperatures. At elevated temperatures, a small concentration (approximately 5 X IOP5, at temperatures near the melting point in copper) of surface adatoms and vacancies are stable due to the lowering of free energy of the surface that occurs from the configurational entropy of the roughened surface. Both treatments conclude, however, that for most metals no large soale disordering (surface melting) of low index planes will occur even at temperatures near their melting points. For the purpose of this treatment, they will be considered as atomically smooth. Vicinal prunes In discussing the origin of inward-dinting cusp in the y-plot, Herring(5@) suggested a structure for the vieinal planes of a low index plane. According to this picture, the vicinal planes should differ in structure from the low index plane only by the presence of certain density of monatomic ledges of height S and average separation
&G
Wu3
where M is the small angle of deviation from the low index plane. The surface energy y of a vicinal plane for small tc is given by
where y0 = surface energy of low index plane (ergslcmz) yI = excess free energy of surface per unit length of ledge (ergs/cm).
VOL.
II,
1963
(III)
PLANE
FIG. 1. Torque terms for (111) cusp of copper, 0 -0FHC e-high purity copper in copper in graphite boat; graphite boat; a-OFHC copper with no graphite in system (3). Taken from Robertson and Shewmon, to be published in Trans. Amer. Inst. Min. (MetaEZ.) Engrs. 224 (1962).
Equation (1) predicts a cusp wit8ha parabolic shape, having
a constant torque
term,
ay 75 acr = x -
yOa.
Experimental determinations of torque terms for nickel(l) and coppeG show that the torque term decreases approximately linearly* with angular deviation from the low index plane, Fig. 1. The magnitude of the decrease, however, is generally much smaller than would be expected from the coefficient yO. This indicates that another second order term in u should appear in equation (1). Robertson and Shewmon(3) have shown how this term can arise, if the steps are randomly distributed, from the probability of forming double ledges-that is, ledges of height 2s.t Their equation may be written:
y=yo+
Fu- [~ + I&]u2
(2)
* In the data of Robertson and Shewmon,‘3) the scatter is too large to conclude that there is any deviation from linearity. Mykura”) indicates some curvature in his plots, but this will be ignored in the present treatment. t Other models also can give rise to a second order term. For example, if ya varies inversely with the single ledge spacing, to a first approx~tion, a relationship similar to equation (2) can be derived. Alternatively, it can be postulated that the ledges are in thermal equilibrium and that two single ledges combine with the formation energy yrbe, where _5? is a length per atom along a ledge. In this case, the COefficient of cxzin equation (2) will be yr/S*M exp [+y&‘/kT]. The treatment given in this paper wa8 carried out only for the statistical model of Robertson and Shewmon because it appears to be the simplest model consistent with the experimental facts. It would not be difficult to modify the present treatment, to t&e into account these other ideas.‘aO’ 3 In equation (2), the density of single ledges a/S has not been modified to take into account the presence of double ledges because the density of double ledges is small in comparison to x/X, especially for small a.
GJOSTEIN:
ADSORPTION
AND
SURFACE
treatment
where is YI = 2YL - YD (positive or negative)
the
interaction
energy
between ledges, yD is the
energy per unit length of double ledge, and M
ENERGY
shows
evident
two
quantities
mentally,
which
to define for later use
may
be
determined
experi-
account
the
&_L
(a2r1=
yoO 8x2
2y,
to this model
for believing
that
+k
(4)
S2MYCW
where yoOrefers to y0 for perfectly
the
heterogeneous
there appears a linear
to be no
extrapolation
curvature
of since
any higher order terms in equation
would give rise to appreciable
(2),
of the torque
term plot at large angles rather than small. For both nickel and copper, however, there are no data at small
angles
information
to verify
a minimum
of a facet.
except
for some
from the equations
These equations
will give
value for the torque terms at the origin
of the cusp.
Torque
these equations taken
this point,
that may be obtained
of equilibrium
from
terms for nickel obtained
agree well with extrapolated
the twin
boundary
method,
from value
according
to Mykura,(l)
but this does not seem to be the case
for
Gjostein(l2)
37 1
copper. 1
MODEL
discussion
it is obvious
estimated
value for -aa to be about Y00 [ 0
structure-that
the
minimum
0.4 for the
(100)
character
is, the defect
of the
sites,
an adsorption
simplicity, such
the
defects
impurity
atoms
considered different double
site
kinks,
and sites arising
sites.
so too, Instead,
adsorption is
shown
although
AG,“,
three
surface,
AG,‘,
in
Fig.
point
3,
with
evident
of these
quantities
bond
be in a state
facets on which measurements
it is necessary
further experimental
is needed.
evidence
an adsorption An
> AG,“,
AG,o,
The relative magnibe determined
from
knowledge
of the nature
of the
in which AG, to
a phenomeno-
can
that an adsorbing
Shewmon. The former value was based on only a few observations due to the difficulty of obtaining Thus,
and
in Fig. 2.
and AGF may be
atom
theory of surface bonding
can be made.
ledge
> AG,’
From
AG,’
lated
a serious discrepancy,
single
respectively.
AG,’
later.
of view, both
either a theoretical
AG,O,
will be
energetically
this order is not essential to the treatment
either greater or less than AG,‘. tude
of
bulk
energy diagram, which depicts this model,
as will become logical
AG,,
of
dislocations
the gas atom
the ledge or double ledge, or from
before this can be considered
adatoms,
ledge sites, as shown schematically
energy,
sake
energy
will be eliminated
with
of sites:
for
from
they
a
that
could serve
surface
To each type of site will be attributed free
In
characterized
Since,
to
vacancies,
as interacting types
being
free energy.
as
as adsorption
each
contribution
have been ignored,
electron
This is much larger than the extrapoof 0.067 obtained by Robertson and
surface.
all the defects that have been mentioned
plane in copper. value
that
must take into
general sense, it would be expected
as adsorption by
clean surface.
term curve to cc = 0 is not valid,
presumably
between
of a surface.
completely
the torque
will
The purpose of this section is to derive a set of adsorption isotherms which reflect (3)
reason
will be made
of gaseous adsorption surface
distance
As will become
sites in this treatment.
the previous
any treatment nature
namely
According
interkink
ADSORPTION
to the mean ledge
959
distances.
later no distinction
kinked and non-kinked
From
direction. At this point, it is convenient
the
be about five interatomic
is the number of possible ledge sites per cm of length along a line normal
that
(I)
with The
presently does not seem to the relative
can be predicted, rely
makes
experiment. magnitude
on experimental
meagre as it may be, regarding
evidence,
this point.
Ledge structure At elevated temperatures, will have a constant
the ledges at equilibrium
mean direction,
but they will
not be atomically straight.@) In analogy to the structure of a low index plane, a ledge will contain a certain density of kinks at temperatures above absolute zero. In contrast to the structure of a low index plane, however, the defect concentration is much larger in the case of the ledge. For example, for copper near it#s melting point, the Burton
et uZ.(~)
of
and therefore
FIG. 2. Schematic representation of interaction of an atom in gas phase with the surface (s), single ledge (I), and double ledge (d) sites on the surface.
ACTA
960
METALLURGICA,
VOL.
11,
1963
different sites on the surface, and are, of course, dependent on the first three reactions. Due to similarity of the first three reactions, the same formalism applies to all of them. It is necessary, therefore, only to derive the adsorption isotherm for the ith site. Consider the Gibbs free energy change SG whieh occurs as a result of transferring an infinitesimal number of atoms, Etnp,from the gas phase to occupy a certain number, 6n,, of the ith site. This free energy change can be written in terms of the chemical potentials, ,u, of the three species participating in reaction (5a), and is given by:
Gas Phoss
‘P FIG. 3. Adsorption energy dirtgram.
That atoms are more strongly adsorbed at ledges than on the surface receives support from the flashfilament-desorption studies of Hiekmott and Ehrlichol) and from low energy eleetron diffr~otion studies,(12) but these investigations do not provide information about the relative magnitudes of AC,“, AC,’ and AC,“. As will be demonstrated later,(*) an analysis of faceting phenomena can give some idea as to the magnitude of AG,‘. The adso~tion isotherm will depend on the mechanisms of adsorption and desorption as well as the type of adsorption layer that is formed, i.e. whether it is mobile or immobile. Since the main features of the present analysis can be illustrated quite well by considering the simplest kind of adsorption process, namely, the localized (immob~e) adsorption of single atoms, this will be considered first, while the complexity arising from dissociation of diatomic molecules will be treated later. It is necessary to develop this latter case also, since there is a strong likelihood that diatomic moleoules do dissociate upon adsorption, although this has been demonstra~d convincingly only for nitrogen on tungsten.ol) As will become evident, the type of adsorption isotherm does not alter the kind of changes that adsorption produces in the y-plot, but only their magnitude. The model just described gives rise to the following equilibria: Gi-X?r?:Ge G+LZG, G+DtiG= Gs*GL
(5a) (5b) (50) (5d)
G,ZG, GLZZGD
OW @f)
G symbolizes the gas atom; S, L, and D represent the surface, single ledge and double ledge sites; the last three reactions represent the equilibria among
where N is Avogadros number. Since the adsorption of a single gas atom destroys one of the ith sites and creates one g-i complex, it follows that
tin,
lingi = 6n, =
At equilibrium, 6G = 0, and equation (6) becomes IQ@6 = PI + I%
(7)
Making use of the thermodynamic relationship pi= & + RT1 n ai,in which pi0 is the chemical potential of the ith species in the standard state and ai is the activity of the ith species, it is possible to rewrite equation (7) in the form a 25
= exp
a@ai
1Pi0+;;-
aq _exp[-;;q
(*)
where ,QO - ,ugO- &’ z AGiQ is the adsorption free energy for the ith site, when the species are in their standard states. The activities aj have the following definitions in terms of surface concentration Ni (atoms/cm2 or atoms/cm), the saturation surface concentration Ni, (atomslom2 or atoms/em), the activity coefficients fsi, fi and f,, the pressure in the gas phase P and the pressure in the standard state PO.
(W
@2=f,
P p
() 0
(90)
GJOSTEIN:
ADSORPTION
AND
From these definitions, equation (8) can be placed in the form*
which defines a characteristic pressure Pi for the ith site given by (II)
SURFACE
ENERGY
expression for y0 is ye - yw =
-kTN,,
’ - p’ps s 0 (1 + P/P,,
=
-kTN,,
In (1 _t P/P,)
dlnp
yzo = --kTfi,,ln
Ideal behavior
(1 + P/PJ + kTflN,9, ln (I + P/P,)
The quantity f,&‘Ji the relation
may be considered, through (12)
as & measure of the lateral interaction energy AGir for this particular site. In general, AGil is a complex quantity and will be some function of the occupation probability pi = NJN,,. The most complete treatment can be worked out for adsorbing atoms which behave ideally. Ideal behavior will be defined as AG,r = 0 for 0 5 2
zs
yI -
I 1. This casewill be treated
for both dissociative and non-di~o~iative adsorption. An attempt will be made, later, to comment on the perturbing influence that AGir # 0 will have on ideal behavior.
+ kTSN,, In (1 + P/P,)
4% -_=
d 1nP
-kTI’,
(dT = 0)
where yk may be yO, yE or yr (from equation (Z)), and where the corresponding surface excess quantities rk are defined by l?,=N,-NN,rNs
(I3a)
rc = N, - XN,
(13b)
I’r = 2I’, -
I’& = 2N, -
equations
(lo),
N& - XN,
(13~)
(12) and (13a), the
* If the ratio P/P, is replaced by the concentmtion ratio C/C,, equation (10) applies equally well to solute adsorption from the crystal interior.
(14~)
where yO,, ylo and yr,, refer to surface, ledge and ledge in~r&ction energies at P = 0, that is, for a perfectly clean surface. Making use of the fact that the total number of all types of sites N, available for adsorption is a constant given by
the final expression can be arrived at by combining equations (2), (14) and (15) y -
y,,,,= -kTN, + i
In (1 + P/P,) yzo -
kTN,,
+ kT(N,, + IIaving obtained equations, giving the degree of occupancy of each type of site as a function of the variables of the system (P, AGiO, T), it is now necessary to relate the occupation probabi~ty to the change in surfaoe free energy, y. This can be done through the use of the Gibbs-Adsorption Isotherm, which may be stated in the form:
(14b)
yro = --BkTN,, In (1 + P/PJ + kTN, In (1 + P/P,)
AG: = - RT In (f,#JJ
(144
remembering that for ideal behavior P, is not a function of N,/N,,. Similarly the expressions for yI and yr are found to be yt -
Combining
961
(I)
-
fiN,J ln (1 +
f’/PJ
PIPJ]
-YOOa2 2 a2
- SrX -
In (1 +
[
‘yr, + kTNd, ln (1 + PIP,)
2kTN,, In (1 + P/P,)
+ kT[XiWN,,
-NT, .,
-
flJf)I In (1 + P/P,)
1
+ -t?-
-
SMln (I + P/P,)]
(16)
-
N,, + Sfl,(l
XSM2
[ETN,,
Some simplification of equation (16) is possible. The third and fourth order terms are present only when P > 0. They have an appreciable magnitude only when ccis large and when there is preferential adsorption at surface sites, P, < fPd, P&. When ot is large, the model used here should be least applicable, particularly without the inclusion of higher order terms in the y vs. M:relationship for P = 0. Accordingly, they will be ignored in deriving expressions for the slope
ACTA
962
METALLURGICA,
and intercept of the torque term plot from equations (3) and (4). In the equations that follow, N,, = N,,.
(17) 2kTN,, +----S2MY,0 f12MY,
23=1+-
X Pn (1 + P./PA + e
[(SM -
i
00
2 ln (1 + P/P,)]
l)(iVz, - SN,)] lnfl -
VOL.
11,
1963
can be viewed as an interaction energy for ith sit,e.
A comparison of equations (10) and (21), for the case where P, = P,,, it shows that at small degrees of occupation equation (21) gives much higher values of N~lNi~. may be derived following Expressions for ylc - ylEO the integration procedure used previously, but for sake of brevity they will not be given here. Moreover, the analog of equation (16), which is very similar to it in appearance will be omitted. If the final result is simplified in the manner discussed in connection with equation (lci), the parameters A and 3 that characterize the y-plot are given by
PIP,) (18)
Adsorption
with dissociation
Allowing the ith site to be represented by Si, a prototype equilibrium for dissociative adsorption can be written as G, + 2X,*
2(G - SJ
+ rg
so that equation (8) becomes --.aoi
2 =
exp
[(~/GO
+
~~2’
-
‘g2’i2
2,Uu,i”)/RT]
= exp (-_G,,‘/RT)
(20)
AG,,O is the adso~tion site for the dissociative case free energy of the itsh when the species are in their standard states. If use is made of the definitions given in equations (9), the adsorption isotherm for the dissociative case can be expressed in the form where 2,~’
-
2pio -
,ugzo =
ENJOY,, -
2 In [l +
SN,) In [I + (PIJ’sa)1’21 (23)
+ 4kTN,, m
oo
(In [l + (fYPdf21
tPI'P,,P21) -tg
00 x [(EM - l)(N,, + SN,)] In [l + (P/f’,,)l~z]
The free energy change (SC result from the transfer of Sng2 diatomic molecules from the gas phase to an, sites on the surface is the same as given in equation (6), except in this case 6?z,z = 2&b, = 2&Z.,,
2~1,
B,=1$-
(19)
(N,, 00
(24) DISCUSSION
According to the above t~atment, there are three major effects that an adsorbing gas will have on the shape of the y-plot in the vicinity of an energy cusp. Each effect comes into prominence as the ambient pressure P approaches the characteristic pressure of the site. In the analysis that follows it will be necessary to treat limiting cases, i.e., cases where one (or two) of the ~harac~risti~ pressures is much larger than the others. Such limiting cases are not without experimental interest. For example, if the AG,O’s are separated by 10 kcal/mole-a figure that is not unreasonable in view of the fact that heats of adsorption
may be as large as 150 kcal/mole-the
characteristic pressures will differ by approximately two orders of magnitude at 1000°K. Consider first the ease where there is enhanced binding at the ledge
where the characteristic pressure P,, of the ith site for the dissociative case can be written as P
id
= L!@- Pli2 exp (AG.ad“/RT) f$ a
Analogously to the previous case, RT
soi_ lnfTi = -
(221
and double ledge sites. Case I:
P, > (Pd, PI) or P,, > (P,,, P,,) As the ambient pressure P is increased from zero the first major change in the y-plot occurs in its as reflected by changes in B. In equations (18) and (24) the terms depending on P, or P,, will
curvature
AG:”
be neglected,
and for convenience
the following
new
GJOSTEIN:
ADSORPTION
AND
SURFACE
ENERGY
963
(I)
variables are defined. 2.0
x=PIP,
Xd = PIP,, 1.5
R=
Pa/P1
With these definitions it is possible to write equations (18) and (24) in the form
t
1.0
a? 2
0.5
0 al" 0 I 0 " -0.5 -1.0
B, = l3,, + 2I3, In
-1.5
1 + x2’2
I
(25b)
II 3 (fi&Pli
The parameters K and K, are a measure of the relative binding energies of the single and double ledge sites. As (Ii, K,) -+ 0 double ledges tend to have the dominant binding energy, which means that the concentration of adsorbed atoms in single ledge sites can be neglected with respect to that in double ledge sites. Figs. 4 and 5 give plots of (B - B,,) /3,and (Bd - 3,,)/2B,,vs. x and xd, for various values of K and K,. For the non-dissociative case, when 0 < K < 0.5, the curves in Fig. 4 show a maximum, which is located at x, = (1 - 2K)lK. For x > x, = (1 - 2K)jK2, (B - B&B, -=c 0. When K > 0.5, (B - &J/B, is never positive and decreases monotonically with increasing x due to the dominance of the single ledge term. In the dissociative case, a maximum at xdm = (1 - 2K,1'2)2/K, appears in the curves in Fig. 5, when 0 < K, < 0.25, and ( Bd - B,,)/2B, < 0 for all
0
5
IO
I5
20 _ 25
30
35
40
45
50
x,j’YPdd FIG. 5. Effect of reduced pressure (X, = P/P& on yz-plot curvature parameter (Rd) for dissociative adsorption. Kd( =P,/P,,) is a measure of the relative adsorption energies of single and double ledge sites. x& when Kd > 0.25. If a comparison is made between dissociative and non-d~sociativ~ cases, by setting x = x~, it can be seen that for a given value of (x, zd), the increase (or decrease, depending on K) in curvature of the y-plot is somewhat greater for non-dissociative adsorption, and that for a given K = Kd value, the maximum occurs earlier for dissociative adsorption. The second major change in the shape of the y-plot occurs in the slope of the energy cusp at its origin, (A, Ad), as is evident from equations (17) and (23). To describe this effect, use is made of the definitions
so that equations (17) and (23) can be rewritten as A=&,-&ln(l
+Kx)
(264
A, = A,,- 2A, In [ 1 + (X,Z#~]
0
I
I
I
,
I
I
I
,
I
IO
20
30
40
50
60
70
60
90
Xf
IO0
P%
pl-3SSUre {x = P/P,) On yz-plot curvature parameter (B), for non-dissociative adsorption. K(= Pa/P,) is a measure of the relative adsorption energies of single and double ledge sites.
FIG.
4.
Effect
Of
lY?dUC%d
(26b) Fig. 6 gives plots of (A - A,,)/A, and (A, - A,)/2A, vs. it: and xcl for various values of K and K,. From equation (26) it is clear that at some value of (x, x& (A, Ad) will become zero and finally negative. This will occur when the logarithmic functions attain a value of A,,/A, (or A,,/2A,). Generally it can be said that when (K,KJ < 1, the effect of adsorption on curvature of the plot. y-plot occurs at much smaller values of (x, z8) than does the effect on the slope at the origin, but that when K - 0.5(orK, -+ 0.25) both effects tend to be present simultaneously. It was mentioned above that it is possible for (A,Ad) to become negative*. When this is the case has discussed surface energy.
* Herring’5)
negative
the
possibility
of having
a
ACTA
964
‘CA,
~ETALLURG
: B
- 1.0
I I
11,
1963
conditions, ledges can bunch together, resulting in a phenomenon known as thermal faceting. The criterion for this to occur will be discussed in an accompsnying paper.(8) A summary of the adsorption effects for Case I is shown schematically in Fig. 8.
0
t 8
VOL.
----_-_
s
N
-2.0 P”0
b
::
a* I 4
4
P -
-3.0
P-
- 4.0 t
1 0
Pd
P,
I 11
20
1
$1 40
x.Por
pd
I 60
x d
11
I
1 80
IO0
-p-m pdd
?&a.
6. Effeot of reduced pressure (X, X,) an the slope of y2-piot at the vertex of the cusp (A, A?). (K, K,) is a mewmro of the relative adsorption energies of single and double ledge sites. The dashed curves represent dissociative adsorption and the solid curves represent non-dissociative adsorption.
an outward pointing cusp will appear in the y-plot. Fig. 7 schematic&~y shows this tradition. There is no ex~riment&l evidence indicting whether or not this effect actually occurs. At first glance, the concept of a negative ledge energy might be rejected, since apparently the ledges could increase in length spont&neously with it decrease in free energy of the system. The total surface energy is not negative, however, and any change in length of the ledges must be subject to the criterion that total surface energy of system be B minimum. TTndercertain thermodynamic
FIU. 7. y,-plots and corresponding Gibbs-Wulff shapes, depicting the trrsnsition from an inward to an outwardpointing cusp.
0 -Ct-
FIG. 8, Schematic ~umm5ry of major effects of adsorption on shape of y vs. cc plots.
Case II: P, < (Pd, PJ or Ptd < (P,,, P,,) A third major change in the y-plot occurs when the ambient pressure P approaches P, (or Pad). This effect which would not be important in Case I, results when there is stronger bonding of the adsorbed atoms on surface sites relative to the ledge sites. This is a complex effect, involving alterations in B, A snd yO. Examination of equations (18) and (24) shows that the important term controlling the magnitude of the curvature effect is, (XiW - l)(Nls + SN,),, which takes on values of -0.086 N,, or -0.17 N,, for the (100) and (111) planes, respectively. Thus, the curvature (B, B,) should decrease with increasing (x, Q) in a manner shown by the curves for (K, KS) = 0 in Figs. 4 and 5, when they are reflected about the zero baseline. Equations (17) and (23) show that the effect on (A, Ad) depends on the term, N,, + SN,, which takes on values of 1.71 N,, and 2.41 N,, for the (100) and (111) plane, respectively. Thus, in contrast to Case I, (A, A,) should increase with increasing (X, X,) for both the (100) and (111) cusp. The form the (A, Ad) vs. (X, X,) relationship can be seen by reflecting the curves in Fig. 6 about the zero baseline. These effects are essentially the reverse of Case I. They can be illustrated by viewing the sequence of curves, from bottom to top, in Fig. 8, as representing increasing pressure.
ADSORPTION
GJOSTEIN:
AND
SURFACE
ENERGY
TABLE 1. Adsorption of oxygen on copper at 1000°C
Heat of adsorption
- 110,000 ealfmole 0,
965
(I)
__I__ Dissociated coverage -_I
Dew point
PO,
Pebar&eristie
Undissociated coverage
- 36°C
3.5 x IO-“2
5 x 10-13
7 X IO-10
2.7 x 10-z
+23*c
3.5 x lo-‘8
5 x 10-13
7 x LO-6
2.7 x 10-a
The effect on y0 can be determined from equation (16). As is evident, y0 will decrease with increasing pressure by an amount kTN, In (1 + P/P,)* for non-dissociative adsorption or by an amount kTN, In [l + (P/PS)1i2] for dissociative adsorption. For the (160) plane in copper ~TN~ = 270 ergs/cm2 at lOOO”C, which means that by, = (0.1, 0.2)~~ when (P/P,, P/P,,) = 1. Since y0 decreases with increasing P, all the vicinal planes will be reduced by the same amount, and thus each point of the y-plot will shrink towards the origin.
The y-plot has been explored for two metals, copper’3) and nickel,(1t2) under entirely different ambient atmospheres. The nickel specimens were annealed at 1000°C in a vacuum of about low5 mm, while the copper specimens were annealed in a hydrogen atmosphere with a dew point of either -36°C or +23”C at temperatures in range 950-1020°C. The torque terms repotid by Robertson and Shewmon(3) for copper annealed under hydrogen are somewhat smaller than those for nickel. This might be viewed as an impurity effect, but as will be shown, there is a strong possibility that this is not the case, and that probably they are intrinsically lower. Dell et aZ.(lQ found the heat of adsorption of oxygen on copper at low coverages to be about Oxygen appears to be the - 110 kcaljmole 0,. most likely adsorbent on copper, since hydrogen has a much smaller, although uncertain, heat of adsorption (it is reported in the range -9.0 to -35.0 kcal./mole H,.) (17~18)Assuming a reasonable value for the entropy of adsorption of hydrogen, it should have a very small free energy of adsorption at elevated ~m~eratu~s (l~O°C~, and thus should have a negligible coverage, unless there is a change in bonding mechanism at these temperatures, which is not reflected in the heat measured at lower temperatures. Table 1 gives the degree of occupancy to be expected for a heat of adsorption of - 110 kcal/mole 0, and * These values c&r&ted by assuming that every lattice
site on the substrate can serve as sn adsorption site. Ledges were assumed to lie along close-packed directions in an f.c.c. lattice. Lscmann et aZ.(‘) were concerned soley with this term; they assumed P and N, to vary in a specific manner with the crystallographic plane under consideration.
entropy of adsorption of -30 e.u.7, under various conditions. It can be seen that the worst possible case is for a dissociated: molecule when the dew point is +23”C, where the degree of occupancy would be 2.7 x 10m3. Since there is no information available which would allow the charac~ristic pressure given in Table 1 to be associated with any particular site, the most unfavorable assumption will be made. This means the characteristic pressure should be equated to Pa,,. If this is true, the change in Bd will be, Bd - B,, = 5.4 x 1O-3 fz,, due to the terms containing P,,. For the (100) cusp in copper, B, is found to be at lOOO*C B, -
2kT Nzs _ 2kT 4 y0 PA!l y0 uf g(1.38 x lo-16)(1.273 x 103) = o 6.
=
(1.7 x 103)(3.61 x 10-8)2
’
For the (ill) cusp of copper, B, is a factor 0.43 smaller than this value. Thus it appears that Robertson and Shewmon essentially de~rmined B,,,. It could be argued, however, that I(AH,,O, AH,,O)I > j AH,“j, and that this larger heat was not detected calorimetrically. If this were the case there should be appreciable contribution to B, from terms containing (P,,, P,,). However, at the lower dew point this contribution should decrease considerably. This means they should have found B to vary with the partial pressure of oxygen whereas no appreciable effect on B was observed when the partial pressure of oxygen was decreased by four orders of magnitude. Similar reasoning can be applied to A,, and again no effect of adsorption seems likely. It must be concluded, then, that adsorbed oxygen did not influence the results of Robertson and Shewmon, and that it seems doubtful that adsorbed hydrogen exerted any influence either, although this latter point must be left open until further information is availab1e.g t The entropy of adsorption was assumed to be the same as entropy of formation Cu,O, namely -30 Cal/mole 0, - “K, since an atom loses essentially the same number of degrees of freedom in both the oxidation and adsorption processes. f It has not been established conclusively whether or not oxygen exists as dissociated atoms on metal surfaces. 3 If it is assumed that B,, W&Rmeasured, it is interesting to note that B, < 1, and therefore it follows that yro, i.e. there is a rep&non between single ledges. The significance of this f&zt is not understood as yet.
ACTA
966
Mykura
found
larger
the
parameters
A and B for nickel than those reported
for copper;
TABLE 2. Comparison
values
METALLURGICA,
of
A and B for nickel and
of parameters copper
11, 1963
VOL.
Deviations from ideality In any real system such as those discussed deviations
from ideal behavior
is independent
of degree of occupancy
and therefore,
of the pressure of the adsorbing
Very little is known theoretically
Nickel
copper
* Actually
so
that it is no longer possible to assume that Pi (or Pid)
A
(111) (100)
above,
can be expected,
gas.
about the deviations
0.17 0.30
0.24t 1.13
that might
not be treated within the scope of this paper.
Experi-
0.11 0.067
0.22 0.14
mentally,
heat
the slope of the experiment
be expected,
of given site,
it is known
adsorption
torque vs. GLplot is
which,
decreases
increasing
coverage
at any rate, could
generally (becomes
that
the
less negative)
of the surface.
The manner
which it decreases is not always simple.
rather than These two quantities are not much and nickel values, where yO/v is nearly unity and A is small. B’ is not constant for the (111) cusp in nickel. This value is for a = 0; B’ increases to 0.8 at cc = 0.3.
viewed
Pi increase
as having
occupancy
of the ifh site.
with
would
retain
the
degree
In this case,
essentially
the
in
This can be
changes in shape, which were described case,
of
with
of
the basic
for the ideal
same
character,
but will occur at higher pressures
(or lower temper-
that the values of A for copper are for a clean surface,
atures).
Given enough information
about the changes
it appears that either Case II-type
on
y-plot
a comparison
place
on nickel,
higher yet,
is given in Table
single
ledge
or that nickel
ledge
energies
principles
energy,
shape of small nickel
by Mykura.
This indicates
have
influenced
been
pressures Without
surface
pressure
it is not
cleanliness
is about
this
of the Gibbs-
by Sundquistc4)
adsorption. systems
The
possible to
the partial pressure of oxygen fraction
of
electron
diffraction
total
pressure
and
from ideal behavior
The phenomenological in the y-plot, has been
developed,
structure
theory
of shape alterations
due to adsorption
based
utilizing
from the gas phase, a model
on three types
that the slope of the y-plot
pressure in a degree
(15).
bulk In
of the adsorbing
decrease
there is enhanced
ledge sites.
binding
pressure
temperature)
of the adsorbing
when
atom at
ledge sites relative to the surface sites (Case I).
The
at
slope of the cusp origin can decrease to a point where
order
to
it becomes negative,
even
the surface
causing a cusp inversion.
sites have
the stronger
of nickel,
(Case II),
the slope of the y-plot
must be a very small
increasing
pressure of the adsorbing
of
sites: Theory
in the vicinity
with increasing
gas (decreasing
surface
NiO
(or probably oxygen?)
The
of a cusp should
of
of adsorption
predicts
1O-5 mm.
it
for a given system.
used by Mykura
set of data.
pressure
temperature,
should be possible to determine the extent of deviation
single ledge and double
to form
with oxidation
with
surface,
to assign
either
required
the
total
the same,
with chemisorbed
the
first
that the results for nickel
8 x 1O-s mm
problems
more critically
from
to decide
check on the partial
of oxygen
eliminate
As
terms than those found
were essentially
a further
system,
1000°C
crystals
by
in the vacuum
and Sundquist
copper.
determinations
much larger torque
takes
has an intrinsically
be calculated
Wulff
of
adsorption
ylO, than
cannot
Some preliminary
each
Since it appears
and thus it is not possible
issue.
indicate
2.
10e5 mm.
The
binding
will increase gas.
forces with
For Case I
the curvature
increases
worthog) show that carbon is also a possible adsorbent
goes through The maximum
a maximum and finally decreases. tends to vanish as the difference in
on nickel, in which case it is to be expected
binding
energy
becomes
large.
adsorption
studies
of
Schlier
and
Farnsthat CO
may be another problem.
There are no other data on the shape of the y-plot for metals,
except
for information
that can be ob-
tained from thermal faceting. This subject cussed in detail in an accompanying paper.(s)
is dis-
t This parenthetical comment was added because normally it would be expected that chemisorbed oxygen would be bonded more strongly than bulk oxide. However, Paravano et aZ.‘l*’ reported a heat of adsorption of - 111 kcal/mole, which is very close to the heat of formation of bulk oxide, - 116 kcal/mole.
with pressure
When
between When
single
critical value, the curvature with increasing curvature involved.
with
and
this difference
pressure. pressure
at first, then
double
ledge
drops below
a
decreases monotonically
Case II gives a decreasing increase;
no minimum
is
References 1. H. MYKURA, Acta Met. 9, 570 (1961). 2. J. M. BLARELY and H. MYKURA, Acta Met. 9. 595 (1961). 3. W. M. ROBERTSON and P. G. SHEWMCIN, Trans. Amer. Inst. Min. (M&U.) Engrs. 224, 804 (1962).
GJOSTEIN:
ADSORPTION
AND
4. B. E. SUNDQUIST. Private communication. 5. C. HERRING, Structure and Properties of Solid Surfaces. (Edited by R. GOMER and C. S. SMITH), p. 5. University of Chicago Press (1953). 6. C. HERRING, Phys. Rev. 82, 87 (1951). 7. R. LACMANN and I. N. STRANSKI, Growth and Perfection of Crystals. (Edited by R. H. DOREMUS, B. W. ROBERTS and D. TURNBULL), p. 427. Wiley, New York (1958). 8. N. A. GJOSTEIN. Acta Met., this issue p. 969. 9. W. K. BURTON, N. CABRERA and F. C. FRANK, Phil. Trans. Roy. Sot. Lond. 243A, 299 (1951). 10. W. W. MULLINS, Acta Met. 7, 746 (1959). 11. T. W. HICKMOTT and G. EHRLICR, J. Phys. Chem. Solids 5, 47 (1958). 12. E. BAUER, Phys. Rev. 123, 1206 (1961). 13. N. A. GJOSTEIN, Acta Met. 7, 812 (1959).
SURFACE
ENERGY
967
(I)
14. 0. D. GONZALEZ and G. PARRAVANO, J. Amer. Chem. Sot. 78, 4533 (1956). 15. F. D. RICHARDSON and J. H. E. JEFFES, J. Iron St. Inst. 160, 261 (1948). 16. R. M. DELL, F. S. STONE and P. F. TILEY, Trans. Fu”aTaday Sot. 49, 201 (1953). 17. T. TAKEUCHI and
M.
SAKAGUCHI,
Japan 39, pp. 177-182 (1957). 18. T. KWAN, Advances in Catalysis 6, 67.
Bull.
Chem. Sot..
Academic
Press,
New York (1954). 19. R. E. SCHILIER and H.
E. FARNSWORTH, Advances in Catalysis 9, 434. Academic Press, New York (1957). 20. The author is indebted to Professor G. M. POUND for discussions on this problem, and to Drs. H. MYKURA and J. HIRTR for communications, which help crystallize some of these ideas.