To quantitatively compare our calculated conductances with those of GNP assemblies, one could incorporate our single-interface model in a multi-time-scale framework. As this would be computationally demanding, we apply a highly simplified geometric model, where the GNPs are assumed to be in a primitive cubic packing and to have no internal resistance (see Supporting Information, Section S2. A more obvious (and serious) shortcoming of ReaxFF for (singly deprotonated) dithiols on single Au(111) surfaces was that, in contrast to monothiols, it led to unlikely chemical reactions involving the free thiol tail groups of dithiol ligands, including S-C bond formation (to a lesser degree, this also happened in sandwich structures). As a pragmatic solution, we set up a Python script capable of detecting S-C-S, S-H, and C-S-C bond formations and excluded the affected structures from further interpretation. On the long run, further reparametrization of ReaxFF or machine-learned force fields would be a more satisfying solution.