24小时热门版块排行榜    

查看: 4989  |  回复: 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 的主题更新
信息提示
请填处理意见