説明

シミュレーション方法、シミュレーションプログラム、半導体製造装置

【課題】紫外線の膜へのダメージ量の正確な予測が行えるようにするシミュレーション方法を提供する。
【解決手段】粒子密度の微分方程式に基づくシミュレーションを行って算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算し、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求め、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから、半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーション方法を実行する。

【発明の詳細な説明】
【技術分野】
【0001】
本技術は、半導体装置の製造における、紫外線によるダメージ量を予測するための、シミュレーション方法及びシミュレーションプログラムに係わる。また、本技術は、紫外線によるダメージ量を予測して、プロセス条件を制御する半導体製造装置に係わる。
【背景技術】
【0002】
半導体装置の製造において、高周波バイアス印加によって生成されるプラズマを利用した加工が広く行われている。
【0003】
プラズマ中では、加工チャンバの構造及びプロセス条件に応じて、粒子同士の衝突やチャンバの壁との相互作用反応が起こり、イオン、ラジカル、光が生成されている。
このうちのイオンが膜へ照射されることで引き起こされる、膜へのダメージ(結晶欠陥)の予測や制御の技術開発については、シミュレーションからプロセス手法、装置構造に至るまで、様々な技術が提案されている。
【0004】
イオンが引き起こすダメージに加えて、近年、プラズマから照射される光、特に紫外線、即ち紫外(UV;ultraviolet)や真空紫外(VUV;vacuum-ultraviolet)による、膜へのダメージが懸念されている。そして、プラズマから照射される紫外線による膜へのダメージに関する基礎研究やモニタリング技術が国内外で活発化してきている。
【0005】
昨今、半導体デバイスが薄膜化し、有機や化合物半導体デバイスの開発が拡大してゆくことにより、紫外領域での発光強度予測、紫外線が引き起こすダメージの定量予測、並びに、制御技術の開発が望まれている。
そして、従来から、半導体装置の製造において、紫外線のモニタ方法や紫外線によるダメージ量予測には、いくつかの方法が提案されている。
【0006】
例えば、ウェハ中やプラズマチャンバの天板に、シリコン膜と酸化膜の積層からなるセンサを設置し、負バイアス印加の下で、センサに入射した紫外線により生成されるホールを誘導電流として検出する方法が提案されている(特許文献1、特許文献2を参照)。
プラズマプロセス中に、この誘導電流の変動を測定することによって、トータル(分光はしていない)の紫外線強度をリアルタイムでモニタリングすることができる。また、この紫外線強度を指標として、ダメージ量を予測する。
【0007】
また、有機物除去のための紫外線洗浄装置において、紫外線ランプから放射される紫外線のモニタリング方法が提案されている(特許文献3を参照)。
洗浄には波長が254nmの紫外線を用いているが、その強い発光強度ゆえに、直接モニタするとセンサがすぐに劣化してしまう問題がある。
そのため、提案されている方法では、直接254nmの波長をモニタせずに、この波長と相関を持ち、かつ、強度の小さい313nmの紫外線強度をリファレンスとしてモニタリングすることにより、本来知りたい波長254nmの強度変化を予測する。
【先行技術文献】
【特許文献】
【0008】
【特許文献1】特開2009−283838号公報
【特許文献2】特開2009−59879号公報
【特許文献3】特開昭62−122130号公報
【発明の概要】
【発明が解決しようとする課題】
【0009】
ところで、紫外線照射による膜へのダメージは、照射される紫外線の波長と対象となる膜のバンドギャップとが関係する。
例えば、シリコン(バンドギャップは1.12eV)であれば、紫外線(エネルギー4eV以上で310nmより短い波長)が入射されると、その紫外線波長全帯で吸収があり、結晶欠陥が引き起こされる。
SiO(8.95eV:139nm)やSiN(5.1eV:243nm)であれば、シリコンよりもバンドギャップが大きいために、310nmより短い波長でも波長帯によって吸収度やダメージ量は大きく異なる。
従って、シリコンよりも上層の膜の構造や厚さ、さらには照射波長によって、シリコンへのダメージ量は大きく異なる。
また、使用するプロセス条件によっても発光スペクトル分布は異なるため、半導体装置の製造への実用という意味では、紫外(UV)から真空紫外(VUV)にかけた分光が重要となる。
そして、実際の半導体装置の製造におけるプラズマプロセス中では、様々な紫外域波長がダメージに関与していると考えられる。
【0010】
しかしながら、特許文献1や特許文献2に記載された方法では、紫外線のトータル量をモニタリングしていることになるが、波長ごとの強度をモニタリングする(分光する)ことはできていない。そのため、ダメージ量の正確な予測・制御が行えない。
また、特許文献3に記載された方法では、特定の2つの波長のみにピンポイントで着目しているため、ダメージ量の正確な予測・制御が行えない。
【0011】
本技術の目的は、紫外線の膜へのダメージ量の正確な予測が行えるようにするシミュレーション方法及びシミュレーションを提供するものである。また、紫外線の膜へのダメージ量を予測して、ダメージ量を制御することができる半導体製造装置を提供するものである。
【課題を解決するための手段】
【0012】
本技術のシミュレーション方法は、半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーション方法である。
まず、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。
そして、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。
さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する。
【0013】
本技術のシミュレーションプログラムは、半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーションプログラムである。
まず、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。
そして、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。
さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する。
そして、これらの機能を情報処理装置に実装して実行させる。
【0014】
本技術の半導体製造装置は、ウェハが配置されるチャンバと、チャンバ内における可視領域の発光スペクトルを検出するセンサと、チャンバ内のプロセス条件を算出する算出部と、チャンバ内のプロセス条件を制御する制御部を備える。
算出部は、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。そして、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、センサで実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の発光スペクトルを予測する。そして、予測した紫外領域の発光スペクトルから紫外線によるダメージ量を予測し、このダメージ量がなるべく少なくなるようにチャンバ内のプロセス条件を算出する。
制御部は、算出部で算出されたプロセス条件となるように、チャンバ内のプロセス条件を制御する。
【0015】
上述の本技術のシミュレーション方法によれば、粒子密度に基づいて計算した、可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測する。これにより、紫外領域の発光スペクトルを正確に予測することができる。
さらに、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測するので、発光する紫外線波長に応じて、発光強度及び紫外線によるダメージ量を正確に予測することができる。
【0016】
上述の本技術のシミュレーションプログラムによれば、粒子密度に基づいて計算した、可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測する。さらに、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する。
これにより、紫外領域の発光スペクトルを正確に予測することができ、発光する紫外線波長に応じて、発光強度及び紫外線によるダメージ量を正確に予測することができる。
【0017】
上述の本技術の半導体製造装置の構成によれば、算出部は、粒子密度に基づいて計算した、可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測する。これにより、紫外領域の発光スペクトルを正確に予測することが可能になる。
そして、算出部は、予測した紫外領域の発光スペクトルから紫外線によるダメージ量を予測する。これにより、発光する紫外線波長に応じて、発光強度及び紫外線によるダメージ量を正確に予測することが可能になる。
さらに、算出部は、このダメージ量がなるべく少なくなるようにチャンバ内のプロセス条件を算出し、制御部は、算出部で算出されたプロセス条件となるように、チャンバ内のプロセス条件を制御する。これにより、ダメージ量がなるべく少なくなるようにチャンバ内のプロセス条件が制御されるので、半導体装置を製造する際の紫外線によるダメージ量を少なくすることが可能になる。
【発明の効果】
【0018】
上述の本技術によれば、紫外線波長に応じて、発光強度及び紫外線によるダメージ量を正確に予測することが可能になるため、効率良く各膜に対する紫外線による変質(ダメージ)を低減させることが可能になる。
従って、本技術により、紫外線による変質(ダメージ)が少なく、良好な特性を有する半導体装置を製造することが可能になる。
【図面の簡単な説明】
【0019】
【図1】本技術に係るシミュレーション方法の基本フロー図である。
【図2】紫外領域の発光スペクトルの定量予測のシミュレーションのフロー図である。
【図3】プラズマ発光シミュレーションの対象となる装置のプラズマチャンバの一部の模式図である。
【図4】図2において、プラズマ発光シミュレーションの発光強度予測からEEDFを求める過程を説明する模式図である。
【図5】ダメージ量データベースの一部の概要である。
【図6】本技術による半導体製造装置の概念図である。
【図7】図6の半導体製造装置において、製造プロセスの条件を補正するフローチャートである。
【図8】第1の実施の形態のシミュレーション方法で使用するパターン構造図である。
【図9】第2の実施の形態のシミュレータの概念図である。
【図10】第3の実施の形態のドライエッチング装置の概略構成図である。
【図11】図10のドライエッチング装置において、製造プロセスの条件を補正するフローチャートである。
【図12】第4の実施の形態のドライエッチング装置の概略構成図(ブロック図)である。
【図13】図12のドライエッチング装置において、製造プロセスの条件を補正するフローチャートである。
【図14】第5の実施の形態のドライエッチング装置において、オーバーエッチングの終点を検出するフローチャートである。
【発明を実施するための形態】
【0020】
以下、本技術を実施するための最良の形態(以下、実施の形態とする)について説明する。
なお、説明は以下の順序で行う。
1.本技術の概要
2.第1の実施の形態(シミュレーション方法)
3.第2の実施の形態(シミュレータ)
4.第3の実施の形態(ドライエッチング装置)
5.第4の実施の形態(ドライエッチング装置)
6.第5の実施の形態(ドライエッチング装置及びエッチング加工方法)
【0021】
<1.本技術の概要>
まず、具体的な実施の形態の説明に先立ち、本技術の概要を説明する。
【0022】
本技術においては、紫外発光を直接モニタリングすることなく、可視発光スペクトルデータに基づいて、紫外発光スペクトルをシミュレーションにより予測する。
また、シミュレーションによる紫外発光スペクトルの予測の結果をインプット情報として、膜へのダメージ量をシミュレーションにより予測する。
さらに、リアルタイムに紫外線ダメージを低減することを目的として、上述した予測方法に基づくソフトウェアを実装した半導体製造装置を構成する。
【0023】
本技術では、下記のように、シミュレーション方法、シミュレーションプログラム、半導体製造装置を構成する。
【0024】
本技術のシミュレーション方法は、半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーション方法である。
まず、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。
そして、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。
この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測する。さらに、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する。
【0025】
本技術のシミュレーションプログラムは、半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーションプログラムである。
まず、まず、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。
そして、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。
この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測する。さらに、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する。
そして、これらの機能をコンピュータに実行させる。
【0026】
本技術の半導体製造装置は、ウェハが配置されるチャンバと、チャンバ内における可視領域の発光スペクトルを検出するセンサと、チャンバ内のプロセス条件を算出する算出部と、チャンバ内のプロセス条件を制御する制御部を備える。
算出部は、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。そして、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、センサで実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の発光スペクトルを予測する。そして、予測した紫外領域の発光スペクトルから紫外線によるダメージ量を予測し、このダメージ量がなるべく少なくなるようにチャンバ内のプロセス条件を算出する。
制御部は、算出部で算出されたプロセス条件となるように、チャンバ内のプロセス条件を制御する。
【0027】
(シミュレーション方法)
本技術に係るシミュレーション方法について、詳細に説明する。
本技術に係るシミュレーション方法の基本フロー図を、図1に示す。
図1に示すフローは、詳細を後述する、シミュレータ(情報処理装置、算出部)により実行される。
図1に示すように、まず、ステップS1において、シミュレータは、可視領域(例えば、300nm−800nm:領域Aと呼ぶ)の発光スペクトルデータをインプットデータとして読み込む。この可視領域の発光スペクトルデータは、実際に、半導体製造装置(加工装置等)のチャンバ内で発光し、OES等のセンサによって検出された、可視領域の発光スペクトルのデータである。
次に、ステップS2において、シミュレータは、プラズマ発光シミュレーションと、レシピ情報とから、紫外領域(例えば、10nm−300nm:領域Bと呼ぶ)の発光スペクトルの定量予測を行う。プラズマ発光シミュレーションの詳細については、後述する。また、レシピ情報は、対象とするプロセス(加工や成膜)における発光種や発光波長群を示す情報である。
次に、ステップS3において、シミュレータは、パターン構造データ、即ち、紫外線が照射される対象であるパターン構造の情報を読み込む。パターン構造データとしては、マスク情報(例えば、GDSファイル情報)と膜厚情報を読み込む。
次に、ステップS4において、シミュレータは、ダメージ量を予測する。ここでは、シミュレータが、ステップS2で予測した紫外領域の発光スペクトルと、読み込んだパターン構造データとから、ダメージ量データベースとの照合を行い、ダメージ量を予測する。
【0028】
さらに、図1のフローの各ステップについて、詳細に説明する。
図1のステップS2における、紫外領域の発光スペクトルの定量予測は、例えば、以下に説明するようにして、行うことができる。
【0029】
紫外領域の発光スペクトルの定量予測のシミュレーションのフロー図を、図2に示す。
図2に示すフローも、詳細を後述する、シミュレータ(情報処理装置、算出部)により実行される。
【0030】
まず、図2の最も左の過程において、プラズマを用いた加工や成膜中に、OES等のセンサで検出した、可視領域のスペクトルデータ101と、プラズマ発光シミュレーションの発光強度予測102を、レシピ情報から得られる代表的な波長λで、比較する。
これにより、詳細を後述するようにして、EEDF(電子エネルギー分布関数;Electron Energy Distribution Function)103を導出する。
プラズマ発光シミュレーションによる発光強度の予測は、文献(特開2011−134927号公報や、Jpn. J. Appl. Phys.,49,(2010),08JD01を参照のこと)に記載されたシミュレーション方法によって、行うことができる。
【0031】
次に、前述した代表的な波長λを発光する、発光種に関する、粒子同士の反応の断面積(σλ)を断面積データベースより照合し、各波長λにおいて、次式のように<σλ’v>104を求める。
【数1】

