24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3002  |  回复: 13

xueht987

木虫 (正式写手)

[求助] CASTEP 能量曲线上升 但是收敛了

采用几何优化——勾选了more中的优化晶胞,结果能量曲线上升,但是收敛,这种计算结果可靠吗?很怀疑!
除了改变ECUT、K点,还有办法吗?这能否说明我建的晶胞就是比较稳定的,无需优化呢?
如果能用的话,我该怎么确定优化之后最终的(稳定)晶胞的能量以及晶胞参数,找最低的那个能量吗?
期待大家的指点,谢谢!


回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueht987

木虫 (正式写手)

自己先顶一下!这个是虚晶近似计算出现的现象,呵呵,欢迎大家来讨论指点!
2楼2011-07-21 22:36:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzl8181

铁杆木虫 (著名写手)

单身浏览汉

【答案】应助回帖

★ ★ ★
franch(金币+3): 谢谢指点 2011-07-22 09:06:20
xueht987(金币+5): 谢谢指点哦,确实是这样 2011-07-22 10:10:56
应该是初始输入的晶格常数过小造成的,这是内聚能偏低,你看一下优化后的晶格和优化前的比一下是不是这样。
版权乜有,盗用不究!
3楼2011-07-22 07:34:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueht987

木虫 (正式写手)

是这样的,优化前的模型我采用的是实验值(8GPa压力下),优化过程中没有加压(0压),优化后晶胞比优化前有所膨胀,晶胞参数变大。请教楼上,这能否理解为:由于能量升高,所以优化后的晶胞相对于优化前是一种次稳定状态,但在0压下就应该是这种存在状态。
4楼2011-07-22 11:17:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzl8181

铁杆木虫 (著名写手)

单身浏览汉

★ ★
zzy870720z(金币+2): 谢谢指教 2011-07-23 16:44:29
‘所以优化后的晶胞相对于优化前是一种次稳定状态’只能说优化前的晶格在8GP时可能是稳定的(这是还要看其声子谱和弹性常量才能确定是否真的稳定);‘0压下就应该是这种存在状态'’0压下应该是优化后的结构,但是否真正稳定也不好说;这两种情况都可能处于亚稳态(这时声子谱或弹性常量不满足稳定条件的状态),当然也可能都处于稳定态。。。
版权乜有,盗用不究!
5楼2011-07-23 06:50:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

引用回帖:
Originally posted by lzl8181 at 2011-07-22 07:34:05:
应该是初始输入的晶格常数过小造成的,这是内聚能偏低,你看一下优化后的晶格和优化前的比一下是不是这样。

能再具体解释一下吗?ms的那张图给的是焓值,而不是能量值啊!结构优化本身就是在寻在焓的极小值,沿着力的方向移动原子和晶格,能量一定是减小的。从正常算法上来看是如果应力的计算是正确的,那么不可能出现能量升高的情况的。
6楼2011-07-25 16:02:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

★ ★
franch(金币+2): 谢谢回帖交流 2011-07-25 20:49:41
楼主的情况感觉是选了对称性优化,或者是加了一些束缚。才可能出现这样的问题。
7楼2011-07-25 16:14:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzl8181

铁杆木虫 (著名写手)

单身浏览汉

★ ★ ★
franch(金币+3): 谢谢回帖交流 2011-07-25 20:50:20
如果在oPa 下它是稳相,而在加压直8GPa 过程中没有相变的话,焓确实是应该升高的。你的体系是否有相变?你可以从零压下逐渐加压再试试。MS给出的是焓还是能量我还真没想过,我一直觉得它是能量,稍后核实一下。H=E+PV,零压下的焓就是能量。
版权乜有,盗用不究!
8楼2011-07-25 20:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueht987

木虫 (正式写手)

引用回帖:
Originally posted by lzl8181 at 2011-07-25 20:12:37:
如果在oPa 下它是稳相,而在加压直8GPa 过程中没有相变的话,焓确实是应该升高的。你的体系是否有相变?你可以从零压下逐渐加压再试试。MS给出的是焓还是能量我还真没想过,我一直觉得它是能量,稍后核实一下。H= ...

谢谢大家的关注!按楼上的试了一下,加压过程中,能量走势和贴出来的类似,往上走,但是8GPa时却往下了。查看总能,发现0-8GPa下,E-P曲线起起伏伏,波浪形的,8GPa这个点能量陡降,和前面的差值比较大。不知何解,请指教
9楼2011-07-26 17:34:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzl8181

铁杆木虫 (著名写手)

单身浏览汉

楼主,不用管了,“结果能量曲线上升,但是收敛”,说明在你的设定收敛参数下它收敛了,结果应该是可靠的,但有可能是亚稳态。“这能否说明我建的晶胞就是比较稳定的,无需优化呢?”根据你的收敛曲线,刚开始时stress并没收敛,说明stress过大,这有可能是原来加压8Gpa引起的,所以MS要调整晶胞参数,使它收敛,这时能量就要发生变化了。
从E-P曲线看不出什么的话(注意加压时最好在前一个加压构型的基础上加,这样不会产生大的突变误差),你再画一下H-P曲线(H=E+PV),看看曲线如何变化,这是要是波浪式变化的话就真有问题了!
版权乜有,盗用不究!
10楼2011-07-26 19:41:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xueht987 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323求调剂(计算机视觉和大模型项目经历) +3 chaoxiicy 2026-03-31 3/150 2026-04-05 10:33 by zhq0425
[考研] 求调剂 +3 电气小神童 2026-04-04 3/150 2026-04-05 10:17 by barlinike
[考研] 求调剂 一志愿西南交通大学085701环境工程 282分 +5 多多爱吃汉堡 2026-04-04 5/250 2026-04-04 22:14 by 啵啵啵0119
[考研] 085601,一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-03-31 5/250 2026-04-04 19:45 by 1753564080
[考研] 求调剂 +6 朔朔话 2026-04-02 7/350 2026-04-04 19:16 by 蓝云思雨
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +4 。。。54 2026-04-03 4/200 2026-04-04 16:09 by dongzh2009
[考研] 一志愿哈尔滨工业大学085600英一数二337分求调剂 +11 lyz0427 2026-04-03 11/550 2026-04-04 15:31 by dongzh2009
[考研] 一志愿沪985,326分求调剂 +3 刘墨墨 2026-04-03 3/150 2026-04-04 11:16 by 悲伤的芋头
[考研] 303求调剂 +9 DLkz1314. 2026-03-30 9/450 2026-04-03 18:34 by ls刘帅
[考研] 275求调剂 +14 waltzh 2026-04-01 14/700 2026-04-03 11:24 by qqq850810
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 22408 266求调剂 +3 masss11222 2026-04-02 3/150 2026-04-02 18:11 by 笔落锦州
[考研] 一志愿上海海洋大学083200食品学硕,求调剂,接受其他专业 +6 what张 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 环境工程 085701,267求调剂 +15 minht 2026-03-29 16/800 2026-04-01 10:13 by li_sujuan99
[考研] 一志愿大连理工大学,机械工程学硕,341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
信息提示
请填处理意见