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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么力没有达到收敛标准，就提示已经达到了收敛精度，并且停止了呢？
413/5返回列表上一页12345下一页
查看: 3789  |  回复: 40
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

c03185

新虫 (小有名气)
引用回帖:
19楼: Originally posted by 贺仪 at 2012-08-05 23:58:21
你把你的ediff提高到1E-7 然后重新写一行
EDIFFG= -0.002
再试试？...

还是那样的，无语了，哎。。。
赞一下 回复此楼
21楼2012-08-06 00:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
送鲜花一朵
引用回帖:
14楼: Originally posted by cenwanglai at 2012-08-05 21:02:08
分析一下你的几个原子最后最后一列的受力数据，你看看。受力是成对的，一个正A，一个负A，两相抵消恰好是0.

为什么不试着把某层原子全部固定呢？

能把POSCAR贴出来不？也可以直接发给我。
cwlzju@zju.edu.cn....

结果还是那样，找不到原因。同样的参数用在其他的晶格结构里边，驰豫原子对应的力是收敛的。在这个晶格结构就是不对。版主爷爷给指条活路吧。。。

INCAR：
SYSTEM = AO
ENCUT = 520
ISTART = 0; ICHARG = 2
ISMEAR = 1; SIGMA = 0.20
NSW = 600; IBRION = 2
POTIM = 0.1
EDIFF = 1E-7
EDIFFG = -5E-2
PREC= Accurate
IALGO=48

POSCAR：
AO(001)_9
   1.00000000000000
     4.2487171130865562   -2.9342791501140022    0.0000000000000010
     0.0000000000000000    4.5431245416429560   -2.4538970306330690
     6.8981356599745016    9.9882204546856510   18.4921082297066661
   O    Ti
   9   6
Selective dynamics
Direct
  0.9790666575233118  0.3333333333333357  0.0505438132934409   F   F   T
  0.6457333241899832  0.6666666666666643  0.2713454927088639   F   F   T
  0.6666666666666643  0.6457333241899832  0.0505438132934409   F   F   T
  0.3333333333333357  0.9790666575233118  0.2713454927088639   F   F   T
  0.3542666758100168  0.0209333424766882  0.0505438132934409   F   F   T
  0.0209333424766882  0.3542666758100168  0.2713454927088639   F   F   T
  0.6876000091433525  0.0000000000000000  0.1609446530011525   F   F   F
  0.0000000000000000  0.6876000091433525  0.1609446530011525   F   F   F
  0.3123999908566475  0.3123999908566475  0.1609446530011525   F   F   F
  0.6666666666666643  0.3333333333333357  0.1168738138540851   F   F   F
  0.3333333333333357  0.6666666666666643  0.3154820665970479   F   F   T
  0.6666666666666643  0.3333333333333357  0.0064072394052569   F   F   T
  0.3333333333333357  0.6666666666666643  0.2050154921482271   F   F   F
  0.0000000000000000  0.0000000000000000  0.0963564137588521   F   F   T
  0.0000000000000000  0.0000000000000000  0.2255328922434601   F   F   T


最终结果:
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.50830     -0.85384      0.11633        -0.671847     -0.681135      0.568972
      4.61545      3.84444      3.38219         0.399076      0.869518     -0.568972
      3.18100      1.48211     -0.65027         1.032169     -0.327554     -0.257660
      3.28815      6.18039      2.61559        -1.032169      0.035882      0.415202
      1.85370     -0.43977      0.88293        -0.399076      0.952574     -0.415202
      1.96085      4.25851      4.14878         0.671847     -0.849285      0.257660
      4.03164     -0.41006      2.97621         0.308830     -0.213287      0.000000
      1.11022      4.73140      1.28891         0.000000      0.330230     -0.178369
      2.43752      2.11015      2.20961        -0.308830     -0.116944      0.178369
      3.63869      0.72555      1.34328        -0.005537     -0.008018     -0.014844
      3.59252      5.20183      4.19812         0.005220      0.007559      0.013995
      2.87663     -0.37788     -0.69961        -0.005220     -0.007559     -0.013995
      2.83046      4.09840      2.15524         0.005537      0.008018      0.014844
      0.66483      0.96265      1.78224         0.011722      0.016973      0.031424
      1.55561      2.25245      4.17017        -0.011722     -0.016973     -0.031424
-----------------------------------------------------------------------------------
    total drift:                                0.000075      0.000108      0.000200
d Force = 0.1619927E-03[ 0.151E-03, 0.173E-03]  d Energy =-0.4911406E-03 0.653E-03
d Force = 0.2019964E+01[ 0.202E+01, 0.202E+01]  d Ewald  = 0.2019964E+01-0.245E-07
--------------------------------------------------------------------------------------------------------
    POTLOK:  cpu time    1.34: real time    0.69
--------------------------------------------------------------------------------------------------------
Steepest descent step on ions:
trial-energy change:    0.000491  1 .order   -0.000162   -0.000173   -0.000151
  (g-gl).g = 0.173E-03      g.g   = 0.173E-03  gl.gl    = 0.000E+00
