24小时热门版块排行榜    

查看: 1096  |  回复: 5

甲斐之虎

银虫 (正式写手)

[求助] 关于计算之后个别分子仍然受力很大请教各位前辈……

问题1:首先上INCAR,EDIFFG设置的是-5E-2,照理说这个标准不高,但是收敛极慢,计算起来都着急,不知道这个INCAR有什么问题没有?
问题2:不知道OUTCAR的这个TOTAL-FORCE怎么看,前三列是原子xyz坐标,后三列是原子在xyz方向受的力吗?
问题3:运行过程的二张截图,一个是最开始的,一个是最后的,可以发现从1F到27F,F的数值一直都是.296E+03,这个应该是原子受到最大的力吧,从头到尾一直都是这个数值,一直到收敛结束,不知道怎么回事……
最后,得到的结果与文献上相比,总是有5%的误差,不知道是不是这里出问题了
关于计算之后个别分子仍然受力很大请教各位前辈……
INCAR.jpg


关于计算之后个别分子仍然受力很大请教各位前辈……-1
OUTCAR-Total-Force.jpg


关于计算之后个别分子仍然受力很大请教各位前辈……-2
运行过程.jpg


关于计算之后个别分子仍然受力很大请教各位前辈……-3
运行过程2.jpg



[ 来自小组 石墨烯 ]
回复此楼

» 猜你喜欢

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

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

别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢交流 2013-10-28 07:36:50
这个收敛速度也跟体系有关吧,不知道先优化了没,第二个问题应该你说的是对的,这个total drift 目前还不知道是什么意思,第三个问题我也等待高手解答,包括DAV:1 这一行整个都需要解释,能力有限啊 。
希望有一天成为大神
2楼2013-10-27 15:19:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)

引用回帖:
2楼: Originally posted by sungjen at 2013-10-27 15:19:05
这个收敛速度也跟体系有关吧,不知道先优化了没,第二个问题应该你说的是对的,这个total drift 目前还不知道是什么意思,第三个问题我也等待高手解答,包括DAV:1 这一行整个都需要解释,能力有限啊 。

多谢参与啊,大侠,我们一起坐等大神出现吧
别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
3楼2013-10-27 16:05:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
甲斐之虎: 金币+30, ★★★★★最佳答案, 這位前輩太?�?�?�?�?姶罅� 2013-11-04 12:07:34
fzx2008: 金币+2, 谢谢指导 2013-11-13 09:30:52
问题1:首先上INCAR,EDIFFG设置的是-5E-2,照理说这个标准不高,但是收敛极慢,计算起来都着急,不知道这个INCAR有什么问题没有?

你指的"慢" 是指 "跑了很多ionic loops"
最簡單最直覺的答案就是 POSCAR "不好"
POSCAR 給的結構不好造成 relaxation 過程需要許多功夫
至於 INCAR 的問題~~不少
我有回信在你別的帖

问题2:不知道OUTCAR的这个TOTAL-FORCE怎么看,前三列是原子xyz坐标,后三列是原子在xyz方向受的力吗?
yes

问题3:运行过程的二张截图,一个是最开始的,一个是最后的,可以发现从1F到27F,F的数值一直都是.296E+03,这个应该是原子受到最大的力吧,从头到尾一直都是这个数值,一直到收敛结束,不知道怎么回事……

NO, F 不是受力
F 是能量, 同一行 E0 也是能量
F 是 free energy, E0 是 total energy
兩者之間差一個 entropy
4楼2013-11-04 02:30:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)

???????:
4?: Originally posted by jpchou at 2013-11-04 02:30:19
????1????????INCAR??EDIFFG???????-5E-2???????????????????????????????????????????????????????INCAR??????????У?

?????"??" ??? "??????ionic loops"
????????X?? ...

???x??????
?????????}????????????Y???D?????????e????????????????????a31(???????е???е????)??-0.126333???@?????????????????????????????????С?5E-2??EDIFFG=-5E-2???????N?????]???_?????????@????????
?????????}???????F???E0?????????????@?????}????^?^????F??????o?B???????ISMEAR=-5???@???]??entropy?????????@?????}
别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
5楼2013-11-04 12:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

★ ★
fzx2008: 金币+2, 谢谢回帖! 2013-11-13 10:52:01
???????:
5?: Originally posted by ?????? at 2013-11-04 12:14:52
???x??????
?????????}????????????Y???D?????????e????????????????????a31(???????е???е????)??-0.126333???@?????????????????????????????????С?5E-2??EDIFFG=-5E-2?????? ...

?????????}????????????Y???D?????????e????????????????????a31(???????е???е????)??-0.126333???@?????????????????????????????????С?5E-2??EDIFFG=-5E-2???????N?????]???_?????????@????????

????
?D????????
????????????~???????????
????
?@?? TOTAL-FORCE data ?????????? ionic loop
?????????? ionic loop ??


?????????}???????F???E0?????????????@?????}????^?^????F??????o?B???????ISMEAR=-5???@???]??entropy?????????@?????}

?? F or E0 ??????
???????
??? case A ?? F
case B ?? E0
?????у????????p
(??????g E0 ... ?]???κ????)

?????]????? MD ~ ???????? entropy ????}
??????????? ISMEAR -5 ??????]???P?B??
6楼2013-11-04 17:24:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 甲斐之虎 的主题更新
信息提示
请填处理意见