This work introduces a systematic and comprehensive computational approach for exploring the configurational space of large atmospheric clusters far beyond the size routinely studied in the literature. We find that parallel ABCluster runs are required to achieve a good guess for the low-energy cluster structures. In addition, we applied the B97-3c method for optimizing the geometries and calculating the vibrational frequencies.