24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3022  |  回复: 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的回帖

春华秋实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的回帖
查看全部 8 个回答

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

金虫 (正式写手)

我的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的回帖
信息提示
请填处理意见