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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 结构优化不收敛
71/1返回列表
查看: 2467  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zxgxbh2008

铜虫 (小有名气)

[求助] 结构优化不收敛

在做硅表面吸附一个有机分子时为什么结构优化总是不收敛呢，是什么问题呢
以下是我的输入文件，请帮我分析一下
# SCF input for VASP
# Note that VASP uses the FIRST occurence of a keyword
# Extra input added by hand
#RWIGS=1.11 0.77 0.32 0.73
LREAL=Auto
NPAR=1
EDIFFG=-0.05
# End of extra input

    PREC = Normal
   ENCUT = 400.000
  IBRION = 2
     NSW = 100
    ISIF = 2
    ALGO = Normal (blocked Davidson)
    NELM = 200
  NELMIN = 2
   EDIFF = 1.0e-04
  EDIFFG = -0.02
VOSKOWN = 1
  NBLOCK = 1
   ISPIN = 1
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
   LWAVE = .FALSE.
  LCHARG = .FALSE.
ADDGRID = .FALSE.
  ISMEAR = -5
   SIGMA = 0.2
   LREAL = .FALSE.
RWIGS = 1.11 0.77 0.32 0.73
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-09-02 19:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)

【答案】应助回帖

★
mazuju028(金币+1): 谢谢交流 2011-09-04 13:48:11
zxgxbh2008(金币+10): 2011-10-10 21:46:01
把NSW调大点 1000左右
赞一下(1人) 回复此楼
相逢是缘！
2楼2011-09-04 10:27:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengnan0501

铜虫 (小有名气)
★
uuv2010(金币+1): 欢迎参加交流 2011-09-06 22:34:18
是电子步不收敛还是离子步不收敛啊？？
赞一下(1人) 回复此楼
3楼2011-09-06 02:05:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wvsw

木虫 (小有名气)

老虎大王

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。正解。 2011-09-07 09:18:12
addgrid=true，ALGO=F;
你是刚开始就把精度设这么高吗？试试逐次提高精度
赞一下(1人) 回复此楼
生于忧患，死于安乐
4楼2011-09-07 08:44:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxgxbh2008

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by mengnan0501 at 2011-09-06 02:05:30:
是电子步不收敛还是离子步不收敛啊？？

就是算完了给提示结果不收敛，什么是电子和离子收敛啊？
赞一下 回复此楼
5楼2011-10-10 21:47:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

repent03

金虫 (小有名气)
★ ★
mazuju028(金币+2): 谢谢交流 2011-10-16 22:12:00
在做scf计算的时候
首先INCAR中友两个参数是关于电子步和离子步的，分别是NELM 和 NSW
计算状况是，电子步先动，然后才是离子步
你可以grep Iteration OUTCAR来看
（）括号中间的是电子步  Iteration 后面数字是离子步
两个的收敛条件预设是10E-4 10E-3
希望对你有帮助
赞一下(2人) 回复此楼
6楼2011-10-16 21:10:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxgxbh2008

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by repent03 at 2011-10-16 21:10:36:
在做scf计算的时候
首先INCAR中友两个参数是关于电子步和离子步的，分别是NELM 和 NSW
计算状况是，电子步先动，然后才是离子步
你可以grep Iteration OUTCAR来看
（）括号中间的是电子步  Iteration 后面数字 ...

POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.79942      0.01772     10.91379         0.002725     -0.002334      0.000618
      2.87011      0.00265      9.83896         0.002108      0.000579      0.003820
      1.44975      0.32883      5.40434        -0.009275      0.002769      0.000619
      1.89524     -0.32483     14.48545         0.014190     -0.016980     -0.010888
      0.86669      1.85559     12.31925        -0.002757     -0.000945      0.000212
      2.89491      2.02419      8.64081         0.010740      0.005914      0.011541
      1.04584      1.95602      7.12113         0.008514     -0.009898     -0.005415
      2.78915      1.74759     13.69920         0.004465     -0.007269     -0.005093
      0.99334      3.91994     11.21190        -0.000643      0.005689      0.015004
      2.98118      3.94592     10.01408         0.006552      0.010938      0.000995
      2.02798      3.92298      6.33231         0.018025     -0.013743     -0.003856
      2.16687      3.55373     15.11104        -0.001723     -0.001289     -0.010739
      1.60983      7.54184     13.18964        -0.004951     -0.002092     -0.001196
      2.87969      5.82854      8.59095         0.004570     -0.001279     -0.007901
      0.99783      5.85381      7.07512         0.000363     -0.011839     -0.005162
      2.98882      5.66578     14.48031        -0.020117     -0.007669     -0.026067
      0.78541      8.06039     11.02779         0.014696      0.007660     -0.003816
      2.74210      7.87728      9.75987         0.004023      0.001639      0.000892
      2.06255      7.74618      6.17885         0.002597     -0.006584     -0.004479
      2.05919      7.14815     16.21688         0.119920      0.128884     -0.007137
     -0.08417      9.98529     12.13861        -0.004977     -0.006067     -0.002880
      2.88159      9.66631      8.26435        -0.010094      0.008291      0.006930
      0.92638      9.68093      6.82296        -0.001072     -0.007392      0.015014
      3.49939      9.14335     13.93827         0.005112      0.004849      0.012649
      0.67753     11.67976     10.67033        -0.002361      0.006425      0.006691
      2.76080     11.57718      9.60403         0.009187     -0.000569      0.004467
      1.92799     11.66782      6.07117         0.006051     -0.003171     -0.009511
      2.61460     11.09928     15.15498        -0.029751     -0.038062     -0.000271
      0.87803     13.41249     12.24410        -0.000018      0.011116     -0.011822
      2.82256     13.53185      8.28444         0.012561     -0.001327      0.011056
      0.98887     13.71150      6.75107         0.011161      0.001746      0.000990
      2.88016     13.01250     13.61096         0.003708      0.008356     -0.006002
      4.83712     -0.03409     11.16295        -0.007920      0.004322      0.003364
      6.56323      0.01794      9.51283         0.002265      0.010666      0.000981
      3.68289      0.18094      6.12939         0.002889      0.005088      0.007516
      7.36525     14.97119     15.30446         0.010848      0.000060     -0.009351
      4.89573      1.84478     12.62419         0.005158      0.003112     -0.000169
      6.63671      2.06457      8.32494        -0.001329      0.014627      0.011012
      4.63102      2.14020      7.02967        -0.004176     -0.000448     -0.006696
      6.84512      1.65776     13.97124        -0.007568      0.015111      0.009145
      4.89706      3.92782     11.41823        -0.004485     -0.012777      0.007589
      6.70216      3.96519      9.81774        -0.003959      0.001678     -0.010973
      4.20192      3.93964      5.52032         0.001377      0.002484     -0.004485
     -0.26253      3.73100     15.06953         0.153676     -0.141692     -0.182662
      4.69289      5.92445     12.79035        -0.033352     -0.031289     -0.027971
      6.59445      5.77268      8.25065        -0.005021     -0.003039     -0.007674
      4.58602      5.82502      6.94821        -0.013478      0.009968      0.013418
      0.92871      4.99210     13.31706        -0.009873      0.011219      0.021944
      4.65449      7.67379     11.07285         0.018349     -0.006302     -0.010247
      6.54089      7.68311      9.59964        -0.003953      0.001072      0.004742
      4.23056      7.61788      5.37210         0.005661      0.006677      0.007342
      0.32647      8.30660     15.00151        -0.048855     -0.016118      0.002843
      5.11027      9.65201     12.21493        -0.002096     -0.000163      0.007527
      6.52719      9.58392      8.05635         0.005906      0.005122      0.004569
      4.59004      9.60553      6.66580         0.005247      0.006909      0.001141
      2.25899     10.74450     12.68542         0.013224      0.006694     -0.001270
      4.72902     11.55761     10.88172         0.001447      0.002421     -0.006773
      6.48143     11.63254      9.24067        -0.009952      0.005700      0.009088
      4.10959     11.51694      5.25718         0.014251     -0.002525      0.000058
      0.32625     10.73089     14.44958         0.013288      0.008683      0.004467
      4.95087     13.27916     12.49089        -0.002672     -0.017045      0.003955
      6.56169     13.55675      7.90114        -0.010890      0.007436     -0.016190
      4.57778     13.40657      6.63870        -0.003631     -0.005576     -0.004300
      6.91248     13.02670     13.89229         0.001234      0.015876     -0.005824
      6.15612      6.76168     13.58812         0.132002      0.114700      0.087658
      6.37731      7.31639     14.80731         0.036158     -0.029339     -0.190193
      3.11454      8.02822     16.82618        -0.103077     -0.089918     -0.036851
      6.89792      6.98439     12.79584        -0.064700     -0.012425      0.050814
      5.70040      7.13363     15.64458        -0.054112     -0.005961      0.062935
      1.54244      6.20583     17.25288        -0.049386     -0.079979      0.044264
      6.51156      4.50042     16.09991        -0.152021      0.128628      0.185995
-----------------------------------------------------------------------------------
    total drift:                               -0.024830     -0.015511      0.006404


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



  FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  ---------------------------------------------------
  free  energy   TOTEN  =      -364.133653 eV

  energy  without entropy=     -364.053056  energy(sigma->0) =     -364.093355

d Force =-0.8499241E-03[-0.289E-02, 0.119E-02]  d Energy =-0.8957913E-03 0.459E-04
d Force =-0.2284160E+01[-0.229E+01,-0.228E+01]  d Ewald  =-0.2284159E+01-0.175E-05

这是我最后一步的计算结果，麻烦你帮我看一下，谢谢
赞一下 回复此楼
7楼2011-10-17 08:34:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zxgxbh2008 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