请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 294|回复: 2

[NWChem] 分数占据数计算FON问题

[复制链接]

115

帖子

0

威望

1394

eV
积分
1509

Level 5 (御坂)

发表于 2018-2-7 17:03:04 | 显示全部楼层 |阅读模式
对分子进行了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是正常。



FON_problem.png

1万

帖子

25

威望

1万

eV
积分
30943

管理员

公社社长

发表于 2018-2-7 17:44:43 | 显示全部楼层
可以再看看比如-0.05的情况
也许是收敛问题,没收敛到对应情况的基态
北京科音自然科学研究中心:http://www.keinsci.com  不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训
思想家公社的门口Blog:http://sobereva.com  发布大量原创计算化学相关博文
Multiwfn量子化学波函数分析程序主页:http://sobereva.com/multiwfn
计算化学公社论坛:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

115

帖子

0

威望

1394

eV
积分
1509

Level 5 (御坂)

 楼主| 发表于 2018-2-9 10:58:09 | 显示全部楼层
应该是收敛的问题,感觉nwchem没Gaussian好用。
最后用关键词convergence rabuck 20 levlon 0.2解决的问题
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-5-24 20:12 , Processed in 0.201782 second(s), 26 queries .

快速回复 返回顶部 返回列表