148 J. Phys. Chem. B, Vol. 106, No. 1, 2002
Rotstein et al.
systems in which polymer degradation and/or cross-linking are
extensive, and those results will be published soon. All the
reactions described here are considered at 90 °C and in the
presence of N2.
The overall chemistry of the system will therefore be a
combination of both pathways described above and will be
determined by the type of the reactive groups of the polymer
in solution and the nature of the metallic species formed. The
overall chemical reaction considered in this work is given by
The process can therefore be modeled by the following system
of nonlinear differential equations:
3
2,33
j-1
∞
1
2
c˘ j )
[Rj-k,k cj-k ck - Q
j-k,k
ck] -
[Rj,k c c - Q cj+k] - δ1,j c1 p
j,k
j
k
∑
∑
k)1
k)1
{
p˘ ) - c p
1
(4)
∞
j)1
for j ) 1, 2, ..., together with initial conditions satisfying ∑
j cj(0) ) 1 and p(0) ) p0. In eq 4, p0 is the initial dimensional
concentration of reactive sites divided by C, Rjk is the
dimensional coefficient of coagulation divided by λ, and Qjk is
the dimensional coefficient of fragmentation divided by λC. All
are nonnegative and dimensionless. Since, as we pointed out
before, a pseudo-equilibrium state is reached as a consequence
of the fragmentation process and the decrease in the mobility
of large clusters, we assume that the concentration of clusters
whose size is larger than some finite N is negligible. This is
also reflected in our simulation, and therefore, the ∞ bound on
the summation in eq 4 can be replaced by N - j. Note that in
that case, the conservation of mass is assured. When the ∞ bound
is left in the summation, conservation of mass may be lost for
t g tc for some nonnegative tc (tc ) gelation time). In that
case, the model will be valid up until t ) tc.
As we pointed out before, and as it is reflected in the model,
the process of metal coagulation and cluster formation (inde-
pendent of the existence of polymer interactions) is a complex
series of reactions, which can be grouped into two categories:
(
m
0)
n[Co (CO) ] + PS f Co + Co - PS*
(3)
x
y
where the asterisk (*) stands for a metal-polymer complex
compound, x ) 1, 2, 4, or 6 and y ) 4, 8, 12, or 16, respectively,
y/x e 4 (according to the coordination chemistry of metal
30
carbonyl complexes ), and m ) j + k in eq 1. In the main
reaction, the clusters are formed from smaller reactions. The
Co-PS* complexes formed during this process are the products
of the secondary reaction described in eq 2, which competes
with the main coagulation-fragmentation process described in
eq 1. The extent of the secondary reaction depends on the
reactivity of the cobalt precursor and on the interaction
coefficient between the reactive cobalt species and the poly-
styrene molecules. Therefore, the rate of the main coagulation
reaction will depend not only on the concentration of the cobalt
precursor and the rate coefficients of coagulation and fragmen-
tation but also on the concentration of the available reactive
sites on the polystyrene in the system that are capable of
interacting with the cobalt precursor.
In this manuscript we present a mathematical model to
qualitatively understand the dynamics of the complex chemical
system under consideration. The model is as simple as possible,
but reflects the main chemical features of the system. To our
knowledge, no attempt to model chemical systems of this kind
has been reported to date.
To translate the physical system considered here into math-
ematical terms, we first scale the time variable by multiplying
it by the coefficient of reaction between cobalt fragments and
polystyrene, λ, a positive constant, and the initial cluster mass
C, also a positive constant. λ in this case is analogous to the
Flory-Huggins interaction parameter,31 and it represents the
strength of the interaction between the metal clusters and the
functional groups of the polymer. Then we define cj as the
concentration of Co clusters of size j divided by C, and similarly
p as the concentration of the reactive sites on the polystyrene,
whose initial value is assumed to be a linear function of the
polymer concentration, divided by C. For the model developed
in this paper, we assumed that the only significant irreversible
interaction between cobalt fragments and polystyrene will occur
for c1 with a coefficient of reaction λ, while for other clusters,
j g 2, the interactions with the polystyrene will be reversible
and weaker.30 Therefore, each c1 particle, which is involved in
this interaction, will be excluded from further participation in
the main coagulation reaction. Moreover, the chemical bonding
between the c1 particles and the polystyrene reduces the number
of reactive sites on the polymer and hence, limits the extent
and rate of the cobalt-styrene interaction. It is important to stress
that the overall chemistry of the system and the final particle
size distribution of the cobalt clusters are dependent on the
relative importance and contribution of the two competing
reactions. These considerations, as we will see later, will be
reflected in the results obtained from our model, which will
illustrate the main features of the dynamics of the system, while
disregarding a more exact and detailed description.
3
2
(a) The growth of clusters either by the formation of atom-
atom bonds between two metal fragments (metal fragment )
c1) and between a small cluster and a metal fragment, or by the
development of surface-surface interactions between two large
clusters; and (b) The decomposition of large metal clusters into
smaller fragments. To properly model these processes, it is
important to construct the appropriate reaction rate expressions
both for the cluster formation step and the cluster dissociation
step (with or without polymer interaction).
The simplest coagulation expressions found in the literature
R
R
that are pertinent for our model are Rjk ) j k (predominantly
34-40
R
R
chemical bond),
interactions),
and Rjk ) j + k (predominantly surface
where in both cases 0 e R e 1. In this model,
34,35
we recognize the fact that there is a decline in particle mobility
through the medium as a function of increased particle size.
Moreover, the presence of the polystyrene molecules in the
reaction solution contributes to the formation of a viscous
41,42
medium,
which will impose an upper bound on the ability
of larger clusters to interact via surface-surface interactions,
and hence these types of interactions will have a negligible
contribution to the overall coagulation process. These assump-
tions are consistent with the fact that in large clusters the fraction
of atoms on the surface is small compared to the total number
of atoms in the cluster, and therefore only this small fraction is
essentially coordinatively unsaturated. In effect, we can view
large clusters as being thermodynamically stable in real time.
This is not the case for small clusters whose surface atoms
represent a considerable fraction of the total metal mass. In this
case, the model has to promote the growth of particles formed
from smaller fragments, because this growth will minimize their
free energy by reducing the number of coordinatively unsatur-
4
3,44
ated surface atoms.
Therefore, the expression that we propose for the coefficients
of coagulation of the cobalt fragments in our model (with or