ここで、f(E)はEEDF、vは電子速度、Eは電子のエネルギーである。
【0032】
次に、得られた<σλ’v>104を用いて、各波長λの紫外発光強度の計算105を行う。
即ち、関係する遷移λ’を考慮することによって、以下の詳細平衡の式から、振動・回転遷移に伴う電子遷移、微細構造遷移、超微細構造遷移による、発光に関する励起準位のpopulationである、[x]λを求める。
【0033】
【数2】

ここで、nは電子密度、σλは反応断面積、Aλは各遷移におけるアインシュタインのA係数である。
【0034】
さらに、この式によって求められた[x]λを用いて、紫外領域の波長λの光の発光強度Iλは、次のように表される。
【0035】
【数3】

ここで、hはプランク係数、νλは対応する周波数、βは逃避率、Sは源泉関数、Pはバックグラウンド関数である。
【0036】
ここで、バックグラウンド成分は無視すると、P=0となる。また、一般的に、半導体製造で用いられているプラズマプロセスでの光学厚さは薄いことから、近似的にβ≒1が成り立つ。すると、上の式は、下記のように表すことができる。
【数4】

【0037】
このようにして、各波長λの紫外発光強度が得られ、紫外領域(領域B)のスペクトル分布が得られる。
そして、得られた紫外領域(領域B)のスペクトル分布と、OES等のセンサで検出した、可視領域(領域A)のスペクトルデータ101を合わせて、紫外領域及び可視領域の発光スペクトル分布106が得られる。
【0038】
次に、プラズマ発光シミュレーションによる発光強度の予測の具体的な手法について、以下に概略を説明する。
プラズマ発光シミュレーションの手法及び条件を説明するために、プラズマ発光シミュレーションの対象となる装置(プラズマを用いて加工や成膜を行う装置)のプラズマチャンバの一部の模式図を、図3に示す。
【0039】
図3に示すように、装置1は、チャンバ基体3から成るプラズマチャンバ2と、プラズマチャンバ2のチャンバ壁9内に形成された空間と、この空間内に設けられた下部電極5Bと、チャンバ基体3に設けられた上部電極5Aとを有する。
このプラズマ発光シミュレーションでは、プラズマチャンバ2内の空間を区画することにより、複数のメッシュ空間2mが設定される。そして、メッシュ空間2m毎に、シミュレーションにより粒子密度nが算出される。
【0040】
チャンバ壁9には、当該チャンバ壁9を区画することにより、複数のメッシュ壁9mが設定される。そして、粒子の付着確率Sが、メッシュ壁9m毎に設定される。
メッシュ壁9m毎に付着確率Sを設定することにより、チャンバ壁9の各部の事情に応じて適切に付着確率Sを設定することができる。例えば、チャンバ壁9において、チャンバ基体3により構成される部分(SiO)と、上部電極5Aにより構成される部分(Si)とで、異なる付着確率Sを設定することができる。
【0041】
メッシュ空間2m及びメッシュ壁9mの形状及び大きさは適宜に設定されてよい。
メッシュ空間2mの(端面の)形状及び大きさとメッシュ壁9mの形状及び大きさは、互いに同一でもよいし、異なっていてもよい。メッシュ空間2m及びメッシュ壁9mの大きさは、例えば、数mmオーダーである。
【0042】
シミュレーションは、2次元シミュレーションでも3次元シミュレーションでもよい。換言すれば、メッシュ空間2mは、2次元的に設定されてもよいし、3次元的に設定されてもよく、メッシュ壁9mは、1次元的に設定されてもよいし、2次元的に設定されてもよい。
【0043】
以下では、メッシュ空間2mの識別番号をpとし、メッシュ空間2m毎に値が変化するパラメータ(pの関数)については、(p)を付すことがある。
また、メッシュ壁9mの識別番号をwとし、メッシュ壁9m毎に値が変化するパラメータ(wの関数)については、(w)を付すことがある。
【0044】
シミュレーションにおいては、メッシュ空間2mと、チャンバ壁9との距離を示すパラメータも考慮される。具体的には、メッシュ空間2mと側壁2aとの(最短)距離L(p)、メッシュ空間2mと天板2bとの(最短)距離R(p)、メッシュ空間2mとメッシュ壁9mとの距離d(w)が考慮される。なお、d(w)もpの関数であるが、pの付加は省略する。距離L(p)、R(p)、d(w)は、例えば、メッシュ空間2mの中央位置からの距離である。また、d(w)は、例えば、メッシュ壁9mの中央位置からの距離である。
【0045】
このように、チャンバ壁9との距離が考慮されることにより、粒子密度nが適切に算出される。特に、メッシュ壁9m毎に、付着確率S(w)及び距離d(w)が設定されることにより、粒子密度nが適切に算出される。
【0046】
粒子密度nの微分方程式は、例えば、以下の数式を用いる。
【0047】
【数5】

