Assessment of lattice energy minimization for the prediction of molecular organic crystal structures

Lattice energy searches for theoretical low-energy crystal forms are presented for 50 small organic molecules, and we compare the experimentally observed crystal forms to these lists of hypothetical polymorphs. For each known crystal, the relative stability is calculated with respect to the global minimum energy structure, and we determine the number of unobserved structures lower in energy than the experimental form. The distributions of these relative energies and their rankings in the predicted lists are used to determine the efficacy of lattice energy minimization in crystal structure prediction. Although a simple form for the interaction energies has been used, the calculations produce almost a third of the known crystals as the global minimum in energy, and approximately a half of the known structures are within 1 kJ/mol of the global minimum. Molecules with no hydrogen-bonding capacity are most likely to be found close to the global minimum in lattice energy, while increasing the number of possible hydrogen-bond donor-acceptor combinations leads to less reliable predictions.