“…The system under study represents a difficult balancing act between dispersion interactions and metal coordination in which pairwise semi-classical approaches like the D-models can lead to misleading results as previously observed, for example, in the study of cryogenic ion vibrational predissociation (CIVP), 133 bond dissociation energies 134 as well as energies, 135 prediction of vibrational spectra 136 and interconversion barriers. 137 We finally note that as shown in this work the inclusion of dispersion corrections, at least via the D3(BJ), D4, and VV10 approaches, although decisive for the description of the interaction energy between the stacked ligandpairs, represents only a necessary but not sufficient condition for an unambiguous prediction of the molecular structures. Therefore, including dispersion corrections of this kind is not always the preferred option when optimizing geometries with DFT.…”