24小时热门版块排行榜    

查看: 2457  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 求机械类四区sci推荐 5+4 迷茫小旷 2024-06-14 6/300 2024-06-20 17:12 by 中外文期刊编辑
[基金申请] 国产期刊影响因子大于12的有20多个个了 +3 babu2015 2024-06-20 3/150 2024-06-20 16:34 by 3115321
[找工作] 高校两个offer选择 +15 cowox2021 2024-06-18 16/800 2024-06-20 15:24 by ase123456
[文学芳草园] 累并快乐着 +21 MYHLD521 2024-06-14 21/1050 2024-06-20 13:34 by 1429524838
[硕博家园] 豫北虫友互识 +11 xuhongli903 2024-06-18 12/600 2024-06-20 13:32 by 1429524838
[基金申请] F03青年基金函评结果 +3 暨阳一只柴 2024-06-19 4/200 2024-06-20 11:15 by 暨阳一只柴
[基金申请] 太卷了 +14 laoyuefubio 2024-06-17 27/1350 2024-06-20 09:52 by htjwqy
[论文投稿] 投稿求助 6+3 小短腿安琪拉 2024-06-19 3/150 2024-06-20 08:59 by nono2009
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中! +10 恐龙爸爸 2024-06-14 10/500 2024-06-20 08:42 by jsguo4127
[论文投稿] 有化学投稿的期刊推荐吗 20+3 荣小撇 2024-06-13 4/200 2024-06-20 08:24 by 1975632963
[海外博后] 墨尔本大学博后offer要不要接 +3 kyxblmm 2024-06-18 3/150 2024-06-19 22:39 by blake1111
[基金申请] F口信息学部拿面上,大概需要什么样的成果 +7 _奋黎_ 2024-06-16 18/900 2024-06-19 20:02 by _奋黎_
[论文投稿] 审稿 +5 香瓜木香 2024-06-19 6/300 2024-06-19 17:44 by xli1984
[基金申请] 基金得中 +4 woaini0218 2024-06-18 4/200 2024-06-19 17:27 by mengzl
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +8 babu2015 2024-06-14 8/400 2024-06-19 15:25 by 风今25
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀? +7 65syn 2024-06-14 11/550 2024-06-19 10:19 by 65syn
[公派出国] CSC德国博后每个月资助多少呀?够用吗 +4 326lhpqk 2024-06-16 7/350 2024-06-19 02:03 by PLHOU
[考博] 2025考博 +8 自强不息a?a 2024-06-15 13/650 2024-06-18 18:12 by 投必得科研顾问
[论文投稿] 审稿问题:为什么荧光激发波长和紫外吸收波长差的大? 10+5 sdawege 2024-06-14 10/500 2024-06-17 18:54 by HH-探针
[论文投稿] 二审返修送审10天了,原来一审的3个审稿人只有2个接受了审稿,会邀请新审稿人么? 50+3 huanpo116 2024-06-15 5/250 2024-06-16 10:27 by bobvan
信息提示
请填处理意见