Novel kinetic trapping in charged colloidal clusters due to self-induced surface charge organization

Interaction between Hydrated Smectite Clay Particles as a Function of Salinity (0-1 M) and Counterion Type (Na, K, Ca)

Swelling clay minerals control the hydrologic and mechanical properties of many soils, sediments, and sedimentary rocks. This important and well-known phenomenon remains challenging to predict because it emerges from complex multiscale couplings between aqueous chemistry and colloidal interaction mechanics in nanoporous clay assemblages, for which predictive models remain elusive. In particular, the predominant theory of colloidal interactions across fluid films, the widely used Derjaguin-Landau-Verwey-Overbeek model, fails to predict the ubiquitous existence of stable swelling states at interparticle distances below 3 nm that are stabilized by specific inter-atomic interactions in overlapping electrical double layers between the charged clay surfaces. Atomistic simulations have the potential to generate detailed insights into the mechanisms of these interactions. Recently, we developed a metadynamics-based molecular dynamics simulation methodology that can predict the free energy of interaction between parallel smectite clay particles in a wide range of interparticle distances (from 0.3 to 3 nm) and salinities (from 0.0 to 1.0 M NaCl). Here, we extend this work by characterizing the sensitivity of interparticle interactions to counterion type (Na, K, Ca). We establish a detailed picture of the free energy of interaction of parallel clay particles across water films as the sum of five interaction mechanisms with different sensitivities to salinity, counterion type, and interparticle distance.

Thermodynamic Signatures of Structural Transitions and Dissociation of Charged Colloidal Clusters: A Parallel Tempering Monte Carlo Study

Computational simulation of colloidal systems make use of empirical interaction potentials that are founded in well-established theory. In this work, we have performed parallel tempering Monte Carlo (PTMC) simulations to calculate heat capacity and to assess structural transitions, which may occur in charged colloidal clusters whose effective interactions are described by a sum of pair potentials with attractive short-range and repulsive long-range components. Previous studies on these systems have shown that the global minimum structure varies from spherical-type shapes for small-size clusters to Bernal spiral and “beaded-necklace” shapes at intermediate and larger sizes, respectively. In order to study both structural transitions and dissociation, we have organized the structures appearing in the PTMC calculations by three sets according to their energy: (i) low-energy structures, including the global minimum; (ii) intermediate-energy “beaded-necklace” motifs; (iii) high-energy linear and branched structures that characterize the dissociative clusters. We observe that, depending on the cluster, either peaks or shoulders on the heat–capacity curve constitute thermodynamics signatures of dissociation and structural transitions. The dissociation occurs at T=0.20 for all studied clusters and it is characterized by the appearance of a significant number of linear structures, while the structural transitions corresponding to unrolling the Bernal spiral are quite dependent on the size of the colloidal system.

Protein-polymer mixtures in the colloid limit: Aggregation, sedimentation, and crystallization

While proteins have been treated as particles with a spherically symmetric interaction, of course in reality, the situation is rather more complex. A simple step toward higher complexity is to treat the proteins as non-spherical particles and that is the approach we pursue here. We investigate the phase behavior of the enhanced green fluorescent protein (eGFP) under the addition of a non-adsorbing polymer, polyethylene glycol. From small angle x-ray scattering, we infer that the eGFP undergoes dimerization and we treat the dimers as spherocylinders with aspect ratio L/D - 1 = 1.05. Despite the complex nature of the proteins, we find that the phase behavior is similar to that of hard spherocylinders with an ideal polymer depletant, exhibiting aggregation and, in a small region of the phase diagram, crystallization. By comparing our measurements of the onset of aggregation with predictions for hard colloids and ideal polymers [S. V. Savenko and M. Dijkstra, J. Chem. Phys. 124, 234902 (2006) and Lo Verso et al., Phys. Rev. E 73, 061407 (2006)], we find good agreement, which suggests that the behavior of the eGFP is consistent with that of hard spherocylinders and ideal polymers.

