24小时热门版块排行榜

返回列表

QuenZ

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 811.5
散金: 5
红花: 1
帖子: 204
在线: 78.7小时
虫号: 1181414
注册: 2011-01-01
性别: MM
专业: 理论和计算化学

[交流] 固定磁矩（NUPDOWN）优化试验

最近跟一个老外老头子询问吸附物质与基底之间的电子交换，被问及在计算时有没有设置NUPDOWN（VASP）参数，他说如果alpha和beta电子数没设好，就不好说计算出来的电子交换数对不对了，让我拿O2去测试一下。
我就把NUPDOWN=-1,0,1,2都算了下，-1是默认值，自动寻找分布电子，0就是alpha=beta，1是alpha-beta=1，2是alpha-beta=2。按道理2应该是最符合O2三线态的情况了。但从能量和态密度上看2还是有点问题，不知道问题出在哪，所以希望共同讨论下。
INCAR：PREC=Accurate
               ISPIN=2
               ENCUT=520
               EDIFFG=-0.02
               ISMEAR=0
               SIGMA=0.05
计算结果：
NUPDOWN = -1,       0,    1,       2
Energy/eV  = -9.859, -8.751, -9.206, -9.837
Mag/uB    = 2,       0,       2,       2
DOS图和能带上电子占据情况需见附件。

0和1从能量上就可以否定了，也列举了NUPDOWN=0时dos及能级占据情况，2s轨道电子被拿到2px和2py的2π*轨道上去了，占据数为0.5，不予考虑了。
主要对比-1和2的结果，发现-1从能量上是要比2的更低的。其它的电子占据及磁矩二者都一致，符合三线态电子分布。唯一不一样的地方是，NUPDOWN=2时，电子排到了高出费米能级5eV的轨道，认为2π*轨道是在费米能以上了。而NUPDOWN=-1时，电子最高分布没有到费米能。我就有些纠结了：2π*轨道就是HOMO轨道吧，那费米能就应该在HOMO附近吧，为什么NUPDOWN=2时HOMO跑到那么高能量的地方去了？那么NUPDOWN=-1时，轨道的能量分布是更可靠的了，所以结论岂不就是用NUPDOWN=-1更靠谱？

回复此楼