计算材料学：分子能带计算与各原子的轨道分析

您可以点此 下载链接 下载实验报告的 PDF 扫描版本，或者继续阅读本页面，其由 AI OCR 技术生成但可能存在识别误差的电子文档，建议以扫描版为准，感谢您的查阅。

一、实验预习

1. 实验目的

（1）以 $\text{CaF}_2$（萤石结构）和 $\text{GaAs}$（闪锌矿结构）为计算对象，在 Materials Studio（CASTEP模块）中完成能带结构（Band Structure）和投影态密度（PDOS）的第一性原理计算，学习电子结构计算的完整流程；

（2）通过对比两种材料的能带结构，识别直接带隙与间接带隙，提取带隙宽度，结合 PDOS 分析各原子轨道对价带和导带的贡献，理解电子态的轨道起源；

（3）通过电子密度差分图（Electron Density Difference）的三维可视化，直观辨别 $\text{CaF}_2$ 的离子键特征与 $\text{GaAs}$ 的共价键特征，理解化学键类型与电子结构的对应关系；

（4）结合附图中的能带拓扑分类（直接/间接/反转带隙三种情形），讨论带隙类型对光电材料应用的影响。

2. 实验原理

2.1 能带结构理论

在周期性势场中，电子的波函数满足 Bloch 定理：

$$ \psi_{n\mathbf{k}}(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}} u_{n\mathbf{k}}(\mathbf{r}) $$

其中 $n$ 为能带指标，$\mathbf{k}$ 为波矢，$u_{n\mathbf{k}}(\mathbf{r})$ 为与晶格同周期的周期函数。对应的本征能量 $\varepsilon_n(\mathbf{k})$ 构成能带，在布里渊区高对称路径上的分布即为能带结构。

Kohn-Sham 方程在倒空间的矩阵形式为：

$$ \sum_{\mathbf{G}'}\left[\frac{\hbar^2}{2m}|\mathbf{k}+\mathbf{G}|^2\delta_{\mathbf{G}\mathbf{G}'} + V_{\text{eff}}(\mathbf{G}-\mathbf{G}')\right]c_{n\mathbf{k}}(\mathbf{G}') = \varepsilon_{n\mathbf{k}}\,c_{n\mathbf{k}}(\mathbf{G}) $$

能带结构计算分为两步：首先做自洽场（SCF）计算得到收敛的电荷密度，再在高对称路径 $\mathbf{k}$ 点上做非自洽计算（Non-SCF）求解本征值。

2.2 带隙的定义与分类

带隙定义为价带顶（VBM）与导带底（CBM）的能量差：

$$ E_g = E_{\text{CBM}} - E_{\text{VBM}} $$

按 VBM 与 CBM 的 $\mathbf{k}$ 空间位置分类（对应附图三种情形）：

情形(a)——直接带隙（Direct Gap）：

$$ \mathbf{k}_{\text{VBM}} = \mathbf{k}_{\text{CBM}} = \mathbf{K}_+ $$

VBM 与 CBM 位于同一 $\mathbf{k}$ 点，电子跃迁无需声子参与，跃迁概率高，对应高效光电材料（如 GaAs）。

情形(b)——间接带隙（Indirect Gap）：

$$ \mathbf{k}_{\text{VBM}} \neq \mathbf{k}_{\text{CBM}}, \quad E_g^{K_+} \neq E_g^{K_-} $$

VBM 与 CBM 位于不同 $\mathbf{k}$ 点，光学跃迁需要声子辅助动量守恒，跃迁概率低（如 Si、Ge）。

情形(c)——带隙反转（Inverted/Topological Gap）：

$$ E_g^{K_+} > 0,\quad E_g^{K_-} < 0 \Rightarrow \text{能带交叉} $$

两个 $\mathbf{k}$ 点的带隙符号相反，出现能带反转，是拓扑绝缘体的标志性特征，表面态受时间反演对称性保护。

2.3 态密度与投影态密度（PDOS）

总态密度（Total DOS）定义为：

