We report constant-pressure molecular dynamics (MD) simulations of binary ${\mathrm{C}}_{60}∕{C}_{n}$ fullerene-mixtures ($n=70$, 76, 84, 96) modeled in terms of a spherically symmetric two-body potential. By starting from a liquid configuration of the system, we cool mixtures down to freezing and beyond, until room temperature is reached, in order to verify the formation of solid solutions, namely, of configurations characterized by a unique crystalline lattice whose sites are randomly occupied by the two component fullerene species. We first explore the entire concentration range of the ${\mathrm{C}}_{60x}∕{\mathrm{C}}_{70(1\ensuremath{-}x)}(0lxl1)$ mixture and find fairly good agreement with experimental data, exhibiting partial reciprocal solubility of the two components into each other with an immiscibility gap at intermediate compositions. In fact, the system we simulate forms substitutional solid solutions over a wide range of concentrations except for $0.3\ensuremath{\leqslant}x\ensuremath{\leqslant}0.5$; over such an interval, it turns out that the initially liquid mixture can be supercooled down to relatively low temperatures, until eventually a glassy phase is formed. The study is then extended to fullerene mixtures of molecular diameter ratio $\ensuremath{\alpha}={\ensuremath{\sigma}}_{{\mathrm{C}}_{60}}∕{\ensuremath{\sigma}}_{{\mathrm{C}}_{n}}$ smaller than in ${\mathrm{C}}_{60}∕{\mathrm{C}}_{70}$ (where $\ensuremath{\alpha}=0.93$), as is the case for ${\mathrm{C}}_{60}∕{\mathrm{C}}_{76}$ $(\ensuremath{\alpha}=0.89)$, ${\mathrm{C}}_{60}∕{\mathrm{C}}_{84}(\ensuremath{\alpha}=0.85)$ and ${\mathrm{C}}_{60}∕{\mathrm{C}}_{96}$ $(\ensuremath{\alpha}=0.79)$. The effect of the size mismatch between the two species is dramatic: The solid immiscibility region rapidly expands even upon a tiny reduction of $\ensuremath{\alpha}$, with formation of an amorphous phase at sufficiently low temperature, as found for the ${\mathrm{C}}_{60}∕{\mathrm{C}}_{70}$ mixture. For the smallest $\ensuremath{\alpha}({\mathrm{C}}_{60}∕{\mathrm{C}}_{96})$ cocrystallization of the two components turns out to be forbidden over the whole concentration axis. A mapping of the MD evidences of the fullerene mixtures' phase behavior onto the phase diagram of binary hard-sphere mixtures (determined by other authors) turns out to be worthwhile and enlightnening. In particular, size ratio effects and the onset of glassy phases emerge in qualitative good agreement with such studies, and with results of phase coexistence calculations in model binary colloidal systems.
Read full abstract