计算化学公社

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

作者
Author: 银河牛仔 时间: 2021-11-28 13:13
标题: 使用VASP进行结构优化时遇到的问题（金属铜吸附的NRR）
最近接触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

作者
Author: 银河牛仔 时间: 2021-11-28 13:19
出现的结果很奇怪，吸附位点的Cu原子被挤到很靠下的位置。有的计算还会出现Error EDDDAV: Call to ZHEGV failed. Returncode = 43 2 48的报错

作者
Author: 夙雪夜归 时间: 2021-11-28 13:39

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

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

作者
Author: 银河牛仔 时间: 2021-11-28 23:40

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

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

作者
Author: wsz 时间: 2021-11-29 03:10
可以试试设置磁矩初猜，磁矩不合理也可能导致这个错误。

作者
Author: 银河牛仔 时间: 2021-11-29 12:38

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

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

作者
Author: wsz 时间: 2021-11-30 11:55

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

建议开自旋设置磁矩

作者
Author: 银河牛仔 时间: 2021-12-1 15:52

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

好的，感谢老师

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