there are many energetically-competitive crystal structures in these landscapes that are not distinguished successfully from the experimentally known forms by the force field energy model. Thus, the FIT+DMA results shown here are not representative of typical results with empirically parameterized, atomic multipole-based force fields, but represent situations where improvements are needed beyond such models.