Statics and dynamics of colloid-polymer mixtures near their critical point of phase separation: A computer simulation study of a continuous Asakura–Oosawa model

Phase behavior and dynamics in a colloid-polymer mixture under spherical confinement

From studies via molecular dynamics simulations, we report results on structure and dynamics in mixtures of active colloids and passive polymers that are confined inside a spherical container with a repulsive boundary. All interactions in the fully passive limit are chosen in such a way that in equilibrium coexistence between colloid-rich and polymer-rich phases occurs. For most part of the studies the chosen compositions give rise to Janus-like structure: nearly one side of the sphere is occupied by the colloids and the rest by the polymers. This partially wet situation mimics approximately a neutral wall in the fully passive scenario. Following the introduction of a velocity-aligning activity to the colloids, the shape of the polymer-rich domain changes to that of an ellipsoid, around the long axis of which the colloid-rich domain attains a macroscopic angular momentum. In the steady state, the orientation of this axis evolves via diffusion, magnitude of which depends upon the strength of activity, but only weakly.

Cross-Correlation Corrected Friction in Generalized Langevin Models: Application to the continuous Asakura-Oosawa Model

In a previous study we proposed a method to parameterize isotropic, configuration independent, non-Markovian generalized Langevin thermostats to achieve dynamic consistency in coarse-grained models. In the current study, by applying the same strategy, we develop coarse-grained implicit solvent models for the continuous Asakura-Oosawa model, which under certain conditions allows to develop very accurate coarse-grained potentials. By developing coarse-grained models for different reference systems with varying parameters, we test the broader applicability of the proposed procedure and demonstrate the relevance of accurate coarse-grained potentials in bottom-up derived dissipative models. We study how different system parameters affect the dynamic representability of the coarse-grained models. In particular we find that the quality of the coarse-grained potential is crucial to correctly model the backscattering effect due to collisions on the coarse-grained scale. In the dynamics of colloid suspensions the hydrodynamic interactions affect the long-time scale dynamics by solvent mediated momentum transfer. These interactions are not explicitly modeled in the presented coarse-grained models, which poses some limitations to the proposed coarse-graining scheme. The Asakura-Oosawa model allows a tuning of system parameters, to gain an improved understanding of these limitations. We also propose three new iterative optimization schemes to fine tune the generalized Langevin thermostat to exactly match the reference velocity-autocorrelation function.

Critical behavior in active lattice models of motility-induced phase separation

Lattice models allow for a computationally efficient investigation of motility-induced phase separation (MIPS) compared to off-lattice systems. Simulations are less demanding, and thus, bigger systems can be accessed with higher accuracy and better statistics. In equilibrium, lattice and off-lattice models with comparable interactions belong to the same universality class. Whether concepts of universality also hold for active particles is still a controversial and open question. Here, we examine two recently proposed active lattice systems that undergo MIPS and investigate numerically their critical behavior. In particular, we examine the claim that these systems and MIPS in general belong to the Ising universality class. We also take a more detailed look on the influence and role of rotational diffusion and active velocity in these systems.

q-Independent Slow Dynamics in Atomic and Molecular Systems

Investigating million-atom systems for very long simulation times, we demonstrate that the collective density-density correlation time (τ_{α}) in simulated supercooled water and silica becomes wave-vector independent (q^{0}) when the probing wavelength is several times larger than the interparticle distance. The q independence of the collective density-density correlation functions, a feature clearly observed in light-scattering studies of some soft-matter systems, is thus a genuine feature of many (but not all) slow-dynamics systems, either atomic, molecular, or colloidal. Indeed, we show that when the dynamics of the density fluctuations includes particle-type diffusion, as in the case of the Lennard-Jones binary-mixture model, the q^{0} regime does not set in and the relaxation time continues to scale as τ_{α}∼q^{-2} even at small q.

Micro-heterogeneity versus clustering in binary mixtures of ethanol with water or alkanes

Snapshots of the difference in complex disorder, with analogy with direct (ethanol–water) and inverse (ethanol–alkanes) emulsions.

How soft repulsion enhances the depletion mechanism

