Simulations of heteroaggregation in a suspension of alumina and silica particles: Effect of dilution

Electrostatic Heteroaggregation: Fundamentals and Applications in Interfacial Engineering

The aggregation of oppositely charged soft materials (particles, surfactants, polyelectrolytes, etc.) that differ in one or more physical or chemical attributes, broadly referred to as electrostatic heteroaggregation, has been an active area of research for several decades now. While electrostatic heteroaggregation (EHA) is relevant to diverse fields such as environmental engineering, food technology, and pharmaceutical formulations, more recently there has been a resurgence to explore various aspects of this phenomenon in the context of interface stabilization and the development of functional materials. In this Feature Article, we provide an overview of the recent contributions of our group to this exciting field with particular emphasis on fundamental studies of electrostatic heteroaggregation between oppositely charged systems in the bulk, at interfaces, and across the bulk/interface. The influence of the size and shape of particles and the surface charge of heteroaggregates on the formation of Pickering emulsions and their utilization in the development of porous ceramics is discussed.

Impact of Charge Regulation on Self-Assembly of Zwitterionic Nanoparticles

Zwitterionic modification of colloids with weak acids and bases represents a promising strategy in creating functional materials with tunable properties and modeling the self-organization of charged proteins. However, accurate incorporation of the dynamic dissociation or association of ionization groups known as charge regulation (CR) is often intractable in theoretical and computational investigations since charge redistribution and configuration need to be evolved self-consistently. Using hybrid Monte Carlo and molecular dynamics simulations, we demonstrate that a dilute suspension of overall charge-neutral zwitterionic Janus nanoparticles shows a conformational transition from an open assembly of string or bundle to compact cluster along with the variation in pH. The behavior under CR is qualitatively different from the commonly employed constant charge condition where the transition is absent. The CR-induced clustering is due to the inhomogeneous and fluctuating charges localized near the equatorial boundary of the Janus particle. These features are enhanced particularly at low salt concentration and high electrostatic coupling strength. Our results indicate the critical role of charge regulation in the spatial self-organization of zwitterionic nanoparticles.

Study of the aggregation behavior of Janus particles by coupling experiments and Brownian dynamics simulations

HYPOTHESIS

New colloids such as inverse patchy particles or Janus particles are considered as smart building blocks in the development of innovative and performant materials. For example, the control of the self-assembly of oxide-based charged Janus particles is interesting for ceramic shaping. Thus, the synthesis of silica based Janus particles as well as a detailed study of their behavior in suspension are presented in this paper.

EXPERIMENTS

Fluorescent silica particles are partially modified in surface by grafting amine groups using a Pickering emulsion route. Zeta potential measurements, sedimentation tests and confocal microscopy observations are carried out to analyze the aggregation of the obtained particles in aqueous suspension as a function of the patch size and of the pH. Brownian dynamics simulations are also performed to better understand the aggregate structures.

FINDINGS

The aggregation of the synthesized silica-based Janus particles can be tuned by modifying the experimental parameters, and elongated or on the contrary more compact structures could be observed. This control of aggregation makes such particles promising to build new ceramic materials.

Shear-Induced Heteroaggregation of Oppositely Charged Colloidal Particles

This paper investigates numerically the shear-induced aggregation of mixed populations of colloidal particles leading to the formation of clusters. Suspensions with different amounts of positively and negatively charged colloidal particles are simulated. To resolve the aggregation kinetics and structural properties of the formed clusters, we resort to a mixed deterministic-stochastic simulation method. The method is built on a combination of a Monte Carlo algorithm to sample a statistically expected sequence of encounter events between the suspended particles and a discrete element method built in the framework of Stokesian dynamics to simulate the encounters in a fully predictive manner. Results reveal a strong influence of the composition of the population on both the aggregation kinetics and the aggregate structure. In particular, we observe a size-stabilization phenomenon taking place in the suspension when the relative concentration of the majority particles lies in the range 80-85%; i.e., starting from primary particles, after a short growth period, we observed a cessation of aggregation. Inspection of the aggregate morphology shows that the formed clusters are composed of few minority particles placed in the inner region, while the aggregate surface is covered by majority particles, acting to provide a shielding effect against further growth.

Computer simulations of heteroaggregation with large size asymmetric colloids

HYPOTHESIS

Hetero-aggregation of inorganic colloids is influenced by numerous parameters, which dictate the suspension properties. When particles are different in size, the suspension can be either stable or unstable according to concentration of components, ionic strength, and pH. Experimentally, understanding the role of each parameter is sometimes difficult because parameters cannot easily be modified independently. Numerical simulations are thus very useful to discriminate between different effects.

SIMULATIONS

Brownian dynamics simulations are used here to study the heteroaggregation of dilute suspensions composed of two populations of colloids with large size asymmetry. Special attention is paid to the effect of small-particle concentration, surface potentials, and ionic strength.

FINDINGS

