24小时热门版块排行榜    

查看: 2424  |  回复: 8

匿名

用户注销 (小有名气)

本帖仅楼主可见
已阅   同方向广播   申请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的键能,而非一般时候采用的总能。
不一定挂在论坛,计算问题问题欢迎留言。
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-16 10/500 2026-03-18 00:01 by 星空星月
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 材料与化工304求B区调剂 +7 邱gl 2026-03-11 8/400 2026-03-17 09:36 by 努力学习赚彩礼
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
信息提示
请填处理意见