Multi-impurity complexes for n-type diamond: a computational study