No Arabic abstract
Monte Carlo simulations of the Asakura-Oosawa (AO) model for colloid-polymer mixtures confined between two parallel repulsive structureless walls are presented and analyzed in the light of current theories on capillary condensation and interface localization transitions. Choosing a polymer to colloid size ratio of q=0.8 and studying ultrathin films in the range of D=3 to D=10 colloid diameters thickness, grand canonical Monte Carlo methods are used; phase transitions are analyzed via finite size scaling, as in previous work on bulk systems and under confinement between identical types of walls. Unlike the latter work, inequivalent walls are used here: while the left wall has a hard-core repulsion for both polymers and colloids, at the right wall an additional square-well repulsion of variable strength acting only on the colloids is present. We study how the phase separation into colloid-rich and colloid-poor phases occurring already in the bulk is modified by such a confinement. When the asymmetry of the wall-colloid interaction increases, the character of the transition smoothly changes from capillary condensation-type to interface localization-type. The critical behavior of these transitions is discussed, as well as the colloid and polymer density profiles across the film in the various phases, and the correlation of interfacial fluctuations in the direction parallel to the confining walls. The experimental observability of these phenomena also is briefly discussed.
We extensively investigated the critical behavior of mixtures of colloids and polymers via the two-component Asakura-Oosawa model and its reduction to a one-component colloidal fluid using accurate theoretical and simulation techniques. In particular the theoretical approach, hierarchical reference theory [Adv. Phys. 44, 211 (1995)], incorporates realistically the effects of long-range fluctuations on phase separation giving exponents which differ strongly from their mean-field values, and are in good agreement with those of the three-dimensional Ising model. Computer simulations combined with finite-size scaling analysis confirm the Ising universality and the accuracy of the theory, although some discrepancy in the location of the critical point between one-component and full-mixture description remains. To assess the limit of the pair-interaction description, we compare one-component and two-component results.
Employing Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene between graphite walls we show that the mass exchange between layers close to the walls is a slow process already in the melt state. For the glass transition of confined polymers this process competes with the slowing down due to packing effects and intramolecular rotation barriers.
We investigated the viscoelastic properties of colloid-polymer mixtures at intermediate colloid volume fraction and varying polymer concentrations, thereby tuning the attractive interactions. Within the examined range of polymer concentrations, the samples ranged from fluids to gels. Already in the liquid phase the viscoelastic properties significantly changed when approaching the gelation boundary, indicating the formation of clusters and transient networks. This is supported by an increasing correlation length of the density fluctuations, observed by static light scattering and microscopy. At the same time, the correlation function determined by dynamic light scattering completely decays, indicating the absence of dynamical arrest. Upon increasing the polymer concentration beyond the gelation boundary, the rheological properties changed qualitatively again, now they are consistent with the formation of colloidal gels. Our experimental results, namely the location of the gelation boundary as well as the elastic (storage) and viscous (loss) moduli, are compared to different theoretical models. These include consideration of the escape time as well as predictions for the viscoelastic moduli based on scaling relations and Mode Coupling Theories (MCT).
We report on the linear viscoelastic properties of mixtures comprising multiarm star (as model soft colloids) and long linear chain homopolymers in a good solvent. In contrast to earlier works, we investigated symmetric mixtures (with a size ratio of 1) and showed that the polymeric and colloidal responses can be decoupled. The adopted experimental protocol involved probing the linear chain dynamics in different star environments. To this end, we studied mixtures with different star mass fraction, which was kept constant while linear chains were added and their entanglement plateau modulus ($G_p$) and terminal relaxation time ($tau_d$) were measured as functions of their concentration. Two distinct scaling regimes were observed for both $G_p$ and $tau_d$: at low linear polymer concentrations, a weak concentration dependence was observed, that became even weaker as the fraction of stars in the mixtures increased into the star glassy regime. On the other hand, at higher linear polymer concentrations, the classical entangled polymer scaling was recovered. Simple scaling arguments show that the threshold crossover concentration between the two regimes corresponds to the maximum osmotic star compression and signals the transition from confined to bulk dynamics. These results provide the needed ingredients to complete the state diagram of soft colloid-polymer mixtures and investigate their dynamics at large polymer-colloid size ratios. They also offer an alternative way to explore aspects of the colloidal glass transition and the polymer dynamics in confinement. Finally, they provide a new avenue to tailor the rheology of soft composites.
As first explained by the classic Asakura-Oosawa (AO) model, effective attractive forces between colloidal particles induced by depletion of nonadsorbing polymers can drive demixing of colloid-polymer mixtures into colloid-rich and colloid-poor phases, with practical relevance for purification of water, stability of foods and pharmaceuticals, and macromolecular crowding in biological cells. By idealizing polymer coils as effective penetrable spheres, the AO model qualitatively captures the influence of polymer depletion on thermodynamic phase behavior of colloidal suspensions. In previous work, we extended the AO model to incorporate aspherical polymer conformations and showed that fluctuating shapes of random-walk coils can significantly modify depletion potentials [W. K. Lim and A. R. Denton, Soft Matter 12, 2247 (2016); J. Chem. Phys. 144, 024904 (2016)]. We further demonstrated that the shapes of polymers in crowded environments depend sensitively on solvent quality [W. J. Davis and A. R. Denton, J. Chem. Phys. 149, 124901 (2018)]. Here we apply Monte Carlo simulation to analyze the influence of solvent quality on depletion potentials in mixtures of hard sphere colloids and nonadsorbing polymer coils, modeled as ellipsoids whose principal radii fluctuate according to random-walk statistics. We consider both self-avoiding and non-self-avoiding random walks, corresponding to polymers in good and theta solvents, respectively. Our simulation results demonstrate that depletion of polymers of equal molecular weight induces much stronger attraction between colloids in good solvents than in theta solvents and confirm that depletion interactions are significantly influenced by aspherical polymer conformations.