Abstract
Organic molecules can be stable in distinct crystalline forms, known as polymorphs, which have significant consequences for industrial applications. Here, we predict the polymorphs of crystalline benzene computationally for an accurate anisotropic model parametrized to reproduce electronic structure calculations. We adapt the basin-hopping global optimization procedure to the case of crystalline unit cells, simultaneously optimizing the molecular coordinates and unit cell parameters to locate multiple low-energy structures from a variety of crystal space groups. We rapidly locate all the well-established experimental polymorphs of benzene, each of which corresponds to a single local energy minimum of the model. Our results show that basin-hopping can be both an efficient and effective tool for polymorphic crystal structure prediction, requiring no a priori experimental knowledge of cell parameters or symmetry.