| 查看: 969 | 回复: 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 有没有人对第一种情况做个解释呢。 |
回复此楼» 猜你喜欢
时间戳又变了8-15
已经有20人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
mazuju028
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 30 (小学生)
- 贵宾: 0.141
- 金币: 2770.4
- 散金: 528
- 红花: 10
- 帖子: 943
- 在线: 419.5小时
- 虫号: 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
- 金币: 2770.4
- 散金: 528
- 红花: 10
- 帖子: 943
- 在线: 419.5小时
- 虫号: 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
- 金币: 2770.4
- 散金: 528
- 红花: 10
- 帖子: 943
- 在线: 419.5小时
- 虫号: 668835
- 注册: 2008-12-05
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2014-09-23 08:24:57
6楼2014-09-23 20:42:53
