计算化学公社

标题: vasp多次优化后的能量均不一样 [打印本页]

作者
Author: yiranfengbai 时间: 2021-8-9 16:41
标题: vasp多次优化后的能量均不一样
优化了不同层数的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.

作者
Author: jiangning198511 时间: 2021-8-11 11:46
有没有对比自旋密度？

作者
Author: yiranfengbai 时间: 2021-8-11 16:54

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

一直都没有考虑到自旋密度的问题

，请教一下，我是要对比两次优化后的自旋密度吗？要作图对比吗

作者
Author: 夙雪夜归 时间: 2021-8-12 10:26
本帖最后由夙雪夜归于 2021-8-12 10:27 编辑

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

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

作者
Author: jiangning198511 时间: 2021-8-12 12:36

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

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

作者
Author: 夙雪夜归 时间: 2021-8-12 12:44

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

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

作者
Author: jiangning198511 时间: 2021-8-12 13:07

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

试试UPDOWN关键词

作者
Author: 夙雪夜归 时间: 2021-8-12 13:38

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

谢谢，我去试试看

作者
Author: yiranfengbai 时间: 2021-8-13 14:49

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

作者
Author: 卡开发发 时间: 2021-8-13 15:13

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

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

作者
Author: yiranfengbai 时间: 2021-8-13 15:18

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

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

作者
Author: 卡开发发 时间: 2021-8-13 15:21

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)