24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 2606  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

伊点若天

金虫 (小有名气)

[求助] Dmol3中参数的设置问题。

在对过渡金属氧化物进行结构优化时,成功优化了其中一个表面,但是对于另外一个表面,总是出现以下错误:
Message: SCF not converging. Choose "Use Smearing" on DMol3 SCF panel
or set "Occupation  Thermal" in the input file
You may also need to change spin or use symmetry
Resubmit DMol3
Message: DMol3 job failed
Error: DMol3 exiting
Use Smearing 选项也有勾选,Occupation Thermal值设置为0.005,use symmetry也有勾选,请教下Spin值的调整问题,怎样才能进行合理的参数设置来进行优化?或者是Spin值设定外的其他问题。
(第一个表面spin设定为0.5成功优化,第二个表面spin值也是0.5,出现上述问题,但不确定是否为spin值的原因)。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)

引用回帖:
3楼: Originally posted by 伊点若天 at 2013-03-24 21:32:11
有没有关于Fermi占据和分数占据的文章,我刚开始做,对这两个东西还不是很懂……...

一搜一大把啊。
一定要精通MaterialsStudios,VASP.
4楼2013-03-25 00:40:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

vasp001

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流,, 2013-03-25 22:48:16
伊点若天: 金币+10, ★★★很有帮助 2013-03-26 08:38:36
可以尝试使用
Occupation Thermal  0.005 ,
如果计算能够收敛的话,再把此值减小,直到,采用feimi占据,因为,要计算体系磁性,尽量使用Fermi占据,不然的话,会导致分数占据的。
一定要精通MaterialsStudios,VASP.
2楼2013-03-23 13:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伊点若天

金虫 (小有名气)

引用回帖:
2楼: Originally posted by vasp001 at 2013-03-23 13:43:13
可以尝试使用
Occupation Thermal  0.005 ,
如果计算能够收敛的话,再把此值减小,直到,采用feimi占据,因为,要计算体系磁性,尽量使用Fermi占据,不然的话,会导致分数占据的。

有没有关于Fermi占据和分数占据的文章,我刚开始做,对这两个东西还不是很懂……
3楼2013-03-24 21:32:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)

★ ★ ★ ★ ★
franch: 金币+5, 谢谢回帖交流,, 2013-03-25 22:48:44
引用回帖:
3楼: Originally posted by 伊点若天 at 2013-03-24 21:32:11
有没有关于Fermi占据和分数占据的文章,我刚开始做,对这两个东西还不是很懂……...

密度泛函理论发展到现在,已经能够相当好的解决很多实际的问题。特别是“B3LYP”由于其在有机化学中的巨大成功几乎已经成为DFT的代名词。DFT的核心就是“交换相关泛函”。如果这个东西的精确形式被确定,那么,也许从此以后量子化学家就会失业了。当然,短时间内这种可能性似乎不存在。这里通过一个有趣的现象:“分数占据数”来讨论一些交换相关泛函的一种系统误差。



我们知道,基于波函数的从头算方法是在3N维Hilbert空间中进行的,而DFT则是在3维的Euclid(真实)空间中进行的,在真实空间中,我们可以考察密度的变化对能量的影响。对于真实的N电子体系,我们能够确定其能量。那么,如果是“分数”电子呢?比如,H2+在两个H无穷远的时候,考虑对称性,每个核实际上是带1/2电子的,那么这时的能量如何计算呢?



利用量子统计的方法,Perdew等人证明,对于电子数N+delta (0


因此,对于精确的密度泛函理论,或者说,精确的交换相关泛函,能量E对电子数N的曲线应为下图1的红色线。很明显,在整数电子处,能量对电子数的一阶微分的不连续的(左,右导数不相等)!这是一个很奇异的性质。也就是说,下面这个关系成立:




并且,由于线性插值性:


其中I和A分别为电离能和电子亲和能。而且,我们发现,E对N是凸的,用化学语言就是说,第一电离能小于第二电离能,第二电离能小于第三电离能,等等。这个结论到目前为止无法从数学上证明,但是在实验和数值计算上没有发现例外。

那么近似的泛函满足这个条件吗?考虑一个模型,假设H2电子数为n(不一定是整数N),则相应的电子密度为:

用这个密度做X-alpha计算,可以得到n belongsto [0,2]的能量曲线,即图1中蓝色线:



图1



很明显,Xα泛函失败了,在N=1处,它的一阶微分是连续的,不存在斜率的跳跃。可以证明,对于包括B3LYP在内的所有现在的泛函,没有一个泛函能满足不连续微分条件!这意味着,对于H2+,有


这是不可能的。这个结果是不符合实际的,这是现代所有交换相关泛函的系统误差之一:它们缺乏在整数电子数处的不可微性。这就是分数占据数(fractional occupation)误差。这使得DFT容易高估离域分子的稳定性。TDDFT常常在计算大共轭体系时出现严重的偏差,原因可能也在于此。
http://prb.aps.org/abstract/PRB/v57/i3/p1427_1
一定要精通MaterialsStudios,VASP.
5楼2013-03-25 00:42:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见