本帖最后由 丁越 于 2021-9-3 19:54 编辑
官网给的这个WO3的例子是立方晶系的，我用origin画了一个四方晶系的WO3能带。大概流程是这样的：
1.把晶体的cif（原胞画，超胞画出来的能带是折叠的）文件上传到seekpath网站，产生高对称点路径，并且下方有cp2k的输出文件模板。
2.cp2k_bs2csv.py转化.bs能带数据文件为.csv文件。举例来说我的四方晶系WO3高对称点路径是Γ—X—M—Γ—Z—R—A—Z|X—R|M—A，因此脚本处理后就会产生9个.csv文件（实际上只是对.bs文件中能量重新排列了一下）。第一个文件对应Γ—X，第二个文件对应X—M。。。。以此类推。.csv文件中，前三列是k点的分数坐标（seekpath网站上倒空间布里渊区信息那一块就有）。举例来说，对于Γ—X，cp2k输入文件中NPOINTS 10，这就相当于从Γ—X之间我们均匀撒10个k点，因此我们可以从第一个.csv文件看到均匀间隔的11个k点，之后的几列是每个k点对应的能带能量E(k)。 3.将这9个.csv文件导入origin同一个worksheet中（file,多个ASCII文件导入，导入设置中‘多文件导入模式’选择‘新建行’ ）。由于所有相邻k点都是均匀间隔的，我们删去前三列k点坐标中任意两列，保留的一列作为X轴，其数值我们按照0，0.05,0.1,0.15......等间距重新设置。这是因为我们要把Γ—X—M—Γ—Z—R—A—Z|X—R|M—A画在一条横坐标轴上，所以只需要考虑相邻高对称点相对间隔并且慢慢往上加间距就行。(gamma-X相距0.5，X-M相距0.5....) 4.能量扣除fermi能级能量后作为Y轴 5.高对称点在X轴位置自己数一下，在X轴上如何标注参考http://blog.wangruixing.cn/2019/07/11/band3/
我也是初学小白，对于能带图处理我的理解不知道有没有错误，有错误的地方还请多多指正。另外还有两个问题没有搞清楚，手动@卡开发发老师,还请老师指教： 1.对于掺杂晶胞，或者缺陷等，我们必须要构建超胞模型，那么这种超胞模型cp2k算能带如何构建高对称点路径，以及后续的能带反折叠处理？我知道vaspkit可以处理这种能带反折叠计算，但是对于cp2k我感觉好难啊，完全无从下手 2.seekpath网站上什么要用它产生的晶胞参数和坐标呢？我看它和原始晶胞参数是有些区别的，但是不知道它为什么这样做。
| 
发表于 Post on 2021-8-27 22:43:17
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author
