Two-dimensional colloidal aggregation: concentration effects.

Abstract

Extensive numerical simulations of diffusion-limited (DLCA) and reaction-limited (RLCA) colloidal aggregation in two dimensions were performed to elucidate the concentration dependence of the cluster fractal dimension and of the different average cluster sizes. Both on-lattice and off-lattice simulations were used to check the independence of our results on the simulational algorithms and on the space structure. The range in concentration studied spanned 2.5 orders of magnitude. In the DLCA case and in the flocculation regime, it was found that the fractal dimension shows a linear-type increase with the concentration phi, following the law: d(f)=d(fo)+aphi(c). For the on-lattice simulations the fractal dimension in the zero concentration limit, d(fo), was 1.451+/-0.002, while for the off-lattice simulations the same quantity took the value 1.445+/-0.003. The prefactor a and exponent c were for the on-lattice simulations equal to 0.633+/-0.021 and 1.046+/-0.032, while for the off-lattice simulations they were 1.005+/-0.059 and 0.999+/-0.045, respectively. For the exponents z and z', defining the increase of the weight-average (S(w)(t)) and number-average (S(n)(t)) cluster sizes as a function of time, we obtained in the DLCA case the laws: z=z(o)+bphi(d) and z'=z'(o)+b'phi(d'). For the on-lattice simulations, z(o), b, and d were equal to 0.593+/-0.008, 0.696+/-0.068, and 0.485+/-0.048, respectively, while for the off-lattice simulations they were 0.595+/-0.005, 0.807+/-0.093, and 0.599+/-0.051. In the case of the exponent z', the quantities z'(o), b', and d' were, for the on-lattice simulations, equal to 0.615+/-0.004, 0.814+/-0.081, and 0.620+/-0.043, respectively, while for the off-lattice algorithm they took the values 0.598+/-0.002, 0.855+/-0.035, and 0.610+/-0.018. In RLCA we have found again that the fractal dimension, in the flocculation regime, shows a similar linear-type increase with the concentration d(f)=d(fo)+aphi(c), with d(fo)=1.560+/-0.004, a=0.342+/-0.039, and c=1.000+/-0.112. In this RLCA case it was not possible to find a straight line in the log-log plots of S(w)(t) and S(n)(t) in the aggregation regime considered, and no exponents z and z' were defined. We argue however that for sufficiently long periods of time the cluster averages should tend to those for DLCA and, therefore, their exponents should coincide with z and z' of the DLCA case. Finally, we present the bell-shaped master curves for the scaling of the cluster size distribution function and their evolution when the concentration increases, for both the DLCA and RLCA cases.