求助做Re掺杂NiFe-LDH体系结构优化时体系原子跑散

wsh1224

计算体系为4×4×1的NiFe-LDH，在对NiFe-LDH进行结构优化收敛后，在CONTCAR基础上掺杂了一个Re原子得到Re掺杂NiFe-LDH体系，并对该体系进行结构优化，但优化时体系总是跑散，猜测是初始体系过于糟糕。
想请教各位老师两个问题：
1.做元素掺杂体系时，应在LDH优化前掺杂，还是在LDH优化后掺杂？
2.在优化LDH后进行Re原子掺杂，但掺杂后优化时结构跑散，这是什么原因，需要怎样调整Re掺杂后的结构以使优化正常进行？

(优化采用很粗糙的参数仍然跑散)
INCAR:
SYSTEM = OPT-1

ISMEAR = 0
SIGMA = 0.1

#ISPIN = 2

IBRION = 2
POTIM = 0.5
NSW = 200

ENCUT = 400
EDIFF = 0.001
EDIFFG = -0.2

IVDW = 11
KPINTS采用gamma111
POSCAR见附件
POSCAR (3.89 KB, 下载次数 Times of downloads: 11)

2023-11-28 13:14 上传 Uploaded

点击下载
Click to download

Exodus

这个体系的Ni和Fe应该是高自旋态，既没有考虑磁矩的变化，k点取样也明显不足

wsh1224

Exodus 发表于 2023-11-28 14:11
这个体系的Ni和Fe应该是高自旋态，既没有考虑磁矩的变化，k点取样也明显不足

老师您好，我最开始优化的时候是加了自旋极化，k点取331，但是Re掺杂后的结构在优化时跑散，在网上查找解决办法时发现有人提到：
“一般来说，针对理论催化的计算，初始结构都是不太合理的。因此一开始采用很粗糙的优化（EDIFF=0.001，EDIFFG=-0.2），很低的K点密度(Gamma)，不考虑自旋就可以了，这样NSW<60的设置就比较好。其它参数可以默认。

经过第一轮优化，就可以进入下一步细致的优化了。就我的经验，EDIFF=1E-4,EDIFFG=-0.05，不考虑自旋，IBRION=2，其它默认，NSW=100;跑完后可以设置IBRION = 1，减小OPTIM（默认为0.5，可以设置0.2）继续优化。
https://blog.sciencenet.cn/blog-567091-895440.html”
我认为跑散是由于初始结构过于不合理所以按照该办法进行分步优化，但在进行第一步粗糙的优化时体系仍然跑散，所以有点摸不着头脑。

Exodus

NiFe-LDH这个结构算的人已经很多了，结构已经是基本确定了的，不会偏差很大。其次用很粗糙的精度只会让原子受力计算不准确，反而会让结构偏离势能面极小值点，我觉得你更应该关注不同化学环境的原子所处的自旋态，否则即使算出来也没什么意义

wsh1224

Exodus 发表于 2023-11-28 14:33
NiFe-LDH这个结构算的人已经很多了，结构已经是基本确定了的，不会偏差很大。其次用很粗糙的精度只会让原子 ...

好的，感谢老师指点，我学习一下

sobereva

帖子标题不得乱用叹号，置顶的新社员必读贴以及http://bbs.keinsci.com/thread-9348-1-1.html都明确说了。这次给你改了，下次乱用叹号直接删帖扣分处理。