24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3017  |  回复: 7

春华秋实2020

金虫 (正式写手)

[求助] VASP计算中total-force为什么等于零

vasp中弛豫石墨烯,但是结果中total-force等于零,不知道是怎么回事,请高手指点
INCAR文件如下:
ISIF=2
ISMEAR=0
SIGMA=0.1
ENCUT=500
IBRION=2
NSW=300
LVTOT=.TRUE.
#EDIFF=1E-5
EDIFFG=-0.05
LREAL=Auto
感觉没什么问题,网上有人说是对称性的问题,但我选的赝势是PAW_PBE,ISYM默认值是2,是考虑对称性的,难道不要考虑对称性。
回复此楼

» 猜你喜欢

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

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

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

171713294

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
春华秋实2020: 金币+2, 有帮助, 多谢 2013-05-15 11:58:18
sunyang1988: 金币+1, 谢谢交流 2013-06-13 00:08:33
应该不考虑其对称性,还有贴出你的outcar里面含有最后总理的那部分看看
坚持,不放弃
2楼2013-05-15 08:14:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.71014      1.23000      9.80000         0.000000      0.000000      0.000000
      1.42028      0.00000      9.80000         0.000000      0.000000      0.000000
      2.84056      0.00000      9.80000         0.000000      0.000000      0.000000
      3.55070     -1.23000      9.80000         0.000000      0.000000      0.000000
      4.97098     -1.23000      9.80000         0.000000      0.000000      0.000000
      5.68112     -2.46000      9.80000         0.000000      0.000000      0.000000
      7.10141     -2.46000      9.80000         0.000000      0.000000      0.000000
      7.81155     -3.69000      9.80000         0.000000      0.000000      0.000000
      9.23183     -3.69000      9.80000         0.000000      0.000000      0.000000
      9.94197     -4.92000      9.80000         0.000000      0.000000      0.000000
      0.71014      3.69000      9.80000         0.000000      0.000000      0.000000
      1.42028      2.46000      9.80000         0.000000      0.000000      0.000000
      2.84056      2.46000      9.80000         0.000000      0.000000      0.000000
      3.55070      1.23000      9.80000         0.000000      0.000000      0.000000
      4.97098      1.23000      9.80000         0.000000      0.000000      0.000000
      5.68112      0.00000      9.80000         0.000000      0.000000      0.000000
      7.10141      0.00000      9.80000         0.000000      0.000000      0.000000
      7.81155     -1.23000      9.80000         0.000000      0.000000      0.000000
      9.23183     -1.23000      9.80000         0.000000      0.000000      0.000000
      9.94197     -2.46000      9.80000         0.000000      0.000000      0.000000
      0.71014      6.15000      9.80000         0.000000      0.000000      0.000000
      1.42028      4.92000      9.80000         0.000000      0.000000      0.000000
      2.84056      4.92000      9.80000         0.000000      0.000000      0.000000
      3.55070      3.69000      9.80000         0.000000      0.000000      0.000000
      4.97098      3.69000      9.80000         0.000000      0.000000      0.000000
      5.68112      2.46000      9.80000         0.000000      0.000000      0.000000
      7.10141      2.46000      9.80000         0.000000      0.000000      0.000000
      7.81155      1.23000      9.80000         0.000000      0.000000      0.000000
      9.23183      1.23000      9.80000         0.000000      0.000000      0.000000
      9.94197      0.00000      9.80000         0.000000      0.000000      0.000000
      0.71014      8.61000      9.80000         0.000000      0.000000      0.000000
      1.42028      7.38000      9.80000         0.000000      0.000000      0.000000
      2.84056      7.38000      9.80000         0.000000      0.000000      0.000000
      3.55070      6.15000      9.80000         0.000000      0.000000      0.000000
      4.97098      6.15000      9.80000         0.000000      0.000000      0.000000
      5.68112      4.92000      9.80000         0.000000      0.000000      0.000000
      7.10141      4.92000      9.80000         0.000000      0.000000      0.000000
      7.81155      3.69000      9.80000         0.000000      0.000000      0.000000
      9.23183      3.69000      9.80000         0.000000      0.000000      0.000000
      9.94197      2.46000      9.80000         0.000000      0.000000      0.000000
      0.71014     11.07000      9.80000         0.000000      0.000000      0.000000
      1.42028      9.84000      9.80000         0.000000      0.000000      0.000000
      2.84056      9.84000      9.80000         0.000000      0.000000      0.000000
      3.55070      8.61000      9.80000         0.000000      0.000000      0.000000
      4.97098      8.61000      9.80000         0.000000      0.000000      0.000000
      5.68112      7.38000      9.80000         0.000000      0.000000      0.000000
      7.10141      7.38000      9.80000         0.000000      0.000000      0.000000
      7.81155      6.15000      9.80000         0.000000      0.000000      0.000000
      9.23183      6.15000      9.80000         0.000000      0.000000      0.000000
      9.94197      4.92000      9.80000         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------
    total drift:                                0.000000      0.000000      0.000000
3楼2013-05-15 11:56:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

我的INCAR是从师兄计算的一个32个原子的石墨烯的文件里直接copy的,他的没问题,
4楼2013-05-15 11:58:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

想到一个可能的原因,我师兄建立的石墨烯的初始结构里健长不是1.42,有1.43的,而我建立的石墨烯结构健长都是1.42,这应该理论上的理想模型了,会不会是因为这个,所以受力为零。
5楼2013-05-15 12:39:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)

引用回帖:
5楼: Originally posted by 春华秋实2020 at 2013-05-15 12:39:10
想到一个可能的原因,我师兄建立的石墨烯的初始结构里健长不是1.42,有1.43的,而我建立的石墨烯结构健长都是1.42,这应该理论上的理想模型了,会不会是因为这个,所以受力为零。

不知道这个total-force 是不是就是每个原子所受的力
希望有一天成为大神
6楼2013-06-11 14:26:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风云火海

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by 春华秋实2020 at 2013-05-15 12:39:10
想到一个可能的原因,我师兄建立的石墨烯的初始结构里健长不是1.42,有1.43的,而我建立的石墨烯结构健长都是1.42,这应该理论上的理想模型了,会不会是因为这个,所以受力为零。

我也遇到了这个问题,但是我觉得应该不是这个原因。你可以试一下。

发自小木虫Android客户端
7楼2019-01-12 07:54:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

8楼2022-06-19 15:12:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 春华秋实2020 的主题更新
信息提示
请填处理意见