24小时热门版块排行榜    

查看: 3917  |  回复: 12

flymice

金虫 (著名写手)

[求助] 关于dmol3的收敛性测试

att,我正在用ms的dmol3计算吸附能,在查阅K点设置相关帖子时,看到很多关于收敛性测试的帖子。castep模块的electronic里有energy cutoff选项,但是在dmol3里面却没有,想问一下,是否dmol3 是不需要进行收敛性测试的?那K点应该怎么取比较合适?
回复此楼

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

ms计算模型及其方法

» 猜你喜欢

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

Believeinyourself,impossibleisnothing.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

uuv2010

荣誉版主 (职业作家)

优秀版主

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
WDD880227(金币+1): 感谢交流提示 2012-03-12 16:14:44
flymice(金币+2): 有帮助, 恩,默认fine的话,有时候并不收敛的,而且有时候计算时间超长,所以必须得选择合适的K点。 2012-03-12 18:00:41
原则上k点和截断能都需要测试,一般的计算中用fine选项一般会满足要求,但是为了有把握,还是测试一下的好
10楼2012-03-12 16:12:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

navyzhang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
WDD880227(金币+1): 感谢回帖交流 2012-03-12 09:37:41
flymice(金币+5): ★★★很有帮助, 谢谢帮助。 2012-03-12 18:00:06
引用回帖:
5楼: Originally posted by flymice at 2012-03-09 17:04:55:
也就是说,必须取多个K点,cutoff不变,每个K点计算一遍?
然后取一个相对合理的K点,变化cutoff,每个cutoff计算一遍?
是否可以这样理解收敛性测试?

此外,K点是a*b*c的形式,怎么增加K点呢?如何在已知 ...

如果是bulk,可以abc三个同时加,不过一般情况下不用加到太大,
你可以参考一些文献,看看你的这个东西别人有没有算过,如果有算过可以不用测试,引用一下就行。
如果真是找不到文献,那就从4*4*4开始吧,这个很多时候凭经验的。
另外表面吸附,你是不是用slab模拟表面,那么abc中应该有一个是1才对,比如slab表面法向方向沿着z轴方向,那么c=1,其他两个测试一下。不过通常表面体系比较大,如果计算资源不是超强的话,建议不要超过4*4*1。
9楼2012-03-12 09:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

navyzhang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
fzx2008(金币+1): 谢谢回帖交流! 2012-03-09 12:22:38
flymice(金币+2): 有帮助 谢谢,不知道测试的具体操作是怎么样的,不知能否截图表示? 2012-03-09 13:02:21
理论基础不一样,electronic里有orbital cutoff选项,可供测试
k点可以测试啊!
2楼2012-03-09 11:55:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flymice

金虫 (著名写手)

引用回帖:
: Originally posted by navyzhang at 2012-03-09 11:55:04:
理论基础不一样,electronic里有orbital cutoff选项,可供测试
k点可以测试啊!

在orbital cutoff那里倒是可以设置orbital cutoff的值,K点也可以在k-points里面设 置。我现在疑惑这个测试究竟是个怎么操作法?依据帖子中的意思,似乎不是把就是个想要的K点,每一组都算一遍?
    我在help文档里面没有找到convergence test 的具体类容,求兄台不吝赐教。
Believeinyourself,impossibleisnothing.
3楼2012-03-09 13:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

navyzhang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
sunyang1988(金币+1): 谢谢交流 2012-03-09 15:22:45
flymice(金币+5): 有帮助 谢谢,但还是有问题不太明白,望讲解。 2012-03-09 17:05:43
k点和cutoff由大到小,选取一定的间隔,找最小的体系,算一遍
然后看变化关系到那个点收敛,就取该点的k和cutoff
(k点和cutoff大到一定程度就对计算结果没有太大影响)
4楼2012-03-09 14:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flymice

金虫 (著名写手)

引用回帖:
4楼: Originally posted by navyzhang at 2012-03-09 14:40:42:
k点和cutoff由大到小,选取一定的间隔,找最小的体系,算一遍
然后看变化关系到那个点收敛,就取该点的k和cutoff
(k点和cutoff大到一定程度就对计算结果没有太大影响)

也就是说,必须取多个K点,cutoff不变,每个K点计算一遍?
然后取一个相对合理的K点,变化cutoff,每个cutoff计算一遍?
是否可以这样理解收敛性测试?

此外,K点是a*b*c的形式,怎么增加K点呢?如何在已知晶格结构的情况下,选择相对合理的点?K点选取与晶格结构有关,高手当然可以自己根据晶格结构推导,但是我们对固体物理不熟悉(勿笑,正在学习中)的人,可以通过这个网站
http://www.cryst.ehu.es/cgi-bin/ ... ph-kv-list?gnum=129
例如这里面的label表下面的值是否都得取,随后连接成路径?

再麻烦你一下,就是在做分子在晶体表面吸附时,是否还需调整K点,还是沿用原来优化晶格时所设定的K点?
Believeinyourself,impossibleisnothing.
5楼2012-03-09 17:04:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天棚元帅

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang(金币+2): 多谢指教 2012-03-10 08:41:02
flymice(金币+3): 有帮助 谢谢。 2012-03-10 10:33:34
1、收敛性测试不仅仅只有K点,对于不同的计算体系还有其它不同需要测试项;

2、在不熟悉的情况下,建议用默认K点设置,再在此基础上简单测试下就行了;

3、多多参考相关的计算文献,看看别人对于类似于你的体系都做了哪些测试,具体又是怎么测试的~~
A man gets what he earns when he earns it.
6楼2012-03-10 08:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flymice

金虫 (著名写手)

引用回帖:
: Originally posted by 天棚元帅 at 2012-03-10 08:38:14:
1、收敛性测试不仅仅只有K点,对于不同的计算体系还有其它不同需要测试项;

2、在不熟悉的情况下,建议用默认K点设置,再在此基础上简单测试下就行了;

3、多多参考相关的计算文献,看看别人对于类似于你的 ...

恩,我现在就是使用默认的K点设置。
Believeinyourself,impossibleisnothing.
7楼2012-03-10 10:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

toushixiao

金虫 (正式写手)

呵呵,在这里半个板凳,学习了!
8楼2012-03-11 13:40:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 flymice 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 86/4300 2026-08-18 00:18 by longsantaizi
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见