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

24小时热门版块排行榜    

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

菲顶顶

铜虫 (正式写手)

[交流] 收敛测试的讨论

我想请教大家，在用Ms寻找适合体系的截断能和K点的时候，到底是怎么做的呢？     
看了很多帖子 都没有说的很详细，在做截断能和k点测试的时候task到底是选择energy还是选择优化，
我以前都是选择energy，然后固定其他的不变 改变截断能，然后把每次计算得出的，castep文件里的final energy记下来 看差不多两个能量差为0.01eV时就选这个截断能（但是看到有的人又说小于0.001eV）k点也是这样
但最近又看到有人说做这种测试的时候，task要选择优化，然后把。castep文件里的 Final Enthalpy记下来，看这个
我也不知道要怎么搞了，有些疑惑，大家可以来具体说下嘛？ 都进来交流一下吧
回复此楼

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

模拟知识汇总 刺客

» 本帖已获得的红花（最新10朵）

1357246

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-11-20 21:18:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jq120

木虫 (正式写手)
★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
菲顶顶: 金币+1 2012-11-21 09:05:39
sunyang1988: 金币+2, 谢谢指点 2012-11-22 21:22:23
菲顶顶: 金币+1 2012-12-09 19:11:46
我不知道你作什么体系，我作的有些体系结构是实验上已经报道的，这个时候我认为几何结构优化已经没有必要了，所需要的就是计算能量以及带隙等。那么我在作收敛测试的时候，当然不用再作结构哦偶优化。
一般而言，我认为对于一个陌生的体系，我们需要作的大概为：
1、结构优化，只有优化好的结构，后面的计算才有意义；
2、能量计算、带隙计算等；
在这些计算过程中，参数的选择当然很重要了，这包括cutoff的选择，K点设置等等，所以有必要作收敛测试。
赞一下 回复此楼
3楼2012-11-20 22:24:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

707462558

银虫 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
菲顶顶: 金币+1 2012-11-21 09:09:09
菲顶顶: 金币+1 2012-12-09 19:11:56
引用回帖:
3楼: Originally posted by jq120 at 2012-11-20 22:24:03
我不知道你作什么体系，我作的有些体系结构是实验上已经报道的，这个时候我认为几何结构优化已经没有必要了，所需要的就是计算能量以及带隙等。那么我在作收敛测试的时候，当然不用再作结构哦偶优化。
一般而言，我 ...

请问 “有些体系结构是实验上已经报道的，这个时候我认为几何结构优化已经没有必要了”不优化 你的体系就可以直接进行后续计算吗
赞一下 回复此楼
4楼2012-11-21 08:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
3楼: Originally posted by jq120 at 2012-11-20 22:24:03
我不知道你作什么体系，我作的有些体系结构是实验上已经报道的，这个时候我认为几何结构优化已经没有必要了，所需要的就是计算能量以及带隙等。那么我在作收敛测试的时候，当然不用再作结构哦偶优化。
一般而言，我 ...

那你的意思就是在做测试的时候只做单点能测试，那你的标准是多少呢？就是两次的能量差是多少呢？我看到有人说把这个总的能量差除以原子数，大概在0.001EV左右就可以了 是这样的吗？我以前就是直接比较总能差值
赞一下 回复此楼
5楼2012-11-21 09:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
4楼: Originally posted by 707462558 at 2012-11-21 08:58:16
请问 “有些体系结构是实验上已经报道的，这个时候我认为几何结构优化已经没有必要了”不优化 你的体系就可以直接进行后续计算吗...

我认为还是要优化下的
赞一下 回复此楼
6楼2012-11-21 09:09:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

707462558

银虫 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
菲顶顶: 金币+1 2012-11-21 16:26:15
菲顶顶: 金币+1 2012-12-09 19:12:01
引用回帖:
6楼: Originally posted by 菲顶顶 at 2012-11-21 09:09:31
我认为还是要优化下的...

