The DFTB theory was combined with the isothermal Brownian-type molecular dynamics (MD) and metadynamics molecular dynamics (MMD) algorithms to perform simulation studies for Au clusters. Two representative DFTB parametrizations were investigated. In one parametrization, the DFTB-A, the Slater–Koster parameters in the DFTB energy function were determined focusing on the ionic repulsive energy part, Erep and the other, the DFTB-B, due attention was paid to the electronic band-structure energy part, Eband. Minimized structures of these two parametrizations were separately applied in MD and MMD simulations to generate unbiased and biased trajectories in collective variable (CV) space, respectively. Here, we found the MD simulations monitored at 300 K manifest fluxional characteristics in planar cluster Au9/DFTB-A, but give no discernible tracts of fluxionality for planar Au8/DFTB-A and Au8/DFTB-B, for nonplanar Au10/DFTB-A and, to some extent, for nonplanar Au9/DFTB-B; they are plausibly being hindered by higher-than kBT energy barriers. Very recent FIR-MPD spectroscopy measurements, however, were reported to have detected at 300 K both the planar and nonplanar neutral Aun clusters in the size range 5 ≤ n ≤ 13. The failure of MD simulations has prompted us to apply the MMD simulation and construct the free energy landscape (FEL) in CV space. Through scrutinizing the FELs of these clusters and their associated structures, we examine the relative importance of Erep/DFTB-A and Eband/DFTB-B in unraveling the covalent-like behavior of valence electrons in Aun. Most important of all, we shall evaluate the DFTB parametrization in MMD strategy through comparing extensively the simulation data recorded with the gas-phase experimental data.