24小时热门版块排行榜    

查看: 5357  |  回复: 24

xylz6188

铁杆木虫 (著名写手)


[交流] CASTEP Geometry Optimization

我用CASTEP优化结构时,出现能量曲线,先先下降,然后又升高,最后优化结束(100次迭代),最后也成功的结束,可是能量不是最低的,这到底是怎么回事???求高手赐教,大家一起学习探讨

优化曲线



收敛曲线



[ 来自科研家族 晶体材料家族 ]
回复此楼

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

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

awmc2008

至尊木虫 (文坛精英)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
uuv2010(金币+1): 欢迎参加交流 2011-11-24 18:27:34
xylz6188(金币+5): 2011-11-25 15:23:34
我估计曲线没有完全显示计算结果。你可以查看.CASTEP文件,结构优化应该不止45步。
9楼2011-11-24 18:15:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

乐石||

木虫 (正式写手)



小木虫(金币+0.5):给个红包,谢谢回帖
xylz6188(金币+5): 2011-11-25 09:17:08
是不是没显示全啊?点击你的结构,然后CASTEP-ANAYSIS-ENERGY EVOLUTION,优化时能量的变化就全显示出来了,嘿嘿
11楼2011-11-24 19:51:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jeeanlee

木虫 (正式写手)


★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
youzhizhe(金币+2): 谢谢交流。 2011-11-24 23:26:47
xylz6188(金币+5): 2011-11-25 09:15:33
引用回帖:
1楼: Originally posted by xylz6188 at 2011-11-24 11:47:58:
我用CASTEP优化结构时,出现能量曲线,先先下降,然后又升高,最后优化结束(100次迭代),最后也成功的结束,可是能量不是最低的,这到底是怎么回事???求高手赐教,大家一起学习探讨
[eimg]ae/97/761930_1322 ...

从你的收敛图上看,应该45步并非最后一步,因为此时你的所有收敛标准都没达到呢。因此45步时的能量(即你所认为的基态能量)并非是最终优化的真实能量值,所以,请你检查确认你的能量和收敛变化曲线已经显示完全了吗?
12楼2011-11-24 20:30:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sarana2010

新虫 (初入文坛)



小木虫(金币+0.5):给个红包,谢谢回帖
xylz6188(金币+5): 那请兄台赐教为谢!!! 2011-12-05 09:34:06
我可以负责任的告诉你,肯定有问题!
22楼2011-12-01 19:21:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)



小木虫(金币+0.5):给个红包,谢谢回帖
xylz6188(金币+5): 2011-12-05 09:33:12
看这个结果好像是没达标就优化完毕了,可以看看最后一步判据收敛是yes还是no来确认收敛的最后情况。
castep有时候并不是3个判据一起满足才完成,有时候2个判据满足就收敛了。具体的原因可以自己看看他的帮助文档关于收敛的说明。
23楼2011-12-03 14:56:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

aaq2800

铁杆木虫 (知名作家)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
xylz6188(金币+5): 2011-11-24 15:04:52
uuv2010(金币+1): 多谢提示 2011-11-24 18:26:21
提高截止能量再试试
2楼2011-11-24 14:28:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

52nana

铁虫 (小有名气)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
uuv2010(金币+1): 欢迎参加交流 2011-11-24 18:26:29
xylz6188(金币+5): 2011-12-01 09:47:52
你的计算最后显示succeed ?
看看有没有warning,估计是有,建议增大迭代次数看看。
3楼2011-11-24 14:53:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
2楼: Originally posted by aaq2800 at 2011-11-24 14:28:15:
提高截止能量再试试

我是通过收敛性测试之后找到平面波截断能的,这应该是没有问题啊
4楼2011-11-24 15:05:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
3楼: Originally posted by 52nana at 2011-11-24 14:53:04:
你的计算最后显示succeed ?
看看有没有warning,估计是有,建议增大迭代次数看看。

我的计算后面显示completed successfully,但是我曾经也增大迭代次数到300,只到100多点就算完了,优化曲线没有显示是能量最低的,而是和上图一样,是不是我的结构不好啊????
5楼2011-11-24 15:09:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
uuv2010(金币+1): 欢迎参加交流 2011-11-24 18:26:41
引用回帖:
4楼: Originally posted by xylz6188 at 2011-11-24 15:05:49:
我是通过收敛性测试之后找到平面波截断能的,这应该是没有问题啊

什么叫通过收敛性测试之后找到平面波截断能的?莫非能量设为100可以收敛,设置400就不能收敛了?
6楼2011-11-24 15:16:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
uuv2010(金币+1): 哪个软件可以用力场优化啊?谢谢指导! 2011-11-24 18:27:27
xylz6188(金币+5): 2011-11-25 15:34:07
可能是初始结构不太好,先clean一下,或者使用力场优化一下(很快的),然后再用castep优化计算
7楼2011-11-24 15:24:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡马西平

金虫 (正式写手)


heh ,不太了解
8楼2011-11-24 15:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘LOVEH

金虫 (小有名气)



