Metastable states and the kinetics of colloid phase separation

Nucleation of melt: From fundamentals to dispersed systems

The most evident aspects of a first order transition of a system from an old to a new phase, are the presence of a discontinuity at the interface between both phases and the thermal effects related to the latent heat exchanged with the surrounding environment. These effects are the result of a sequence of events promoted by thermodynamic conditions persisting over the equilibrium in a metastable state. The breakdown of metastability is promoted by infinitesimal energy fluctuations resulting in the germination of clusters of the new phase that can grow to a critical size (nucleus) and then develop or vanish. Examples of these sequences are common in various technological fields such as combustion, food processing, pharmaceutical manufacturing, condensation, and phase change heat transfer, etc. This work aims to highlight a logical path that leads the readers from the fundamental phenomenology to the most intricated aspects of the nucleation within dispersed systems such as oil-in-water emulsions. Differences between the homogeneous and heterogeneous mechanisms are, under the light of the Classical Nucleation Theory (CNT), presented in bulk and confined systems until defining a minimum confinement size. By collecting insights coming from a rich scientific literature mostly focused on the stability of emulsified systems, the discussion is then on the aspects related to the surface related mechanisms. Two main aspects are then considered: a) the wettability of the nucleating cluster by the surrounding melt; b) the affinity between the adsorbed layer, where a surfactant is located, and the oil melt phase (mainly n-alkanes and triacylglycerols with different moieties). In cases where nucleation is dominating over the dewetting of the nucleus, the contact angle can be considered as a constant value. The affinity in terms of molecular features between the surfactant and the oil phase can promote the template effect. Several factors seem to play a role in this interaction such as the thermal characteristics of the surfactant and comparable dimensions between the molecule (or fractions) of the dispersed compound and the tail of the surfactant.

Effects of non-pairwise repulsion on nanoparticle assembly

Electrostatic interactions provide a convenient way to modulate interactions between nanoparticles, colloids, and biomolecules because they can be adjusted by the solution pH or salt concentration. While the presence of salt provides an easy method to control the net interparticle interaction, the nonlinearities arising from electrostatic screening make it difficult to quantify the strength of the interaction. In particular, when charged particles assemble into clusters or aggregates, nonlinear effects render the interactions strongly non-pairwise. Here, we report Brownian dynamics simulations to investigate the effect that the non-pairwise nature of electrostatic interactions has on nanoparticle assembly. We compare these simulations to a system in which the electrostatics are modeled by a strictly pairwise Yukawa potential. We find that both systems show a narrow range in parameter space where the particles form well-ordered crystals. Bordering this range are regions where the net interactions are too weak to stabilize aggregated structures or strong enough that the system becomes kinetically trapped in a gel. The non-pairwise potential differs from the pairwise system in the appearance of an amorphous state for strongly charged particles. This state appears because the many-body electrostatic interactions limit the maximum density achievable in an assembly.

Vitrification and gelation in sticky spheres

Glasses and gels are the two dynamically arrested, disordered states of matter. Despite their importance, their similarities and differences remain elusive, especially at high density, where until now it has been impossible to distinguish them. We identify dynamical and structural signatures which distinguish the gel and glass transitions in a colloidal model system of hard and "sticky" spheres. It has been suggested that "spinodal" gelation is initiated by gas-liquid viscoelastic phase separation to a bicontinuous network and the resulting densification leads to vitrification of the colloid-rich phase, but whether this phase has sufficient density for arrest is unclear [M. A. Miller and D. Frenkel, Phys. Rev. Lett. 90, 135702 (2003) and P. J. Lu et al., Nature 435, 499-504 (2008)]. Moreover alternative mechanisms for arrest involving percolation have been proposed [A. P. R. Eberle et al., Phys. Rev. Lett. 106, 105704 (2011)]. Here we resolve these outstanding questions, beginning by determining the phase diagram. This, along with demonstrating that percolation plays no role in controlling the dynamics of our system, enables us to confirm spinodal decomposition as the mechanism for gelation. We are then able to show that gels can be formed even at much higher densities than previously supposed, at least to a volume fraction of ϕ = 0.59. Far from being networks, these gels apparently resemble glasses but are still clearly distinguished by the "discontinuous" nature of the transition and the resulting rapid solidification, which leads to the formation of inhomogeneous (with small voids) and far-from-equilibrium local structures. This is markedly different from the glass transition, whose continuous nature leads to the formation of homogeneous and locally equilibrated structures. We further reveal that the onset of the attractive glass transition in the form of a supercooled liquid is in fact interrupted by gelation. Our findings provide a general thermodynamic, dynamic, and structural basis upon which we can distinguish gelation from vitrification.

