The lyotropic phase behavior of cuboidal particles was investigated via Monte Carlo simulations. Hard cubes were approximated by suitably shaped clusters of hard spheres. Changes in concentration and structure of the system were monitored as a function of osmotic pressure P* (imposed in an isobaric ensemble). As expected, an isotropic phase prevailed at low concentrations (low P*) and a crystalline phase formed at high concentrations (high P*). A third distinct phase was also observed for an intermediate range of concentrations (approximately marked by breaks in the P* versus concentration curve). The structure of this mesophase was characterized both visually and analytically by calculating radial distribution functions and order parameters. It was found that such a mesophase exhibits orientational ordering along three axes (cubatic order) but significant translational disorder, thus having a structure clearly distinct from both isotropic and crystalline phases. (C) 2004 American Institute of Physics.