By extending the nonequilibrium potential refinement algorithm and lattice switch method to the semigrand ensemble, the semigrand potentials of the fcc and hcp structures of polydisperse hard-sphere crystals are calculated with the bias sampling scheme. The result shows that the fcc structure is more stable than the hcp structure for polydisperse hard-sphere crystals below the terminal polydispersity.
We compute the equilibrium concentration of stacking faults and point defects in polydisperse hard-sphere crystals. We find that, while the concentration of stacking faults remains similar to that of monodisperse hard sphere crystals, the concentration of vacancies decreases by about a factor two. Most strikingly, the concentration of interstitials in the maximally polydisperse crystal may be some six orders of magnitude larger than in a monodisperse crystal. We show that this dramatic increase in interstitial concentration is due to the increased probability of finding small particles and that the small-particle tail of the particle size distribution is crucial for the interstitial concentration in a colloidal crystal.
An approach to obtain the structural properties of additive binary hard-sphere mixtures is presented. Such an approach, which is a nontrivial generalization of the one recently used for monocomponent hard-sphere fluids [S. Pieprzyk, A. C. Branka, and D. M. Heyes, Phys. Rev. E 95, 062104 (2017)], combines accurate molecular-dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation. A comparison of the direct correlation functions obtained with the present scheme with those derived from theoretical results stemming from the Percus-Yevick (PY) closure and the so-called rational-function approximation (RFA) is performed. The density dependence of the leading poles of the Fourier transforms of the total correlation functions and the decay of the pair correlation functions of the mixtures are also addressed and compared to the predictions of the two theoretical approximations. A very good overall agreement between the results of the present scheme and those of the RFA is found, thus suggesting that the latter (which is an improvement over the PY approximation) can safely be used to predict reasonably well the long-range behavior, including the structural crossover, of the correlation functions of additive binary hard-sphere mixtures.
The solid-solid coexistence of a polydisperse hard sphere system is studied by using the Monte Carlo simulation. The results show that for large enough polydispersity the solid-solid coexistence state is more stable than the single-phase solid. The two coexisting solids have different composition distributions but the same crystal structure. Moreover, there is evidence that the solid-solid transition terminates in a critical point as in the case of the fluid-fluid transition.
Colloidal systems observed in video microscopy are often analysed using the displacements correlation matrix of particle positions. In non-thermal systems, the inverse of this matrix can be interpreted as a pair-interaction potential between particles. If the system is thermally agitated, however, only an effective interaction is accessible from the correlation matrix. We show how this effective interaction differs from the non-thermal case by comparing with high-statistics numerical data from hard-sphere crystals.
We report numerical calculations of the concentration of interstitials in hard-sphere crystals. We find that, in a three-dimensional fcc hard-sphere crystal at the melting point, the concentration of interstitials is 2 * 10^-8. This is some three orders of magnitude lower than the concentration of vacancies. A simple, analytical estimate yields a value that is in fair agreement with the numerical results.