(57)【要約】【課題】結晶性を有する部材に複数回イオン注入をする場合のイオン分布を求める際の誤差を少なくする。【解決手段】イオン分布特定ステップ３０は、既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定する。イオン分布算出ステップ３１は、第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出する。差分計算ステップ３２は、イオン分布算出ステップ３１によって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、イオン分布特定ステップ３０によって特定された既存のイオン分布との差分を計算する。イオン分布計算ステップ３３は、差分計算ステップ３２によって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算する。【選択図】　　　　図１
2004022965.tif 000002
【特許請求の範囲】
【請求項１】
結晶性を有する部材に複数回に分けてイオンを注入する場合に当該部材内部のイオン分布を計算するイオン分布計算方法において、
既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定するイオン分布特定ステップと、
第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出するイオン分布算出ステップと、
前記イオン分布算出ステップによって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、前記イオン分布特定ステップによって特定された既存のイオン分布との差分を計算する差分計算ステップと、
前記差分計算ステップによって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算するイオン分布計算ステップと、
を有することを特徴とするイオン分布計算方法。
【請求項２】
前記複数回のイオン注入は、全ドーズ量をＮ回（Ｎ≧ｎ）に均等に分割して注入することを特徴とする請求項１記載のイオン分布計算方法。
【請求項３】
前記複数回のイオン注入は、全ドーズ量を任意の回数に分割して注入することを特徴とする請求項１又は２記載のイオン分布計算方法。
【請求項４】
前記イオン注入は、複数の種類が異なるイオンを注入することを特徴とする請求項１、２又は３記載のイオン分布計算方法。
【請求項５】
前記イオン分布特定ステップは、注入パスに沿って既存のイオン分布を特定し、
前記イオン分布計算ステップは、前記注入パスに沿って第（ｎ＋１）回目のイオン分布を計算する、
ことを特徴とする請求項１、２、３又は４記載のイオン分布計算方法。
【請求項６】
前記イオン分布は、アモルファス領域と、チャネリング領域とのそれぞれについてイオンの分布であることを特徴とする請求項１、２、３、４又は５記載のイオン分布計算方法。
【請求項７】
前記結晶性を有する部材は、半導体基板であることを特徴とする請求項１、２、３、４、５又は６記載のイオン分布計算方法。
【請求項８】
結晶性を有する部材に複数回に分けてイオンを注入する場合に当該部材内部のイオン分布をコンピュータに計算させるコンピュータ読み取り可能なプログラムにおいて、
コンピュータを、
既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定するイオン分布特定ステップ、
第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出するイオン分布算出ステップ、
前記イオン分布算出ステップによって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、前記イオン分布特定ステップによって特定された既存のイオン分布との差分を計算する差分計算ステップ、
前記差分計算ステップによって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算するイオン分布計算ステップ、
として機能させるコンピュータ読み取り可能なプログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明はイオン分布計算方法およびプログラムに関し、特に、結晶性を有する部材に複数回に分けてイオンを注入する場合に当該部材内部のイオン分布を計算するイオン分布計算方法およびそのような計算をコンピュータに実行させるプログラムに関する。
【０００２】
【従来の技術】
半導体装置の製造プロセスにおいては、半導体基板に、例えば、砒素や硼素などのイオンを注入し、部材の電気的性質を変更することが一般的に行われている。
【０００３】
ＭＯＳ−ＦＥＴ（ｍｅｔａｌ　ｏｘｉｄｅ　ｓｅｍｉｃｏｎｄｕｃｔｏｒ−Ｆｉｅｌｄ　Ｅｆｆｅｃｔ　Ｔｒａｎｓｉｓｔｏｒ）の製造プロセスにおいては、Ｅｘｔｅｎｓｉｏｎ，Ｓｏｕｒｃｅ／Ｄｒａｉｎ領域へ４方向からイオンを注入する場合がある。
【０００４】
図９は、このようなイオン注入の態様を説明するための図である。この図に示すように、ＭＯＳ−ＦＥＴの製造プロセスにおいて、イオン注入をする際には、図の左右方向（図中矢印で示す方向）および前後方向の４方向（図示せず）からイオンを注入し、ｎ＋領域を形成する。
【０００５】
このようなイオン注入を行う際には、注入されたイオンが所望の分布となるようにイオンのエネルギーと、ドーズ量とを決定する必要があることから、パーソナルコンピュータやワークステーション（以下単にコンピュータと称する）を用いてイオン分布をシミュレーションすることが従来から行われてきた。
【０００６】
ところで、結晶構造を有する部材にイオンを注入する場合、イオンが衝突した領域では、結晶構造が壊れてアモルファス（ａｍｏｒｐｈｏｕｓ）構造に転化することが知られている。
【０００７】
結晶構造が壊れてしまうと、それよりも奥にはイオンが到達しにくくなることから、アモルファス領域よりも奥の領域ではドーズ量を増加してもイオンの濃度が増加しなくなってしまう。
【０００８】
図１０は、このような現象を説明するための図である。この図の横軸は、イオンが注入される部材の表面からの深さを示し、また、縦軸は注入されたイオンの濃度を示している。また、各曲線は、ドーズ量を変化させた場合の結果を示している。
【０００９】
この図に示すように、ドーズ量が増加するに従って、グラフのピーク部分（アモルファス領域に対応する部分）は略比例して増加しているのに対し、それよりも奥に位置する部分（チャネリング領域）はドーズ量が増加してもほとんど変化していない。
【００１０】
【発明が解決しようとする課題】
従来において、イオン注入を複数回行う場合のシミュレーションについては、１回分のイオン注入を単純に注入の回数分だけ累積加算して求めるべき結果としていた。
【００１１】
アモルファス領域については、前述のように線形性を有するので重ね合わせの定理が成立し、得られる結果も実際の値に近い値が得られていた。
しかし、チャネリング領域については、線形性を有しないことから重ね合わせの定理が成立せず、単純に加算しただけでは得られる結果が実際の値から大きく乖離したものとなっていた。
【００１２】
本発明はこのような点に鑑みてなされたものであり、結晶構造を有する部材に複数回に分けてイオン注入を行う場合でも実際の値に近いシミュレーション結果を得ることが可能なイオン分布計算方法およびそのようなイオン分布計算方法を実行するプログラムを提供することを目的とする。
【００１３】
【課題を解決するための手段】
本発明では上記課題を解決するために、図１に示す、結晶性を有する部材に複数回に分けてイオンを注入する場合に当該部材内部のイオン分布を計算するイオン分布計算方法において、既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定するイオン分布特定ステップ３０と、第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出するイオン分布算出ステップ３１と、前記イオン分布算出ステップ３１によって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、前記イオン分布特定ステップ３０によって特定された既存のイオン分布との差分を計算する差分計算ステップ３２と、前記差分計算ステップ３２によって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算するイオン分布計算ステップ３３と、を有することを特徴とするイオン分布計算方法が提供される。
【００１４】
ここで、イオン分布特定ステップ３０は、既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定する。イオン分布算出ステップ３１は、第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出する。差分計算ステップ３２は、イオン分布算出ステップ３１によって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、イオン分布特定ステップ３０によって特定された既存のイオン分布との差分を計算する。イオン分布計算ステップ３３は、差分計算ステップ３２によって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算する。
【００１５】
また、本発明では上記課題を解決するために、結晶性を有する部材に複数回に分けてイオンを注入する場合に当該部材内部のイオン分布をコンピュータに計算させるコンピュータ読み取り可能なプログラムにおいて、コンピュータを、既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定するイオン分布特定ステップ、第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出するイオン分布算出ステップ、前記イオン分布算出ステップによって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、前記イオン分布特定ステップによって特定された既存のイオン分布との差分を計算する差分計算ステップ、前記差分計算ステップによって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算するイオン分布計算ステップ、として機能させるコンピュータ読み取り可能なプログラムが提供される。
【００１６】
ここで、イオン分布特定ステップ３０は、既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定する。イオン分布算出ステップ３１は、第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出する。差分計算ステップ３２は、イオン分布算出ステップ３１によって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、イオン分布特定ステップ３０によって特定された既存のイオン分布との差分を計算する。イオン分布計算ステップ３３は、差分計算ステップ３２によって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算する。
【００１７】
【発明の実施の形態】
以下、本発明の実施の形態を図面を参照して説明する。
図１は、本発明の動作原理を説明する原理図である。この図に示すように、本発明のイオン分布計算方法は、イオン分布特定ステップ３０、イオン分布算出ステップ３１、差分計算ステップ３２およびイオン分布計算ステップ３３を有している。
【００１８】
ここで、イオン分布特定ステップ３０は、既に実行されたｎ回（ｎは自然数）のイオン注入によって生じた既存のイオン分布を特定する。
イオン分布算出ステップ３１は、第（ｎ＋１）回目のイオン注入によって生じると想定されるイオン分布を算出する。
【００１９】
差分計算ステップ３２は、イオン分布算出ステップ３１によって算出された第（ｎ＋１）回目のイオン注入によるイオン分布と、イオン分布特定ステップ３０によって特定された既存のイオン分布との差分を計算する。
【００２０】
イオン分布計算ステップ３３は、差分計算ステップ３２によって求められた差分から第（ｎ＋１）回目のドーズ量を求めて第（ｎ＋１）回目のイオン分布を計算する。
【００２１】
次に、以上の原理図の動作について説明する。なお、以下では、合計で５回のイオン注入を行う場合であって、第３回目のイオン注入を求める場合を例に挙げて説明する。
【００２２】
先ず、イオン分布特定ステップ３０は、前回、即ち、第２回目までのイオン注入によって形成されたイオン分布の状態を特定し、特定された情報をイオン分布算出ステップ３１に供給する。
【００２３】
イオン分布算出ステップ３１は、次回、即ち、第３回目のイオン注入によって生じると想定されるイオン分布を算出し、差分計算ステップ３２に供給する。
差分計算ステップ３２は、イオン分布算出ステップ３１によって算出された第３回目のイオン注入によるイオン分布と、イオン分布特定ステップ３０によって特定された第２回目までのイオン注入によって形成されたイオン分布との差分を計算し、イオン分布計算ステップ３３に供給する。
【００２４】
イオン分布計算ステップ３３は、差分計算ステップ３２から供給された差分のドーズ量を、イオン注入した際に得られるイオン分布を計算する。そして、イオン分布計算ステップ３３は、得られた結果を、第３回目のイオン注入によって得られたイオン分布の計算結果として出力する。
【００２５】
以上に説明したように、本発明のイオン分布計算方法によれば、既存のイオン分布と、次のイオン注入によって生じると予想されるイオン分布の差分に基づいて注入イオンのドーズ量を計算し、得られたドーズ量に基づいて次のイオン注入によるイオン分布の結果を計算するようにしたので、従来の方法のように１回のイオン注入の結果を複数倍にする場合に比較して、チャネリング領域のイオン分布を正確にシミュレーションすることが可能になる。
【００２６】
次に、本発明の実施の形態について説明する。
図２は、本発明の実施の形態の構成例を示す図である。この図に示すように、本発明のイオン分布計算方法を実施するコンピュータ５０は、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）５０ａ、ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）５０ｂ、ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）５０ｃ、ＨＤＤ（Ｈａｒｄ　Ｄｉｓｋ　Ｄｒｉｖｅ）５０ｄ、ＧＣ（ＧｒａｐｈｉｃｓＣａｒｄ）５０ｅ、Ｉ／Ｆ（Ｉｎｔｅｒｆａｃｅ）５０ｆおよびバス５０ｇによって構成されており、その外部には表示装置５１および入力装置５２が接続されている。
【００２７】
ここで、ＣＰＵ５０ａは、ＨＤＤ５０ｄに格納されているプログラムに従って各種演算処理を実行するとともに、装置の各部を制御する。
ＲＯＭ５０ｂは、ＣＰＵ５０ａが実行する基本的なプログラムやデータを保持している。
【００２８】
ＲＡＭ５０ｃは、ＣＰＵ５０ａが実行対象とするプログラムやデータを一時的に格納する。
ＨＤＤ５０ｄは、ＣＰＵ５０ａが実行するプログラムやデータを格納している。また、ＣＰＵ５０ａの演算処理によって生成されたデータを格納する。
【００２９】
ＧＣ５０ｅは、ＣＰＵ５０ａから供給された描画データに基づいて描画処理を実行し、得られた画像を映像信号に変換して表示装置５１に出力する。
Ｉ／Ｆ５０ｆは、入力装置５２から出力されたデータの表現形式を内部の表現形式に変換して入力する。
【００３０】
バス５０ｇは、ＣＰＵ５０ａ、ＲＯＭ５０ｂ、ＲＡＭ５０ｃ、ＨＤＤ５０ｄ、ＧＣ５０ｅおよびＩ／Ｆ５０ｆを相互に接続し、これらの間でデータの授受を可能にする。
【００３１】
表示装置５１は、例えば、ＬＣＤ（Ｌｉｑｕｉｄ　Ｃｒｙｓｔａｌ　Ｄｉｓｐｌａｙ）またはＣＲＴ（Ｃａｔｈｏｄｅ　Ｒａｙ　Ｔｕｂｅ）ディスプレイによって構成されており、ＧＣ５０ｅから出力された映像信号を表示出力する。
【００３２】
入力装置５２は、例えば、マウスやキーボードによって構成されており、ユーザの操作に応じたデータを生成して出力する。
次に、以上の実施の形態の動作について説明する。
【００３３】
以下では、先ず、イオン分布の計算方法について説明した後、図２に示す実施の形態の動作について説明する。
ドーズ量Φにより半導体基板に対してイオン注入を行った場合における濃度分布Ｎ（ｘ）は、以下の式によって表される。
【００３４】
【数１】
2004022965.tif 000003
【００３５】
ここで、ｘは、半導体基板の表面からの距離（深さ）を、ｒはドーズレシオを、ｎａはアモルファス領域の濃度を、ｎｃはチャネリング領域の濃度を示している。換言すると、本発明の計算方法では、図３に示すように、イオン分布をアモルファス領域におけるイオン分布ｎａと、チャネリング領域におけるイオン分布ｎｃとに分割して計算し、得られた結果を合計することにより、全体のイオン分布を得る方法を採用している。
【００３６】
なお、式（１）は、等価的に以下のように表すことができる。
【００３７】
【数２】
2004022965.tif 000004
【００３８】
なお、Φｃｈａｎはチャネリング領域へのドーズ量を示している。ここで、ｒおよびΦｃｈａｎはドーズ量Φに依存することから、これらをそれぞれｒ（Φ）およびΦｃｈａｎと表すことにする。すると、両者は、以下のような関係を有することになる。
【００３９】
【数３】
2004022965.tif 000005
【００４０】
ここで、全ドーズ量Φをｎ回に均等に分割して注入する場合について考える。この場合、各回において注入されるドーズ量はΦ／ｎとなる。式（２）に式（３）を代入すると、各回におけるイオン分布は以下のようになる。
【００４１】
【数４】
2004022965.tif 000006
【００４２】
ところで、従来の方法においては、ドーズ量がΦ／ｎのイオン注入をｎ回繰り返す場合には、ドーズ量がΦ／ｎである場合のイオン分布を求めこれを単純にｎ回分だけ加算することにより計算結果を得ていた。即ち、次のような式によってｎ回分のイオン分布を得ていた。
【００４３】
【数５】
2004022965.tif 000007
【００４４】
しかし、このような方法では、図３に示すｎｃの領域が単純にｎ倍にされてしまうため、正確なシミュレーション結果を得ることができない。
そこで、本発明の実施の形態では、先ず、図４に示すように、ドーズ量Φとチャネリング領域のドーズ量Φｃｈａｎとの関係を求める。そして、第（ｉ＋１）番目のイオン分布を求める場合には、第ｉ番目のイオン分布と第（ｉ＋１）番目のイオン分布の差分を計算し、得られた差分をチャネリング領域へ注入されるドーズ量としてイオン分布を計算する。
【００４５】
即ち、第（ｉ＋１）回目のイオン注入により新たに追加されるイオンの濃度分布ΔＮｉ（ｘ）は、以下の式によって表される。
【００４６】
【数６】
2004022965.tif 000008
【００４７】
従って、第（ｉ＋１）回目の注入は、ドーズ量Φ／ｎであり、また、チャネリング領域のドーズ量が以下の式で示す注入分布であると擬制することができる。
【００４８】
【数７】
2004022965.tif 000009
【００４９】
このようにして、各回目におけるイオン注入の結果に応じて変化するイオン分布を計算することが可能になる。
次に、以上の実施の形態の動作を具体的に説明する。図５は、図２に示す実施の形態の動作を説明するためのフローチャートである。この処理が開始されると、以下のステップが実行される。
【００５０】
ステップＳ１０：
ＣＰＵ５０ａは、入力装置５２を介してイオンの注入条件を入力する。なお、イオンの注入条件としては、例えば、全イオンのドーズ量Φ、ドーズレシオｒおよび図４に示すドーズ量とチャネリング領域のドーズ量Φｃｈａｎとの関係を示すデータ等を入力する。
【００５１】
ステップＳ１１：
ＣＰＵ５０ａは、既存のイオン注入分布のドーズ分布を評価する。即ち、式（４）に示すＮｉ（ｘ）を求める。
【００５２】
ステップＳ１２：
ＣＰＵ５０ａは、今回注入するドーズ量を用いて差分ドーズを評価する。即ち、式（６）を用いて差分ドーズΔＮｉ（ｘ）を計算する。
【００５３】
ステップＳ１３：
ＣＰＵ５０ａは、ステップＳ１２で求めた差分ドーズを、ステップＳ１１で求めた既存のイオン分布に加算することにより、第（ｉ＋１）番目の注入によるイオン分布を求める。
【００５４】
ステップＳ１４：
ＣＰＵ５０ａは、イオン注入処理を繰り返す場合には、ステップＳ１１に戻って前述の場合と同様の処理を繰り返し、それ以外の場合にはステップＳ１５に進む。
【００５５】
ステップＳ１５：
ＣＰＵ５０ａは、最終的なイオン分布を発生し、対応する描画命令を生成してＧＣ５０ｅに供給する。
【００５６】
ステップＳ１６：
ＧＣ５０ｅは、ステップＳ１５においてＣＰＵ５０ａから供給された計算結果を表示装置５１に表示させる。
【００５７】
以上の処理によれば、チャネリング領域に誤差を生じることなく、イオン分布を正確に計算することが可能になる。
なお、比較のために従来の計算方法の処理の流れについても説明する。
【００５８】
図６は、従来の計算方法の処理の流れを説明するためのフローチャートである。この図に示す処理が開始されると、以下のステップが実行される。
ステップＳ３０：
ＣＰＵ５０ａは、入力装置５２を介してイオンの注入条件を入力する。
【００５９】
ステップＳ３１：
ＣＰＵ５０ａは、イオン注入分布を発生する。即ち、式（５）の何れかの式に示す１回分のイオン注入分布を求める。
【００６０】
ステップＳ３２：
ＣＰＵ５０ａは、ステップＳ３１で得られたイオン分布を既存のイオン注入分布に加算する。
【００６１】
ステップＳ３３：
ＣＰＵ５０ａは、イオン注入処理を繰り返す場合には、ステップＳ３１に戻って前述の場合と同様の処理を繰り返し、それ以外の場合にはステップＳ３４に進む。
【００６２】
ステップＳ３４：
ＣＰＵ５０ａは、最終の分布を発生する。
ステップＳ３５：
ＣＰＵ５０ａは、計算結果に対応する描画命令をＧＣ５０ｅに供給する。その結果、計算結果を示す画面が表示装置５１に表示される。
【００６３】
図７は、本発明の計算方法と、従来の計算方法によってそれぞれ得られた濃度分布を示すグラフである。このグラフの横軸は、半導体基板の表面からの深さを示しており、縦軸は濃度を示している。また、破線（ｏｌｄ　ｄ−ｒａｔｉｏｎ）は、従来の方法によるグラフを示しており、点線（ｎｅｗ　ｄ−ｒａｔｉｏｎ）は、本発明の方法によるグラフを示している。
【００６４】
このグラフから分かるように、従来の方法では、本発明の方法と比較して、深い位置にある領域（チャネリング領域）の値が大きくなっている。これは、前述のようにチャネリング領域の誤差に起因するものである。
【００６５】
図８は、本発明の計算方法と従来の計算方法のそれぞれによって得られたイオン分布に基づいてシミュレーションされたトランジスタのゲート長と、閾値との関係を示すグラフである。
【００６６】
このグラフの横軸は、トランジスタのゲート長を示している。また、縦軸は、トランジスタの閾値を示している。また、破線（ｏｌｄ　ｄ−ｒａｔｉｏｎ）は、従来の方法によるグラフを示しており、点線（ｎｅｗ　ｄ−ｒａｔｉｏｎ）は、本発明の方法によるグラフを示している。
【００６７】
この図から分かるように、従来の計算方法と、本発明の計算方法の誤差は、ゲート長が短くなるほど顕著となっている。近年では、半導体装置上に形成される素子の微細化が進んでいるので、誤差の影響はますます顕著となるものと考えられるため、本発明の計算方法の利用価値が増大するものと考えられる。
【００６８】
なお、以上に説明した実施の形態では、ドーズ量を均等にｎ分割して注入する場合について説明したが、例えば、任意に分割して計算することも可能である。そのような方法について以下に説明する。
【００６９】
即ち、ドーズ量Φ１でイオン注入を行った後に、ドーズ量Φ２でイオン注入を行った場合、あるべきイオン分布は以下の式によって表される。
【００７０】
【数８】
2004022965.tif 000010
【００７１】
ここで、第２回目のイオン注入後のイオン分布は、以下の式によって表される。
【００７２】
【数９】
2004022965.tif 000011
【００７３】
従って、以上の式を用いて、前述の場合と同様に計算をすることにより、任意に分割した場合の結果を得ることができる。
また、以上の実施の形態では、同一のイオンを注入する場合を例に挙げて説明したが複数のイオンを注入する場合に対しても本発明を適用することが可能である。そのような方法を以下に説明する。
【００７４】
例えば、ドーズ量Φｘにより所定のイオンを注入した後、ドーズ量Φにより異種のイオンを注入する場合を考える。この場合、イオン分布は以下の式で表される。
【００７５】
【数１０】
2004022965.tif 000012
【００７６】
従って、複数のイオンを使用してイオン注入を行う場合であっても、これらのイオンによるチャネルドーズの分布が分かれば、前述の場合と同様にイオン分布を計算することが可能になる。
【００７７】
また、以上の実施の形態では、イオン濃度を深さのパラメータｘのみを用いて表すようにしたが、イオンの注入方向が異なれば同一の深さにおいても濃度が異なる場合が生じ得る。その場合には、深さ方向に直交する平面の位置を示すパラメータ（例えば、ｙ）を用いてイオン濃度を示すことも可能である。
【００７８】
なお、そのような場合には、イオンを注入するパスに沿ってイオン濃度の分布を計算するようにすればよい。
最後に、上記の処理機能は、コンピュータによって実現することができる。その場合、イオン分布計算方法を実施する装置が有すべき機能の処理内容を記述したプログラムが提供される。そのプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリなどがある。磁気記録装置には、ハードディスク装置（ＨＤＤ）、フレキシブルディスク（ＦＤ）、磁気テープなどがある。光ディスクには、ＤＶＤ（Ｄｉｇｉｔａｌ　Ｖｅｒｓａｔｉｌｅ　Ｄｉｓｋ）、ＤＶＤ−ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）、ＣＤ−ＲＯＭ（Ｃｏｍｐａｃｔ　Ｄｉｓｋ　Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）、ＣＤ−Ｒ（Ｒｅｃｏｒｄａｂｌｅ）／ＲＷ（ＲｅＷｒｉｔａｂｌｅ）などがある。光磁気記録媒体には、ＭＯ（Ｍａｇｎｅｔｏ−Ｏｐｔｉｃａｌ　ｄｉｓｋ）などがある。
【００７９】
プログラムを流通させる場合には、たとえば、そのプログラムが記録されたＤＶＤ、ＣＤ−ＲＯＭなどの可搬型記録媒体が販売される。また、プログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することもできる。
【００８０】
プログラムを実行するコンピュータは、たとえば、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、自己の記憶装置に格納する。そして、コンピュータは、自己の記憶装置からプログラムを読み取り、プログラムに従った処理を実行する。なお、コンピュータは、可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することもできる。また、コンピュータは、サーバコンピュータからプログラムが転送される毎に、逐次、受け取ったプログラムに従った処理を実行することもできる。
【００８１】
【発明の効果】
以上説明したように本発明では、既存のイオン分布を特定し、次のイオン注入によって生じるイオン分布を算出し、既存のイオン分布と次のイオン注入によって生じるイオン分布との差分を求め、この差分に基づいて求められるドーズ量に応じてイオン分布を計算するようにしたので、チャネリング領域において発生する誤差を少なくすることが可能になる。
【００８２】
また、以上説明したように本発明では、コンピュータに、既存のイオン分布を特定させ、次のイオン注入によって生じるイオン分布を算出させ、既存のイオン分布と次のイオン注入によって生じるイオン分布との差分を求めさせ、この差分に基づいて求められるドーズ量に応じてイオン分布を計算させるようにしたので、イオン分布を求める際に、チャネリング領域において発生する誤差を少なくすることが可能になる。
【図面の簡単な説明】
【図１】本発明の動作原理を説明するための原理図である。
【図２】本発明の実施の形態の構成例を示す図である。
【図３】本発明の計算方法の原理を説明するための図である。
【図４】ドーズ量とチャネルドーズ量との関係を示す図である。
【図５】本発明の計算方法の処理の流れを説明するためのフローチャートである。
【図６】従来の計算方法の処理の流れを説明するためのフローチャートである。
【図７】本発明の計算方法と、従来の計算方法によって得られたイオン分布を示す図である。
【図８】本発明の計算方法と、従来の計算方法によって得られたゲート長と閾値との関係を示す図である。
【図９】イオンの注入の様子を示す図である。
【図１０】イオン分布のドーズ量との関係を示す図である。
【符号の説明】
３０　イオン分布特定ステップ
３１　イオン分布計算ステップ
３２　差分計算ステップ
３３　イオン分布計算ステップ
５０　コンピュータ
５０ａ　ＣＰＵ
５０ｂ　ＲＯＭ
５０ｃ　ＲＡＭ
５０ｄ　ＨＤＤ
５０ｅ　ＧＣ
５０ｆ　Ｉ／Ｆ
５０ｇ　バス
５１　表示装置
５２　入力装置
【図１】
2004022965.tif 000013
【図２】
2004022965.tif 000014
【図３】
2004022965.tif 000015
【図４】
2004022965.tif 000016
【図５】
2004022965.tif 000017
【図６】
2004022965.tif 000018
【図７】
2004022965.tif 000019
【図８】
2004022965.tif 000020
【図９】
2004022965.tif 000021
【図１０】
2004022965.tif 000022
【手続補正書】
【提出日】平成15年7月3日(2003.7.3)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００４０
【補正方法】変更
【補正の内容】
【００４０】
ここで、全ドーズ量Φをｎ回に均等に分割して注入する場合について考える。この場合、各回において注入されるドーズ量はΦ／ｎとなる。式（２）のΦにΦ／ｎを代入すると、各回におけるイオン分布は以下のようになる。


【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００６３
【補正方法】変更
【補正の内容】
【００６３】
図７は、本発明の計算方法と、従来の計算方法によってそれぞれ得られた濃度分布を示すグラフである。このグラフの横軸は、半導体基板の表面からの深さを示しており、縦軸は濃度を示している。また、破線（ｏｌｄ　ｄ−ｒａｔｉｏ）は、従来の方法によるグラフを示しており、点線（ｎｅｗ　ｄ−ｒａｔｉｏ）は、本発明の方法によるグラフを示している。


【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００６６
【補正方法】変更
【補正の内容】
【００６６】
このグラフの横軸は、トランジスタのゲート長を示している。また、縦軸は、トランジスタの閾値を示している。また、点線（ｏｌｄ　ｄ−ｒａｔｉｏ）は、従来の方法によるグラフを示しており、破線（ｎｅｗ　ｄ−ｒａｔｉｏ）は、本発明の方法によるグラフを示している。