| 查看: 429 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 405963512 的 4 个金币 | |||
[交流]
【求助】关于晶体优化的困难
|
|||
|
大家好，我说一下，自己的情况，我做的是有机晶体的能带计算，晶胞内原子也比较多，比如这次，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 ] |
回复此楼» 猜你喜欢
277求调剂 数一104分
已经有13人回复
-
304求调剂
已经有10人回复
-
求调剂
已经有5人回复
-
336求调剂，一志愿中科大
已经有6人回复
-
071000生物学，一志愿深圳大学296分，求调剂
已经有3人回复
-
一志愿北京化工085600 310分求调剂
已经有19人回复
-
274求调剂求调剂
已经有7人回复
-
277工科求调剂
已经有10人回复
-
283分求调剂
已经有10人回复
-
295求调剂
已经有16人回复
calos818
木虫 (著名写手)
新人菜鸟
- 应助: 11 (小学生)
- 金币: 3125.7
- 散金: 77
- 红花: 3
- 帖子: 1201
- 在线: 730.8小时
- 虫号: 698507
- 注册: 2009-02-08
- 专业: 凝聚态物性 II ：电子结构
2楼2010-02-06 20:18:12
3楼2010-02-06 22:26:57
ice_rain
铁杆木虫 (正式写手)
- 应助: 34 (小学生)
- 贵宾: 0.871
- 金币: 5590.5
- 散金: 60
- 红花: 6
- 帖子: 848
- 在线: 963.6小时
- 虫号: 283377
- 注册: 2006-10-08
- 专业: 凝聚态物性 II ：电子结构
4楼2010-02-06 22:30:01
5楼2010-02-07 16:14:34
