CyRhmU.jpeg
查看: 4080  |  回复: 20

stractor

金虫 (著名写手)


[交流] 【求助】如何看力

在VASP中,优化一般以力为收敛标准,怎么看力呢?
我在OUTCAR中发现了很多关于力的,比如
FORCE on cell =-STRESS in cart. coord.  units (eV/reduce length):
  Direction    X        Y        Z        XY       YZ       ZX
  --------------------------------------------------------------------------------------
  Alpha Z    228.09    228.09    228.09

。。。。。。

Total       -3.54     -4.54     13.70     -4.58     -0.30      0.42
  in kB       -5.41     -6.94     20.92     -6.99     -0.46      0.64
  external pressure =        2.86 kB  Pullay stress =        0.00 kB

。。。。。。

FORCES acting on ions
    electron-ion (+dipol)            ewald-force                    non-local-force                 convergence-correction
-----------------------------------------------------------------------------------------------
   -.397E+00 0.126E-01 -.195E+00   0.345E+00 0.378E-02 0.191E+00   0.492E-01 0.508E-02 0.287E-02   0.149E-02 0.121E-02 -.427E-03
   -.279E+00 0.448E+00 0.232E+00   0.268E+00 -.385E+00 -.192E+00   0.151E-01 -.731E-01 -.463E-01   0.142E-02 -.120E-02 -.535E-03

。。。。。。

0.720E+00 0.667E+00 0.351E+01   -.124E+00 -.874E+00 0.151E+02   0.326E-01 0.278E-01 -.409E-01   0.250E-03 0.143E-03 -.281E-03
-----------------------------------------------------------------------------------------------
   0.963E+00 -.298E+00 0.194E+02   0.185E-12 -.108E-13 -.799E-13   -.970E+00 0.301E+00 -.194E+02   0.211E-02 0.110E-02 -.395E-02


POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
     10.06599     10.06719      0.01673        -0.001720      0.022612     -0.001236
     10.06703      5.03749     10.34262         0.005402     -0.011136     -0.007025
。。。。。。
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

fanchen021

铁杆木虫 (正式写手)


stractor(金币+1):但 OUTCAR 里面没有这一项 2010-11-25 11:22:39
看作用在原子上的最大力
FORCES: max atom
2楼2010-11-25 09:46:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)


stractor(金币+1): 2010-11-26 12:19:49
你的EDIFFG=-1E-3吗?
3楼2010-11-25 10:54:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
4楼2010-11-25 10:57:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
stractor(金币+2):谢谢交流。可能需要在http://theory.cm.utexas.edu/vtsttools/bader/下载编译一些脚本。 2010-11-26 12:20:23
引用回帖:
Originally posted by fanchen021 at 2010-11-25 09:46:29:
看作用在原子上的最大力
FORCES: max atom

我在结构优化计算的OUTCAR中查过了,没有这一项了呢。
5楼2010-11-25 11:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)



youzhizhe(金币+1):谢谢参与交流。 2010-11-25 13:32:38
引用回帖:
Originally posted by fzuvivi at 2010-11-25 10:54:44:
你的EDIFFG=-1E-3吗?

EDIFFG=-0.03。
个人感觉是
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
     10.06491     10.06976      0.01770         0.010282      0.007834      0.002618

......................
不知道是不是?

[ Last edited by stractor on 2010-11-25 at 11:24 ]
6楼2010-11-25 11:21:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


引用回帖:
Originally posted by cenwanglai at 2010-11-25 11:08:45:


我在结构优化计算的OUTCAR中查过了,没有这一项了呢。

不会吧,我的所有OUTCAR里都有这一项啊
     3.51061      2.02685      5.67094         0.000000      0.000000      0.018799
      5.85101      2.02685      0.00000         0.000000      0.000000      0.715101
      4.68081      4.05370      1.65492         0.000000      0.000000     -0.774121
      3.51061      2.02685      3.30983         0.000000      0.000000     -0.118627
      5.85101      2.02685      4.90216         0.000000      0.000000      0.016457
