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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 关于dmol3的收敛性测试
131/2返回列表12下一页
查看: 3128  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

flymice

金虫 (著名写手)

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

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

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

ms计算模型及其方法

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
Believeinyourself,impossibleisnothing.
1楼 2012-03-09 10:46:42
已阅   回复此楼   关注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选项一般会满足要求，但是为了有把握，还是测试一下的好
赞一下(2人) 回复此楼
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。
赞一下(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 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