$$ g(E) = \frac{1}{N}\sum_{n,\mathbf{k}} \delta(E - \varepsilon_{n\mathbf{k}}) $$

实际计算中用 Gaussian 展宽代替 $\delta$ 函数：

$$ g(E) = \frac{1}{N}\sum_{n,\mathbf{k}} \frac{1}{\sigma\sqrt{2\pi}}\exp\left(-\frac{(E-\varepsilon_{n\mathbf{k}})^2}{2\sigma^2}\right) $$

投影态密度（PDOS） 进一步将各 $\mathbf{k}$ 点本征态投影到原子轨道基上：

$$ g_{\alpha lm}(E) = \sum_{n,\mathbf{k}} |\langle\phi_{\alpha lm}|\psi_{n\mathbf{k}}\rangle|^2\, \delta(E-\varepsilon_{n\mathbf{k}}) $$

其中 $\phi_{\alpha lm}$ 为第 $\alpha$ 个原子的 $|l,m\rangle$ 轨道。PDOS 直接揭示各原子、各轨道对特定能量范围电子态的贡献，是判断化学键性质的核心工具。

2.4 电子密度差分图

电子密度差分定义为：

$$ \Delta\rho(\mathbf{r}) = \rho_{\text{crystal}}(\mathbf{r}) - \sum_{\alpha}\rho_{\text{atom},\alpha}(\mathbf{r}) $$

即晶体电子密度减去孤立原子电子密度之和。$\Delta\rho > 0$（电荷积累）区域颜色偏蓝，$\Delta\rho < 0$（电荷亏损）区域颜色偏红，通过等值面（Isosurface）的空间分布，直观判断键合类型：

$$ \begin{cases} \Delta\rho\ \text{积累在键区（两原子间）} & \Rightarrow \text{共价键（方向性）} \\ \Delta\rho\ \text{积累在阴离子球形区域} & \Rightarrow \text{离子键（球形对称性）} \end{cases} $$

2.5 离子键与共价键的电子结构特征对比

二、实验过程

1. 实验过程记录

（1）晶体结构建立与弛豫

本实验计算两种材料：

$\text{CaF}_2$（萤石型，空间群 $Fm\bar{3}m$）：
原胞含3个原子（1个Ca位于 $4a$ $(0,0,0)$，2个F位于 $8c$ $(0.25, 0.25, 0.25)$），初始晶格常数 $a = 5.463$ Å。从实验截图可见，CASTEP Electronic Options 中 Selected potentials 设置为：F采用 F\_00PBE.usp（截断能330 eV），Ca采用 Ca\_00PBE.usp（截断能280 eV），整体截断能设为 400 eV（取最高要求），K点网格 $4\times4\times4$，Ultrasoft赝势，Koelling-Harmon相对论修正，Representation选择 Reciprocal space。

$\text{GaAs}$（闪锌矿型，空间群 $F\bar{4}3m$）：
原胞含2个原子（Ga位于 $4a$ $(0,0,0)$，As位于 $4c$ $(0.25,0.25,0.25)$），初始晶格常数 $a = 5.653$ Å。CASTEP Electronic Options 中 Selected potentials：Ga采用 Ga\_00PBE.usp（截断能255 eV），As采用 As\_00PBE.usp（截断能300 eV），整体截断能400 eV，K点 $4\times4\times4$。

两种材料均先进行几何优化（Relaxed CaF2 CASTEP GeomOpt、Relaxed GaAs CASTEP GeomOpt），再以弛豫后结构开展能带计算。

（2）能带结构计算设置

CASTEP Calculation 中 Properties 选项卡设置如下（截图所示）：

勾选：Band structure（能带），Density of states（态密度，勾选 Calculate PDOS），Energy range: 10.0 eV，k-point set: Medium（$6\times6\times6$，比SCF步更密以获得平滑能带），Display DOS: Total，Function: DOS，Scissors: 0.0 eV。

不勾选：Core level spectroscopy，Electron density difference（单独分析步骤），Electron localization function。