We investigate binary mixtures of large colloids interacting through soft potentials with small, ideal depletants. We show that softness has a dramatic effect on the resulting colloid-colloid effective potential when the depletant-to-colloid size ratio q is small, with significant consequences on the colloidal phase behaviour. We provide an exact relationship that allows us to obtain the effective pair potential for any type of colloid-depletant interaction in the case of ideal depletants, without having to rely on complicated and expensive full-mixture simulations. We also show that soft repulsion among depletants further enhances the tendency of colloids to aggregate. Our theoretical and numerical results demonstrate that--in the limit of small q--soft mixtures cannot be mapped onto hard systems and hence soft depletion is not a mere extension of the widely used Asakura-Oosawa potential.

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)].

Phase behavior of active swimmers in depletants: molecular dynamics and integral equation theory

We study the structure and phase behavior of a binary mixture where one of the components is self-propelling in nature. The interparticle interactions in the system are taken from the Asakura-Oosawa model for colloid-polymer mixtures for which the phase diagram is known. In the current model version, the colloid particles are made active using the Vicsek model for self-propelling particles. The resultant active system is studied by molecular dynamics methods and integral equation theory. Both methods produce results consistent with each other and demonstrate that the Vicsek model-based activity facilitates phase separation, thus, broadening the coexistence region.

Hydrodynamic mechanisms of spinodal decomposition in confined colloid-polymer mixtures: a multiparticle collision dynamics study

A multiscale model for a colloid-polymer mixture is developed. The colloids are described as point particles interacting with each other and with the polymers with strongly repulsive potentials, while polymers interact with each other with a softer potential. The fluid in the suspension is taken into account by the multiparticle collision dynamics method (MPC). Considering a slit geometry where the suspension is confined between parallel repulsive walls, different possibilities for the hydrodynamic boundary conditions (b.c.) at the walls (slip versus stick) are treated. Quenching experiments are considered, where the system volume is suddenly reduced (keeping the density of the solvent fluid constant, while the colloid and polymer particle numbers are kept constant) and thus an initially homogeneous system is quenched deeply into the miscibility gap, where it is unstable. For various relative concentrations of colloids and polymers, the time evolution of the growing colloid-rich and polymer-rich domains are studied by molecular dynamics simulation, taking hydrodynamic effects mediated by the solvent into account via MPC. It is found that the domain size [script-l](d)(t) grows with time t as [script-l](d)(t) [proportionality] t(1/3) for stick and (at late stages) as [script-l](d)(t) [proportionality] t(2/3) for slip b.c., while break-up of percolating structures can cause a transient "arrest" of growth. While these findings apply for films that are 5-10 colloid diameters wide, for ultrathin films (1.5 colloid diameters wide) a regime with [script-l](d)(t) [proportionality] t(1/2) is also identified for rather shallow quenches.

Controlling the wetting properties of the Asakura-Oosawa model and applications to spherical confinement

We demonstrate for the Asakura-Oosawa model and an extension of this model that uses continuous rather than hard potentials, how wetting properties at walls can be easily controlled. By increasing the interaction range of the repulsive wall potential acting on the colloids (while keeping the polymer-wall interactions constant) polymers begin to substitute colloids at walls and the system can be driven from complete wetting of colloids via partial wetting to complete wetting of polymers. As an application, we discuss the morphology and wetting behavior of colloid-polymer mixtures in spherical confinement. We apply the recently developed 'ensemble switch method' where the Hamiltonian is extended to a combination of a system with walls and of a system without walls to calculate the surface excess free energies of colloid-rich and polymer-rich phases. The contact angle then is inferred from Young's equation.

Colloids and polymers in random colloidal matrices: demixing under good-solvent conditions

We consider a simplified coarse-grained model for colloid-polymer mixtures, in which polymers are represented as monoatomic molecules interacting by means of pair potentials. We use it to study polymer-colloid segregation in the presence of a quenched matrix of colloidal hard spheres. We fix the polymer-to-colloid size ratio to 0.8 and consider matrices such that the fraction f of the volume that is not accessible to the colloids due to the matrix is equal to 40%. As in the Asakura-Oosawa-Vrij (AOV) case, we find that binodal curves in the polymer and colloid volume-fraction plane have a small dependence on disorder. As for the position of the critical point, the behavior differs from that observed in the AOV case: While the critical colloid volume fraction is essentially the same in the bulk and in the presence of the matrix, the polymer volume fraction at criticality increases as f increases. At variance with the AOV case, no capillary colloid condensation or evaporation is generically observed.

Spinodal decomposition of polymer solutions: molecular dynamics simulations of the two-dimensional case