【0048】
この数式は、上述したように、メッシュ空間2m毎に用いられる。n(i,t)及びn(j,t)は、種類i及びjの粒子の時間tにおける粒子密度である。粒子の種類i(j)は、例えば、ガスが水素の場合、H、H、H、H、H、eである。以下では、i,jに代えて、H、H等を( )内に示すことがある。
は反応レート、τは排気特性時間、τは壁消失特性時間である。vは粒子速度、Dは拡散定数、Pは圧力、μは移動度、kはボルツマン定数、Tは電子温度、eは素電荷である。
【0049】
この数式の第1行目の式において、第1項は、粒子同士の衝突反応による粒子の生成及び消失による密度変化を示している。第2項は、粒子のプラズマチャンバ2からの排気による密度変化を示している。第3項は、粒子とチャンバ壁9との相互作用及び粒子のウェハでの消失による密度変化を示している。
【0050】
数式の第1行目及び第2行目の式に示されるように、メッシュ空間2mとメッシュ壁9mとの距離d(w)の値が小さくなるほど粒子の消失が多くなるように、距離d(w)が組み込まれている。
【0051】
次に、算出した粒子密度n(i)に基づいて、発光に寄与する粒子Hの粒子密度n(H)を算出する。なお、Hは、HやH等の粒子の種類(上記の数式における種類i,j)とは別個の概念であり、粒子密度n(H)は、上記の数式では求められない。
具体的には、メッシュ空間2m毎に、数式で求められた粒子密度n(i)をパラメータとして含む、粒子密度n(H)の微分方程式をシミュレーションにより解くことにより、粒子密度n(H)を算出する。
【0052】
次に、メッシュ空間2m毎の粒子密度n(H)に基づいて、メッシュ空間2m毎の発光強度I(p)を算出する。
そして、図3においてハッチングして示すように、メッシュ空間2m毎の発光強度I(p)を所定の視線方向y1において積分し、プラズマチャンバ2を視線方向y1に見たときの発光強度Iを算出する。
このように、視線方向y1における発光強度Iを算出することにより、OES10により検出される視線方向y1における発光強度Iを再現することができる。
なお、視線方向y1は、水平方向、鉛直方向、これらに傾斜する方向、メッシュ空間2mの配列方向、又は、当該配列方向に傾斜する方向等、適宜な方向に設定されてよい。また、積分される空間の、視線方向y1における断面積も、適宜に設定されてよい。積分される空間の、視線方向y1に直交する方向における位置は、適宜に設定されてよく、また、複数位置において視線方向y1における発光強度Iが算出されてもよい。
【0053】
上述した手法のプラズマ発光シミュレーションによって、各波長におけるシミュレーション発光強度Iλが得られる。即ち、図2のプラズマ発光シミュレーションの発光強度予測102が得られる。
また、図2において、プラズマ発光シミュレーションの発光強度予測102からEEDF103を求める過程を説明する模式図を、図4に示す。
この発光強度予測102の各波長におけるシミュレーション発光強度Iλは、図4に示すように、装置構造パラメータα、プロセス条件パラメータβ、電子温度Tの関数として表される。
ここで、例えば、電子温度Tを0〜10eVの範囲で変化させて得られる発光強度Iλと発光強度の実測値IOλを代表的な波長λ(レシピ情報から、比較対象の発光種を判断する)で比較する。そして、発光強度が所望の誤差範囲内(例えば、±10%)を満たすような電子温度Tを探す。その電子温度Tを用いてMaxwell分布等を仮定して、EEDF(電子エネルギー分布関数)103を導出する。
このようにして、図2に示したEEDF103が導出されるので、このEEDF103を用いて、図2に示した過程を経て、紫外領域のスペクトル分布を得ることができる。
【0054】
図1のステップS4における、ダメージ量の予測は、ダメージ量データベースを利用する。このダメージ量データベースは、膜種(Si,SiO,SiN,有機膜、化合物)に対して、特定波長ごとに、パターン構造、ドーズ量、プロセス条件(ガス種、流量、圧力、パワー)の各データと、ダメージ量(結晶欠陥)のデータとの組み合わせから成る。
このダメージ量データベースの一部の概要を、図5に示す。図5では、Siと有機膜に対して、それぞれ、特定の波長1、波長2、波長3ごとに、パターン構造、ドーズ量、プロセス条件のデータとダメージ量のデータとの組み合わせを示している。
なお、図5では、波長1に1組のパターン構造、ドーズ量、プロセス条件のデータとダメージ量のデータとの組み合わせを記載しているが、実際のダメージ量データベースでは、波長1、波長2、波長3のそれぞれに多数の組の組み合わせがある。
ダメージ量データベースは、あらかじめ実測によって作成しても良いし、MD(Molecular Dynamics)等による数値シミュレーション等で作成しても良い。
【0055】
図5に示したような構成のダメージ量データベースを使用することにより、紫外線の波長及び強度と、プロセス条件・パターン構造・開口率との関連付けを行って、ダメージ量を予測することができる。
これにより、予測したダメージ量に基づいて、紫外線によるダメージをできるだけ少なく抑制する、プロセス条件やデバイス構造(又はパターンレイアウト)の最適化を定量的に行うことが可能になる。
【0056】
そして、シミュレータは、計算で得られた紫外領域の発光スペクトルと、図1のステップS3で読み込んだ、紫外線が照射される実際のパターン構造の情報とを用いて、図5に示したようなダメージ量データベースと照合し、スプライン補間する。これにより、紫外線ダメージ量を予測する。
このようにして紫外線ダメージ量を予測して、図1に示す処理が終了する。
【0057】
上述した本技術のシミュレーション方法によれば、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。そして、対象とする製造プロセスにおける発光種及び発光波長の情報(レシピ情報)を参照して、計算した可視領域の各波長の発光強度と、センサ実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求める。さらに、この電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測している。これにより、紫外領域の発光スペクトルを正確に予測することができる。
さらに、予測して得られた紫外領域の発光スペクトルから、紫外線によるダメージ量を予測するので、プロセス中に発光する紫外線波長に応じて、発光強度及び紫外線によるダメージ量を正確に予測することができる。
【0058】
(半導体製造装置)
本技術による半導体製造装置の概念図を、図6に示す。また、本技術による半導体製造装置のプロセス補正動作のフローチャートを、図7に示す。
図6に示す半導体製造装置11は、プラズマチャンバ12、可視領域発光モニタであるOES、シミュレーションを実行するためのシミュレーションプログラム(ソフトウェア)を実装したシミュレータ14、制御部15を有して、構成されている。シミュレータ14は、コンピュータ等の情報処理装置であり、本技術に係る「算出部」に該当する。
【0059】
プラズマチャンバ12内には、下部に電極13が設けられている。この電極13上に、ウェハ16を載置して、ウェハ16に対して成膜や膜の加工が行われる。
なお、図示しないが、プラズマチャンバ12内の上部にも電極(上部電極)を設けており、電極(下部電極)13と上部電極との間に電界をかけて、プラズマを発生させた状態で、ウェハ16に対して成膜や膜の加工が行われる。
【0060】
シミュレータ14に実装させるシミュレーションプログラム(ソフトウェア)は、オンライン及び/又はオフラインで動作可能な構成とする。
そして、シミュレータ14は、図6の右側の図に示すように、OESにより取得したモニタリングデータから、紫外領域のスペクトルを予測し、ダメージ量を予測して、最適化されたプロセス条件を探し出す、操作を実行する。
なお、シミュレータ14に実装させるシミュレーションプログラム(ソフトウェア)を、ダメージ量の予測の後で、さらに、破線で囲まれた、形状変動予測を行って、最適化されたプロセス条件を探し出す構成としても良い。
【0061】
制御部15は、シミュレータ14において実行されたシミュレーションで得られた、プラズマチャンバ12のプロセス条件の補正を、プラズマチャンバ12へフィードバックするように制御を行う。
なお、可視領域発光モニタとしては、図6に示したOES以外の構成も採用することが可能である。
【0062】
そして、図6に示した半導体製造装置11では、例えば、図7に示すフローチャートに従って、製造プロセスの条件を補正する。
図7に示すように、まず、ステップS31において、プラズマ加工中に、プラズマチャンバ12に備え付けられたOESによって、可視領域の発光スペクトルデータ(例えば、波長300−800nm)を逐次取得する。
次に、ステップS32において、シミュレータ14は、前述したプラズマ発光シミュレーションと、レシピ情報から、紫外領域の発光スペクトルを予測する。
次に、ステップS33において、シミュレータ14は、紫外線が照射される対象である、パターン構造のデータ(マスク情報と膜厚情報等)を読み込む。
次に、ステップS34において、シミュレータ14は、前述したシミュレーション方法により、紫外線ダメージ量を予測する。ここでは、シミュレータ14が、ステップS32で予測した紫外領域の発光スペクトルと、読み込んだパターン構造データとから、ダメージ量データベースとの照合を行い、ダメージ量を予測する。
次に、ステップS35において、シミュレータ14は、ダメージ量データベースから、レシピ条件に近い範囲内(例えば、±50%)でダメージを低減するプロセス条件(圧力、流量、パワー)を選定する。
次に、図6及び図7において破線で囲われた、形状変動予測を行わない場合には、そのままステップS37に進み、シミュレータ14は、選定したプロセス条件の補正パラメータを、制御部15に受け渡す。
【0063】
これに対して、図6及び図7において破線で囲われた、形状変動予測を行う場合には、ステップS35からステップS36に進んで、シミュレータ14が形状変動予測を行う。具体的には、シミュレータ14が、選定されたプロセス条件及び、ウェハ開口率・Semi-Local開口率・パターン構造を用いて、形状シミュレーションを行うことにより、形状変動を予測する。形状シミュレーションは、例えば、特開2009−152269号公報に記載された方法を採用することができる。
予測の結果、変動が許容範囲内(例えば、±10%以内)であれば、ステップS37に進み、シミュレータ14は、選定したプロセス条件の補正パラメータを、制御部15に受け渡す。
一方、変動が許容範囲を超える場合には、ステップS35に戻り、シミュレータ14が、ダメージを低減することができ、かつ、前回選定したプロセス条件とは別のプロセス条件を選定する。そして、ステップS36に進み、シミュレータ14が、選定したプロセス条件に対して、形状変動予測を行う、
これを繰り返すことで、シミュレータ14は、スペックを満たす形状と紫外線ダメージ量の低減を行うことのできる最適プロセス条件を探し出す。
プロセス条件が決定したら、シミュレータ14は、決定したプロセス条件を、制御部15を通じてプラズマチャンバ12側に送り、プロセス条件の補正をリアルタイムに行う。
【0064】
このように、本技術による半導体製造装置では、加工中のプラズマ状態やチャンバ壁状態を反映する、可視領域発光の情報を用いて紫外線スペクトルの予測を行うため、これまでの手法よりも、正確に紫外線ダメージ量を制御することが可能となる。これにより、半導体デバイスの電気特性、即ち、トランジスタI−V特性、イメージャの白点・暗電流、太陽電池の変換効率、レーザーの出力特性等の向上を行うことが可能になる。
また、形状変動予測も行う場合には、プロセス条件をさらに最適化することが可能となる。
【0065】
なお、本技術による半導体製造装置は、プラズマを用いる加工であれば、ドライエッチング、CVD(Chemical Vapor Deposition)、PVD(Physical Vapor Deposition)等、加工の種類を問わず、適用することができる。
【0066】
<2.第1の実施の形態(シミュレーション方法)>
続いて、具体的な実施の形態を説明する。
まず、第1の実施の形態のシミュレーション方法を説明する。
本実施の形態のシミュレーション方法は、OESによるモニタリングデータを用いて、紫外発光スペクトル及びダメージ量のシミュレーションを行う。
【0067】
本実施の形態のシミュレーション方法で使用する、パターン構造図を、図8に示す。
図8に示すパターン構造は、裏面照射型イメージセンサ(画素サイズ1.5μm)のフォトダイオードPDからOCL(On Chip Lens)までの各層である。
フォトダイオードPDは厚さ3000nm、その上の層間膜21は厚さ1500nm、その上の3色R,G,Bのカラーフィルター22は厚さ700nm、その上のOCL23は厚さ700nmである。
図8は、このパターン構造を作製する際の、OCL23のドライエッチング加工の際に生じる真空紫外線VUVが、OCL23、カラーフィルター22、層間膜21中を減衰・透過し、フォトダイオードPDへダメージを与える概念図である。ここでは、真空紫外線VUVは、波長10〜200nmとする。
本実施の形態では、この状況で真空紫外線VUVがフォトダイオードPDへ与えるダメージ量をシミュレーション予測する。
OCL23の加工に用いる加工装置の構造情報とプロセス条件を、以下に示す。
【0068】
(加工装置)
CCP(Capacitive Coupled Plasma)ドライエッチング装置
周波数:上部=60MHz、下部=2MHz
チャンバ直径:100cm
上下電極ギャップ:40mm
可視光モニタ装置:OES
【0069】
(プロセス条件)
ガス種類:CF
ガス流量:150sccm
圧力:60mTorr
上部印加パワー:1000W
エッチング時間:15分
【0070】
可視光モニタ装置であるOES(波長帯域:200−800nm、波長分解能:0.5nm)は、加工装置のチャンバの壁のポート(窓材は、例えば、サファイア)越しに設置する。
このOESにより、エッチング中にプラズマから発せられる発光スペクトルを、逐次データ取得する。
本実施の形態では、上述した加工装置の情報、構造情報(OCL23、カラーフィルター22、層間膜21の各層の厚さ)、上述したプロセス条件の情報、OESにより取得した可視発光スペクトル情報を、入力情報とする。
電子エネルギー分布を導出するために必要なスペクトル波長としては、例えば、321nm(CFラジカル起因)、240nm(CFラジカル起因)、686nm(Fラジカル起因)を用いることができる。ただし、強度が大きければ、波長はこの限りではない。
そして、これらの入力情報と、図1に示した本技術に係るダメージ量の予測方法を用いる。
これにより、本実施の形態では、波長200nm以下の紫外領域のスペクトル、並びにフォトダイオードPDの表層に形成されるシリコンの結晶欠陥量をシミュレーション予測することができる。
【0071】
上述の本実施の形態のシミュレーション方法によれば、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した粒子密度に基づいて、可視領域の各波長の発光強度を計算する。そして、レシピ情報(対象とする製造プロセスにおける発光種及び発光波長の情報)を参照して、計算した可視領域の各波長の発光強度と、実際にOESにより取得した可視領域の発光スペクトルとを入力情報として用いている。これらの入力情報から、本技術に係るダメージ量の予測方法を用いて、紫外領域の発光スペクトルを予測し、予測して得られた紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する。これにより、プロセス中に発光する紫外線波長に応じて、発光強度及び紫外線によるダメージ量を正確に予測することができる。
【0072】
なお、本実施の形態ではCCPドライエッチング装置に本技術を適用したが、本技術に係るシミュレーション方法を適用する加工装置は、ICP(Inductive Coupled Plasma)装置やECR(Electron Cyclotron Resonance)装置等、種類は問わない。
また、本技術に係るシミュレーション方法を適用する加工工程は、図8に示したパターン構造のOCL23のドライエッチング以外にも、トランジスタ、やSTI(Shallow Trench Isolation)、サイドウォール、コンタクトホール等、種類は問わない。
さらにまた、本技術に係るシミュレーション方法を適用するプロセス条件についても、ターゲット膜を加工できる条件であれば、特に限定されない。
【0073】
<3.第2の実施の形態(シミュレータ)>
続いて、第2の実施の形態のシミュレータを説明する。
本実施の形態のシミュレータ(情報処理装置、算出部)は、本技術に係るダメージ量を予測するシミュレーションの機能を実行させるコンピュータプログラム(ソフトウェア)を実装しており、本技術に係るシミュレーション方法を実行する。
【0074】
第2の実施の形態のシミュレータの概念図を、図9に示す。
図9に示すシミュレータ30は、入力部31、紫外スペクトル演算部32、プロセス条件別の発光の波長データベース33、ダメージ量演算部34、ダメージ量データベース35、出力部36から構成されている。
【0075】
シミュレータ30は、エッチング装置等の加工装置40と信号やデータをやり取りする。
そして、シミュレータ30内では、図9の矢印に従って処理が進む。
【0076】
入力部31には、入力情報として各種のデータDATAが入力される。
この各種のデータDATAとしては、パターン構造データ(例えば、GDSデータ)、膜厚データ、装置パラメータ、プロセスパラメータ・可視発光スペクトルデータ、等が挙げられる。
【0077】
紫外スペクトル演算部32では、入力部31に入力されたデータDATAを用い、このデータDATAに基づいて、プロセス条件別の発光の波長データベース33を参照して、対応する波長データを取得し、紫外領域の発光スペクトルを求める演算を行う。プロセス条件別の発光の波長データベース33は、プロセス条件毎に、それぞれ波長群をデータとして持っている。
ダメージ量演算部34では、紫外スペクトル演算部32で得た紫外領域の発光スペクトルを用い、この紫外領域の発光スペクトル及び入力されたデータDATAに基づいてダメージ量データベース35を参照し、ダメージ量を求める演算を行う。
出力部36では、求めたダメージ量をシミュレータ30から外部に出力する。
【0078】
本実施の形態のシミュレータに実装させるコンピュータプログラム(ソフトウェア)は、好ましくは、GUI(グラフィカルユーザーインターフェイス)から操作を行うことによって、データ入力から計算の実行、計算経緯や結果の表示がなされる構成とする。GUIとしては、例えば、tcl/tkやMotif等で構築されたGUIを採用することができるが、これらに限らない。
また、本実施の形態のシミュレータに実装させるコンピュータプログラム(ソフトウェア)の動作プラットフォームは、OS(Operation System)の種類を問わず、各種OSに適用することが可能である。
さらに、本実施の形態の本実施の形態のシミュレータに実装させるコンピュータプログラム(ソフトウェア)は、オフラインでの動作でもよいし、オンラインでシミュレータ30と接続されていても良い。
【0079】
<4.第3の実施の形態(ドライエッチング装置)>
次に、第3の実施の形態のドライエッチング装置を説明する。
本実施の形態のドライエッチング装置は、紫外線によるダメージ量を最適化するように、プロセス条件を自動制御する構成である。
第3の実施の形態のドライエッチング装置の概略構成図(ブロック図)を、図10に示す。
【0080】
図10に示すドライエッチング装置41は、プラズマチャンバ42、OES、シミュレーションを実行するためのシミュレーションプログラム(ソフトウェア)を実装したシミュレータ44、制御部45を有して、構成されている。シミュレータ44は、コンピュータ等の情報処理装置であり、本技術に係る「算出部」に該当する。
【0081】
プラズマチャンバ42内には、下部に電極43が設けられている。この電極43上に、ウェハ46を載置して、ウェハ46に対してエッチング加工が行われる。
なお、図示しないが、プラズマチャンバ42内の上部にも電極(上部電極)を設けており、電極(下部電極)43と上部電極との間に電界をかけて、プラズマを発生させた状態で、ウェハ46に対してエッチング加工が行われる。
OESは、プラズマチャンバ42の壁のポート窓越しに設置されている。
OESは、例えば、波形分解能が0.5nmであり、波長帯域が300nm−800nmである構成とすることができる。
【0082】
本実施の形態のドライエッチング装置41において、エッチングの際の条件としては、例えば、上部電極の周波数を60MHz、下部電極の周波数を13.56MHz、プラズマチャンバ42の直径を100cm、上下の電極のギャップを40mmとする。
【0083】
シミュレータ44に実装させるシミュレーションプログラム(ソフトウェア)は、オンライン及び/又はオフラインで動作可能な構成とする。
そして、シミュレータ44は、図10の右側の図に示すように、OESにより取得したモニタリングデータから、紫外領域のスペクトルを予測し、ダメージ量を予測して、最適化されたプロセス条件を探し出す、操作を実行する。
【0084】
制御部45は、シミュレーションで得られた、プラズマチャンバ42のプロセス条件の補正を、プラズマチャンバ42へフィードバックするように制御を行う。
なお、可視領域発光モニタとしては、図10に示したOES以外の構成も採用することが可能である。
【0085】
そして、図10に示したドライエッチング装置41では、例えば、図11に示すフローチャートに従って、製造プロセスの条件を補正する。
図11に示すように、まず、ステップS41において、ドライエッチング加工中に、プラズマチャンバ42に備え付けられたOESによって、例えば1秒毎に、可視領域の発光スペクトルデータ(例えば、波長300−800nm)を取得する。
次に、ステップS42において、シミュレータ44は、前述したプラズマ発光シミュレーションと、レシピ情報から、紫外領域の発光スペクトルを予測する。
次に、ステップS43において、シミュレータ44は、紫外線が照射される対象である、パターン構造のデータ(マスク情報と膜厚情報等)を読み込む。
次に、ステップS44において、シミュレータ44は、前述したシミュレーション方法により、紫外線ダメージ量を予測する。ここでは、シミュレータ44が、予測した紫外領域の発光スペクトルと、読み込んだパターン構造データとから、ダメージ量データベースとの照合を行い、ダメージ量を予測する。
次に、ステップS45において、シミュレータ44は、ダメージ量データベースから、レシピ条件に近い範囲内(例えば、±50%)でダメージを低減するように、プロセス条件(圧力、流量、パワー)を選定する。
次に、ステップS46において、シミュレータ44は、選定したプロセス条件の補正パラメータを、制御部45に受け渡す。
これらのステップを、エッチングの終了まで、逐次繰り返すことにより、紫外線ダメージ量を低減したドライエッチングが可能になる。
【0086】
そして、ステップS47において、シミュレータ44は、エッチング時間が終了したかを判断する。
エッチング時間が終了している場合には、エッチング処理を終了させる。
エッチング時間が終了していない場合には、ステップS41に戻り、OESによって、可視領域の発光スペクトルデータを取得する。
このようにして、紫外線ダメージ量を低減したドライエッチングが可能となる。
【0087】
なお、本実施の形態において、ドライエッチング装置は、CCP装置、ICP装置、ECR装置、のいずれも使用することが可能である。
また、CVD(化学的気相成長法)装置やPVD(物理的気相成長法)装置等のプラズマを用いた半導体加工装置にも、本実施の形態を適用することができる。
【0088】
上述の本実施の形態のドライエッチング装置41によれば、シミュレータ44において、紫外領域の発光スペクトルの予測とダメージ量の予測が行われ、ダメージ量を最適化するプロセス条件が選定される。そして、制御部45において、シミュレーションで得られた、プラズマチャンバ42のプロセス条件の補正を、プラズマチャンバ42へフィードバックするように制御が行われる。
これにより、ダメージ量を最適化する(なるべく少なくなる)ようにプロセスチャンバ42のプロセス条件が制御されるので、ドライエッチングの際の紫外線によるダメージ量を少なくすることが可能になる。
従って、ドライエッチングの際の紫外線によるダメージが少なく、良好な特性を有する半導体装置を製造することが可能になる。
【0089】
<5.第4の実施の形態(ドライエッチング装置)>
次に、第4の実施の形態のドライエッチング装置を説明する。
本実施の形態のドライエッチング装置は、形状変動も考慮して、紫外線によるダメージ量を最適化するように、プロセス条件を自動制御する構成である。
第4の実施の形態のドライエッチング装置の概略構成図(ブロック図)を、図12に示す。
【0090】
図12に示すドライエッチング装置51は、プラズマチャンバ52、OES、シミュレーションを実行するためのシミュレーションプログラム(ソフトウェア)を実装したシミュレータ54、制御部55を有して、構成されている。シミュレータ54は、コンピュータ等の情報処理装置であり、本技術に係る「算出部」に該当する。
【0091】
プラズマチャンバ52内には、下部に電極53が設けられている。この電極53上に、ウェハ56を載置して、ウェハ56に対してエッチング加工が行われる。
なお、図示しないが、プラズマチャンバ52内の上部にも電極(上部電極)を設けており、電極(下部電極)53と上部電極との間に電界をかけて、プラズマを発生させた状態で、ウェハ56に対してエッチング加工が行われる。
OESは、プラズマチャンバ52の壁のポート窓越しに設置されている。
OESは、例えば、波形分解能が0.5nmであり、波長帯域が300nm−800nmである構成とすることができる。
【0092】
本実施の形態のドライエッチング装置51において、エッチングの際の条件としては、例えば、上部電極の周波数を60MHz、下部電極の周波数を13.56MHz、プラズマチャンバ52の直径を100cm、上下の電極のギャップを40mmとする。
【0093】
シミュレータ54に実装させるシミュレーションプログラム(ソフトウェア)は、オンライン及び/又はオフラインで動作可能な構成とする。
そして、シミュレータ54は、図12の右側の図に示すように、OESにより取得したモニタリングデータから、紫外領域のスペクトルを予測し、ダメージ量を予測する、操作を実行する。
シミュレータ54は、ダメージ量を予測した後に、さらに、形状変動を予測して、最適化されたプロセス条件を探し出す、操作を実行する。
【0094】
制御部55は、シミュレーションで得られた、プラズマチャンバ52のプロセス条件の補正を、プラズマチャンバ52へフィードバックするように制御を行う。
なお、可視領域発光モニタとしては、図12に示したOES以外の構成も採用することが可能である。
【0095】
そして、図12に示したドライエッチング装置51では、例えば、図13に示すフローチャートに従って、製造プロセスの条件を補正する。
図13に示すように、まず、ステップS51において、ドライエッチング加工中に、プラズマチャンバ52に備え付けられたOESによって、例えば1秒毎に、可視領域の発光スペクトルデータ(例えば、波長300−800nm)を取得する。
次に、ステップS52において、シミュレータ54は、前述したプラズマ発光シミュレーションと、レシピ情報から、紫外領域の発光スペクトルを予測する。
次に、ステップS53において、シミュレータ54は、紫外線が照射される対象である、パターン構造のデータ(マスク情報と膜厚情報等)を読み込む。
次に、ステップS54において、シミュレータ54は、前述したシミュレーション方法により、紫外線ダメージ量を予測する。ここでは、シミュレータ54が、予測した紫外領域の発光スペクトルと、読み込んだパターン構造データとから、ダメージ量データベースとの照合を行い、ダメージ量を予測する。
次に、ステップS55において、シミュレータ54は、ダメージ量データベースから、レシピ条件に近い範囲内(例えば、±50%)でダメージを低減するように、プロセス条件(圧力、流量、パワー)を選定する。
【0096】
続いて、ステップS56に進んで、シミュレータ54が形状変動予測を行う。具体的には、シミュレータ54が、選定されたプロセス条件及び、ウェハ開口率・Semi-Local開口率・パターン構造を用いて、形状シミュレーションを行うことにより、形状変動を予測する。なお、Semi-Local開口率は、注目するパターン周辺の数十mm領域(例えば、30mm。プロセス条件によって異なり、粒子の平均自由行程の数倍程度の値である。)における、レジストマスクの開口率である。
形状シミュレーションは、例えば、特開2009−152269号公報に記載された方法を採用することができる。
予測の結果、形状変動が許容範囲内であれば、ステップS57に進み、シミュレータ54は、選定したプロセス条件の補正パラメータを、制御部55に受け渡す。
一方、形状変動が許容範囲を超える場合には、ステップS55に戻り、シミュレータ54が、ダメージを低減することができ、かつ、前回選定したプロセス条件とは別のプロセス条件を選定する。そして、ステップS56に進み、シミュレータ54は、選定したプロセス条件に対して、形状変動予測を行う、
プロセス条件が決定したら、シミュレータ54は、決定したプロセス条件を、制御部55を通じてプラズマチャンバ52側に送り、プロセス条件の補正をリアルタイムに行う。
これらのステップを、エッチングの終了まで、逐次繰り返すことにより、スペックを満たす形状の制御と、紫外線ダメージ量の低減とが可能になる。
【0097】
次に、ステップS58において、シミュレータ54は、エッチング時間が終了したかを判断する。
エッチング時間が終了している場合には、エッチング処理を終了させる。
エッチング時間が終了していない場合には、ステップS51に戻り、OESによって、可視領域の発光スペクトルデータを取得する。
このようにして、形状変動を抑制し、かつ紫外線ダメージ量を低減したドライエッチングが可能となる。
【0098】
なお、本実施の形態において、ドライエッチング装置は、CCP装置、ICP装置、ECR装置、のいずれも使用することが可能である。
また、CVD(化学的気相成長法)装置やPVD(物理的気相成長法)装置等のプラズマを用いた半導体加工装置にも、本実施の形態を適用することができる。
【0099】
上述の本実施の形態のドライエッチング装置51によれば、シミュレータ54において、紫外領域の発光スペクトルの予測とダメージ量の予測が行われ、ダメージ量を最適化するプロセス条件が選定される。そして、制御部55において、シミュレーションで得られた、プラズマチャンバ52のプロセス条件の補正を、プラズマチャンバ52へフィードバックするように制御が行われる。
これにより、ダメージ量を最適化する(なるべく少なくなる)ようにプロセスチャンバ52のプロセス条件が制御されるので、ドライエッチングの際の紫外線によるダメージ量を少なくすることが可能になる。
従って、ドライエッチングの際の紫外線によるダメージが少なく、良好な特性を有する半導体装置を製造することが可能になる。
【0100】
また、上述の本実施の形態のドライエッチング装置51によれば、さらに、シミュレータ54において、紫外領域の発光スペクトルの予測とダメージ量の予測の間に、形状変動の予測が行われる。
従って、ドライエッチングの際の形状の変動が少なく、良好な特性を有する半導体装置を安定して歩留まり良く製造することが可能になる。
【0101】
<6.第5の実施の形態(ドライエッチング装置及びエッチング加工方法)>
次に、第5の実施の形態のドライエッチング装置を説明する。また、このドライエッチング装置を用いたエッチング加工方法を説明する。
本実施の形態のドライエッチング装置は、オーバーエッチングの終点の検出を、ダメージ量のシミュレーションで制御する構成である。
【0102】
ドライエッチング装置の各部の構成は、図6、図10、図12に示した装置と同様とすることができる。
【0103】
本実施の形態のドライエッチング装置では、例えば、図14に示すフローチャートに従って、オーバーエッチングの終点を検出する。
図14に示すように、ステップS61において、ドライエッチング装置は、メインエッチングを開始する。
メインエッチングが進行していくと、ステップS62において、ドライエッチング装置は、メインエッチングの終点検出を行う。
メインエッチングの終点は、例えば、前述したOESによって可視領域の発光をモニタリングすることによって、検出することができる。
【0104】
メインエッチングの終点を検出した後に、ステップS63において、ドライエッチング装置は、オーバーエッチングを開始する。
次に、ステップS64において、ドライエッチング装置は、紫外線によるダメージ量を予測する。ダメージ量の予測は、図1に示したシミュレーションフローに従い行う。
【0105】
続いて、ステップS65において、ドライエッチング装置は、予め設定されたダメージ量に到達したかどうかを判断する。
設定されたダメージ量に到達している場合には、ステップS66に進み、ドライエッチング装置は、オーバーエッチングの終点検出を行う。
設定されたダメージ量に到達していない場合には、ステップS64に戻り、ドライエッチング装置は、オーバーエッチングを続行する。そして、再びダメージ量の予測を行う。
このようにして、オーバーエッチングの終点を検出することができる。
図14に示したオーバーエッチングの終点を検出するフローチャートは、図6、図10、図12に示した装置のシミュレータ14,44,54に実装されたコンピュータプログラム(ソフトウェア)によって、実行することができる。
【0106】
ここで、本実施の形態のエッチングの終点検出を、半導体装置の製造プロセスにおけるエッチングに適用する場合の具体例を示す。
この具体例では、シリコン基板上のSiOが厚さ2nm、その上のポリシリコンが厚さ150nm、その上に形成するレジストマスクが厚さ300nmであり、レジストマスクを用いたエッチングによって、ポリシリコンをゲートのパターンにパターニングする。
初期のレジストマスクの線幅は、100nmとする。
この具体例において、加工に用いる装置とプロセス条件を、以下に示す。
【0107】
(加工装置)
CCP(Capacitive Coupled Plasma)ドライエッチング装置
周波数:上部=60MHz、下部=13.56MHz
チャンバ直径:100cm
上下電極ギャップ:60mm
可視光モニタ装置:OES
【0108】
(プロセス条件)
ガス種類:HBr/O
ガス流量:500sccm/5sccm
圧力:30mTorr
上部印加パワー:150W
下部印加パワー:60W
エッチング時間:メインエッチング100秒+オーバーエッチング
【0109】
ポリシリコンのメインエッチングの終点検出は、従来から行われているように、OES(波長帯域:300−800nm、波長分解能:0.5nm)により、可視領域の発光を指標にして、その発光強度の変動によって行う。本具体例では、例えば、SiBrの波長526nm発光を指標にする。
また、オーバーエッチングについては、紫外線ダメージシミュレーションで予測されるシリコン基板へのダメージ量が許容値(例えば、シリコン結晶欠陥密度が1017cm)を超えたときに、終点検出とする。
また、電子エネルギー分布を導出するための代表スペクトル波長は、例えば、484nm(COラジカル起因)、777nm(Oラジカル起因)、656nm(Hラジカル起因)を用いることができる。ただし、強度が大きければ、波長はこの限りではない。
これにより、シリコン基板へのダメージ量を可視化しながら、ポリシリコンのオーバーエッチングを行うことが可能となり、シリコン基板へのダメージを低減することができる。
【0110】
上述の本実施の形態のエッチング加工方法によれば、紫外線によるダメージ量を指標にして、オーバーエッチングの終点検出を行うので、リアルタイムにダメージ量を確認しながらダメージを低減したプラズマ加工が可能になる。
【0111】
なお、上述した第5の実施の形態のエッチング方法は、加工条件やその時のリファレンスとする可視領域発光を最適化することで、ゲートエッチングに限らず、STIやサイドウォール、コンタクトホール等のエッチングにも適用可能である。
また、第5の実施の形態のエッチング方法は、有機膜や化合物半導体(GaNやSiGe等)といった、その他の膜のオーバーエッチング終点検出にも適用できる。
また、第5の実施の形態の終点検出を、第3の実施の形態のエッチング装置や第4の実施の形態のエッチング装置に具備させても構わない。
【0112】
以上、上述の実施の形態は本技術の好ましい様態を示すものである。
本技術の技術的範囲は、上述した実施の形態に限定されるものではなく、紫外線以外の波長領域への応用も同様に可能である。さらに、ドライエッチングに限定されるものではなく、プラズマを利用する半導体製造プロセス及び半導体製造装置に対して広く適用できるものである。
【0113】
なお、本開示は以下のような構成も取ることができる。
(1)半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーション方法であって、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した前記粒子密度に基づいて、可視領域の各波長の発光強度を計算し、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求め、前記電子エネルギー分布関数と、前記発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測するシミュレーション方法。
(2)予測した紫外領域の発光スペクトルと、作製するパターン構造のデータから、紫外線によるダメージ量を予測する前記(1)に記載のシミュレーション方法。
(3)膜の材料、紫外線の波長、パターン構造、紫外線のドーズ量、プロセス条件と、ダメージ量を関連付けたダメージ量データベースを予め作成しておき、紫外線によるダメージ量を予測する際に、前記ダメージ量データベースと照合を行う、前記(2)に記載のシミュレーション方法。
(4)半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーションプログラムであって、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した前記粒子密度に基づいて、可視領域の各波長の発光強度を計算し、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求め、前記電子エネルギー分布関数と、前記発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する、機能を情報処理装置に実装して実行させるシミュレーションプログラム。
(5)ウェハが配置されるチャンバと、前記チャンバ内における可視領域の発光スペクトルを検出するセンサと、粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した前記粒子密度に基づいて、可視領域の各波長の発光強度を計算し、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、前記センサで実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求め、前記電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから紫外線によるダメージ量を予測し、前記ダメージ量がなるべく少なくなるように前記チャンバ内のプロセス条件を算出する算出部と、前記算出部で算出されたプロセス条件となるように、前記チャンバ内のプロセス条件を制御する制御部を備えた半導体製造装置。
(6)前記算出部は、紫外線によるダメージ量を予測し、前記チャンバ内のプロセス条件を算出した後に、形状シミュレーションによってプロセス中の形状変動を予測し、形状変動が許容範囲を超える場合には、先に算出したプロセス条件とは別のプロセス条件を算出する、前記(5)に記載の半導体製造装置。
(7)前記算出部は、さらに、メインエッチング後のオーバーエッチングの終点検出を、予測した紫外線によるダメージ量が、予め設定したダメージ量に到達したことにより行う、前記(5)又は前記(6)に記載の半導体製造装置。
【符号の説明】
【0114】
11 半導体製造装置、12,42,52 プラズマチャンバ、13,43,53 電極、14,44,54 シミュレータ、15,45,55 制御部、16,46,56 ウェハ、21 層間膜、22 カラーフィルター、23 OCL、30 シミュレータ、31 入力部、32 紫外スペクトル演算部、34 ダメージ量演算部、35 ダメージ量データベース、36 出力部、40 加工装置、41,51 ドライエッチング装置、PD フォトダイオード、VUV 真空紫外線

