[论文解读] A Novel Approach to Spherical Harmonics Expansion for Electron Transport in Semiconductors

S.F. Liotta, Armando Majorana|arXiv (Cornell University)|Mar 28, 2000

Gas Dynamics and Kinetic Theory参考文献 8被引用 5

一句话总结

本文提出了一种新颖的球谐函数展开方法，用于半导体中的玻尔兹曼输运方程，实现了在抛物带和非抛物带近似下电子输运的精确建模。通过在散射算子中引入能量截断并推导渐近解，该方法在高电场下与数值结果高度一致，同时保持了物理一致性，并相比全动能模拟显著降低了计算成本。

ABSTRACT

A set of equations is derived from the Boltzmann kinetic equation describing charge transport in semiconductors. The unknowns of these equations depend on the space-time coordinates and the electron energy. The non-parabolic and parabolic band approximation are treated in detail. In these cases, the set of equations is equivalent to those obtained in the spherical-harmonics expansion. Stationary and homogeneous solutions are explicitly treated. In order to solve numerically the equations, a cut in the energy range is introduced. The modified model maintains the physical characteristics of the original equations. The solution of an asymptotic equation is found and compared to the numerical solutions.

研究动机与目标

• 开发一种计算高效且物理准确的模型，用于亚微米半导体器件中的电子输运，超越漂移-扩散和流体动力学近似。
• 将球谐函数展开技术扩展至非抛物带结构（如Kane模型），同时保持与玻尔兹曼方程的一致性。
• 通过引入物理解释的能量截断，解决碰撞算子中无限能量范围带来的数值挑战。
• 推导并验证高电场区域的渐近解，实现与数值解的对比，并提升计算效率。

提出的方法

• 从玻尔兹曼方程出发，采用一般形式的电子能量，推导出一组偏微分-差分方程，适用于抛物带和非抛物带。
• 对分布函数应用球谐函数展开，将动能方程转化为能量和角度变量上的矩方程系统。
• 在散射率核中引入有限能量截断，使问题在数值上可解，同时保持关键物理特性。
• 使用能量区间离散化方案对所得系统进行数值求解，边界条件确保质量守恒和物理一致性。
• 通过解析近似推导大电场下的渐近解，特别针对速度分布函数。
• 通过匹配分布函数在能量范围内的积分，将渐近解与数值结果进行验证，调整归一化常数以确保一致性。

实验结果

研究问题

• RQ1 如何将球谐函数展开推广以处理电子输运模型中的非抛物带结构？
• RQ2 在碰撞算子中引入能量截断对输运模型的物理保真度和数值稳定性有何影响？
• RQ3 所推导的高电场渐近解与矩方程的数值解匹配程度如何？
• RQ4 该模型在预测饱和速度方面的准确性如何，尤其与蒙特卡罗模拟和实验数据相比？
• RQ5 渐近解在何种条件下成为可靠近似？其与全数值解在不同电场下的表现如何比较？

主要发现

• 该模型成功将球谐函数展开推广至非抛物带，且在抛物带和Kane带近似下与现有框架保持等价。
• 在抛物带情况下，饱和速度达到 $1.27 \times 10^5~{}\text{m/sec}$，略低于实验值 $1.0 \times 10^5~{}\text{m/sec}$，表明仍有改进空间。
• 在非抛物带情况下，电场为 $10^5~{}\text{V/cm}$ 时速度达到 $0.95 \times 10^5~{}\text{m/sec}$，与实验结果的吻合优于抛物带模型。
• 渐近解在高电场下与数值结果高度一致，能量分布和速度分布仅存在微小偏差。
• 能量截断策略保持了物理一致性，通过渐近解与数值解之间分布函数积分的匹配得到验证。
• 在 $w = \alpha$ 处观察到速度分布函数导数的不连续性，归因于方程中的 $\sqrt{w - \alpha}$ 项，表明其与能带非抛物性存在物理关联。