Real space analysis of colloidal gels: triumphs, challenges and future directions

Colloidal gels constitute an important class of materials found in many contexts and with a wide range of applications. Yet as matter far from equilibrium, gels exhibit a variety of time-dependent behaviours, which can be perplexing, such as an increase in strength prior to catastrophic failure. Remarkably, such complex phenomena are faithfully captured by an extremely simple model-'sticky spheres'. Here we review progress in our understanding of colloidal gels made through the use of real space analysis and particle resolved studies. We consider the challenges of obtaining a suitable experimental system where the refractive index and density of the colloidal particles is matched to that of the solvent. We review work to obtain a particle-level mechanism for rigidity in gels and the evolution of our understanding of time-dependent behaviour, from early-time aggregation to ageing, before considering the response of colloidal gels to deformation and then move on to more complex systems of anisotropic particles and mixtures. Finally we note some more exotic materials with similar properties.

Inverse design of equilibrium cluster fluids applied to a physically informed model

Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. We then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.

Self-Assembly of Colloidal Molecules Based on Host-Guest Chemistry and Geometric Constraints

The preparation of colloidal molecules (CMs), that is, clusters of colloids with a defined aggregation number and configuration, is of continued and significant interest in colloid chemistry and materials science and numerous interactions have been utilized to drive their (self-)assembly. However, only very few reports are available on the assembly of CMs based on host-guest chemistry. In this paper, we investigate the assembly of like-charged silica particles into well-defined, core-satellite AB_n-type CMs in water, mediated by host-guest interactions and geometric constraints. Exploiting the inherent dynamics of noncovalent attraction and making use of a soft polymer shell to enhance multivalent host-guest interactions, we successfully synthesized AB₃, AB₄, and AB₆ CMs by selecting the appropriate size ratio of satellite to core particles.

Opposed flow focusing: evidence of a second order jetting transition

We propose a novel microfluidic "opposed-flow" geometry in which the continuous fluid phase is fed into a junction in a direction opposite to the dispersed phase. This pulls out the dispersed phase into a micron-sized jet, which decays into micron-sized droplets. As the driving pressure is tuned to a critical value, the jet radius vanishes as a power law down to sizes below 1 μm. By contrast, the conventional "coflowing" junction leads to a first order jetting transition, in which the jet disappears at a finite radius of several μm, to give way to a "dripping" state, resulting in much larger droplets. We demonstrate the effectiveness of our method by producing the first microfluidic silicone oil emulsions with a sub micron particle radius, and utilize these droplets to produce colloidal clusters.

Hydrodynamic simulations of charge-regulation effects in colloidal suspensions

Self-organization of charged soft matter is of crucial importance in biology. However, it is an extremely complex phenomenon due to dynamical couplings between the hydrodynamic flow, ions, and charges of soft matter. For colloidal suspensions, the coupling between the former two has already been studied by numerical simulations, with the colloid surface charge being fixed. However, the self-organization of colloids and/or the application of an external electric field make the electrostatic environment of each colloid inhomogeneous in both space and time. Thus, this leads to inhomogenisation of the surface charge of each colloid under ionisation equilibrium conditions. This effect is known as "charge regulation" and is of great importance in various electrostatic and electrokinetic phenomena not only in colloid suspensions but also in solutions of biomolecules. However, there has so far been no success in taking the charge regulation effect into account in numerical simulations of colloidal electrokinetics. Here, we extend the fluid particle dynamics (FPD) method to incorporate the charge regulation effect. We present a theoretical formulation of the method and its application to two types of problems, where charge regulation plays an important role: (i) cluster formation of colloid particles and (ii) a single colloid particle under an external field. By these examples, we show not only the importance of considering charge-regulation effects in the self-organization of charged systems but also the applicability of our simulation method to more complex problems of charged soft matter systems.

Hunting mermaids in real space: known knowns, known unknowns and unknown unknowns

