The Li + ion was allowed to optimize its position along the z-coordinate only. At the optimized position, z min , the single-point energy of the aromatic molecule and Li + complex, E Li + +mol (x, y, z min ), was calculated. The binding energy landscape in DFT, ∆E DFT (x, y, z min ), is defined as follows,