Computer simulations of phase transitions of bulk and confined colloid–polymer systems

Grand canonical Monte Carlo, histogram reweighting and finite-size scaling methods are used to determine the phase transitions of bulk (three-dimensional) and confined (quasi-two-dimensional) neutral colloid–polymer systems. The colloids are modeled as hard spheres and the polymer molecules as hard chains, and the only attractive forces are effective ones induced by depletion effects. In contrast to the predictions of mean field and other approximate theories, the nature of the coexistence phases is found to not depend solely on the polymer-to-colloid size ratio, q, but on the colloid diameter, the polymer radius of gyration, and the polymer monomer size. The threshold values of q for the onset of liquid–liquid phase separation differ significantly from earlier predictions, and depend strongly on the dimensionality of space. Extrapolation to the “protein limit” of very small colloid and very long polymer indicates that immiscibility persists at this limit in three dimensions, while it does not always do so for confined systems.