Formation of porous crystals via viscoelastic phase separation

Viscoelastic phase separation of colloidal suspensions can be interrupted to form gels either by glass transition or by crystallization. With a new confocal microscopy protocol, we follow the entire kinetics of phase separation, from homogeneous phase to different arrested states. For the first time in experiments, our results unveil a novel crystallization pathway to sponge-like porous crystal structures. In the early stages, we show that nucleation requires a structural reorganization of the liquid phase, called stress-driven ageing. Once nucleation starts, we observe that crystallization follows three different routes: direct crystallization of the liquid phase, the Bergeron process, and Ostwald ripening. Nucleation starts inside the reorganized network, but crystals grow past it by direct condensation of the gas phase on their surface, driving liquid evaporation, and producing a network structure different from the original phase separation pattern. We argue that similar crystal-gel states can be formed in monatomic and molecular systems if the liquid phase is slow enough to induce viscoelastic phase separation, but fast enough to prevent immediate vitrification. This provides a novel pathway to form nanoporous crystals of metals and semiconductors without dealloying, which may be important for catalytic, optical, sensing, and filtration applications.

The relationship between local density and bond-orientational order during crystallization of the Gaussian core model

Whether nucleation is triggered by density or by bond-orientational order is one of the most hotly debated issues in recent investigations of the crystallization process. Here, we present a numerical study of the relationship between them for soft particles within the isothermal-isobaric ensemble. We compress the system and thus obtain the fluid-solid transition. By investigating locally dense-packed particles and particles with a relatively high bond-orientational order in the compressing process, we find a sharp increase of the spatial correlations for both densely packed particles and highly bond-orientational ordered particles at the phase transition point, which provide new characterization methods for the liquid-crystal transition. We also find that it is the bond-orientational order rather than density that triggers the nucleation process. The relationship between the local density and the bond-orientational order parameter is strongly affected by the characterization methods used. The local bond order parameter (q6) shows clear correlation with the local density (ρ) in the fluid stage, while the coarse-grained form (q[combining macron]6) does not correlate with ρ at all, owing to the comparable spatial scales of q6 and ρ. Nevertheless, q[combining macron]6 shows an obvious advantage in distinguishing between solid and liquid particles in our work. These results may elevate our understanding of the mechanism of the crystallization process.

Crystallization and arrest mechanisms of model colloids

We performed dynamic simulations of spheres with short-range attractive interactions for many values of interaction strength and range. Fast crystallization occurs in a localized region of this parameter space, but the character of crystallization pathways is not uniform within this region. Pathways range from one-step, in which a crystal nucleates directly from a gas, to two-step, in which substantial liquid-like clusters form and only subsequently become crystalline. Crystallization can fail because of slow nucleation from either gas or liquid, or because of dynamic arrest caused by strong interactions. Arrested states are characterized by the formation of networks of face-sharing tetrahedra that can be detected by a local common neighbor analysis.

Optimization of crystal nucleation close to a metastable fluid-fluid phase transition

The presence of a metastable fluid-fluid critical point is thought to dramatically influence the crystallization pathway, increasing the nucleation rate by many orders of magnitude over the predictions of classical nucleation theory. We use molecular dynamics simulations to study the kinetics of crystallization in the vicinity of this metastable critical point and throughout the metastable fluid-fluid phase diagram. To quantitatively understand how the fluid-fluid phase separation affects the crystal nucleation, we evaluate accurately the kinetics and reconstruct the thermodynamic free-energy landscape of crystal formation. Contrary to expectations, we find no special advantage of the proximity of the metastable critical point on the crystallization rates. However, we find that the ultrafast formation of a dense liquid phase causes the crystallization to accelerate both near the metastable critical point and almost everywhere below the fluid-fluid spinodal line. These results unveil three different scenarios for crystallization that could guide the optimization of the process in experiments

