vasp多次优化后的能量均不一样

yiranfengbai

优化了不同层数的Co3O4模型，收敛后cp CONTCAR POSCAR做单点能，发现与优化的能量相差较大(优化：-82.5392，单点：-81.4278)，在第一次优化的基础上cp CONTCAR POSCAR再次优化，得到的能量也与第一次优化的能量相差较大。
修改了INCAR中的部分参数重新优化，发现每次优化的能量都与单点的能量相差1eV左右，这是为什么啊？
SYSTEM = 111-A-1
ISTART = 0
ICHARG = 2
PREC = Accurate
LREAL = Auto
GGA = PE
ALGO = Fast

ENCUT = 450
IBRION = 2
NSW = 400
ISIF = 2

ISMEAR = 0
SIGMA = 0.1
EDIFF = 1E-4
EDIFFG = -0.03
POTIM = 0.1
NELM = 300
NELMDL = 6

#ISPIN
ISPIN = 2
MAGMOM = -3 3 13*0

#LDA+U
LDAU = .TURE.
LDAUTYPE = 2
LDAUL = 2 -1
LDAUU = 3.3 0
LDAUJ = 0 0
LMAXMIX = 4

AMIX = 0.2
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001

LCHARG = .F.
LWAVE = .F.

jiangning198511

有没有对比自旋密度？

yiranfengbai

jiangning198511 发表于 2021-8-11 11:46
有没有对比自旋密度？

一直都没有考虑到自旋密度的问题，请教一下，我是要对比两次优化后的自旋密度吗？要作图对比吗

夙雪夜归

yiranfengbai 发表于 2021-8-11 16:54
一直都没有考虑到自旋密度的问题，请教一下，我是要对比两次优化后的自旋密度吗？要作图对比吗

在INCAR里设置LORBIT=11，OUTCAR里每个离子步结束会输出magnetization (x)，这一信息可以作为自旋密度的参考[img]
@jiangning198511老师您指的是这个吗

jiangning198511

夙雪夜归 发表于 2021-8-12 10:26
在INCAR里设置LORBIT=11，OUTCAR里每个离子步结束会输出magnetization (x)，这一信息可以作为自旋密度的 ...

是的 估计每次都收敛到不同的自旋态了

夙雪夜归

jiangning198511 发表于 2021-8-12 12:36
是的 估计每次都收敛到不同的自旋态了

老师您好，请问VASP有没有办法人为控制体系是开壳层单重态还是三重态呢，如果有的话该如何设置参数呢。体系是一个不饱和的表面B-O键和一个OOH自由基，类似双自由基

jiangning198511

夙雪夜归 发表于 2021-8-12 12:44
老师您好，请问VASP有没有办法人为控制体系是开壳层单重态还是三重态呢，如果有的话该如何设置参数呢。体 ...

试试UPDOWN关键词

夙雪夜归

jiangning198511 发表于 2021-8-12 13:07
试试UPDOWN关键词

谢谢，我去试试看

yiranfengbai

jiangning198511 发表于 2021-8-12 12:36
是的 估计每次都收敛到不同的自旋态了

老师，查看OUTCAR发现自旋态果然不对，请问该怎么办？
magnetization (x)

# of ion       s       p       d       tot
------------------------------------------
    1       -0.015  -0.042  -2.579  -2.636
    2       -0.014  -0.047  -2.570  -2.631
    3       -0.011  -0.042  -2.500  -2.553
    4        0.014   0.047   2.570   2.631
    5        0.015   0.042   2.579   2.636
    6        0.011   0.042   2.500   2.553
    7       -0.011  -0.056  -1.832  -1.900
    8        0.011   0.056   1.832   1.900
    9       -0.000   0.000   0.005   0.005
   10       -0.000  -0.000  -0.005  -0.005
   11        0.000  -0.000  -0.005  -0.005
   12       -0.005  -0.002   0.120   0.113
   13        0.005   0.002  -0.120  -0.113
   14        0.000   0.000   0.005   0.005
   15       -0.000  -0.000  -0.000  -0.000
   16       -0.012  -0.063  -1.898  -1.973
   17        0.012   0.063   1.898   1.973
   18       -0.011  -0.055  -1.818  -1.884
   19        0.000   0.000   0.005   0.005
   20        0.011   0.055   1.818   1.884
   21       -0.000  -0.000  -0.005  -0.005
   22       -0.004  -0.024   0.000  -0.029
   23       -0.004  -0.030   0.000  -0.035

卡开发发

yiranfengbai 发表于 2021-8-13 14:49
老师，查看OUTCAR发现自旋态果然不对，请问该怎么办？
magnetization (x)

用MAGMOM初猜一下，以及多cp CONCAR POSCAR多优化几遍。

yiranfengbai

卡开发发 发表于 2021-8-13 15:13
用MAGMOM初猜一下，以及多cp CONCAR POSCAR多优化几遍。

好的老师，我试试先优化几遍

卡开发发

yiranfengbai 发表于 2021-8-13 15:18
好的老师，我试试先优化几遍

MAGMOM一定得设置下，这个直接决定了初始猜测的电子分布情况。