24小时热门版块排行榜    

查看: 3665  |  回复: 40

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 ]
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 的主题更新
信息提示
请填处理意见