Fe单原子掺杂石墨烯结构优化问题求助

小张要努力

大家好，采用VASP优化Fe单原子掺杂石墨烯结构时发现，优化完Fe原子的z坐标与碳原子z坐标一致，并不是向文献中，是向外凸起的情况。而且计算出的磁矩也与文献不一致（文献为0，我计算出1点多）。请教一下大家
（1）利用MS建立掺杂模型：先import graphene --Make P1-删除一层——超胞为4 4 1——将晶格z轴改为20——利用Vesta导出为vasp格式。
（2）INCAR文件为：
ISTART=0   
ICHARG=2   
PREC=N

ISIF=2（这个2和3都试了，结果没啥变化）
IBRION=2
NSW=300

EDIFFG=-0.02  

ALGO=N
ENCUT=500   
ISPIN=2  
MAGMON=31*0 1*2 （这句存在和不存在情况也都试了，磁矩和结构依旧不正确）
EDIFF=1E-5
NELM=200


ISMEAR=0
SIGMA=0.05
LREAL=Auto
IVDW = 11

K点文件：
Automic generation
0
Monhkorst-Pack
7 7 1
0.0 0.0 0.0

结构文件：

优化后的结构文件：

（3）文献为：The adsorption characteristics of As2O3,Pb0, PbO and PbCl2 on single atom iron adsorbent with graphene-based substrates
文献结果为：

白皛白

可能是在建模的时候就严格建在一个平面上，这种情况下很难直接通过结构优化形成弯曲或者凸起的结构，可以在POSCAR中把铁原子的位置往垂直于石墨烯平面的方向移动0.5~1埃再优化一下试试。

小张要努力

白皛白 发表于 2024-11-7 20:43
可能是在建模的时候就严格建在一个平面上，这种情况下很难直接通过结构优化形成弯曲或者凸起的结构，可以在 ...

好的，十分感谢您的建议，我试试

后到咸阳斯蒂庞

小张要努力 发表于 2024-11-8 09:37
好的，十分感谢您的建议，我试试

因为结构优化只能保证构型的能量收敛到势能面上的极小值，而不是全局最小值，所以初始构型会优化结果产生显著影响，不妨可以按照楼上白皛白老师说的调整一下Fe原子在空间上的距离，例如分别取0.5、1、1.5的距离进行优化，最后选择能量最低的构型用于后续计算。

小张要努力

后到咸阳斯蒂庞 发表于 2024-11-12 10:40
因为结构优化只能保证构型的能量收敛到势能面上的极小值，而不是全局最小值，所以初始构型会优化结果产生 ...

十分感谢您的回复后面根据文献调整的初始优化结构，确实可以了！