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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 其他 » Doml3中的smearing设置问题
91/1返回列表
查看: 2923  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

匿名

用户注销 (小有名气)
本帖仅楼主可见

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

第一性原理学习

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-27 21:51:31
已阅   同方向广播   申请QC强帖   回复此楼   编辑   查看我的主页

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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的键能，而非一般时候采用的总能。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
2楼2014-06-28 00:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
3楼2014-06-28 11:07:31
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
有这种可能性。不过你搜索完过渡态是否检查虚频和优化过渡态？另外你的体系是表面？还是孤立分子？这两种在smearing的处理可能都会有差异。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
4楼2014-06-28 12:19:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
5楼2014-06-28 14:26:20
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
5楼: Originally posted by glddoyss at 2014-06-28 14:26:20
表面，我的设置是看了一些文献之后弄的。...

金属表面还是非金属表面？
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
6楼2014-06-28 14:38:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
7楼2014-06-28 16:51:47
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

linabeiruoyi

新虫 (小有名气)
引用回帖:
4楼: Originally posted by 卡开发发 at 2014-06-28 12:19:57
有这种可能性。不过你搜索完过渡态是否检查虚频和优化过渡态？另外你的体系是表面？还是孤立分子？这两种在smearing的处理可能都会有差异。

请问反应能能垒是负值说明什么？这种反应非常容易发生？
赞一下 回复此楼
8楼2016-07-12 16:25:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
8楼: Originally posted by linabeiruoyi at 2016-07-12 16:25:51
请问反应能能垒是负值说明什么？这种反应非常容易发生？...

如果是自由能能垒的话这都不算是过渡态了，也有可能是搜索的结构不合理啊，比如不太物理的固定，或者smearing不合理，这些都有可能导致能垒负值。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
9楼2016-07-12 17:55:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 glddoyss 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 06:13 by xzcv0Tehd8lQ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 05:54 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 05:47 by xzcv0Tehd8lQ
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 05:41 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
关闭
请填处理意见
关闭