For the parameterization of the potential, we fitted the binding energy landscape obtained from MS, ∆E MS (x, y, z min ) (see Eq. 4), with the reference binding energy obtained from DFT, ∆E DFT (x, y, z min ) (see Eq. 3), by tuning some of the interactions between carbon/sulfur from the aromatic molecule and Li + .