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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于晶体优化的困难
51/1返回列表
查看: 400  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 405963512 的 4 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

405963512

铜虫 (正式写手)

[交流] 【求助】关于晶体优化的困难

大家好，我说一下，自己的情况，我做的是有机晶体的能带计算，晶胞内原子也比较多，比如这次，114个（并且应该存在金属间的弱相互作用）。INCAR如下：
SYSTEM=RINQOV
ENCUT=500
ISTART=0;ICHARG=2
EDIFF=1e-5;EDIFFG=-0.001
ISMEAR=-5;SIGMA=0.10
PREC=Accurate
LREAL=AUTO;IALGO=48
NSW=100;IBRION=2;ISIF=2
但是优化起来结果是这样：
    N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.763084401788E+03   -0.76308E+03   -0.23897E+02  4736   0.288E+01    0.759E+00
RMM:   2    -0.770146065305E+03   -0.70617E+01   -0.12625E+01  4922   0.796E+00    0.234E+01
RMM:   3    -0.763160670548E+03    0.69854E+01   -0.66953E+00  4968   0.540E+00    0.826E+00
RMM:   4    -0.761916499352E+03    0.12442E+01   -0.10953E+00  5136   0.185E+00    0.460E+00
RMM:   5    -0.762446243620E+03   -0.52974E+00   -0.16318E-01  4855   0.760E-01    0.794E+00
RMM:   6    -0.762232743839E+03    0.21350E+00   -0.15162E-01  4999   0.675E-01    0.735E+00
RMM:   7    -0.762248549758E+03   -0.15806E-01   -0.13834E-01  4768   0.706E-01    0.755E+00
RMM:   8    -0.762057779684E+03    0.19077E+00   -0.35203E-02  4826   0.443E-01    0.653E+00
RMM:   9    -0.762211969469E+03   -0.15419E+00   -0.30147E-02  4839   0.401E-01    0.761E+00
RMM:  10    -0.761999504985E+03    0.21246E+00   -0.25964E-02  4762   0.392E-01    0.651E+00
RMM:  11    -0.762169670911E+03   -0.17017E+00   -0.21575E-02  4772   0.373E-01    0.761E+00
RMM:  12    -0.762003600561E+03    0.16607E+00   -0.21977E-02  4769   0.368E-01    0.654E+00
RMM:  13    -0.762135883903E+03   -0.13228E+00   -0.21323E-02  4750   0.371E-01    0.756E+00
RMM:  14    -0.762011013261E+03    0.12487E+00   -0.21609E-02  4763   0.370E-01    0.657E+00
RMM:  15    -0.762121079598E+03   -0.11007E+00   -0.21543E-02  4783   0.373E-01    0.754E+00
RMM:  16    -0.762011353686E+03    0.10973E+00   -0.21623E-02  4751   0.371E-01    0.658E+00
RMM:  17    -0.762110773639E+03   -0.99420E-01   -0.21603E-02  4769   0.374E-01    0.752E+00
RMM:  18    -0.762011372508E+03    0.99401E-01   -0.21683E-02  4769   0.372E-01    0.658E+00
RMM:  19    -0.762103022471E+03   -0.91650E-01   -0.21666E-02  4763   0.375E-01    0.751E+00
RMM:  20    -0.762011607510E+03    0.91415E-01   -0.21744E-02  4761   0.373E-01    0.658E+00
RMM:  21    -0.762096706132E+03   -0.85099E-01   -0.21729E-02  4786   0.376E-01    0.750E+00
RMM:  22    -0.762012171324E+03    0.84535E-01   -0.21812E-02  4755   0.374E-01    0.658E+00
RMM:  23    -0.762091313366E+03   -0.79142E-01   -0.21781E-02  4769   0.376E-01    0.749E+00
RMM:  24    -0.762012945985E+03    0.78367E-01   -0.21889E-02  4781   0.374E-01    0.658E+00
RMM:  25    -0.762086615040E+03   -0.73669E-01   -0.21834E-02  4764   0.377E-01    0.748E+00
RMM:  26    -0.762013865929E+03    0.72749E-01   -0.21941E-02  4760   0.375E-01    0.659E+00
RMM:  27    -0.762082396477E+03   -0.68531E-01   -0.21886E-02  4792   0.377E-01    0.747E+00
RMM:  28    -0.762014884063E+03    0.67512E-01   -0.21999E-02  4752   0.375E-01    0.659E+00
RMM:  29    -0.762078552598E+03   -0.63669E-01   -0.21920E-02  4771   0.378E-01    0.747E+00
RMM:  30    -0.762015965279E+03    0.62587E-01   -0.22052E-02  4786   0.376E-01    0.659E+00
RMM:  31    -0.762075021687E+03   -0.59056E-01   -0.21965E-02  4761   0.378E-01    0.746E+00
RMM:  32    -0.762017092387E+03    0.57929E-01   -0.22088E-02  4765   0.376E-01    0.659E+00
RMM:  33    -0.762071752003E+03   -0.54660E-01   -0.22008E-02  4789   0.378E-01    0.746E+00
RMM:  34    -0.762018222886E+03    0.53529E-01   -0.22129E-02  4760   0.376E-01    0.660E+00
RMM:  35    -0.762068692903E+03   -0.50470E-01   -0.22036E-02  4771   0.378E-01    0.745E+00
RMM:  36    -0.762019363113E+03    0.49330E-01   -0.22169E-02  4785   0.377E-01    0.660E+00
RMM:  37    -0.762065815196E+03   -0.46452E-01   -0.22073E-02  4764   0.378E-01    0.745E+00
RMM:  38    -0.762020498774E+03    0.45316E-01   -0.22194E-02  4768   0.377E-01    0.660E+00
RMM:  39    -0.762063104633E+03   -0.42606E-01   -0.22108E-02  4790   0.378E-01    0.744E+00
RMM:  40    -0.762021629829E+03    0.41475E-01   -0.22224E-02  4755   0.377E-01    0.660E+00
RMM:  41    -0.762060515406E+03   -0.38886E-01   -0.22130E-02  4783   0.378E-01    0.744E+00
RMM:  42    -0.762022757288E+03    0.37758E-01   -0.22258E-02  4782   0.377E-01    0.661E+00
RMM:  43    -0.762058059447E+03   -0.35302E-01   -0.22159E-02  4761   0.378E-01    0.743E+00
RMM:  44    -0.762023858202E+03    0.34201E-01   -0.22276E-02  4774   0.378E-01    0.661E+00
RMM:  45    -0.762055719519E+03   -0.31861E-01   -0.22190E-02  4792   0.379E-01    0.743E+00
RMM:  46    -0.762024958392E+03    0.30761E-01   -0.22300E-02  4760   0.378E-01    0.661E+00
RMM:  47    -0.762055759273E+03   -0.30801E-01   -0.22173E-02  4774   0.379E-01    0.743E+00
RMM:  48    -0.762028211859E+03    0.27547E-01   -0.22217E-02  4786   0.377E-01    0.662E+00
RMM:  49    -0.762052441269E+03   -0.24229E-01   -0.22204E-02  4767   0.379E-01    0.742E+00
RMM:  50    -0.762030938571E+03    0.21503E-01   -0.22202E-02  4769   0.377E-01    0.662E+00
RMM:  51    -0.762049432743E+03   -0.18494E-01   -0.22227E-02  4795   0.379E-01    0.742E+00
RMM:  52    -0.762022999932E+03    0.26433E-01   -0.22544E-02  4759   0.379E-01    0.660E+00
RMM:  53    -0.762048151006E+03   -0.25151E-01   -0.22481E-02  4778   0.381E-01    0.742E+00
RMM:  54    -0.762024317265E+03    0.23834E-01   -0.22606E-02  4785   0.380E-01    0.661E+00
RMM:  55    -0.762045481329E+03   -0.21164E-01   -0.22499E-02  4758   0.381E-01    0.741E+00
RMM:  56    -0.762025509814E+03    0.19972E-01   -0.22613E-02  4776   0.380E-01    0.661E+00
RMM:  57    -0.762042722358E+03   -0.17213E-01   -0.22559E-02  4783   0.381E-01    0.741E+00
RMM:  58    -0.762026589221E+03    0.16133E-01   -0.22637E-02  4760   0.380E-01    0.661E+00
RMM:  59    -0.762040544381E+03   -0.13955E-01   -0.22596E-02  4774   0.381E-01    0.740E+00
RMM:  60    -0.762027666326E+03    0.12878E-01   -0.22649E-02  4786   0.380E-01
   3 F= -.76202767E+03 E0= -.76202767E+03  d E =0.360103E+00
