Effects of colloid polydispersity on the phase behavior of colloid-polymer mixtures

Deterministic Formation and Growth of Dendritic Crystals of Attractive Colloids

Dendrites are ubiquitous crystals produced in supersaturated solutions and supercooled melts, but considerably less is known about their formation and growth kinetics. Here, the key factors are explored that dictate dendrite formation and growth, utilizing experimental colloidal models in which the particles act as molecules with Mie potential. Depletion attraction is employed to colloids and manipulate their strength to control supersaturation. Dendrites are predominantly produced under conditions of low supersaturation, where the separation between crystals is large due to slow nucleation. The dendrites do not emerge directly from nuclei. Instead, isotropic grains, initially produced from nuclei, morph into polygons. Arms then sprout from the vertices of these polygons, eventually giving rise to dendrites. Triggering this polygon-to-dendrite transformation requires a high diffusional flux. This necessitates a prolonged diffusion time to maintain a steep concentration gradient in the surrounding environment even after the transformation from circular grains to polygons.

Multiphase Coexistence in Binary Hard Colloidal Mixtures: Predictions from a Simple Algebraic Theory

A general theoretical framework is proposed to quantify the thermodynamic properties of multicomponent hard colloidal mixtures. This framework is used to predict the phase behavior of mixtures of rods with spheres and rods with plates taking into account (liquid) crystal phases of both components. We demonstrate a rich and complex range of phase behaviors featuring a large variety of different multiphase coexistence regions, including two five-phase coexistence regions for hard rod/sphere mixtures, and even a six-phase equilibrium for hard rod/plate dispersions. The various multiphase coexistences featured in a particular mixture are in line with a recently proposed generalized phase rule and can be tuned through subtle variations of the particle shape and size ratio. Our approach qualitatively accounts for certain multiphase equilibria observed in rod/plate mixtures of clay colloids and will be a useful guide in tuning the phase behavior of shape-disperse mixtures in general.

Sedimentation path theory for mass-polydisperse colloidal systems

Both polydispersity and the presence of a gravitational field are inherent to essentially any colloidal experiment. While several theoretical works have focused on the effect of polydispersity on the bulk phase behavior of a colloidal system, little is known about the effect of a gravitational field on a polydisperse colloidal suspension. We extend here the sedimentation path theory to study sedimentation-diffusion-equilibrium of a mass-polydisperse colloidal system: the particles possess different buoyant masses but they are otherwise identical. The model helps to understand the interplay between gravity and polydispersity on sedimentation experiments. Since the theory can be applied to any parent distribution of buoyant masses, it can also be used to study the sedimentation of monodisperse colloidal systems. We find that mass-polydispersity has a strong influence in colloidal systems near density matching for which the bare density of the colloidal particles equals the solvent density. To illustrate the theory, we study crystallization in sedimentation-diffusion-equilibrium of a suspension of mass-polydisperse hard spheres.

Phase stability of colloidal mixtures of spheres and rods

We determined the phase boundaries of aqueous mixtures containing colloidal rod-like fd-viruses and polystyrene spheres using diffusing-wave spectroscopy and compared the results with free volume theory predictions. Excluded volume interactions in mixtures of colloidal rods and spheres lead to mediated depletion interactions. The strength and range of this attractive interaction depend on the concentrations of the particles, the length L and diameter D of the rods, and the radius R of the spheres. At strong enough attraction, this depletion interaction leads to phase separation. We experimentally determined the rod and sphere concentrations where these phase transitions occur by systematically varying the size ratios L/R and D/R and the aspect ratio L/D. This was done by using spheres with different radii and modifying the effective diameter of the rods through either the ionic strength of the buffer or anchoring a polymeric brush to the surface of the rods. The observed phase transitions were from a binary fluid to a colloidal gas/liquid phase coexistence that occurred already at very low concentrations due to the depletion efficiency of highly anisotropic rods. The experimentally measured phase transitions were compared to phase boundaries obtained using free volume theory (FVT), a well established theory for calculating the phase behavior of colloidal particles mixed with depletants. We find good correspondence between the experimental phase transitions and the theoretical FVT model where the excluded volume of the rod-like depletants was explicitly accounted for in both the reservoir and the system.

