能带中心计算

d-带中心描述符广泛应用在多相催化研究领域，计算公式如下：

\[\epsilon_d = \frac{\int_{-\infty}^{\infty} \epsilon\rho _d(\epsilon)d\epsilon}{\int_{-\infty}^{\infty} \rho _d(\epsilon)d\epsilon}\]

使用 GVasp ，你可以使用下述命令很容易的计算能带中心：

gvasp band-center -j center.json

center.json 示例如下：

{
  "pos_file": "CONTCAR_dos",
  "dos_file": "DOSCAR_dos",
  "atoms": "C",
  "orbitals": "p",
  "xlim": [-8.76, 11.55]
}

center.json 文件中各参数的含义与 dos-plot 任务类似。

重要

If your LORBIT in INCAR is not 11 or 12, please add “LORBIT”: 10 in the center.json.