We review efforts to realise so-called mermaid (or short-ranged attraction/long ranged repulsion) interactions in 3d real space. The repulsive and attractive contributions to these interactions in charged colloids and colloid-polymer mixtures, may be accurately realised, by comparing particle-resolved studies with colloids to computer simulation. However, when we review work where these interactions have been combined, despite early indications of behaviour consistent with predictions, closer analysis reveals that in the non-aqueous systems used for particle-resolved studies, the idea of summing the attractive and repulsive components leads to wild deviations with experiment. We suggest that the origin lies in the weak ion dissociation in these systems with low dielectric constant solvents. Ultimately this leads even to non-centro-symmetric interactions and a new level of complexity in these systems.

Temperature- and field-induced structural transitions in magnetic colloidal clusters

Magnetic colloidal clusters can form chain, ring, and more compact structures depending on their size. In the present investigation we examine the combined effects of temperature and external magnetic field on these configurations by means of extensive Monte Carlo simulations and a dedicated analysis based on inherent structures. Various thermodynamical, geometric, and magnetic properties are calculated and altogether provide evidence for possibly multiple structural transitions at low external magnetic field. Temperature effects are found to overcome the ordering effect of the external field, the melted stated being associated with low magnetization and a greater compactness. Tentative phase diagrams are proposed for selected sizes.

Thermodynamic signatures and cluster properties of self-assembly in systems with competing interactions

Colloidal particles, amphiphiles and functionalized nanoparticles are examples of systems that frequently exhibit short-range attraction coupled with long-range repulsion. We vary the ratio of attraction and repulsion in a simple isotropic model with competing interactions, using molecular simulations, and observe significant differences in the properties of the self-assembled clusters that form. We report conditions that lead to the self-assembly of clusters of a preferred size, accompanied by a change in the slope of the pressure with respect to density, similar to micelles formed by amphiphilic molecules. We also report conditions where repulsion dominates, clusters of a preferred size form and the pressure vs. density slope is unaffected by self-assembly. We investigate cluster structure by calculating the size distributions, free colloid density, cluster shape and density profiles. The system dynamics are characterized by cluster life-times. We do not find qualitative differences in structure or dynamics of the clusters, regardless the pressure behavior. Therefore, thermodynamic and structural quantities are required to classify the different clustering characteristics that are observable in systems with competing interactions. Our results have implications in terms of development of design principles for stable cluster self-assembly.

Multistep Heterogeneous Nucleation in Binary Mixtures of Charged Colloidal Spheres

Nucleation plays a decisive role in determining the crystal structure and size distribution; however, understanding of the fundamentals of nucleation is quite limited. In particular, it is unclear whether a nucleus forms spontaneously from solution via a single- or multiple-step process. Here we show how a binary mixture of charged colloidal spheres nucleates heterogeneously on a flat substrate by means of Bragg microscopy, laser diffraction, and laser microscopy. In contrast with the conventional one-step and two-step nucleation mechanisms, a novel pathway of multistep heterogeneous nucleation under certain experimental conditions is highlighted by four steps: initial homogeneous fluid → prenucleation clusters → preordered prenucleation clusters → intermediate ordered phase → final crystal. It is expected that the obtained results would be helpful in rationalizing the rich phase behavior exhibited by the binary mixture systems and in developing better and broadly applicable nucleation models as well as in designing defect-free single-crystal alloys.

Structure of colloidal gels at intermediate concentrations: the role of competing interactions

Colloidal gels formed by colloid-polymer mixtures with an intermediate volume fraction (ϕ_c ≈ 0.4) are investigated by confocal microscopy. In addition, we have performed Monte Carlo simulations based on a simple effective pair potential that includes a short-range attractive contribution representing depletion interactions, and a longer-range repulsive contribution describing the electrostatic interactions due to the presence of residual charges. Despite neglecting non-equilibrium effects, experiments and simulations yield similar gel structures, characterised by, e.g., the pair, angular and bond distribution functions. We find that the structure hardly depends on the strength of the attraction if the electrostatic contribution is fixed, but changes significantly if the electrostatic screening is changed. This delicate balance between attractions and repulsions, which we quantify by the second virial coefficient, also determines the location of the gelation boundary.

