搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

a-Si:H薄膜中SiyHx结构组态的原子模拟研究

翟世铭 廖黄盛 周耐根 黄海宾 周浪

HE KAI-FEN, A. SALAT. WINDING NUMBER BIFURCATION AND SUPERSTRUCTURE IN NONLINEAR DRIFT WAVES DRIVEN BY A TRAVELING WAVE. Acta Phys. Sin., 1990, 39(2): 204-211. doi: 10.7498/aps.39.204
Citation: HE KAI-FEN, A. SALAT. WINDING NUMBER BIFURCATION AND SUPERSTRUCTURE IN NONLINEAR DRIFT WAVES DRIVEN BY A TRAVELING WAVE. Acta Phys. Sin., 1990, 39(2): 204-211. doi: 10.7498/aps.39.204

a-Si:H薄膜中SiyHx结构组态的原子模拟研究

翟世铭, 廖黄盛, 周耐根, 黄海宾, 周浪

Atomic simulation of SiyHx structure configuration in a-Si:H thin films

Zhai Shi-Ming, Liao Huang-Sheng, Zhou Nai-Gen, Huang Hai-Bin, Zhou Lang
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 氢化非晶硅薄膜(a-Si:H)中SiyHx结构组态对薄膜应用性能有重要影响, 然而现有的分析测试手段难以对其进行深入细致的研究. 本文运用分子动力学方法模拟分析了a-Si:H/c-Si薄膜中SiyHx结构组态, 以及衬底温度对其含量的影响; 并进一步运用第一性原理方法计算了各SiyHx组态中的Si-H键能. 结果发现a-Si:H薄膜中SiyHx结构可以归纳为六种组态. 三类为以化学键结合的SiHx组态, 包括SiH, SiH2和SiH3; 另外三类为以物理键结合的HSiy组态, 包括HSi2(s), HSi2(l)和HSi3. 键能结果反映出六种组态的稳定性由高到低的顺序为SiH > SiH2 > SiH3 > HSi2(s) > HSi2(l) > HSi3. HSiy组态中Si—H键能在太阳光中的可见光和红外线的能量范围内, 阳光照射引起HSiy组态中的Si-H物理键断裂, 是非晶硅薄膜电池产生S-W (Steabler-Wronski)效应的主要机理. 另外, 薄膜沉积过程中衬底温度的升高将导致各类SiyHx组态含量大幅降低.
    The hydrogenated amorphous silicon (a-Si:H) film is the core structure of hetero junction with intrinsic thin layer solar cell. Its quality determinates the photoelectric conversion efficiency of this solar cell directly. The configuration of SiyHx is an important structure characteristic of a-Si:H films, and it can influence on the quality of a-Si:H thin films and their application properties. However, it is difficult to study them in depth and detail by the existing analytical and testing methods. In this paper, the structure configuration of SiyHx in a-Si:H /c-Si thin films and the effect of substrate temperature on its content have been simulated and analyzed by molecular dynamics method. A modified Tersoff potential developed by Murty was used to calculate the inter-atomic forces. The results showed that the SiyHx structure in a-Si:H thin films can be summarized into six configurations. Three traditional configurations, including SiH, SiH2 and SiH3, can be referred to as SiHx configurations.The other three nove configurations, including HSi2(s), HSi2(l) and HSi3, can be referred to as HSiy configurations. The main differences between the configurations of HSi2(l) and HSi2(s) are the longer Si—H bonds and bigger bond angle in HSi2(l) configuration than those in HSi2(s) configuration. All of the Si-H bonds in SiHx configurations are strong chemical bonds, while the Si—H bonds in HSiyconfigurations are weak physical bonds. The further calculations of the Si-H bond energies in six configurations have been carried out by the first principle method. According the bond energies results, we can deduce that the order of the stability of six configurations from high to low is SiH > SiH2 > SiH3 > HSi2(s) > HSi2(l) > HSi3. Comparing the Si—H bond energies of the six configurations with the solar energy, it is found that the Si-H bond energy in HSiy configuration is in the range of visible and infrared light in solar light. Si—H physical bonds are easy to fracture in HSiy configuration caused by solar light. This may be the main mechanism of producing Steabler-Wronski (S-W) effect in amorphous silicon thin film cells. In addition, the rise of substrate temperature in the deposition process of a-Si:H films will lead to a significant decrease in the configuration content of all kinds of SiyHx configurations.
      通信作者: 周耐根, ngzhou@ncu.edu.cn ; 黄海宾, haibinhuang@ncu.edu.cn
    • 基金项目: 国家级-氢化非晶硅薄膜局域结构特征及其形成的分子动力学模拟研究(51561022)
      Corresponding author: Zhou Nai-Gen, ngzhou@ncu.edu.cn ; Huang Hai-Bin, haibinhuang@ncu.edu.cn

    氢化非晶硅薄膜(a-Si:H)具有良好的光电特性. 在太阳能电池领域, 其与晶体硅形成的异质结作为基础的非晶硅/晶体硅太阳能电池保持了晶体硅太阳能电池的最高转换效率的世界纪录26.7%(单面结构)[1]和24.7%(双面结构)[2], 与传统的晶体硅太阳能电池相比, 非晶硅/晶体硅太阳能电池还具有优良的温度特性[3,4]、较高的开路电压、高稳定性[5]等优势, 被业内公认为是下一代适宜规模生产、市场竞争力强的太阳能电池.

    但我们对a-Si:H的认识仍不够深入, 其为非晶结构, 其结构和性能受制备过程中的工艺参数, 如衬底温度[6-8]、入射动能[9]、入射频率[10,11]等影响显著, 且其结构在后继使用和热处理过程中会发生变化, 典型的如光致衰减效应(Steabler-Wronski效应, S-W效应)[12]、热处理提高本征非晶硅层的钝化效果等[13,14]. 普遍认为薄膜中特殊的Si—H键中弱键断裂、迁移等机制是上述现象的可能原因[15]. 但在实验研究方面受限于制样和观测手段, 只能对薄膜的整体结构进行宏观统计上的分析, 难以得到进一步的微观结构信息, 更无法对不同工艺条件下制备的薄膜结构进行详细分类, 深入理解其性能随工艺演变的本质. 分子动力学模拟基于经典力学, 可追踪薄膜每个原子在相空间的运动规律和轨迹, 以及便于统计分析原子结构包括成键信息等优点, 是研究薄膜沉积生长和结构分析的有力工具. 本课题组前期利用分子动力学模拟研究了沉积参数对a-Si:H薄膜中悬挂键含量、H含量、晶化率和SiHx键合等微观结构的影响, 给出了本征非晶硅钝化晶体硅的钝化效果在适度的温度下热处理可获得改善的合理解释[6,9]. 本文在该研究基础上进一步模拟分析了a-Si:H/c-Si薄膜中SiyHx组态的种类、结构、稳定性, 及衬底温度对其含量的影响, 以期对非晶硅薄膜的结构和性能随工艺的演变规律有更进一步的认识.

    图1为a-Si:H/c-Si (001)薄膜生长的分子动力学模型示意图. 衬底尺寸为15a × 15a × 5.3a(a为硅的晶格常数), 共约9900个原子. 在Z方向, 衬底自下而上被分为三层: 1)底部为固定原子层, 以防止沉积过程中衬底模型受入射原子的冲击而飘移; 2)中间部分为恒温层, 在该层中采用Langevin方法控制其温度, 即为衬底温度, 防止由于沉积过程中入射原子的动能导致衬底温度的升高; 3)顶部为自由层, 在该部分的原子处于自由运动状态, 不对其施加其他约束条件以避免恒温算法可能引入的影响. 在XY方向上运用周期性边界条件, 用以模拟薄膜在无限大的表面上生长. 原子间相互作用采用由Murty和Atwater[16]在Tersoff势[17]的基础上修改得到的Murty-Tersoff势函数来描述, 时间步长设为0.001 ps. 现有研究表明, SiH3基团是获得高质量薄膜的主要自由基团, 因此仅采用SiH3基团进行沉积. 沉积时SiH3基团在衬底上方一定高度的空间区域内随机产生, 再垂直入射到衬底表面. 结构分析时, 体系将快冷至约0.01 K, 以消除热振动的影响.

    图 1 晶体硅衬底表面沉积生长a-Si:H/c-Si薄膜的模型结构示意图\r\nFig. 1. Schematic diagram of the model structure about a-Si:H thin film deposited on the surface of crystalline silicon substrate.
    图 1  晶体硅衬底表面沉积生长a-Si:H/c-Si薄膜的模型结构示意图
    Fig. 1.  Schematic diagram of the model structure about a-Si:H thin film deposited on the surface of crystalline silicon substrate.

    SiyHx组态结构的键合状态、键长、键角均采用Mercury软件[18]进行分析. Mercury作为一款可视化软件, 可以提供体系所有原子成键的键长、键角信息. SiyHx组态中Si—H键的键能采用VASP软件[19]进行计算, 计算模型中Si—H键的键长键角均采用各SiyHx组态统计得到的平均值.

    在不同沉积参数下模拟沉积a-Si:H, 得到的a-Si:H/c-Si薄膜原子结构基本类似. 图2示出了模拟衬底温度为500 K时沉积得到a-Si:H/c-Si薄膜原子结构. 图中原子分为两部分, 底部为衬底层, 原子排列规则; 顶部为沉积生长的a-Si:H薄膜层, 原子排列混乱无序, 二者界面处在一定范围内Si原子从规则排列渐变为无序排列. 说明在模拟考察的沉积生长参数下, 均能形成非晶结构的薄膜, 没有出现结晶现象. 其中棕色原子和较小的红色原子分别为沉积的硅原子和氢原子, 可以发现氢原子主要分布在薄膜层的内部, 但因其半径小, 少部分可扩散进入衬底层. 图2(b)为薄膜层的径向分布函数. 图中径向分布函数曲线的第一峰为尖锐的高峰, 后续第二峰、第三峰非常弱, 说明该结构短程有序长程无序, 为典型的非晶结构的径向分布函数曲线, 进一步说明了薄膜为非晶结构.

    图 2 衬底温度为500 K时, 入射约8000个基团得到的a-Si:H/c-Si薄膜的原子结构示意图(a) 和薄膜层的RDF图(b)\r\nFig. 2. Atomic structure of a-Si:H/c-Si film after depositing about 8000 groups with the substrate temperature of 500 K (a) and RDF curve of the film layer (b).
    图 2  衬底温度为500 K时, 入射约8000个基团得到的a-Si:H/c-Si薄膜的原子结构示意图(a) 和薄膜层的RDF图(b)
    Fig. 2.  Atomic structure of a-Si:H/c-Si film after depositing about 8000 groups with the substrate temperature of 500 K (a) and RDF curve of the film layer (b).

    为研究H原子在a-Si:H/c-Si薄膜中的存在状态, 我们统计了Si-H键的键长分布, 如图3所示. 在图中Si-H键分布有两个明显的主峰, 以及之后两部分较矮的峰群, 分别处于1.42—1.49 Å, 1.49—1.57 Å, 1.57—1.69 Å和1.69—1.74 Å. 这可能对应于不同的SiyHx组态.

    图 3 a-Si:H薄膜中Si-H键的键长分布 (沉积条件: 衬底温度为500 K, 入射动能为1.45 eV, 入射频率为1 ps/个基团)\r\nFig. 3. The bond length distribution of the Si-H bond in the a-Si:H film (Deposition conditions: substrate temperature is 500 K, the incident kinetic energy is 1.45 eV, and the incident frequency is 1 ps/group).
    图 3  a-Si:H薄膜中Si-H键的键长分布 (沉积条件: 衬底温度为500 K, 入射动能为1.45 eV, 入射频率为1 ps/个基团)
    Fig. 3.  The bond length distribution of the Si-H bond in the a-Si:H film (Deposition conditions: substrate temperature is 500 K, the incident kinetic energy is 1.45 eV, and the incident frequency is 1 ps/group).

    为进一步分析a-Si:H/c-Si薄膜中的SiyHx组态, 利用Mercury软件对模拟沉积生长得到的a-Si:H薄膜中的Si—H的键合状态、键长、键角进行了统计分析, 发现H原子在连接不同数量的Si原子时, 其Si—H键的键长、键角呈现出不同的特征. 根据传统价键理论, 在Si—H结构中, 单个H原子通常只与一个Si原子形成化学键, 结合成三种组态: 分别为SiH, SiH2和SiH3, 记为SiHx组态. 这三类结构很容易鉴别出来, 其结构示于图4(a)图4(c)中. 这三种组态中Si-H键的键长相近, 集中分布在1.46 Å附近, 对应于键长分布图(图3)中的第一峰, 在a-Si:H/c-Si薄膜中各SiyHx组态中长键角所占的比率最大. 除了上述三种化学键结合的SiHx组态外, 我们还发现了H与周围的Si原子形成物理键, 以氢桥Si-H-Si(HSi2)或HSi3的形式出现, 记为HSiy. a-Si:H薄膜中氢桥的存在曾被报道过[20,21], 但对于氢桥的结构特征分析不够明晰. 我们进一步分析还发现HSi2结构中存在两种具有不同键长键角特征的组态, 一种组态的Si—H键长处在1.49—1.57 Å范围, 键角集中分布在123°—137°之间, 因Si—H键长较短, 记为HSi2(s), 如图4(d)所示. 另一种组态Si—H键分布在1.5 —1.68 Å范围, Si—H—Si键角大于140°, 记为HSi2(l), 如图4(e)所示. 不难发现HSi2(s)和HSi2(l)两种组态对应于键长分布图(图3)中的第二峰和第三峰群, 其中HSi2(s)的含量明显高于HSi2(l). 图4(f)示出的是HSi3结构组态. 它的Si-H键在所有组态中最长, 大于1.7 Å. 同时它的Si-H-Si键角最小, 小于118°. HSi3组态对应于键长分布图(图3)中的最后一个小峰, 在所有组态中占比最少.

    图 4 a-Si:H薄膜中SiyHx组态的结构示意图 (a) SiH; (b) SiH2; (c) SiH3; (d) HSi2(s); (e) HSi2(l); (f) HSi3\r\nFig. 4. The structure Schematic diagram of SiyHx in the a-Si:H film: (a) SiH; (b) SiH2; (c) SiH3; (d) HSi2(s); (e) HSi2(l); (f) HSi3.
    图 4  a-Si:H薄膜中SiyHx组态的结构示意图 (a) SiH; (b) SiH2; (c) SiH3; (d) HSi2(s); (e) HSi2(l); (f) HSi3
    Fig. 4.  The structure Schematic diagram of SiyHx in the a-Si:H film: (a) SiH; (b) SiH2; (c) SiH3; (d) HSi2(s); (e) HSi2(l); (f) HSi3.

    为分析以上六种SiyHx组态的稳定性, 我们采用第一性原理方法计算了它们的Si-H键能, 计算时各组态模型中的键长键角均采用统计得到的平均值. 同时将这些键能与太阳光的光子能量进行对比, 结果示于图5中. 通过键能的高低可以发现它们的稳定性由高到低的顺序为: SiH > SiH2 > SiH3 > HSi2(s) > HSi2(l) > HSi3. 只有SiH键的键能超出太阳光能量区间, 而SiH2, SiH3, HSi2(s), HSi2(l), HSi3组态的键能都处于该能量区间内. 但在实际非晶硅薄膜电池结构中, 非晶硅表面还有一层铟锡氧化物半导体透明导电膜(ITO)或加上一层玻璃层. 实验研究表明, ITO层或玻璃层几乎可以吸收掉太阳光中的全部紫外线[22]. 也就是说, 只有可见光和红外线可以真正照射到a-Si:H薄膜上. 从能量上正好对应于Si-H物理键的三种组态结构HSiy.

    图 5 各SiyHx组态中Si-H键能与太阳光能量对比图\r\nFig. 5. Comparisons of Si-H bond energy in each SiyHx configuration with solar energy
    图 5  各SiyHx组态中Si-H键能与太阳光能量对比图
    Fig. 5.  Comparisons of Si-H bond energy in each SiyHx configuration with solar energy

    可以想象一下太阳光照射到电池上的场景: 首先紫外线被ITO或玻璃层全部吸收, 只有可见光和红外线透过ITO层到达a-Si:H薄膜表面. 它们因能量太低, 不会引起SiHx组态中Si—H键断键, 也就不会产生自由H原子; 但它们可与三种HSiy组态中Si—H物理键产生交互作用, 引起某个Si-H物理键断裂. 这种物理键的断裂, 只会让H原子位置产生偏移, 远离原成键Si原子, 偏向其余的Si原子, 同时引起Si原子电荷中心的变化. 对于断键的Si原子来说, 断键将增强其价电子的不饱和度, 类似于产生了一个强度较弱的Si悬挂键, 这将导致电池性能的下降. 因为光照断键的H原子仍被原HSiy组态中的其他Si原子约束, 不易发生长距离的扩散, 所以在没有光照的条件下, 经过结构弛豫, 原来的断裂的Si—H键又可能重新成键, 恢复电池的使用性能. 这个现象与非晶硅薄膜电池应用过程中的S-W效应的演化规律非常符合. 对于S-W效应起因, 人们已经做了大量的工作, 并提出了多种微观模型试图解释它, 但至今仍没有形成统一的观点. 而根据前面的分析结果, 我们认为光照引起HSiy组态中的Si—H物理键断裂, 才是S-W效应的主要机理.

    为了探讨沉积生长条件对SiyHx组态的影响规律, 我们选择衬底温度这一对薄膜结构影响较大的条件进行考察. 控制模拟参数: 衬底温度分别为300, 400, 500, 600和700 K, 基团入射频率为1个/ps, 入射动能为1.45 eV, 在Si(100)面上沉积约8000个基团. 然后统计了薄膜中悬挂键、SiHx和HSiy含量随衬底温度的变化, 如图6(a)所示. 该图以及后续图中各组态的含量都以其中的最大值为参考进行了归一化处理. 在不同衬底温度下, SiHx含量始终远大于HSiy的含量. 随温度升高, 薄膜中的SiHx和HSiy含量呈下降趋势, 这是因为衬底温度越高, 原子热振动加剧, Si—H弱键断裂形成氢气离开薄膜概率增加, 导致H原子含量降低的结果. 在500 K温度以下, 随温度升高, 悬挂键含量迅速降低; 500 K后, 下降趋势变得平缓. 这是因为虽然温度升高导致Si—H弱键断裂, H原子溢出, 会形成更多的Si悬挂键; 但高温有利于Si悬挂键之间也发生重组, 形成Si—Si强键, 使得Si悬挂键含量降低, 后者的作用大于前者, 总体上使得薄膜中Si悬挂键含量的降低.

    图 6 不同衬底温度下沉积生长的a-Si:H/c-Si薄膜中的悬挂键、SiHx和HSiy相对含量(a), SiH, SiH2和SiH3的相对含量(b)和HSi2(s), HSi2(l)和HSi3相对含量(c)\r\nFig. 6. Relative contents of dangling bonds, SiHx and HSiy(a), relative contents of SiH, SiH2 and SiH3 (b), and relative content of HSi2(s), HSi2 (l) and HSi3 (c) in a-Si:H/c-Si films deposited with different substrate temperatures.
    图 6  不同衬底温度下沉积生长的a-Si:H/c-Si薄膜中的悬挂键、SiHx和HSiy相对含量(a), SiH, SiH2和SiH3的相对含量(b)和HSi2(s), HSi2(l)和HSi3相对含量(c)
    Fig. 6.  Relative contents of dangling bonds, SiHx and HSiy(a), relative contents of SiH, SiH2 and SiH3 (b), and relative content of HSi2(s), HSi2 (l) and HSi3 (c) in a-Si:H/c-Si films deposited with different substrate temperatures.

    进一步分析了SiHx的三种组态和HSiy三种组态含量随衬底温度的变化, 分别示于图6(b)图6(c)中. 从图6(b)中可发现, SiHx键合方式始终以SiH键为主, SiH2含量次之, SiH3含量最少. SiHx键含量随衬底温度的升高而降低. 同时SiH3含量随衬底温度的升高也有大幅的减少. 而HSiy中主要为HSi2(s)组态, HSi2(l)和HSi3含量很少. 随衬底温度的升高, 三种组态含量呈现明显的下降趋势. HSi2(l)在400 K之前含量迅速降低, 之后下降趋势变缓. HSi3在500 K后含量降为0, 这可能因为HSi3键能最低, 在高温下最不稳定.

    本文采用分子动力学方法, 模拟了在不同沉积条件下得到a-Si:H/c-Si薄膜, 并对其中SiyHx组态进行了结构分析, 发现在a-Si:H/c-Si薄膜中除常见的以化学键结合的SiHx组态外, 还存在以物理键结合的HSiy组态. 根据H原子连接Si原子的数量、Si—H键长及键角特征, 我们把SiyHx组态分为SiH, SiH2, SiH3, HSi2(s), HSi2(l)和HSi3六种类型, 前三种为SiHx组态, 后三种为HSiy组态. 通过第一性原理计算各组态中的Si—H键能发现, 它们的稳定性由高到低的顺序为: SiH > SiH2 > SiH3 > HSi2(s) > HSi2(l) > HSi3. 且HSiy组态中Si—H键能在太阳光中的可见光和红外线的能量范围内. 阳光照射引起HSiy组态中的Si—H物理键断裂, 可能是非晶硅薄膜电池产生S-W效应的主要机理. 另外, 薄膜沉积过程中衬底温度的升高将导致各类SiyHx组态含量的大幅降低.

    [1]

    Yoshikawa K, Kawasaki H, Yoshida W, Irie T, Konishi K, Nakano K, Uto T, Adachi D, Kanematsu M, Uzu H, Yamamoto K 2017 Nat. Energy 2 17032Google Scholar

    [2]

    Liu J, Huang S H, He L 2015 J. Semicond. 36 4Google Scholar

    [3]

    Mishima T, Taguchi M, Sakata H, Maruyama E 2011 Sol. Energy Mater. Sol. Cells 95 18Google Scholar

    [4]

    Li Z, Zhang X W, Han G R 2010 Phys. Status Solidi A 207 144Google Scholar

    [5]

    Masuko K, Shigematsu M, Hashiguchi T, Fujishima D, Kai M, Yoshimura N, Yamaguchi T, Ichihashi Y, Mishima T, Matsubara N, Yamanishi T, Takahama T, Taguchi M, Maruyama E, Okamoto S 2014 IEEE. J. Photovoltaics 4 1433Google Scholar

    [6]

    Luo Y R, Gong H Y, Zhou N G, Huang H B, Zhou L 2018 Appl. Phys. A 124 18Google Scholar

    [7]

    Illiberi A, Creatore M, Kessels W M M, de Sanden M 2010 Phys. Status SolidiR 4 206Google Scholar

    [8]

    Bronsveld P C P, Mates T, Fejfar A, Kocka J, Rath J K, Schropp R E I 2010 Phys. Status Solidi A 207 525Google Scholar

    [9]

    Luo Y, Zhou N, Gong H, Huang H, Lang Z 2018 IOP Conf. Ser.: Mater. Sci. Eng. 284 012006Google Scholar

    [10]

    Hou G F, Fan Q H, Liao X B, Chen C Y, Xiang X B, Deng X M 2011 J. Vac. Sci. Technol., A 29 061201Google Scholar

    [11]

    Andujar J L, Bertran E, Canillas A, Roch C, Morenza J L 1991 J. Vac. Sci. Technol., A 9 2216Google Scholar

    [12]

    Staebler D L, Wronski C R 1977 Appl. Phys. Lett. 31 292Google Scholar

    [13]

    Zhang D, Tavakoliyaraki A, Wu Y, van Swaaij R, Zeman M 2011 Energy Procedia 8 207Google Scholar

    [14]

    Alnuaimi A, Islam K, Nayfeh A 2013 Sol. Energy 98 236Google Scholar

    [15]

    秦国刚, 孔光临 1988 半导体学报 01 103Google Scholar

    Qin G G, Kong G L 1988 J. Semicond. 01 103Google Scholar

    [16]

    Murty M V R, Atwater H A 1995 Phys. Rev. B 51 4889Google Scholar

    [17]

    Tersoff J 1988 Phys. Rev. B 37 6991Google Scholar

    [18]

    Macrae C F, Bruno I J, Chisholm J A, Edgington P R, Wood P A 2008 J. Appl. Crystallogr. 41 466Google Scholar

    [19]

    Hafner J, Kresse G 1997 Properties of Complex Inorganic Solids (New York: Plenum Press) pp69–99

    [20]

    Sriraman S, Agarwal S, Aydil E S, Maroudas D 2002 Nature 418 62Google Scholar

    [21]

    Robertson J 2000 J. Non-Cryst. Solids 266 79Google Scholar

    [22]

    Kim H, Horwitz J S, Kushto G, Piqué A, Kafafi Z H, Gilmore C M, Chrisey D B 2000 J. Appl. Phys. 88 6021Google Scholar

  • 图 1  晶体硅衬底表面沉积生长a-Si:H/c-Si薄膜的模型结构示意图

    Fig. 1.  Schematic diagram of the model structure about a-Si:H thin film deposited on the surface of crystalline silicon substrate.

    图 2  衬底温度为500 K时, 入射约8000个基团得到的a-Si:H/c-Si薄膜的原子结构示意图(a) 和薄膜层的RDF图(b)

    Fig. 2.  Atomic structure of a-Si:H/c-Si film after depositing about 8000 groups with the substrate temperature of 500 K (a) and RDF curve of the film layer (b).

    图 3  a-Si:H薄膜中Si-H键的键长分布 (沉积条件: 衬底温度为500 K, 入射动能为1.45 eV, 入射频率为1 ps/个基团)

    Fig. 3.  The bond length distribution of the Si-H bond in the a-Si:H film (Deposition conditions: substrate temperature is 500 K, the incident kinetic energy is 1.45 eV, and the incident frequency is 1 ps/group).

    图 4  a-Si:H薄膜中SiyHx组态的结构示意图 (a) SiH; (b) SiH2; (c) SiH3; (d) HSi2(s); (e) HSi2(l); (f) HSi3

    Fig. 4.  The structure Schematic diagram of SiyHx in the a-Si:H film: (a) SiH; (b) SiH2; (c) SiH3; (d) HSi2(s); (e) HSi2(l); (f) HSi3.

    图 5  各SiyHx组态中Si-H键能与太阳光能量对比图

    Fig. 5.  Comparisons of Si-H bond energy in each SiyHx configuration with solar energy

    图 6  不同衬底温度下沉积生长的a-Si:H/c-Si薄膜中的悬挂键、SiHx和HSiy相对含量(a), SiH, SiH2和SiH3的相对含量(b)和HSi2(s), HSi2(l)和HSi3相对含量(c)

    Fig. 6.  Relative contents of dangling bonds, SiHx and HSiy(a), relative contents of SiH, SiH2 and SiH3 (b), and relative content of HSi2(s), HSi2 (l) and HSi3 (c) in a-Si:H/c-Si films deposited with different substrate temperatures.

  • [1]

    Yoshikawa K, Kawasaki H, Yoshida W, Irie T, Konishi K, Nakano K, Uto T, Adachi D, Kanematsu M, Uzu H, Yamamoto K 2017 Nat. Energy 2 17032Google Scholar

    [2]

    Liu J, Huang S H, He L 2015 J. Semicond. 36 4Google Scholar

    [3]

    Mishima T, Taguchi M, Sakata H, Maruyama E 2011 Sol. Energy Mater. Sol. Cells 95 18Google Scholar

    [4]

    Li Z, Zhang X W, Han G R 2010 Phys. Status Solidi A 207 144Google Scholar

    [5]

    Masuko K, Shigematsu M, Hashiguchi T, Fujishima D, Kai M, Yoshimura N, Yamaguchi T, Ichihashi Y, Mishima T, Matsubara N, Yamanishi T, Takahama T, Taguchi M, Maruyama E, Okamoto S 2014 IEEE. J. Photovoltaics 4 1433Google Scholar

    [6]

    Luo Y R, Gong H Y, Zhou N G, Huang H B, Zhou L 2018 Appl. Phys. A 124 18Google Scholar

    [7]

    Illiberi A, Creatore M, Kessels W M M, de Sanden M 2010 Phys. Status SolidiR 4 206Google Scholar

    [8]

    Bronsveld P C P, Mates T, Fejfar A, Kocka J, Rath J K, Schropp R E I 2010 Phys. Status Solidi A 207 525Google Scholar

    [9]

    Luo Y, Zhou N, Gong H, Huang H, Lang Z 2018 IOP Conf. Ser.: Mater. Sci. Eng. 284 012006Google Scholar

    [10]

    Hou G F, Fan Q H, Liao X B, Chen C Y, Xiang X B, Deng X M 2011 J. Vac. Sci. Technol., A 29 061201Google Scholar

    [11]

    Andujar J L, Bertran E, Canillas A, Roch C, Morenza J L 1991 J. Vac. Sci. Technol., A 9 2216Google Scholar

    [12]

    Staebler D L, Wronski C R 1977 Appl. Phys. Lett. 31 292Google Scholar

    [13]

    Zhang D, Tavakoliyaraki A, Wu Y, van Swaaij R, Zeman M 2011 Energy Procedia 8 207Google Scholar

    [14]

    Alnuaimi A, Islam K, Nayfeh A 2013 Sol. Energy 98 236Google Scholar

    [15]

    秦国刚, 孔光临 1988 半导体学报 01 103Google Scholar

    Qin G G, Kong G L 1988 J. Semicond. 01 103Google Scholar

    [16]

    Murty M V R, Atwater H A 1995 Phys. Rev. B 51 4889Google Scholar

    [17]

    Tersoff J 1988 Phys. Rev. B 37 6991Google Scholar

    [18]

    Macrae C F, Bruno I J, Chisholm J A, Edgington P R, Wood P A 2008 J. Appl. Crystallogr. 41 466Google Scholar

    [19]

    Hafner J, Kresse G 1997 Properties of Complex Inorganic Solids (New York: Plenum Press) pp69–99

    [20]

    Sriraman S, Agarwal S, Aydil E S, Maroudas D 2002 Nature 418 62Google Scholar

    [21]

    Robertson J 2000 J. Non-Cryst. Solids 266 79Google Scholar

    [22]

    Kim H, Horwitz J S, Kushto G, Piqué A, Kafafi Z H, Gilmore C M, Chrisey D B 2000 J. Appl. Phys. 88 6021Google Scholar

  • [1] 周明锦, 侯氢, 潘荣剑, 吴璐, 付宝勤. 锆铌合金的特殊准随机结构模型的分子动力学研究. 物理学报, 2021, 70(3): 033103. doi: 10.7498/aps.70.20201407
    [2] 王艳, 徐进良, 李文, 刘欢. 超临界Lennard-Jones流体结构特性分子动力学研究. 物理学报, 2020, 69(7): 070201. doi: 10.7498/aps.69.20191591
    [3] 李锐, 刘腾, 陈翔, 陈思聪, 符义红, 刘琳. 界面结构对Cu/Ni多层膜纳米压痕特性影响的分子动力学模拟. 物理学报, 2018, 67(19): 190202. doi: 10.7498/aps.67.20180958
    [4] 张忠强, 李冲, 刘汉伦, 葛道晗, 程广贵, 丁建宁. 石墨烯碳纳米管复合结构渗透特性的分子动力学研究. 物理学报, 2018, 67(5): 056102. doi: 10.7498/aps.67.20172424
    [5] 惠治鑫, 贺鹏飞, 戴瑛, 吴艾辉. 硅功能化石墨烯热导率的分子动力学模拟. 物理学报, 2014, 63(7): 074401. doi: 10.7498/aps.63.074401
    [6] 张程宾, 程启坤, 陈永平. 分形结构纳米复合材料热导率的分子动力学模拟研究. 物理学报, 2014, 63(23): 236601. doi: 10.7498/aps.63.236601
    [7] 周耐根, 胡秋发, 许文祥, 李克, 周浪. 硅熔化特性的分子动力学模拟–-不同势函数的对比研究. 物理学报, 2013, 62(14): 146401. doi: 10.7498/aps.62.146401
    [8] 曹万强, 舒明飞. 弛豫铁电体的键能与配位数模型. 物理学报, 2013, 62(1): 017701. doi: 10.7498/aps.62.017701
    [9] 肖红星, 龙冲生. UO2 晶体中低密勒指数晶面表面能的分子动力学模拟. 物理学报, 2013, 62(10): 103104. doi: 10.7498/aps.62.103104
    [10] 周化光, 林鑫, 王猛, 黄卫东. Cu固液界面能的分子动力学计算. 物理学报, 2013, 62(5): 056803. doi: 10.7498/aps.62.056803
    [11] 柯川, 赵成利, 苟富均, 赵勇. 分子动力学模拟H原子与Si的表面相互作用. 物理学报, 2013, 62(16): 165203. doi: 10.7498/aps.62.165203
    [12] 兰惠清, 徐藏. 掺硅类金刚石薄膜摩擦过程的分子动力学模拟. 物理学报, 2012, 61(13): 133101. doi: 10.7498/aps.61.133101
    [13] 马颖. 非晶态石英的变电荷分子动力学模拟. 物理学报, 2011, 60(2): 026101. doi: 10.7498/aps.60.026101
    [14] 马文, 祝文军, 陈开果, 经福谦. 晶界对纳米多晶铝中冲击波阵面结构影响的分子动力学研究. 物理学报, 2011, 60(1): 016107. doi: 10.7498/aps.60.016107
    [15] 刘建廷, 段海明. 不同势下铱团簇结构和熔化行为的分子动力学模拟. 物理学报, 2009, 58(7): 4826-4834. doi: 10.7498/aps.58.4826
    [16] 廖乃镘, 李 伟, 蒋亚东, 匡跃军, 祁康成, 李世彬, 吴志明. 椭偏透射法测量氢化非晶硅薄膜厚度和光学参数. 物理学报, 2008, 57(3): 1542-1547. doi: 10.7498/aps.57.1542
    [17] 周宗荣, 王 宇, 夏源明. γ-TiAl金属间化合物面缺陷能的分子动力学研究. 物理学报, 2007, 56(3): 1526-1531. doi: 10.7498/aps.56.1526
    [18] 张 林, 王绍青, 叶恒强. 大角度Cu晶界在升温、急冷条件下晶界结构的分子动力学研究. 物理学报, 2004, 53(8): 2497-2502. doi: 10.7498/aps.53.2497
    [19] 罗 志, 林璇英, 林舜辉, 余楚迎, 林揆训, 余云鹏, 谭伟锋. 氢化非晶硅薄膜中氢含量及键合模式的红外分析. 物理学报, 2003, 52(1): 169-174. doi: 10.7498/aps.52.169
    [20] 梁海弋, 王秀喜, 吴恒安, 王宇. 纳米多晶铜微观结构的分子动力学模拟. 物理学报, 2002, 51(10): 2308-2314. doi: 10.7498/aps.51.2308
计量
  • 文章访问数:  9998
  • PDF下载量:  144
出版历程
  • 收稿日期:  2019-08-22
  • 修回日期:  2020-02-06
  • 刊出日期:  2020-04-05

/

返回文章
返回