计算化学公社

标题: 分数占据数计算FON问题 [打印本页]

作者
Author: superrice 时间: 2018-2-7 17:03
标题: 分数占据数计算FON问题
对分子进行了FON的计算，计算都正常收敛。但是对缺电子体系的结果如下，其中p代表正电荷，p0.50代表(N-0.5)体系：

p0.00.out:       Total DFT energy = -2479.653546498976
p0.10.out:       Total DFT energy = -2479.518121053243
p0.20.out:       Total DFT energy = -2479.502377444411
p0.30.out:       Total DFT energy = -2479.486829518009
p0.40.out:       Total DFT energy = -2479.471498662635
p0.50.out:       Total DFT energy = -2479.456375104017
p0.60.out:       Total DFT energy = -2479.441466023836
p0.70.out:       Total DFT energy = -2479.426773457255
p0.80.out:       Total DFT energy = -2479.412296466813
p0.90.out:       Total DFT energy = -2479.398077776328
p1.00.out:       Total DFT energy = -2479.384068279371

发现p0.00和p0.10之间的差值非常大，而p0.10到p1.00之间的差值很小。导致E(Delta_N)-E()对Delta_N的绘图时导致如图问题，请问这是nwchem的数值误差吗，该如何解决，谢谢！。
图中问题泛函是LC-wPBE_2，而LC-wPBE_1是正常。

作者
Author: sobereva 时间: 2018-2-7 17:44
可以再看看比如-0.05的情况
也许是收敛问题，没收敛到对应情况的基态

作者
Author: superrice 时间: 2018-2-9 10:58
应该是收敛的问题，感觉nwchem没Gaussian好用。
最后用关键词convergence rabuck 20 levlon 0.2解决的问题

欢迎光临计算化学公社 (http://bbs.keinsci.com/)