Polymer nanocomposite capsules formed by droplet extraction: spontaneous stratification and tailored dissolution

We report the formation of polymeric and nanocomposite capsules via droplet solvent extraction, focusing on the interplay between solvent exchange and removal, demixing and directional solidification kinetics. We investigate a model system of sodium poly(styrene sulfonate), NaPSS and silica nanoparticles in aqueous solution, whose phase behaviour is experimentally measured, and examine a series of selective extraction solvents (toluene, butyl acetate, ethyl acetate and methyl ethyl ketone), ranging from 0.04 to 11% v/v water solubility. Tuning the rate of solvent exchange is shown to provide an effective means of decoupling demixing and solidification timescales, and thereby tunes the internal microstructure of the capsule, including hollow, microporous, core-shell, and bicontinuous morphologies. In turn, these determine the capsule dissolution mechanism and kinetics, ranging from single to pulsed release profiles of nanoparticle clusters (at intermediate solubilities), to minimal dissolution (at either extremes). These findings provide facile design and assembly strategies for functional capsules with time-varying release profiles.

Phase behavior of colloid-polymer depletion mixtures with unary or binary depletants

Adding depletants to a colloidal suspension induces an attractive interparticle interaction that can be tuned to obtain desired structures or to probe phase behavior. When the depletant is not uniform in size, however, both the range and strength of the attraction become difficult to predict and hence control. We investigated the effects of depletant bidispersity on the non-equilibrium phase behavior of colloid-polymer mixtures. We added unary or binary mixtures of polystyrene as the depletant to suspensions of charged poly(methyl methacrylate) particles. The structure and dynamics of the particles were compared over three sets of samples with various mixtures of two different polystyrenes whose size varied by an order of magnitude. The structure and dynamics were nearly independent of depletant dispersity if the polymer concentration was represented as a sum of normalized concentrations of each species. Near the transition region between a fluid of clusters and an interconnected gel at intermediate volume fractions, partitioning of polymers in a binary mixture into colloid-rich and polymer-rich phase leads to a slightly different gelation pathway.

Two- and three-phase equilibria of polydisperse Yukawa hard-sphere fluids confined in random porous media: high temperature approximation and scaled particle theory

We have studied the phase behavior of polydisperse Yukawa hard-sphere fluid confined in random porous media using extension and combination of high temperature approximation and scaled particle theory. The porous media are represented by the matrix of randomly placed hard-sphere obstacles. Due to the confinement, polydispersity effects are substantially enhanced. At an intermediate degree of fluid polydispersity and low density of the matrix, we observe two-phase coexistence with two critical points, and cloud and shadow curves forming closed loops of ellipsoidal shape. With the increase of the matrix density and the constant degree of polydispersity, these two critical points merge and disappear, and at lower temperatures the system fractionates into three coexisting phases. A similar phase behavior was observed in the absence of the porous media caused, however, by the increase of the polydispersity.

Phase behaviour of colloids plus weakly adhesive polymers

The phase behaviour of a colloidal dispersion mediated by weakly adhesive polymers is considered. The polymers are depleted but are weakly adhesive and hence comprise a non-zero polymer concentration at the colloid's surface, in contrast to the classical assumption in depletion theories involving a zero polymer concentration at the surface. The theory is composed of a generalized free-volume theory for colloid-polymer mixtures and a self-consistent mean-field theory for polymers at surfaces. It is found that the weak adhesion of the polymers shifts the phase stability of the colloid-polymer mixtures to higher polymer concentrations as compared to assuming a full depletion effect. The predicted phase diagrams employing the new theory are consistent with experiments on mixtures of silica spheres coated with stearyl alcohol and polydimethylsiloxane in cyclohexane and with Monte Carlo simulation results.

Equilibrium fluid-crystal interfacial free energy of bcc-crystallizing aqueous suspensions of polydisperse charged spheres