不优化的画建模出来的体系可能是很粗糙  我的就是手画的  很粗糙的
赞一下 回复此楼
8楼2012-11-21 09:54:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
今天看了很多帖子  应该市选择单点能吧，然后就是优化我认为是必须要做的，测试的话精度差不多在0.001我感觉就可以了。是今天又看了侯老师的书，有点明白了了

还有就是顺便问大家一个问题 掺杂后优化结构总是不收敛 请问原因是什么呢？
赞一下 回复此楼
9楼2012-11-21 16:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
还有想问大家 我建立好超胞以后 k点是按比例增加 还是减少啊
赞一下 回复此楼
10楼2012-11-22 16:59:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jq120

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by 菲顶顶 at 2012-11-22 16:59:32
还有想问大家 我建立好超胞以后 k点是按比例增加 还是减少啊

看下固体物理或者其他的书。
实空间和倒空间存在一定关联。
建立超胞以后，K点设置相应可以减少的。超胞比较大，仅用Gamma点就可以了。
赞一下 回复此楼
11楼2012-12-12 12:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Brilliantly

捐助贵宾 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主你现在知道收敛的具体步骤了吗？求具体步骤，谢谢了
赞一下 回复此楼
12楼2013-10-02 12:29:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangyangwlhx

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
9楼: Originally posted by 菲顶顶 at 2012-11-21 16:26:06
今天看了很多帖子  应该市选择单点能吧，然后就是优化我认为是必须要做的，测试的话精度差不多在0.001我感觉就可以了。是今天又看了侯老师的书，有点明白了了

还有就是顺便问大家一个问题 掺杂后优化结构总是不 ...

楼主  您好 我在测E-CUT时 也是像您前面帖子上说的
我以前都是选择energy，然后固定其他的不变 改变截断能，然后把每次计算得出的，castep文件里的final energy记下来 看差不多两个能量差为0.01eV时就选这个截断能（但是看到有的人又说小于0.001eV）k点也是这样
请问这样到底对不对 CASTEP算的能量 相差能达到0.001？ 我计算下来最接近的也相差0.01  还有我在做收敛性测试时 在final energy 下面出现 energies not corrected for finite basis set  这是什么原因 谢谢您
我们的常规思维是不是就是先测e-cut和kpoint 然后结构优化 然后计算和分析   诚心求教  麻烦您了
赞一下 回复此楼
13楼2013-10-07 11:08:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
13楼: Originally posted by yangyangwlhx at 2013-10-07 11:08:13
楼主  您好 我在测E-CUT时 也是像您前面帖子上说的
我以前都是选择energy，然后固定其他的不变 改变截断能，然后把每次计算得出的，castep文件里的final energy记下来 看差不多两个能量差为0.01eV时就选这个截断 ...

是的 即说的是对的 关于那个energies not corrected for finite basis set你在小木虫上搜一下 有讲解的 我的也出现过 不过没什么大的影响 好像没事
赞一下 回复此楼
14楼2013-10-07 14:55:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangyangwlhx

木虫 (正式写手)
引用回帖:
14楼: Originally posted by 菲顶顶 at 2013-10-07 14:55:05
是的 即说的是对的 关于那个energies not corrected for finite basis set你在小木虫上搜一下 有讲解的 我的也出现过 不过没什么大的影响 好像没事...

好的 谢谢
回复此楼
15楼2013-10-07 20:27:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
13572462楼
2012-11-20 21:51   回复  
送鲜花一朵
菲顶顶: 金币+1 2012-11-21 09:53:06
菲顶顶: 金币+3 2012-12-09 19:12:12
菲顶顶7楼
2012-11-21 09:53   回复  
引用回帖:
2楼: Originally posted by 1357246 at 2012-11-20 21:51:40

相关版块跳转 我要订阅楼主 菲顶顶 的主题更新
151/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 04:21 by Equinoxhua
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:25 by HajghKOMynmz
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:14 by HajghKOMynmz
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:08 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 3/150 2026-11-12 08:26 by HajghKOMynmz
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:25 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:03 by 9fOTqj074tIg
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 06:03 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
信息提示
关闭
请填处理意见
关闭