minimal model of the graphite surface and its corrugation was realized with two carbon layers. The bottom layer was frozen during the simulations to avoid the graphite surface to bend due to the small thickness of the simulated surface. Conversely, the atoms of the top graphitic layer were allowed to change their nuclear coordinate during molecular dynamics simulations. In all simulations, conducted by employing 3D periodic boundary conditions, the vertical side of the box, corresponding to the direction normal to the graphite surface, was fixed to 10 nm to allow enough vacuum space between the periodic graphite surface replicas.