24小时热门版块排行榜    

查看: 2777  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

匿名

用户注销 (小有名气)

本帖仅楼主可见

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

第一性原理学习

» 猜你喜欢

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

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

匿名

用户注销 (小有名气)

本帖仅楼主可见
3楼2014-06-28 11:07:31
已阅   申请QC强帖   回复此楼   编辑   查看我的主页
查看全部 9 个回答

卡开发发

专家顾问 (著名写手)

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的回帖

卡开发发

专家顾问 (著名写手)

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强帖   回复此楼   编辑   查看我的主页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +12 我4大白菜 2026-08-20 12/600 2026-08-20 20:42 by 6543yes
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见