| 查看: 1319 | 回复: 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. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
回复此楼» 猜你喜欢
312求调剂
已经有12人回复
-
生物学学硕，一志愿湖南大学，初试成绩338
已经有3人回复
-
学硕274求调剂
已经有3人回复
-
279求调剂
已经有6人回复
-
085701环境工程，267求调剂
已经有3人回复
-
中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培
已经有5人回复
-
274求调剂
已经有20人回复
-
265求调剂
已经有16人回复
-
材料与化工304求B区调剂
已经有6人回复
-
321求调剂
已经有5人回复
jmlv
金虫 (著名写手)
- 应助: 50 (小学生)
- 金币: 4640.2
- 散金: 20
- 红花: 6
- 帖子: 1357
- 在线: 129小时
- 虫号: 83527
- 注册: 2005-07-30
- 性别: GG
- 专业: 催化化学
6楼2015-11-26 16:06:18
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
