  The melting transition of a system of hard cubes is studied numerically both
in the case of freely rotating cubes and when there is a fixed orientation of
the particles (parallel cubes). It is shown that freelly rotating cubes melt
through a first-order transition, whereas parallel cubes have a continuous
transition in which positional order is lost but bond-orientational order
remains finite. This is interpreted in terms of a defect-mediated theory of
melting
