To understand how the electric field affects the adsorption geometries of isolated molecules, we start from the geometries obtained in earlier work by Egger et al. in absence of an external electric field. There, 11 distinct adsorption geometries were obtained by performing geometry optimizations from four different high-symmetry points (top, bridge, hcp-hollow and fcc-hollow) with the TCNE molecule in various orientations. Going beyond the method used by Egger et al., we replaced the (then state-of-the-art) vdW Surf dispersion correction with the more recent MBD-NL dispersion correction, which has been previously shown to yield adsorption energies in excellent agreement with experiment. Furthermore, we modelled the metal slab using 9 layers in order to obtain highly accurate adsorption energies. With these new settings, we re-optimized these geometries in a 6 × 6 super cell until the remaining forces fell below 0.01 eV Å -1