Analytical phase diagrams for colloids and non-adsorbing polymer

We review the free-volume theory (FVT) of Lekkerkerker et al. [Europhys. Lett. 20 (1992) 5591 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) 755] 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 H 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 approximate to xi, where the De Gennes blob size (correlation length) xi scales as xi similar to phi(-gamma), with gamma = 0.77 for good solvents and I for a theta solvent. In this limit Pi = Pi(sd) similar to phi(3 gamma). We now apply the following additional modifications: iv) Pi = Pi(0) + Pi(sd), where Pi(sd) = A phi(3 gamma); the prefactor A is known from renormalization group theory. This simple additivity describes the crossover for the osmotic pressure from the dilute limit to the semidilute limit excellently. (v) delta(-2) = delta(-2)(0) + xi(-2), where delta(0)approximate to R is the dilute limit for the depletion thickness delta. This equation describes the crossover in length scales from delta(0) (dilute) to xi(semidilute). 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 R diagrams phi q(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. (C) 2008 Elsevier B.V. All rights reserved.