24小时热门版块排行榜    

查看: 2432  |  回复: 14

菲顶顶

铜虫 (正式写手)


[交流] 收敛测试的讨论

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

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

模拟知识汇总 刺客

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注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

相关版块跳转 我要订阅楼主 菲顶顶 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 h4CP7TrQR8Lg 2026-08-22 3/150 2026-08-23 16:30 by lXgU0RvNTDC9
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 15:18 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 14:30 by KM5EcsNQRBPn
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见