Numerical study of cluster formation in binary charged colloids

Cluster formation of oppositely charged colloidal particles is studied numerically. A simple Brownian dynamics method with a screened-Coulomb (Yukawa) potential is employed for numerical simulations. An equilibrium phase which consists of clusters and unassociated particles is obtained. It is shown that the equilibrium association number of clusters and their shapes are determined by charge numbers and charge ratio of the binary particles. The phase diagram of cluster formation for various charge numbers and their ratios is obtained. A simple relation between the association number and the charge ratio is found. It is demonstrated that in the case of high charge ratio the cluster takes a multilayer structure which is highly symmetric. It is also pointed out that the cluster-particle interaction changes dynamically in the cluster formation process, which is involved in the selection of final cluster structure.

A Detailed Study on the Low-Energy Structures of Charged Colloidal Clusters

The target of this investigation is the systematic characterization of the low-energy structures of charged colloidal clusters that may be important to understand the self-assembling process of biomolecules. The aggregation of charged colloidal particles is governed by the attractive short-ranged Morse potential and the Yukawa repulsive tail to describe the long-range charge effect. A global optimization strategy, based on our own evolutionary algorithm, was adopted to discover the low-energy structures of colloidal clusters composed of up to 20 particles. A detailed analysis of the low-energy structures involving charged particles shows that the appearance of the Bernal spiral as the most stable motif occurs, first, at N = 6, but it is favored for larger clusters (N ≥ 13); for 6 ≤ N ≤ 12, there is a competition between the spiral (which is favored for higher charges) and more spherical-like structures. Finally, we study binary clusters composed by two sets of differently charged colloidal particles. Although a great diversity of low-energy structures is observed (especially for aggregates with one of the components in excess), the global minimum is disputed by three structural motifs depending on the composition of the cluster and, in some cases, on the range of the Morse potential.

Spontaneous Formation of Eutectic Crystal Structures in Binary and Ternary Charged Colloids due to Depletion Attraction

Crystallization of colloids has extensively been studied for past few decades as models to study phase transition in general. Recently, complex crystal structures in multi-component colloids, including alloy and eutectic structures, have attracted considerable attention. However, the fabrication of 2D area-filling colloidal eutectics has not been reported till date. Here, we report formation of eutectic structures in binary and ternary aqueous colloids due to depletion attraction. We used charged particles + linear polyelectrolyte systems, in which the interparticle interaction could be represented as a sum of the electrostatic, depletion, and van der Waals forces. The interaction was tunable at a lengthscale accessible to direct observation by optical microscopy. The eutectic structures were formed because of interplay of crystallization of constituent components and accompanying fractionation. An observed binary phase diagram, defined by a mixing ratio and inverse area fraction of the particles, was analogous to that for atomic and molecular eutectic systems. This new method also allows the adjustment of both the number and wavelengths of Bragg diffraction peaks. Furthermore, these eutectic structures could be immobilized in polymer gel to produce self-standing materials. The present findings will be useful in the design of the optical properties of colloidal crystals.

Equilibrium Phase Behavior of a Continuous-Space Microphase Former

Periodic microphases universally emerge in systems for which short-range interparticle attraction is frustrated by long-range repulsion. The morphological richness of these phases makes them desirable material targets, but our relatively coarse understanding of even simple models hinders controlling their assembly. We report here the solution of the equilibrium phase behavior of a microscopic microphase former through specialized Monte Carlo simulations. The results for cluster crystal, cylindrical, double gyroid, and lamellar ordering qualitatively agree with a Landau-type free energy description and reveal the nontrivial interplay between cluster, gel, and microphase formation.

