1. 引言
1960年Clough首次引用“有限单元”这一术语,提出一种应用于离散系统的标准求解方法 [1] 。随后数学计算方法的融入,有限元法吸取了差分法对求解域进行离散处理的启示,又继承了里兹法选择试探函数的合理计算方式。20世纪70年代后,有限元法的发展步入工学领域,其“有限单元”的概念来源于工程上的“直接类比”法,是连续介质问题数值解法中最活跃的分支,也成为各类高精度器件设计和性能分析的可靠依据。
有限元法的核心思想是“离散化”和“分片插值”,以简单逼近复杂。当划分的区域足够小、单元数目足够多时,有限单元解收敛于问题的精确解。早期的有限元法的应用集中在各种孤立的物理效应上,但是现实世界的物理现象不会独立依存。随着计算科学的迅速发展,灵巧、简洁、快速的算法和强劲的硬件配置为多物理场耦合打开高速发展的大门。基于有限元法的COMSOL Multiphysics软件可针对大型且复杂的工程问题进行高效地求解并快速产生精确的结果 [2] 。它的特色在于丰富的多物理场耦合,本质上是建模偏微分方程组,任意连接与耦合多物理场方程 [3] 。
近十年来,金属/半导体复合结构的优化设计与改造是光催化领域的研究热点之一。大量的研究工作聚焦于热电子的迁移和肖特基势垒的协同效应,优化复合结构并提高体系的光催化活性。随着界面热电子运动过程的不断探索,热电子和光生电子的悖向运动和往复循环使得金属/半导体间的界面电子流动过程复杂。D. M. Kolb课题组报道了一种热电子直接经界面迁移注入到半导体的隧穿迁移途径。金属/半导体界面间的肖特基势垒处于极小的尺度 [4] ,量子隧道无法被忽视,由此电子的自身能量小于势垒高度时依旧贯穿势垒。Elsharif等 [5] 对热电子的隧穿给出了形象化的比拟,他们将低能态的热电子比作光波,将势垒比作一堵墙,光波在触碰到墙壁时并不会立即消失,而是穿过墙壁。绝缘层的添加分隔金属和半导体,避免热电子的回流和光生电子的悖向运动,有效调节热电子的单向注入效率。然而,入射光、电磁场、半导体器件等多物理场的交织,亟需借助模拟仿真手段来窥探金属/半导体界面的电子流动过程。当前的模拟集中于稳态局域电场的可视化表征,局域电场的分布和强度展示了复合模型的电子产生量情况,但是电场的模拟不能直观地显示电子的迁移轨迹。借助COMSOL软件的电场和半导体场耦合,一种新型的直观表示界面电子运动路径的可视化表征方法被提出。然而,实现可视化表征的前提是构建适用于金属/绝缘层/半导体三明治结构的异质界面的光电子传输的数学物理模型。
基于绝缘中间层的引入和单向隧穿电子路径的设想,本文建立稳态电场和载流子浓度的研究步骤,每个研究步骤分别求解各自的物理场,随即将电场的变量输入到载流子浓度场的相应位置,构建半导体光电子物理场模型,在此模型的基础上加入绝缘层和隧穿的约束边界条件,实现单向传输性能。
2. 半导体光电子物理场耦合过程及隧穿传输的边界条件
金属/半导体热电子器件受光照后,金属内的自由电子“热化”并弛豫,半导体内的光生电子吸收光子能量从价带跃迁至导带。给与光照条件才能使得自由电子的集体运动和电磁场相互作用产生共谐振荡,才能激发半导体内的带间跃迁。为诠释电子的激发现象,半导体光电子学-频域接口将波动光学模块与半导体模块串接耦合,模拟电磁波与半导体间的相互作用。光电子物理场耦合给半导体添加了光跃迁条件,增强了半导体的受激或自发生成速率,体现了半导体材料吸收或发射光引起的电磁场增强或损耗。耦合作用使得光场强度从电磁波界面传递到半导体界面的光跃迁特性中,从而半导体中吸收引起的复折射率或介电函数的变化再次被传回电磁波界面。在模拟中光跃迁条件被添加来体现耦合特性。
2.1. 在光照下半导体内的载流子分布
传播的电磁波随入射方向衰减,出现振荡电场,两个量子态之间随之发生跃迁,相应地产生受激吸收。对于半导体而言,受激吸收过程表现为价带的一个电子吸收一个光子运动至导带,促使光的相干吸收。受激发的电子在导带内迅速平衡,导带的电子产生率G可用基于电子费米能级的费米-狄拉克函数 [6] 来描述:
(1)
式(1)中
——k空间的态密度,通常取
。
——k空间内电子从受激态到激发态的平均矩阵元;
和
可简化k形式,带入式(1)内,整理可得:
(2)
将积分变换为球极坐标,改变积分变量得到:
(3)
将式(3)简化为:
(4)
因为受激吸收率与受激自发发射率直接相关,辐射场的能量在跃迁能量范围内是均匀分布的,所以给定此情况下受激发射和吸收的现象模型 [7] ,如下:
(5)
(6)
其中E——光子能量,
;
是电子从受激态到激发态的吸收过程;
是电子从激发态到受激态的发射过程;
、
——与半导体材料相关的常数;
——单位体积单位能量的平均光子数(J−1∙m−3);
自发发射率不直接依赖于辐射场,基于此,其模型为:
(7)
式(7)中
也是与半导体材料相关的常数。
当半导体与辐射场处于平衡态时,发射速率与吸收速率也相应达到平衡。在此平衡点下,受激吸收、受激发射与自发发射速率的关系为:
(8)
(9)
设定(5)和(6)中
,将其变换带入式(9)中得到:
(10)
电磁场增强引起的高温使得受激辐射超过自发辐射,
值变大,由此
。将其带入式(10)重新整理得到:
(11)
式(11)与普朗克黑体辐射光谱的形式类似,类比可得:
(12)
通常情况下
,c是光在真空中的速度,引入系数C、频率为
的单色线谱对
建立吸收速率模型 [8] 为:
(13)
式(13)中
——光子能量E处的光子密度;
建立分布光谱就需要对一系列线谱求和,设光子驱动能量为
,光子驱动密度
,求和公式可用积分转换为以下形式:
(14)
将式(14)与(5)类比可得:
(15)
对全部频率下的线谱积分得到总电子产生率为:
(16)
类比式(16)和式(4)可得出自发发射寿命与矩阵元的关系:
(17)
而对于常用的n型半导体材料而言,自发发射寿命为常数,将电子占有率等公式带入,整理最终得到:
(18)
当热平衡时,载流子的复合速率与载流子的产生率相等,即为:
(19)
式(19)中电子与空穴的复合速率的求解采用Shockley-Read-Hall模型 [9] ,定义为:
(20)
式(20)中
和
为空穴与电子的寿命;
(21)
(22)
(23)
(24)
式(21)到(24)中
和
——电子与空穴的简并度因子;
和
——导带与价带的有效态密度;
和
——带隙和带隙缩小量;
——陷阱能级。
2.2. 金属接触边界条件的添加
半导体物理场内的金属接触条件被用于表征金属与半导体接触形成的势垒,在真空条件下,此势垒为理想肖特基势垒,利用Crowell和Sze [10] 提出的肖特基接触的简化模型计算。常温下,半导体材料的载流子浓度分布遵从于经典的Boltzamn统计分布,具有非简并性 [11] ,更适用于简化模型的算法。
金属球与半导体颗粒的接触视为载流子的点源,因此电子与空穴对的复合选择为表面复合机制:
(25)
(26)
式(25)和(26)中
——半导体域的外向法线;
和
——半导体内电子和空穴的复合速度;
和
——准平衡态电子和空穴密度(不改变局域能带结构)。
准平衡态情况下,金属的费米能级与半导体在边界处的费米能级相等,由此载流子密度为 [12] :
(27)
(28)
2.3. 薄绝缘栅边界条件的添加
薄绝缘栅边界条件被用于建模金属与半导体之间的绝缘层,绝缘层的厚度极薄时,电场垂直于金属表面的同时也近似垂直于绝缘层/半导体界面。在此条件下,绝缘层/半导体界面处的电位移场
可以定义为:
(29)
式(29)中
——栅极上的势;
——绝缘层/半导体界面的势;
——薄绝缘层厚度;
——薄绝缘体的相对介电函数。
在施加外界偏压V0下的栅极电压为:
(30)
流入薄绝缘栅的电子和空穴的电流必须为零,因此Neumann边界条件 [13] 被应用:
(31)
(32)
2.4. Wentzel-Kramers-Brillouin隧穿方程
金属/半导体异质结界面的电荷流动采用热电子发射模型来表征 [14] ,目的是维持异质结面的电荷连续性。当电子不被约束任意流动时,穿过金属/半导体异质结的电流为:
(33)
(34)
(35)
(36)
式(33)到(36)中
——金属内的外法向电子电流;
——远离半导体的外法向电子电流;
和
——金属的电子与空穴的复合速度;
和
——半导体的电子与空穴的复合速度。
金属/半导体异质结构内的电子与空穴的复合速度被定义为:
(37)
(38)
(39)
(40)
式(37)到(38)中有效Richardson’s系数给出如下:
(41)
(42)
(43)
(44)
式(41)至式(44)中
——金属内电子的质量;
——金属内空穴的质量;
——半导体内电子的质量;
——半导体内空穴的质量。
金属向半导体传递热电子的方式分为两种:电子隧穿和越过势垒。对于金属/绝缘层/半导体三明治结构,加入中间绝缘层的目的是为热电子搭建隧穿通道,因此在模拟中越过势垒这一迁移方式被略去,假设在中间界面的电子传递方式为隧穿。当载流子的能垒厚度等于或小于衰减长度时,量子隧穿效应将发挥重要作用。通过引入Wentzel-Kramers-Brillouin (WKB)隧道模型来模拟异质结或肖特基接触的边界条件,形象化表征异质结构界面的隧穿行为。WKB模型适用于计算任意形状势垒的晶界隧穿电流,考虑的是隧穿电子的动能近似接近于势垒顶部 [15] 。WKB近似假设提出隧穿电流会向热离子电流增加一个分数因子
[16] ,应用于金属/半导体模型内,金属表示为材料1,半导体表示为材料2,在式(33)和式(34)给出的电子与空穴电流密度中引入额外的
,即量子隧穿注入半导体的电子与空穴电流密度近似于:
(45)
(46)
(47)
(48)
这其中分数因子
被定义为:
(49)
式(49)中l——电场线;
Vx——能量轴X;
h——普朗克常数;
——在选择的势垒域中两个相反边界1和2之间的电场线的线积分。
能量积分的限制条件给出如下:
(50)
(51)
式(50)中
和
——跨过势垒的两个相反边界1和2处的电场
值,1和2间通过电场线连接,在选择的势垒域外计算。
在势垒的域选择范围内取式(51)的最大函数,势垒变量
为:
(52)
3. 结论
可视化金属/半导体界面的电子迁移的表征体系的建立,可改良表征电子生成、聚集和运动的模拟手段的单一性。可视化表征体系建立的前提是厘清独立的物理场间的相互连接和作用。本文基于COMSOL软件的多物理场特性,耦合电磁场和半导体,将入射光连结致使半导体内光跃迁行为发生,对基于COMSOL软件中波动光学模块、半导体模块的基础性麦克斯韦方程、漂移–扩散方程,代入金属接触和薄绝缘栅等边界条件,以此推导适用于金属/半导体模型和多物理场耦合条件下的偏微分方程组,建立电场强度、载流子浓度的数学物理模型。
基金项目
江苏省高等学校基础科学(自然科学)研究重大项目(项目编号:22KJA460008)。