Molecular dynamics study of coagulation in silica-nanocolloid-water-NaCl systems based on the atomistic model

In the present work, large scale molecular dynamics (MD) simulations of nanocolloidal silica in aqueous NaCl solutions were performed using a fully atomistic model to study the microscopic structures and dynamics of the systems that lead to aggregation or gelation. Our attention is focused on the self-organizations that occur in the structures of the colloidal silica and water for various concentrations of NaCl. As the salt concentration increased, coagulation developed through the direct bonding of SiO4 units. The trend was explained by the systematic changes in the pair correlation functions related to the barrier height in the potential of mean force [J. G. Kirkwood, J. Chem. Phys., 1935, 3, 300]. Network structures of silica were visualised, and their fractal dimensions were examined by computing the running coordination numbers of Si-Si pairs and also by the analysis of two dimensional images. The calculated dimension by the former method was comparable to the experimental observations for the aggregation of silica colloids, and at longer length scales, super-aggregation was evident in the gelation process. The result from the 2D images is found to be insensitive to the differences in the structure. Clear changes in both the structure and mobility of the water were observed as the NaCl concentration increased, suggesting the importance of the solvent structures to these processes, although such a feature is lacking in the conventional models and most simulations of colloids.

Phase diagrams and kinetics of phase transitions in protein solutions

The phase behavior of proteins is of interest for fundamental and practical reasons. The nucleation of new phases is one of the last major unresolved problems of nature. The formation of protein condensed phases (crystals, polymers, and other solid aggregates, as well as dense liquids and gels) underlies pathological conditions, plays a crucial role in the biological function of the respective protein, or is an essential part of laboratory and industrial processes. In this review, we focus on phase transitions of proteins in their properly folded state. We first summarize the recently acquired understanding of physical processes underlying the phase diagrams of the protein solutions and the thermodynamics of protein phase transitions. Then we review recent findings on the kinetics of nucleation of dense liquid droplets and crystals. We explore the transition from nucleation to spinodal decomposition for liquid-liquid separation and introduce the new concept of solution-to-crystal spinodal. We review the two-step mechanism of protein crystal nucleation, in which mesoscopic metastable protein clusters serve as precursors to the ordered crystal nuclei. The concepts and mechanisms reviewed here provide powerful tools for control of the nucleation process by varying the solution thermodynamic parameters.

Kinetics of nanochain formation in a simplified model of amelogenin biomacromolecules

We show that the kinetics of nanochain formation of amelogenin molecules is well described by a combination of translational and rotational diffusion of a simplified anisotropic bipolar model consisting of hydrophobic spherical colloid particles and a point charge located on each particle surface. The colloid particles interact via a standard depletion attraction whereas the point charges interact through a screened Coulomb repulsion. We study the kinetics via a Brownian dynamics simulation of both translational and rotational motions and show that the anisotropy brought in by the charge dramatically affects the kinetic pathway of cluster formation and our simple model captures the main features of the experimental observations.

Nucleation

Crystallization starts with nucleation and control of nucleation is crucial for the control of the number, size, perfection, polymorphism and other characteristics of crystalline materials. This is particularly true for crystallization in solution, which is an essential part of processes in the chemical and pharmaceutical industries and a major step in physiological and pathological phenomena. There have been significant recent advances in the understanding of the mechanism of nucleation of crystals in solution. The foremost of these are the two-step mechanism of nucleation and the notion of the solution-crystal spinodal. According to the two-step mechanism, the crystalline nucleus appears inside pre-existing metastable clusters of size several hundred nanometers, which consist of dense liquid and are suspended in the solution. While initially proposed for protein crystals, the applicability of this mechanism has been demonstrated for small molecule organic materials, colloids, polymers, and biominerals. This mechanism helps to explain several long-standing puzzles of crystal nucleation in solution: nucleation rates which are many orders of magnitude lower than theoretical predictions, the significance of the dense protein liquid, and others. At high supersaturations typical of most crystallizing systems, the generation of crystal embryos occurs in the spinodal regime, where the nucleation barrier is negligible. The solution-crystal spinodal helps to understand the role of heterogeneous substrates in nucleation and the selection of crystalline polymorphs. Importantly, these ideas provide powerful tools for control of the nucleation process by varying the solution thermodynamic parameters.