【特許請求の範囲】
【請求項1】
半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーション方法であって、
粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、
算出した前記粒子密度に基づいて、可視領域の各波長の発光強度を計算し、
対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求め、
前記電子エネルギー分布関数と、前記発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測し、
予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する
シミュレーション方法。
【請求項2】
予測した紫外領域の発光スペクトルと、作製するパターン構造のデータから、紫外線によるダメージ量を予測する請求項1に記載のシミュレーション方法。
【請求項3】
膜の材料、紫外線の波長、パターン構造、紫外線のドーズ量、プロセス条件と、ダメージ量を関連付けたダメージ量データベースを予め作成しておき、紫外線によるダメージ量を予測する際に、前記ダメージ量データベースと照合を行う、請求項2に記載のシミュレーション方法。
【請求項4】
半導体装置の製造における、紫外線によるダメージ量を予測するシミュレーションプログラムであって、
粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した前記粒子密度に基づいて、可視領域の各波長の発光強度を計算し、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求め、前記電子エネルギー分布関数と、前記発光種に関する反応断面積とを用いて、紫外領域の各波長の発光強度を計算することにより、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから、紫外線によるダメージ量を予測する、機能を情報処理装置に実装して実行させる
シミュレーションプログラム。
【請求項5】
ウェハが配置されるチャンバと、
前記チャンバ内における可視領域の発光スペクトルを検出するセンサと、
粒子密度の微分方程式に基づくシミュレーションを行って、粒子密度を算出し、算出した前記粒子密度に基づいて、可視領域の各波長の発光強度を計算し、対象とする製造プロセスにおける発光種及び発光波長の情報を参照して、計算した可視領域の各波長の発光強度と、前記センサで実際に検出された可視領域の発光スペクトルとを比較して、電子エネルギー分布関数を求め、前記電子エネルギー分布関数と、発光種に関する反応断面積とを用いて、紫外領域の発光スペクトルを予測し、予測した紫外領域の発光スペクトルから紫外線によるダメージ量を予測し、前記ダメージ量がなるべく少なくなるように前記チャンバ内のプロセス条件を算出する算出部と、
前記算出部で算出されたプロセス条件となるように、前記チャンバ内のプロセス条件を制御する制御部を備えた
半導体製造装置。
【請求項6】
前記算出部は、紫外線によるダメージ量を予測し、前記チャンバ内のプロセス条件を算出した後に、形状シミュレーションによってプロセス中の形状変動を予測し、形状変動が許容範囲を超える場合には、先に算出したプロセス条件とは別のプロセス条件を算出する、請求項5に記載の半導体製造装置。
【請求項7】
前記算出部は、さらに、メインエッチング後のオーバーエッチングの終点検出を、予測した紫外線によるダメージ量が、予め設定したダメージ量に到達したことにより行う、請求項5に記載の半導体製造装置。

【図1】
image rotate

【図5】
image rotate

【図6】
image rotate

【図7】
image rotate

【図8】
image rotate

【図9】
image rotate

【図10】
image rotate

【図11】
image rotate

【図12】
image rotate

【図13】
image rotate

【図14】
image rotate

【図2】
image rotate

【図3】
image rotate

【図4】
image rotate


【公開番号】特開2013−115354(P2013−115354A)
【公開日】平成25年6月10日(2013.6.10)
【国際特許分類】
【出願番号】特願2011−262417(P2011−262417)
【出願日】平成23年11月30日(2011.11.30)
【出願人】(000002185)ソニー株式会社 (34,172)
【Fターム(参考)】