Nevertheless, the results obtained for all our models are in good agreement with the literature. By performing AIMD and using the CSHE method, Bramley et al. 111 found values of 4.80 eV (rVV10) and 4.65 eV (PBE). When using the WF method, the authors found slightly higher values of 5.00 eV and 4.83 eV, further highlighting how the treatment of χ S can affect the workfunctions. Similarly, Le et al. 110 reported a value of 4.7 eV for the workfunction at the solvated Pt(111) interface using AIMD with the PBE+D3 functional.