Melting of Hard Cubes

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