We have studied the structures and stabilities of Au6 to explore the origin of the large discrepancy between relative energies obtained from the density functional theory (DFT) and ab initio correlated levels of theory. The MP2 methods significantly overestimate the stability of the non-planar isomer when the double-$zeta$ polarization quality of basis sets, such as LANL2DZ+1f and CEP31G+1f, are used. However, we show that such preference for the non-planar structure at the MP2 level mainly originates from the large basis set superposition error.