-----------------------------------------------------------------------------------
    total drift:                                0.000000     -0.000001      0.012219


--------------------------------------------------------------------------------------------------------


  FORCES: max atom, RMS     0.026838    0.012964
  FORCE total and by dimension    0.073335    0.026838
OPT: skip step - force has converged
7楼2010-11-25 11:55:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by fanchen021 at 2010-11-25 11:55:09:


不会吧,我的所有OUTCAR里都有这一项啊
     3.51061      2.02685      5.67094         0.000000      0.000000      0.018799
      5.85101      2.02685      0.00000         0.000000      0.000000 ...

什么版本的vasp?

vasp4.6.xx?
vasp5.xx?
8楼2010-11-25 12:59:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


引用回帖:
Originally posted by cenwanglai at 2010-11-25 12:59:54:


什么版本的vasp?

vasp4.6.xx?
vasp5.xx?

vasp 4.6.36
9楼2010-11-25 13:10:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by stractor at 2010-11-25 11:21:51:

EDIFFG=-0.03。
个人感觉是
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
     10.0 ...

是。呵呵。不过我没高清楚,力的正负是什么意思?
10楼2010-11-25 13:31:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)


引用回帖:
Originally posted by fanchen021 at 2010-11-25 11:55:09:


不会吧,我的所有OUTCAR里都有这一项啊
     3.51061      2.02685      5.67094         0.000000      0.000000      0.018799
      5.85101      2.02685      0.00000         0.000000      0.000000 ...

我的5.2版本的,也没有你说的这一行啊 “  FORCES: max atom, RMS     0.026838    0.012964”
  FORCE total and by dimension    0.073335    0.026838
OPT: skip step - force has converged

我的只是在最后写一个total drift
11楼2010-11-25 14:56:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


vasp 5.2.2 我也用了,我优化的OUTCAR里也是有那一行的。
12楼2010-11-25 16:38:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)



youzhizhe(金币+1):谢谢提示。 2010-11-25 23:54:57
stractor(金币+2):谢谢提示 2010-11-26 12:17:08
我一般看收敛都是用vef.pl 这个脚本看的(vtstscripts里的)。这个脚本里就是在OUTCAR提取原子的最大受力的
13楼2010-11-25 16:42:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao72379

木虫 (正式写手)


stractor(金币+2):谢谢交流 2010-11-26 12:21:37
呵呵 这个问题我原来也是不懂  用了半年vasp也没关心着问题 主要功夫是看能量去了 嘿嘿 我给个正解吧(自认为正确的 哈哈)
在OUTCAR里是能找到如下信息的(4.6版)
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
     2.32720      9.15829      7.18156         0.002637     -0.004234      0.000292
     8.62320      9.15829      7.18156         0.002637     -0.004234      0.000292
     0.75579     11.92243      7.16200        -0.001375      0.001914     -0.000926
     7.05179     11.92243      7.16200        -0.001375      0.001914     -0.000926
     0.75649     17.47136      7.12124        -0.000949      0.002844      0.000262
     7.05249     17.47136      7.12124        -0.000949      0.002844      0.000262
     5.47484      9.15824      7.18139         0.001707     -0.003726     -0.000903
每排后3个数据就是力,一般要所有的力数据收敛到你INCAR设定的标准 relax就会结束了
14楼2010-11-25 16:54:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)


引用回帖:
Originally posted by fanchen021 at 2010-11-25 16:42:55:
我一般看收敛都是用vef.pl 这个脚本看的(vtstscripts里的)。这个脚本里就是在OUTCAR提取原子的最大受力的

我找到了这个脚本,但怎么用,在当前作业目录下采用 ./vef.pl 命令吗?
15楼2010-11-26 12:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


stractor(金币+1):谢谢解答,是不是要编译什么,我的不会出现你这样的结果,fe.dda文件什么内容都没有 2010-11-27 22:44:57
引用回帖:
Originally posted by stractor at 2010-11-26 12:19:24:

