When modeling extremely large cluster systems, we must accept a decrease in the applied level of theory. However, we still need to ensure that the applied methods yield reliable results. In the recent work by Engsvang et al., we benchmarked the cluster structures and binding energies for the (SA) n (AM) n clusters, with n = 1-6. Compared to the benchmark geometries, optimized at ωB97X-D/6-31++G(d,p), GFN1-xTB method performed well, providing similar cluster structures (RMSD = 0.31 Å) and thermal contributions to the free energy (MAE = 2.0 kcal/mol). It was also found that B97-3c yielded good agreement in binding energies with the benchmark DLPNO-CCSD(T 0 )/aug-cc-pVTZ level, with an MAE of 2.1 kcal/mol. Here, we extend this analysis to the SA-DMA and SA-AM-DMA systems, as well as assess clusters with one more acid or base molecule in the clusters. In addition, compared to our previous study, we also test how the empirically corrected PBEh-3c, B97-3c, and r 2 SCAN-3c functionals perform for obtaining the geometries.