直拉单晶硅的制备 掺杂
直拉单晶硅的制备硅、锗等单晶制备，就是要实现由多晶到单晶的转变，即原子由液相的随机排列直接转变为有序阵列;由不对称结构转变为对称结构但这种转变不是整体效应，而是通过固液界面的移动而逐渐完成的为实 现上述转化过程，多晶硅就要经过由固态到熔融态，然后又由熔融态硅到固态晶体硅的转变这就是从熔 体硅中生长单晶硅所遵循的途径从熔体中生长硅单晶的方法，目前应用最广泛的主要有两种：有坩埚直 拉法和无坩埚悬浮区熔法在讨论这两种制备方法之前，还应讨论在制备单晶过程中必不可少的一些准备 工序包括掺杂剂的选择、坩埚的选择、籽晶的制备等，分别介绍如下：一、掺杂在制备硅、锗单晶时，通常要加入一定数量杂质元素（即掺杂）加入的杂质元素决定了被掺杂半导 体的导电类型、电阻率、少子寿命等电学性能掺杂元素的选择必须以掺杂过程方便为准，又能获得良好 的电学性能和良好晶体完整性为前提1掺杂元素的选择（1） 根据导电类型和电阻率的要求选择掺杂元素制备N型硅、锗单晶，必须选择V族元素（如P、As、Sb、Bi）；制备P型硅、锗单晶必须选择皿族元素（如B、Al、Ga、In、Ti）杂质元素在硅、锗晶体中含量的多少决定了硅、锗单晶的电阻率。
电阻率不仅与杂质浓度有关，而且与载流子的迁移率有关当 杂质浓度较大时，杂质对载流子的散射作用，可使载流子的迁移率大大降低，从而影响材料的导电能力 考虑到以上因素，从理论上计算了电阻率与杂质浓度的关系曲线，如图9-5所示在生产工艺上按电阻率的高低分档掺杂有三档：轻掺杂（适用于大功率整流级单晶）、中掺杂（适 用于晶体管级单晶）、重掺杂（适用于外延衬底级单晶）2） 根据杂质元素在硅、锗中溶解度选择掺杂元素各种杂质元素在硅、锗中溶解度相差颇大例如， 采用大溶解度的杂质，可以达到重掺杂的目的，又不会使杂质元素在晶体中析出影响晶体性能下表列出了常用掺杂元素在硅、锗单晶生长时掺入量的极限，超过了极限量，单晶生长不能进行3） 根据分凝系数选择掺杂元素前面已介绍了几种掺杂元素，在硅、锗熔体生长时的平衡分凝系数 的值可见，分凝系数小的掺杂元素，在熔体生长晶体时很难从熔体进入晶体因此要使K0小的杂质进 行重掺杂显然是不适宜的所以要根据分凝系数来选择掺杂元素是很重要的条件之一4） 根据杂质元素在硅、锗中扩散系数选择掺杂元素考虑到整个半导体器件的热稳定性和在半导体 器件制造工艺，特别是高温工艺如扩散、外延等工艺中，常常要求硅、锗单晶中掺杂元素的扩散系数要小 些好。
否则当高温扩散制作器件时，衬底的杂质也同时以反扩散方式进入外延层中，影响杂质的再分布， 对器件电性能不利在硅中，属于快扩散杂质元素的有：H、Li、Na、Cu、Fe、K、Au、He、Ag、S；属 于慢扩散杂质元素的有：Al、P、B、Ga、Ti、Sb、As、Bi等5） 根据杂质元素的蒸发常数选择掺杂元素由于单晶制备有时在真空条件下进行，故在晶体生长过 程中气液相之间处于非平衡态，杂质的挥发变得更为容易杂质元素的挥发往往会造成熔体中杂质量的变 化，这是掺杂过程中应加以考虑的因素在真空条件下杂质元素的挥发程度通常用蒸发触表示因为在 一定温度和真空条件下，单位时间内从熔体中蒸发出的杂质数量N与熔体表面积A、以及熔体中的杂质浓 度CL成正比，当浓度较低时，上述参量间的关系可用下式表示：N=EACl （9-10）式中比例常数E称蒸发常数因此可知，E的数值即单位时间从单位表面积上蒸发出去的杂质数因 此，E的数值大小可以表征杂质蒸发速率的快慢程度各种杂质元素在硅熔体中的蒸发常数E的数值由表 8给出掺杂元素的选择除上述四方面的因素外，还要考虑掺杂元素原子半径的影响杂质元素的原子半径与 硅、锗原子半径之差是影响晶体完整性的重要因素之一，应尽量选用与硅、锗原子半径相似的杂质元素作 为掺杂剂。
二、掺杂方式按掺杂剂的形式可分为元素掺杂和合金掺杂1) 元素掺杂即直接将纯杂质元素加入硅中它适于制备电阻率10-2~10-3Q«cm的重掺杂硅单晶2) 母合金掺杂 是将掺杂元素与硅先做成母合金(例如，硅锑合金、硅硼合金)，根据母合金含的 杂质量相应的加入母合金量3) 中子辐照掺杂(NTD)按上述两种掺杂方式掺杂，由于杂质分凝、蒸发、温度分布等因素的影 响使得杂质均匀性很难作好，尤其是高阻材料更困难近年来发展了一种中子嬗变掺杂技术可大大提高 硅单晶电阻率的均匀性中子嬗变掺杂的主要特征是：杂质不是从外部掺进原材料中去，而是被掺材料本身的一部分原子在辐 照过程中直接转变而来对硅材料而言，中子嬗变掺杂技术的基本原理简述如下：自然硅中有三种稳定的同位素，即285!、2%!、30Si，它们所占比例分别为92.28%、4.67%和3.05%这些同位素在Si中呈均匀分布将未掺杂的硅单晶材料放入核反应堆中辐照30 Si俘获热中子后嬗变成不稳定的同位素31 Si，并放出了丫射线；31 Si的半衰期为2.62小时，释放出B射线后变成稳定的磷(31P )这样就相当于在硅中引进了均匀分布的施主杂质 磷由于中子嬗变是在固相低温下进行的，而且中子在硅中射程为90〜100cm，所以它们的杂质可以参杂 很均匀、很准确。
被辐照的单晶应放在反应堆均匀中子流的范围内，在辐照完成后，须存放一定时间，经 测定放射性后才能投入安全使用辐射中，产生了大量的晶格缺陷，使得硅的电阻率高达105Q・cm，经 过900°C、3小时退火，电阻率即可达稳定值三、掺杂的计算掺入杂质的量，一般应考虑以下几个因素：①原材料中的杂质含量；②杂质的分凝效应；③如果是在 真空系统中拉制单晶，还应考虑杂质在真空中的蒸发效应；④另外还应考虑到在生长过程中坩埚或系统内 杂质的污染考虑了以上几个影响因素后才能正确计算该加入单晶内的杂质量然而，以上这些影响因素 的大小是随材料的生长工艺而变动的，必须针对具体问题具体分析下面对两种掺杂形式分别计算如下：(1) 掺杂杂质元素重量的计算在只有分凝效应影响时，对于直拉硅或锗单晶掺杂量的计算因直拉 晶体生长过程实际上是一个正常凝固过程，如果原材料较纯，则材料的电阻率P与杂质浓度CS有如下关 系(9-11)式中U为电子迁移率(对n型半导体)或空穴迁移率(对p型半导体)q为电子电荷根据正常凝固 的杂质分布表达式可知，对于拉出单晶的某一位置g的电阻率和熔体里的杂质浓度C0有如下关系：P= 1/q •日• k • C0(1 - g)-(1-k) (9-12)(9-12)式表示电阻率沿晶锭长度方向之变化。
假如我们所用的原材料是高纯硅或高纯锗，高纯原材料本身的杂质相对于掺入的杂质量可以忽略不 计；又假定杂质并不蒸发，在生长过程中也没有污染，这时(9-12)式中的C0就是m_ C M A _ 1 M A叭 0 d N 一 q—Kp(1 - g)-(1-k) d N ③⑶0 0式中d为硅或锗的密度，N0为阿伏伽德罗常数，A为杂质元素的原子量，A/N0为每个杂质原子的质量在计算时，需要有确定的g值（即相对位置），也就是说要知道所要求的电阻率的位置假定我们需 要生长某一电阻率的单晶，可以取X=1/2来计算；如果我们要获得某一电阻率范围例如1〜3Q・cm的单晶， 则可以取g=0 （头部）处的P=3Q・cm来计算下面举两个例子来具体说明例题I 要求拉出在g=1/2处，P=1Q・cm的锗单晶50g，所用原材料是区熔提纯的高纯锗，问需要 掺入杂质砷的质量为多少？（已知砷的分凝系数K=0.02,砷原子量为74, U用电子迁移率为4000cm2/v -s， N0=6.03X 1023/mol）解：用（9-12）和（9-13）式分别可求出C0和m为C0^5X1016cm-3m^0.057mg例题II若要求X=1/2处，P=1Q・cm的n型硅单晶50g，问需要掺入砷杂质多少？解：经计算得到C0^4X 1015cm-3m^0.01mg由该例可看出，掺入的杂质量是很少的，如此小的量在称重方面不易准确，会产生很大误差。
因此除 拉制重掺杂单晶外，一般将掺的杂质与原材料先做合金，再以合金的方式掺入，这样就可更准确的控制掺 杂量常用的母合金P-Si、B-Si、Ge-Sb、Ge-Ga等合金2）采用中间合金（母合金）掺杂时对母合金重量的计算①拉制P型硅或锗单晶时掺入母合金质量的计算：若掺入方式是掺入中间合金，比如硅硼合金的电阻 率为P］，对应的杂质浓度CH，若在硅中投入中间合金重为W］，因中间合金中的杂质含量相对硅的比例来 看一般是很小的可以把合金的密度看作和硅的密度相同，故可以硅的密度dds代替则投入的合金体积 为W1/ dds，那么投入多晶硅熔体中的中间合金的杂质原子总数为WJ d合）・CH设投入多晶熔硅原材料的 质量为M，那么熔硅体积为M/ dds，（忽略合金小体积），则掺入的杂质在熔体中的浓度CS导出如下根据杂质在合金中的总数与在硅熔体中总数不变得因为d合=dds，M+W^M (因W1很小可略)所以(9-14)(9-15)W _ C^M _1 C Kp q (1 - g)-(1-k) C在实际拉高阻单晶过程中，拉制出的单晶，其单晶头部的杂质浓度Cq（头）为:SWCS (头)=KCL=K M1CH±Cf式中第一项为投入的杂质在单晶头部凝固出来的杂质浓度，第二项Cf为多晶硅料（原材料）中和坩埚 污染引入的杂质对单晶头部杂质的净贡献。
在拉制P型单晶时，若多晶硅料含有杂质也是P型的，则它起 着补充P型杂质含量的作用此时（9-15）式的Cf应取正号；如果多晶硅中含有N型杂质，则它对投入 的P型杂质起补偿作用，则Cf应取负号因此，要拉制电阻率上限P上（对应CS（头））的P型单晶，应投入含杂质浓度为CH中间合金的重量W1是对P型多晶硅为WM [Cs - Cf]i K CH(9-16a)对N型多晶硅为WM [Cs (^) + Cf]1 KCH(9-16b)我们应用上式可计算投入中间合金的重量W］②拉制N型硅单晶时掺入母合金重量的计算：N型硅单晶通常是在真空条件下进行的，又因从中间合 金的投入到放肩生长完成，需要一段时间，在这段时间内熔硅中掺入的杂质元素必会有部分挥发，因此还 应考虑蒸发作用的影响显然杂质蒸发量与时间、熔体表面积、熔体中杂质浓度成正比因此，熔体中因 杂质挥发使杂质浓度的降低量dN可用下式表示：dN = - EAC ° - dt（9-17）式中E为蒸发常数，A^为熔体表面积，C0为熔体中杂质的初始浓度，dt为杂质浓度降低dN所需时间，C0与熔体中杂质数目的关系式为:(9-18)V为熔体体积，则c0在数值上为单位熔体中杂质数目。
将式（9-18）代入（9-17）得dN = - EAsN V上式积分后（9-19）所以EAInN In = 一所s t |tN 0 V 0N = N户((-EA^ / V)・ t]（9-20）当t=0时，即中间合金投入尚未开始挥发的情况下，此时熔体中的杂质数目N为:(9-21)即N等于中间合金中的杂质数熔体的体积可用熔体重M和密度d表示为:(9-22)将（9-21）、（9-22）式代入（9-20）式得(9-23)N = Kc exp[(-EA d/M)・t] = Kc exp[(-EA /或・t]d H s d H s式（9-23）描述了熔体中杂质浓度随时间推移而减小的变化规律9-24)前面已叙述过，已知开始结晶时单晶体头部的杂质浓度CS（P上）为：cs（P 上）=KC土 Cf式中CS（P上）表示单晶头部的电阻率最高，作为电阻率的上限P上所对应的杂质浓度CS（P上）式中右 边第一项为掺入杂质在单晶头部凝固出来的杂质浓度;第二项Cf为多晶硅及坩埚对单晶头部杂质浓度的净 贝献将（9-23）式代入（9-24）式得(9-25)(9-26)(9-27a)(9-27b)NCs(p 上)=K y 土 Cf也可写成:Cs(p 上)=K上式，即表示单晶头部杂质浓度CS（P上）随生长时间而变化的规律。
若多晶硅为N型，则Cf取正号，得W = M [ C(「U 一 * ]e -号1 K CH若多晶硅为P型，则Cf取负号，得W = M[ Cs(P 上)+ Cf ]e-伊1 K CH上式即为拉制单晶硅时掺入母合金重的计算式A0/V=1/2的熔体中的杂质浓度从CS降到CS /e所经历的时间记为t蒸发因为1/e=0.368，所以t蒸发表示 熔体中的杂质浓度经过时间t后蒸发掉63.2%剩下36.8%所需的时间，定义为t蒸发若知道杂质的蒸发常数 E就可算出t蒸发常见硅中杂质的t蒸发列于表9 寮由表9看出，一些杂质如Sb、As、In在硅中蒸发很快，只要几分钟就蒸发了大部分因此，掺有快 蒸发溶质的熔体，不宜在真空中生长掺杂单晶，但可在保护气氛下拉制3）直拉法拉制硅单晶的气氛有正压氩气、减压氩气（也有少数用氨气）和真空三种在不同的气 氛下，掺杂剂的蒸发情况不同掺杂剂算时，必须考虑它的影响下面我们分别讨论不同的气氛下的掺杂 计算A在正压或减压氩（氮）气气氛下拉晶的掺杂计算生长集成电路和分立器件用的大直径中、低阻晶体时，普遍采用减压氩气气氛只有部分直径较小的 晶体，还在正压氩气气氛下进行除重掺杂外，在两种气氛下拉晶，都不用纯元素，而用合金共熔法掺杂， 因为纯元素量太小且物理化学性质与硅熔体也相差太远。
在氮气气氛下拉晶的掺杂计算与氩气气氛下拉晶 相同在掺杂计算上忽略蒸发的影响CZ法生长是正常凝固过程，掺杂剂的轴向分布遵从Pfann关系式为简化计算，选取晶（1） 杂质的挥发；（2） 石英坩埚引入的杂质；（3） 多晶硅中的初始杂质浓度；（4）合金电阻率的不均匀度这些假设基本符合实际情况合金中的掺杂元素分凝进入晶体，达到与目标电阻率相应的掺杂剂浓度， 即有（a + b）Ct = kC b式中，a为原料多晶硅的重量，单位gb为掺入合金的重量，单位gCm为合金中的掺杂剂的浓度， 单位cm-3k为有效分凝系数Ct为晶体肩部大约锭长10%位置处目标电阻率对应的杂质浓度，单位cm-3 因而掺入合金重量由下式计算：C 1 a kC - Ct在低压氩气气氛下，拉晶速度为1mm/min时，在肩部为锭长10%位置上，几种常用掺杂元素在硅中 的有效分凝系数为：磷KPQ0.406；硼Kbq0.91；锑Ksb^0.052o按上式计算掺杂量，能够比较准确的控制晶体的电阻率在必要时，根据实际拉晶情况略加修正B当考虑石英坩埚沾污和多晶硅中的初始杂质浓度时，有下列计算公式：1当考虑多晶硅的初始杂质浓度C0时，有下列两种情况：（1）用p （或n）型多晶硅生长n （或p）型晶体C=心-kC b C+kC0 at a + b 0 kC - C - kC0(2) 用n (或p)型多晶硅生长n (或p)型晶体C _ k(C0a + C/)t a + bb £^S akC - Ct2当同时考虑石英坩埚沾污CB和多晶硅的初始杂质浓度C0时，有下列两种情况:（1）用p （或n）型多晶硅生长n （或p）型晶体b C + kC。
气 akC - C — ：C0 + 气④（2）用n （或p）型多晶硅生长n （或p）型晶体b Ct- kC0 一 CB akC - C + CbC _ k号 + Cmb) + Ct a+b B其中CB的符号，生长n型晶体取负号，生长p型晶体取正号3石英坩埚引入硅中的杂质浓度CB的计算熔硅与石英作用，使坩埚中的硼、铝、钙、氧及金属杂质引入硅中铝、钙、氧和金属或者不是浅能 级或者分凝系数很小因此，这里主要考虑B的沾污，其浓度CB由下式表示：式中，k方为B在Si中的有效分凝系数，rB为石英坩埚的沾污率，一般为（2〜10）X1010cm-2・S-1，SL为熔硅与石英坩埚的接触面积，cm2； VL为熔硅体积，cm3； t为开始熔化到晶体收肩所经历的时间，单 位So4在真空下拉晶的掺杂计算在真空下拉晶必须考虑掺杂剂的蒸发掺杂元素的蒸发系数以Sb、As最大，P、Al较小，B很小元 素的蒸发使其浓度发生变化可由下式表示：一— _ A 一dC = - EC —l dtL式中，E为掺杂元素的蒸发系数，Al为熔体蒸发面积，单位cm2, CL为晶体放肩时熔体中的掺杂元素 的浓度，单位cm-3其他符号意义同上对上式积分，得 C = C exp（ - EADt）L LO —L式中，CLO为熔体中初始杂质浓度，单位cm-3。
因CL0 =孚，故:C =~^exp （ —lt）L又由（a + b）C _ bCd^J =可式中，dL是熔体的密度，为2.5g/cm3, dS为晶体的密度，为2.33g/cm3, CS0为晶体中的平均杂质浓度， 单位cm-3由于dL^dS，Cm>>CL上式可简化为b = ―^a = — S0 a exp （- E^~l t）C m kCm 匕式中CS0的计算分下列四种情况1、 用p型多晶制备n型单晶时，C^ 0 = C + kC 0 + kCB2、 用n型多晶制备n型单晶时，C^0 = C - kC0 + kCB3、 用p型多晶制备p型单晶时，C^0 = Ct - kC0 - kCB4、 用n型多晶制备p型单晶时，Cs0 = C + kC0 - kCB5掺杂结果的修正在掺杂过程中，下列因素可能与掺杂计算不一致：（1）多晶硅的杂质测试不准确，（2）多晶中存在着 对导电性起作用的各种杂质，（3）多晶中有明显的杂质补偿，（4）母合金测算不准确，（5）石英坩埚中主 要杂质含量不一致，（6）蒸发系数和干果的沾污率选择不合适，（7）拉晶中从开始熔化到收肩的时间控制 不严格，（8）拉晶工艺的其他不稳定因素。
因此，拉出晶体的电阻率与掺杂计算有偏离，可在下一次拉晶中进行修正修正式可用:式中，场为上次拉晶时掺入的合金重量，P］为上次拉晶时晶体的电阻率，b2为下次拉晶时应掺入的 合金重量；p2为下次拉晶时晶体的电阻率四、单晶体中杂质浓度的均匀性（电阻率的均匀性）及其控制电阻率的均匀性是单晶制备过程中的一个重要问题半导体器件要求单晶电阻率很均匀，但在拉制单 晶过程中，影响电阻率均匀性因素很多例如：在拉晶过程中，生成的单晶中的杂质，是用熔体中得来的， 所以固态中杂质的分布均匀与否直接与熔体中的杂质分布有关一方面，由于液态中的杂质分布是否均匀 是由搅拌决定，因此固态中的电阻率均匀与否也和动态的搅拌有关另一方面，由于单晶生长过程同时存 在一个分凝过程，所以杂质的分布必然受到分凝作用的影响分凝作用使得在拉晶过程中，熔体中的杂质 浓度不断增加（或减少），因此拉出的晶体中的杂质浓度也不断变化又由于固态中杂质浓度cs与有效分 凝系数Ke,以及液态中杂质浓度CL（B）有关；而有效分凝系数Ke又与拉速f、转速w以及搅拌情况有关； 所以拉速、转速以及温度分布等因素都通过杂质分凝的影响而影响到电阻率的分布再有，在整个拉晶过 程中由于杂质蒸发也对电阻率的均匀性产生影响。
而且这些影响因素之间又有着紧密的联系为了便于叙 述起见，我们分纵向电阻率的均匀性和横向电阻率的均匀性两个问题分别加以说明1纵向电阻率均匀性的控制直拉单晶是个正常凝固的过程如果只考虑杂质分凝效应的影响，且K0<1时，在拉晶的过程中，杂 质不断由固态向液态集聚，这就使坩埚里的熔硅中的杂质浓度愈来愈大，因而单晶尾部的杂质浓度比头部 大得多正常凝固使晶体电阻率分布由大到小，其原因主要是由于拉晶过程中坩埚里的熔体不断减少；同时由 于杂质分凝作用，也会造成熔体内杂质浓度cl（b）的不断增加，从而造成晶体中杂质分布不均匀所以， 只要设法使坩埚中熔体的杂质浓度CL（B）和有效分凝系数Ke的乘积不变，即可使分凝出来的晶体中的 杂质浓度CS不变这样既可保证晶体电阻率均匀要达到纵向电阻率比较均匀的目的，下面介绍坩埚里 熔体杂质浓度CL（B）和Ke的乘积不变的几种方法：（1）改变拉速法从有效分凝系数的定义出发，有CS=KeCL（B）式中CS为固相中的杂质浓度；CL（B）为液相中的杂质浓度（只要熔体的体积远大于边界层的体积，CL⑻可以用熔体的平均浓度Cl来表示），Ke为有效分凝系数，它与f、以的关系可由下式表示:K 0 八 、 ,f・8 、K + (1 - K ) exp(- ^c )K0为平衡分凝系数，f为晶体生长速率，D为杂质在熔体中的扩散系数，5C为杂质积累层厚度（相当 于溶质边界层厚度），它与搅拌有关，其关系式为5 = 1.6 ID 1/3 门 1/63-1/2C由上述三个表示式看出，有效分凝系数与晶体生长速率（拉速）f有关。
当拉速非常慢时，在通常拉 晶工艺中拉速是很慢的，而且满足f • 5C/D<<1，因此Ke^K0这样就保证了 CL （B）^Cl不变和Ke^K0 的乘积不变从而保证了纵向电阻率的均匀性2） 悬浮坩埚法 悬浮坩埚法的装置如图9-6所示在大石墨坩埚中放一个小石墨坩埚，在小坩埚下 开一个连通孔（见图），当锗熔化后，借助熔锗的浮力和重力平衡原理，小坩埚就浮置于熔体中，这时将 杂质掺入小坩埚中，拉晶在小坩埚内进行来拉时小坩埚充满熔锗，且保持内外液面一样高，随着单晶的 长大，小坩埚内的熔锗逐渐减少，内液面随之下降，外面的锗熔液就要通过连通孔进入小坩埚，使之保持 内外液面一样高，这样就保证了小坩埚内锗熔体体积不变小坩埚中熔锗的杂质浓度C广M/V，M为熔锗 中的杂质总量，V为熔锗的体积因为一般选用的杂质的分凝系数比较小，如表中的铟的K值为0.001， 锑的K值为0.004因此对于分凝系数远小于1的杂质来说，分凝到单晶锗中的杂质量是非常少的另外， 杂质放在内坩埚中，外坩埚仅放纯锗熔体，杂质在内坩埚不易通过小孔流出，因而内坩埚中熔体锗中的杂 质总量M可看作是不变的因而CL就近似不变如果拉晶速率较慢，且拉晶条件不变，即K值也不变。
那么K与CL的乘积就基本不变因而晶体中的杂质浓度Cs=KCl在拉晶过程中保持不变这样就可以保 证纵向电阻率基本上是均匀的实践证明，用此法所拉制的锗单晶有70%以上部分的电阻率是均匀的而 一般直拉法拉出的单晶的均匀部分，最高也不过1/3，由此可见浮置坩埚法的优点3） 连续加料法这种方法和通常直拉单晶法基本相同，也是在单晶炉内进行，只是多加一个装置 这个装置的作用是使在拉出单晶的过程中连续不断地熔入与之等量的锗原料，以保持熔锗的量不变，则 CL基本上不变如果在拉晶过程中保证拉晶条件不变，那么纵向电阻率将是很均匀的该方法的设备如图 9-7所示有图看出，这个附加装置与拉晶杆连动使它以相同速度将锗锭通入坩埚中，最简单的装置，可 以用定滑轮改变速度方向的特点来达到，滑轮能人工的上下做适当调节当然这里有一个要求，那就是拉 出的单晶截面积要与熔入的锗锭的截面积相等，这样才能保证熔入量与拉出量相等，这需要使生长单晶处 与熔入锗锭处的温度分布情况尽可能接近，实际上应使拉晶处与熔处靠近，在拉出单晶上加冷却圈，在熔 入锗锭处加辅助加热器这种不断加料法拉出的单晶，经测试表明，其纵向电阻率分布的确是较一般正常凝固法拉出的单晶来 的均匀。
4） 双层坩埚法双层坩埚的结构是一个小坩埚外再套一个大埚其结构如图9-8所示所需的杂质 放在小埚内，在小埚内进行拉晶，随着晶体生长，小埚内的熔体慢慢减少，大埚内的熔体借助连通孔缓慢 流入小埚，使小埚内的熔体量基本上不变，杂质（对分凝系数小于1的杂质而言）浓度几乎不变，从而使 电阻率均匀双层坩埚的第二个优点是可以减少原材料的损耗因为所剩尾料还可以继续留在坩埚内使用第三个 优点是，操作方便四是对减少位错密度有利因此，该法是目前应用较广泛的一种方法在一定程度上， 该法克服了上述几种方法的缺点，保留了它们的优点上面介绍几种方法都未考虑蒸发效应，都属获得纵向均匀电阻率的直拉法下面介绍在直拉单晶过程中，同时考虑杂质分凝和蒸发效应是获得纵向均匀电阻率的条件，对于实际 工作是由指导意义的前面已经说过，如果在真空下拉晶，应考虑杂质的蒸发效应这时一方面分凝效应会使K<1的杂质富 集于熔体；另一方面蒸发效应将使杂质不断地从熔体逸出最后熔体内杂质将随时间而变化，其变化的程 度取决于这两个效应之和我们假定拉晶过程的情况如图9-9所示设欲拉总长为L的晶体，某时刻已拉出长为y的晶体，y/L=g 为在某时刻t已拉出的单晶占全部拉光熔体形成的单晶总长度L的分数。
设晶体截面积为AS，且在拉晶过 程中不变，熔体表面积为A也不变某时刻熔体体积V（g）可表示为V（g）= LAs（1 -g） （9-28）根据单位时间由蒸发引起熔体中杂质减少量应等于杂质的蒸发量，可写成：EAC =-V （g）（%） （9-29）L dt蒸发式中（V）蒸发表示单位时间由于蒸发引起熔体中杂质浓度的变化’可表示为:（虬）二（、如dt蒸发 dg）*） = （土）f dt dg L（9-30）f为拉晶速率，把（9-28）和（9-29）代入（9-30）式得（1 -幺）（导））发=-*C （9-31）dt 蒸发 fA LS另一方面，由于分凝作用，在单位时间内由于分凝作用引起熔体中杂质的变化量，应等于结晶时熔体 变为同体积晶体由分凝作用而进入熔体的杂质量（或者说熔体变为等量晶体后使熔体中杂质的改变量） 即：V（g）（dCL） =（1 - K）dyA Cdt 分凝 s L 没有电磁搅拌但由于熔体对流，起到搅拌作用另外，籽晶和坩埚反向旋转，也起到了搅拌作用在拉晶 过程中，由于杂质分凝，使得固液交界面附近熔体一侧有杂质积累（对K<1的杂质），且拉速愈快积累愈 多，若搅拌不均匀，则界面上杂质分布不均。
如果对熔体搅拌非常充分，使杂质刚在固液交界面上积累就 立刻被搅合均匀，那么拉出的单晶截面上的电阻率就趋于均匀了由此可见，造成横向电阻率不均匀的因 素之一就搅拌不均匀或搅拌不充分所致因此，为了控制横向均匀性，在拉晶过程中必须进行充分搅拌， 尤其是机械搅拌；另一个措施是适当降低拉速在拉晶过程中，处于固液交界面附近晶体径向各点的杂质 向交界面附近熔体一侧分凝出的杂质的速度（或反扩散的速度）是不等的从而造成了晶体径向杂质浓度 分布不均匀，随着拉速的降低，晶体径向各处的杂质有足够的时间向熔体进行扩散，会使晶体径向杂质分 布得到改善一般拉速的控制还与所需单晶直径的粗细有关，生长大直径单晶时，拉速过快，不但径向电 阻率均匀性变差，而且晶体直径的控制也很困难所以拉晶时适当降低拉速可使横向电阻率均匀性得到改 善或写成（1一 G（务）分凝=（1一 K）C（9-32）把（9-31）与（9-32）式相加，即得到蒸发与分凝影响同时存在时对晶体中杂质净贡献的方程式:（1 - g）（虬）+ （1 - g）（虬）=[1 - K - ^A]Cdt蒸发 dt分凝fAS或写成（1 - g）（E） = [1 - K - EA]Cdg蒸发 fA lS（9-33）（9-33）式为一次微分方程，解上述方程得:EAK + 1C = KC 0（1 - g）K+fAs'式中c0为开始时熔体中的杂质浓度，在有蒸发影响存在时蒸发常数愈大 之比A/A1愈大，拉速愈小（即拉晶时间愈长），则蒸发作用的影响愈大。
9-34）熔体表面积与单晶截面r EA如果K + TA - 1 = 0，即拉速为：f = — fAS （1- K）气S（9-35）CS=KC0也就是说，在这样的生长条件下（即上述的拉速下），直拉单晶纵向杂质分布基本上是均匀的在实 际生长中因为公式推导时假设的某些常量有出入，所以上述结果只能当做近似考虑2横向电阻率均匀性的控制横向电阻率分布不均匀，是指拉出的单晶横截面上杂质浓度分布不均匀一般来说，造成横向不均匀 性的原因比起纵向不均匀性的原因复杂，其影响因素较多不过概括起来主要有两个方面：其一是搅拌不 均；其二是固液交界面不平坦1）搅拌不均匀的影响已知直拉单晶过程是一个正常凝固过程，若熔体中杂质分布不均匀，必然使 凝固的单晶中的杂质分布也不均匀所以在拉晶过程中要进行充分搅拌，使熔体中杂质分布均匀是很重要 的要使固液界面熔体一侧杂质分布均匀，可采取两方面措施，一是进行充分搅拌，另一是适当降低拉速 一般用高频感应加热时产生涡流，同时其电磁搅拌作用，使杂质分布趋于均匀；在用电热炉加热时，虽然（2）界面不平坦的影响造成晶体横向电阻率不均匀的第二个原因是固液界面不是一个平面，如图 9-10所示固液交界面的形状和炉温分布，拉晶过程中由于加热系统和散热情况等因素的影响，使得晶体生 长面不是一个与生长方向垂直的平面，而是与等温面形状相同的面，且等电阻率面也是同样形状。
这是因 为分凝作用在同一固液界面上杂质浓度分布均匀但硅片不是沿固液界面切割，而是垂直生长方向切割 因此，同一平面上将包括多个固液界面造成的电阻率分布不均匀，即在同一水平上中心部分电阻率偏低， 四周高，如图（c）所示；若界面凹向熔体，则电阻率中心部分高，四周低关于生长界面的形状，与单 晶炉的热场分布有关而热场的分布与单晶炉系统石墨加热器、保温系统、支撑系统、石墨坩埚托、石墨 坩埚、拉晶杆等的传热情况有关所有这些热传导情况（包括加热与冷却）综合起来，究竟使晶体生长面 凸向熔体还是凹向熔体，对于不同的加热系统应作不同的具体分析一般情况下，将一根单晶沿纵向剖开， 用界面染色法可见到晶体径向结晶前沿的情况，如图9-11所示放肩处的结晶前沿呈上凸形，晶体的中部 结晶前沿呈平直形，而尾部则呈上凹形测其径向各点的电阻率会发现头部及尾部电阻率分布是不均匀的 而在晶体的中部结晶前沿平坦处，径向电阻率均匀性较好由此可见，平整的结晶前沿有利于形成均匀的 横向电阻率分布为了控制横向电阻率的均匀性，人们在实践中提出了如下一些措施：① 在拉晶过程中，搅拌要均匀，其简单办法是适当加大籽晶和坩埚转速，使固液交界面上杂质分布均 匀，对横向电阻率的均匀性有利。
② 要有适宜的温场使轴向温度梯度足够大，使径向温度梯度足够小因为只有较大的纵向温度梯度 （dT/dZ），才能及时传导出晶体生长时释放的潜热，维持晶体继续生长；而只有径向温度梯度（dT/dr）足够小， 才能尽量使固液交界平坦有利于改善电阻率的晶向均匀性和减小晶体内部热应力，提高晶体的完整性③ 采用偏心籽晶杆随着投料量的增大，所使用的坩埚直径也随之增加，使埚内熔体所受到的机械搅 拌作用相对减弱，造成熔体中的杂质分布不均匀因此，曾经在一些单晶炉的设计上，改变了过去传统方 法，将籽晶杆做成与坩埚轴不同心（并不意味着不平行），两轴杆可偏离数毫米，如图9-12所示这样就 增大了机械搅拌作用，使晶体径向电阻率的均匀性得到提高采用这种方法同时可降低晶体的位错但采 用这种方法时，应注意，由于坩埚轴与籽晶杆的不同心，拉晶时，不宜使晶体生长过粗，尤其在长到尾部 时更应如此具体偏离多少视不同的设备及所需的晶体直径而定这种偏心设备，目前已很少采用④ 采用双层坩埚前面已介绍过，该法不但对改善纵向电阻率均匀性有利，对改善横向电阻率的均匀 性也很有利。
