Self-consistent phonon theory of the crystallization and elasticity of attractive hard spheres

We propose an Einstein-solid, self-consistent phonon theory for the crystal phase of hard spheres that interact via short-range attractions. The approach is first tested against the known behavior of hard spheres, and then applied to homogeneous particles that interact via short-range square well attractions and the Baxter adhesive hard sphere model. Given the crystal symmetry, packing fraction, and strength and range of attractive interactions, an effective harmonic potential experienced by a particle confined to its Wigner-Seitz cell and corresponding mean square vibrational amplitude are self-consistently calculated. The crystal free energy is then computed and, using separate information about the fluid phase free energy, phase diagrams constructed, including a first-order solid-solid phase transition and its associated critical point. The simple theory qualitatively captures all the many distinctive features of the phase diagram (critical and triple point, crystal-fluid re-entrancy, low-density coexistence curve) as a function of attraction range, and overall is in good semi-quantitative agreement with simulation. Knowledge of the particle localization length allows the crystal shear modulus to be estimated based on elementary ideas. Excellent predictions are obtained for the hard sphere crystal. Expanded and condensed face-centered cubic crystals are found to have qualitatively different elastic responses to varying attraction strength or temperature. As temperature increases, the expanded entropic solid stiffens, while the energy-controlled, fully-bonded dense solid softens.