Monitoring early stages of silver particle formation in a polymer solution by in situ and time resolved small angle X-ray scattering

Silver particles have been prepared by reduction of silver nitrate with ascorbic acid in acidic aqueous solution containing a low concentration of a commercial polynaphthalene sulfonate polymer (Daxad 19) as dispersant agent. The reduction has been induced and controlled by the slow addition of ascorbic acid at a fixed rate; in this way, we were able to monitor the formation of a silver crystalline colloidal dispersion by in situ and time resolved Small Angle X-ray Scattering measurements. Modeling the scattering intensity with interacting spherical particles in a polymer-Ag like-fractal template allowed us to distinguish different stages involving liquid-like ordered cluster nucleation, cluster growth up to primary particle formation and particle coalescence. Between primary particle formation and particle coalescence, we observed the occurrence of a transient phase of core-shell type structures having primary particles as stable cores in expanding shells built by the organic polymer. We discuss these results in a twofold perspective pertaining both to technology, relative to controlled fabrication of metal nanoparticles and to basic chemical physics, dealing with non standard stepwise crystallization from solutions.

The two-step mechanism of nucleation of crystals in solution

The formation of crystalline nanoparticles starts with nucleation and control of nucleation is crucial for the control of the number, size, perfection, polymorph modification and other characteristics of particles. Recently, there have been significant advances in the understanding of the mechanism of nucleation of crystals in solution. The most significant of these is the two-step mechanism of nucleation, according to which the crystalline nucleus appears inside pre-existing metastable clusters of size several hundred nanometers, which consist of dense liquid and are suspended in the solution. While initially proposed for protein crystals, the applicability of this mechanism has been demonstrated for small-molecule organic and inorganic materials, colloids, and biominerals. This mechanism helps to explain several long-standing puzzles of crystal nucleation in solution: nucleation rates which are many orders of magnitude lower than theoretical predictions, nucleation kinetic dependencies with steady or receding parts at increasing supersaturation, the role of heterogeneous substrates for polymorph selection, the significance of the dense protein liquid, and others. More importantly, this mechanism provides powerful tools for control of the nucleation process by varying the solution thermodynamic parameters so that the volume occupied by the dense liquid shrinks or expands.

Key role of hydrodynamic interactions in colloidal gelation

Colloidal gelation is caused by the formation of a percolated network of colloidal particles suspended in a liquid. Thus far the major transport process leading to gelation has been believed to be the brownian diffusion of particles. Contrary to this common belief, we reveal by numerical simulations that many-body hydrodynamic interactions between colloidal particles also play an essential role in gelation: They significantly promote gelation, or lower the colloid volume fraction threshold for percolation, as compared to their absence. We find that the incompressible nature of a liquid component and the resulting self-organization of hydrodynamic flow with a transverse (rotational) character are responsible for this enhancement of network-forming ability.

Nucleation of crystals from solution: classical and two-step models

