Effect of three-body interactions on the phase behavior of charge-stabilized colloidal suspensions

We study numerically the effect of attractive triplet interactions on the phase behavior of suspensions of highly charged colloidal particles at low salinity. In our computer simulations, we employ the pair and triplet potentials that were obtained from a numerical Poisson-Boltzmann study [Phys. Rev. E 66, 011402 (2002)]]. On the basis of free energy calculations, we determine the phase diagram of an aqueous suspension of identical spheres of diameter sigma=32 nm and charge Z=80 as a function of colloid concentration and salinity, both for the purely pairwise additive system and for the system with pair and triplet interactions. The main effect of including the triplet interactions is a destabilization of the body-centered-cubic (bcc) crystal phase in favor of the face-centered-cubic (fcc) crystal phase. As a consequence the phase diagram features the coexistence of a rather dilute fluid with an almost-close-packed fcc phase at low salinity and bcc-fcc coexistence with a big density jump at intermediate salinity. The triplet attractions do not affect the phase behavior at sufficiently high salinity; under these conditions the system is well described by the pairwise potential.