| 查看: 823 | 回复: 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 有没有人对第一种情况做个解释呢。 |
回复此楼» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有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
10