1.Department of Physics, Shanghai University, Shanghai 200444, China 2.National Laboratory for Infrared Physics, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, Shanghai 200083, China
Fund Project:Project supported by the National Natural Science Foundation of China (Grant Nos. 61474073, 61874069)
Received Date:05 March 2021
Accepted Date:06 May 2021
Available Online:07 June 2021
Published Online:05 October 2021
Abstract:The cross-sectional (110) surface of In0.53Ga0.47As/InP hetero-structure grown by molecular beam epitaxy on an InP (001) substrate is characterized by the cross-sectional scanning tunneling microscopy (XSTM). The cleaved (110) surface across the interface between the In0.53Ga0.47As layer and InP layer is atomically flat but displays slight different image contrast between the two neighbor regions. The scanning tunneling spectroscopy (STS) is used to measure the current/voltage (I-V) spectra. The I-V data of the InGaAs surface and InP (110) surface show the different characteristics. The voltage range of zero-current plateau (apparent band gap) in the I-V spectrum of InP displays the values close to its energy band gaps whereas the plateau ranges in the spectra of In0.53Ga0.47As are by contrast generally 50% larger than the energy band gap of In0.53Ga0.47As. The above phenomenon implies the different physical pictures on the tunneling of two surfaces. In the case of InP, the flat band model is feasible since the band edge states existing in the InP (110) surface can prevent the surface from being affected by the tip –induced band bending (TIBB) effect. In contrast, the TIBB effect must be taken into account to explain the I-V spectra of the In0.53Ga0.47As (110) surface. A statistical analysis of the I-V data of In0.53Ga0.47As reveals that the width of current plateau in the I-V spectrum is generally between 1.05 eV and 1.20 eV and the current onset points (turn-points) with the plateau for the different spectra are slightly different from each other. We are able to explain quantitatively the above features based on the three-dimensional TIBB model given by Feenstra (2003 J.Vac. Sci. Technol. B21 2080). Our calculation reveals that the parameter of density of surface states (DOSS) is a sensitive parameter responsible for the I-V features mentioned above. According to an appropriate assignment of the value of DOSS, which is generally taken in the scope of (0.8–3.0) × 1012 (cm2·eV)–1, we well predict both the width and the onset points of the current-plateau. Moreover, the model also reproduces the line-shapes of the I-V spectra measured on In0.53Ga0.47As. Keywords:semiconductor surface/ scanning tunneling microscopy/ scanning tunneling spectroscopy/ tip-induced band bending
全文HTML
--> --> -->
3.结果与讨论图2(a)为InGaAs/InP异质结(110)解理面的STM形貌结果. 形貌图中标出了晶体生长的(001)方向, 可以看到沿该方向在界面两侧呈现出较为明显的衬度变化. 左侧颜色较暗淡的为InP缓冲层, 右侧较明亮的为InGaAs层. 本工作中重点分析这两种外延层的STS测量结果. 在STM实验中, I-V谱是最为常用的STS测量方式. 对于半导体样品, I-V谱在V = 0 V附近所呈现的零电流平台宽度由半导体能带带隙决定. 图2(b)和图2(c)分别给出了InP和InGaAs的典型测量结果(为了突出零电流平台附近谱线特征, 图中仅显示了 ± 1.2 V范围的谱线). 需要指出的是, 虽然很多文献在进行I-V数据分析时习惯对I-V数据作微分处理, 进而以dI/dV谱或归一化微分隧道谱的形式来表示[18], 但上述数据处理过程一般要对原始I-V数据进行必要的平滑处理以降低数值微分所放大的噪声干扰. 该数据处理也会对谱线特征结构(如峰位或拐点)的准确指认有一定影响. 为了避免此类误差, 本文的数据分析基于未做任何处理的原始I-V谱线数据[19]. 另外, 为了尽可能减小数据指认的不确定性以及局域缺陷态叠加产生的影响, 对每一材料的I-V数据分析基于在表面不同空间随机取样十几组数据获得, 并采用如图2(d)和图2(e)所示的表观带隙能量分布直方图的形式(能量间隔0.05 eV)归纳获得统计性测量结果. 图 2 (a) InGaAs/InP异质结(110)解理面的STM形貌; (b)和(c)为InP和InGaAs(110)面的典型的I-V隧道谱给出的表观带隙测量结果; (d)和(e)分别为在InP和InGaAs(110)表面分别随机采样13和17条I-V谱线给出的表观带隙直方图统计结果 Figure2. (a) STM topography of the cleaved (110) surface of InGaAs/InP hetero-structure; (b) and (c) are the typical I-V tunneling spectra of InP and InGaAs (110) surface, respectively; (d) and (e) are the statistical histograms of the apparent tunneling gaps with the sampling of 13 and 17 spectra on the (110) surfaces, respectively.
图2(b)为InP表面给出的一典型I-V谱线, 其零电流平台宽度(隧道谱表观带隙)约为1.11 eV, 还注意到I-V曲线零电流平台正负电压侧对应的拐点位置分别为–0.59和0.52 eV(从平带模型来说它们分别对应价带顶和导带底能量位置, 见图3). InP隧道谱表观带隙略小于其材料带隙1.35 eV. 实验中在InP表面随机采集了13条谱线, 总体上其谱线特征和图2(b)所示均非常类似. 根据这些谱线计算出的零电流平台宽度均值为1.098 eV, 均方差为0.035 eV, 进而可利用正态分布函数更直观地反映带隙能量的统计分布, 如图2(d)中与直方图叠加的蓝色虚线所示. 图 3 平带模型的STM隧道结能带示意图 (a)正偏压; (b)负偏压; TIBB模型下的隧穿能带图 (c)正偏压; (d)负偏压.EC, EV和EF分别是样品的导带、价带和费米能级, EF, tip是针尖的费米能级, Edefect表示缺陷态能级. Figure3. Schematic energy band diagram of STM tunneling junction with flat band model: (a) Positive sample bias; (b) negative sample bias; schematic energy band diagram of tunneling junction with the TIBB model: (c) positive sample bias; (d) negative sample bias. EC, EV and EF are the conduction band, valence band and Fermi level of the sample respectively, EF, tip is the Fermi level of the tip, and Edefect is the defect level.
图2(c)是测量给出的InGaAs(110)面的典型I-V谱线. 隧道谱表观带隙为1.13 eV, 价带顶和导带底能量分别为–0.67和0.46 eV. 图2(e)为用直方图展示的基于17条测量谱线所绘出的表观带隙直方图的统计分布结果. 其均值为1.108 eV, 均方差为0.038 eV. 值得注意的是, 前面InP的隧道谱表观带隙较接近其材料带隙1.35 eV, 但In0.53Ga0.47As表观带隙明显大于材料带隙0.74 eV; 测量值超过材料带隙的50%. 这一不同特征, 意味着这两种外延层的隧穿模型机制是不同的. 对InP的测量结果, 可基于图3(a)和3(b)所示的平带模型能带图进行解释. 平带模型一般适用于材料高掺杂或有高的DOSS引起的费米能级钉扎情形. 平带模型条件下, 样品正偏压下针尖费米能级高于导带底能量EC (如图3(a)所示eV ≥ EC, e为单位电荷), 电子由针尖隧穿至样品. 如定义eVC = EC, 则平带模型下VC对应图3(a)中零电流平台的右侧拐点. 当偏压如图3(b)所示偏压为负时, 针尖费米能级对准价带顶, 发生电子由样品到针尖的隧穿, VV对应I-V谱零电流平台的左拐点. 这样有e(VC + |VV|) = Eg, 即隧道谱的表观带隙应等于材料的带隙Eg. 不过图2(b)和图2(d)的实验测量结果显示, 平均表观带隙约1.10 eV比InP的带隙1.35 eV (室温)要小一些. 这一现象可解释为: 已有关于InP表面态的实验研究发现[20,21], InP(110)表面很容易因氧吸附而产生在距离导带底0.1—0.3 eV的能量区间产生大量深能级, 而本文I-V隧道谱所测量的InP(110)面的表观带隙比实际带隙数值上小约0.2—0.3 eV, 和上述氧吸附产生的深能级态恰在基本相同的数值范围. 由此看来InP表观带隙的窄化或源自上述深能级表面态对隧穿电流的贡献. 对此还要补充说明的是, 虽然前面提到本样品解理是在氮气手套箱保护氛围下进行的, 包括样品到STM真空腔的传递也是通过氮气手套箱与STM进样室对接来避免样品在空气中的暴露, 但是该制备条件尚不能达到对环境分子吸附充分抑制的超高真空条件原位解理所能提供的效果. 对于洁净半导体表面来说, 表面敏感的气体分子的吸附会是迅速的过程, 所以在我们的样品制备条件下InP表面仍存在产生氧吸附表面态的机制. 不过对同一气体分子的吸附敏感性也会因材料而异. 同一样品中InGaAs(110)解理面的STS结果则展示出不同于InP的特点. 如图2(c)的隧道谱显示, 其零电流平台宽度(表观带隙)比材料带隙大出很多, 意味着平带模型对应的物理图像是不适用的. 故对InGaAs隧道谱的解释必须应用TIBB模型而非平带模型机制[8-12]. 先从图3(c)和图3(d)的隧穿能带图就TIBB效应对带隙测量的影响给予定性说明. 当针尖和半导体表面施加电压时, 电场会感应表面处原来电中性的电荷载体(电子、空穴、电离的施主和受主)形成空间电荷的重新分布而对外加电场产生屏蔽作用, 导致外加电势部分偏压降落在半导体上产生所谓TIBB效应: 正偏压下如图3(c)能带向上弯曲; 负偏压下如图3(d)所示能带向下弯曲. 它必然导致大于材料带隙的隧道谱表观带隙, 即e(VC + |VV|) > Eg. TIBB效应下的I-V谱能量宽度可基于电势微分方程的求解给出定量结果[22,23]. 图3是基于隧穿的一维能带图的说明. 2003年Feenstra[8]进一步考虑STM针尖形状的影响, 发展了更符合STM隧穿情形的三维能带弯曲模型. 对于InGaAs计算所用到的参数选取如下: 针尖曲率半径R = 30 nm, 张角θ = 90°, 针尖-样品间距s = 0.75 nm, 针尖(钨)功函数采用${\phi _m}$=4.74 eV[4]. 基于该模型的表面电势分析就可以预言TIBB效应存在时I-V谱的零电流平台宽度(表观带隙)的大小. 如图4所示, 计算出的InGaAs表面电势能随所加偏压的变化曲线如图中实线所示[24]. 图4的横坐标为针尖-样品的偏压, 纵坐标为以样品费米能级作为零参考点的能量变化值; 图中两组实线分别为表面处导带和价带电势随所加偏压的变化情况, 除了浅蓝色实线为未考虑表面态修正情形, 红、橙、浅绿和蓝色实线分别为考虑不同DOSS修正后情形, 稍后讨论. 另外, 图中黑色虚线是针尖费米能级随所加偏压改变而产生的与样品费米能级差的变化. 它与表面处价带或导带的交点为产生隧道电流的阈值点, 即I-V谱线上零电流平台的左、右拐点电压值(分别记为V–onset和V+onset). 二者差值ΔV给出的表观带隙取值(eΔV). 图 4 三维TIBB模型给出的InGaAs表面处带边电势能随所加样品偏压的变化关系, 不同颜色的实线对应于不同DOSS, 对应箭头所指为该态密度下计算得到的表观带隙, 优化的DOSS对应计算结果如红色实线所示 Figure4. Variation of surface potential of InGaAs with the sample voltage calculated by the 3D TIBB model. The solid lines of different colors correspond to different DOSS, and the corresponding arrows indicate the calculated apparent tunneling gaps at the DOSS. The optimized calculation result is shown by the red line.
Calculations with the density of surface states/(cm2·eV)–1
0
1.25
–0.89
0.36
1.0 × 1011
1.29
–0.89
0.40
1.0 × 1012
1.21
–0.75
0.46
2.0 × 1012
1.11
–0.66
0.45
1.0 × 1013
0.90
–0.49
0.41
表1InGaAs隧道谱特征参量的实验数据与TIBB模型计算的对比 Table1.Comparison of experimental data and TIBB model calculations on the characteristic parameters of InGaAs tunneling spectra.
注意到Feenstra[11,12]在TIBB模型中引入过DOSS这一唯像参量修正. 模型中该修正对于能量位于禁带中的表面态分布采取如下简化处理, 即中性能级位于禁带中央, 禁带中的施主和受主能态分别在中性能级上下均匀分布. 而DOSS大小的取值目前文献中一般在1011—1013 (cm2·eV)–1量级[25]. 为了考察DOSS取值对I-V谱的影响, 选取1.0 × 1011—1.0 × 1013 (cm2·eV)–1态密度范围内几组典型取值进行计算. 各组取值下的电势能变化如图4所示. 如前所述, 每组的价带和导带带边与针尖费米能量(黑色虚线)交点间的电压间隔预言了相应DOSS下I-V谱的零电流平台宽度(STS表观带隙), 如图4中不同颜色水平箭号线所指示, 具体数值见表1. 为了与实验比较, 图中实验测量给出的表观带隙用黑色水平箭头标明. 如图4上方所标示, 当DOSS取为2.0 × 1012 (cm2·eV)–1时所给出的计算结果(红色线)与实验最为一致: 此时不仅表观带隙值最接近实验结果, 而且对零电流平台起、止能量位置的预测也与实验吻合最好. 另外, 通过比较不同DOSS下的计算结果, 观察到DOSS增大时表观带隙展宽效应变小并趋近材料带隙. 例如当取较高DOSS (1.0 × 1013 (cm2·eV)–1)时, 表观带隙值为0.90 eV已经很接近材料带隙值(室温0.74 eV). 事实上, 对上述模型结果的更充分的验证还可以通过与InGaAs的I-V谱的理论计算结果比较来获得[11,12]. 图2(c)的实验测量谱线与计算谱线放到一起的比对如图5(a)所示. 其中, 红色虚线为计算结果. 图5(b)则是测量谱线(黑线)与计算谱线在更宽测量范围即(–1.6 V, 1.6 V)区间的比对. 可以看到, 无论从图5(a)所展示的细节特征还是图5(b)对整条谱线的线形趋势特征的预言, 计算结果均与实验吻合得很好. 图 5 (a) InGaAs的I-V谱的计算(红色虚线)结果与实验谱线(黑色实线)的对比; (b) 图(a)谱线在整个(–1.6 V, 1.6 V)电压范围对计算/实验对比结果的呈现 Figure5. (a) Comparison of the I-V spectra between the calculated (red dashed line) and experimental (black solid line); (b) comparison of the spectra within the whole voltage range of (–1.6 V, 1.6 V).
分析认为, 以图5中特定的I-V谱线进行的讨论说明是代表性的. 结合图2(e)的测量统计结果, I-V谱线的零电流平台宽度及其电压起止位置等参数对于来自表面不同位置的测量谱线会略有差异, 但是采用同样方法对其他谱线进行测试计算, 发现只要适当调节DOSS取值, 可以获得对相应实验谱线较满意的拟合结果. 其道理也不难解释, 因为STS给出的是局域电子态信息, 不同位置的I-V线形难免也存在局域性的表面态的调控, I-V谱线形特征自然也携带着具有空间局域特征的表面态信息. 作为进一步的说明, 在图6(a)—(d)中又给出了InGaAs解理面上另外4组测量谱线(黑色)与计算(红色虚线)放到一起比对的结果, 可以看到二者的吻合均令人满意. 图 6 (a)—(d) InGaAs的四组I-V谱实验谱线(黑色实线)与模型计算(红色虚线)的比对, 各图中分别注明了计算时采用的表面态密度(DOSS)数值 Figure6. (a)–(d) Comparison of four groups of experimental I-V spectra (black solid lines) with the corresponding calculated I-V spectra (red dashed line). The DOSS for each curve is indicated, respectively.