The interfacial free energy is a central quantity in crystallization from the metastable melt. In suspensions of charged colloidal spheres, nucleation and growth kinetics can be accurately measured from optical experiments. In previous work, from these data effective nonequilibrium values for the interfacial free energy between the emerging bcc nuclei and the adjacent melt in dependence on the chemical potential difference between melt phase and crystal phase were derived using classical nucleation theory (CNT). A strictly linear increase of the interfacial free energy was observed as a function of increased metastability. Here, we further analyze these data for five aqueous suspensions of charged spheres and one binary mixture. We utilize a simple extrapolation scheme and interpret our findings in view of Turnbull's empirical rule. This enables us to present the first systematic experimental estimates for a reduced interfacial free energy, σ(0,bcc), between the bcc-crystal phase and the coexisting equilibrium fluid. Values obtained for σ(0,bcc) are on the order of a few k(B)T. Their values are not correlated to any of the electrostatic interaction parameters but rather show a systematic decrease with increasing size polydispersity and a lower value for the mixture as compared to the pure components. At the same time, σ(0) also shows an approximately linear correlation to the entropy of freezing. The equilibrium interfacial free energy of strictly monodisperse charged spheres may therefore be still greater.

Gelation in mixtures of polymers and bidisperse colloids

We investigated the effects of varying the volume fraction of large particles (r) on the linear rheology and microstructure of mixtures of polymers and bidisperse colloids, in which the ratio of the small and large particle diameters was α=0.31 or α=0.45. Suspensions formulated at a total volume fraction of ϕ_{T}=0.15 and a constant concentration of polymer in the free volume c/c^{*}≈0.7 contained solid-like gels for small r and fluids or fluids of clusters at large r. The solid-like rheology and microstructure of these suspensions changed little with r when r was small, and fluidized only when r>0.8. By contrast, dense suspensions with ϕ_{T}=0.40 and α=0.31 contained solid-like gels at all concentrations of large particles and exhibited only modest rheological and microstructural changes upon varying the volume fraction of large particles. These results suggest that the effect of particle-size dispersity on the properties of colloid-polymer mixtures are asymmetric in particle size and are most pronounced near a gelation boundary.

Crystallization kinetics of colloidal binary mixtures with depletion attraction

In this work the crystallization kinetics of colloidal binary mixtures with attractive interaction potential (Asakura-Oosawa) has been addressed. Parameters such as fraction of crystals, linear crystal dimension and crystal packing have been quantified in order to understand how the crystal formation is driven in terms of the depth of the attractive potential and the composition of the binary mixture (described by the number ratio). It was found that inside the eutectic triangle, crystallization is mainly governed by nucleation and the crystal packing is close to the close-packing of hard spheres. Moving out from the eutectic triangle towards small component results in the crystallization of small spheres. Enrichment of the eutectic mixture with large component results in the crystallization of both large and small spheres, however, the kinetics are completely different from those of the eutectic composition. Crosslinked polystyrene microgels with nearly hard sphere interactions were used as model systems. Attraction was introduced by addition of linear polystyrene. The time evolution of crystallization has been followed by static light scattering.

Crystallization kinetics of colloidal model suspensions: recent achievements and new perspectives

Colloidal model systems allow studying crystallization kinetics under fairly ideal conditions, with rather well-characterized pair interactions and minimized external influences. In complementary approaches experiment, analytic theory and simulation have been employed to study colloidal solidification in great detail. These studies were based on advanced optical methods, careful system characterization and sophisticated numerical methods. Over the last decade, both the effects of the type, strength and range of the pair-interaction between the colloidal particles and those of the colloid-specific polydispersity have been addressed in a quantitative way. Key parameters of crystallization have been derived and compared to those of metal systems. These systematic investigations significantly contributed to an enhanced understanding of the crystallization processes in general. Further, new fundamental questions have arisen and (partially) been solved over the last decade: including, for example, a two-step nucleation mechanism in homogeneous nucleation, choice of the crystallization pathway, or the subtle interplay of boundary conditions in heterogeneous nucleation. On the other hand, via the application of both gradients and external fields the competition between different nucleation and growth modes can be controlled and the resulting microstructure be influenced. The present review attempts to cover the interesting developments that have occurred since the turn of the millennium and to identify important novel trends, with particular focus on experimental aspects.

