Doml3中的smearing设置问题 - 量子化学 - 其他 - 小木虫论坛-学术科研互动平台

版块导航: 正在加载中...

客户端APP下载

论文辅导

申博辅导

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (3960)
>文献求助 (408)
>虫友互识 (345)
>导师招生 (323)
>博后之家 (166)
>考博 (155)
>硕博家园 (154)
>休闲灌水 (141)
>论文投稿 (118)
>基金申请 (90)
>招聘信息布告栏 (55)
>考研 (54)
>教师之家 (52)
>公派出国 (43)
>绿色求助(高悬赏) (38)
>找工作 (37)

返回列表

匿名

用户注销 (小有名气)

应助: 3 (幼儿园)
金币: 1282.6
散金: 853
红花: 1
帖子: 172
在线: 85小时
虫号: 0
注册: 2012-03-25
专业: 凝聚态物性 II ：电子结构

本帖仅楼主可见

» 收录本帖的淘帖专辑推荐

第一性原理学习

» 猜你喜欢

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

1楼2014-06-27 21:51:31

已阅同方向广播申请QC强帖回复此楼编辑查看我的主页

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

专家经验: +224
应助: 47 (小学生)
金币: 3480.8
散金: 14
红花: 223
帖子: 2477
在线: 1253.7小时
虫号: 1369423
注册: 2011-08-16
性别: GG
专业: 金属材料的磨损与磨蚀
管辖: 第一性原理

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gmy1990: 金币+3, thanks 2014-06-28 06:13:47
glddoyss(zhou2009代发): 金币+8 2014-06-29 20:10:04

说说我自己的观点，可能不一定完全正确，望各位前辈指正：
1、原则上热拖尾（或称为展宽、分数占据）对于孤立分子不应该考虑，热拖尾计算得到的能量和一些其他性质都是不准确的。
2、smearing来自于有限温度的方法，其中一种如Fermi-Dirac展宽，smearing的值=Boltzmann常量*温度。当0K时，电子在分子自旋轨道的占据只能是0或1，若前线轨道较为密集，则迭代过程占据、非占据轨道发生混合，能量振荡。若计算中考虑Fermi-Dirac展宽，电子就会按照这个F-D分布函数将一部分电子排到混合的原本的非占据态上，此时能量振荡减小。问题在于scf中得到的能量是电子的自由能，自由能和能量的关系F=E-TS，由于采用了“有限温度”，熵也不是0，导致能量计算不准确。
3、不知道你指的步长是？对于几何优化限制位移的步长可通过界面下的MaxstepSize进行调整。

Dmol3遇到难收敛的情况看场合，如果体系有磁性，那么需要合理设定每个原子的初始自旋，保证结构、体系电荷大致合理，这是第一步。然后这些都没问题先尝试减小charge mix，如果convergence整体呈下降趋势但较慢可适当增加DIIS size，除了这些Dmol3也可以像Gaussian那样才有能级提升的方法，将计算文件保存后将occupied这一项（原来可能是thermal或Fermi）修改为level_shift xxx（xxx为能级提升的值，单位为Ha），手头上没MS，有些地方记得不是很清楚，建议看看help。如果上述还不行，可适当采用smearing做优化，得到结构后再降低直至去除smearing进行优化。如果实在去除smearing无法进行计算，smeaing很小的时候，此时温熵项可能也很小，比如1meV/atom几乎可忽略，此时应读取的能量是外推到0K的键能，而非一般时候采用的总能。