我找到了这个脚本,但怎么用,在当前作业目录下采用 ./vef.pl 命令吗?

是的。比如,我在我的一个优化好的目录下输入vef.pl 后,会出现,
    0           1.59483500          -216.043454                    0
    1           1.50700800          -216.063328            -0.019874
    2           1.39912000          -216.202084             -0.15863
    3           1.35256700          -216.214994             -0.17154
    4           0.97341100          -216.346551            -0.303097
    5           0.90999800          -216.357609            -0.314155
    6           0.74630000          -216.456104             -0.41265
    7           0.70556300          -216.463665            -0.420211
    8           0.60424300          -216.549739            -0.506285
    9           0.56631000          -216.556446            -0.512992
   10           0.34305100          -216.586833            -0.543379
   11           0.31455300          -216.591458            -0.548004
   12           0.55732000          -216.624470            -0.581016
   13           0.52821900          -216.628386            -0.584932
   14           0.28016700          -216.658036            -0.614582
   15           0.25022700          -216.661423            -0.617969
   16           0.27186400          -216.670846            -0.627392
   17           0.22158600          -216.673514             -0.63006
   18           0.15353700          -216.679686            -0.636232
   19           0.12014000          -216.681291            -0.637837
   20           0.20970900          -216.685369            -0.641915
   21           0.13338200          -216.686432            -0.642978
   22           0.08766300          -216.687047            -0.643593
   23           0.06567300          -216.687499            -0.644045
   24           0.07568400          -216.687493            -0.644039
   25           0.04350000          -216.688126            -0.644672
   26           0.07302700          -216.688421            -0.644967
   27           0.04070300          -216.688689            -0.645235
   28           0.06882400          -216.688672            -0.645218
   29           0.05552600          -216.689060            -0.645606
   30           0.05419000          -216.689335            -0.645881
   31           0.04710200          -216.689607            -0.646153
   32           0.03033700          -216.689594             -0.64614
   33           0.03087700          -216.689797            -0.646343
   34           0.03965800          -216.689809            -0.646355
   35           0.04456100          -216.689733            -0.646279
   36           0.02284800          -216.689808            -0.646354

[ Last edited by fanchen021 on 2010-11-26 at 14:21 ]
16楼2010-11-26 14:19:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)


★ ★
zzy870720z(金币+2):谢谢指教 2010-11-26 15:44:55
stractor(金币+1):估计是我没编译那里面的一些东西。看了上面的解释,说要重新编译VASP,对吗? 2010-11-27 22:46:27
其中第一列是优化步数,第二列是最大原子受力,小于收敛标准优化就结束了,第三列是总能,第四列是当前步总能和第一步总能的差值,可以看总能的变化趋势。

总的来说,UT的这个脚本看力还是很方便的
17楼2010-11-26 14:24:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1051358187

金虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
16楼: Originally posted by fanchen021 at 2010-11-26 14:19:37
是的。比如,我在我的一个优化好的目录下输入vef.pl 后,会出现,
    0           1.59483500          -216.043454                    0
    1           1.50700800          -216.063328            -0 ...

请教一下,为什么我用了这个脚本什么都没有输出来呢,输出一个空文件

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

18楼2013-05-29 21:30:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asiawan

铁虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
引用回帖:
18楼: Originally posted by 1051358187 at 2013-05-29 21:30:18
请教一下,为什么我用了这个脚本什么都没有输出来呢,输出一个空文件...

你好,请问你解决了vef.pl这个脚本输出空文件的问题了吗?我遇到了相同的问题,可否讲解一下
19楼2015-10-27 19:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asiawan

铁虫 (初入文坛)


送红花一朵
引用回帖:
18楼: Originally posted by 1051358187 at 2013-05-29 21:30:18
请教一下,为什么我用了这个脚本什么都没有输出来呢,输出一个空文件...

你好,请问你解决了vef.pl这个脚本输出空文件的问题了吗?我遇到了相同的问题,可否讲解一下
20楼2015-10-27 19:23:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

21楼2021-01-09 20:53:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 stractor 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见