Crystallization is vital to many processes occurring in nature and in the chemical, pharmaceutical, and food industries. Notably, crystallization is an attractive isolation step for manufacturing because this single process combines both particle formation and purification. Almost all of the products based on fine chemicals, such as dyes, explosives, and photographic materials, require crystallization in their manufacture, and more than 90% of all pharmaceutical products contain bioactive drug substances and excipients in the crystalline solid state. Hence control over the crystallization process allows manufacturers to obtain products with desired and reproducible properties. We judge the quality of a crystalline product based on four main properties: size, purity, morphology, and crystal structure. The pharmaceutical industry in particular requires production of the desired crystal form (polymorph) to assure the bioavailability and stability of the drug substance. In solution crystallization, nucleation plays a decisive role in determining the crystal structure and size distribution. Therefore, understanding the fundamentals of nucleation is crucial to achieve control over these properties. Because of its analytical simplicity, researchers have widely applied classical nucleation theory to solution crystallization. However, a number of differences between theoretical predictions and experimental results suggest that nucleation of solids from solution does not proceed via the classical pathway but follows more complex routes. In this Account, we discuss the shortcomings of classical nucleation theory and review studies contributing to the development of the modern two-step model. In the two-step model that was initially proposed for protein crystallization, a sufficient-sized cluster of solute molecules forms first, followed by reorganization of that cluster into an ordered structure. In recent experimental and theoretical studies, we and other researchers have demonstrated the applicability of the two-step mechanism to both macromolecules and small organic molecules, suggesting that this mechanism may underlie most crystallization processes from solutions. Because we have observed an increase in the organization time of appropriate lattice structures with greater molecular complexity, we propose that organization is the rate-determining step. Further development of a clearer picture of nucleation may help determine the optimum conditions necessary for the effective organization within the clusters. In addition, greater understanding of these processes may lead to the design of auxiliaries that can increase the rate of nucleation and avoid the formation of undesired solid forms, allowing researchers to obtain the final product in a timely and reproducible manner.

Experimental evidence for two-step nucleation in colloidal crystallization

We investigated the freezing of colloidal spheres in two dimensions with single-particle resolution. Using micron-size, charge-stabilized polystyrene spheres with a temperature-dependent depletion attraction induced by surfactant micelles, we supercooled an initially amorphous (gaslike) system. Particle motions were monitored as crystallization proceeded. At low concentrations, freezing occurred in a single step in a manner consistent with classical nucleation theory. In other samples two-step nucleation was found, in which amorphous clusters grew to approximately 30 particles, then rapidly crystallized. Measured free energies show the role of metastable gas-liquid coexistence, which also enhanced the rate of nucleation following deeper quenches.

Shear effects on crystal nucleation in colloidal suspensions

Extensive two-dimensional Langevin dynamics simulations are used to determine the effect of steady shear flows on the crystal nucleation kinetics of charge stabilized colloids and colloids whose pair potential possess an attractive shallow well of a few k_{B}T 's (attractive colloids). Results show that in both types of systems small amounts of shear speeds up the crystallization process and enhances the quality of the growing crystal significantly. Moderate shear rates, on the other hand, destroy the ordering in the system. The very high shear rate regime where a reentering transition to the ordered state could exist is not considered in this work. In addition to the crystal nucleation phenomena, the analysis of the transport properties and the characterization of the steady state regime under shear are performed.

Fractal aggregates in protein crystal nucleation

Monte Carlo simulations of homogeneous nucleation for a protein model with an exceedingly short-ranged attractive potential yielded a nonconventional crystal nucleation mechanism, which proceeds by the formation of fractal, low-dimensional aggregates followed by a concurrent collapse and increase of the crystallinity of these aggregates to become compact ordered nuclei. This result corroborates a recently proposed two-step mechanism for protein crystal nucleation from solution.

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: II. Nucleation in the metastable liquid immiscibility region

In the second part of our paper, we address crystal nucleation in the metastable liquid miscibility region of eutectic systems that is always present, though experimentally often inaccessible. While this situation resembles the one seen in single component crystal nucleation in the presence of a metastable vapor-liquid critical point addressed in previous works, it is more complex because of the fact that here two crystal phases of significantly different compositions may nucleate. Accordingly, at a fixed temperature below the critical point, six different types of nuclei may form: two liquid-liquid nuclei: two solid-liquid nuclei; and two types of composite nuclei, in which the crystalline core has a liquid "skirt," whose composition falls in between the compositions of the solid and the initial liquid phases, in addition to nuclei with concentric alternating composition shells of prohibitively high free energy. We discuss crystalline phase selection via exploring/identifying the possible pathways for crystal nucleation.

Stability diagram for dense suspensions of model colloidal Al2O3 particles in shear flow

