优化收敛问题

lujinhui396

69 F= -.35908044E+03 E0= -.35908391E+03  d E =-.662334E-03  mag=    39.9793
trial-energy change:   -0.000662  1 .order   -0.000591   -0.000771   -0.000410
step:   0.8858(harm=  0.8858)  dis= 0.01514  next Energy=  -359.080598 (dE=-0.823E-03)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.359080483068E+03   -0.20475E-04   -0.23540E-01  4588   0.756E-01    0.189E-01
RMM:   2    -0.359080799125E+03   -0.31606E-03   -0.30718E-03  4491   0.132E-01    0.146E-01
RMM:   3    -0.359080830041E+03   -0.30917E-04   -0.22250E-04  2683   0.407E-02
  70 F= -.35908083E+03 E0= -.35908411E+03  d E =-.105527E-02  mag=    39.9753
curvature:  -8.32 expect dE=-0.239E-01 dE for cont linesearch -0.776E-05
trial: gam= 1.53799 g(F)=  0.287E-02 g(S)=  0.000E+00 ort = 0.180E-03 (trialstep = 0.207E+00)
search vector abs. value=  0.860E-01
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.359081266983E+03   -0.46786E-03   -0.11848E-01  4612   0.555E-01    0.146E-01
RMM:   2    -0.359081450893E+03   -0.18391E-03   -0.18601E-03  4441   0.104E-01    0.127E-01
RMM:   3    -0.359081440733E+03    0.10160E-04   -0.16456E-04  2650   0.353E-02
  71 F= -.35908144E+03 E0= -.35908486E+03  d E =-.610692E-03  mag=    39.9783
trial-energy change:   -0.000611  1 .order   -0.000524   -0.000650   -0.000398
step:   0.5327(harm=  0.5327)  dis= 0.01471  next Energy=  -359.081668 (dE=-0.838E-03)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.359081535863E+03   -0.84971E-04   -0.29542E-01  4572   0.810E-01    0.192E-01
RMM:   2    -0.359081891368E+03   -0.35550E-03   -0.37504E-03  4481   0.139E-01    0.148E-01
RMM:   3    -0.359081902138E+03   -0.10770E-04   -0.23106E-04  2661   0.383E-02
  72 F= -.35908190E+03 E0= -.35908526E+03  d E =-.107210E-02  mag=    39.9796
curvature:  -7.28 expect dE=-0.158E-01 dE for cont linesearch -0.508E-05
trial: gam= 0.85670 g(F)=  0.217E-02 g(S)=  0.000E+00 ort = 0.245E-03 (trialstep = 0.272E+00)
search vector abs. value=  0.657E-01
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.359082314221E+03   -0.42285E-03   -0.15806E-01  4652   0.567E-01    0.136E-01
RMM:   2    -0.359082493583E+03   -0.17936E-03   -0.19240E-03  4279   0.956E-02    0.112E-01
RMM:   3    -0.359082506508E+03   -0.12926E-04   -0.15480E-04  2624   0.336E-02
  73 F= -.35908251E+03 E0= -.35908575E+03  d E =-.604371E-03  mag=    39.9765
trial-energy change:   -0.000604  1 .order   -0.000501   -0.000647   -0.000355
step:   0.6026(harm=  0.6026)  dis= 0.01526  next Energy=  -359.082618 (dE=-0.716E-03)



这是我优化的结果，发现电子步很快就收敛，而离子步都70多步还没收敛，请问怎么解决呢？

nuclessenwin

70步还没收敛，恐怕初始构型不太好吧

nuclessenwin

INCAR贴上来看看

jiewei

离子难以优化，是不是自选改了？改了自旋优化就难。楼上说得对，上incar

lujinhui396

SYSTEM=Y-100-LITTLE-NMAG-C2H4
ENCUT=550
ISTART=0
ICHARG=2
VOSKOWN =1
PREC=Accurate
ISPIN=2
ISMEAR=1
SIGMA=0.1
ISIF=2
IBRION=2
POTIM=0.5
EDIFF=0.0001
EDIFFG=-0.02
NSW=100
NELM=100
ISYM=0
ALGO=FAST
LREAL=Auto
LCHARG = .F.
LMAXMIX = 4
AMIX=0.2
BMIX=0.0001
AMIX_MAG=0.8
BMIX_MAG=0.0001
LWAVE = .F.
LDAU= .TRUE.
LDAUTYPE= 2
LDAUL=2 -1 -1 -1
LDAUU=3.5 0.0 0.0 0.0
LDAUJ=0.0 0.0 0.0 0.0
LDAUPRINT=2
NPAR=4
NSIM=1
ADDGRID=.TRUE.