As a generic model system for phase separation in polymer solutions, a coarse-grained model for hexadecane/carbon dioxide mixtures has been studied in two-dimensional geometry. Both the phase diagram in equilibrium (obtained from a finite size scaling analysis of Monte Carlo data) and the kinetics of state changes caused by pressure jumps (studied by large scale molecular dynamics simulations) are presented. The results are compared to previous work where the same model was studied in three-dimensional geometry and under confinement in slit geometry. For deep quenches the characteristic length scale ℓ(t) of the formed domains grows with time t according to a power law close to [Formula: see text]. Since in this problem both the polymer density ρ(p) and the solvent density ρ(s) matter, the time evolution of the density distribution P(L)(ρ(p),ρ(s),t) in L × L subboxes of the system is also analyzed. It is found that in the first stage of phase separation the system separates locally into low density carbon dioxide regions that contain no polymers and regions of high density polymer melt that are supersaturated with this solvent. The further coarsening proceeds via the growth of domains of rather irregular shapes. A brief comparison of our findings with results of other models is given.

Crowding of polymer coils and demixing in nanoparticle-polymer mixtures

The Asakura-Oosawa-Vrij (AOV) model of colloid-polymer mixtures idealises nonadsorbing polymers as effective spheres that are fixed in size and impenetrable to hard particles. Real polymer coils, however, are intrinsically polydisperse in size (radius of gyration) and may be penetrated by smaller particles. Crowding by nanoparticles can affect the size distribution of polymer coils, thereby modifying effective depletion interactions and thermodynamic stability. To analyse the influence of crowding on polymer conformations and demixing phase behaviour, we adapt the AOV model to mixtures of nanoparticles and ideal, penetrable polymer coils that can vary in size. We perform Gibbs ensemble Monte Carlo simulations, including trial nanoparticle-polymer overlaps and variations in the radius of gyration. Results are compared with predictions of free-volume theory. Simulation and theory consistently predict that ideal polymers are compressed by nanoparticles, and that compressibility and penetrability stabilise nanoparticle-polymer mixtures.

Confined binary two-dimensional colloidal crystals: Monte Carlo simulation of crack formation

Binary mixtures (A, B) of colloidal particles of different sizes in two dimensions may form crystals with square lattice structure (the A-particles occupying the white sites and the B-particles the black sites of a checkerboard). Confining such a system by two parallel 'walls' a distance D apart, long-range order in the direction parallel to the walls is stabilized by 'corrugated walls' that are commensurate with the lattice structure but destabilized by structureless 'hard walls', even if there is no misfit between the strip width D and the crystal lattice spacing. The crossover to quasi-one-dimensional behavior is studied by Monte Carlo simulations, analyzing Lindemann parameters and displacement correlation functions. When D is reduced and thus a misfit created, the stress in the crystal increases up to a critical value, at which the stress jumps to much smaller values due to the formation of an (almost periodic) crack pattern. These cracks typically have a width of several particle diameters, and are mostly disordered, although sometimes small domains with hexagonal order can be identified. At very large misfits, glass-like structures appear. We discuss various methods to characterize order and disorder in such systems.

A combined molecular dynamics and Monte Carlo study of the approach towards phase separation in colloid-polymer mixtures

A coarse-grained model for colloid-polymer mixtures is investigated where both colloids and polymer coils are represented as point-like particles interacting with spherically symmetric effective potentials. Colloid-colloid and colloid-polymer interactions are described by Weeks-Chandler-Andersen potentials, while the polymer-polymer interaction is very soft, of strength k(B)T/2 for maximum polymer-polymer overlap. This model can be efficiently simulated both by Monte Carlo and molecular dynamics methods, and its phase diagram closely resembles that of the well-known Asakura-Oosawa model. The static and dynamic properties of the model are presented for systems at critical colloid density, varying the polymer density in the one-phase region. Applying Lees-Edwards boundary conditions, colloid-polymer mixtures exposed to shear deformation are considered, and the resulting anisotropy of correlations is studied. Whereas for the considered shear rate, [Formula: see text], radial distribution functions and static structure factors indicate only small structural changes under shear, an appropriate projection of these correlation functions onto spherical harmonics is presented that allows us to directly quantify the structural anisotropies. However, the considered shear rate is probably not high enough to see anisotropies in static structure factors at small wavenumbers that have been predicted by Onuki and Kawasaki (1979 Ann. Phys. 121 456) for the critical behavior of systems under shear. The anomalous dependence of the polymer's self-diffusion constant on polymer density is referred to the clustering of the colloid particles when approaching the critical point.