24小时热门版块排行榜    

查看: 3701  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 22:19 by BqkHsOAGISWU
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 PnssoZglMj3v 2026-08-19 4/200 2026-08-19 22:16 by BqkHsOAGISWU
[基金申请] 今天基金会出结果吗?20260819 +12 kkkl_v 2026-08-19 12/600 2026-08-19 22:12 by 小布吉岛
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +54 医学老男孩 2026-08-13 121/6050 2026-08-19 20:28 by 医学老男孩
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 5/250 2026-08-19 17:42 by oCMhZAsgNFoO
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
信息提示
请填处理意见