Because
of their anisotropic electron distribution and electron
deficiency, halonium ions are unusually strong halogen-bond donors
that form strong and directional three-center, four-electron halogen
bonds. These halogen bonds have received considerable attention owing
to their applicability in supramolecular and synthetic chemistry and
have been intensely studied using spectroscopic and crystallographic
techniques over the past decade. Their computational treatment faces
different challenges to those of conventional weak and neutral halogen
bonds. Literature studies have used a variety of wave functions and
DFT functionals for prediction of their geometries and NMR chemical
shifts, however, without any systematic evaluation of the accuracy
of these methods being available. In order to provide guidance for
future studies, we present the assessment of the accuracy of 12 common
DFT functionals along with the Hartree–Fock (HF) and the second-order
Møller–Plesset perturbation theory (MP2) methods, selected
from an initial set of 36 prescreened functionals, for the prediction
of
1
H,
13
C, and
15
N NMR chemical
shifts of [N–X–N]
+
halogen-bond complexes,
where X = F, Cl, Br, and I. Using a benchmark set of 14 complexes,
providing 170 high-quality experimental chemical shifts, we show that
the choice of the DFT functional is more important than that of the
basis set. The M06 functional in combination with the aug-cc-pVTZ
basis set is demonstrated to provide the overall most accurate NMR
chemical shifts, whereas LC-ωPBE, ωB97X-D, LC-TPSS, CAM-B3LYP,
and B3LYP to show acceptable performance. Our results are expected
to provide a guideline to facilitate future developments and applications
of the [N–X–N]
+
halogen bond.