应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何看力
81/1返回列表
查看: 4227  |  回复: 20
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
。。。。。。
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-11-25 09:37:04
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖

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提取原子的最大受力的
赞一下(1人) 回复此楼
13楼2010-11-25 16:42:55
已阅   回复此楼   关注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的这个脚本看力还是很方便的
赞一下(2人) 回复此楼
17楼2010-11-26 14:24:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 stractor 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