curvature:  -0.12 expect dE=-0.376E+00 dE for cont linesearch -0.157E+00

我设置了100步，但是他60步就停了，再有，看样子能量差值基本保持不变，知识上下摆动，
所以请大家帮忙看下，有什么解决办法没有，

谢谢大家了！

[ Last edited by 405963512 on 2010-2-6 at 19:59 ]
回复此楼

» 猜你喜欢
1楼 2010-02-06 19:58:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice_rain

铁杆木虫 (正式写手)
★
aylayl08(金币+1):感谢交流 2010-02-08 09:14
引用回帖:
Originally posted by 405963512 at 2010-02-06 19:58:10:
大家好，我说一下，自己的情况，我做的是有机晶体的能带计算，晶胞内原子也比较多，比如这次，114个（并且应该存在金属间的弱相互作用）。INCAR如下：
SYSTEM=RINQOV
ENCUT=500
ISTART=0;ICHARG=2
EDIFF=1e-5 ...

你先好好看看什么电子步，什么离子步吧。
赞一下(1人) 回复此楼
4楼2010-02-06 22:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

calos818

木虫 (著名写手)

新人菜鸟
★
aylayl08(金币+1):感谢交流 2010-02-08 09:13
提高ENCUT，降低NSW 呵呵，我也不懂，但帮你顶个贴吧呵呵
赞一下(1人) 回复此楼
2楼2010-02-06 20:18:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)
★ ★ ★
405963512(金币+1): 2010-02-07 16:14
aylayl08(金币+3):感谢交流 2010-02-08 09:14
ENCUT=500，应该差不多了吧！VASP的计算不是设置的参数越多越好，VASP有很多默认的参数，比如你的PREC=Accurate就已经可以接受了。在有，你设置的ISIF=2
我忘记不计算说明了，这个最好是小心设置，还有NSW=100，也差不多吧，一般的relaxation计算，只要不是和平衡态差的太多，基本上还是可以的，VASP会在你能量误差达到精度的时候停止的。但是你说的:震荡,肯定不好！还是查查吧！
赞一下(1人) 回复此楼
3楼2010-02-06 22:26:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

405963512

铜虫 (正式写手)
引用回帖:
Originally posted by ice_rain at 2010-02-06 22:30:01:

你先好好看看什么电子步，什么离子步吧。

你能解释一下，你是什么意思吗
赞一下 回复此楼
5楼2010-02-07 16:14:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