Phase diagram of mixtures of colloids and polymers in the thermal crossover from good to θ solvent

We determine the phase diagram of mixtures of spherical colloids and neutral nonadsorbing polymers in the thermal crossover region between the θ point and the good-solvent regime. We use the generalized free-volume theory, which takes into account the polymer-concentration dependence of the depletion thickness and of the polymer compressibility. This approach turns out to be quite accurate as long as q = Rg/Rc ≲ 1 (Rg is the radius of gyration of the polymer and Rc is the colloid radius). We find that, close to the θ point, the phase diagram is not very sensitive to solvent quality, while, close to the good-solvent region, changes of the solvent quality modify significantly the position of the critical point and of the binodals. We also analyze the phase behavior of aqueous solutions of charged colloids and polymers, using the approach proposed by Fortini et al. [J. Phys.: Condens. Matter 17, 7783 (2005)].

Spinodal fractionation in a polydisperse square-well fluid

Using kinetic Monte Carlo simulation, we model gas-liquid spinodal decomposition in a size-polydisperse square well fluid, representing a "near-monodisperse" colloidal dispersion. We find that fractionation (demixing) of particle sizes between the phases begins asserting itself shortly after the onset of phase ordering. Strikingly, the direction of size fractionation can be reversed by a seemingly trivial choice between two interparticle potentials which, in the monodisperse case, are identical--we rationalize this in terms of a perturbative, equilibrium theory of polydispersity. Furthermore, our quantitative results show that kinetic Monte Carlo simulation can provide detailed insight into the role of fractionation in real colloidal systems.

Polydispersity effects in colloid-polymer mixtures

We study phase separation and transient gelation experimentally in a mixture consisting of polydisperse colloids (polydispersity: ≈ 6%) and non-adsorbing polymers, where the ratio of the average size of the polymer to that of the colloid is ≈ 0.062. Unlike what has been reported previously for mixtures with somewhat lower colloid polydispersity (≈ 5%), the addition of polymers does not expand the fluid-solid coexistence region. Instead, we find a region of fluid-solid coexistence which has an approximately constant width but an unexpected re-entrant shape. We detect the presence of a metastable gas-liquid binodal, which gives rise to two-stepped crystallization kinetics that can be rationalized as the effect of fractionation. Finally, we find that the separation into multiple coexisting solid phases at high colloid volume fractions predicted by equilibrium statistical mechanics is kinetically suppressed before the system reaches dynamical arrest.

Attractive emulsion droplets probe the phase diagram of jammed granular matter

It remains an open question whether statistical mechanics approaches apply to random packings of athermal particles. Although a jamming phase diagram has recently been proposed for hard spheres with varying friction, here we use a frictionless emulsion system in the presence of depletion forces to sample the available phase space of packing configurations. Using confocal microscopy, we access their packing microstructure and test the theoretical assumptions. As a function of attraction, our packing protocol under gravity leads to well-defined jammed structures in which global density initially increases above random close packing and subsequently decreases monotonically. Microscopically, the fluctuations in parameters describing each particle, such as the coordination number, number of neighbors, and local packing fraction, are for all attractions in excellent agreement with a local stochastic model, indicating that long-range correlations are not important. Furthermore, the distributions of local cell volumes can be collapsed onto a universal curve using the predicted k-gamma distribution, in which the shape parameter k is fixed by the polydispersity while the effect of attraction is captured by rescaling the average cell volume. Within the Edwards statistical mechanics framework, this result measures the decrease in compactivity with global density, which represents a direct experimental test of a jamming phase diagram in athermal systems. The success of these theoretical tools in describing yet another class of materials gives support to the much-debated statistical physics of jammed granular matter.

Analytical phase diagrams for colloids and non-adsorbing polymer