（3）电子密度差分图计算

能带计算完成后，在 CASTEP Analysis 工具中选择 Electron density，Results file 选择对应 .castep 文件，Density field 选择 Total density，勾选 View isosurface on import，导入后在 Display Style 中调整 Isovalue：

• $\text{CaF}_2$：Isovalue = 1.0768，Type: accessible，Color by mapped field（CASTEP density difference），显示 $\text{Ca}^{2+}$ 周围电荷亏损（红色）与 $\text{F}^-$ 周围电荷积累（等值面封闭球形）
• $\text{GaAs}$：Isovalue = 0.019662，蓝色区域为 Ga-As 键区电荷积累，红色为亏损区，呈现出明显方向性

（4）PDOS分析

在 CASTEP Analysis → Density of states 中，Atom Selection 选择各元素（Ca、F 分别选择；Ga、As 分别选择），生成各原子的投影分波态密度曲线，由 Python 脚本整合输出全范围 PDOS 图（图3、图4）。

三、分析讨论

1. 实验结果与分析

（1）$\text{CaF}_2$ 能带结构与PDOS分析

能带结构：

CASTEP 计算得到的 $\text{CaF}_2$ 能带图沿高对称路径 W-L-G-X-W-K 展示。主要特征如下：

价带顶（VBM）位于 $\Gamma$ 点（G点），导带底（CBM）位于 W 点（或 X 点附近），因此 $\text{CaF}_2$ 为间接带隙绝缘体，计算带隙为：

$$ \boxed{E_g(\text{CaF}_2) = 7.101\ \text{eV}} $$

实验值约为 $11.8$ eV，GGA-PBE 计算值严重低估（约低估40%），这是 GGA 对带隙系统性低估的典型体现（自相互作用误差，导带混入了过多 Ca 的 $3d$ 空轨道）。

能带整体特征：价带区域（$-20$ 至 $0$ eV）带较平坦，色散小，说明 F $2p$ 态局域性强；导带色散较大；价带底约 $-35$ eV 处有一条极平的孤立带，对应 F $2s$ 轨道的强局域态。

PDOS分析（CaF2 Projected Density of States）：

离子键特征分析：

PDOS 图中 Ca 与 F 的峰几乎不重叠，Ca 的峰集中在导带（未占据），F 的峰集中在价带（已占据），说明 Ca 的价电子（$4s^2$）已完全转移至 F，形成 Ca$^{2+}$ 和 F$^-$。这是离子键最典型的电子结构特征。

电子密度差分图（Isovalue = 1.0768）印证了这一判断：等值面以 Ca 原子为中心呈圆形亏损，以 F 原子为中心呈球形积累，无方向性聚集，符合离子键的球形对称电荷分布。

$$ \text{Ca}\ (4s^2) + \text{F}\ (2p^5) \rightarrow \text{Ca}^{2+} + 2\text{F}^- \quad \text{（两个F得到共2个电子）} $$

（2）GaAs 能带结构与PDOS分析

能带结构：

GaAs 能带图同样沿 W-L-G-X-W-K 路径展示。关键特征：

VBM 与 CBM 均位于 $\Gamma$ 点，因此 GaAs 为直接带隙半导体，计算带隙为：

$$ \boxed{E_g(\text{GaAs}) = 0.139\ \text{eV}} $$

实验值为 1.42 eV，GGA-PBE 严重低估（约低估90%），GaAs 的带隙低估问题比 $\text{CaF}_2$ 更严重，原因在于 Ga 的 $3d$ 半芯态与 As $4s$ 轨道的 $p$-$d$ 排斥效应在 GGA 中描述不足，若采用 GGA+U 或 HSE06 杂化泛函可将带隙修正至 $\sim1.1$ eV，更接近实验值。

能带色散特征：价带顶附近有明显抛物线色散，表明空穴有效质量较小；导带底 $\Gamma$ 点的色散也较大（电子有效质量小），两者共同决定了 GaAs 优异的载流子迁移率。