In Al2O3 suspensions, depending on the experimental conditions, very different microstructures can be found, comprising fluidlike suspensions, a repulsive structure, and a clustered microstructure. For technical processing in ceramics, the knowledge of the microstructure is of importance, since it essentially determines the stability of a workpiece to be produced. To enlighten this topic, we investigate these suspensions under shear by means of simulations. We observe cluster formation on two different length scales: the distance of nearest neighbors and on the length scale of the system size. We find that the clustering behavior does not depend on the length scale of observation. If interparticle interactions are not attractive the particles form layers in the shear flow. The results are summarized in a stability diagram.

Gel to glass transition in simulation of a valence-limited colloidal system

We numerically study a simple model for thermoreversible colloidal gelation in which particles can form reversible bonds with a predefined maximum number of neighbors. We focus on three and four maximally coordinated particles, since in these two cases the low valency makes it possible to probe, in equilibrium, slow dynamics down to very low temperatures T. By studying a large region of T and packing fraction phi we are able to estimate both the location of the liquid-gas phase separation spinodal and the locus of dynamic arrest, where the system is trapped in a disordered nonergodic state. We find that there are two distinct arrest lines for the system: a glass line at high packing fraction, and a gel line at low phi and T. The former is rather vertical (phi controlled), while the latter is rather horizontal (T controlled) in the phi-T plane. Dynamics on approaching the glass line along isotherms exhibit a power-law dependence on phi, while dynamics along isochores follow an activated (Arrhenius) dependence. The gel has clearly distinct properties from those of both a repulsive and an attractive glass. A gel to glass crossover occurs in a fairly narrow range in phi along low-T isotherms, seen most strikingly in the behavior of the nonergodicity factor. Interestingly, we detect the presence of anomalous dynamics, such as subdiffusive behavior for the mean squared displacement and logarithmic decay for the density correlation functions in the region where the gel dynamics interferes with the glass dynamics.

Flocculation and percolation in reversible cluster-cluster aggregation

Off-lattice dynamic Monte-Carlo simulations were done of reversible cluster-cluster aggregation for spheres that form rigid bonds at contact. The equilibrium properties were found to be determined by the life time of encounters between two particles (te). te is a function not only of the probability to form or break a bond, but also of the elementary step size of the Brownian motion of the particles. In the flocculation regime the fractal dimension of the clusters is df=2.0 and the size distribution has a power law decay with exponent tau=1.5. At larger values of te transient gels are formed. Close to the percolation threshold the clusters have a fractal dimension df=2.7 and the power law exponent of the size distribution is tau=2.1. The transition between flocculation and percolation occurs at a characteristic weight average aggregation number that decreases with increasing volume fraction.

Thermodynamically stable dispersions induced by depletion interactions

When small particles are added to a colloidal dispersion of large particles, a depletion interaction between large particles occurs because the small ones are depleted from the gaps between the former particles. In the present paper, a cell model is employed to examine the behavior of a dispersion of large particles immersed in an electrolyte solution containing small particles. In this model, each cell consists of one large particle in its center and an associated atmosphere. Double-layer, van der Waals, and depletion interactions, as well as entropic effects, have been taken into account. When the change of the free energy with respect to that of the electrolyte solution is negative (and this happens in most cases), the dispersions of large particles are stable from a thermodynamic point of view. With increasing volume fraction of the small particles, the free energy change becomes more negative. The formation of gels observed experimentally in concentrated emulsions is explained through the formation of a thermodynamically stable dispersion.

Model for reversible colloidal gelation

We report a numerical study, covering a wide range of packing fraction Phi and temperature T, for a system of particles interacting via a square well potential supplemented by an additional constraint on the maximum number n(max) of bonded interactions. We show that, when n(max)<6, the liquid-gas coexistence region shrinks, giving access to regions of low Phi where dynamics can be followed down to low T without an intervening phase separation. We characterize these arrested states at low densities (gel states) in terms of structure and dynamical slowing down, pointing out features which are very different from the standard glassy states observed at high Phi values.

Structure factor scaling in colloidal phase separation

