24小时热门版块排行榜    

CyRhmU.jpeg
查看: 943  |  回复: 2
当前主题已经存档。

vegalew

金虫 (小有名气)

[交流] 【求助】十万火急求救:CASTEP的优化时候SCF的smearing有什么讲究

在CASTEP模块里,Electronic options里面有四个选项,分别是basis,SCF,K-points,Potential,我在SCF设置的时候碰到困难了,在SCF选项的下面有个Orbital occupancy部分里面的smearing值是可以调的,默认一般是0.1,我把值变小了对非金属,氧化钛没有问题,但是单独计算Pt的时候边成0.01根本就不收敛.

这次我做新任务的时候又碰到问题了,我把值写成0.03,算72个原子的氧化钛的表面可以30步左右收敛,用了32个小时,但是在表面加了一个Pt原子和以前的72个氧化钛原子一起算,算了100个小时了,优化了50步还没停下来 我看的新里那叫一个急啊.一般优化算到多少步不收敛就可以停机了?这种贵金属和半导体一起算的时候smearing选择有什么标准么?

十万火急求救

[ Last edited by freshgirl on 2009-6-28 at 22:01 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★
zzgyb(金币+3,VIP+0):谢谢你的参与,欢迎继续关注计算模拟版!
个人认为,smearing参数是虚拟的"电子温度", 是人为引入的一种算法,使得费米能级之上的电子也可以被占据,在费米能级附近的电子态会出现分数占据.这样做的原因是加速收敛,尤其是费米能级电子结构比较复杂(如金属体系)的情况;但若准确DFT计算的话应该选择smearing=0,而如果直接这样选择对于一些体系会不收敛, 因此有文献提到几何优化时应该先选择大的smearing参数,优化之后再使用小的参数继续优化,直到逐渐降低为0
2楼2008-02-07 17:48:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
zzgyb(金币+3,VIP+0):谢谢你的参与,欢迎继续关注计算模拟版!
smearing可以用于scf或几何优化的中间过程,帮助收敛。但它本身是没有物理意义的,很多程序的manual都不建议用它得到最终结果。

对于不收敛的体系,可以改为开壳层(有些程序如crystal、adf-band还可以指定电子占据方式,但我不知道castep是否支持),或者用level shift。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2008-02-08 11:29:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vegalew 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见