| 查看: 2358 | 回复: 9 | |||||
nkleof银虫 (正式写手)
|
[交流]
【原创】ICHARG=11和+U的计算 已有4人参与
|
|
vasp计算DOS有一种方法，就是以较低k网格达到收敛，再以很高的k网格通过ICHARG=11读取自洽CHGCAR进行DOS计算。 我昨天进行了这样的计算，但是在计算中进行了+U，如下是INCAR PREC = Accurate LREAL = A IALGO = 48 NPAR = 8 NSIM = 8 ISTART = 1 ICHARG = 11 NELMIN = 4 LORBIT = 12 EDIFF = 1E-5 ISMEAR = -5 EMIN = -20 EMAX = 20 NEDOS = 11000 ISPIN = 2 LDAU = .TRUE. LDAUTYPE = 2 LDAUL = -1 2 -1 LDAUU = 0 2 0 LDAUJ = 0 0 0 LDAUPRINT = 2 其中自洽收敛是以11 11 11网格，非自洽DOS是以21 21 21网格，另外也进行了21 21 21网格的自洽计算进行对比。 如下是三次的DOS 图，分别为k11网格自洽的DOS，在k11网格自洽基础上k21的DOS结果，k21网格自洽的DOS    从图上可以看到，在k11自洽基础上进行k21的ICHARG=11的DOS结果和其他两个有点不同，费米面以上的部分有所变化。 所以，我猜测可能是因为是INCAR里面的+U设置在后来的非自洽计算中重复起了作用，使得实际作用的U值成了4eV。即原来以U=2进行了自洽计算，得到的CHGCAR已经包含2eV 的Coulomb repulsion，后来进行非自洽计算时INCAR中的+U设置又起了作用，再加上了2eV的Coulomb repulsion，所以出来的DOS结果和原来U=2时的不同。 希望有高手来证实一下，或者解释一下k21的非自洽结果与原自洽结果的DOS不同的原因。 [ Last edited by nkleof on 2010-4-29 at 10:51 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 材料计算模拟实用技巧 | VASP |
» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
静态计算能收敛，同样的参数弛豫不能收敛，这是为什么？
已经有4人回复
- vasp5.2计算半导体Si光学性质的总结
已经有106人回复
- 态密度图跟文献里面的怎么都不一样，求救
已经有36人回复
- ISPIN=1时态密度有为负的区域，为嘛呢？
已经有6人回复
- 【求助成功】DOS计算ICHARG=11出现警告，请高手帮忙
已经有9人回复
- 【求助】VASP计算体系性能没有f能带电子
已经有24人回复
- 【求助】态密度困惑
已经有8人回复
- 【求助】关于拟合得到体弹性模量&自旋极化&LDA+U等问题，望大家指点指点
已经有4人回复
- 【求助】vasp静态计算能量值错误
已经有21人回复
- 【求助】LDA+U优化问题
已经有10人回复
- 【求助】VASP进行DFT+HF计算的问题？
已经有7人回复
- 【求助】计算单个Fe原子，无法收敛，请教
已经有19人回复
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
2楼2010-04-29 20:14:05
锐利的碎片
木虫 (正式写手)
star watcher
- 1ST强帖: 1
- 应助: 136 (高中生)
- 金币: 3637.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构
3楼2010-04-29 20:32:49
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
4楼2010-04-29 23:00:53
锐利的碎片
木虫 (正式写手)
star watcher
- 1ST强帖: 1
- 应助: 136 (高中生)
- 金币: 3637.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构
5楼2010-04-29 23:03:42
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
6楼2010-04-29 23:59:59
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
7楼2010-04-30 08:54:07
wujiandongde
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 4265.4
- 散金: 30
- 红花: 5
- 帖子: 182
- 在线: 704.8小时
- 虫号: 941604
- 注册: 2010-01-13
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2011-04-11 12:12:43
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
9楼2011-04-15 09:49:56
wujiandongde
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 4265.4
- 散金: 30
- 红花: 5
- 帖子: 182
- 在线: 704.8小时
- 虫号: 941604
- 注册: 2010-01-13
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2011-04-15 14:28:28
