Multiwfn (http://sobereva.com/multiwfn) is able to generate Hirshfeld surface, which can then be plotted by VMD. Please check examples in Section 4.12.5 and 4.12.6 of Multiwfn manual. There is also a video tutorial showing the whole plotting process: https://youtu.be/Vy9pLQt4aDs