图便宜的话，B97-3c做优化，用直接得到的（不用额外再算单点）molden文件作为下文脚本的输入文件绘制表面静电势图
使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图（含视频演示）
http://sobereva.com/443（http://bbs.keinsci.com/thread-11080-1-1.html）

http://bbs.keinsci.com/thread-13560-1-1.html

但我一般使用ORCA计算RESP电荷的懒人脚本进行优化+单点计算，并且注释这一行rm SP.* SP_* Nval.txt，以保留我想要的含有波函数信息的文件。