On the DFT ground state of crystalline bromine and iodine

We report on an erroneous ground state within common density functional theory (DFT) methods for the solid elements bromine and iodine. Phonon computations at the GGA level for both molecular crystals yield imaginary vibrational modes, erroneously indicating dynamic instability-that fact alone could easily pass as a computational artefact, but these imaginary modes lead to energetically more favorable and dynamically stable structures, made up of infinite monoatomic chains. In contrast, meta-GGA and hybrid functionals yield the correct energetic order for bromine, while for iodine, most global hybrids do not improve the GGA result significantly. The qualitatively correct answer, in both cases, is given by the long-range corrected hybrid LC-ωPBE, the Minnesota functionals M06L and M06, and by periodic Hartree-Fock and MP2 theory. This poor performance of economic DFT functionals should be kept in mind, for example, during global structure optimizations of systems with significant contributions from halogen bonds.