| 查看: 1907 | 回复: 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人回复
感谢感谢 
9楼2014-11-27 04:09:15
2楼2014-11-25 16:05:27
3楼2014-11-26 01:20:55
4楼2014-11-26 03:40:40