The dynamical scaling hypothesis for the structure factor, S (q) , in depletion-driven colloidal phase separation is studied by carrying out Brownian dynamics simulations. A true dynamical scaling is observed for shallow quenches into the two-phase coexistence region. In such a quench, compact clusters nucleate and grow with time and there is only one characteristic length scale in the system after an initial transient period. Scaling is satisfied beyond this initial period. In contrast, deep quenches lead to fractal cluster growth, and the system is controlled by two characteristic lengths that evolve differently in time [Huang, Oh, and Sorensen (HOS), Phys. Rev. E 57, 875 (1998)]. True dynamical scaling thus cannot be expected to hold. However, an apparent scaling for the structure factor is observed over some period of time when these two characteristic length scales become comparable to each other. We compare our simulation results for the total structure factor to theoretical predictions by HOS by writing it as a product of cluster-cluster and the averaged single-cluster structure factors, each with its own characteristic length.

Effect of bond lifetime on the dynamics of a short-range attractive colloidal system

We perform molecular dynamics simulations of short-range attractive colloid particles modeled by a narrow (3% of the hard sphere diameter) square well potential of unit depth. We compare the dynamics of systems with the same thermodynamics but different bond lifetimes, by adding to the square well potential a thin barrier at the edge of the attractive well. For permanent bonds, the relaxation time tau diverges as the packing fraction phi approaches a threshold related to percolation, while for short-lived bonds, the phi dependence of tau is more typical of a glassy system. At intermediate bond lifetimes, the phi dependence of tau is driven by percolation at low phi , but then crosses over to glassy behavior at higher phi . We also study the wave vector dependence of the percolation dynamics.

Kinetics of phase transformations in depletion-driven colloids

We present results from a detailed numerical study of the kinetics of phase transformations in a model two-dimensional depletion-driven colloidal system. Transition from a single, dispersed phase to a two-phase coexistence of monomers and clusters is obtained as the depth of the interaction potential among the colloidal particles is changed. Increasing the well depth further, fractal clusters are observed in the simulation. These fractal clusters have a hybrid structure in the sense that they show hexagonal closed-packed crystalline ordering at short length scales and a ramified fractal nature at larger length scales. For sufficiently deep potential wells, the diffusion-limited cluster-cluster aggregation model is recovered in terms of the large-scale fractal dimension Df of the clusters, the kinetic exponent z, and the scaling form of the cluster size distribution. For shallower well depths inside the two-phase coexistence region, simulation results for the kinetics of cluster growth are compared with intermediate-stage phase separation in binary mixtures. In the single-phase region, growth kinetics agree well with a mean-field aggregation-fragmentation model of Sorensen, Zhang, and Taylor.

Lateral heterogeneity of photosystems in thylakoid membranes studied by Brownian dynamics simulations

The aggregation and segregation of photosystems in higher plant thylakoid membranes as stromal cation-induced phenomena are studied by the Brownian dynamics method. A theoretical model of photosystems lateral movement within the membrane plane is developed, assuming their pairwise effective potential interaction in aqueous and lipid media and their diffusion. Along with the screened electrostatic repulsive interaction the model accounts for the van der Waals-type, elastic, and lipid-induced attractive forces between photosystems of different sizes and charges. Simulations with a priori estimated parameters demonstrate that all three studied repulsion-attraction alternatives might favor the local segregation of photosystems under physiologically reasonable conditions. However, only the lipid-induced potential combined with the size-corrected screened Coulomb interaction provides the segregated configurations with photosystems II localized in the central part of the grana-size simulation cell and photosystems I occupying its margins, as observed experimentally. Mapping of thermodynamic states reveals that the coexistence curves between isotropic and aggregated phases are the sigmoidlike functions regardless of the effective potential type. It correlates with measurements of the chlorophyll content of thylakoid fragments. Also the universality of the phase curves characterizes the aggregation and segregation of photosystems as order-disorder phase transitions with the Debye radius as a governing parameter.

Brownian dynamics simulations of aging colloidal gels