PDOS分析（GaAs Projected Density of States）：

共价键特征分析：

与 $\text{CaF}_2$ 形成鲜明对比：GaAs 的 Ga 与 As 的 PDOS 峰在价带区间（$-10 \sim 0$ eV）大量能量重叠，说明两种原子的轨道发生了显著的杂化混合，是共价键的直接体现。

GaAs 为 $sp^3$ 杂化，Ga（$4s^24p^1$）和 As（$4s^24p^3$）各贡献 4 个价电子（共8个），形成 4 个 $sp^3$ 杂化轨道，每个 Ga-As 键由两个原子各贡献 1 个电子共享。

电子密度差分图（Isovalue = 0.019662）直观显示：电荷积累（蓝色等值面）主要分布在 Ga-As 键轴方向，呈明显方向性，且等值面呈哑铃形（类 $p$ 轨道特征），这与共价键的方向性和电子共享特征完全吻合。

注意 GaAs 的 Isovalue 远小于 $\text{CaF}_2$（0.0197 vs 1.0768），说明 GaAs 的电荷转移量小得多，共价成分更强，但因 Ga 和 As 的电负性差异（$\chi_{\text{As}} - \chi_{\text{Ga}} \approx 0.4$），仍具有一定的极性共价键成分（介于纯共价和离子之间）。

（3）两种材料的综合对比

2. 问题提出与讨论

问题：附图中三种能带拓扑情形（a）（b）（c）各对应什么物理体系，对光电器件设计有何影响？

附图展示了在两个不等价 $\mathbf{K}$ 点（$K_+$ 和 $K_-$）处能带的三种拓扑关系：

情形(a)——两点均为正带隙直接型：

$$ E_g^{K_+} > 0,\quad E_g^{K_-} > 0, \quad \mathbf{k}_{\text{VBM}} = \mathbf{k}_{\text{CBM}} $$

两个谷均为直接带隙，代表材料如本实验的 GaAs（$\Gamma$ 点直接带隙）以及六方氮化硼（hBN）。对应高效光电跃迁，无需声子参与，是激光器和 LED 的理想材料。

情形(b)——两点带隙不等，存在间接成分：

$$ E_g^{K_+} \neq E_g^{K_-}, \quad \mathbf{k}_{\text{VBM}} \neq \mathbf{k}_{\text{CBM}} $$

对应 Si、Ge 等间接带隙半导体，以及本实验 $\text{CaF}_2$ 的带隙类型（VBM 在 $\Gamma$，CBM 在 W/X）。光学跃迁需声子辅助，效率低，不适合发光器件，但适合光电探测器（对各频率光响应更宽）。

情形(c)——能带反转（拓扑相变点）：

$$ E_g^{K_+} > 0,\quad E_g^{K_-} < 0 \Rightarrow \text{能带在}\ K_-\ \text{处交叉} $$

这是拓扑绝缘体（Topological Insulator, TI）的标志性能带特征。当体系从情形(a)过渡到情形(c)时，经历拓扑相变，体带隙在某一点关闭并重新打开，带隙符号反转，导致受时间反演对称性 $\mathcal{T}$ 保护的表面/边界态出现：

$$ \mathcal{T}^2 = -1 \Rightarrow \text{Kramers简并保护，表面态无法被无磁散射体打开带隙} $$

代表体系如 $\text{Bi}_2\text{Se}_3$、$\text{HgTe/CdTe}$ 量子阱。对器件设计的意义在于表面态的导电性受拓扑保护，适用于低耗散量子器件和自旋电子学应用。

三种情形与本实验材料的关联：

CaF₂ → 情形(b)：宽带隙间接型
        绝缘体性质，光学透明窗口，
        不适合光电器件

GaAs  → 情形(a)：直接带隙
        高效光-电转换，
        激光器/LED/太阳能电池首选

若对GaAs施加应变或合金化 → 可能出现情形(c)
        即拓扑相变，产生拓扑表面态

作者：GARFIELDTOM
邮箱：coolerxde@gt.ac.cn