| 查看: 1173 | 回复: 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. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
回复此楼» 猜你喜欢
三无产品还有机会吗
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有7人回复
-
压汞仪和BET测气凝胶孔隙率
已经有4人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
-
投稿精细化工
已经有6人回复
jmlv
金虫 (著名写手)
- 应助: 50 (小学生)
- 金币: 4640.2
- 散金: 20
- 红花: 6
- 帖子: 1357
- 在线: 129小时
- 虫号: 83527
- 注册: 2005-07-30
- 性别: GG
- 专业: 催化化学
4楼2015-11-26 10:08:29
wangsoft
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 3085.6
- 红花: 8
- 帖子: 311
- 在线: 104.6小时
- 虫号: 1470573
- 注册: 2011-11-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2015-11-25 10:03:22
WanderingHeart
铁杆木虫 (著名写手)
- 应助: 90 (初中生)
- 金币: 13889.8
- 散金: 5
- 红花: 18
- 帖子: 1885
- 在线: 467.1小时
- 虫号: 257464
- 注册: 2006-06-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2015-11-25 12:45:20
WanderingHeart
铁杆木虫 (著名写手)
- 应助: 90 (初中生)
- 金币: 13889.8
- 散金: 5
- 红花: 18
- 帖子: 1885
- 在线: 467.1小时
- 虫号: 257464
- 注册: 2006-06-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2015-11-26 11:06:26