lujinhui396

nuclessenwin 发表于 2014-11-21 21:10
70步还没收敛，恐怕初始构型不太好吧

进行了初步优化后，在进一步优化的。

lujinhui396

jiewei 发表于 2014-11-21 22:42
离子难以优化，是不是自选改了？改了自旋优化就难。楼上说得对，上incar

因为有磁性的钴，所以打开自旋优化的

卡开发发

现在不知道结构及对应的磁性设置是不是有问题，先检查看看这个。
如果结构及磁性没问题，不知道KPOINTS和SIGMA是否充分，如果这两个不充分导致的结果就是力计算不准确，会造成不收敛。
如果还是不行，就得考虑调整POTIM，或者切换离子步的代数。

目前我也只能考虑到这么多。

lujinhui396

卡开发发 发表于 2014-11-22 20:41
现在不知道结构及对应的磁性设置是不是有问题，先检查看看这个。
如果结构及磁性没问题，不知道KPOINTS和S ...

谢谢。我的体系是四氧化三钴，体相的基态是反铁磁。但是我切面以后，发现对SLAB模型，我只打开自旋，即顺磁计算，能量更低。若保持体相的反铁磁，则能量高。
另外，我的是底层固定的slab，在做K点测试的时候，使用这个底层固定的slab模型吗？我的K点取的是SLAB模型没固定，测试的K点

卡开发发

lujinhui396 发表于 2014-11-24 08:52
谢谢。我的体系是四氧化三钴，体相的基态是反铁磁。但是我切面以后，发现对SLAB模型，我只打开自旋，即顺 ...

KPOINTS测试基于单点即可，这样锁定与否就没有影响了，但优化后最好进行重新验证。
我不知道你怎么保持表面反铁磁的，反铁磁需要开启自旋极化并指定每个原子初始磁矩。
还有ISMEAR=1？你的KPOINTS在2x2x1以下？

lujinhui396

卡开发发 发表于 2014-11-24 10:59
KPOINTS测试基于单点即可，这样锁定与否就没有影响了，但优化后最好进行重新验证。
我不知道你怎么保持 ...

我保持反铁磁，就是设置了原子的MAGMOM，并且正负相见，在上面贴出的INCAR里没有。
对于ISMEAR=1的疑问，是因为我的体系是半导体吗？我的K点是2 2 1 。对于这个问题，我询问了原作者,他说VASP manual is not right for all cases. MP is the best choice for the setting of the partial occupancies of a surface. We also did some tests with other methods such as Gaussian smearing. I reckon even the second order MP for your calculations.
另外，我用Ismear=0，sigma=0.05 计算，得出的能量：ISMEAR=0：-327.317EV ;ISMEAR=1：-327.514877 EV
（2）请问，一般，体相是反铁磁，切的面是否要保持这个磁序呢？有的文献里说，保持。但是我计算出来的能量，却是不保持的更低。

卡开发发

lujinhui396 发表于 2014-11-24 11:51
我保持反铁磁，就是设置了原子的MAGMOM，并且正负相见，在上面贴出的INCAR里没有。
对于ISMEAR=1的疑问 ...

（1）没，k=2x2x1去ismear=1倒是没啥太大争议，因为ismear>0适合体系较大k较少的情况。
（2）我做的一些表面也有磁序稍有变化的，但是好像从反铁磁变成铁磁的情况我没见过，不知道提高层数是否有帮助，如果层数不充分也可能导致一些问题。

这个问题看来涉及的因素比较多。。。

lujinhui396

卡开发发 发表于 2014-11-24 12:38
（1）没，k=2x2x1去ismear=1倒是没啥太大争议，因为ismear>0适合体系较大k较少的情况。
（2）我做的一些 ...

好的。我再找找问题，非常感谢

lujinhui396

卡开发发 发表于 2014-11-24 12:38
（1）没，k=2x2x1去ismear=1倒是没啥太大争议，因为ismear>0适合体系较大k较少的情况。
（2）我做的一些 ...

在一篇文献里找到。一般采用加U计算Co3O4,但是这篇文献指出，加U计算时，体相基态是铁磁的更稳定，而不是反铁磁。

卡开发发

lujinhui396 发表于 2014-11-24 17:18
在一篇文献里找到。一般采用加U计算Co3O4,但是这篇文献指出，加U计算时，体相基态是铁磁的更稳定，而不是 ...

这里我没明白你的意思，你计算的结果体相到底铁磁性稳定还是反铁磁稳定？表面呢？