使用VASP进行结构优化时遇到的问题（金属铜吸附的NRR）

银河牛仔 · 发表于 Post on 2021-11-28 13:13:56

最近接触VASP计算，在以Cu金属slab模型为基底计算NRR过程时，结构优化的结果总是不尽人意，以下是相关参数，向各位前辈请教

k点设置3x3x1,slab真空层15A，固定下面两层，放开上面两层和吸附的原子，进行结构优化，
ISTART=0
ICHARG=2
ISPIN=2
PREC=N
LREAL=.F.

ALGO=V
NELM=80
EDIFF=1E-5
ENCUT=400

LDIPOL=.TRUE.
IDIPOL=3

IBRION=2
NSW=200
ISIF=2
EDIFFG=-0.05

ISMEAR=1
SIGMA=0.2

银河牛仔 · 发表于 Post on 2021-11-28 13:19:28

出现的结果很奇怪，吸附位点的Cu原子被挤到很靠下的位置。有的计算还会出现Error EDDDAV: Call to ZHEGV failed. Returncode = 43 2 48的报错

夙雪夜归 · 发表于 Post on 2021-11-28 13:39:11

银河牛仔发表于 2021-11-28 13:19
出现的结果很奇怪，吸附位点的Cu原子被挤到很靠下的位置。有的计算还会出现Error EDDDAV: Call to ZHEGV fa ...

algo=V换成normal或者fast试试看还会不会报错

银河牛仔 · 发表于 Post on 2021-11-28 23:40:09

夙雪夜归发表于 2021-11-28 13:39
algo=V换成normal或者fast试试看还会不会报错

之前用的Normal就是这样，后来用very fast还是不行，吸附位点的Cu原子会被挤下去很多。会不会是自旋极化要关掉？

wsz · 发表于 Post on 2021-11-29 03:10:23

可以试试设置磁矩初猜，磁矩不合理也可能导致这个错误。

银河牛仔 · 发表于 Post on 2021-11-29 12:38:32

wsz 发表于 2021-11-29 03:10
可以试试设置磁矩初猜，磁矩不合理也可能导致这个错误。

纯铜的体系要不要把自旋关掉好一点

wsz · 发表于 Post on 2021-11-30 11:55:32

银河牛仔发表于 2021-11-29 12:38
纯铜的体系要不要把自旋关掉好一点

建议开自旋设置磁矩

银河牛仔 · 发表于 Post on 2021-12-1 15:52:08

wsz 发表于 2021-11-30 11:55
建议开自旋设置磁矩

好的，感谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 使用VASP进行结构优化时遇到的问题（金属铜吸附的NRR）

浏览过的版块