| 查看: 1901 | 回复: 8 | |||
[交流]
求助：vasp计算rutile TiO2的电子DOS时，为何单胞和supercell的结果不一样？ 已有2人参与
|
|
我一般都是用单胞来计算电子dos。这次因为想看一个特定的晶格distortion对rutile的影响，而这个distortion需要在一个2*2*2的supercell实现，所以我就尝试用一个supercell来算dos了。 优化的单胞展开成2*2*2的supercell（这个supercell没有任何distortion，我是想看一看单胞和8倍的supercell得到的东西是不是一样），开始计算，想着结果应该和单胞出来的没有什么区别啊（我觉得dos的值应该大约是8倍，除此之外分布应该没啥区别），但是奇怪的是supercell产生的dos很明显的少了一个高能级的conduction band peak，我不知道是咋回事。比较单胞和supercell的区别，只有kpoint的密度低了一倍，这应该是没有问题的，以下是我对这个超胞的计算文件 （1）先以4*4*4的 K point 优化结构，进而得到chgcar INCAR: SYSTEM = rutile Relax ISTART = 0 ICHARG = 2 ENCUT = 500 PREC = Accurate NELMIN = 4 # do a minimum of four electronic steps NELM = 40 # maximum of 40 electronic steps EDIFF = 1E-8 EDIFFG = -1E-6 # high accuracy ADDGRID = .TRUE. IBRION = -1 ISMEAR = 0 SIGMA = 0.05 LREAL = .False. LWAVE = .False. LPLANE = .TRUE. NPAR = 8 # twice the number of nodes ------------------------------------------------------------------------------ (2) 拷贝chgcar到新directory，增加KPOINT 到 11*11*11 INCAR SYSTEM = rutile Relax ISTART = 0 ICHARG = 11 ENCUT = 500 PREC = Accurate NELMIN = 4 # do a minimum of four electronic steps NELM = 40 # maximum of 40 electronic steps EDIFF = 1E-8 EDIFFG = -1E-6 # high accuracy ADDGRID = .TRUE. IBRION = -1 ISMEAR = 0 SIGMA = 0.05 LORBIT = 10 #EMIN = -15 #EMAX = 15 #NEDOS = 1000 LREAL = .False. LWAVE = .False. LPLANE = .TRUE. NPAR = 8 # twice the number of nodes ----------------------------------------------------------------------- 单胞的计算和上面一样，除了KPOINT 分别是8*8*8 和 21*21*21，我也改过ISMEAR 到-1，木有啥变化。两个系统内力都收敛的很小，是优化的很好的结构。二者的能量是7.999倍，非常好的8倍关系。附件中给出了元胞和supercell的比较，很明显，最上面的一个peak在超胞计算中不见了 本人菜鸟，实在是想不出来什么地方不对了，望高人指点，谢谢啦!  dos_rutile.png |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP计算TiO2表面能问题 求助
已经有8人回复
- VASP计算TiO2能带竟然没有带隙
已经有3人回复
- VESTA画电荷密度图的疑问
已经有20人回复
- TiO2（101）晶面的优化
已经有13人回复
- 【转帖】用Material Studio构建复杂模型及转换成poscar的2种方法
已经有56人回复
2楼2014-11-25 16:05:27
3楼2014-11-26 01:20:55
4楼2014-11-26 03:40:40
失败的路人
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2824.7
- 散金: 20
- 红花: 1
- 帖子: 362
- 在线: 292.9小时
- 虫号: 2506020
- 注册: 2013-06-13
- 性别: GG
- 专业: 理论和计算化学
5楼2014-11-26 08:19:41
|
遇见大神了  果然是专家，一语点破玄机，我发现1个单胞的时候vasp by default的nbands = 32, 但是8个超胞的情况下nbands = 216，分明小于8倍的关系，结果导致上层空带没有被考虑进去，所以少了一个Peak。。。我把nbands提高到256，也就是正好8倍，结果就一样了 这个问题我是不可能想到的，我以为限定了能量区间，该算的都能算了呢 这个问题vasp online manul说的好像不太对啊，上面说如果非spin polarized的计算，default的nbands = nelect/2+nion/2，我用的ti_pv所以是10个电子，o是6个，一共44个电子/cell，这个怎么也算不出nbands=32啊？ |
6楼2014-11-26 13:01:07
7楼2014-11-26 16:07:57
8楼2014-11-26 16:09:18
感谢感谢 9楼2014-11-27 04:09:15
100