“…47 We employed the triple-z Slater basis set, plus two polarization functions (STO-TZ2P) for valence electrons, within the generalized gradient approximation (GGA) according to the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional 48,49 because of its reliable performance at reasonable computational cost for similar clusters, leading to accurate geometrical parameters in good agreement to experimental results. [50][51][52][53] London dispersion correction to DFT via Grimme's empirical dispersion correction with the Becke-Johnson damping functions, D3(BJ), was employed for all the calculations. 54,55 A similar level of theory has been employed in related medium-sized gold nanoclusters.…”