Are high-<i>Z</i>′ polymorphs metastable? Insight from pharmaceutical polymorphs

High-<i>Z</i>' polymorphs are hypothesized to be kinetically trapped products and hence less stable than their low-<i>Z</i>' counterparts (<i>Z</i>' = number of symmetry-independent molecules in a crystal structure). Although the existence of such a symmetry-stability relationship would be of fundamental importance to the structural chemistry of molecular crystals, it has yet to be supported by experimental or computational evidence. A major challenge in analyzing high-<i>Z</i>' crystal structures is the large number of atoms, which makes accurate quantum chemical evaluation of their lattice cohesive energies (LCEs) computationally intensive. A systematic test of this hypothesis in drug polymorphs is made feasible by the CE-B3LYP method, which uses pairwise summation of interaction energies for LCE estimation. Here, we have analyzed 15 drugs (49 polymorphs) of low- and high-<i>Z</i>' crystal forms in terms of lattice cohesive energies, Kitaigorodskii packing indices (KPI), in-crystal molecular volumes, and crystal densities. Our results show no direct relation between <i>Z</i>' and the stability of drug polymorphs.