砷化镓晶体力学特性的各向异性计算
Anisotropic Calculation of Mechanical Property of GaAs Crystal
通讯作者: 姜 晨, 教授. E-mail:jc_bati@163.com
收稿日期: 2020-08-28 修回日期: 2020-10-15 网络出版日期: 2021-06-20
基金资助: |
|
Corresponding authors: JIANG Chen, professor. E-mail:jc_bati@163.com
Received: 2020-08-28 Revised: 2020-10-15 Online: 2021-06-20
Fund supported: |
|
砷化镓因其良好的光电特性被广泛应用于电子与半导体领域, 为推动砷化镓解理加工技术, 对砷化镓材料力学特性的各向异性进行计算并分析。本研究对砷化镓各个晶面之间的夹角、面间距、原子的密度等结构参数进行计算, 基于广义胡克定律结合压痕实验, 分析砷化镓材料表层弹性模量、泊松比、剪切模量、硬度、断裂韧性等力学特性在{100}晶面沿不同晶向力学性能的变化规律。结果表明: 砷化镓不同晶面间结构参数的不同是导致砷化镓力学特性呈现各向异性的主要原因; 砷化镓在{100}晶面上弹性模量、泊松比、剪切模量的各向异性均呈现出周期性变化, 且{100}晶面的剪切模量为恒值59.4 GPa; 砷化镓{100}晶面硬度的各向异性变化幅度较小, 断裂韧性变化幅度较大, 最小值为0.304 MPa·m1/2, 位于<110>晶向, 确定<110>晶向是裂纹最容易扩展的晶向。
关键词:
Gallium arsenide possesses excellent photoelectric property which have been widely used in electron and semiconductor industries. To promote the cleavage manufacturing technology of gallium arsenide, the anisotropy of mechanical property of gallium arsenide is calculated in this work. The structural parameters, such as angle between the crystal planes, the spacing between planes and density of atoms, were calculated. The change rule of elastic modulus, Poisson’s ratio, shear modulus, hardness and fracture toughness of monocrystalline gallium arsenide on crystal {100} plane with crystal orientation was analyzed by generalized Hooke’s law and indentation tests. The results show that difference of structural parameters between planes is the main reason for the anisotropy of mechanical properties of GaAs. With variation of crystal orientation, the elastic modulus and Poisson’s ratio on the {100} crystal plane appear periodic changes, while shear modulus is identically equal to 59.4 GPa. Anisotropy of GaAs {100} crystal plane hardness changes slightly while fracture toughness changes greatly with a minimum value of 0.304 MPa·m1/2 in <110> crystal direction which is sensitive to the crack propagation.
Keywords:
本文引用格式
董康佳, 姜晨, 任绍彬, 郎小虎, 高睿, 叶卉.
DONG Kangjia, JIANG Chen, REN Shaobin, LANG Xiaohu, GAO Rui, YE Hui.
目前, 国内外针对砷化镓材料力学性能及其结构特点展开了系列研究。Chen等[7]通过不同载荷下的压痕试验, 对砷化镓{100}晶面的<100>晶向和<110>晶向进行测试, 得到{100}晶面两个晶向的硬度、弹性模量、断裂韧性等相关参数; Xu等[8]使用波氏压头进行纳米压痕试验, 对砷化镓{100}晶面同一晶向施加不同载荷, 计算出不同载荷下的弹性模量、维氏硬度、断裂韧性, 得到下压载荷与力学参数的规律, 同时对砷化镓的损伤机理进行了说明; Ponraj等[9]使用不同压头在较高温度下对砷化镓进行压痕实验, 得到单一晶向的弹性模量、硬度、残余应力以及压头下压载荷与下压深度的关系, 同时采用电子显微镜对压痕形貌进行观察; Hjort等[10]通过对立方晶体结构的分析, 绘制出砷化镓的晶体结构, 对{111}晶面的特点进行详细介绍, 同时对弹性模量, 断裂韧性进行相关的计算, 得到单个晶面上弹性模量和断裂韧性沿着不同晶向的变化规律。
为推动砷化镓解理制造技术, 完整揭示砷化镓力学特性的各向异性, 本研究基于砷化镓的结构特点, 分析砷化镓力学特性各向异性产生的原因, 并且结合广义胡克定律和压痕试验, 计算出砷化镓弹性模量、泊松比、剪切模量、断裂韧性等力学特性, 得到砷化镓力学特性在{100}晶面沿不同方向上的变化规律, 为进一步提升砷化镓解理加工工艺技术提供指导依据。
1 砷化镓的结构特点
1.1 砷化镓的晶胞结构
砷化镓(GaAs)作为一种III-V族半导体材料, 属于闪锌矿晶格结构, 单个晶胞结构可以看作是两个面心立方晶格组成, 一个砷晶格和一个镓晶格。若以任一镓原子为空间坐标原点建立坐标系A, a为两相邻顶角原子中心距离, 其四个砷原子的空间坐标分别为(1/4a, 1/4a, 1/4a), (3/4a, 3/4a, 1/4a), (3/4a, 1/4a, 3/4a), (1/4a, 3/4a, 3/4a), 8个镓原子位于在面心立方晶格的8个顶角, 其余6个镓原子位于正六面体的6个面的中心处。砷化镓的各向异性主要表现在如图1所示的三个晶面, 其中镓原子为蓝色, 砷原子为红色。
图1
图1
砷化镓结构及主要晶面示意图
Fig. 1
Structure and main crystal plane of GaAs
(a) GaAs crystal structure; (b) GaAs (100) plane; (c) GaAs (110) plane; (d) GaAs (111) plane
1.2 砷化镓主要晶面的结构参数
由坐标系A可知, {100}晶面的面间距为a, {110}晶面的面间距为$\sqrt{2}/2a$, {111}晶面的面间距借助于面方程得到。相邻两个{111}晶面的面方程分别为式(1)和式(2):
进而得到{111}晶面的面间距为$\sqrt{2}/3a$。不同晶面之间的夹角是通过面方程计算出其法向量然后求得夹角。{100}晶面的面方程为式(3)和{110}晶面的面方程为式(4):
{100}晶面、{110}晶面、{111}晶面对应面的空间法向量分别为: (1, 0, 0)、(1, -1, 0)、(1, 1, 1), 由此求出{100}晶面与{110}晶面的夹角为45°, {100}晶面与{111}晶面的夹角为54.74°, {110}晶面与{111}晶面的夹角为90°。
砷化镓主要晶面原子分布如图2所示, 若以单位面积内原子的占有个数来描述其原子分布的密度, 则{110}晶面原子密度>{111}晶面原子密度>{100}晶面原子密度。正是由于各个晶面间的面间距、面夹角、原子密度的不同导致了砷化镓力学特性呈现出各向异性。
图2
图2
砷化镓各晶面原子分布结构图
Fig. 2
Atomic distribution structure of each crystal plane of GaAs
(a) GaAs (100) plane; (b) GaAs (110) plane; (c) GaAs (111) plane
2 砷化镓的力学性能
2.1 砷化镓的弹性模量
由广义胡克定律可知弹性范围内应力与应变之间呈现线性关系, 可以得出式(5)和式(6)[9], 其中σ为应力, ε为应变, C为刚度矩阵, S为C的可逆矩阵称为柔度矩阵。表达式如下:
因砷化镓为立方晶体结构, 具有高度的对称性, 其中${{S}_{11}}$、${{S}_{12}}$、${{S}_{44}}$为柔度系数。可以将上述式(6)转化为:
如图3所示, 在空间中建立直角坐标系B, 三个相互垂直的坐标轴分别为X, Y, Z, 假设l为砷化镓空间晶格内任意方向, ${{l}_{1}}$、${{l}_{2}}$、${{l}_{3}}$分别为<100>晶向、<010>晶向、<001>晶向与l方向上的余弦, 其中α是空间向量l与Z轴所形成的夹角, 而β是l在XY平面的投影与X轴的夹角, m为向量l终点且垂直于l向量的平面(此平面记作剪切面)上的任意向量, γ为向量m与向量l所在延长线在剪切面上投影所形成的夹角。利用柔度矩阵得到砷化镓空间晶格内任意方向上的弹性模量只与向量l有关, 可以表示为式(8)[11,12,13,14,15,16,17]:
图3
其中:
且$l_{1}^{2}+l_{2}^{2}+l_{3}^{2}=1$,$\ {{l}_{1}}$、${{l}_{2}}$、${{l}_{3}}\in [0,\ 1]$。根据表1中砷化镓的材料参数得到其空间单晶格的弹性模量以及{100}晶面的弹性模量。
Mechanical parameters of GaAs | S11/ (×10-12, Pa-1) | S12/ (×10-12, Pa-1) | S44/ (×10-12, Pa-1) |
---|---|---|---|
Value | 11.7 | 3.7 | 16.8 |
由上述公式可绘制出图4(a)砷化镓空间单晶格的弹性模量以及图4(b)砷化镓{100}晶面沿不同晶向的弹性模量, 进而得到砷化镓空间单晶格弹性模量最大值在<110>晶向处与面心立方晶格的8个顶角所对应, 为141.16 GPa, 最小值位于<100>晶向, 为85.26 GPa。砷化镓{100}晶面的弹性模量呈现周期性变化, 每90°一个周期, 一共4个周期, 若取$<\bar{1}10>$晶向到<110>晶向为一个周期, 其弹性模量的变化是沿着$<\bar{1}10>$晶向到<100>晶向递减, <100>晶向到<110>晶向递增, 且单个周期也存在对称性, 本周期沿着<100>晶向的弹性模量完全对称。
图4
图4
砷化镓弹性模量
Fig. 4
Elastic modulus of gallium arsenide
(a) Single lattice elastic modulus of GaAs; (b) Elastic modulus of GaAs (100) plane
2.2 砷化镓的泊松比
其中:
并且式(15)、式(16)、式(17)满足此式:$m_{1}^{2}+m_{2}^{2}+m_{3}^{2}=1$,${{m}_{1}}$、${{m}_{2}}$、${{m}_{3}}\in [0,\ 1]$。
由上述公式可绘制出图5(a)砷化镓空间单晶格的泊松比以及图5(b)砷化镓{100}晶面沿不同晶向的泊松比, 进而得到砷化镓空间单晶格泊松比最大值在<100>晶向处, 值为0.44, 最小值位于<110>晶向, 值为0.19。由图5(b)可知砷化镓{100}晶面的泊松比沿着不同晶向的值呈现出周期性变化, 共存在4个周期, 若取图中0°~270°为一个周期, 即沿着<110>晶向到$<\bar{1}10>$晶向的变化。单个周期内泊松比沿着<110>晶向到<100>晶向递增, 在<100>晶向上出现最大值, 然后沿着<100>晶向到$<\bar{1}10>$晶向递减, 在$<\bar{1}10>$晶向出现最小值, 且单个周期内沿着<100>晶向呈对称性。
图5
图5
砷化镓泊松比
Fig. 5
Poisson’s ratio of gallium arsenide
(a) Lattice Poisson’s ratio of GaAs; (b) Poisson’s ratio of GaAs (100) plane
2.3 砷化镓的剪切模量
图6
图6
砷化镓剪切模量
Fig. 6
Shear modulus of gallium arsenide
(a) Shear modulus of GaAs spatial lattice; (b) Shear modulus of GaAs (100) plane
3 砷化镓压痕试验
3.1 维氏压痕试验
压痕试验是通过压头垂直压入测试材料, 通过分析压头的下压载荷与在样品上形成的压痕特征, 从而计算得出材料的弹性模量、硬度、断裂韧性等相关力学特性。本次压痕试验, 采用HXS-1000Ak维氏压头显微硬度计, 试验所用压头为标准维氏金刚石压头, 压痕特征的观察也在该设备上进行。进行压痕实验时, 首先对{100}晶面为抛光面的砷化镓晶圆片预处理, 使用金刚石刀将砷化镓晶圆片按标定好的晶向进行切割, 方便压痕试验时确定材料的晶向。其中第一次取向为{100}晶面的<110>晶向, 之后逆时针每间隔15°进行一次取向, 共进行7次取向, 然后使用1.961 N的压力进行加载。
3.2 试验结果分析
砷化镓在相同压力下、不同晶向的压痕形状相同, 但是产生的裂纹长度、压痕的宽度不同, 由此可知不同晶向下的硬度和断裂韧性有差异(图7)。
图7
当试验结果满足c/b≥3时, 断裂韧性${{K}_{\text{IC}}}$的值较为准确。其中${{H}_{v}}$为维氏硬度, $F$为施加载荷, b为压痕半对角线长度、c为压痕选取角度对应晶向上的裂纹长度。根据式(20)和式(21)对砷化镓{100}晶面单个周期内的维氏硬度和断裂韧性进行求解。由图8可知: 砷化镓{100}晶面的维氏硬度波动范围较小, 在5.688~6.207 GPa之间, 均值为6.004 GPa, 在<100>晶向上存在最小值5.688 GPa, 在<110>晶向上存在最大值6.207 GPa。砷化镓{100}晶面的断裂韧性波动范围较大, 在0.304~0.602 MPa·m1/2, 均值为0.433 MPa·m1/2, <110>晶向到$<\bar{1}10>$晶向的变化规律为先增大后减小, 在<110>晶向存在最小值0.304 MPa·m1/2, 即{100}晶面上裂纹沿<110>晶向最易扩展。表2为已有相关研究和本次研究结果的汇总比较。
图8
图8
砷化镓维氏硬度和断裂韧性
Fig. 8
Vickers hardness and fracture toughness of gallium arsenide
Data source | Hardness/GPa | KIC/(MPa·m1/2) | Comments |
---|---|---|---|
Chen, et al[7] | 5.27-5.45 | 0.482-0.579 | <110> and <100> crystal direction measured by different loads |
Xu, et al[8] | 6-10 | - | Measured by different loads |
Ponraj, et al[9] | 8.2-10.4 | - | Measured by different loads |
Hjort, et al[10] | 7 | 0.44 | Average value |
This work | 5.69-6.21 | 0.30-0.60 | Measured by single loads and different crystal orientations |
4 结论
本研究基于砷化镓结构特点, 结合广义胡克定律和压痕试验, 研究了砷化镓的力学特性的各向异性, 结论如下:
1) {100}晶面的面间距为a, <110>晶面的面间距为$\sqrt{2}/2a$, <111>晶面的面间距为$\sqrt{3}/3a$, <110>晶面的原子密度大于{111}晶面的原子密度大于{100}晶面的原子密度; {100}晶面之间的距离大于{110}晶面的距离大于{111}晶面之间的距离; {100}晶面于{110}晶面的夹角为45°, {100}晶面与{111}晶面的夹角为54.74°, {110}晶面与{111}晶面的夹角为90°, 上述参数的不同导致了力学特性的各向异性。
2) 砷化镓空间弹性模量和泊松比以及剪切模量存在高度对称性, 弹性模量最大值为141.16 GPa, 最小值为85.26 GPa, 相差接近39.6%; 剪切模量的最大值为59.4 GPa, 最小值为32.5 GPa, 相差接近45.3%; 泊松比最大值为0.44, 最小值为0.19, 相差接近56.81%, {100}晶面砷化镓弹性模量和泊松比沿着不同晶向呈现4个周期分布, {100}晶面剪切模量沿着不同晶向没有变化, 恒值为59.4 GPa。
3) 砷化镓{100}晶面断裂韧性具有明显各向异性, 最大值与最小值相差49.54%, 在<110>晶向具有最小断裂韧性, 裂纹最容易扩展。
参考文献
Dicing of gallium arsenide high performance laser diodes for industrial applications: part I: scratching operation
Dicing of gallium- arsenide high performance laser diodes for industrial applications: Part II. cleavage operation
Cleavage fracture of brittle semiconductors from the nanometre to the centimetre scale
Estimation of energy savings when adopting ultrasonic vibration-assisted magnetic compound fluid polishing
Estimation of energy and time savings in optical glass manufacturing when using ultrasonic vibration- assisted grinding
Experimental investigation of influence of scratch features on GaAs cleavage plane during cleavage processing using a scratching capability index
Fundamental study of ductile- regime diamond turning of single crystal gallium arsenide
Mechanical behavior of undoped n-Type GaAs under the indentation of berkovich and Flat- Tip indenters
Studies of nanoindentation and residual stress analysis of Ge/GaAs epilayers
Gallium arsenide as a mechanical material
What is the Young's modulus of silicon?
On anisotropic elastic materials for which Young’s modulus E(n) is independent of n or the shear modulus G(n, m) is independent of n and m
The directional dependence of elastic stiffness and compliance shear coefficients and shear moduli in cubic materials
Calculated elastic constants for stress problems associated with semiconductor devices
The mechanical properties and elastic anisotropies of cubic Ni3Al from first principles calculations
Elastic properties of diamond-like phases based on carbon nanotubes
Anisotropy of elasticity of a Ni base single crystal superalloy
Indentation response of single-crystalline GaAs in the nano-, micro-, and macroregime
Analysis of anisotropic mechanical properties of monocrystalline silicon callow
/
〈 |
|
〉 |