g(Force)  = 0.173E-03   g(Stress)= 0.000E+00 ortho     = 0.000E+00
gamma     =   0.00000
trial     =   1.00000
opt step  =   4.00000  (harmonic =   7.63264) maximal distance =0.00162217
next E    =  -127.127400   (d E  =  -0.00066)
--------------------------------------------------------------------------------------------------------
reached required accuracy - stopping structural energy minimisation
writing wavefunctions
赞一下 回复此楼
22楼2012-08-06 01:04:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xyzzz986

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2012-08-07 09:20:47
c03185: 金币+10, ★有帮助 2012-08-15 07:32:30
1、换个优化的方法，比如IBRION=1，或者改变POTIM的值应该可以。

2、氧化物一般都不好算，特别是极化比较厉害的。

3、算TiO2表面还要加自旋的。
赞一下 回复此楼
23楼2012-08-06 02:00:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
引用回帖:
23楼: Originally posted by xyzzz986 at 2012-08-06 02:00:30
1、换个优化的方法，比如IBRION=1，或者改变POTIM的值应该可以。

2、氧化物一般都不好算，特别是极化比较厉害的。

3、算TiO2表面还要加自旋的。

谢谢，请问如何判断加不加自旋啊？
赞一下 回复此楼
24楼2012-08-06 04:46:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者
我也出现过这样的问题，知道的大侠赶紧帮帮啊！！！！！！！！！
赞一下 回复此楼
臾兮，臾兮，何时是归期！
25楼2012-08-06 09:04:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
★ ★
liliangfang: 金币+2, 谢谢指教 2012-08-07 09:20:59
怀疑问题不是出在算法上，而是作业输入文件有可能写错了，建议认真检查输入参数和输入文件
如果书写正确的话，不应该出现这样的问题
EDIFFG is negative it has a different meaning: In this case the relaxation will stop if all forces are smaller than |EDIFFG |
如果可以的话，可以将INCAR,POSCAR和KPOINTS的内容都贴上来，让大家一起找找原因

[ Last edited by uuv2010 on 2012-8-6 at 11:33 ]
赞一下(1人) 回复此楼
26楼2012-08-06 11:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xyzzz986

木虫 (著名写手)

【答案】应助回帖

★ ★
liliangfang: 金币+2, 谢谢指教 2012-08-07 09:21:09
1、我记得jcp上面曾经出过一个计算的专刊，大概是2006-2008年之间的事儿。里面有篇关于TiO2计算的综述，该大牛提到过TiO2的表面应该是有自旋的，但绝大多数作TiO2表面计算的人都没有考虑这样的效应。

2、你考虑自旋作个自洽计算，不考虑自旋作个自洽计算，然后比较一下。TiO2中Ti都是6配位的，O都是3配位的，做成表面必然有悬挂键，对于Ti，那些悬挂键可能是单电子。
赞一下 回复此楼
27楼2012-08-06 11:48:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2012-08-07 09:21:18
尝试检查OUTCAR里写的EDIFFG是不是确实是你需要的值。如果不是，用vim手打一遍INCAR吧...


emacs爱好者莫喷俺啊...俺只是路过
赞一下 回复此楼
28楼2012-08-06 19:57:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
引用回帖:
14楼: Originally posted by cenwanglai at 2012-08-05 21:02:08
分析一下你的几个原子最后最后一列的受力数据，你看看。受力是成对的，一个正A，一个负A，两相抵消恰好是0.

为什么不试着把某层原子全部固定呢？

能把POSCAR贴出来不？也可以直接发给我。
cwlzju@zju.edu.cn....

请问您确定一个正的一个负的相互抵消就是对的吗？
貌似我不管怎么修改参数，最终结果都是一个正，对应的一个负值。
赞一下 回复此楼
29楼2012-08-07 16:09:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
c03185: 金币+10, ★有帮助 2012-08-15 07:33:10
liliangfang: 金币+1, 谢谢交流 2012-08-16 21:46:52
引用回帖:
29楼: Originally posted by c03185 at 2012-08-07 16:09:24
请问您确定一个正的一个负的相互抵消就是对的吗？
貌似我不管怎么修改参数，最终结果都是一个正，对应的一个负值。...

我觉得不对~

理由不充分~
赞一下 回复此楼
30楼2012-08-07 17:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 c03185 的主题更新
413/5返回列表上一页12345下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂，求导师收 +6 天天好运来上岸� 2026-03-24 6/300 2026-03-24 17:19 by dick_runner
[考研] 收08调剂生 +4 komorebi69 2026-03-18 4/200 2026-03-24 16:53 by 不知名小学生
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 328求调剂，英语六级551，有科研经历 +7 生物工程调剂 2026-03-17 12/600 2026-03-23 18:18 by YMU施老师
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
信息提示
关闭
请填处理意见
关闭