nwchem输出的.cube文件

lianghua飘

请教sob老师，想做一个电子密度差图（双电荷减去单电荷），用nwchem产生了电子密度的cube文件，按照您手册的步骤，在用multiwfn算电子密度差图时，屏幕没有电子密度这一选项，请问这是怎么回事呢，谢谢。

sobereva

cube文件里不含波函数信息，显然没法计算电子密度，获得密度差。
cube里仅有格点数据和原子坐标信息，顶多通过两个含有密度格点数据的cube文件用multiwfn主功能13里的格点运算功能求差。

lianghua飘

sobereva 发表于 2018-1-12 12:38
cube文件里不含波函数信息，显然没法计算电子密度，获得密度差。
cube里仅有格点数据和原子坐标信息，顶多 ...

谢谢sob老师，nwchem如何产生波函数，从而计算密度差呢，谢谢

sobereva

lianghua飘 发表于 2018-1-12 14:28
谢谢sob老师，nwchem如何产生波函数，从而计算密度差呢，谢谢

目前没有简单的办法。
目前最新版NWChem (6.8)虽然支持输出.molden文件，但是内容严重不标准，虽然能被Multiwfn读取，但结果是错的
利用JANPA程序自带的工具也能把nwchem的波函数信息转换成.molden文件形式，但用起来比较麻烦，而且也未必能被Multiwfn正常载入。
简单来说建议换个程序。