Self-assembly of colloidal magnetic particles: energy landscapes and structural transitions

The self-assembly of colloidal magnetic particles is of particular interest for the rich variety of structures it produces and the potential for these systems to be reconfigurable.

Equilibrium cluster fluids: pair interactions via inverse design

Inverse methods of statistical mechanics are becoming productive tools in the design of materials with specific microstructures or properties. While initial studies have focused on solid-state design targets (e.g., assembly of colloidal superlattices), one can alternatively design fluid states with desired morphologies. This work addresses the latter and demonstrates how a simple iterative Boltzmann inversion strategy can be used to determine the isotropic pair potential that reproduces the radial distribution function of a fluid of amorphous clusters with prescribed size. The inverse designed pair potential of this "ideal" cluster fluid, with its broad attractive well and narrow repulsive barrier at larger separations, is qualitatively different from the so-called SALR form most commonly associated with equilibrium cluster formation in colloids, which features short-range attractive (SA) and long-range repulsive (LR) contributions. These differences reflect alternative mechanisms for promoting cluster formation with an isotropic pair potential, and they in turn produce structured fluids with qualitatively different static and dynamic properties. Specifically, equilibrium simulations show that the amorphous clusters resulting from the inverse designed potentials display more uniformity in size and shape, and they also show greater spatial and temporal resolution than those resulting from SALR interactions.

Controlled Clustering in Binary Charged Colloids by Adsorption of Ionic Surfactants

We report on the controlled clustering of oppositely charged colloidal particles by the adsorption of ionic surfactants, which tunes charge numbers Z of particles. In particular, we studied the heteroclustering of submicron-sized polystyrene (PS) and silica particles, both of which are negatively charged, in the presence of cetylpyridinium chloride (CPC), a cationic surfactant. The surfactant concentration Csurf was selected below the critical micelle concentration. As CPC molecules were adsorbed, Z values of the PS and silica particles decreased, inverting to positive when Csurf exceeded the isoelectric point Ciep. Hydrophobic PS particles exhibited much lower Ciep than hydrophilic silica particles. At Csurf valuess between their Ciep values, the particles were oppositely charged, and clustering was enabled. To explain the clustering behavior, we investigated adsorption isotherms of the CPC and screened-Coulomb-type pair potential. Expected applications of the present findings are the control of colloidal associations and construction of various particle types into heterogeneous colloidal clusters.

Structure and dynamics of model colloidal clusters with short-range attractions

We examine the structure and dynamics of small isolated N-particle clusters interacting via short-ranged Morse potentials. "Ideally prepared ensembles" obtained via exact enumeration studies of sticky hard-sphere packings serve as reference states allowing us to identify key statistical-geometrical properties and to quantitatively characterize how nonequilibrium ensembles prepared by thermal quenches at different rates T[over ̇] differ from their equilibrium counterparts. Studies of equilibrium dynamics show nontrivial temperature dependence: nonexponential relaxation indicates both glassy dynamics and differing stabilities of degenerate clusters with different structures. Our results should be useful for extending recent experimental studies of small colloidal clusters to examine both equilibrium relaxation dynamics at fixed T and a variety of nonequilibrium phenomena.

Free-energy landscapes of granular clusters grown by magnetic interaction

We experimentally study the aggregation of small clusters made of non-Brownian dipolar beads in a vibro-fluidized system. The particles are paramagnetic spheres that add around a fixed magnetic seed inside a granular gas of glass beads. We observe that under appropriate physical conditions symmetric and asymmetric cluster configurations are created and, as the number of particles increases, the aggregation time obeys a power law. We use an ensemble statistics to evaluate the free-energies and entropies landscapes of the granular clusters. The correspondence between such landscapes shows that, even if the system is of macroscopic scale and not in strict equilibrium, our approach to understand the relationship between the cluster structures and the interactions that create them is reliable.