植成楊,甘志銀,潘建秋
華中科技大學微系統研究中心,湖北武漢 430074
GaN 材料具有良好的電學特性[1],如寬帶隙(3.39 eV)、高電子遷移率(室溫 1 000 cm 2 /V·s)、高擊穿電壓(3×106 V/cm)等,其優良的特性,誘人的應用前景和巨大的市場潛力,引來各國激烈的研究熱潮[2]。
氣相淀積生長法易于控制薄膜的厚度、組分和摻雜,是目前制備GaN 薄膜的主流方法,其中主要有金屬有機物化學氣相沉積(MOCVD)法和氫化物氣相外延(HVPE)法。
GaN 的MOCVD 法工藝是液態的TMGa 由載氣(H2 或N2)攜帶進腔體與另一路由載氣攜帶的NH3 在溫度1100℃左右的襯底上反應沉積得到[3]。MOCVD 系統的設計思想,通常要考慮系統密封性,流量、溫度控制要精確,組分變換要迅速,系統要緊湊等。
HVPE 技術具有設備簡單、成本低、生長速度快等優點,可以生長均勻、大尺寸GaN 厚膜,作為進一步用MOCVD 生長器件結構的襯底。HVPE 法中Ga 源先在800℃左右的溫區與HCl 反應生成GaCl,GaCl 再與NH3 混合在1 050℃左右的高溫區反應生成GaN[4-5],因此對應的設備需提供兩個溫區。
對于以上兩種主流的GaN 薄膜制備方法,本文提出一種在結構上可兼容兩種工藝的新式實驗型GaN 外延設備。在充分考慮了氣體密度、導熱系數、比熱容、粘性系數等參數與溫度的關系, 建立了二維和三維熱流耦合有限元模型,對襯底溫度、氣相溫度、流場分布等進行了系統分析。
氣相淀積生長設備構架存在共性,其都包含反應腔體、加熱系統、氣體輸運系統、冷卻系統等。需用理論和數據指導新設備腔體的設計,并合理搭建相匹配的加熱系統、氣體輸運系統、冷卻系統等。本文提出的新型設備如圖1 所示。加熱系統位于襯底基座下方,采用電阻片輻射加熱方式,側壁和底座為水冷結構。襯底上方布置可拆裝的陣列型鎵源支撐,并在襯底基座和鎵源支撐之間設置可拆裝的輔助加熱系統,以調節襯底基座及鎵源支撐的溫度。加熱片下方布置多道熱屏蔽板,提高加熱效率。進氣方式采用多氣孔進口方式,可以使反應腔內氣流均勻性得到保證。襯底基座采用石墨制造。

圖1 MOCVD 與HVPE 相兼容的腔體結構示意圖
在 comsol multiphysics 中構建三維溫度場有限元模型。該模型中設定反應腔體側壁與下底盤為水冷恒溫體。能量交換方式考慮加熱系統與襯底基座的輻射換熱, 襯底基座內部的熱傳導,輔助加熱系統與襯底基座及鎵源支撐的輻射換熱,鎵源支撐與襯底基座的輻射換熱,氣體和襯底基座的對流換熱,加熱片與熱屏蔽板的輻射換熱,冷卻水與腔體內壁及底座的對流換熱。各材料的熱物性參數取值于文獻[6]。
本文充分考慮了氣體密度、導熱系數、比熱容、粘性系數等參數與溫度的關系, 計算模型中采用的氫氣熱物性參數與溫度有如下關系(P=2.5kPa):

根據流體力學傳熱學理論, 反應氣體的狀態由以下方程描述:

Qs:熱源項;
Qr:輻射換熱項;其中面1 和面2 之間的輻射換熱可表示為下面的式子:

Xi,j:表示面i 到面j 的角系數;
ε :表面發射率;
Eb:黑體輻射力。
式 ( 5)為 連 續 性 方 程, (6)、(7)、(8)為 動 量 守 恒Navier-Stokes 方程, ( 9)為能量方程。式中 u、v、w 分別為氣體x、y、z 三分方向上的速度分量,C p 、 μT、T、λ、ρ 分別為氣體的比熱、粘度、溫度、導熱系數、密度。所有固體壁面均設置無滑移邊界條件即 u = v =w=0,氣體入口設置速度條件,出口邊界采用壓力邊界條件,電阻片施加功率邊界條件,反應室水冷壁面溫度為Tw = 340 K。構建的熱流耦合的有限元模型如圖2 所示。在該模型中, 氣體狀態為連續不可壓縮的層流氣體[7]。

