| 查看: 1317 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
jmlv金虫 (著名写手)
|
[求助]
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. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
回复此楼» 猜你喜欢
292求调剂
已经有3人回复
-
340求调剂
已经有5人回复
-
生物学学硕求调剂
已经有10人回复
-
上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境）
已经有4人回复
-
材料学求调剂
已经有6人回复
-
303求调剂
已经有5人回复
-
一志愿武理085500机械专业总分300求调剂
已经有7人回复
-
考研调剂
已经有4人回复
-
281求调剂
已经有4人回复
-
0805 316求调剂
已经有6人回复
aaq2800
铁杆木虫 (知名作家)
- 应助: 51 (初中生)
- 金币: 27903.5
- 散金: 1522
- 红花: 29
- 沙发: 1
- 帖子: 9462
- 在线: 2189.1小时
- 虫号: 598401
- 注册: 2008-09-10
- 性别: GG
- 专业: 原子和分子物理
8楼2015-11-27 14:12:34
wangsoft
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 3085.6
- 红花: 8
- 帖子: 311
- 在线: 104.6小时
- 虫号: 1470573
- 注册: 2011-11-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2015-11-25 10:03:22
WanderingHeart
铁杆木虫 (著名写手)
- 应助: 90 (初中生)
- 金币: 13890.8
- 散金: 5
- 红花: 18
- 帖子: 1885
- 在线: 467.1小时
- 虫号: 257464
- 注册: 2006-06-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2015-11-25 12:45:20
jmlv
金虫 (著名写手)
- 应助: 50 (小学生)
- 金币: 4640.2
- 散金: 20
- 红花: 6
- 帖子: 1357
- 在线: 129小时
- 虫号: 83527
- 注册: 2005-07-30
- 性别: GG
- 专业: 催化化学
4楼2015-11-26 10:08:29
20