A large body of experimental work has established that athermal colloid/polymer mixtures undergo a sequence of transitions from a disordered fluid state to a colloidal crystal to a second disordered phase with increasing polymer concentration. These transitions are driven by polymer-mediated interparticle attraction, which is a function of both the polymer density and size. It has been posited that the disordered state at high polymer density is a consequence of strong interparticle attractions that kinetically inhibit the formation of the colloidal crystal, i.e., the formation of a non-equilibrium gel phase interferes with crystallization. Here we use molecular dynamics simulations and density functional theory on polymers and nanoparticles (NPs) of comparable size and show that the crystal-disordered phase coexistence at high polymer density for sufficiently long chains corresponds to an equilibrium thermodynamic phase transition. While the crystal is, indeed, stabilized at intermediate polymer density by polymer-induced intercolloid attractions, it is destabilized at higher densities because long chains lose significant configurational entropy when they are forced to occupy all of the crystal voids. Our results are in quantitative agreement with existing experimental data and show that, at least in the nanoparticle limit of sufficiently small colloidal particles, the crystal phase only has a modest range of thermodynamic stability.npj Computational Materials (2017) 3:3 ; doi:10.1038/s41524-016-0005-8
INTRODUCTIONOver five decades of experimental, theoretical and simulation work 1-13 have led to the knowledge that polymer-induced depletion attractions between colloids in solution can cause them to crystallize even at low polymer concentrations. For large colloids, both face centered cubic (FCC) and hexagonally closed packed (HCP) structures tend to form, since the free energy difference between these two polymorphs is very small. 14 In the nanoparticle limit, when the colloids and chains become comparable in size, we 13 recently showed that colloids preferentially crystallize into the HCP structure in the presence of a dilute concentration of long enough polymers, in preference to the FCC morphology. This preference is driven by the difference in the free energy cost of confining individual polymer chains within the voids of the two different colloidal crystal morphologies. Specifically, polymers prefer to occupy the octahedral voids (OV), which are~6 times larger in volume than the tetrahedral voids (TVs). The HCP structure features connected OVs, while in the FCC OVs are completely surrounded by TVs. Therefore, because a sufficiently long polymer chain must spread across multiple connected voids, it has significantly lower free energy in a HCP crystal than in the FCC analog. Further increasing the polymer density ϕ p , however, experimentally leads to disordered structures, and the formation of apparently kinetically arrested states (diffusion limited or reaction limited aggregations) that has generally been...