圖2 三維熱流耦合有限元模型
本文提出的HVPE 結構部分是可以進行拆裝的,把這部分結構去掉,調整襯底基座與進氣口的相對位置,可以得到適合MOCVD 工藝的設備構架。
鑒于結構的對稱性,本文中的MOCVD 構架的有限元模擬分析可采用二維軸對稱模型,在保證模擬準確性的同時,減少網格數目,提高計算速率。其二維軸對稱模型如圖3 所示。圖中是有限元計算完的后處理溫度分布及流場分布圖。

圖3 二維軸對稱有限元模型溫度及流場分布圖
通過調整加熱電阻片的布置格局,及功率分配,可以使得襯底基座溫度在有效的半徑內達到非常好的均勻性。對于襯底基座外圍區域,低溫區域對基座的角系數占較大的比例,所以襯底基座外圍區域功率損耗較大。為保證襯底基座表面的溫度均勻性,在外圍區域設置獨立的加熱源,提供更大的加熱功率。

圖4 襯底基座徑向溫度分布圖
提取基座表面徑向的溫度分布,如圖4 所示,可以看出,在基座半徑50mm 范圍內,溫差在±1℃以內。在此半徑范圍內,可以放置3 片2 英寸的襯底。
另外,從反應腔內流線的分布可以看出,在襯底基座上方,氣體處于層流狀態,這將為外延層的生長提供了非常有利的穩定流場環境。
在上述二維軸對稱模型確定的加熱片格局后,確定了加熱系統布置結構。HVPE 法中Ga 源先在850℃左右的溫區與HCl反應生成GaCl,GaCl 再與NH3 混合在1 050℃左右的高溫區反應生成GaN。在原有MOCVD 結構基礎上,在襯底上方增加可拆裝的鎵源支撐,可組裝成適合HVPE 工藝的設備構架。在鎵源支撐于襯底之間可布置輔助加熱陣列,便于調整鎵源支撐所需的合適溫度和合理的安放高度。
考慮模型的對稱關系,選取完整模型的四分之一進行計算,在兩個對稱平面設置對稱邊界條件,以減少有限元模型的網格數目,提高運算效率。計算得到溫度分布云圖如圖5 所示。

圖5 三維有限元模型溫度及流場分布圖
通過合理布置鎵源支撐的位置,以及調節加熱功率,可使得其溫度在800℃~850℃之間。提取鎵源支撐所在高度截面可得到溫度分布曲線,如圖6 所示。

圖6 鎵源支撐溫度分布
圖中波峰近似于直線,是鎵源支撐體的溫度,其中間下凹部分是腔體中的氣體流動區域。
1)對于襯底基座的上表面,不同區域的功率損耗也不同。為保證表面溫度均勻性,可設置不同的加熱源,提供相匹配的加熱功率;
2)設計可拆裝的HVPE 結構,在共享MOCVD 構架的基礎上作出適當的調整,使得一臺設備適用于不同的工藝方法,并且在多物理場耦合有限元模型上得到了理論論證;
3)在進行有限元建模時,應注意建模策略,根據模型幾何特點,可簡化成二維軸對稱模型或者四分之一的幾何模型。如此,可在保證模型準確性的前提下,減少有限元模型的網格數目,提高運算效率。
[1]鄧志杰,鄭安生.半導體材料[M].北京:化學工業出版社,2004.
[2]李寶珠.寬禁帶半導體材料技術[J].電子工業專用設備,2010(187):05-10.
[3]孫玉芹.用于固態照明的非極性a面GaN薄膜的MOCVD生長及表征[D].華中科技大學,2011.
[4]Akinori K,Miho M,et al.Journal of Crystal Growth,246(3-4),(2002):230-236.
[5]Wei T B,Hu Q,Duan R F,et al.Journal of CrystalGrowth,2009,311(17):4153-4157.
[6]趙鎮南.傳熱學[M].北京:高等教育出版社,2002:492-509.
[7]李建國,劉實,李依依,等.熱絲化學氣相沉積金剛石薄膜空間場的數值分析[J].金屬學報,2005,41(4):437-443.