24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jmlv

金虫 (著名写手)

应助: 50 (小学生)
金币: 4640.2
散金: 20
红花: 6
帖子: 1357
在线: 129小时
虫号: 83527
注册: 2005-07-30
性别: GG
专业: 催化化学

[求助] Ni团簇的磁矩为什么总是波动？已有1人参与

我创建了一个气态游离(free standing)、79个原子的截角八面体Ni团簇，对其进行自旋极化计算是，出现咄咄怪事：

第一步，对从体相截出的Ni79团簇进行完全弛豫，收敛门槛设为EDIFF=1E-5，得到磁矩83.10，总能量-121.724 eV。

第二步，对第一步得到的优化构型进行完全弛豫，收敛门槛仍设为EDIFF=1E-5，按理说应该得到同样的磁矩和总能量，但事与愿违：磁矩82.64，总能量-137.479 eV。

总能量差别如此之大，令人无法信任也无法使用这种结果呀。有达人碰到过类似的问题吗？求解！谢谢。

附上我的INCAR参数：
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ENCUT  = 400.00 eV
PREC = Normal
LREAL  = Auto

EDIFF  = 1E-5
NSW = 1000
IBRION = 2
ISIF = 2

ISMEAR = 2
SIGMA  = 0.2

ISPIN  = 2             # spin polarized

# DFT+U settings
LDAU = .TRUE.    # switch on the +U corrections
LDAUTYPE = 2       # Dudarev's approach
LDAUL  = 2       # -1 no U; 1 for p electrons; 2 d and 3 f.
LDAUU  = 6.0    # U values
LDAUJ  = 0.0    # do not specify them for Dudarev's approach

# output controls
LCHARG = .FALSE.
LWAVE  = .FALSE.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

回复此楼