24小时热门版块排行榜    

查看: 5135  |  回复: 11

yangyangwlhx

木虫 (正式写手)

少校

[求助] CASTEP几何结构优化的问题 虽然收敛了 感觉不是很对 请大家帮忙找一下原因 已有1人参与

我的体系用castep几何结构优化时,从图中看只有energy和stress这俩条线达到收敛,而force和displacement这俩条 线  压根就只有直线  。从几何结构优化的输出文件中看,也没有force 收敛的情况,我就想搞清楚 我的问题出在哪里了 是不是几何结构优化时 漏了什么参数 我先建了模型 然后做ECUT和kpoint的测试 再几何结构优化  顺便想问一下 有没有使用LDA+U方法的前辈 在计算过程中 加U 是在几何结构优化的时候 就加 还是在几何结构优化的时候不加,等优化结束以后,在计算能带和态密度以及其他性质时加U进行计算   请前辈指点 小弟新手  也没有人可问  只好在这里提一些粗浅的问题 麻烦大家了
CASTEP几何结构优化的问题 虽然收敛了 感觉不是很对 请大家帮忙找一下原因
DWJA8SHT$8V%8Y~E%WZ65WL.jpg


CASTEP几何结构优化的问题 虽然收敛了 感觉不是很对 请大家帮忙找一下原因-1
WJK8`FO9MKK@S0E9X`VBG}C.jpg


