Colloids and Surfaces A: Physicochem. Eng. Aspects 436 (2013) 325–332 Contents lists available at ScienceDirect Colloids and Surfaces A: Physicochemical and Engineering Aspects jo ur nal ho me page: www.elsevier.com/locate/colsurfa A population balance equation model to predict regimes of controlled nanoparticle aggregation Anand K. Atmuri a , Michael A. Henson a , Surita R. Bhatia a,b,c,∗ a Department of Chemical Engineering, University of Massachusetts, Amherst, MA 01003, United States b Department of Chemistry, Stony Brook University, Stony Brook, NY 11794, United States c Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, NY 11793, United States h i g h l i g h t s • Formation of stable finite-sized clus- ters observed in charged colloidal dispersions. • A population balance equation (PBE) model was developed to capture this phenomenon. • Model successfully predicts aggrega- tion regimes and final aggregate size. g r a p h i c a l a b s t r a c t a r t i c l e i n f o Article history: Received 15 February 2013 Received in revised form 27 June 2013 Accepted 3 July 2013 Available online 12 July 2013 Keywords: Cluster Aggregation Aggregate Processing Theory Modeling a b s t r a c t Forming stable clusters or aggregates of nanoparticles is of interest in a number of emerging applications. While formation of unstable fractal aggregates and flocs has been well-studied with both experiments and theory, the conditions that lead to stable, finite-sized clusters is not as well understood. Here, we present an integrated experimental and modeling study to explore aggregation in concentrated attrac- tive colloidal suspensions. A population balance equation (PBE) model is used to predict the aggregation dynamics of quiescent colloidal suspensions. A DLVO (Derjaguin–Landau–Verwey–Overbeek) type poten- tial is used to describe the interparticle potential, with attractive interactions arising from van der Waals forces and long-range repulsive interactions caused by electrostatics. The PBE model includes a full cal- culation of stability ratio variations as a function of aggregate size, such that the energy barrier increases with increasing size. As the ionic strength is decreased, the model predicts three regimes of behavior: uncontrolled aggregation into large flocs, controlled aggregation into stable clusters, and no aggrega- tion. The model is tested experimentally using latex particles at different salt concentrations and particle concentrations. When the Hamaker constant and surface potential are fit to aggregate size measure- ments collected at one salt concentration, the model accurately predicts the final mean aggregate size and regimes of aggregation at other salt concentrations and the same particle concentration. This result suggests that van der Waals and electrostatic forces are the dominant particle interactions in determining the final aggregate state. The mean aggregate size and aggregation regimes at different particle concen- trations could be accurately predicted by adjusting the surface potential. This parameter adjustment is consistent with the expectation that increasing colloid weight fractions cause aggregates to have a more fractal nature and hence have a lower effective repulsion. However, the model predicts much faster aggregation rates than what are observed experimentally. This discrepancy may be due to hydrodynamic effects or another slow dynamical process which is not accounted for in the model. Nevertheless, this study presents the first PBE model that can successfully predict stable aggregate size and aggregation regimes of charged colloidal particles over a range of salt concentrations and particle concentrations. © 2013 Elsevier B.V. All rights reserved. ∗ Corresponding author at: Department of Chemistry, Stony Brook University, Stony Brook, NY 111794, United States. Tel.: +1 631 632 7788. E-mail addresses: surita.bhatia@stonybrook.edu, sbhatia@ecs.umass.edu (S.R. Bhatia). 0927-7757/$ – see front matter © 2013 Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.colsurfa.2013.07.002