“…Density functional theory (DFT) calculations of this study done with the Gaussian 09 program package [8]. A selection of functionals were used, namely the GGA (Generalized Gradient Approximation) functional PW91 [9] and LC-BLYP [10,11,12,13] (long range corrected BLYP), meta-GGA functional M06-D3 [14] (M06 with the Grimme empirical dispersion correction D3 [15]) and the hybrid functionals B3LYP (Becke 1993 and Lee-Yang-Parr) [16,17] and CAM-B3LYP [18] (long range corrected version of B3LYP), with the triple-ζ basis set 6-311G (d,p) [20,21,22]. The LC-BLYP optimised gas phase structures were used to further conduct a natural bond orbital (NBO) analysis (using the NBO 3.1 module [23] in Gaussian 09), as well as an electronic density analysis (using Bader's quantum theory of atoms in molecules (QTAIM) [24,25,26], as implemented in ADF2017 [27,28,29]).…”