| 查看: 796 | 回复: 5 | ||
mazuju028木虫 (正式写手)
|
[求助]
使用NUPDOWN的怪异现象 已有1人参与
|
|
最近对同一个体系，分别作了两次自洽计算， 第一次计算使用以下参数： ISPIN = 2 NUPDOWN = 2 第二次计算： ISPIN = 2 #NUPDOWN = 2 可以简单理解为第一次作了个三重态计算，第二次应该就是普通的单重态（基态）计算， 这个体系共有960个电子，所以理论上说第一种情形的积分态密度在费米能级处的值应该为向上481个，向下479个， 可是，从DOSCAR中得到这样的数据： -2.008 0.5450E+01 0.1612E+02 0.9596E+03 0.9577E+03 费米能级为-2.0019eV，可以看到向上的电子为959.6个，向下的电子为957.7个，两者相减为2，满足NUPDOWN的设置，可是电子数相加岂不是两倍总电子数了？ 对于第二次不加NUPDOWN, 计算结果就正常了，向上和向下的电子分别为480个。 DOSCAR显示如下： -1.226 0.2460E-08 0.2461E-08 0.4800E+03 0.4800E+03 费米能级为-1.2259eV 有没有人对第一种情况做个解释呢。 |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
mazuju028
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 30 (小学生)
- 贵宾: 0.141
- 金币: 2794.4
- 散金: 528
- 红花: 10
- 帖子: 942
- 在线: 419.4小时
- 虫号: 668835
- 注册: 2008-12-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2014-09-22 09:09:45
3楼2014-09-22 10:56:34
mazuju028
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 30 (小学生)
- 贵宾: 0.141
- 金币: 2794.4
- 散金: 528
- 红花: 10
- 帖子: 942
- 在线: 419.4小时
- 虫号: 668835
- 注册: 2008-12-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
请问你可以把你的INCAR贴出来，我看看吗？ 我的INCAR是这样的： SYSTEM = TiO ENCUT = 400 ISTART =0; ICHARG = 2 ISMEAR = 0 SIGMA = 0.2 #NSW = 200 #IBRION = 2; POTIM = 0.2 #ISIF = 2 #EDIFF = 1E-4; EDIFFG = -2E-2 EDIFF = 1E-6 PREC = Accurate LREAL = Auto IALGO = 38 NPAR = 4 #GGA+U LDAU= T LDAUTYPE = 2 LDAUL = 1 2 LDAUU = 5 4.2 LDAUJ = 0 0 LASPH = .T. LMAXMIX = 4 #spin ISPIN = 2 NUPDOWN = 2 |
4楼2014-09-22 14:21:48
mazuju028
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 30 (小学生)
- 贵宾: 0.141
- 金币: 2794.4
- 散金: 528
- 红花: 10
- 帖子: 942
- 在线: 419.4小时
- 虫号: 668835
- 注册: 2008-12-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2014-09-23 08:24:57
6楼2014-09-23 20:42:53