  We describe a Monte Carlo procedure which allows sampling of the disjoint
configuration spaces associated with crystalline and fluid phases, within a
single simulation. The method utilises biased sampling techniques to enhance
the probabilities of gateway states (in each phase) which are such that a
global switch (to the other phase) can be implemented. Equilibrium
freezing-point parameters can be determined directly; statistical uncertainties
prescribed transparently; and finite-size effects quantified systematically.
The method is potentially quite general; we apply it to the freezing of hard
spheres.