The simulation results show that hetero-aggregation can be tuned by modifying these different parameters, and that the resulting aggregate structures depend more on the surface properties of small particles than on those of large particles. The simulations shed light on a further parameter crucially influencing hetero-aggregation, i.e. the mobility of small particles when adsorbed on large ones. The present results rationalize numerous experimental observations reported in the literature and can be used as reference to explain future experimental observations.

Aggregation and fragmentation in liquids with dispersed nanoparticles

Nanoparticle-induced aggregation and fragmentation phenomena in liquid media are investigated by applying a model of preferential attachment of dispersing molecules to randomly chosen nanoparticles and larger particles, each containing a single nanoparticle.

Nanoparticle-Induced Gelation of Bimodal Slurries with Highly Size-Asymmetric Particles: Effect of Surface Chemistry and Concentration

A systematic study has been performed to investigate the effect of surface potential of nanoparticles on the rheological behavior of bimodal suspensions, using a model system consisting of polystyrene latex (primary size ∼530 nm) and alumina-coated silica (primary size ∼12 nm) particles. The surface potential of small particles was tuned by varying the solution pH, causing them to be repulsive to each other, attractive to each other, and oppositely charged to the large particles, while the large particles remained electrostatically stabilized. We found that the rheological properties could be dramatically changed from viscous to gel-like depending on the surface potential and concentration of small particles. A colloidal gel was induced by small particles when the small particles had the opposite charge to the large particles and a volume fraction of 10(-4) < ϕsmall < 10(-3), and when the small particles were attractive to each other above a critical threshold, ϕsmall > 10(-4). Cryo-SEM distinguished the gel structures to be either short bridging gels produced by oppositely charged small particles or long bridging gels or dense gels produced by attractive small particles. On the basis of this rheological behavior and microstructure, we prepared a phase diagram of highly size-asymmetric bimodal colloids with respect to the surface chemistry and concentration of small particles.

Numerical and experimental study of suspensions containing carbon blacks used as conductive additives in composite electrodes for lithium batteries

Suspensions of carbon blacks and spherical carbon particles are studied experimentally and numerically to understand the role of the particle shape on the tendency to percolation. Two commercial carbon blacks and one lab-synthesized spherical carbon are used. The percolation thresholds in suspensions are experimentally determined by two complementary methods: impedance spectroscopy and rheology. Brownian dynamics simulations are performed to explain the experimental results taking into account the fractal shape of the aggregates in the carbon blacks. The results of Brownian dynamics simulations are in good agreement with the experimental results and allow one to explain the experimental behavior of suspensions.

Structure and rheology of colloidal particle gels: insight from computer simulation

A particle gel is a network of aggregated colloidal particles with soft solid-like mechanical properties. Its structural and rheological properties, and the kinetics of its formation, are dependent on the sizes and shapes of the constituent particles, the volume fraction of the particles, and the nature of the interactions between the particles before, during and after gelation. Particle gels may be permanent or transient depending on whether the colloidal forces between the aggregating particles lead to irreversible bonding or weak reversible interactions. With short-range reversible interactions, network formation is typically associated with phase separation or kinetic arrest due to particle crowding. Much existing knowledge has been derived from computer simulations of idealized model systems containing spherical particles interacting with well-defined pair potentials. The status of current progress is reviewed here by summarizing the underlying methodology and key findings from a range of simulation approaches: Monte Carlo, molecular dynamics, Brownian dynamics, Stokesian dynamics, dissipative particle dynamics, multiparticle collision dynamics, and fluid particle dynamics. Then it is described how the technique of Brownian dynamics simulation, in particular, has provided detailed insight into how different kinds of bonding and weak reversible interactions can affect the aggregate fractal structure, the percolation behaviour, and the small-deformation rheological properties of network-forming colloidal systems. A significant ongoing development has been the establishment and testing of efficient algorithms that are able to capture the subtle dynamic structuring effects that arise from effects of interparticle hydrodynamic interactions. This has led to an appreciation recently of the potentially important role of these particle-particle hydrodynamic effects in controlling the evolving morphology of simulated colloidal aggregates and in defining the location of the sol-gel phase boundary.

Kinetically driven ordered phase formation in binary colloidal crystals

The aggregation of binary colloids of the same size and balanced charges is studied by Brownian dynamics simulations for dilute suspensions. It is shown that, under appropriate conditions, the formation of colloidal crystals is dominated by kinetic effects leading to the growth of well-ordered crystallites of the sodium-chloride (NaCl) bulk phase. These crystallites form with very high probability even when the cesium-chloride (CsCl) phase is more stable thermodynamically. Global optimization searches show that this result is not related to the most favorable structures of small clusters, which are either amorphous or of the CsCl structure. The formation of the NaCl phase is related to the specific kinetics of the crystallization process, which takes place by a two-step mechanism. In this mechanism, dense fluid aggregates form at first and then crystallization follows. It is shown that the type of short-range order in these dense fluid aggregates determines which phase is finally formed in the crystallites. The role of hydrodynamic effects in the aggregation process is analyzed by stochastic rotation dynamics - molecular dynamics simulations, and we find that these effects do not play a major role in the formation of the crystallites.