The aging of colloidal gels is investigated using very long duration Brownian dynamics simulations. The Asakura-Oosawa description of the depletion interaction is used to model a simple colloid polymer mixture. Several regimes are identified during gel formation. The intermediate scattering function displays a double decay characteristic of systems where some kinetic processes are frozen. The beta relaxation at short times is explained in terms of the Krall-Weitz model for the decorrelation due to the elastic modes present. The alpha relaxation at long times is well described by a stretched exponential, showing a wide spectrum of relaxation times for which the q dependence is tau(alpha)=q(-2.2), lower than for diffusion. For the shortest waiting times, a combination of two stretched exponentials is used, suggesting a bimodal distribution. The extracted relaxation times vary with waiting time as tau(alpha)=tau(0.66)(w), more slowly than in the simple aging case. The real space displacements are found to be strongly non-Gaussian, correlated in space and time. We were unable to find clear evidence that the gel aging was driven by internal stresses. Rather, we hypothesize that in this case of weakly interacting gels, the aging behavior arises due to the thermal diffusion of strands, constrained by the percolating network, which ruptures discontinuously. Although the mechanisms differ, the similarity of some of the results to aging of glasses is striking.

Relation between aggregation and phase separation: three-dimensional Monte Carlo simulations

We study phase separation of particles in solution using Monte Carlo simulations of reversible aggregation on a cubic lattice. Two stages of the phase separation can be clearly distinguished: initial random aggregation and subsequent densification. Step one leads to a distribution of fractal aggregates close to the binodal and to a temporary gel for large attractive interaction. Step two leads to isolated spherical dense domains close to the binodal and branched wormlike strands for large attractive interactions. The transition between the two types of structure is gradual and there is no clear feature that shows the existence of a spinodal. The first stage of the phase separation is metastable very close to the binodal or at very large interaction energy. In the latter case, the second step can be viewed as an aging process of the gel formed in the first step.

Liquid-liquid separation in solutions of normal and sickle cell hemoglobin

We show that in solutions of human hemoglobin (Hb)--oxy- and deoxy-Hb A or S--of near-physiological pH, ionic strength, and Hb concentration, liquid-liquid phase separation occurs reversibly and reproducibly at temperatures between 35 and 40 degrees C. In solutions of deoxy-HbS, we demonstrate that the dense liquid droplets facilitate the nucleation of HbS polymers, whose formation is the primary pathogenic event for sickle cell anemia. In view of recent results that shifts of the liquid-liquid separation phase boundary can be achieved by nontoxic additives at molar concentrations up to 30 times lower than the protein concentrations, these findings open new avenues for the inhibition of the HbS polymerization.

Insights into phase transition kinetics from colloid science

Colloids display intriguing transitions between gas, liquid, solid and liquid crystalline phases. Such phase transitions are ubiquitous in nature and have been studied for decades. However, the predictions of phase diagrams are not always realized; systems often become undercooled, supersaturated, or trapped in gel-like states. In many cases the end products strongly depend on the starting position in the phase diagram and discrepancies between predictions and actual observations are due to the intricacies of the dynamics of phase transitions. Colloid science aims to understand the underlying mechanisms of these transitions. Important advances have been made, for example, with new imaging techniques that allow direct observation of individual colloidal particles undergoing phase transitions, revealing some of the secrets of the complex pathways involved.

Classification of ordering kinetics in three-phase systems

Though equations of motion containing transport coefficients are required to quantitatively predict the phase-ordering dynamics of any given system, a great deal can be gleaned just from the shape of the free-energy landscape. We demonstrate how to extract the most information concerning phase-ordering phenomenology from a knowledge of a system's free-energy function, or phase diagram. Many putative pathways to equilibrium can be ruled out on the grounds of the second law of thermodynamics. In some parts of the phase diagram, these considerations are sufficient to completely determine the phase-ordering process without ever having to calculate a transport coefficient, even when three phases are present. The results include a large number of regions of the phase diagram with distinct phase-ordering kinetics, and some surprisingly elaborate routes to the equilibrium state. A process is found whereby a crystalline condensation nucleus becomes coated with a shell of gas, buffering it from a majority metastable liquid phase. Our results, based on thermodynamic arguments, are supported by numerical solution of model B, which describes diffusive phase-ordering kinetics. Some of our predictions are tested against experimental observations of colloid-polymer mixtures, described in more detail in the preceding paper [F. Renth, W. C. K. Poon, and R. M. L. Evans, Phys. Rev. E 64, 031402 (2001)]. A compact notation is developed to represent intricate phase-ordering pathways.