| 查看: 2501 | 回复: 4 | |||||
cenwanglai荣誉版主 (知名作家)
老和山猥琐派九段
|
[交流]
【求助完毕】IBZKPT,KPOINTS与DOS计算
|
|
对于边长10A以上的金属氧化立方体Bulk体系，做DOS计算。一般遵循这样的路线 （1）ISMEAR =0 Gamma点结构优化 （2）增加K点，ISMEAR = 0 自洽计算电荷密度 （3）增加k点，ISMEAR = -5 非自洽计算DOSCAR 第（3）中，偷懒的K的设置是Auto,Monkhorst-Pack方式。但是常常出现问题，比如数字的奇偶性，是否Gamma Centered 等等。更常见的问题是由于K点不是很好的按照四面体的方式分布的，所以手册上推荐明确手动列出四面体方法的K点来，给出的方法是根据POSCAR和Auto generated KPOINTS生成IBZKPT的Fortran程序kpoints。kpoints不会编译，官方提供的是一个有makefile的几个模块组成的程序包。自己手写又没学好能带结构k空间理论。 有两个问题： (1)合理的计算DOS的KPOINTS应该是怎样的呢？包含所有高对称点的（Gamma点和，边，面，棱等周期性边界上的特殊点），按照四面体方式划分简约布里渊区的方式来取样？ (2)下面的这个方法行不行呢？ 在（3）中先简单的跑两步，就会生成IBZKPT文件，然后断掉任务，cp IBZKPT KPOINTS，然后重新运行。这样是不是要好一些呢？ [ Last edited by cenwanglai on 2010-9-1 at 17:10 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 | VASP |
» 猜你喜欢
范进中举一文的中心思想
已经有8人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有6人回复
-
我面上完蛋了
已经有13人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
人气不行了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp 计算表面DOS与文献对应不上 请求大家帮助.......
已经有20人回复
- 【求助】这样的dos合适吗？
已经有26人回复
- 【求助】VERY BAD NEWS! internal error in subroutine IBZKPT
已经有13人回复
- 【求助】关于ISMEAR
已经有12人回复
- 【求助】Hse06计算时KPOINTS出错
已经有30人回复
- 【重点讨论】计算能带时如何确定高对称点和k-points path?
已经有40人回复
- 【交流】从一个实际例子解读IBZKPT
已经有4人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
- 【第一性原理活动贴】VASP 官方手册学习活动
已经有56人回复
- 【求助】计算dos时很慢，而且WAVECAR很大
已经有21人回复
- 【转帖】关于从vasp4.5、4.6计算光学性质（组合帖）
已经有33人回复
★ ★
zzy870720z(金币+2):谢谢交流，呵呵 2010-08-23 19:52:12
cenwanglai(金币+2):谢谢回复，呵呵。 2010-08-24 15:43:11
zzy870720z(金币+2):谢谢交流，呵呵 2010-08-23 19:52:12
cenwanglai(金币+2):谢谢回复，呵呵。 2010-08-24 15:43:11
|
本帖内容被屏蔽 |
2楼2010-08-23 18:48:39
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
|
可能是我的问题没有表述清楚。 计算DOS有两种KPOINTS设置方式.一种是Auto M-P方式，另外一种是四面体方式，即列出k点坐标，然后给出四面体列表和权重。 手册上说以四面体方式给出的KPOINTS有两种方式生成，一种是根据POSCAR和Auto M-P generated KPOINTS 运行程序kpoint得到,另外是运行vasp，ISMEAR=-5的时候会生成IBZPKT，这个就是按照四面体方式划分第一布里渊区的。 问题是，ISMEAR=-5的时候，Auto M-P的KPOINTS既然就可以了，那为什么还要运行kpoint得到四面体方式的kpoints (explicit list of all coordinates and of the connecton tables for the tetrahedra)? |
3楼2010-08-24 15:42:11
★ ★ ★
zzy870720z(金币+1):谢谢指教，期待你更好的表现 2010-08-24 16:19:27
aylayl08(金币+2):谢谢提示 2010-08-24 16:26:54
cenwanglai(金币+1):谢谢回复，呵呵。我自觉手册读得还是比较认真的。看来还不行。 2010-08-24 17:59:34
zzy870720z(金币+1):谢谢指教，期待你更好的表现 2010-08-24 16:19:27
aylayl08(金币+2):谢谢提示 2010-08-24 16:26:54
cenwanglai(金币+1):谢谢回复，呵呵。我自觉手册读得还是比较认真的。看来还不行。 2010-08-24 17:59:34
|
本帖内容被屏蔽 |
4楼2010-08-24 16:08:52
江山8806
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1175.9
- 散金: 5
- 红花: 1
- 帖子: 177
- 在线: 175小时
- 虫号: 769154
- 注册: 2009-05-13
- 性别: MM
- 专业: 污染控制化学
5楼2010-08-24 17:17:38
20