小木虫(金币+0.5):给个红包,谢谢回帖
xylz6188(金币+5): 2011-11-25 15:26:02
我的也遇到过这种问题,不过优化到110多步,能量没有楼主的起伏大,后果是算的时候没有算出结果。
10楼2011-11-24 18:33:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
6楼: Originally posted by aaq2800 at 2011-11-24 15:16:33:
什么叫通过收敛性测试之后找到平面波截断能的?莫非能量设为100可以收敛,设置400就不能收敛了?

我通过计算得出,平面波截断能再300时是最优的,所以平面波阶段能应该没有问题,应该是我的结构不稳定、不够好导致呢!兄台想请教一下是不是每建一个周期性模型都必须切面|添加真空层啊,我是做石墨烯的,我没有切面|添加真空层,谢谢了
13楼2011-11-25 09:14:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
12楼: Originally posted by jeeanlee at 2011-11-24 20:30:24:
从你的收敛图上看,应该45步并非最后一步,因为此时你的所有收敛标准都没达到呢。因此45步时的能量(即你所认为的基态能量)并非是最终优化的真实能量值,所以,请你检查确认你的能量和收敛变化曲线已经显示完全 ...

从收敛曲线看,确实没有收敛,所以我感觉很困惑啊
14楼2011-11-25 09:16:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
11楼: Originally posted by 乐石|| at 2011-11-24 19:51:43:
是不是没显示全啊?点击你的结构,然后CASTEP-ANAYSIS-ENERGY EVOLUTION,优化时能量的变化就全显示出来了,嘿嘿

谢谢兄台,我试试
15楼2011-11-25 09:17:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

乐石||

木虫 (正式写手)


引用回帖:
15楼: Originally posted by xylz6188 at 2011-11-25 09:17:21:
谢谢兄台,我试试

嘿嘿~
16楼2011-11-25 10:20:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyh987

金虫 (小有名气)



小木虫(金币+0.5):给个红包,谢谢回帖
请问楼主如何进行收敛性测试找到平面截断能啊,我也一直遇到收敛的问题,能不能跟我说说具体步骤啊
17楼2011-11-25 10:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
9楼: Originally posted by awmc2008 at 2011-11-24 18:15:09:
我估计曲线没有完全显示计算结果。你可以查看.CASTEP文件,结构优化应该不止45步。

你说的对头,可是我的能量曲线实在是不好,能量忽高忽低的,有点反常啊,我这样也有优化结束呢,可是就是出现能量不是最低的构型
18楼2011-11-25 15:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
10楼: Originally posted by 刘LOVEH at 2011-11-24 18:33:08:
我的也遇到过这种问题,不过优化到110多步,能量没有楼主的起伏大,后果是算的时候没有算出结果。

请问兄台你出现这样的情况最后是怎么解决的呢
19楼2011-11-25 15:26:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
17楼: Originally posted by fyh987 at 2011-11-25 10:21:56:
请问楼主如何进行收敛性测试找到平面截断能啊,我也一直遇到收敛的问题,能不能跟我说说具体步骤啊

我是这样做的:先选取较大的截断能,比如800eV,然后K点从1*1*1逐渐增大到比如15*15*1,看能量怎么变,然后选取较大的k点,比如15*15*1,然后截断能从200eV逐渐增大到比如800eV,看你能量如何变化,从而选取你所需要的平面波截断能和k点。我是这样做的
20楼2011-11-25 15:31:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
7楼: Originally posted by 贺仪 at 2011-11-24 15:24:15:
可能是初始结构不太好,先clean一下,或者使用力场优化一下(很快的),然后再用castep优化计算

兄台我们想到一块了,我当时也是觉得我的初始结构不好,找了好多种方法来构建模型,几种建模方法都出现同样的问题,后来就真的不知道怎么解决了,现在看到你的点评,我来试试,谢谢了兄台
21楼2011-11-25 15:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xylz6188

铁杆木虫 (著名写手)


引用回帖:
23楼: Originally posted by redskywei at 2011-12-03 14:56:34:
看这个结果好像是没达标就优化完毕了,可以看看最后一步判据收敛是yes还是no来确认收敛的最后情况。
castep有时候并不是3个判据一起满足才完成,有时候2个判据满足就收敛了。具体的原因可以自己看看他的帮助文档 ...

好的,多谢指点!
24楼2011-12-05 09:33:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

415999915

铁虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
21楼: Originally posted by xylz6188 at 2011-11-25 15:37:45:
兄台我们想到一块了,我当时也是觉得我的初始结构不好,找了好多种方法来构建模型,几种建模方法都出现同样的问题,后来就真的不知道怎么解决了,现在看到你的点评,我来试试,谢谢了兄台

我也遇到了和楼主同样的问题 就是收敛曲线波动很大,而且出现上升的趋势,不知道楼主的问题是怎样解决的啊?
还有 有人说是先clean 一下,在用力场优化,这两个步骤是怎么操作的呢?
这个问题烦了很久啦,希望楼主不吝赐教,不胜感谢 ~~~~
25楼2012-05-16 10:58:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xylz6188 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
信息提示
请填处理意见