| 查看: 1605 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有8人回复
-
0703化学调剂
已经有3人回复
-
材料，纺织，生物（0856、0710），化学招生啦
已经有9人回复
-
一志愿中国海洋大学，生物学，301分，求调剂
已经有3人回复
-
280求调剂
已经有4人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有9人回复
-
工科材料085601 279求调剂
已经有6人回复
-
070300化学319求调剂
已经有5人回复
-
生物学071000 329分求调剂
已经有3人回复
-
0703化学调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP计算TiO2表面能问题 求助
已经有8人回复
- VASP计算TiO2能带竟然没有带隙
已经有3人回复
- VESTA画电荷密度图的疑问
已经有20人回复
- TiO2（101）晶面的优化
已经有13人回复
- 【转帖】用Material Studio构建复杂模型及转换成poscar的2种方法
已经有56人回复
7楼2014-11-26 16:07:57
2楼2014-11-25 16:05:27
3楼2014-11-26 01:20:55
4楼2014-11-26 03:40:40
50