absoluteju 发表于 2019-3-23 16:02
Multiwfn版本是3.6版本，文件是这样产生的。先从materials studio里保存pbd文件，然后用Multiwfn转换成.x ...

跟文件大小没直接关系。不应该显示什么unknown

计算流程看起来没问题

3.6目前根本没有发布。对于3.6(dev)一定要看清楚你用的具体是哪天更新的，和官网上最后更新日期对照

absoluteju 发表于 2019-3-23 14:16
Sobereva老师，计算得到fch文件后，参考http://sobereva.com/443做静电势图，但是没有输出ES ...

如果用的不是最新Multiwfn版本，更新到最新再试
如果还有问题，把fch文件具体产生过程都告诉我。如果fch压缩后不超过5MB，上传我看看

sobereva 发表于 2019-3-21 11:32
和带电状态没直接关系。记得MS里只优化氢原子，不要优化重原子。

你目前的收敛精度还明显不够

感谢老师！

absoluteju 发表于 2019-3-20 19:45
结构是从剑桥数据库找的，在MS上优化过切割出来的，一共有4个分子，两个分子各带一个正电荷，另外两个各 ...

和带电状态没直接关系。记得MS里只优化氢原子，不要优化重原子。

你目前的收敛精度还明显不够
倒是可以考虑用xtb算一下得到.molden文件。xtb的速度极快，也容易SCF收敛，至少给你个定性正确的静电势分布应该是够的
然后把.molden用Multiwfn转换为.fch文件（为了能利用cubegen快速算静电势），再按照博文里的做法去绘图

sobereva 发表于 2019-3-20 19:36
结构怎么来的？如果初始结构可能有不合理性，可以先用分子力场优化下再用DFT搞

结构是从剑桥数据库找的，在MS上优化过切割出来的，一共有4个分子，两个分子各带一个正电荷，另外两个各带一个负电荷，不知道是不是这个原因，如果用Multiwfn+VMD画静电势图，不收敛的计算得到的gbw文件可以用吗

结构怎么来的？如果初始结构可能有不合理性，可以先用分子力场优化下再用DFT搞

plus 发表于 2019-3-19 20:50
纯GGA泛函加上KDIIS比较好用，有收敛的趋势可以加大迭代次数

我试了增加迭代次数，但是显示输入文件有问题，能帮忙写一下吗

纯GGA泛函加上KDIIS比较好用，有收敛的趋势可以加大迭代次数

理论方法改成B97-3c再试。虽然跟这个问题没直接关系，但这个级别比起BLYP-gCP-D3更好，而且耗时相仿佛，说不定恰好也没不收敛问题

如果不改理论方法，把%scf去掉，加上veryslowconv再试

liyuanhe211 发表于 2019-3-19 19:52
倒数第二步的能量变化上百Hartree，感觉（不太确定）是类似程序Bug发散了？

这个是输出结果，最后波动很大

倒数第二步的能量变化上百Hartree，感觉（不太确定）是类似程序Bug发散了？