CASTEP几何结构优化的问题 虽然收敛了 感觉不是很对 请大家帮忙找一下原因-2
2CIOK_QKPVI{$V9`}N9~})K.jpg
回复此楼

» 猜你喜欢

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

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-03-24 17:13:41
第一步计算表明你的force和displacement已经远远小于收敛标准,程序就不考虑这个了。从图上也可以看出来,你的Max force、Max displacement其实再图的最低端(-6的位置),之所以这样是因为他们已经低于程序默认的限度。这个结构没有问题,如果你不相信,我建议你提升scf、displacement和force的收敛标准,但是我自己的经验来说这是无用功,换句话说,其实已经收敛了。

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

不一定挂在论坛,计算问题问题欢迎留言。
2楼2013-11-30 23:33:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流 2013-12-02 22:12:20
引用回帖:
3楼: Originally posted by yangyangwlhx at 2013-12-01 12:20:30
非常感谢您的指点  请问您计算的体系中 有没有使用用LDA+u的方法 U是在什么阶段加 是在几何结构优化时加U 还是在结构优化后 计算性质时加U  谢谢...

老实说,我一直用的是castep,也没怎么使用过+U,因为可能我做的体系有的时候不太适合+U,根本算不动。不过前两天看到有人说+U对结构影响不大。这个我不能确定,但是有一点,我们知道Hubbard U显然是包含在Hamilton量里面的那么对能量有影响,构造出来的势能面多多少少有区别。退一万步说,我们把U值划分到Exc中,不同的U其实就如同不同的Exc,做出来的结构多多少少会有些不同的吧。如果有必要你可以做个测试看看,比如U值0和某个值优化的结构去比较一下晶格常数。我觉得+U来进行几何优化肯定是一件相当耗时耗力的事情,所以往往做一些简化的处理。

这个是个人愚见,望各位指正。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2013-12-01 15:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-03-24 17:13:56
引用回帖:
7楼: Originally posted by yangyangwlhx at 2013-12-01 20:53:56
加U 计算量猛然增大  比较费劲...

所以绝大部分人似乎都是直接优化,结构靠谱了再去+U。类似的案例也有,我之前听我师兄说他看到有人先用Dmol3优化,再CASTEP计算单点能,同样我也见过PBE优化HSE计算能带的。我觉得有两点要考虑,对于发文章要说服审稿人,从科研角度我觉得怎么说明这种“近似”的合理性比较重要。不过castep确实很慢。
不一定挂在论坛,计算问题问题欢迎留言。
8楼2013-12-01 22:39:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yangyangwlhx

木虫 (正式写手)

少校

送红花一朵
引用回帖:
2楼: Originally posted by 卡开发发 at 2013-11-30 23:33:52
第一步计算表明你的force和displacement已经远远小于收敛标准,程序就不考虑这个了。从图上也可以看出来,你的Max force、Max displacement其实再图的最低端(-6的位置),之所以这样是因为他们已经低于程序默认的 ...

非常感谢您的指点  请问您计算的体系中 有没有使用用LDA+u的方法 U是在什么阶段加 是在几何结构优化时加U 还是在结构优化后 计算性质时加U  谢谢
孤狼
3楼2013-12-01 12:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangyangwlhx

木虫 (正式写手)

少校

引用回帖:
4楼: Originally posted by 卡开发发 at 2013-12-01 15:20:30
老实说,我一直用的是castep,也没怎么使用过+U,因为可能我做的体系有的时候不太适合+U,根本算不动。不过前两天看到有人说+U对结构影响不大。这个我不能确定,但是有一点,我们知道Hubbard U显然是包含在Hamilto ...

谢谢楼主 关于加U  在什么时候加 众说纷纭 不知道谁的对 真的 有的说加U对晶格常数影响很大 有的说 没影响  我在小木虫上 看了好多关于LDA+U的贴  最终还是确定不下来 在那个阶段加U 。。。。。。。。。。。。。。
孤狼
5楼2013-12-01 17:43:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流 2013-12-02 22:12:37
引用回帖:
5楼: Originally posted by yangyangwlhx at 2013-12-01 17:43:32
谢谢楼主 关于加U  在什么时候加 众说纷纭 不知道谁的对 真的 有的说加U对晶格常数影响很大 有的说 没影响  我在小木虫上 看了好多关于LDA+U的贴  最终还是确定不下来 在那个阶段加U 。。。。。。。。。。。。。。...

个人觉得一开始就加原理上应该不会错,计算量上如何是个问题
不一定挂在论坛,计算问题问题欢迎留言。
6楼2013-12-01 20:20:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangyangwlhx

木虫 (正式写手)

少校

引用回帖:
6楼: Originally posted by 卡开发发 at 2013-12-01 20:20:17
个人觉得一开始就加原理上应该不会错,计算量上如何是个问题...

加U 计算量猛然增大  比较费劲
孤狼
7楼2013-12-01 20:53:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangyangwlhx

木虫 (正式写手)

少校

引用回帖:
8楼: Originally posted by 卡开发发 at 2013-12-01 22:39:11
所以绝大部分人似乎都是直接优化,结构靠谱了再去+U。类似的案例也有,我之前听我师兄说他看到有人先用Dmol3优化,再CASTEP计算单点能,同样我也见过PBE优化HSE计算能带的。我觉得有两点要考虑,对于发文章要说服审 ...

CASTEP 之前的版本上 在结构优化时是不能加U 只能在结构优化完了 再去加U  但是5.5和6.0在结构优化时 是可以加U 的  这就让人比较困惑  不过 我看一本有关材料计算的书上说 计算的流程是 先输入计算模型和一些物理参数-------判断晶体结构是否为稳定结构-----计算晶体的物理性质-----输出计算结果  如果按照这个流程   有一种说法就是不妥当的 那就是 在几何结构优化时加U并且计算材料的一些性质,我觉得应该是在几何结构优化的时候 加U 结构优化出来后 再计算晶体的物理性质
孤狼
9楼2013-12-02 09:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 谢谢回帖交流 2013-12-02 22:12:48
引用回帖:
9楼: Originally posted by yangyangwlhx at 2013-12-02 09:35:28
CASTEP 之前的版本上 在结构优化时是不能加U 只能在结构优化完了 再去加U  但是5.5和6.0在结构优化时 是可以加U 的  这就让人比较困惑  不过 我看一本有关材料计算的书上说 计算的流程是 先输入计算模型和一些物理 ...

后面两者没什么区别,其实无非就是多加了一下自己的判断,如果测试工作做得到位这些也都无所谓的,当然也可以参考别人的文献。
不一定挂在论坛,计算问题问题欢迎留言。
10楼2013-12-02 13:09:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangyangwlhx 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 8/400 2026-03-17 16:57 by ruiyingmiao
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 285化工学硕求调剂(081700) +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
请填处理意见