We review the free-volume theory (FVT) of Lekkerkerker et al. [Europhys. Lett. 20 (1992) 559] for the phase behavior of colloids in the presence of non-adsorbing polymer and we extend this theory in several aspects: (i) We take the solvent into account as a separate component and show that the natural thermodynamic parameter for the polymer properties is the insertion work Pi(v), where Pi is the osmotic pressure of the (external) polymer solution and v the volume of a colloid particle. (ii) Curvature effects are included along the lines of Aarts et al. [J. Phys.: Condens. Matt. 14 (2002) 7551] but we find accurate simple power laws which simplify the mathematical procedure considerably. (iii) We find analytical forms for the first, second, and third derivatives of the grand potential, needed for the calculation of the colloid chemical potential, the pressure, gas-liquid critical points and the critical endpoint (cep), where the (stable) critical line ends and then coincides with the triple point. This cep determines the boundary condition for a stable liquid. We first apply these modifications to the so-called colloid limit, where the size ratio q(R)=R/a between the radius of gyration R of the polymer and the particle radius a is small. In this limit the binodal polymer concentrations are below overlap: the depletion thickness delta is nearly equal to R, and Pi can be approximated by the ideal (van't Hoff) law Pi=Pi(0)=phi/N, where phi is the polymer volume fraction and N the number of segments per chain. The results are close to those of the original Lekkerkerker theory. However, our analysis enables very simple analytical expressions for the polymer and colloid concentrations in the critical and triple points and along the binodals as a function of q(R). Also the position of the cep is found analytically. In order to make the model applicable to higher size ratio's q(R) (including the so-called protein limit where q(R)>1) further extensions are needed. We introduce the size ratio q=delta/a, where the depletion thickness delta is no longer of order R. In the protein limit the binodal concentrations are above overlap. In such semidilute solutions delta approximately xi, where the De Gennes blob size (correlation length) xi scales as xi approximately phi(-gamma), with gamma=0.77 for good solvents and gamma=1 for a theta solvent. In this limit Pi=Pi(sd) approximately phi(3gamma). We now apply the following additional modifications: With these latter two modifications we obtain again a fully analytical model with simple equations for critical and triple points as a function of q(R). In the protein limit the binodal polymer concentrations scale as q(R)(1/gamma), and phase diagrams phiq(R)(-1/gamma) versus the colloid concentration eta become universal (i.e., independent of the size ratio q(R)). The predictions of this generalized free-volume theory (GFVT) are in excellent agreement with experiment and with computer simulations, not only for the colloid limit but also for the protein limit (and the crossover between these limits). The q(R)(1/gamma) scaling is accurately reproduced by both simulations and other theoretical models. The liquid window is the region between phi(c) (critical point) and phi(t) (triple point). In terms of the ratio phi(t)/phi(c) the liquid window extends from 1 in the cep (here phi(t)-phi(c)=0) to 2.2 in the protein limit. Hence, the liquid window is narrow: it covers at most a factor 2.2 in (external) polymer concentration.

Phase equilibria in polydisperse nonadditive hard-sphere systems

Colloidal particles naturally exhibit a size polydispersity that can greatly influence their phase behavior in solution. Nonadditive hard-sphere (NAHS) mixtures are simple and well-suited model systems to represent phase transitions in colloid systems. Here, we propose an analytical equation of state (EOS) for NAHS fluid mixtures, which can be straightforwardly applied to polydisperse systems. For positive values of the nonadditivity parameter Delta the model gives accurate predictions of the simulated fluid-fluid coexistence curves and compressibility factors. NPT Monte Carlo simulations of the mixing properties of the NAHS symmetric binary mixture with Delta>0 are reported. It is shown that the enthalpy of mixing is largely positive and overcomes the positive entropy of mixing when the pressure is increased, leading to a fluid-fluid phase transition with a lower critical solution pressure. Phase equilibria in polydisperse systems are predicted with the model by using the density moment formalism [P. Sollich, Adv. Chem. Phys. 116, 265 (2001)]. We present predictions of the cloud and shadow curves for polydisperse NAHS systems composed of monodisperse spheres and polydisperse colloid particles. A fixed nonadditivity parameter Delta > 0 is assumed between the monodisperse and polydisperse spheres, and a Schulz distribution is used to represent the size polydispersity. Polydispersity is found to increase the extent of the immiscibility region. The predicted cloud and shadow curves depend dramatically on the upper cutoff diameter sigmac of the Schulz distribution, and three-phase equilibria can occur for large values of sigmac.

Colloid-polymer mixtures in the protein limit

This review discusses the structure and phase behaviour of mixtures of colloidal particles and non-adsorbing polymers in the protein limit of large polymers and small colloids. The vast majority of work on colloid-polymer mixtures has been concerned with the colloid limit of large colloidal particles and small polymer chains. In this regime, the diameter of the colloidal particles, , is larger than the characteristic size of the polymer-taken as twice their radius of gyration, . The opposite limit, of size ratios , is called the protein limit due to the common practice of adding polymer to protein solutions in order to aid protein crystallisation. Theoretical predictions for systems in the protein limit are considered briefly and then the main focus is on recent experimental studies of mixtures in the protein limit.

Phase behavior of weakly polydisperse sticky hard spheres: Perturbation theory for the Percus-Yevick solution

We study the effects of size polydispersity on the gas-liquid phase behavior of mixtures of sticky hard spheres. To achieve this, the system of coupled quadratic equations for the contact values of the partial cavity functions of the Percus-Yevick solution [R. J. Baxter, J. Chem. Phys. 49, 2770 (1968)] is solved within a perturbation expansion in the polydispersity, i.e., the normalized width of the size distribution. This allows us to make predictions for various thermodynamic quantities which can be tested against numerical simulations and experiments. In particular, we determine the leading order effects of size polydispersity on the cloud curve delimiting the region of two-phase coexistence and on the associated shadow curve; we also study the extent of size fractionation between the coexisting phases. Different choices for the size dependence of the adhesion strengths are examined carefully; the Asakura-Oosawa model [J. Chem. Phys. 22, 1255 (1954)] of a mixture of polydisperse colloids and small polymers is studied as a specific example.

Influence of polydispersity on the critical parameters of an effective-potential model for asymmetric hard-sphere mixtures

We report a Monte Carlo simulation study of the properties of highly asymmetric binary hard-sphere mixtures. This system is treated within an effective fluid approximation in which the large particles interact through a depletion potential [R. Roth, Phys. Rev. E 62 5360 (2000)] designed to capture the effects of a virtual sea of small particles. We generalize this depletion potential to include the effects of explicit size dispersity in the large particles and consider the case in which the particle diameters are distributed according to a Schulz form having a degree of polydispersity 14%. The resulting alteration (with respect to the monodisperse limit) of the metastable fluid-fluid critical point parameters is determined for two values of the ratio of the diameters of the small and large particles: q(triple bond)sigma(s)/(-)sigma(b)=0.1 and q=0.05. We find that the inclusion of polydispersity moves the critical point to lower reservoir volume fractions of the small particles and high volume fractions of the large ones. The estimated critical point parameters are found to be in good agreement with those predicted by a generalized corresponding states argument which provides a link to the known critical adhesion parameter of the adhesive hard-sphere model. Finite-size scaling estimates of the cluster percolation line in the one phase fluid region indicate that inclusion of polydispersity moves the critical point deeper into the percolating regime. This suggests that phase separation is more likely to be preempted by dynamical arrest in polydisperse systems.

Experimental phase diagram of a model colloid-polymer mixture in the protein limit

Sterically stabilized silica nanoparticles were synthesized, and turbidity measurements confirmed that they behaved as hard spheres in cyclohexane. Poly(isoprene) was added to give mixtures in the protein limit with a polymer coil/colloid radius ratio of 4.8. Their phase behavior under good solvent conditions was studied experimentally. The critical colloid volume fraction was phi = 0.13, whereas recent simulations (Bolhuis, P. G.; Meijer, E. J.; Louis, A. A. Phys. Rev. Lett. 2003, 90, 068304) predicted phi = 0.24. This difference is ascribed to the fact that many systems showing good solvent scaling behavior of the polymer still have a Flory-Huggins parameter close to 0.5, for instance, chi = 0.45 in this work.