24小时热门版块排行榜    

查看: 3576  |  回复: 17

xiaobai216

新虫 (初入文坛)

[交流] 【求助】收敛测试怎么做??? 已有7人参与

请教大家收敛测试具体该怎么做?可以说下步骤吗,比如确定Ecut的,谢谢
回复此楼

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

第一原理资料汇编 模拟知识汇总

» 猜你喜欢

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

第一性原理新手...非常非常新
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylli

木虫 (小有名气)

作两个图

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢指点 2010-09-27 22:42:17
逐步增加Ecut,观测总能随着Ecut增加的变化,直到总能基本不变,此时的Ecut为合理的Ecut.

此外,还要测试k点选取,方法类似. 即:绘制随着k点增加,总能的改变.

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

2楼2010-09-27 21:48:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaobai216

新虫 (初入文坛)

引用回帖:
Originally posted by ylli at 2010-09-27 21:48:05:
逐步增加Ecut,观测总能随着Ecut增加的变化,直到总能基本不变,此时的Ecut为合理的Ecut.

此外,还要测试k点选取,方法类似. 即:绘制随着k点增加,总能的改变.

那就是要建立一个固定模型,固定住其他参数改变Ecut来进行结构优化得到体系总能?我这么做过,但得到的结果很费解,随着E的增大,Total Energy 出现了一个峰值,这是怎么回事呢
第一性原理新手...非常非常新
3楼2010-09-27 22:20:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaobai216

新虫 (初入文坛)

引用回帖:
Originally posted by ylli at 2010-09-27 21:48:05:
逐步增加Ecut,观测总能随着Ecut增加的变化,直到总能基本不变,此时的Ecut为合理的Ecut.

此外,还要测试k点选取,方法类似. 即:绘制随着k点增加,总能的改变.

再麻烦问您下,收敛测试是不是在结构优化的过程中逐步改变E?
第一性原理新手...非常非常新
4楼2010-09-27 22:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylli

木虫 (小有名气)

在结构优化中改变Ecut

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
cenwanglai(金币+1):谢谢回复! 2010-09-28 08:16:45
有峰值有可能,但是应该还会下来,也就是最后的能量是最低的,不然可能你的结构不是局域最小(第一性原理中的结构优化不是全局优化).
5楼2010-09-28 07:44:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-11-09 13:50:20
全局优化=分子动力学优化+第一性原理优化?
6楼2010-09-28 08:17:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limin2007

铁杆木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by ylli at 2010-09-27 21:48:05:
逐步增加Ecut,观测总能随着Ecut增加的变化,直到总能基本不变,此时的Ecut为合理的Ecut.

此外,还要测试k点选取,方法类似. 即:绘制随着k点增加,总能的改变.

我增加截断能,总能量一直还在增加,还没有出现最小的数值,是不是还需要一直怎加截断能计算?你说的绘制图形,是不是每一次的截断能和对应的计算能量表示出来,在orgi里面绘制?
7楼2010-11-08 12:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ym23(金币+2):感谢指导 2010-11-08 16:33:57
引用回帖:
Originally posted by limin2007 at 2010-11-08 12:03:38:

我增加截断能,总能量一直还在增加,还没有出现最小的数值,是不是还需要一直怎加截断能计算?你说的绘制图形,是不是每一次的截断能和对应的计算能量表示出来,在orgi里面绘制?

不用试最小值,而是基本不变,比如小数点后面四五位在变化等等,当然由你自己的判别标准就可以
增加Ecut不会有最小值的,理论上肯定是Ecut越大能量应该最小,计算就需要考虑计算时间和计算精度的平衡了
直观上肯定是图形方便判断,E-Ecut的能量曲线,当Ecut增加到一定时,E就基本不再变化了就可以了
另外至于用什么画图,我的感觉是测试的话基本都差不多,如果文章中需要的话可以考虑用origin
博学、审问、慎思、明辨、笃学
8楼2010-11-08 14:45:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limin2007

铁杆木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by zzy870720z at 2010-11-08 14:45:01:


不用试最小值,而是基本不变,比如小数点后面四五位在变化等等,当然由你自己的判别标准就可以
增加Ecut不会有最小值的,理论上肯定是Ecut越大能量应该最小,计算就需要考虑计算时间和计算精度的平衡了
直观 ...

谢谢你,我选择一个截断能计算后,能量收敛图形也趋于一条直线了,但是我增加截断能后在计算,castep文件里面的总能量还是减小, 按照您的意思增加截断能总能量是不变才最好 是这样吗?
那我的问题是不是还需改变截断能继续计算?

[ Last edited by limin2007 on 2010-11-8 at 15:44 ]
9楼2010-11-08 15:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-11-09 00:14:14
引用回帖:
Originally posted by limin2007 at 2010-11-08 15:28:16:

谢谢你,我选择一个截断能计算后,能量收敛图形也趋于一条直线了,但是我增加截断能后在计算,castep文件里面的总能量增大了, 按照您的意思增加截断能总能量是不变才最好 是这样吗?
那我的问题是不是还需改变 ...

看最后的能量,最初来Etot vs Ecut 曲线在判断比较好
总能量应该是随着Ecut的增加减小才对

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

博学、审问、慎思、明辨、笃学
10楼2010-11-08 15:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaobai216 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:33 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ph7bxbho8f 2026-07-09 3/150 2026-07-10 05:53 by j4frgf6zcw
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[基金申请] 祈祷青基必中 50+3 hadengry 2026-07-09 9/450 2026-07-09 17:37 by 文杰猪猪
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[基金申请] 关于如何从代码看上不上会 +11 且听虎啸 2026-07-09 14/700 2026-07-09 15:27 by libeny
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评?应用基础研究和基础研究,希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[基金申请] 今年E04面上 +15 布布和一二 2026-07-05 21/1050 2026-07-09 10:25 by 学员fOzRO9
[硕博家园] 大龄残疾硕士的一点执念 +14 3776737374 2026-07-06 22/1100 2026-07-09 06:07 by lovechizi
[基金申请] 生命口会评 +4 wuke100666 2026-07-07 4/200 2026-07-08 16:38 by 逍遥不够自在
[教师之家] 网络举报学术不端成为了男女之间情杀的手段了,真荒唐。。。 +3 zju2000 2026-07-04 3/150 2026-07-07 21:00 by hyx552200
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联,生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
[论文投稿] 刑法学论文投稿求助 +3 高源678 2026-07-04 5/250 2026-07-06 12:20 by 平天尺
[论文投稿] Journal of Environmental Chemical Engineering +3 雪花天涯 2026-07-05 3/150 2026-07-06 09:00 by bobvan
[基金申请] 最后一年,祈求好运 +3 pyw1992 2026-07-06 3/150 2026-07-06 08:53 by 苏_frp
信息提示
请填处理意见