BoltzTraP2 石墨烯二维等能图
用实际石墨烯均匀网格能带运行二维入口,说明化学势、能量窗、BZ边界和二维单位。
阅读本篇目录
fermisurface_2d适用于垂直于z方向的准二维体系:输入k点的z分量为0,晶胞的平面内基矢与z方向分开。本例使用已有真实石墨烯15×15×1电子结构。
下载完整BoltzTraP2实例包。包内有已有集群Si与石墨烯DFT数据、明确标注来源的上游格式例子、32个操作条件、原始输出、独立提取与绘图脚本,不含POTCAR。
1. 一次完成插值和两张图#
在新目录:
bash fermi2d.sh
ls work/graphene-fermi-electrons
脚本先生成石墨烯m5插值,再分别画µ−E_F(DFT)=−0.30和+0.30 eV的情况。能量参照是原DFT费米能,不自动等同于任意模型的Dirac点。
2. 底层入口#
在包含自己 bands.bt2的目录中:
btp2 fermisurface_2d bands.bt2 0.3
有图形显示环境时是带滑块的窗口;集群无显示环境时,请用配套 fermi2d.sh或 run.sh入口保存PNG。配套渲染器保持原始能带和等能数据,仅适配现代Matplotlib轮廓裁剪接口,并让初始滑块标注与请求的化学势一致。
3. 读图#


黑色六边形是第一布里渊区。当前原生绘图带有有限能量宽度,默认在µ附近±0.05 eV形成区域;不要把有宽度的区域当成无限精细的一条费米线。
边界上的多片等能区域可能通过周期边界相互等价,不能仅数图上的六个角就断言有六个独立谷。
4. 改化学势#
需要自己尝试更多值时,重新解压到新目录,使用原生交互窗口或在独立输出目录调用配套渲染入口。fermisurface_2d默认化学势单位eV;不要把插值能量窗参数使用的Ha直接填到这里。
5. 二维电导另有厚度约定#
本页展示等能区域,没有把含真空的三维超胞体积当作二维材料的固定实体体积。若进一步计算输运,需要说明采用有效厚度还是片电导;任意改变真空层后仍把三维体电导直接互比,会引入人为体积依赖。
默认环境未启用依赖VTK的三维 fermisurface入口;这与已经运行的二维入口是两回事。