JPH05262593A - 単結晶の酸素濃度予測方法及び単結晶の引上方法 - Google Patents
単結晶の酸素濃度予測方法及び単結晶の引上方法Info
- Publication number
- JPH05262593A JPH05262593A JP34718792A JP34718792A JPH05262593A JP H05262593 A JPH05262593 A JP H05262593A JP 34718792 A JP34718792 A JP 34718792A JP 34718792 A JP34718792 A JP 34718792A JP H05262593 A JPH05262593 A JP H05262593A
- Authority
- JP
- Japan
- Prior art keywords
- oxygen concentration
- melt
- crystal
- pulling
- single crystal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Landscapes
- Crystals, And After-Treatments Of Crystals (AREA)
Abstract
(57)【要約】
【目的】結晶成長中に結晶中の酸素濃度を予測する。イ
ンゴット全長に亘る酸素濃度の予測を可能とし、評価の
時間短縮を図る。この酸素濃度予測方法を基に、酸素濃
度予測値を引上パラメータにフィードバックすることに
よって、結晶中酸素濃度の精度よい制御を行う。 【構成】CZ法による単結晶の引上げにおいて、引上げ
中のルツボ回転数、融液レベル、結晶固化率及び融液表
面温度を測定し、これらの値からモデル式によって引上
げ単結晶中の酸素濃度を求める。
ンゴット全長に亘る酸素濃度の予測を可能とし、評価の
時間短縮を図る。この酸素濃度予測方法を基に、酸素濃
度予測値を引上パラメータにフィードバックすることに
よって、結晶中酸素濃度の精度よい制御を行う。 【構成】CZ法による単結晶の引上げにおいて、引上げ
中のルツボ回転数、融液レベル、結晶固化率及び融液表
面温度を測定し、これらの値からモデル式によって引上
げ単結晶中の酸素濃度を求める。
Description
【0001】
【産業上の利用分野】本発明は、CZ法により製造され
る単結晶中の酸素濃度を予測する方法及びその予測値に
基づいて単結晶を引上げる方法に関する。
る単結晶中の酸素濃度を予測する方法及びその予測値に
基づいて単結晶を引上げる方法に関する。
【0002】
【従来の技術】従来、結晶の酸素濃度は次のような手順
で評価されていた。 (a)単結晶引き上げ終了後、単結晶を円筒状に加工す
る。 (b)加工が終わった単結晶を切断機に取付け、サンプ
ルを切断する。 (c)サンプルをエッチングする。 (d)サンプルを研磨する。 (e)サンプルを洗浄する。 (f)サンプルをFTIR(フーリエ変換型赤外分光
法)で測定する。
で評価されていた。 (a)単結晶引き上げ終了後、単結晶を円筒状に加工す
る。 (b)加工が終わった単結晶を切断機に取付け、サンプ
ルを切断する。 (c)サンプルをエッチングする。 (d)サンプルを研磨する。 (e)サンプルを洗浄する。 (f)サンプルをFTIR(フーリエ変換型赤外分光
法)で測定する。
【0003】また、シリコン結晶中酸素濃度の予測モデ
ルとしては、結晶中に取り込まれる酸素量を解析したモ
デルがある。(例えば、平田、干川:日本結晶成長学会
誌Vol.15No.2 1988、P.A.Rama
chandran J.Electrochem.Se
c.Vol.137,No.10,1990)。しか
し、実操業における種々の引上条件での結晶成長に対し
て、定量的に酸素濃度を予測した報告はない。
ルとしては、結晶中に取り込まれる酸素量を解析したモ
デルがある。(例えば、平田、干川:日本結晶成長学会
誌Vol.15No.2 1988、P.A.Rama
chandran J.Electrochem.Se
c.Vol.137,No.10,1990)。しか
し、実操業における種々の引上条件での結晶成長に対し
て、定量的に酸素濃度を予測した報告はない。
【0004】
【発明が解決しようとする課題】サンプルを取り出して
評価する方法では、ウェハーに加工してからの評価とな
るため、結果が得られるまでに時間がかかるという問題
がある。さらに、全数検査でなく抜き取り検査になるか
ら、規格はずれを見落とす危険性もある。従来の酸素濃
度予測モデルは、実操業における種々の引上条件での結
晶成長に対して、インゴットの軸方向酸素濃度分布を定
量的に予測する目的では、使用できなかった。この理由
は、モデル中で使用しているルツボ−融液界面濃度境膜
層厚(δx),雰囲気−融液界面濃度境膜層厚(δ
y),各界面温度(T1 ),(T2 )等のパラメータデ
ータと実際の引き上げ機で測定可能なルツボ回転数(C
R),メルトレベル(ML)(融液表面とヒーター上端
との差),メルト表面温度(T(Pyro))などのプロセス
データとの関係が分かっていなかったことによる。さら
に、従来文献では定常状態を仮定して解析しているの
で、時間変化を扱うことが出来なかった。つまり、実際
の引上においてはインゴット酸素濃度は軸方向で変化し
ているが、この酸素濃度変化を従来のモデルでは予測す
ることができなかった。
評価する方法では、ウェハーに加工してからの評価とな
るため、結果が得られるまでに時間がかかるという問題
がある。さらに、全数検査でなく抜き取り検査になるか
ら、規格はずれを見落とす危険性もある。従来の酸素濃
度予測モデルは、実操業における種々の引上条件での結
晶成長に対して、インゴットの軸方向酸素濃度分布を定
量的に予測する目的では、使用できなかった。この理由
は、モデル中で使用しているルツボ−融液界面濃度境膜
層厚(δx),雰囲気−融液界面濃度境膜層厚(δ
y),各界面温度(T1 ),(T2 )等のパラメータデ
ータと実際の引き上げ機で測定可能なルツボ回転数(C
R),メルトレベル(ML)(融液表面とヒーター上端
との差),メルト表面温度(T(Pyro))などのプロセス
データとの関係が分かっていなかったことによる。さら
に、従来文献では定常状態を仮定して解析しているの
で、時間変化を扱うことが出来なかった。つまり、実際
の引上においてはインゴット酸素濃度は軸方向で変化し
ているが、この酸素濃度変化を従来のモデルでは予測す
ることができなかった。
【0005】本発明の目的は、これらの問題点を解決し
た結晶中酸素濃度予測方法を提供することにある。また
本発明の別の目的は、この予測に基づいて引上げ結晶中
の酸素濃度を制御する方法を提供するものである。
た結晶中酸素濃度予測方法を提供することにある。また
本発明の別の目的は、この予測に基づいて引上げ結晶中
の酸素濃度を制御する方法を提供するものである。
【0006】
【課題を解決するための手段】本発明は、CZ法による
単結晶の引上げにおいて、引上げ中のルツボ回転数、融
液レベル、結晶固化率の操業条件と、融液表面温度との
関係を用いて、引上げ時に測定した実際の融液表面温度
から、引上げ単結晶中の酸素濃度を高精度に求めること
を特徴とする単結晶の酸素濃度予測方法である。
単結晶の引上げにおいて、引上げ中のルツボ回転数、融
液レベル、結晶固化率の操業条件と、融液表面温度との
関係を用いて、引上げ時に測定した実際の融液表面温度
から、引上げ単結晶中の酸素濃度を高精度に求めること
を特徴とする単結晶の酸素濃度予測方法である。
【0007】また、上記方法により求めた引上げ単結晶
中の酸素濃度により、ルツボ回転数及び融液レベルの操
業条件を結晶固化率に応じて調整し、同時にヒータ温度
を変化させて融液表面温度を調整しながら単結晶を引上
げることを特徴とする単結晶の引上方法である。これを
さらに具体的に説明すると、次のとおりである。すなわ
ち、CZ法によって石英ルツボ内のシリコン融液からシ
リコン単結晶を引き上げる際に、予めルツボ回転数(C
R),メルトレベル(ML),融液表面温度
(T(Pyro))と、濃度境膜層厚(δx),(δy),ル
ツボ−融液界面温度(T1 ),雰囲気−融液界面温度
(T2 )との関係を固化率(g)の関数の形で求めてお
く。実際の引き上げ中に測定可能なルツボ回転数(C
R),メルトレベル(ML),融液表面温度
(T(Pyro))を、非定常のまま解いた酸素濃度予測式に
代入し、インゴットの酸素濃度を予測する。さらに、そ
の予測を基にルツボ回転数(CR),メルトレベル(M
L),融液表面温度(T(Pyro))を制御することによっ
て、結晶中酸素濃度を制御する。
中の酸素濃度により、ルツボ回転数及び融液レベルの操
業条件を結晶固化率に応じて調整し、同時にヒータ温度
を変化させて融液表面温度を調整しながら単結晶を引上
げることを特徴とする単結晶の引上方法である。これを
さらに具体的に説明すると、次のとおりである。すなわ
ち、CZ法によって石英ルツボ内のシリコン融液からシ
リコン単結晶を引き上げる際に、予めルツボ回転数(C
R),メルトレベル(ML),融液表面温度
(T(Pyro))と、濃度境膜層厚(δx),(δy),ル
ツボ−融液界面温度(T1 ),雰囲気−融液界面温度
(T2 )との関係を固化率(g)の関数の形で求めてお
く。実際の引き上げ中に測定可能なルツボ回転数(C
R),メルトレベル(ML),融液表面温度
(T(Pyro))を、非定常のまま解いた酸素濃度予測式に
代入し、インゴットの酸素濃度を予測する。さらに、そ
の予測を基にルツボ回転数(CR),メルトレベル(M
L),融液表面温度(T(Pyro))を制御することによっ
て、結晶中酸素濃度を制御する。
【0008】
【作用】本発明によれば、実際の引き上げ機で測定可能
なルツボ回転数(CR),メルトレベル(ML),結晶
の固化率(g)などのプロセスデータを用い、さらに時
間変化(t)を扱うことができるモデルであるので、成
長中の結晶中酸素濃度(Cc)を予測することができ
る。
なルツボ回転数(CR),メルトレベル(ML),結晶
の固化率(g)などのプロセスデータを用い、さらに時
間変化(t)を扱うことができるモデルであるので、成
長中の結晶中酸素濃度(Cc)を予測することができ
る。
【0009】CZ法により製造されるシリコン結晶中の
酸素は、石英ルツボの溶解により融液に混入する。この
溶解反応はルツボ−融液界面での酸素の拡散速度が律束
であると言われている。融液に混入した酸素は大部分が
蒸発し、一部分がインゴット中に取り込まれる。ここ
で、石英ルツボの溶解量をX、雰囲気への蒸発量をY、
結晶へ取り込まれる量をZ、融液内の酸素量をQとする
と、結晶中への酸素混入機構を定量的に解析する上での
基本式は、次のようになる。
酸素は、石英ルツボの溶解により融液に混入する。この
溶解反応はルツボ−融液界面での酸素の拡散速度が律束
であると言われている。融液に混入した酸素は大部分が
蒸発し、一部分がインゴット中に取り込まれる。ここ
で、石英ルツボの溶解量をX、雰囲気への蒸発量をY、
結晶へ取り込まれる量をZ、融液内の酸素量をQとする
と、結晶中への酸素混入機構を定量的に解析する上での
基本式は、次のようになる。
【0010】 dQ/dt=X−Y−Z (1) 各界面における酸素の移動を考慮すると、融液内の酸素
量の保存式として次式が得られる。 dQ/dt=As・D・(Cs−Cm)/δx −Aa・D・(Cm−Ca)/δy −Ac・K・Vg・Cm (2) t:時間(sec) Q:融液中の酸素量(atoms) Ca:雰囲気−融液界面の酸素濃度(atoms/cm
3 ) Cs:ルツボ−融液界面の酸素濃度(atoms/cm
3 ) Cm:融液中の平均酸素濃度(atoms/cm3 ) K:偏析係数 D:酸素の拡散係数(cm2 /sec) As:融液−ルツボ界面積(cm2 ) Aa:融液−雰囲気界面積 Ac:融液−結晶界面積(cm2 ) Vg:結晶の成長速度(cm/sec) 結晶中の酸素濃度は融液中の酸素濃度を用いて Cc=K・Cm (3) ただし、Cc:結晶中酸素濃度(atoms/cm3 )
と表される。
量の保存式として次式が得られる。 dQ/dt=As・D・(Cs−Cm)/δx −Aa・D・(Cm−Ca)/δy −Ac・K・Vg・Cm (2) t:時間(sec) Q:融液中の酸素量(atoms) Ca:雰囲気−融液界面の酸素濃度(atoms/cm
3 ) Cs:ルツボ−融液界面の酸素濃度(atoms/cm
3 ) Cm:融液中の平均酸素濃度(atoms/cm3 ) K:偏析係数 D:酸素の拡散係数(cm2 /sec) As:融液−ルツボ界面積(cm2 ) Aa:融液−雰囲気界面積 Ac:融液−結晶界面積(cm2 ) Vg:結晶の成長速度(cm/sec) 結晶中の酸素濃度は融液中の酸素濃度を用いて Cc=K・Cm (3) ただし、Cc:結晶中酸素濃度(atoms/cm3 )
と表される。
【0011】ここまでは、従来の酸素濃度予測モデルで
ある。ここで本発明では、各々の量に対して時間変化お
よび、結晶の固化率に対する変化を考える。上式で融液
内の酸素量(Q)は引き上げた結晶の固化率(g)を用
いて、次のように表すことができる。なお、式中の
(t)は時間変化を考慮に入れるという意味である。
ある。ここで本発明では、各々の量に対して時間変化お
よび、結晶の固化率に対する変化を考える。上式で融液
内の酸素量(Q)は引き上げた結晶の固化率(g)を用
いて、次のように表すことができる。なお、式中の
(t)は時間変化を考慮に入れるという意味である。
【0012】 Q(t)=Cm(t)・Mo ・[1−g(t)] (4) Mo :総融液量(cm3 ) g:固化率(×100%) さらに固化率g(t)は次のようになる。
【0013】
【数1】
【0014】そして、本酸素濃度予測法では、ルツボ−
融液界面の酸素濃度として、シリコンの酸素溶解度の温
度依存性を表わす式を用いた。(平田、干川(NTT
LSI研究所)応用物理学会公演予稿集CZ法シリコン
融液中の酸素溶解度とその温度依存性)。これにより、
引き上げ中の変化を取り込んだ。これは、次式のようで
ある。
融液界面の酸素濃度として、シリコンの酸素溶解度の温
度依存性を表わす式を用いた。(平田、干川(NTT
LSI研究所)応用物理学会公演予稿集CZ法シリコン
融液中の酸素溶解度とその温度依存性)。これにより、
引き上げ中の変化を取り込んだ。これは、次式のようで
ある。
【0015】 Cs(t)=4.0×1023exp(−2.0×104 /T1 (t)) (6) T1 :ルツボ−融液界面温度(K) さらに本酸素濃度予測法では、雰囲気−融液界面の酸素
濃度として、雰囲気界面温度(T2 (t))とSiOの
分圧(P(t))により界面濃度を表す式を用いた。こ
れにより、引き上げ中の変化を取り込んだ。これは、次
式のようである。
濃度として、雰囲気界面温度(T2 (t))とSiOの
分圧(P(t))により界面濃度を表す式を用いた。こ
れにより、引き上げ中の変化を取り込んだ。これは、次
式のようである。
【0016】 Ca(t)=5.0×1020P(t)exp(17290 /T2 (t)−11.1) (7) P:SiOの分圧 T2 :雰囲気−融液界面温度(K) 上記(7)式は、T.Carlberg:J.Elec
trochem.Soe.”SOLID−STATES
CIENCE AND TECHNOLOGY”Vo
l.133.No.9(1986)による。
trochem.Soe.”SOLID−STATES
CIENCE AND TECHNOLOGY”Vo
l.133.No.9(1986)による。
【0017】基本式(2)、(3)を非定常のまま扱
い、さらに時間の関数でなく、固化率の関数として表し
たものが以下の式である。ここで、各界面の酸素濃度C
s(g),Ca(g)は式(6)、(7)を用い温度の
関数として表わされる。
い、さらに時間の関数でなく、固化率の関数として表し
たものが以下の式である。ここで、各界面の酸素濃度C
s(g),Ca(g)は式(6)、(7)を用い温度の
関数として表わされる。
【0018】
【数2】
【0019】CM :総融液重量(kg) dl :融液の密度(kg/cm3 ) ds :結晶の密度(kg/cm3 ) a,b:次式で定義される定数 Vg=aL+b Vg:結晶の成長速度(cm/sec) L:結晶の長さ(cm) ここで、ルツボ−融液界面の温度T1 (t)、雰囲気−
融液界面の温度T2 (t)は、標準条件での引上げにお
ける融液内温度分布の測定から、例えば下の関数で決定
する。なお、この関数形は、このような簡単な形でも、
使用可能なことを示すものであり、式(8)(9)のよ
うに、固化率(g)・ルツボ回転数(CR)・メルトレ
ベル(ML)等引き上げ条件によりルツボ内融液温度分
布を表わすものであれば、精度を増すために式を複雑に
することも可能である。
融液界面の温度T2 (t)は、標準条件での引上げにお
ける融液内温度分布の測定から、例えば下の関数で決定
する。なお、この関数形は、このような簡単な形でも、
使用可能なことを示すものであり、式(8)(9)のよ
うに、固化率(g)・ルツボ回転数(CR)・メルトレ
ベル(ML)等引き上げ条件によりルツボ内融液温度分
布を表わすものであれば、精度を増すために式を複雑に
することも可能である。
【0020】 T1 =(a1 +a2 g)+(a3 +a4 CR) +(a5 +a6 ML)+T(PYRO) (8) T2 =(b1 +b2 g)+(b3 +b4 CR) +(b5 +b6 ML)+T(PYRO) (9) また、上記Cc(g)の式を用いて結晶中の酸素濃度を
計算した結果、結晶中酸素濃度は融液−ルツボ界面と融
液−雰囲気界面の濃度境膜層厚の比(δx/δy)に依
存することがわかった。この計算の際、各界面における
濃度境膜層厚の比はフィッティングパラメータとする。
このフィッティングパラメータδx/δyは他の条件が
同一であればルツボ回転数により一義的に決まる関数に
より決定され、さらに、固化率に対しても一義的に決ま
る関数により決定できることを見出した。
計算した結果、結晶中酸素濃度は融液−ルツボ界面と融
液−雰囲気界面の濃度境膜層厚の比(δx/δy)に依
存することがわかった。この計算の際、各界面における
濃度境膜層厚の比はフィッティングパラメータとする。
このフィッティングパラメータδx/δyは他の条件が
同一であればルツボ回転数により一義的に決まる関数に
より決定され、さらに、固化率に対しても一義的に決ま
る関数により決定できることを見出した。
【0021】拡散層厚δx,δyは、以下のように求め
る。融液ルツボ界面の濃度境膜層厚は無限円盤下の濃度
境膜層厚を与える実験式によりつぎのように与えられ
る。 δx=1.6×D1/3 ・(4×10-3)1/6 ×ω-1/2 (10) ここで、 D:拡散定数 ω:無限円盤の回転数(rps) である。ただし、CZ法によるシリコン単結晶成長の場
合、回転するのは無限円盤でなく、有限の範囲にあるル
ツボである。そこで、実際に結晶成長を行い、その結晶
の評価結果から、無限円盤の回転数(ω)と有限範囲ル
ツボの回転数(CR)の対応関係を導いた結果、例えば
一次式などの関数で決定されることを見出した。c1 〜
c4 をCZ炉により異なる係数として、下式のように表
わすことができる。
る。融液ルツボ界面の濃度境膜層厚は無限円盤下の濃度
境膜層厚を与える実験式によりつぎのように与えられ
る。 δx=1.6×D1/3 ・(4×10-3)1/6 ×ω-1/2 (10) ここで、 D:拡散定数 ω:無限円盤の回転数(rps) である。ただし、CZ法によるシリコン単結晶成長の場
合、回転するのは無限円盤でなく、有限の範囲にあるル
ツボである。そこで、実際に結晶成長を行い、その結晶
の評価結果から、無限円盤の回転数(ω)と有限範囲ル
ツボの回転数(CR)の対応関係を導いた結果、例えば
一次式などの関数で決定されることを見出した。c1 〜
c4 をCZ炉により異なる係数として、下式のように表
わすことができる。
【0022】 ω=(c1 +c2 CR)+(c3 +c4 g) (11) よって、融液−ルツボ界面の濃度境膜層厚δxは下式の
ように、ルツボ回転数と固化率の一次式の−1/2乗の
関数により決定される。 δx=1.6×D1/3 ・(4×10-3)1/6 ×{(c1 +c2 CR)+(c3 +c4 g)}-1/2 (12) 図1は(12)式から決定されるδxとCRとの関係の
一例を示したものである。この図よりルツボ回転数が低
い領域を除けば、(12)式の関係を図2のように直線
で近似しても誤差は少ないことがわかる。従って、実用
的な引上げ条件の範囲では、融液ルツボ界面の濃度境膜
層厚δxは下式のように、ルツボ回転数と固化率の一次
式により決定することも可能である。
ように、ルツボ回転数と固化率の一次式の−1/2乗の
関数により決定される。 δx=1.6×D1/3 ・(4×10-3)1/6 ×{(c1 +c2 CR)+(c3 +c4 g)}-1/2 (12) 図1は(12)式から決定されるδxとCRとの関係の
一例を示したものである。この図よりルツボ回転数が低
い領域を除けば、(12)式の関係を図2のように直線
で近似しても誤差は少ないことがわかる。従って、実用
的な引上げ条件の範囲では、融液ルツボ界面の濃度境膜
層厚δxは下式のように、ルツボ回転数と固化率の一次
式により決定することも可能である。
【0023】 δx=(c5 +c6 CR)+(c7 +c8 g) (13) 融液−雰囲気界面の濃度境膜層厚さδyの変化はルツボ
回転数(CR)・固化率(g)の変化に対して変化が小
さいため本方法では一定と仮定したが、厚さの変化を固
化率やルツボ回転数・結晶回転数・メルトレベル・雰囲
気のSiOの分圧・等による関数として決定することも
可能である。
回転数(CR)・固化率(g)の変化に対して変化が小
さいため本方法では一定と仮定したが、厚さの変化を固
化率やルツボ回転数・結晶回転数・メルトレベル・雰囲
気のSiOの分圧・等による関数として決定することも
可能である。
【0024】結晶中酸素濃度は融液−ルツボ界面と融液
−雰囲気界面の濃度境膜層厚の比(δx/δy)に依存
するため、融液−雰囲気界面の濃度境膜層厚δyを一定
値にしてもフィッテイング可能であり、酸素濃度予測・
制御可能である。 δx:融液−ルツボ界面の濃度境膜層厚(cm) δy:融液−雰囲気界面の濃度境膜層厚(cm) δx/δy:融液−坩堝界面と融液−雰囲気界面の濃度
境膜層厚の比 g:固化率(×100%) CR:坩堝回転数(rpm) ML:メルトレベル(mm)(ヒータと融液表面の相対
的高さ) T(Pyro):放射温度計測定による融液表面温度(K) D:シリコン融液中での酸素の拡散定数 係数a1 〜a6 、b1 〜b6 はもちろんのこと係数c1
〜c6 もCZ炉によって異なるので、後記実施例に示す
ように、標準引き上げの結果から求める。
−雰囲気界面の濃度境膜層厚の比(δx/δy)に依存
するため、融液−雰囲気界面の濃度境膜層厚δyを一定
値にしてもフィッテイング可能であり、酸素濃度予測・
制御可能である。 δx:融液−ルツボ界面の濃度境膜層厚(cm) δy:融液−雰囲気界面の濃度境膜層厚(cm) δx/δy:融液−坩堝界面と融液−雰囲気界面の濃度
境膜層厚の比 g:固化率(×100%) CR:坩堝回転数(rpm) ML:メルトレベル(mm)(ヒータと融液表面の相対
的高さ) T(Pyro):放射温度計測定による融液表面温度(K) D:シリコン融液中での酸素の拡散定数 係数a1 〜a6 、b1 〜b6 はもちろんのこと係数c1
〜c6 もCZ炉によって異なるので、後記実施例に示す
ように、標準引き上げの結果から求める。
【0025】なお、融液−結晶界面積Acは結晶成長中
の実測値を用いる。融液−ルツボ,融液−雰囲気界面の
接触面積As,Aaは使用するルツボの形に対応させ
て、各々固化率の関数として表わし、本酸素濃度予測式
に用いる。またシリコン融液中での酸素の拡散定数は種
々の値が報告されているが、本酸素濃度予測法では例え
ば次の値を用いる。
の実測値を用いる。融液−ルツボ,融液−雰囲気界面の
接触面積As,Aaは使用するルツボの形に対応させ
て、各々固化率の関数として表わし、本酸素濃度予測式
に用いる。またシリコン融液中での酸素の拡散定数は種
々の値が報告されているが、本酸素濃度予測法では例え
ば次の値を用いる。
【0026】 D=5.0×10-4cm2 /sec (14) 以上のように、各界面の温度式(8)、(9)及び拡散
層厚比の推定式(12)または(13)を用いて、T
1 ,T2 ,δxを求め(δyは一定値を用いる)、これ
を酸素濃度予測式に代入することにより、結晶成長条件
と融液表面温度の変化から結晶中酸素濃度のインゴット
軸方向分布(固化率変化)を予測することができる。
層厚比の推定式(12)または(13)を用いて、T
1 ,T2 ,δxを求め(δyは一定値を用いる)、これ
を酸素濃度予測式に代入することにより、結晶成長条件
と融液表面温度の変化から結晶中酸素濃度のインゴット
軸方向分布(固化率変化)を予測することができる。
【0027】
【実施例】以下、本発明を具体的に説明する。式(8)
〜(13)での係数決定のための標準引き上げ条件とし
て、例えば、チャージ量45kg、ルツボ径16イン
チ、結晶径6インチに対して、ルツボ回転数4rpm、
14rpm、メルトレベル(融液表面とヒーター上端の
差)100mm、30mmの4条件での引き上げ結果か
ら、フィッティングにより各係数を決定した。今回フィ
ッティングにより決定した各パラメータは次の通りであ
る。
〜(13)での係数決定のための標準引き上げ条件とし
て、例えば、チャージ量45kg、ルツボ径16イン
チ、結晶径6インチに対して、ルツボ回転数4rpm、
14rpm、メルトレベル(融液表面とヒーター上端の
差)100mm、30mmの4条件での引き上げ結果か
ら、フィッティングにより各係数を決定した。今回フィ
ッティングにより決定した各パラメータは次の通りであ
る。
【0028】結晶中酸素濃度は融液−ルツボ界面と融液
−雰囲気界面の濃度境膜層厚の比(δx/δy)に依存
するので融液−雰囲気界面の濃度境層厚δyはある一定
値でかまわない。本実施例では例えば次の値を用いる。 δy=0.01(通常のδyのオーダー) この時 δx=−0.0117g−0.0016CR+0.0461 T1 =T(Pyro)+3.00g+1.10CR−0.33ML+31.0 T2 =T(Pyro)+5.00g+0.80CR−0.12ML+3.0 この係数の値は、本モデルにより酸素濃度を予測制御す
ることを繰り返すとき、フィッティングを続けること
で、精度を増すことが可能となる。
−雰囲気界面の濃度境膜層厚の比(δx/δy)に依存
するので融液−雰囲気界面の濃度境層厚δyはある一定
値でかまわない。本実施例では例えば次の値を用いる。 δy=0.01(通常のδyのオーダー) この時 δx=−0.0117g−0.0016CR+0.0461 T1 =T(Pyro)+3.00g+1.10CR−0.33ML+31.0 T2 =T(Pyro)+5.00g+0.80CR−0.12ML+3.0 この係数の値は、本モデルにより酸素濃度を予測制御す
ることを繰り返すとき、フィッティングを続けること
で、精度を増すことが可能となる。
【0029】この決定された係数をもとに、結晶中酸素
濃度を予測した例を以下に示す。本実施例では、CZ炉
にて16インチルツボで6インチ結晶をメルトレベル7
0mm,ルツボ回転数6rpm及び12rpmの2条件
で結晶成長させた。また、ルツボ回転数6rpmで結晶
成長させたときの、放射温度計で測定した融液表面温度
は図3である。ルツボ回転数12rpmで結晶成長させ
たときの融液表面温度測定値は図示しなかったが図3に
示したのと同様に測定を行い、本酸素濃度予測法に用い
た。これらの条件により、引き上げた結晶の酸素濃度を
FTIRによって評価したものと、本酸素濃度予測法で
計算したものと比較したものが図4である。このよう
に、ルツボ回転数6rpm,12rpmいずれの場合も
結晶の酸素濃度を精度よく予測できていることが分る。
濃度を予測した例を以下に示す。本実施例では、CZ炉
にて16インチルツボで6インチ結晶をメルトレベル7
0mm,ルツボ回転数6rpm及び12rpmの2条件
で結晶成長させた。また、ルツボ回転数6rpmで結晶
成長させたときの、放射温度計で測定した融液表面温度
は図3である。ルツボ回転数12rpmで結晶成長させ
たときの融液表面温度測定値は図示しなかったが図3に
示したのと同様に測定を行い、本酸素濃度予測法に用い
た。これらの条件により、引き上げた結晶の酸素濃度を
FTIRによって評価したものと、本酸素濃度予測法で
計算したものと比較したものが図4である。このよう
に、ルツボ回転数6rpm,12rpmいずれの場合も
結晶の酸素濃度を精度よく予測できていることが分る。
【0030】また、以上により結晶中酸素濃度の予測が
可能となったことから、この予測をもとに、(ルツボ回
転数(CR),メルトレベル(ML),融液表面温度
(T(P yro))などの引上パラメータを制御することによ
って、結晶中酸素濃度の制御も可能となった。図6は結
晶成長中にリアルタイムで酸素濃度を予測・制御する方
法の、概念図を表すものである。この方法を基に、結晶
中酸素濃度を一定に制御したものが図5である。本図の
ように、酸素濃度を目標値の±0.5ppma以内に制
御することが可能となった。
可能となったことから、この予測をもとに、(ルツボ回
転数(CR),メルトレベル(ML),融液表面温度
(T(P yro))などの引上パラメータを制御することによ
って、結晶中酸素濃度の制御も可能となった。図6は結
晶成長中にリアルタイムで酸素濃度を予測・制御する方
法の、概念図を表すものである。この方法を基に、結晶
中酸素濃度を一定に制御したものが図5である。本図の
ように、酸素濃度を目標値の±0.5ppma以内に制
御することが可能となった。
【0031】
【発明の効果】本発明により、結晶成長時に成長を行な
っているその場で、結晶中酸素濃度を予測することが可
能となり、初めて行う引き上げ条件でも酸素濃度を予測
することが可能となった。さらに、引き上げ中の結晶変
形抑制等のためのヒーター温度変化に対する酸素濃度変
化も予測可能となった。したがって、 (a)インゴット全長に亘る酸素濃度の予測が可能とな
り、スペック外れの危険が減った。 (b)結晶成長中に酸素濃度予測が可能となり評価の時
間短縮が出来る。
っているその場で、結晶中酸素濃度を予測することが可
能となり、初めて行う引き上げ条件でも酸素濃度を予測
することが可能となった。さらに、引き上げ中の結晶変
形抑制等のためのヒーター温度変化に対する酸素濃度変
化も予測可能となった。したがって、 (a)インゴット全長に亘る酸素濃度の予測が可能とな
り、スペック外れの危険が減った。 (b)結晶成長中に酸素濃度予測が可能となり評価の時
間短縮が出来る。
【0032】この酸素濃度予測法を基に、酸素濃度予測
値を引上パラメータにフィードバックすることによっ
て、結晶中酸素濃度の精度よい制御が可能となった。
値を引上パラメータにフィードバックすることによっ
て、結晶中酸素濃度の精度よい制御が可能となった。
【図1】δxとCRとの関係の一例を示すグラフであ
る。
る。
【図2】図1のグラフを直線で近似したグラフである。
【図3】融液表面温度と固化率との関係を示すグラフで
ある。
ある。
【図4】各条件での酸素濃度(FTIR)と本発明での
計算値の比較を示すグラフである。
計算値の比較を示すグラフである。
【図5】実施例の実験結果を示すグラフである。
【図6】リアルタイム酸素濃度予測、制御のフローチャ
ートである。
ートである。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 金田 洋 千葉市中央区川崎町1番地 川崎製鉄株式 会社技術研究本部内 (72)発明者 村上 雅宏 千葉市中央区川崎町1番地 川崎製鉄株式 会社技術研究本部内
Claims (2)
- 【請求項1】 CZ法による単結晶の引上げにおいて、
引上げ中のルツボ回転数、融液レベル、結晶固化率の操
業条件と、融液表面温度との関係を用いて、引上げ時に
測定した実際の融液表面温度から、引上げ単結晶中の酸
素濃度を高精度に求めることを特徴とする単結晶の酸素
濃度予測方法。 - 【請求項2】 請求項1により求めた引上げ単結晶中の
酸素濃度により、ルツボ回転数及び融液レベルの操業条
件を結晶固化率に応じて調整し、同時にヒータ温度を変
化させて融液表面温度を調整しながら単結晶を引上げる
ことを特徴とする単結晶の引上方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3-346259 | 1991-12-27 | ||
JP34625991 | 1991-12-27 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH05262593A true JPH05262593A (ja) | 1993-10-12 |
Family
ID=18382192
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP34718792A Withdrawn JPH05262593A (ja) | 1991-12-27 | 1992-12-25 | 単結晶の酸素濃度予測方法及び単結晶の引上方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH05262593A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11228286A (ja) * | 1998-02-13 | 1999-08-24 | Shin Etsu Handotai Co Ltd | 単結晶製造方法 |
EP1734157A1 (en) * | 2005-06-15 | 2006-12-20 | Siltronic AG | Production process of silicon single crystal |
-
1992
- 1992-12-25 JP JP34718792A patent/JPH05262593A/ja not_active Withdrawn
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11228286A (ja) * | 1998-02-13 | 1999-08-24 | Shin Etsu Handotai Co Ltd | 単結晶製造方法 |
EP1734157A1 (en) * | 2005-06-15 | 2006-12-20 | Siltronic AG | Production process of silicon single crystal |
US7582160B2 (en) | 2005-06-15 | 2009-09-01 | Siltronic Ag | Silicone single crystal production process |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8641822B2 (en) | Method and apparatus for controlling diameter of a silicon crystal ingot in a growth process | |
JP4380537B2 (ja) | シリコン単結晶を製造する方法 | |
JP4833747B2 (ja) | 結晶成長されたインゴットの品質評価方法 | |
JP4808917B2 (ja) | 易揮発性不純物をドープしたシリコンからなる単結晶を製造する方法およびこの種の単結晶およびこれから製造した半導体ウェーハ | |
JP3097838B2 (ja) | 結晶成長の制御方法および装置 | |
CN107109687A (zh) | 能够控制锭界面形状的单晶生长***和方法 | |
JPH03159987A (ja) | 単結晶の比抵抗コントロール方法 | |
JP2008280211A (ja) | 単結晶の製造方法 | |
JPH0777994B2 (ja) | 単結晶の酸素濃度コントロール方法及び装置 | |
JP4955237B2 (ja) | 無駄時間をもつ時変系制御対象のための制御システム及び方法 | |
US7582160B2 (en) | Silicone single crystal production process | |
JPH05262593A (ja) | 単結晶の酸素濃度予測方法及び単結晶の引上方法 | |
JP2750935B2 (ja) | 分子線制御方法及び分子線結晶成長装置 | |
US5176787A (en) | Method and apparatus for measuring the diameter of a silicon single crystal | |
JP2787042B2 (ja) | 単結晶引き上げ方法 | |
JP3041990B2 (ja) | 単結晶育成方法及び装置 | |
JP2579761B2 (ja) | 単結晶直径の制御方法 | |
WO1989000211A1 (en) | Growth of semiconductor single crystals | |
JPH01313385A (ja) | 半導体単結晶の直径制御方法 | |
KR20050021382A (ko) | 실리콘 단결정 제조방법 | |
JPH08259380A (ja) | シリコン結晶成長方法 | |
JPH0416436B2 (ja) | ||
JP2007099531A (ja) | 分子線源セルの制御システム、その制御方法および分子線エピタキシャル装置 | |
JPH06316481A (ja) | 結晶引上方法 | |
JP4407539B2 (ja) | 単結晶インゴットの引上げ速度のシミュレーション方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20000307 |