  Based on simplifications of previous numerical calculations [Graf and
L\"{o}wen, Phys. Rev. E \textbf{59}, 1932 (1999)], we propose algebraic free
energy expressions for the smectic-A liquid crystal phase and the crystal
phases of hard spherocylinders. Quantitative agreement with simulations is
found for the resulting equations of state. The free energy expressions can be
used to straightforwardly compute the full phase behavior for all aspect ratios
and to provide a suitable benchmark for exploring how attractive interrod
interactions mediate the phase stability through perturbation approaches such
as free-volume or van der Waals theory.
