We have performed extensive numerical simulations of diffusion-limited (DLCA) and reaction-limited (RLCA) colloid aggregation to obtain the dependence on concentration of several structural and dynamical quantities, among them the fractal dimension of the clusters before gelation, the average cluster sizes, and the scaling of the cluster size distribution function. A range in volume fraction φ spanning two and a half decades was used for this study. For DLCA, a square root type of increase of the fractal dimension with concentration from its zero-concentration value was found: ${\mathit{d}}_{\mathit{f}}$=${\mathit{d}}_{\mathit{f}}^0$+a${\mathrm{\phi }}^{\mathrm{\beta }}$, with ${\mathit{d}}_{\mathit{f}}^0$=1.80±0.01, a=0.91±0.03, and β=0.51±0.02. For RLCA the same type of behavior was found, this time with ${\mathit{d}}_{\mathit{f}}^0$=2.10±0.01, a=0.47±0.03, and β=0.66±0.08. In the case of DLCA, the exponent z that defines the power law increase of the weight-average cluster size (${\mathit{S}}_{\mathit{w}}$) with time also increases as a square root type with concentration: z=${\mathit{z}}^0$+b${\mathrm{\phi }}^{\mathrm{\alpha }}$, with ${\mathit{z}}^0$=1.07±0.06, b=3.09±0.22, and α=0.55±0.03, while the exponent z′ that describes the power law increase of the number-average cluster size (${\mathit{S}}_{\mathit{n}}$) with time follows the same law: z′=${\mathit{z}}^{\prime 0}$+b′${\mathrm{\phi }}^{\mathrm{\alpha }\prime }$, now with ${\mathit{z}}^{\prime 0}$=1.05±0.04, b′=3.41±0.24, and α′=0.46±0.02. We have also found that the cluster size distribution function scales as ${\mathit{N}}_{\mathit{s}}$(t)≊${\mathit{N}}_0$${\mathit{S}}_{\mathit{w}}^{\mathrm{-}2}$f(s/${\mathit{S}}_{\mathit{w}}$), where ${\mathit{N}}_0$ is the number of initial colloidal particles and f is a concentration-dependent function displaying an asymmetric bell shape in the limit of zero concentration. For RLCA, we found an exponential increase of the average cluster sizes for a substantial range of the aggregation time: ${\mathit{S}}_{\mathit{w}}$∼${\mathit{e}}^{\mathit{p}\mathrm{\phi }\mathit{t}}$ and ${\mathit{S}}_{\mathit{n}}$∼${\mathit{e}}^{\mathit{q}\mathrm{\phi }\mathit{t}}$, with p≊2q. For longer times the behavior departs from the exponential increase and, in the case of ${\mathit{S}}_{\mathit{w}}$ for low concentration, it crosses over to a power law increase. In the RLCA case the scaling is as in DLCA where now a power law decay of the function f defines the exponent τ, f(x)∼${\mathit{x}}^{\mathrm{-}\mathrm{\tau }}$g(x), with g(x) decaying exponentially fast for x>1. A slight dependence of the exponent τ on concentration was computed around to the value τ=1.5. © 1996 The American Physical Society.

• Received 15 March 1996

DOI:https://doi.org/10.1103/PhysRevE.54.5456