计算化学公社

标题: 求助表面模型SPIN DENSITY 出现问题 [打印本页]

作者
Author: 夙雪夜归 时间: 2021-7-20 09:32
标题: 求助表面模型SPIN DENSITY 出现问题

各位老师同学大家好，

学生在建立了一个1x1的磷酸硼表面模型之后对表面进行了饱和处理，得到了表面是羟基的一个slab模型（如图1）。我尝试的表面反应中有一步是羟基脱氢，我想得到脱去羟基后的slab模型的能量，于是直接去掉一个上表面的羟基氢进行结构优化，得到优化后的结构CONTCAR-1。但是结构优化后计算SPIN DENSITY的时候发现上下表面的羟基氧均出现了不同程度的自旋密度（如图2），在INCAR里添加了LORBIT=11后得到了vasp输出的magnetization （x）（如图3），其中显示的“magnetization density ”信息和SPINDENSITY图显示的信息一致；VASP版本5.3.5，KPOINTS设置的是6x6x6，Gamma

另外除了使用给出的结构文件进行了计算，还尝试添加z轴的偶极矫正（LDIPOL=.TRUE.，IDIPOL=3），设置RWIGS，增加真空层厚度，移动slab在盒子中z方向的位置，得到的magnetization（x）结果均显示上下表面都出现了“magnetization density ”

我尝试将下表面的羟基换成了H原子（如图4），则得到的结果中下表面没有出现“magnetization density ”（如图5）

(, 下载次数 Times of downloads: 0)

(, 下载次数 Times of downloads: 23) (, 下载次数 Times of downloads: 36) (, 下载次数 Times of downloads: 24) (, 下载次数 Times of downloads: 19) (, 下载次数 Times of downloads: 18)

目前的问题主要有三个：

1.是否存在INCAR设置不合理或者结构不合理的问题，从而导致了这样不合理的SPIN DENSITY 结果

2.下表面羟基直接换成H的做法是否可取

3.请问为什么会出现这样的结果呢

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