STM calculation with VASP

First you need to have your geometry and DOS calculations done. Once that is done you need to generate a PARCHG file. This is done by adding this into the INCAR:

partial charge densities:
  LPARD = .TRUE.
  LSEPK = .FALSE.
  LSEPB = .FALSE.
  NBMOD = -3
  EINT = -1.50 0

Once this is done, you can plot the STM image with p4v:

Activate the "View structure" to be sure that you are above the surface and move the "Tip position" cursor to change the height of the tip.