(57)【特許請求の範囲】
【請求項１】
電子線を１度に露光できる領域をサブフィールドとし前記電子線の偏向によって前記サブフィールドを被加工基板上の所定位置に順次露光することで所定パターンを形成するパターン形成方法において、
前記所定パターンの全領域を一定の大きさを有する複数個の第１区分領域に分割し、一つの前記第１区分領域での前記電子線が透過するパターン面積が前記一つの第１区分領域全体に占める割合をパターン面積密度として前記所定パターン領域全体のパターン面積密度を前記第１区分領域のそれぞれに対応して構成される２次元の配列データとする第１の面積密度マップ作成工程と、
前記第１区分領域の幾つかのパターン面積密度を用いて計算された値を第２区分領域のパターン面積密度とする第２の面積密度マップを作成する工程と、
前記サブフィールドを露光する時の電子線の偏向位置で定まる電子線のボケ量を前記第１又は前記第２の面積密度マップと前記電子線の電流密度により修正した前記電子線のボケ量と、前記第１及び前記第２の面積密度マップとは異なる第３の面積密度マップを用いて前記電子線の前記被加工基板内での後方散乱の影響を考慮することによって電子線露光・現像後のパターン形状の変化を予測してパターン変形を施した転写用マスクを用いて、前記サブフィールドを露光する時に想定した前記サブフィールドの偏向位置と電流密度を概ね一致させて露光する工程とを有することを特長とするパターン形成方法。
【請求項２】
請求項１記載のパターン形成方法において、
前記パターン変形の変形量は電子線露光によりレジスト中に形成された潜像の現像工程あるいは現像後のレジストパターンを用いて基板加工を行なう工程での寸法変化量を加味したパターン変形量を前記第１、第２及び第３の面積密度マップとはさらに異なる第４の面積密度マップを用いて予め予測して修正したものであることを特徴とするパターン形成方法。
【請求項３】
請求項１記載のパターン形成方法において、
前記第２の区分領域は一つの前記サブフィールド全体であって、前記転写用マスクのパターン形成に用いられるパターンデータは、前記電子線のボケを前記一つのサブフィールド全体の電流に依存したクーロン効果による前記電子線のボケを前記一つのサブフィールドの面積密度と前記一つのサブフィールドの位置での電子光学系の収差とより予測してその形状を変形したものを含むことを特徴とするパターン形成方法。
【請求項４】
請求項３記載のパターン形成方法において、
前記一つのサブフィールドの位置での前記電子線のボケを前記一つのサブフィールドの偏向中心からの距離と前記一つのサブフィールドを透過する電流値の関数で表された近似式を用いて求めることを特徴とするパターン形成方法。
【請求項５】
請求項１記載のパターン形成方法において、
前記パターンデータは、前記サブフィールド内の局所的な電子線のボケの分布を前記第２の面積密度マップを用いて予測しその形状を修正したパターンであることを特徴とするパターン形成方法。
【請求項６】
請求項１記載のパターン形成方法における面積密度マップを作成する工程において、前記第２の面積密度マップの各前記第２の区分領域内の要素は考慮対象の第２の区分領域の周辺の第２の区分領域が当該第２の区分領域に与える影響を取り込んで求められた要素であることを特徴とするパターン形成方法。
【請求項７】
請求項１記載のパターン形成方法において、パターン形状を修正したことによる前記第１の面積密度マップの前記第１の区分領域内のパターン面積密度の変化量が予め定められた値より大きい場合、修正後のパターン形状により作成された面積密度マップに基づきパターン形状の修正を行なうことを特徴とするパターン形成方法。
【請求項８】
請求項７記載のパターン形成方法において、パターン形状を修正した後パターン形状の修正に伴う各パターンの面積変化量を算出し該面積変化量より前記第１の面積密度マップの前記第１の区分領域内の各要素の変化量を求め該要素の変化量が予め定められた値より大きい場合、該要素の変化量を前記第１の区分領域に加算することで修正面積密度マップを形成し該修正面積密度マップを基にパターン形状の修正を行なうことを特徴とするパターン形成方法。
【請求項９】
請求項７記載のパターン形成方法において、面積密度を求める各区分領域内の各パターンの周辺長を保存しパターン形状の修正量とパターンの周辺長よりパターン面積の変化量を計算することを特徴とするパターン形成方法。
【請求項１０】
請求項１記載のパターン形成方法において、パターンの修正量を前記第３の面積密度マップより計算される後方散乱のレベルと前記第２の面積密度マップあるいは複数の面積密度マップより計算される電子線のボケ量および修正対象パターン辺の属性とパターン変形量を関係付けられた補正テーブルまたは近似式を用いて決定することを特徴とするパターン形成方法。
【請求項１１】
請求項１記載のパターン形成方法において、前記修正対象パターン辺の属性は長辺長、短辺長またはパターン幅のいずれかを含むことを特徴とするパターン形成方法。
【請求項１２】
請求項１記載のパターン形成方法において、前記所定のパターンは複数のマスクを通過した電子線露光の重ね合せにより形成され、各々の前記マスクに対応したパターンのマスク別面積密度マップを形成する工程と、前記マスク別面積密度マップの同一位置の区分領域内要素を加算することにより合成面積密度マップを形成する工程とを含み、前記電子線のボケの算出には前記マスク別面積密度マップを用い、前記電子線の後方散乱の影響の算出には前記合成面積密度マップを用いることを特徴とするパターン形成方法。
【請求項１３】
請求項１記載のパターン形成方法において、前記所定のパターンは複数のマスクを通過した電子線露光の重ね合せにより形成され、各々の前記マスクに対応したパターンのマスク別面積密度マップを形成する工程と、各々の前記マスクに対応したパターンが一体となったパターンを用いて一体面積密度マップを形成する工程とを含み、前記電子線のボケの算出には前記マスク別面積密度マップを用い、前記電子線の後方散乱の影響の算出には前記一体面積密度マップを用いることを特徴とするパターン形成方法。
【請求項１４】
半導体回路の１層の回路パターンの全領域内に第１の大きさを有する複数個の第１の区分領域を設定し、該区分領域内のパターン面積を計算し第１の面積密度マップを作成する工程と、前記第１の区分領域とは異なる第２の区分領域を設定し第２の区分領域内のパターン面積を既に求められた第１の面積密度マップより計算し第２の面積密度マップを作成する工程と、前記第２の区分領域とは異なる第３の区分領域を設定し第３の区分領域内のパターン面積を既に求められた第１の面積密度マップより計算し第３の面積密度マップを作成する工程と、前記第２の面積密度マップを用いて露光用電子線の後方散乱の影響を考慮しかつ第３の面積密度マップを用いて前記露光用電子線のボケの影響を考慮した時のパターン変形量を予め予測してパターン形状を修正する工程とを含む工程により作成したパターンデータを用いて形成した電子線転写用マスクを用いて、レジストの塗布された半導体基板に電子線を照射し所望の回路パターンを形成することを特徴とする半導体装置の製造方法。
【請求項１５】
請求項１４記載の半導体装置の製造方法において、前記露光用電子線を１度に露光できる領域をサブフィールドとした時、１層の回路パターンをサブフィールド単位に分割し１サブフィールド内のパターンを複数の相補サブフィールドに分割する工程とを含み前記露光用電子線の後方散乱の影響は１サブフィールドおよび該サブフィールドに近接するサブフィールド内の全体のパターンを考慮した面積密度マップを使用し、前記露光用電子線のボケは分割された相補サブフィールド単位に作成された面積密度マップを使用して算出したものであって、電子線転写用マスクはサブフィールド単位および相補パターン単位に異なる位置に形成されており、前記露光用電子線の照射に当たっては同一位置の相補サブフィールドを前記半導体基板上の同一位置に重ねて露光し所望の回路パターンを形成することを特徴とする半導体装置の製造方法。
【請求項１６】
請求項１５記載の半導体装置の製造方法において、半導体装置のパターンデータの設計段階で使用された階層構造の少なくとも１部を保持した状態で該パターンデータを前記相補パターンへの分割を行った後に、前記サブフィールドへの分割および近接効果補正のための前記パターン変形を行なうことを特徴とする半導体装置の製造方法。
【請求項１７】
請求項１４記載の半導体装置の製造方法において電子線を１度に露光できる領域をサブフィールドとした時、１層の回路パターンを複数の相補パターンに分割し、各々の相補パターンを同一分割位置を有するサブフィールド単位に分割する工程とを含み電子線の後方散乱の影響は考慮対象サブフィールドおよび該サブフィールドに近接するサブフィールド内の全相補パターンを考慮した面積密度マップを使用し、電子線のボケは分割された相補サブフィールド単位に作成された面積密度マップを使用して算出したものであって、電子線転写用マスクはサブフィールド単位及び相補パターン単位に異なる位置に形成されており、電子線の照射に当たっては同一位置の相補サブフィールドを前記半導体基板上の同一位置に重ねて露光し所望の回路パターンを形成することを特徴とする半導体装置の製造方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は転写用マスクと電子線を用いた微細なパターン形成方法とそれを用いた半導体装置の製造方法に関する。
【０００２】
【従来の技術】
電子線によるパターン形成法ではレジスト膜が表面に形成された基板に照射された電子線が基板内で散乱し、広い範囲に渡って再びパターンを形成すべきレジスト膜の方に戻ってくるため、電子線の堆積エネルギが近傍にあるパターンの寸法や配置によって不均一な分布となる現象がある。この結果、仕上がりレジスト寸法は設計寸法からずれる。これは近接効果と呼ばれている。近接効果の及ぶ範囲は電子線の加速電圧によって異なるが例えば１００ｋＶの加速電圧で照射した電子はシリコン基板内で照射点から半径３０μｍ以上の範囲に渡って影響を与え、パターン寸法を変化させる。従って高精度な露光を行なうためには近接効果を適切に補正する必要がある。
【０００３】
また、電子線描画法は小パターンを逐次露光していくためスループットが低いという問題があったが、最近大面積転写用マスクを用いスループットを大きく改善した方式の実用化が進められている。この方式は一般的に電子線転写リソグラフィ（Ｅｌｅｃｔｒｏｎ  Ｐｒｏｊｅｃｔｉｏｎ  Ｌｉｔｈｏｇｒａｐｈｙ：以下ＥＰＬと略す）と呼ばれている。このＥＰＬでは一括で露光する範囲をサブフィールドと呼ぶが、サブフィールドは約２５０μｍ角である。近接効果の及ぶ範囲は前述の例では半径３０μｍ程度であるのでパターンの配置によってサブフィールド内での近接効果は不均一に現れる。このため従来転写用マスクを用いない電子線直接描画法で一般的であった露光パターン毎の露光量を変化させて最適なパターン寸法を得ようとする方式を使用することができない。
【０００４】
ＥＰＬ方式での近接効果を補正する手法の１つは補助露光法である。この手法は露光パターンの不均一性に基づく電子線の後方散乱による堆積エネルギ分布の不均一性を補助露光によって均一にする方式である。露光パターンが疎な領域で後方散乱による影響が少なくなるため疎な領域に補助露光を行なう。より正確に言えば露光パターンと白黒反転したパターンを電子線のボケ量をおよそ後方散乱の範囲に等しくし適切な露光量で露光することで達成される。この手法は近接効果のための複雑な計算を必要としないところが最大の利点である。しかしながら本来の露光パターンの露光の他に補助露光を実行する必要があるため、補助露光用のマスクが必要となり、また２度露光のため実質的なスループットが低下するという問題があった。さらに補助露光は露光パターン以外の領域に露光を行なうことになるためパターン部と非パターン部の堆積エネルギコントラストの低下を招くといった問題もあった。
【０００５】
ＥＰＬ方式での異なる近接効果補正方式としてパターン変形法がある。これは予め近接効果によって露光後にパターンが変形するのを予想してマスクパターンを変形させておくものである。露光後近接効果によって寸法が変化した結果仕上がり寸法が設計寸法になるようにマスクパターンの変形量を設定するものである。この方式では一旦適正なマスクを作成すれば、露光時は近接効果を考慮することなく露光を実行できるという利点がある。また２度露光が不要であるのでスループットを低下させるという問題もない。一方、近接効果によるパターン変形を予測したマスクパターンの変形量の計算は複雑であり長時間の計算機処理を要するという欠点があった。このため例えばジャーナル・オブ・ヴァキューム・サイエンス・アンド・テクノロジー（Ｊｏｕｒｎａｌ  ｏｆ  Ｖａｃｕｕｍ  Ｓｃｉｅｎｃｅ  ａｎｄ  Ｔｅｃｈｎｏｌｏｇｙ  Ｖｏｌ．Ｂ９  ｐｐ．３０４８−３０５３  Ｎｏ．６  Ｎｏｖ／Ｄｅｃ  １９９１）  に開示されているように同じ特徴を持つパターンをグループ化して図形変形量を決定するといういわゆるルールベースによる補正法が用いられた。また特開平３−２２５８１６号公報においてパターン面積密度マップを用いて露光量補正によって近接効果を補正する手法が開示されている。同様にパターン面積密度マップを用いて図形変形量を決定する方式がジャパニーズ・ジャーナル・オブ・アプライド・フィジックス（Ｊａｐａｎｅｓｅ  Ｊｏｕｒｎａｌ  ｏｆ  Ａｐｐｌｉｅｄ  Ｐｈｙｓｉｃｓ  Ｖｏｌ．３７  ｐｐ．６７６７−６７７３  Ｎｏ．１２Ｂ  Ｄｅｃ．  １９９８）に示されている。
【０００６】
【発明が解決しようとする課題】
しかしながら、上記従来法には以下のような欠点があった。先ず第１に大面積を一括で露光する方式では多数の電子が電子光学系に滞在するため空間電荷による電子線のボケ量が変化する問題がある。これは一般にクーロン効果と呼ばれている。前述したサブフィールド単位でトータル電流が露光パターンによって異なるためクーロン効果によってサブフィールド単位でボケの量が異なる。更に、実際にはクーロン効果はサブフィールド内であっても必ずしも一定ではなくパターンの粗密によって局所的にボケ量が変化している（ここではこれを局所クーロン効果と呼ぶ）。従来方式では局所クーロン効果まで考慮した近接効果補正は行われていなかった。
【０００７】
第２に、近接効果の補正には光学系の収差とクーロン効果で決まる電子線のボケ量と後方散乱を考慮してパターン変形量を決定する必要があるがパターン変形によって他のパターンにも影響を与えることになるため、一般的には連立方程式の解として相互に矛盾のない変形量を決定する必要がある。これには今日の複雑なＬＳＩパターンにおいては多大な計算時間を必要とする問題があった。
【０００８】
また、連立方程式を解く代わりにパターン変形に伴うパターン相互の影響を繰り返し演算によって収束させる手段が考えられる。この場合パターン変形量は面積密度と電子線のボケ量の関数で表されるため１元方程式を解くことで求められる。しかし方程式の解として変形量を決定したパターン変形と面積計算を繰り返すことになるためやはり計算時間がかかるのは避けられない。
【０００９】
第３に、相補分割、サブフィールド分割、近接効果補正の演算は複雑なパターン演算を膨大なＬＳＩデータに対して行なう必要がある。通常、ＬＳＩのレイアウトデータは階層的なデータとして処理されることによってそのデータ量は数１０分の１以下に圧縮されている。しかし上記従来法ではサブフィールド毎にフラット化されたデータに対してこれらの処理を行なうため、上記階層化の効果を活用することができず処理時間およびデータハンドリング上で大きな問題があった。
【００１０】
本発明の目的はＥＰＬ方式の電子線露光法において、電子線のボケを含めた前方散乱、後方散乱、クーロン効果のいずれにも対応可能な近接効果補正方法とこれを用いたパターン形成方法および半導体装置の製造方法を提供することにある。
【００１１】
本発明の他の目的はＥＰＬ方式の電子線露光法において複雑なＬＳＩパターンにおいても高速で精密な近接効果補正を可能とした処理によるパターン形成方法を提供することにある。
【００１２】
本発明の他の目的はＥＰＬ方式の電子線露光法において高速なデータ処理と圧縮されたデータの扱いによって半導体装置を低コストで製造する方法を提供することにある。
【００１３】
【課題を解決するための手段】
本願明細書に開示される発明のうち、代表的なものの概要は以下の通りである。
【００１４】
露光パターンを後方散乱、クーロン効果、プロセス起因の寸法誤差のそれぞれの影響範囲を考慮した区分領域に分割し、区分領域内のパターン面積占有率（パターン面積密度）を保存しておき、寸法補正に必要なパターン変形量を電子線の偏向位置と電子線の電流密度で定まるボケ量とそれぞれのパターン面積密度の関数としパターン寸法・属性に応じて決定し、電子線露光に際しては想定した電子線の偏向位置と電流密度と概ね一致させて露光する。ここでサブフィールド全体あるいは露光パターン全体の区分領域内のパターン面積密度を配列状データとして保存したものを面積密度マップと呼ぶ。上記以外にレジストの現像過程や後のエッチング工程において局所的なパターンの粗密依存性が発生する場合にはその影響範囲に応じた区分領域を設定しその区分領域内のパターン面積密度の関数としてパターン変形量を決定する。これら異なる大きさの複数の区分領域内のパターン面積密度は基準となる区分領域（最も小さい区分領域が望ましい）を設定しておきパターンデータから面積を計算する時は基準区分領域単位で行い、他の大きさの区分領域の面積密度は既に求められた基準領域内の面積密度から計算することが望ましい。この理由は最小の区分領域を持つ面積密度マップから、より大きいサイズの区分領域内の面積密度マップを作成することは単純な四則演算で可能であるため直接パターンデータから面積を計算するよりもはるかに高速に計算することができるためである。
【００１５】
近接効果やクーロン効果あるいはエッチングなどのプロセスの影響は特定の領域内のみのパターンの影響を受けるものと距離の関数として他の領域のパターンの影響を受けるものとがある。例えばクーロン効果はサブフィールド内のみのパターンの影響を受けるものであるが、近接効果は特定の区分領域内だけでなく周辺の領域からの影響を距離の関数として受け、異なるサブフィールドのパターンの影響も受けることになる。これら距離依存の影響は前記基準区分領域で求められたパターン面積密度を、例えば近接効果であれば堆積エネルギの分布関数で畳み込み積分、またはこれと同様の処理を行なうことで取り込むことができる。
【００１６】
ＥＰＬ方式の電子線露光方式ではいくつかの理由によりウエハ上の同一領域に複数のマスクを通過した電子線を重ね合わせて所望のパターンを露光することがある。その理由の１つは、転写用マスクとして、電子線が通過できる所望の形状の開孔部を有するステンシルマスクを用いた場合、１枚のマスクでは実現できないドーナッツ状のパターンや機械的強度が弱いため２枚に分割する必要のあるパターンが存在することである。他の理由はパターン密度が非常に高いパターンデータにおいてはクーロン効果が強く影響し、所望の解像性が得られない場合があり、２枚のマスクに分割して露光する必要性が出てくるからである。これら２枚（あるいは必要によりそれ以上）のマスクは相補マスクと呼ばれている。相補マスクを用いた露光において、クーロン効果は個別のマスクによる露光で影響を受けるが、近接効果は相補マスク全体のパターンの影響を受けることになる。この場合には相補マスク毎の面積密度マップとそれらを合成した合成面積密度マップを作成してそれぞれの影響を考慮してパターン変形を行なうことにより高精度なパターンを得ることが可能となる。
【００１７】
また、パターン変形に伴う各区分領域のパターン面積密度の変化を前記基準区分領域内で再計算し、パターン面積密度の変化量が予め定められた値以下になるまで繰り返す。このとき、面積計算を行なう区分領域内に各パターンの周辺長を保存しておくことにより、パターン変形量（各辺のシフト量）と周辺長からパターン面積の変化量が計算できる。変形後のパターンから面積を計算しその変化量を求めることは時間のかかる計算となるが、周辺長より面積の変化量を計算することは数値計算のみで実行できるため繰り返し演算の大幅な高速化が実現できる。
【００１８】
また、階層的に設計されたレイアウトデータに対して、できるだけ設計階層構造を保持したまま相補分割または相補分割と近接効果補正を行なう。このためには、サブフィールド分割は相補分割を行った後に行なうことが望ましい。
【００１９】
【発明の実施の形態】
（実施例１）
以下、本発明の実施例を図１を用いて説明する。図１は本発明のパターン形成方法の流れを示す図である。四角で囲まれた部分はデータ処理の機能ブロックである。先ずＬＳＩパターンデータ１０１が入力として与えられる。次に１０２で示すように電子線露光装置で定まるサブフィールド（電子線を一度に転写できる最大の大きさを持つ領域でウエハ上換算でおよそ２５０μｍ角）単位にパターンデータを分割する。またＬＳＩの層によっては必要に応じて相補パターンに分割を行なう。このサブフィールド分割と相補パターン分割の処理順序は逆であっても構わない。次にパターンデータを面積密度を計算する単位となる基準区分領域に分割し（１０３）、それぞれの区分領域内の個々のパターンの面積を計算し積算する（１０４）。区分領域の面積に占めるパターン面積の総和の比率を面積密度と呼ぶことにする。ＬＳＩパターン全体あるいはサブフィールド単位で処理を行なう場合はサブフィールド全体でそれぞれの基準区分領域について面積密度を計算しこれを２次元の配列データとしたものを基準面積密度マップと呼び保存しておく（１０５）。基準区分領域の大きさは１〜１０μｍ程度が望ましく、後にパターン変形のために用いるいくつかの面積密度マップの区分領域より小さいことが必要である。また、この基準区分領域は後述する局所クーロン効果、近接効果、プロセス要因の補正のための区分領域の中で最も小さい区分領域でもって代用することも可能である。
【００２０】
初めて基準面積密度マップが作成された場合には１０６で示される判定処理は無条件でＮの処理に進む。次に、パターンの寸法精度劣化をもたらすいくつかの要因毎に区分領域の大きさの異なるいくつかの面積密度マップを作成する（１０７）。本実施例ではサブフィールド（ＳＦ）全体のクーロン効果（１０８）。局所クーロン効果（１０９）、近接効果（１１０）、エッチング等のプロセス要因（１１１）の４種類の面積密度マップを作成する。それぞれの区分領域の大きさは、１０８で２５０μｍ角、１０９で２０μｍ角、１１０で１０μｍ角、１１１で３０μｍ角とした。これら４つの区分領域内の面積密度は前述した基準面積密度マップの要素の単純な加算平均処理だけで計算することができるため高速な処理が可能となる。１０８のＳＦクーロン効果は個別のサブフィールド内のパターン面積率と電子線の電流密度のみでクーロン効果が決まるため基準面積密度マップより計算された面積密度マップをそのままで保存する。一方、その他の３種の要因は考慮対象の区分領域内のパターンだけでなく周辺の区分領域のパターンであっても考慮対象のパターン寸法に対して距離に依存した影響を与える。そこでそれぞれの影響の程度を表す分布関数を用いて周辺区分領域パターンの影響取り込み計算を行った結果をそれぞれの面積密度マップとして保存する。寸法補正に必要なパターン変形量は電子線のボケ量にも大きく依存するため、電子線のボケ量は電子光学系で定まるサブフィールド偏向位置での収差とクーロン効果を考慮して決定する。上述した複数の面積密度マップの各要素とパターンの変形量を関係付けられた補正テーブルあるいは近似式１１２を用いて個々のパターンの変形を行なう（１１３）。
【００２１】
上述した処理によってパターン変形を行った結果、各区分領域内のパターン面積密度は変化し、従ってパターン変形量に誤差が発生することになる。これを補正するため変形後のパターンを用いて再び基準区分領域内のパターン面積を計算し（１０４）、面積密度マップ１０５を再度作成する。次に、最も新しい面積密度マップの各要素と１つ前に作成した面積密度マップの各要素との差（変化量）が予め定めた規定値より大きいか否かを判定し（１０６）、大きければ上記処理を繰り返す。もし規定値以下であればパターン変形の処理は終了であり、次工程であるマスク用パターンデータの作成１１４に進む。ここでは電子光学系の縮小倍率（例えば１／４）に合わせたパターンデータの拡大処理、鏡面反転処理等のパターン処理を行ないマスク描画機のデータフォーマットに合わせたデータ変換処理を行なう。次に、前記マスク描画データを用いてマスク描画・レジスト現像の後、基板エッチング等のマスク加工処理を経てＥＰＬマスクを作製する（１１５）。上記ＥＰＬマスクを用いて電子線露光を実行する際にはサブフィールド偏向位置と電子線の電流密度をパターン変形を実施した時に想定したサブフィールド偏向位置と電流密度に概ね一致させて行ない（１１６）所望のパターン形成を行なう。
【００２２】
（実施例２）
図２及び図３は相補データに分割を行なう場合の処理の詳細を示したものである。図２は面積密度マップを求めるまでの手法を示したものである。チップサイズ２０ｍｍ角のＬＳＩのレイアウト設計が完了したＣＡＤデータ２０１は最初に２５０μｍ角のサブフィールド２０２に分割される。分割後のデータはＳＦ１，１〜ＳＦ８０，８０の６４００個のデータとなる。次にそれぞれのサブフィールドを２個の相補データ２０３と２０４、例えばＳＦ１，１ではＳＦ１，１ＡとＳＦ１，１Ｂに分割する。これはそれぞれの相補データ内ではドーナッツ状のパターンを含まないように、またパターンの長手寸法／短手寸法（アスペクト比）が３０を越えない大きさに分割している。ここまでの処理を１つの機能ブロックとして中間データファイル２０５に出力した。中間データファイルの形式は以後の処理に用いる機能ブロックと整合のとれるものが良い。一般的にＬＳＩデータの授受に用いられるストリームファイル形式であると異なるソフトウエアを用いた処理にも汎用性が高い。ただし、上記の処理と以後の処理を連続して行なうにはソフトウエア固有のデータ形式で保存することが処理速度向上の上で望ましい。
【００２３】
次に上記の処理でサブフィールド分割および相補データ分割されたデータは近接効果等補正処理（クーロン効果、近接効果、プロセスに起因する寸法誤差補正処理をここでは近接効果等補正処理と略して呼ぶことにする）のためのコンピュータＰＣ０（２０６）に順次入力される。ＰＣ０にはＰＣ１からＰＣ１０の１０台のコンピュータ（２０７）が接続されており入力されたデータをサブフィールド単位で個別のコンピュータに分配する。個別コンピュータはサブフィールド毎に相補分割されたデータを５μｍ角の基準区分領域に分割しそれぞれの基準区分領域内のパターン面積密度を計算し基準面積密度マップｓｆｉ，ｊ＿ａ、ｓｆｉ，ｊ＿ｂ（２０８、２０９）を作成する。次に局所クーロン効果の補正のため２５μｍ角の区分領域Ｃ１を設定する。局所クーロン効果の影響が距離ｘ離れた点で関数ｆ１（ｘ）としてＣ１中心での局所クーロン効果面積密度マップの要素ｓｆｉ，ｊ＿ｃ１ａ（Ｋ，Ｌ）、ｓｆｉ，ｊ＿ｃｌｂ（Ｋ，Ｌ）は（式１）で求められる。
【００２４】
【式１】
0003686367.tif 000002
また同様にしてサブフィールド全体でのクーロン効果は２５０μｍ角を１つの区分領域としてパターン面積密度を計算すればよいので（式２）で求められる。
ここでｎはサブフィールド内の基準区分領域のｘ方向およびｙ方向個数である。個別のコンピュータはサブフィールド単位を扱うのでｓｆｉ，ｊ−ｃ２ａ、ｓｆｉ，ｊ−ｃ２ｂはそれぞれ１つの値を持つだけとなる。
【００２５】
【式２】
0003686367.tif 000003
これらのデータをコンピュータＰＣ０（２０６）に返しＬＳＩチップ全体の大きさの基準面積密度マップＭＡＰ＿Ａ、ＭＡＰ＿Ｂ（２１０、２１１）を作成する。同様に局所クーロン効果補正用面積密度マップＭＡＰ＿Ｃ１Ａ、ＭＡＰ＿Ｃ１Ｂ（２１２、２１３）とサブフィールドクーロン効果補正用面積密度マップＭＡＰ＿Ｃ２Ａ、ＭＡＰ＿Ｃ２Ｂ（２１４、２１５）を作成する。
【００２６】
一方後方散乱はサブフィールドを越えて影響を与えるためＬＳＩチップ全体で考慮する必要がある。このため１０μｍ角の区分領域を有する後方散乱補正用面積密度マップＭＡＰ＿ＰＲＯＸ（２１６）を基準面積密度マップＭＡＰ＿Ａ、ＭＡＰ＿Ｂ（２１０、２１１）より作成する。ＭＡＰ＿ＰＲＯＸの各区分領域内の要素は先ず基準面積密度マップＭＡＰ＿ＡとＭＡＰ＿ＢとでのＬＳＩチップ上同じ位置の各要素を（式３）のように加算する。
【００２７】
【式３】
0003686367.tif 000004
次に周辺区分領域からの影響を取り込む。後方散乱の分布関数をｆ２（ｘ）とすれば後方散乱の面積密度マップの要素は（式４）で求められる。
【００２８】
【式４】
0003686367.tif 000005
ここでｆ２（ｘ）はガウス分布で良く近似ができ通常（式５）のように記述される。
【００２９】
【式５】
0003686367.tif 000006
ここでβｂは点入射の電子線による後方散乱のエネルギ堆積量が入射点の１／ｅになる距離で後方散乱係数と呼ばれるものである。
【００３０】
次に作成された面積密度マップを用いてパターン寸法の変形を行なう手順を図３を用いて説明する。ここではプロセス起因の寸法誤差は無視できるとしてサブフィールドクーロン効果、局所クーロン効果、後方散乱による近接効果を補正する。従って図３（ａ）に示すように３種類の面積密度マップを使用する。また相補パターンの内の一方についてのみ説明するが他方についても同様の処理が適用できる。それぞれの面積密度マップの配列をサブフィールドクーロン効果マップでは（ｉ，ｊ）、局所クーロン効果マップでは（ｋ，ｌ）、後方散乱マップでは（ｍ，ｎ）と表記する。近接効果による寸法シフトは電子線のボケ（レジスト中での前方散乱を含む）と後方散乱によって発生する。従って本実施例の場合、区分領域のサイズが最も小さい後方散乱の面積密度マップの単位毎に電子線のボケ量と後方散乱の堆積エネルギ量を求める必要がある。電子線のボケの原因の１つは電子光学系で定まる収差とサブフィールド単位のクーロン効果である。光学系の収差はサブフィールドの偏向中心からの座標（ｘ，ｙ）で定まる。また電子線の電流密度をｉｂとすればサブフィールド単位の電流は当該サブフィールドの面積密度ｓｆｉ，ｊ＿ｃ２ａ  を用いてｉｂ＊ｓｆｉ，ｊ＿ｃ２ａで計算される。この時の電子線のボケ量が図３（ｂ）のようにβｆ１であったとする。さらに（ｉ，ｊ）サブフィールド内での局所クーロン効果は５０μｍ角内の面積密度ｓｆｉ，ｊ＿ｃｌａ（ｋ，ｌ）と電子線電流密度ｉｂよりｉｂ＊ｓｆｉ，ｊ＿ｃｌａ（ｋ，ｌ）の依存した値βｆ２として与えられる。ここで注意すべきは面積密度ｓｆｉ，ｊ＿ｃｌａ（ｋ，ｌ）は単に（ｋ，ｌ）区分領域内のみで決定した面積密度ではなく周辺区分領域の影響と取り込んだ面積密度であるという点である。また電子線のボケ量βｆ１、βｆ２はそれぞれの電流と光学系により定まる固有の値である。局所クーロン効果の１区分領域内での電子線のトータルのボケ量を、
【００３１】
【式６】
0003686367.tif 000007
の（式６）より求める。電子線のボケ量には０．１μｍ程度以下の領域で影響を及ぼすプロセス上の要因を加えることによってその影響を補正することも可能である。例えば化学増幅レジストの酸拡散の影響は電子線のボケとは独立の現象であるので（式６）に加算して等価的に電子線のボケが増大したものと見なすことができる。次に後方散乱の影響は（ｍ，ｎ）区分領域内の面積密度ＭＡＰ＿ＰＲＯＸｉ，ｊ（ｍ，ｎ）と電子線の基板材料からの反射の程度を表すηの積η＊ＭＡＰ＿ＰＲＯＸｉ，ｊ（ｍ，ｎ）に比例する。ここでηは正確には点入射の電子線が前方散乱でレジストに堆積するエネルギを１とした時に基板からの後方散乱によりレジストに堆積するエネルギの値として定義される。
【００３２】
図３（ｃ）に示したラインパターンとコンタクトホールについて寸法変形量は以下のようにして求められる。まず、パターン幅Ｗと寸法変化量Δｗを一般化の
【００３３】
【式７】
0003686367.tif 000008
ため電子線のボケ量βｆで正規化し、さらに式の簡単化のため後方散乱の当該区分領域での面積密度を単にαで表す。すなわち、
とすると、ラインパターンとコンタクトホールの場合の寸法変化量Δｗはそれぞれ以下（式８）より求められる。
【００３４】
【式８】
0003686367.tif 000009
ここでＥｓは基準エネルギで例えば１：１ライン／スペースパターンが設計値通りに仕上がる露光量でのパターンエッジでの堆積エネルギを表すものである。
【００３５】
それぞれの区分領域内のパターンについてパターン変形量を求め、パターン変形を実施する。異なる局所クーロン効果区分領域単位で電子線のボケ量が異なり、異なる後方散乱区分領域毎に後方散乱のレベルが異なる。従って同一の後方散乱区分領域内でのみ同一寸法パターンに同一のパターン変形量が与えられる。
【００３６】
以下、実施例１で説明したごとくパターン変形に伴う面積密度マップの変化を繰り返し計算で収束させ最終的なパターン変形量を決定する。
【００３７】
本実施例では面積密度の周辺領域の取り込み処理や寸法変化量の計算方法の一つの例を示したもので計算方法は同等のものであれば上述した式に限定するものではない。また、複数のコンピュータを使った並列計算での例を示したが、マルチＣＰＵを有するコンピュータを使った並列計算も同等に可能であり、処理データ量によっては１台のコンピュータのみで処理することが可能である。さらにまた、寸法変化量を決定するために（式８）の方程式の解を用いたが、これはパターン幅や面積密度に応じた補正テーブルを用意してテーブルを参照して必要な寸法変化量およびライン端のシフト量を求めることも可能である。
（実施例３）
図４はさらに他の実施例を示すものでパターン変形に伴う面積密度の変化による誤差を収束計算によって低減し、さらに高速に収束計算を実行する方式について説明したものである。本実施例では基準区分領域内のパターンを面積密度マップを用いて変形する工程のみを述べる。
【００３８】
図４に示す大きさａの基準区分領域４０１内にパターン４０２〜４０６の５つの矩形パターンが存在している。それぞれの矩形パターンの短辺をｗ１〜ｗ５とし長辺をｈ１〜ｈ５（図中に示していない）とする。但しパターン４０２と４０３は辺ｓ１を共有しており、パターン４０５は隣接する基準区分領域に跨って存在している。パターン４０５が基準区分領域の境界線で切られる線を仮想的な辺ｓ２として定義する。この時パターン変形を実施する前の基準区分領域内の面積密度α０は図中４０７で示す以下（式９）で与えられる。
【００３９】
【式９】
0003686367.tif 000010
ここで面積密度の計算においてパターン４０５は仮想辺ｓ２で切断されたものとして扱う。また各パターン毎の周辺長を４０８に保存する。この時、周辺長の計算においては共有辺ｓ１と仮想辺ｓ２は除外する。実際の面積密度の初期値α０と周辺長ｌｉの計算は図形演算によって行なう。次にパターン変形に伴う面積密度の変化を４０９で示す（式１０）により計算する。
【００４０】
【式１０】
0003686367.tif 000011
ただし、第１回目の計算では面積密度の変化量は０である。次に個々の基準区分領域で計算した面積密度マップを基に実施例２で示したと同様に電子線のボケ量の計算のための面積密度マップ１、４１０と後方散乱量の計算のための面積密度マップ２、４１１を作成する。これらの複数の面積密度マップより各区分領域内の個々のパターン変形量Δｗｉを求める（４１２）。この時、パターン変形に伴うパターン面積密度の変化は図中４１３のように（式１１）によって求めることができる。
【００４１】
【式１１】
0003686367.tif 000012
（式１１）では辺の移動に伴うパターンの頂点部分での僅かな面積の変化を無視しているが、より正確にはパターン毎に頂点の数も基準区分領域内で保存しておいてΔｗｉ２／ａ２＊（頂点数）を（式１１）に加算（太くする場合）、減算（細くする場合）することが望ましい。  パターンの変形量はパターン変形を行なう前の面積密度により決定しているため（式１１）で求めた面積密度の変化が大きい場合にはパターン変形量に誤差が生じることになる。このため図４で示した変形の妥当性の判定を（式１２）により行なう（４１４）。
【００４２】
【式１２】
0003686367.tif 000013
この結果でパターン面積密度の変化量が予め設定した許容値より大きい場合、すなわち、４１４において
（判定Ｎの場合）には面積密度の再計算（４０９）、面積密度マップの作成（４１０、４１１）、図形変形量の決定（４１２）および面積密度の変化量（４１３）を求めた上で４１４の判定を再度実行する。これを繰り返し４１４が許容値以下になった段階（判定Ｙの場合）でパターン変形量を適正値と決定し実際のパターン変形４１５を行なう。
【００４３】
一般的にはパターン変形量に伴うパターン面積の変化はパターン変形を実施した後のパターンに対して面積を計算し変形前後での面積差を計算する必要がある。この場合には各パターンの辺の移動処理を行ない、さらにその結果のパターンについて面積を計算する必要があるため処理時間がかかるという問題がある。本実施例では各パターンの周辺長を保存しておくことによって、元パターンの辺の移動処理を行なわずパターン変形に伴う面積密度の変化量を直接計算することが可能となった。このため単純な配列数値の計算のみでパターン変形量が適正であったかどうかの判定が可能である。パターン変形量が適正であることが判明した変形量に対して実際のパターンの辺の移動による変形を行なえばよいので時間のかかるパターン処理は１度だけですむという大きな利点がある。
（実施例４）
本実施例ではいくつかの計算処理をフィッティングより求めた近似式により高速化する手法を図５および図６を用いて説明する。
【００４４】
電子線のボケ量は電子光学系の収差とクーロン効果によるボケによって決まることが知られている。光学系の収差は電子線の電流値と偏向位置によって異なるため、個々のサブフィールドの電子線ボケ量はサブフィールド偏向位置、電子線の電流密度、サブフィールド内の開口面積によってそれぞれに異なることになる。この個々のサブフィールドの電子線ボケ量をシミュレーションによって求めることは多大な時間を要する。
【００４５】
本実施例の方法においては図５（ａ）に示すようにサブフィールド偏向可能な領域５０１に複数の代表点５０２を設ける。次にサブフィールド電流がＩ１である時、代表点Ｐ（ｘ，ｙ）での電子線のボケをシミュレーションにより求める。これを設定した代表点全てに対して行なう。またサブフィールド電流がＩ２である時も同様にして代表点での電子線ボケ量をシミュレーションによって求める。サブフィールド電流がＩ１である時の各代表点の位置（ｘ，ｙ）と電子線ボケ量βｆ０＿Ｉ１の関係を多項式によってフィッティングする。例えば５次の多項式によって近似すると（式１３）に示すような形式となる。
【００４６】
【式１３】
0003686367.tif 000014
これをサブフィールド電流がＩ２の場合にも同様に行ない電子線ボケ量βｆ０＿Ｉ２の近似式を求める。個別サブフィールド（（ｉ，ｊ）位置とする）のボケ量を求めるにはまず偏向位置（ｘ，ｙ）でのβｆ０＿Ｉ１またはβｆ０＿Ｉ２は（式１３）に座標を代入することによって求められる。（ｉ，ｊ）サブフィールドの面積密度は図５（ｂ）の面積密度マップ１（５０３）によって求められるのでこれをα１（ｉ，ｊ）とし、電子線の電流密度をｉｂとするとサブフィールド電流はα１（ｉ，ｊ）＊ｉｂとなる。電子線ボケ量が図５（ｃ）に示したようにサブフィールド電流に対して直線的であるとすれば（式１４）を用いて個別サブフィールドの電子線ボケ量を求めることができる。
【００４７】
【式１４】
0003686367.tif 000015
ここでは簡単のために（式１４）を１次多項式でフィッティングしたが、シミュレーションで求めるボケ量をさらに異なる電流値で行ない高次の多項式でフィッティングすることでより高精度な近似ができる。
【００４８】
次にサブフィールド内で分布を持つ局所クーロン効果についても同様に図５（ｄ）の面積密度マップ２（５０４）の面積密度α２（ｋ，ｌ）より
【００４９】
【式１５】
0003686367.tif 000016
として求めることができる。ただしここでの電子線ボケ量はサブフールド内電流分布による電子線ボケの変化量として予め求めておく。これによって電子光学系の収差、サブフィールドクーロン効果、サブフィールド内局所クーロン効果の全てを考慮した電子線ボケ量βｆは（式１６）によって求めることができる。
【００５０】
【式１６】
0003686367.tif 000017
（式１６）のβｆにより規格化したパターン寸法ｗと規格化したパターン寸法変化量Δｗを用いると（式８）の方程式は（式１７）のようになる。
【００５１】
【式１７】
0003686367.tif 000018
（式１７）よりΔｗを求めることは困難ではないがＬＳＩの個々のパターンについて実行すると多大な時間が必要となる。そこで本実施例では（式１７）の解を多項式で近似することを行なう。後方散乱による基板からの反射率であるηは基板材料が同じであれば基本的には一定であるので、ここではシリコン基板上に加速電圧１００ｋＶで電子線描画した場合を仮定しη＝０．６とする。規格化パターン寸法ｗをパラメータとして（式１７）をΔｗに対して解いた結果を図示したものが図６（ａ）である。面積密度αとパターン変形量Δｗの関係は次の３次の多項式（式１８）
【００５２】
【式１８】
0003686367.tif 000019
で非常によくフィッティングができる。また規格化パターン寸法が大きくなるとパターン寸法に依存せず１本の近似曲線に収斂することも分かる。さらに（式１８）３次式の各項の係数ａ３、ａ２、ａ１、ａ０と規格化パターン寸法ｗの関係を表したものが図６（ｂ）である。この結果もまたおよそ０．８＜ｗ＜３．２の範囲で３次の多項式によってよく表すことができる。３．２＜ｗの領域ではほぼ一定値である。すなわちパターン寸法に依存しないということが言える。
【００５３】
図６の結果を用いれば個別のパターン毎に方程式を解く必要はなく３次多項式に数値を代入するだけでパターン形状の変化量を求めることができ、処理の大幅な高速化が可能となる。なお本実施例ではラインパターンについて説明したがホールパターンや矩形パターンについても多項式近似によって同様な高速化が可能である。また、電子光学系の収差、サブフィールドクーロン効果、サブフィールド内局所クーロン効果の全てを考慮した電子線ボケ量を（式１６）の和の形で求めたが局所クーロン効果によるボケの寄与分の算出法によっては二乗和の形で求めることもできる。
（実施例５）
本発明を論理ＬＳＩのパターン露光に適用した例を図７を用いて説明する。図７はＬＳＩの第１配線層の製造工程に適用した場合の例である。図７（ａ）はＬＳＩパターンをサブフィールドに分割した時の一部分を示すもので７０１の領域が２５０μｍ角に分割された１つのサブフィールド領域である。図７（ａ）では実際のパターンそのものではなくパターン群を矩形として表してある。その中の一部の領域７０２および７０３の詳細パターンを図７（ｂ）に示す。この２つの領域ではほとんどのパターンが本ＬＳＩパターンの最小寸法である７０ｎｍで構成されている。パターン７０４は配線が７０ｎｍの等間隔ライン／スペースで配置された領域の１つのパターンを示す、またパターン７０５は疎なパターン領域での７０ｎｍのパターンを示す。
【００５４】
本実施例では後方散乱による近接効果補正のための区分領域を１０μｍとして面積密度マップを作成した（図には示していない）。またサブフィールド内の局所クーロン効果を補正するために５０μｍ角の区分領域を設定し、面積密度マップを作成した。この内サブフィールド７０１について面積密度の分布を３次元グラフで表現したものを図７（ｃ）に示す。図の見やすさのために図７（ａ）のＸ，Ｙ方向と図７（ｃ）では逆になっていることに注意が必要である。図７（ｃ）の面積密度マップの作成に当たっては隣接する区分領域内パターンの影響を取り込むために下記（式１９）のように隣接する９個の区分領域で重み付け平滑化処理を行なっている。
【００５５】
【式１９】
0003686367.tif 000020
ここで、α０、α１はそれぞれ平滑化前および平滑化後の面積密度である。また、ｇ（ｌ，ｍ）は重み係数で、
【００５６】
【式２０】
0003686367.tif 000021
の値を用いた。
【００５７】
近接効果、サブフィールドクーロン効果およびサブフィールド内局所クーロン効果を考慮に入れたマスク寸法変化量を図７（ｄ）に示す。この結果より設計寸法が７０ｎｍであるのでパターン７０４は５２．５ｎｍ、パターン７０５は１１０ｎｍである。４倍転写用マスクとして、それぞれ２１０ｎｍ、４４０ｎｍのパターンを形成し電子線露光を行なったところ図７（ｅ）の条件Ａの結果を得た。寸法測定の再現性±２ｎｍの範囲で７０ｎｍのパターンを得ることが可能であった。一方、図７（ｅ）の条件Ｂは近接効果とサブフィールドクーロン効果のみを考慮したものでサブフィールド内の局所クーロン効果を無視して寸法補正を行なった実験的な結果である。この場合にはパターン密度のやや高いパターン７０４で＋５ｎｍの寸法誤差が認められた。このように３種類の面積密度マップを使用してマスクパターン寸法を変化させることによってサブフィールド内の不均一なパターン配置による寸法誤差を高精度に補正することが可能となった。
【００５８】
本実施例では配線パターンのように比較的パターン密度の高い層について説明したが、論理ＬＳＩにおいてはスルーホール層（あるいはコンタクトホール層）では元々のパターン面積密度が低くサブフィールド内の局所クーロン効果の影響は軽微であるため近接効果とサブフィールドクーロン効果を補正する処理を行なえば十分である。また論理ＬＳＩのゲート層も同様にややパターン密度が低いためサブフィールド内の局所クーロン効果を無視することができる場合が多い。但しゲートパターンに要求される寸法精度は一般にＬＳＩの他の層より厳しいので要求精度によってはやはり局所クーロン効果の補正が必要となる。局所クーロン効果の補正の必要性は局所面積密度の許容最大値を設定しておきこの許容値を超える区分領域があるか否かで判断すればよい。また相補マスクに分割する必要がある場合には分割されたそれぞれのマスクパターン内での面積密度で判断する必要がある。
（実施例６）
図８はＬＳＩの設計段階で使用される階層的なデータ構造をできるだけ活かした形で近接効果等の補正データ処理を行う方式について説明したものである。図８（ａ）に示すＬＳＩデータの最上位階層のセルＡＢＣ（８０１）とし、その下位にＡＡ、ＢＢ、ＣＣ等がある。さらにＡＡ階層の下にはＡ、Ａの下にはａ、ｂ、ｃという階層構造を形成している。
【００５９】
この階層構造を保持したままでデータ処理するためには、サブフィールド分割はできるだけ処理の後工程で行なうことが望ましい。このため、本実施例では図８（ｂ）のようなデータ処理フローとした。まず相補パターン分割の必要のある場合にはＬＳＩデータＡＢＣ（８０１）を最初に２つの相補パターンＡＢＣ−１とＡＢＣ−２に分割する。本実施例では相補パターン分割のために汎用的なＬＳＩのデザインルールチェック（ＤＲＣ）用ソフトウエアを使用した。ＡＢＣ−１とＡＢＣ−２はＡＢＣのほとんどの階層構造が保持されている。
【００６０】
次に分割されたパターンＡＢＣ−１内のセル間のパターンに重なりがあるものについて一体化して重なりを除去する。その後１０μｍの基準区分領域を設定し区分領域内のパターン面積を計算し基準面積密度マップＭＡＰ−Ｂ１を作成する。またパターンＡＢＣ−２についても同様に処理して基準面積密度マップＭＡＰ−Ｂ２を作成する。基準区分領域の大きさを後方散乱補正用の区分領域の１０μｍと一致させたため、後方散乱面積密度マップＭＡＰ−ＰＲＸはＭＡＰ−Ｂ１とＭＡＰ−Ｂ２のマップの同一位置での区分領域の面積密度の値を単に加算するだけで作成が可能である。
【００６１】
次に、基準面積密度マップＭＡＰ−Ｂ１を元に区分領域が３０μｍの局所クーロン効果補正用面積密度マップＭＡＰ−ＣＬ１および区分領域が２５０μｍのサブフィールドクーロン効果補正用面積密度マップＭＡＰ−ＳＦ１を作成する。またもう一方の相補パターンついても同様にして面積密度マップＭＡＰ−ＣＬ２およびＭＡＰ−ＳＦ２を作成する。
【００６２】
次にパターンデータＡＢＣ−１（あるいはＡＢＣ−２）を２５０μｍ角のサブフィールド単位に分割する場合にはサブフィールド内にて階層が閉じているセルについてはそのままの形で階層を保存する。また、下位階層のセルがサブフィールド間に跨っている場合にはサブフィールド内で階層が閉じた複数のセルとして構成し直す。このようにサブフィールド内に階層を保持したパターンデータＡＢＣ−１−ＳＦ（あるいはＡＢＣ−２−ＳＦ）を作成する。これらのパターンデータと面積密度マップＭＡＰ−ＰＲＸ、ＭＡＰ−ＣＬ１、ＭＡＰ−ＳＦ１（あるいはＭＡＰ−ＣＬ２、ＭＡＰ−ＳＦ２）より補正演算処理８０２を行なうには実施例１から実施例５で述べられた手法が用いられる。補正演算の結果相補マスク用データＡＢＣ−Ｍ１とＡＢＣ−Ｍ２が作られる。
【００６３】
この段階ではほとんどの階層は展開されフラットなデータ構造となる。この理由は同一のセルであってもその配置位置によって寸法補正量が異なるためデータ圧縮ができなくなるためである。しかしながら、ＬＳＩの層によっては例えばホールパターンあるいはＬＳＩの論理ゲートの部分ではサブフィールド内の局所的なパターン密度が低い領域もありサブフィールド内の領域において一定の寸法補正量となり階層を維持することも可能である。この段階でパターンデータＡＢＣ−Ｍ１およびＡＢＣ−Ｍ２は直接マスク描画用電子線描画装置のフォーマットに出力すると効率がよい。８０６のデータを用いて作成された転写用相補マスクにより同一位置のサブフィールドはウエハ上の同一位置に転写されるように電子線露光８０７が実行される。
【００６４】
以上のように相補パターンデータへの分割をサブフィールド分割の前に行なうことでデータの階層構造をできるだけ保持したままデータ処理ができるので処理速度の高速化、処理データ量の低減に効果がある。もちろんサブフィールド分割後相補パターン分割を行なってもデータの階層を保持することは可能である。この時は元々のＬＳＩデータの構造によってはフラットな構造に展開するセルが多くなることがありうる。
【００６５】
上記のデータ処理においてサブフィールドサイズは２５０μｍ角であるので２０ｍｍ角のＬＳＩチップでは６４００個のサブフィールドデータが形成される。この時パターンデータＡＢＣ−１−ＳＦ（あるいはＡＢＣ−２−ＳＦ）は異なる６４００個のファイルとして保存してもよいが、全体を１つのファイルとしてサブフィールドデータを６４００個の下位階層のセルとしておくとデータの入出力の際に好都合である。
（実施例７）
次に、具体的な半導体集積回路装置の製造例を説明する。ここでは、例えばゲートアレイ、スタンダードセル等のようなセミカスタム方式で製造される半導体集積回路装置、半導体基板上にカスタムＩ／Ｏ（Ｉｎｐｕｔ／Ｏｕｔｐｕｔ）回路、カスタム論理回路またはＩ／Ｆ（インターフェイス）制御回路を有する半導体集積回路装置の製造方法に本発明を適用する場合を説明する。
【００６６】
図９（ａ）は、半導体集積回路装置における論理素子の一部を示す平面図である。ここでは２個のｎＭＩＳＱｎと、２個のｐＭＩＳＱｐとから構成されている。ｎＭＩＳＱｎは、半導体基板に形成されたｐ型ウエル領域ＰＷの表面のｎ型半導体領域９１１ｎ上に、ｐＭＩＳＱｐは、ｎ型ウエル領域ＮＷの表面のｐ型半導体領域９１１ｐ上に、それぞれ形成されている。ゲート電極９１２Ａは、ｎＭＩＳＱｎおよびｐＭＩＳＱｐに共有となっている。ゲート電極９１２Ａは、例えば低抵抗多結晶シリコンの単体膜、低抵抗多結晶シリコン膜の上部にシリサイド層を設けたポリサイド構造、低抵抗多結晶シリコン膜上に窒化タングステン等のようなバリア膜を介してタングステン等のような金属膜を堆積してなるポリメタル構造あるいは絶縁膜に掘られた溝内に窒化チタン等のようなバリア膜を堆積し、さらにその上に銅等のような金属膜を埋め込むことで形成されたダマシンゲート電極構造で構成されている。ゲート電極９１２Ａの下方の半導体基板部分はチャネル領域となる。
【００６７】
配線９１３Ａは、例えば高電位（例えば３．３Ｖまたは１．８Ｖ程度）側の電源配線であり、コンタクトホールＣＮＴを通じて２個のｐＭＩＳＱｐのｐ型半導体領域９１１ｐと電気的に接続されている。また、配線９１３Ｂは、例えば低電位（例えば０Ｖ程度）側の電源配線であり、コンタクトホールＣＮＴを通じて１個のｎＭＩＳＱｎのｎ型半導体領域９１１ｎと電気的に接続されている。配線９１３Ｃは、２入力ＮＡＮＤゲート回路の入力配線であり、コンタクトホールＣＮＴを通じてゲート電極９１２Ａの幅広部分で接触し電気的に接続されている。配線９１３Ｄは、コンタクトホールＣＮＴを通じてｎ型半導体領域９１１ｎおよびｐ型半導体領域９１１ｐの両方に電気的に接続されている。配線９１４Ａは、コンタクトホールＣＮＴを通じて配線９１３Ｄと電気的に接続されている。
【００６８】
ここでは、図９（ａ）の回路を製作する工程の内で本発明のパターン形成方法による電子線転写方式を用いた工程、光リソグラフィを用いた工程の使い分けを図９（ｂ）に示した。ｐ型あるいはｎ型のウエル領域の形成工程では微細な寸法を必要としないため通常の光リソグラフィによってパターン形成を行った。素子分離層の形成においては基準面積密度マップの他に、後方散乱マップ（ＭＡＰ１）、サブフィールド単位のクーロン効果マップ（ＭＡＰ２）およびサブフィールド内局所クーロン効果マップ（ＭＡＰ３）を作製してこれらの誤差補正を行なった２枚の相補マスクを使用してパターン形成を行なった。次にシリコン基板をエッチングにより溝形成を行なった後、シリコン酸化膜を埋め込みＣＭＰ（Ｃｈｅｍｉｃａｌ  Ｍｅｃｈａｎｉｃａｌ  Ｐｏｌｉｓｈｉｎｇ）による平坦化を行なった。次にゲート形成工程では上述のＭＡＰ１、ＭＡＰ２、ＭＡＰ３の他にドライエッチングでのローディング効果による寸法ばらつきを補正するためプロセスマップ（ＭＡＰ４）を用いてパターン形状の補正を行なった。またクーロン効果を低減するため２枚の相補マスクに分割し電子線露光時に重ね合わせることで所定のパターンを得ている。
【００６９】
コンタクト穴の形成工程では相補マスクへの分割は行なわずＭＡＰ１およびＭＡＰ２のみを使用したパターン変形を行なった。これは後の工程であるスルーホール１の形成工程でも同じである。金属配線１の工程では比較的パターン密度が高い領域があることと、長い配線もあるためステンシルマスクの機械強度確保のために２枚の相補マスクに分割している。また、補正のためのパターン変形にはＭＡＰ１、ＭＡＰ２およびＭＡＰ３を使用した。この後の工程は金属配線２、スルーホール２、金属配線３、スルーホール３、金属配線４、スルーホール４、金属配線５、スルーホール５、プロテクション工程があるが、これらはすべてゲート電極パターンに比べてパターン寸法が大きく、光リソグラフィによってパターン形成した。
【００７０】
また、上記は半導体装置の製造方法に関するものであったが、特にそれに制限されることはない。例えば液晶ディスプレイ等、エネルギ線を照射することによって感応基板上にパターン転写を行なうものであれば同様な手法を用いることによって応用可能である。
（実施例８）
本実施例では、マスク補正に用いるＤＡシステムの一例を図１０を用いて説明する。
【００７１】
まず、所望の半導体回路パターンを含む設計レイアウトデータを入力し、これを必要に応じて相補マスクパターンデータＰ１、Ｐ２へ分割する。なお、このとき併せて、分割辺図形ＰＤ、及びサブフィールド境界線を新たに図形データとして生成することが好ましい。但し、この段階ではサブフィールドへのデータの分割は行なわない。上記パターンデータに対して、前記実施例に述べたのとほぼ同様の方法により複数種類の区分領域に基づく面積密度マップ（ここでは後方散乱マップＤ１（ｘ，ｙ）、クーロン効果マップＤ２（ｘ，ｙ）、ウエハープロセスＤ３（ｘ，ｙ）を用いた）を生成する。ｘ，ｙはレイアウトパターンの座標である。
【００７２】
次に、上記設計パターンＰ、又は相補マスクパターンＰ１、Ｐ２の含むパターンの各辺を適宜分割し、補正対象辺及び補正点を抽出、番号付けする。補正点は、補正辺のほぼ中央に設定する。又、必要に応じて各補正対象辺の属性の場合分けを行なう。辺の属するパターンの寸法（電子線透過部の幅）、辺の長さ、隣接パターンまでの距離、辺のパターンにおける機能（長辺、ライン端、ホール等）等に応じて、全ての辺が分類される。これらの分類は、現在広く行われている光近接効果補正で用いられている図形演算による図形抽出を用いた。なお、必ずしも上記の全項目に対して場合分けを行なう必要はない。例えば、パターン寸法と辺の機能だけの分類で十分な精度が得られる場合もある。各補正対象辺（辺番号ｉ＝１，２，．．．）に対して、補正点座標、辺の長さ、辺の属性の対応表を記憶装置に記憶する。
【００７３】
上記第ｉ補正点における３種類の面積密度Ｄ１ｉ、Ｄ２ｉ、Ｄ３ｉを、上記面積密度マップＤ１（ｘ，ｙ）、Ｄ２（ｘ，ｙ）、Ｄ３（ｘ，ｙ）から内挿近似して算出する。ここで、Ｄ１ｉの値はパターンと第ｉ補正点座標（ｘｉ，ｙｉ）を中心とする第１の分布関数ｆ１（ｘ−ｘｉ，  ｙ−ｙｉ）の畳み込み積分の近似値、即ち、おおよそ上記補正点を中心とする半径ｒ１（後方散乱の範囲）内の後方散乱に寄与するパターンの有効面積密度を表している。ここに、関数ｆ１は、実施例２の式５に示したものと同様のものである。同様に、Ｄ２ｉ、Ｄ３ｉは、各々上記補正点を中心とする局所クーロン効果作用範囲、及びウエハエッチング時ローディング効果作用範囲内の面積密度を表す。
【００７４】
次に、各補正点における補正量、即ち補正辺の移動量を、Ｄ１ｉ、Ｄ２ｉ、Ｄ３ｉから決定する。ここでは、補正対象辺の属性及び、Ｄ１、Ｄ２、Ｄ３に応じて最適補正量ＰＢをあらかじめ計算し、補正ルールテーブルとしてあらかじめ準備しておき、各補正点毎に上記ルールテーブルを参照することにより補正量を決定した。ルールテーブルの生成方法については、後で述べる。又、補正に伴う面積変化量ｄＳｉを、補正辺毎に（対象辺の長さＬｉ  ｘ  補正量ＰＢｉ）により計算する。以上求めた補正量、さらに必要に応じて面積変化量と補正点（補正辺）との対応表を生成し記憶装置に記憶する。
【００７５】
補正に伴う面積又は面積密度の変化量が十分小さい場合には、前記パターンデータ上で各補正対象辺を上記補正量だけ移動させる。辺データの移動は通常の図形演算処理により行なう。辺の末端処理等については、通常の光近接効果補正で行われているのと同様にして、不要な凹凸が生じない様に適宜処理することが好ましい。
【００７６】
各補正辺移動による面積変化ｄＳｉにともなう各面積密度マップの変化分ｄＤｊ（ｘ，ｙ）  （ｊ＝１，２，３）  は次のようにして計算することができる。但し、ここでは上記変化分は連続関数として表現される。
ｄＤｊ（ｘ，ｙ）＝ΣｄＳｉ・ｆｊ（ｘ−ｘｉ，ｙ−ｙｉ）  （式２１）
（和は全ての、又は（ｘ，ｙ）近傍の補正点に対してとる。）
補正後の各補正点の面積密度は、
Ｄ’ｊｉ＝Ｄｊｉ＋ｄＤｊ（ｘｉ，ｙｉ）        （式２２）
により与えられる。補正に伴う面積又は面積密度の変化量が無視できない場合には、上式により各補正点に計算し直した複数の面積密度を用いて、再度前記補正ルールテーブルを参照して補正量を再設定する。以下上記課程を補正に伴う面積又は面積密度の変化量が十分小さくなるまで繰り返す。この様にして各補正対象辺の補正量を最終的に決定した後、図形演算処理により前記パターンデータ上で各補正対象辺を上記補正量だけ移動させる。
【００７７】
その後、サブフィールド分割を行った後、パターンデータをＥＢ描画データとして出力し、これを用いて電子線投影露光装置用ステンシルマスクを描画、製作し、上記ステンシルマスクを用いて半導体基板上のレジスト膜を電子線露光、現像して所望の半導体回路パターンを形成した。なお、上記ステンシルマスクを電子線描画により作成する際にもマスク描画用近接効果補正が必要である。ここでは、上記マスク描画における後方散乱近接効果補正に用いる面積密度マップもＤ１、Ｄ２、Ｄ３と併せて前記基準面積密度マップから作成し、これに基づいて描画を行った。これにより補正通りの寸法をもつマスクパターンを形成することができた。
【００７８】
本実施例によれば、パターンの属性に応じた最適補正量を図形演算によるパターン分類及び高速テーブル参照により極めて迅速に決定できる。又、補正対象辺の位置に応じた最適補正量を比較的空間的に粗い面積密度マップから高精度に内挿近似可能である。また、各辺の長さと補正量から面積密度の変化分が簡単に計算できるので、実際に図形変形を行うことなく最適補正量の収束計算が可能である。
【００７９】
次に、上記補正テーブルの決め方について、ここでは簡単のためラインパターンを例として（図１１を用いて）説明する。幅Ｗのラインパターンマスクを透過した電子線の透過率分布をＵ（ｘ）とする。ここではＵ（ｘ）として幅Ｗの単純なステップ関数を仮定する。電子ビームのパターンに依存しない（クーロン効果を除いた前方散乱による）ぼけの分布関数をＢｆ（ｘ）、クーロン効果による（パターン面積密度に依存する）ぼけの分布関数をＢｃ（ｘ）、後方散乱によるかぶり電子線強度をＶとすると、ウエハ上の電子線強度分布Ｉ（ｘ）は、およそ次式で与えられる。（但し（＊）は畳み込み積分）
0003686367.tif 000022
但し、ここで上記２つのぼけを各々広がりβｆ、βｃのガウス分布で仮定した。分布Ｉ（ｘ）は、Ｗ，  βｆ，  βｃ，  Ｖの関数であることに注意する。従って、ここで現像後のレジストパターンのエッジがＩ（ｘ）  ＝  Ｉｃｒｉｔとなるｘ＝Ｘｅｄｇｅで決まると仮定すると、ＸｅｄｇｅもＷ，  βｆ，  βｃ，Ｖの関数となる。簡単のため、β＝（βｆ２＋βｃ２）０．５とすると以下の通り。
【００８０】
Ｘｅｄｇｅ＝ｆ（Ｗ，β，Ｖ）              （式２４）
設計ラインエッジ位置をＸｄｅｓｉｇｎとすると、設計値からのズレはＸｅｄｇｅ−Ｘｄｅｓｉｇｎで表されるので、このずれを補正するのに必要な最適バイアス量（辺の移動量）ＰＢは次式で求められる。
【００８１】
ＰＢ  ＝‐（Ｘｅｄｇｅ−Ｘｄｅｓｉｇｎ）／ＭＥＦ    （式２５）
ここでＭＥＦはいわゆるマスクエラーファクターと呼ばれる量で、マスク寸法変化に対する実際のレジスト寸法変化の割合を表し、次式で求まる。（マスク上の寸法はウエハ上に縮小された時の寸法に換算する。）
0003686367.tif 000023
以上より、
0003686367.tif 000024
上記バイアスＰＢ  をあらかじめ様々なＷ、β、Ｖに対して計算し３次元テーブル化した。本実施例ではＭＥＦを上式により計算してバイアスＰＢを求めたが、簡単のため
ＰＢ＝−（Ｘｅｄｇｅ−Ｘｄｅｓｉｇｎ）            （式２８）
で近似してもよい。又、同様のテーブルは、ホールパターン、ライン端に対して作成してもよい。さらに、上記Ｉ（ｘ）に対する式２３中の  Ｕ（ｘ）  を様々に設定することにより、近傍に隣接パターンが存在する場合の効果をとりこむことができる。このようにして、前記補正点抽出において行った辺の属性分類に応じてテーブルの使い分けを行うことにより、極めて高精度の補正が可能となる。
（実施例９）
ここでは補正点毎に、電子線投影像を簡易シミュレーションして寸法変動量を求め、これによりバイアス補正量を算出する方法について説明する。パターンの微細化に伴い開口パターン間のスペースがクーロン効果によるビームぼけ量に近づくにつれ、クーロン効果に依存して近傍パターンの状況が寸法に影響を与えるようになり、この効果をより正確に考慮する必要が生じる。これらの影響は実施例８における補正テーブルの場合分け（特に辺の属性の場合分け）をより木目細かく行うことによっても対応できるが、補正テーブルが膨大になってしまうという問題点がある。
【００８２】
そこで、本実施例では、各補正点毎に該補正点近傍で電子線強度分布のシミュレーションを行ない、補正対象エッジの位置ズレを求めその最適補正量を計算した。電子線強度分布Ｉ（ｘ，ｙ）は、式２３のＵ（ｘ）を、マスクの開口パターンの２次元分布Ｕ（ｘ，ｙ）に置きかえることにより計算される。ここでβ及びＶの値はは実施例８同様にして複数の面積密度マップに基づいて補正点毎に計算した値、または該補正点の属する各マップの区分領域の値を用いる。
【００８３】
又、レジスト内の酸やクエンチャの拡散の影響や、現像等に関する各種の公知の計算モデルを加えてもよい。シミュレーションされた電子線強度分布から最適補正量を求める方法については、いわゆるモデルベースもしくはシミュレーションベースの光近接効果補正として広く知られる技術で用いられている手法を用いることができる。
【００８４】
補正のためにはＩ（ｘ，ｙ）を全領域で計算する必要はなく、例えば補正点（ｘｉ，ｙｉ）を中心とする近傍３点もしくは５点だけでも十分である。本実施例が、従来光近接効果補正と異なるのは、補正の元となる投影像の物理的意味が異なることに加え、投影像が複数の面積密度マップに依存し、同一パターンであっても露光領域内の位置により異なることにある。
【００８５】
以上のように、本発明では複数の誤差要因による影響をパターン面積密度マップを利用することで高速高精度に補正する方式について説明した。実施例は本発明のいくつかの例でありこれに制限されるものではない。パターンの寸法精度を劣化させる要因としては電子線転写時のみならず、マスクを製作する工程においても露光方式、基板構造によって誤差を発生させる可能性がある。
【００８６】
マスク描画において電子線描画を用いることが多いためこの段階でも近接効果の影響を受けることになる。その場合にはさらに異なる面積密度マップを作成することによって補正することが可能となる。この場合電子線転写時の縮小倍率を考慮した近接効果補正が必要なことは言うまでもない。
【００８７】
また電子線の露光照明系では光学方式と異なり蝿の目レンズによる照度均一化を行なうことが困難であるため、電子線の一括の露光範囲（サブフィールド）内での照度均一性が問題にあることがある。このような場合においてもサブフィールド内での不均一性を面積密度マップと同等の２次元マップを作成してパターン変形を行なうことができる。
【００８８】
【発明の効果】
以上説明したように、本発明の方法では露光パターンを、後方散乱、クーロン効果、プロセス起因の寸法誤差のそれぞれの影響範囲を考慮した区分領域に分割して複数の面積密度マップを作成し、区分領域内の各パターンの補正に必要な寸法変化量を面積密度の関数として決定することによって複数の効果を取り込んだマスクを作成することが可能となる。この結果、従来補正処理が困難であった複数の要因による寸法誤差を同時に補正することができる。また、面積密度マップを用いることで補正に必要なデータ処理時間を大幅に短縮でき高精度な露光パターンを得ることができる。
【００８９】
またパターン変形に伴うパターン面積の変化をパターンの周辺長より計算するすることによって面積密度マップの再計算を高速に行なうことができるためさらに処理の高速化が可能となる。
【００９０】
また、各種誤差要因による寸法誤差を補正する時に近似関数を用いることによって最適な寸法変化量を単純な数値代入による計算で求めることができるため高速な処理が可能となる。
【００９１】
さらに、ＬＳＩパターンの階層構造を保持したまま演算を行なうことによって扱うデータ量を大幅に削減できるため使用する計算機が能力の低いものでも可能となり、また高速な処理が実現できる。
【００９２】
これら処理の高速化は最適な寸法変化量を求める際の繰り返し演算において、より精密に実行することができるため、今後のＬＳＩで要求される厳しい寸法精度を満足することができＬＳＩの性能向上に大きく貢献できるものである。
【図面の簡単な説明】
【図１】本発明に係るパターン形成方法を説明するためのフロー図である。
【図２】本発明に係るパターン形成方法において、特にパターンデータを複数の相補パターンに分割し、面積密度マップを求める方法を説明するための図である。
【図３】本発明に係るパターン形成方法において、特に面積密度マップを用いてパターン寸法の変形を行なう方法を説明する図で（ａ）はサブフィールド、局所クーロン効果及び後方散乱を求めるための面積密度マップに対するそれぞれの区分領域サイズ及び位置表記法を示す図、（ｂ）は各位置におけるサブフィールド電流、局所クーロン効果、面積密度及びそれに対応する電子線のボケ量、後方散乱エネルギ量を示す図、（ｃ）はラインパターン及びコンタクトホールの寸法変形量を示す図である。
【図４】本発明に係るパターン形成方法において、特にパターン変形を行う方法を具体的なパターンで説明するための図である。
【図５】本発明に係るパターン形成方法において、特に電子線のボケ量を簡便に求める方法を説明するための図で、（ａ）はサブフィールド変更可能な領域を示す図、（ｂ）はクーロン効果に対する面積密度マップを示す概念図、（ｃ）はサブフィールド電流値に対するサブフィールドクーロン効果による電子線のボケを示す図、（ｄ）はサブフィールド内局所クーロン効果に対する面積密度マップを示す概念図、（ｅ）は電流値に対する局所的クーロン効果による電子線のボケを示す図である。
【図６】本発明に係るパターン形成方法において、特にパターン寸法の変化量を簡便に求める方法を説明するための図で、（ａ）はパターン面積密度に対する規格化したパターン変形量を示す図、（ｂ）は規格化パターン寸法に対する近似多項式の係数示す図である。
【図７】本発明に係るパターン形成方法をＬＳＩの第１配線層の工程に適用した例を示す図で、（ａ）ＬＳＩパターンをサブフィールドに分割したときの概略平面図、（ｂ）は具体的な配線パターンに例を示す平面図、（ｃ）サブフィールド内での面積密度分布を３次元グラフィックスで示す図、（ｄ）はパターン位置に対応するマスク寸法変化量を示す図、（ｅ）はパターン位置に対するレジストパターン寸法を示す図である。
【図８】本発明に係るパターン形成方法において、特に近接効果の補正データ処理を説明するための図で、（ａ）はＬＳＩパターンデータの構成を示す概念図、（ｂ）はデータ処理フローを示す図である。
【図９】本発明に係るパターン形成方法をＬＳＩ製造に適用した例を示す図で、（ａ）ＬＳＩパターンの概略平面図、（ｂ）は各種パタン形成工程における本願発明に係るステンシルマスクを適用したＥＰＬ及び光リソグラフィの使い分を説明するための図である。
【図１０】本発明に係るパターン形成方法において、特にマスク補正に用いるＤＡシステムを説明するためのフロー図である。
【図１１】本発明の別の実施例で用いる補正テーブルの作成方法の説明のための模式図。
【符号の説明】
１０１…ＬＳＩパターンデータの入力工程、１０２…サブフィールドおよび相補パターンに分割する工程、１０３…基準区分領域に分割する工程、１０４…区分領域内の個々のパターンの面積を計算し積算する工程、１０５…面積密度マップを作成する工程、１０６…面積密度の変化量が規定値以下か否かを判定する工程、１０７…複数の面積密度マップの全体、１０８…１サブフィールド（ＳＦ）内の全パターンの面積を計算して形成した面積密度マップ、１０９…サブフィールド内の局所的な領域毎にパターン面積を計算して形成した面積密度マップ、１１０…電子線の後方散乱の影響範囲でパターン面積を計算して形成した面積密度マップ、１１１…プロセス起因の誤差の影響範囲でパターン面積を計算して形成した面積密度マップ、１１２…複数の面積密度マップの各要素とパターンの変形量が関係付けられた補正テーブルあるいは近似式、１１３…補正のためのパターン変形工程、１１４…マスク用パターンデータの作成工程、１１５…ＥＰＬマスクを作製する工程、１１６…電子線露光工程、２０１…ＬＳＩのＣＡＤデータ、２０２…サブフィールドに分割されたＬＳＩデータ、２０３…２個の相補パターンに分割されたサブフィールドデータの一方、２０４…２個の相補パターンに分割されたサブフィールドデータの他方、２０５…データ授受のための中間データファイル、２０６…データ管理と複数のコンピュータ制御のためのコンピュータ、２０７…データを並列処理するための複数のコンピュータ、２０８…１サブフィールドの基準面積密度マップ（相補パターンの一方）、２０９…１サブフィールドの基準面積密度マップ（相補パターンの他方）、２１０…ＬＳＩチップ全体の基準面積密度マップ（相補パターンの一方）、２１１…ＬＳＩチップ全体の基準面積密度マップ（相補パターンの他方）、２１２…局所クーロン効果補正用面積密度マップ（相補パターンの一方）、２１３…局所クーロン効果補正用面積密度マップ（相補パターンの他方）、２１４…サブフィールドクーロン効果補正用面積密度マップ（相補パターンの一方）、２１５…サブフィールドクーロン効果補正用面積密度マップ（相補パターンの他方）、２１６…ＬＳＩチップ全体の後方散乱補正用面積密度マップ、４０１…基準区分領域、４０２〜４０６…矩形パターン、４０７…基準区分領域内の面積密度を求める工程、４０８…基準区分領域内の各パターンの周辺長を保存したテーブル、４０９…パターン変形後の基準区分領域内の面積密度を求める工程、４１０…ボケ量の計算のための面積密度マップ、４１１…後方散乱量の計算のための面積密度マップ、４１２…個々のパターン変形量を求める工程、４１３…パターン変形に伴うパターン面積密度の変化量を求める工程、４１４…パターン変形量の妥当性の判定を行なう工程、４１５…適正なパターン変形量によりパターン変形を行なう工程、５０１…サブフィールド偏向可能な領域、５０２…電子線のボケ量をシミュレーションにより求めるための代表点、５０３…サブフィールドクーロン効果補正用の面積密度マップ、５０４…サブフィールド内局所クーロン効果補正用の面積密度マップ、７０１…ＬＳＩパターン内の１サブフィールド領域、７０２…サブフィールド内の一部の領域およびその拡大図、７０３…サブフィールド内の他の一部の領域およびその拡大図、７０４…等間隔ライン／スペースの領域にある１つの７０ｎｍラインパターン、７０５…疎なパターン領域での７０ｎｍのラインパターン、８０１…ＬＳＩデータの最上位階層のセル、８０２…ＬＳＩパターンデータ、８０３…近接効果等の補正演算工程、８０５…近接効果等の補正処理後のＬＳＩパターンデータ（相補マスクの一方）、８０６…近接効果等の補正処理後のＬＳＩパターンデータ（相補マスクの他方）、８０７…２枚の相補マスクを用いた電子線露光工程、９１１ｎ…ｎ型半導体領域、９１１ｐ…ｐ型半導体領域、９１２Ａ〜９１２Ｂ…ゲート電極、９１３Ａ〜９１３Ｄ…金属配線１、ＡＢＣ−１…相補パターンに分割された一方のパターン、ＡＢＣ−２…相補パターンに分割された他方のパターン、ＭＡＰ−Ｂ１…基準面積密度マップ（相補パターンの一方）、ＭＡＰ−Ｂ２…基準面積密度マップ（相補パターンの他方）、ＭＡＰ−ＰＲＸ…後方散乱補正用面積密度マップ、ＭＡＰ−ＣＬ１…局所クーロン効果補正用面積密度マップ（相補パターンの一方）、ＭＡＰ−ＣＬ２…局所クーロン効果補正用面積密度マップ（相補パターンの他方）、ＡＢＣ−１−ＳＦ…サブフィールド単位に分割されたパターンデータ（相補パターンの一方）、ＡＢＣ−２−ＳＦ…サブフィールド単位に分割されたパターンデータ（相補パターンの他方）、ＭＡＰ−ＳＦ１…サブフィールドクーロン効果補正用面積密度マップ（相補パターンの一方）、ＭＡＰ−ＳＦ２…サブフィールドクーロン効果補正用面積密度マップ（相補パターンの他方）、ＡＢＣ−Ｍ１…相補パターンの一方のマスク製作用パターンデータ、ＡＢＣ−Ｍ２…相補パターンの他方のマスク製作用パターンデータ。
【図１】
0003686367.tif 000025
【図２】
0003686367.tif 000026
【図３】
0003686367.tif 000027
【図４】
0003686367.tif 000028
【図５】
0003686367.tif 000029
【図６】
0003686367.tif 000030
【図７】
0003686367.tif 000031
【図８】
0003686367.tif 000032
【図９】
0003686367.tif 000033
【図１０】
0003686367.tif 000034
【図１１】
0003686367.tif 000035