[论文解读] Simple remedy for spurious states in discrete k(dot)p models of semiconductor structures

William R. Frensley, Raja N. Mir|arXiv (Cornell University)|Dec 22, 2014

Spectral Theory in Mathematical Physics被引用 4

一句话总结

该论文提出用一阶迎风有限差分格式替代离散k·p模型中的标准中心差分离散化，以消除由二阶格式中非单调色散关系引起的虚假态——这些虚假态是数值误差的产物。该方法保持了物理对称性，提高了对连续能带结构的保真度，确保了绝对稳定性，并在不损失精度的前提下消除了高k值下的振荡解。

ABSTRACT

The spurious states found in numerical implementations of envelope function models for semiconductor heterostructures and nanostructures are artifacts of the use of the centered-difference formula. They are readily removed by employing a first-order difference scheme. The technique produces none of the loss of accuracy that is commonly feared from lower-order formulations; the fidelity to the continuum band structure is actually improved. Moreover the stability is absolute, thus no adjustments to the model parameters are required. These properties are demonstrated in an explicit calculation of a 3-band model of a semiconductor superlattice.

研究动机与目标

• 解决半导体异质结构与纳米结构数值k·p模型中持续存在的虚假态问题。
• 挑战传统观点，即中心差分是k·p离散化中保证精度与厄米特性的必要条件。
• 证明一阶非对称差分可消除虚假解，同时保持物理对称性与稳定性。
• 表明一阶格式可产生单调的能带色散关系，避免了二阶中心差分的反常行为。
• 提供一种鲁棒、无需参数调节的替代方法，避免现有方法依赖精细调参或复杂修改来抑制非物理态。

提出的方法

• 用一阶迎风算子 $\hat{K}_L$ 和 $\hat{K}_R$ 替代标准的中心差分梯度算子 $\hat{K}_C$，其定义为 $\hat{K}_L f_n = -i/\Delta (f_n - f_{n-1})$ 和 $\hat{K}_R f_n = -i/\Delta (f_{n+1} - f_n)$。
• 利用恒等式 $\hat{K}_L \hat{K}_R = \hat{K}_R \hat{K}_L = \hat{L}$，其中 $\hat{L}$ 为离散拉普拉斯算子，确保二阶动能项的一致性。
• 通过 $\hat{H}_1 = \begin{bmatrix} E_C & \hat{K}_R P \\ P \hat{K}_L & E_V \end{bmatrix}$ 构造离散哈密顿量，替代标准的基于 $\hat{K}_C$ 的公式。
• 分析所得色散关系 $K_{L,R}(k)$，其为复数且避免了中心差分 $K_C(k) = \frac{1}{\Delta} \sin(k\Delta)$ 在 $k = \pm \pi/\Delta$ 处的虚假零点。
• 通过2带和3带k·p模型对比一阶与二阶格式的能带结构，表明一阶格式具有单调性且无虚假态。
• 在3带超晶格模型上验证该方法，确认波函数无高k值振荡，且能带保持单调。

实验结果

研究问题

• RQ1 为何在使用中心差分近似的离散k·p模型中会出现虚假态？
• RQ2 一阶差分格式是否可在不损害精度或对称性的情况下消除虚假态？
• RQ3 一阶与中心差分梯度算子的色散关系在能带结构上的影响有何不同？
• RQ4 一阶格式在不同半导体材料与能带模型中是否均保持稳定与准确？
• RQ5 一阶方法能否在避免非对称离散化伪影的同时保持物理对称性？

主要发现

• 一阶迎风差分格式消除了由中心差分格式中非单调色散关系引起的虚假高k态。
• 一阶格式的色散关系 $|K_{L,R}| = \frac{1}{\Delta} \sqrt{2 - 2\cos(k\Delta)}$ 是单调的，且比中心差分的 $K_C(k) = \frac{1}{\Delta} \sin(k\Delta)$ 更接近真实的 $k$-依赖关系。
• 一阶格式在2带与3带模型中均产生单调能带，而中心差分格式则表现出反常（非单调）行为。
• 使用一阶格式计算的波函数无高k值振荡，而中心差分格式的波函数除在能隙外均表现出此类伪影。
• 该方法无条件稳定，无需调节模型参数，而其他方法通常依赖参数调优来抑制虚假态。
• 一阶格式保持了物理对称性与对连续能带结构的保真度，在介观半导体模型背景下相比二阶格式具有更高的精度。