(57)【要約】【課題】  高速の処理速度つまりは計算スピードでかつ精度よくパターン分割後の各矩形からの蓄積エネルギーを求めることのできる描画方法の提供。【解決手段】  近接効果を補正した電子ビームの最適照射量を求めるために用いられる、試料上に描画するパターンを分割した各矩形からの蓄積エネルギーの算出方法に関し、予め複数の各要素毎のエネルギー分布からなる離散的な関数マップを求めておき、この関数マップに従って各矩形からの蓄積エネルギーを求めるようにした。これによれば、予め求めておいた関数マップつまりは各要素のエネルギー分布に従って簡単な計算を行うだけでよいことから、任意の点における各矩形からの蓄積エネルギーを高速につまりは時間をかけることなく求めることができるようになる。また、個々の要素のグリッド間隔を狭く設定した場合であっても、短時間で精度のよい矩形からの蓄積エネルギーを求めることができる。【選択図】        図３
2014029888.tif 000002
【特許請求の範囲】
【請求項１】
試料上に描画する回路パターンを分割する任意の矩形毎に、近接効果による前記パターンの寸法変化を低減するのに最適な電子ビームの照射量を求めるため、電子ビームの照射対象とされた第１の矩形を基準とする所定範囲内にある前記第１矩形以外の第２の矩形からの実質的な電子ビームの照射量である蓄積エネルギーを前記第２の矩形毎に求める描画方法であって、
所定位置に電子ビームを照射した際における複数の各要素毎のエネルギー分布からなる離散的な関数マップを求めるステップと、
前記求めた離散的な関数マップに従って前記第２の矩形からの蓄積エネルギーを求めるステップと
を備える描画方法。
【請求項２】
前記離散的な関数マップを求めるステップは、離散的な関数マップを、
2014029888.tif 000020
（ただし、ｓ；各要素間のグリッド間隔、Ｎ；Ｒ／ｓ（ＲはＰＳＦ定義域））
により求めることを特徴とする請求項１に記載の描画方法。
【請求項３】
前記第２の矩形からの蓄積エネルギーを求めるステップは、第２の矩形からの蓄積エネルギーを、
2014029888.tif 000021
（ただし、x；第１矩形上の任意の点のＸ座標位置、y；第１矩形上の任意の点のＹ座標位置、xi；関数マップ上の各要素のｘ座標位置、yj；関数マップ上の各要素のｙ座標位置、ｓ；各要素間のグリッド間隔、Ｎ；Ｒ／ｓ（ＲはＰＳＦ定義域））
に基づいて求めることを特徴とする請求項２に記載の描画方法。
【請求項４】
前記第２の矩形からの蓄積エネルギーを求めるステップは、前記第２の矩形に従って任意の大きさからなる複数の第３の矩形の組み合わせを特定するステップと、離散的な関数マップに従って前記特定した複数の第３の矩形毎に蓄積エネルギーを求めるステップと、該求めた第３の矩形からの蓄積エネルギーを前記組み合わせ態様に従って加減算するステップとを含んでなることを特徴とする請求項３に記載の描画方法。
【請求項５】
試料上に描画する回路パターンを分割する任意の矩形毎に、近接効果による前記パターンの寸法変化を低減するのに最適な電子ビームの照射量を求めるため、電子ビームの照射対象とされた第１の矩形を基準とする所定範囲内にある前記第１矩形以外の第２の矩形からの実質的な電子ビームの照射量である蓄積エネルギーを前記第２の矩形毎に求める描画装置であって、
所定位置に電子ビームを照射した際における複数の各要素毎のエネルギー分布からなる離散的な関数マップを求める手段と、
前記求めた離散的な関数マップに従って前記第２の矩形からの蓄積エネルギーを求める算出手段と
を備える描画装置。
【請求項６】
請求項１乃至４のいずれかに記載の描画方法における各ステップをコンピュータで実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電子ビームを用いてマスク基盤などの試料上に所望の回路パターンを描画するために、近接効果による回路パターンの寸法変化を低減するのに最適な電子ビームの照射量を求める描画方法及び描画装置並びにプログラムに関する。特に、最適な照射量を求めるために考慮する必要がある前記パターンを分割した各矩形からの蓄積エネルギーの算出に関し、高速にかつ精度よく各矩形からの蓄積エネルギーを求める技術に関する。
【背景技術】
【０００２】
近年、ＬＳＩなどの半導体集積回路のより一層の高集積化に伴い、半導体デバイスに要求される回路寸法及び回路線幅が年々微細化されてきている。これらの半導体デバイスへ所望の寸法及び線幅の回路パターンを形成するには高精度の原画パターン（レチクル又はマスクとも呼ぶ）が必要とされており、こうした原画パターンを生成するものとして、例えばレジスト膜が塗布された金属基盤（マスク基盤など）の試料上に電子ビーム（電子線）等を照射して原画パターンを描画する、所謂リソグラフィ技術を適用した描画装置などが従来から知られている。
【０００３】
ところで、マスク基盤等の試料上に電子ビームが照射されると、試料上に形成されるレジストパターンの寸法を変動させる近接効果と呼ばれる影響が現われる。これは、照射された電子がレジストや金属基盤などと衝突することにより生ずる前方散乱電子や後方散乱電子によって、試料上の意図していない箇所にもビームが照射されてしまい、結果として主に回路パターンの疎密に応じてレジストパターンの線幅などに変化を生じさせてしまう現象であり、従来の描画装置においては回路のより一層の微細化に伴ってこうした近接効果の影響がより顕在化するといった不都合があった。
【０００４】
そこで、前記近接効果を補正する有効な方法の一つとして、近接効果によるレジストパターンの線幅変化などを低減するのに最適な電子ビームの照射量（最適照射量あるいは最適ドーズ量などと呼ばれる）を、回路パターンの疎密に応じて決定する照射量補正法が知られている。具体的には、回路パターンが密の箇所では実質的な電子ビームの照射量が超過するので照射時間を減らすように制御する一方で、回路パターンが疎の箇所では実質的な電子ビームの照射量が不足するので照射時間を増すように制御する。こうすることによって、近接効果によるレジストパターンの線幅変化を低減するようになっている。このような電子ビームの最適照射量を決定する方法の一例を挙げるとすれば、例えば下記に示す非特許文献１に記載の技術などがある。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】M.Parikh,J.Appl.Phys50 (1979),p4371-4383
【０００６】
上記非特許文献１に記載の方法は、電子ビームの照射量（ドーズ量）と試料上に形成される回路パターンの予め決められた位置（これを評価点と呼ぶ）における感光量との関係を数１に示す行列式で表現しておき、この行列の逆行列を求めることによって、前記各位置での電子ビームの最適照射量を求めるという方法である（自己整合法あるいは行列法などと呼ばれる方法である）。この数１に示される行列式（ＦＤ＝Ｅ）は、照射量を設定する複数の矩形（試料上に形成する回路パターンを構成する各種図形をビーム照射単位に分割したもの）それぞれからの影響を加味し、数２及び数３に示すようなエネルギーの散乱の式Ｆijを用いて「すべての評価点での蓄積エネルギーが等しい」という目的条件を方程式で表したものである。
【数１】
2014029888.tif 000003
ただし、数１はｎ個の評価点、ｍ個の矩形があるときの方程式を例に示したものである。Ｄiは矩形ｉのドーズ量であり、Ｅthreasholdは各評価点での蓄積エネルギーの目標値（共通の値：定数）である。
【０００７】
上記相関関数行列の要素Ｆij（≡＝１，・・・ｍ、ｊ＝１，・・・ｎ）は、例えば次に示す数２に従って計算される。
【数２】
2014029888.tif 000004
上記相関関数行列の要素Ｆijは、矩形ｉからの評価点ｊでのエネルギーの散乱の作用の大きさ（蓄積エネルギー）を示す。
【０００８】
ここで、例えば図９左図に示されるような四辺が座標（l,b,r,t）それぞれの位置にある長方形（具体的には描画対象とされたパターン分割後の矩形）の周囲の任意の点（x,y）における蓄積エネルギーｅ（x,y）は、試料上のある一点に電子ビームを照射したときの最終的なエネルギーの分布を示すＰＳＦ関数を積分したものにより求めることができる（数３参照）。また、数３に示されるＰＳＦ関数がガウス分布で近似できるならば、この数３は数４に示されるようなガウス分布の積としてｅｒｆ関数（誤差関数；Error Function）の積により表現することができる。
【数３】
2014029888.tif 000005
【数４】
2014029888.tif 000006
【０００９】
上記数４に示したガウス分布の組み合わせで表現される近似式は、積分が容易であるために速く計算できるという利点がある。しかし、最近の半導体デバイスに要求される回路寸法及び回路線幅の微細化のより一層の進行に伴い、上記したようなガウス分布による近似ではその精度が不十分になってきた（すなわち、ＰＳＦ関数がガウス分布で近似できない）。
【００１０】
ところで、ＰＳＦ関数がガウス分布で近似できない場合には、数４に示した近似式を使えないことになる（ただし、このことは自己整合法を利用する場合に限らない）。そこで、ＰＳＦ関数がガウス分布で近似できない場合には、図９右図に示すように四辺が座標（l,b,r,t）にそれぞれ位置する長方形（描画対象の矩形）の一面を多数の点（座標（Xi,Yj））の集合と看做して、当該長方形の周囲の任意の点（x,y）における蓄積エネルギーｅ（x,y）を数５に従って計算すればよい。
【数５】
2014029888.tif 000007
ここで、ｓは隣接する点と点との間隔（グリッド間隔）、ｗは長方形（描画対象の矩形）の幅、ｈは長方形（描画対象の矩形）の高さである（図９右図参照）。
【発明の概要】
【発明が解決しようとする課題】
【００１１】
上記した数５に従って任意の点（x,y）における蓄積エネルギーｅ（x,y）を計算する方法の利点は、個々の点のグリッド間隔ｓを狭く（細かく）に設定すれば計算精度が上がるという点にある。したがって、微細化された回路パターンを試料上に描画する必要がある現状において、微細化の進行にあわせて補正精度のよい電子ビームの最適照射量を得るためには、グリッド間隔ｓを狭く設定すればよいことは理解できる。
【００１２】
しかし、上記数５に従って任意の点（x,y）における蓄積エネルギーｅ（x,y）を計算する方法の短所は、グリッド間隔ｓを狭く設定すると計算にかかる時間が増えてしまう点にある。したがって、昨今の回路パターンの微細化に応じた十分な補正精度を持つ電子ビームの最適照射量を得るためにグリッド間隔ｓを狭く設定する必要があるような場合には、前記電子ビームの最適照射量を求める際に用いる蓄積エネルギーの算出に関してその計算に莫大な時間がかかることになり、時間の制約の点から特にＬＳＩ用パターンなどの高集積化された半導体デバイスに用いることは難しい、という問題があった。
【００１３】
本発明は上述の点に鑑みてなされたもので、近接効果による回路パターンの寸法変化を低減するのに最適な電子ビームの照射量を求めるために用いられるパターン分割後の各矩形からの蓄積エネルギーの算出に関し、高速の処理速度つまりは計算スピードでかつ精度よく前記各矩形からの蓄積エネルギーを求めることのできる描画方法及び描画装置並びにプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１４】
本発明に係る描画方法は、試料上に描画する回路パターンを分割する任意の矩形毎に、近接効果による前記パターンの寸法変化を低減するのに最適な電子ビームの照射量を求めるため、電子ビームの照射対象とされた第１の矩形を基準とする所定範囲内にある前記第１矩形以外の第２の矩形からの実質的な電子ビームの照射量である蓄積エネルギーを前記第２の矩形毎に求める描画方法であって、所定位置に電子ビームを照射した際における複数の各要素毎のエネルギー分布からなる離散的な関数マップを求めるステップと、前記求めた離散的な関数マップに従って前記第２の矩形からの蓄積エネルギーを求めるステップとを備える。
【００１５】
本発明の好ましい実施例として、前記離散的な関数マップを求めるステップは、離散的な関数マップを、
2014029888.tif 000008
（ただし、ｓ；各要素間のグリッド間隔、Ｎ；Ｒ／ｓ（ＲはＰＳＦ定義域））
により求めることを特徴とする。
【００１６】
また、前記第２の矩形からの蓄積エネルギーを求めるステップは、第２の矩形からの蓄積エネルギーを、
2014029888.tif 000009
（ただし、x；第１矩形上の任意の点のＸ座標位置、y；第１矩形上の任意の点のＹ座標位置、xi；関数マップ上の各要素のｘ座標位置、yj；関数マップ上の各要素のｙ座標位置、ｓ；各要素間のグリッド間隔、Ｎ；Ｒ／ｓ（ＲはＰＳＦ定義域））
に基づいて求めることを特徴とする。
【００１７】
さらに、前記第２の矩形からの蓄積エネルギーを求めるステップは、前記第２の矩形に従って任意の大きさからなる複数の第３の矩形の組み合わせを特定するステップと、離散的な関数マップに従って前記特定した複数の第３の矩形毎に蓄積エネルギーを求めるステップと、該求めた第３の矩形からの蓄積エネルギーを前記組み合わせ態様に従って加減算するステップとを含んでなることを特徴とする。
【００１８】
こうすることによれば、予め求めておいた関数マップつまりは各要素のエネルギー分布に従って簡単な計算を行うだけでよいことから、パターン分割後の各矩形からの蓄積エネルギーを高速につまりは時間をかけることなく求めることができるようになる。また、これによれば、個々の要素のグリッド間隔を狭く設定した場合であっても、短時間で精度のよい各矩形からの蓄積エネルギーを求めることができる。
【００１９】
本発明は方法の発明として構成し実施することができるのみならず、装置の発明として構成し実施することができる。また、本発明は、コンピュータまたはＤＳＰ等のプロセッサのプログラムの形態で実施することができるし、そのようなプログラムを記憶した記憶媒体の形態で実施することもできる。
【発明の効果】
【００２０】
本発明によれば、予め複数の各要素毎のエネルギー分布からなる離散的な関数マップを求めておき、この関数マップに従って各矩形からの蓄積エネルギーを求めるようにしたことから、これにより高速の処理速度つまりは計算スピードでかつ精度よくパターン分割後の各矩形からの蓄積エネルギーを求めることができる、という効果を奏する。
【図面の簡単な説明】
【００２１】
【図１】本発明に係る描画方法を適用した描画装置の全体構成の一実施例を示す概念図である。
【図２】最適照射量計算部における矩形ごとの電子ビームの最適照射量算出処理の一実施例を示すフローチャートである。
【図３】蓄積エネルギー計算部での各評価点における蓄積エネルギー算出処理の一実施例を示すフローチャートである。
【図４】近接効果補正により最適照射量を算出するチップ範囲の一例を示す概念図である。
【図５】離散的な「関数マップ」の一例を示す概念図である。
【図６】ＰＳＦ定義域を決めるシミュレーション結果の一例を示す概念図である。
【図７】離散的な関数マップを用いて任意の点における矩形からの蓄積エネルギーを求める手順を説明するための概念図である。
【図８】ＰＳＦ定義域外のエネルギー分布を参照する必要がある場合に用いる蓄積エネルギーの式の一例を示す概念図である。
【図９】離散化による蓄積エネルギーを計算する従来例を示す概念図である。
【発明を実施するための形態】
【００２２】
以下、この発明の実施の形態を添付図面に従って詳細に説明する。
【００２３】
図１は、本発明に係る描画方法を適用した描画装置の全体構成の一実施例を示す概念図である。ここに示された描画装置は一例として電子ビーム描画装置を示したものであって、１０は試料室、１１はターゲット（試料）、１２は試料台、２０は電子光学鏡筒、２１は電子銃、２２ａ〜２２ｅは各種レンズ系、２３〜２６は各種偏向系、２７ａはブランキング板、２７ｂ，２７ｃはビーム成形用アパーチャマスクを示している。また、３１は試料台駆動回路部、３２はレーザ側長系、３３は偏向制御回路部、３４はブランキング制御回路部、３５は可変成形ビーム寸法制御回路部、３６はバッファメモリ及び制御回路部、３７は制御計算部、３８は最適照射量計算部、４０はＣＡＤシステムを示している。
【００２４】
図１に示した電子ビーム描画装置の動作を簡単に説明すると、電子銃２１から照射された電子ビームはブランキング用偏向器２３によりオン／オフされる。本装置はこの際の電子ビーム照射時間の長短を調整することによって、試料台１２に載置されたターゲット１１への照射位置に応じた電子ビームの照射量を変化させることができるようになっている。ブランキング板２７ａを通過した電子ビームは、ビーム成形用偏向器２４及びビーム成形用アパーチャマスク２７ｂ，２７ｃにより矩形ビームに成形されると共に、その矩形の寸法が可変される。そして、この矩形状に成形された電子ビームは走査用偏向器２５，２６によりターゲット１１上で偏向走査され、このビーム走査によりターゲット１１上に所望のパターンが描画されるものとなっている。すなわち、ターゲット１１上に描画される所望のパターンは複数の矩形（第１の矩形）の組み合わせに分割され、個々の矩形（第１の矩形）ごとに電子ビームが照射される。
【００２５】
前記矩形ごとの電子ビームの最適照射量（つまりは電子ビーム照射時間の長短）は、ＣＡＤシステム４０で作成された電子ビーム露光用の元データに基づき最適照射量計算部３８によって算出される。最適照射量計算部３８は例えばＣＰＵ、ＲＯＭ、ＲＡＭ等を含んでなるコンピュータであって、近接効果を補正して矩形ごとに前記電子ビームの最適照射量を算出する。ここでは、共役勾配法を利用して上記数１で表せるような行列式（ここでは便宜上、Ａｘ＝ｂと記載する）を解くことによって、矩形ごとの電子ビームの最適照射量を算出する最適照射量計算部３８を例に示す。
【００２６】
最適照射量を算出するために共役勾配法を利用すると、行列式Ａｘ＝ｂは相関関数行列（Ａ）とベクトル（ｘ）の掛け算の繰り返しにより解くことができる。そこで、共役勾配法による行列式（Ａｘ＝ｂ）を解くための処理手順について以下に説明する。まず、初期値としてｒk＝ｂ−Ａｘk、ｐk＝ｒk（ｋ＝０）とする。これらｐkとｒkは中間変数ベクトルであって、ｋは繰り返し計算の回数を示す。次に、以下に示す数８〜数１３までの計算を所定の繰り返し条件を満たすまで繰り返し行うことによって行列式Ａｘ＝ｂを解きｘを求める。なお、以下に示す各数式内におけるドットつきの括弧表記（・）は内積を表す。
【数８】
2014029888.tif 000010
【数９】
2014029888.tif 000011
【数１０】
2014029888.tif 000012
【数１１】
2014029888.tif 000013
【数１２】
2014029888.tif 000014
【数１３】
2014029888.tif 000015
【００２７】
ところで、上記計算手順の過程で現われるＡｐkの「Ａ」はｍ×ｎの大きさの行列であり、「ｐk」はｍの大きさのベクトルである（数１参照）。したがって、通常の方法でＡｐkを計算するにはｍ×ｍ回の掛け算が必要となり、しかも繰り返し計算ごとにそうした計算が行われるために、当該処理の計算量を削減することが処理時間を高速化するために重要であることが理解できる。また、「Ａ」の要素数は矩形数の２乗になることに鑑みれば、上記計算手順をそのまま用いたのでは多大な計算時間が必要とされて都合が悪い。
【００２８】
そこで、行列式（Ａｘ＝ｂ）を高速に解くために、本実施形態にかかる最適照射量計算部３８は蓄積エネルギー計算部３９を有する。この蓄積エネルギー計算部３９では、上記計算手順の過程で現われるＡｐkが各矩形のドーズ量が「ｐk」であるときの各評価点での蓄積エネルギーに等しいことに注目して、Ａｐkの値を（後述のように１つの評価点を有する）矩形要素ごとに１つの蓄積エネルギーで管理する。すなわち、先の行列［Ａ］は系の振る舞いを表すものであるので、行列式Ａｘ＝ｂとは各矩形のドーズ量が［ｘ］であるときの各評価点での蓄積エネルギーが［ｂ］になっていると言うことを意味する。したがって、蓄積エネルギー計算部３９では各評価点における蓄積エネルギーを計算すればよく、評価点の数及び矩形の数に比例して増減する多数の要素からなる巨大な行列Ａを管理してＡｐkの計算をわざわざ行う必要がない。
【００２９】
ここで、蓄積エネルギー計算部３９では、予め計算しておいた離散的な関数マップ（後述する図５参照）を用いて複数の四半平面（第３の矩形に相当）からのエネルギー分布をそれぞれ求めておき、それらを組み合わせて任意の矩形（第２の矩形）からのエネルギー分布に従う任意の点（詳しくは第１の矩形における評価点）における蓄積エネルギーを求めることのできるようにしている（後述する図７参照）。こうすることによって、前記矩形（第２の矩形）ごとの電子ビームの最適照射量の算出をより高速に処理することのできるようにしている。具体的な計算手順については後述する。
【００３０】
以下では、別途示した図４〜図７を適宜に参照しながら、矩形（第２の矩形）ごとの電子ビーム最適照射量（最適ドーズ量）の算出手順について説明する。図２は、図１に示した最適照射量計算部３８における矩形ごとの電子ビーム最適照射量（最適ドーズ量）の算出処理の一実施例を示すフローチャートである。
【００３１】
ステップＳ１は、ＣＡＤシステム４０から取得した電子ビーム露光用の元データに基づき特定される複数の図形Ｚ１〜Ｚ７（回路パターンを構成する図形；図４参照）それぞれを、適切に複数の矩形に分割する。例えば、図４に示すように、各図形Ｚ１〜Ｚ７を後方散乱径の１／１０程度の大きさの複数の矩形に分割する。この際に、図４において示される図形Ｚ１，Ｚ３，Ｚ４，Ｚ５などのようにもともとの形状が四角形以外である図形である場合には、前記サイズの複数の矩形を組み合わせてなるように当該図形を分割する。ここに示す例では、図形Ｚ１が３個の矩形に、図形Ｚ３が２個の矩形に、図形Ｚ４が４個の矩形に、図形Ｚ５が２個の矩形にそれぞれ分割されている。このようにして、各図形Ｚ１〜Ｚ７は大きさが異なる１乃至複数の矩形に分割される。こうした図形Ｚ１〜Ｚ７を１乃至複数の矩形に分割する方法は公知のどのようなものであってもよいことから、ここでの詳しい説明を省略する。
【００３２】
ステップＳ２は、前記分割した各矩形につき１つのエネルギー評価点を決める。図４に示すように、通常のPECでは各矩形を形成する最も長い外周部分の一片の中央位置をエネルギー評価点に決める（図中において黒丸で示す）。この点、各矩形を形成する各片すべての中央位置をエネルギー評価点に決定する（つまり４つの評価点がある）、従来知られた自己整合法とは異なる。なお、Gray Scale PECである場合には、各矩形の中心位置をエネルギー評価点に決めるとよい。
【００３３】
ステップＳ３は、数１に示されるような行列式Ａｘ＝ｂを解くための初期化処理として、各矩形の（初期）ドーズ量ｘ0に初期値「０」を、ｐ0及びｒ0（＝ｂ−Ａｘ0）に適当なエネルギー目標値を設定する。ここで、ｂはｍ元の列ベクトルである。また、上記したようにｐkとｒkは中間変数ベクトルであり、ｋは繰り返し計算の回数を示すものである。前記エネルギー目標値として、通常のPECである場合には各矩形が属する図形Ｚ１〜Ｚ７外周部分のエネルギー目標値を、Gray Scale PECである場合には各矩形が属する図形Ｚ１〜Ｚ７（を含むレイヤ）の蓄積エネルギー目標値をそれぞれ設定する。
【００３４】
ステップＳ４は、図９に示すように面を多数の点の集合と看做して離散的な計算を行うために、離散化グリッド間隔ｓ及びＰＳＦ定義域を決定する。離散化グリッド間隔ｓは、ＰＳＦ関数の変化に対して十分に小さい値に決定する。他方、ＰＳＦ定義域は後述する「関数マップ」（図５参照）の定義域を定義するためのものであって、例えば図６に示されるようなＰＳＦ関数（近似式）を用いたシミュレーション結果に従って決定される。具体的には、エネルギー値が所定値以上である範囲（電子ビーム入射点からの距離）をＰＳＦ定義域として決定する。あるいは、シミュレーション結果に従ってユーザが任意の範囲をＰＳＦ定義域として決定するようにしてもよい。
【００３５】
上記のシミュレーションを行うために用いられるＰＳＦ関数（近似式）としては、例えば数１４に示すようなダブルガウシアン近似式などがある。ここでは、単にＰＳＦ定義域を決めるだけであるので、ＰＳＦ関数をダブルガウシアン近似式などで近似してよい。
【数１４】
2014029888.tif 000016
【００３６】
この数１４に示されるＰＳＦ関数は、試料上のある一点に電子ビームを照射したときの最終的なエネルギーの分布を示し、ｒは電子ビーム入射点からの距離、Ｃは定数、ηは電子ビームの前方散乱によるレジストの感光量と後方散乱によるレジストの感光量との比（近接効果補正係数）、α及びβはそれぞれが加速電圧に応じて決まる前方散乱の広がり（前方散乱径）又は後方散乱の広がり（後方散乱径）を表す所定値である。前記α及びβの各値は、例えば加速電圧が20KeVである場合に（27nm,2um）、50KeVである場合に（30nm,10um）、100KeVである場合に（10nm,32um）である。
【００３７】
なお、ＰＳＦ定義域を決めるシミュレーションを行うために用いるＰＳＦ関数（近似式）としては上記したダブルガウシアン近似式に限らず、公知の他の近似式を用いてよいことは言うまでもない。また、例えば主に電子ビームの加速電圧や基盤の材質に従って、シミュレーションに用いるＰＳＦ関数（近似式）を決めるようにしてもよい。
【００３８】
ステップＳ５は、離散的な「関数マップＥ」の初期化を実行する。ここでは、あらかじめ２Ｎ×２Ｎの大きさの行列［Ｅ］を確保しておき、関数マップ上の任意の点（座標（I,J））に電子ビームを照射した際におけるエネルギー分布（Ｅij）を当該行列の要素毎に初期化する計算を数１５に従って行う。
【数１５】
2014029888.tif 000017
上記Ｎは、Ｎ＝Ｒ（ＰＳＦ定義域）／ｓ（離散化グリッド間隔）により求められる。なお、この数１５で用いるＰＳＦ関数は、ＰＳＦ定義域を決定する際に用いたのと同じ近似式（例えば数１４に示したダブルガウシアン近似式など）を用いればよい。また、上記エネルギー分布（Ｅij）を初期化する際には、個々の点のスポット配置や密度等を考慮してエネルギー計算を行うようにしてよい。
【００３９】
図５に、離散的な「関数マップＥ」の一例を示す。この図５では、ＤIJ＝１となる要素を黒丸で示し、ＤIJ＝０となる要素を白丸で示している。この離散的な「関数マップＥ」は２Ｒ×２Ｒ（ＲはＰＳＦ定義域）の大きさからなる面を点（要素）の集合と看做し、各要素に電子ビームを照射したときのＤIJ＝１となる要素のみを含んでなる四半平面（図中において座標軸ｘｙの右上の範囲）からのエネルギー分布Ｅijを個々の要素毎に表したものである。
【００４０】
ステップＳ６は、行列式（Ａｘ＝ｂ）を共役勾配法により解く処理手順に現れるＡｐk（数８参照）を蓄積エネルギー計算部３９で計算する。上述したように、蓄積エネルギー計算部３９で上記Ａｐkを計算するには、前記決定したエネルギー評価点における蓄積エネルギーｑkを算出すればよい。蓄積エネルギーｑkは、レジストが電子の衝突によって得たエネルギー（分布）である。
【００４１】
ここで、蓄積エネルギー計算部３９での各評価点における蓄積エネルギーｑkの算出手順について、図３を用いて説明する。図３は、蓄積エネルギー計算部Ａでの各評価点における蓄積エネルギー算出処理の一実施例を示すフローチャートである。
【００４２】
ステップＳ２１は、１つのエネルギー評価点（電子ビームの照射対象と仮定された分割後の矩形；第１の矩形に相当）を特定する。ステップＳ２２は、前記特定した１つのエネルギー評価点に影響を及ぼしうる所定範囲（影響範囲）に含まれる、前記図形Ｚ１〜Ｚ７を分割した矩形のうち１つの矩形（第２の矩形に相当）を処理対象に決定する。ここで、前記特定した１つのエネルギー評価点に影響を及ぼしうる範囲である「影響範囲」としては、例えば矩形（第１の矩形）の外周部を前方散乱径の４倍程度の長さだけ広げた範囲、あるいは500nm程度の範囲など任意の範囲に決定してよい。その際には、前方散乱による影響や後方散乱による影響をどの程度加味するかに応じて「影響範囲」を決定するとよい。
【００４３】
図４では一例として、図形Ｚ２（第１の矩形）のエネルギー評価点Ｈ２が特定されたと仮定した場合における、前記エネルギー評価点Ｈ２に影響を及ぼしうる「影響範囲」を点線で示した四角枠で便宜的に表している。そして、この「影響範囲」に含まれる図形Ｚ１を分割した矩形Ｚ１ａ及び矩形Ｚ１ｂのいずれかが順に、処理対象の矩形（第２の矩形）として決定されることになる。
【００４４】
ステップＳ２３は、前記決定した矩形（第２の矩形）を形成する複数の任意の大きさからなる四半平面（第３の矩形）の組み合わせを特定する。ステップＳ２４は、離散的な関数マップＥ（図５参照）を用いて前記特定した複数の四半平面それぞれからの蓄積エネルギーを求める。ステップＳ２５は、前記求めた複数の四半平面それぞれからの蓄積エネルギーに基づき、前記特定した１つのエネルギー評価点における前記決定した矩形からの蓄積エネルギーを求める。
【００４５】
ここで、上記したステップＳ２３〜ステップＳ２５に従う１つのエネルギー評価点における矩形からの蓄積エネルギーを求める手順について、具体的に説明する。例えば、図４に示したエネルギー評価点Ｈ２における矩形Ｚ１ａからの蓄積エネルギーｅ（x,y）を求める場合を例に、図７を用いて説明する。ここでは、矩形Ｚ１ａの四辺が座標（l,b,r,t）にそれぞれ位置するものと仮定する。図７は、離散的な関数マップを用いて任意の点における矩形からの蓄積エネルギーを求める手順を説明するための概念図である。
【００４６】
ここで、上記した離散的な関数マップＥ（図５参照）を用いて任意の点（x,y）における蓄積エネルギーF（x,y）を求める式を数１６に示す。
【数１６】
2014029888.tif 000018
【００４７】
xiは関数マップＥ上の各要素のｘ座標位置、yjは関数マップＥ上の各要素のｙ座標位置を示す。この数１６から理解できるように、任意の点（x,y）における蓄積エネルギーF（x,y）は、離散的な関数マップＥにおいて前記任意の点（x,y）を囲む周囲４つの点（要素）の各エネルギー分布（Ｅij）に基づき双一次補間によって求められる。
【００４８】
まず、上記数１６に従って図中×印で示した座標位置（x,y）にあるエネルギー評価点Ｈ２の蓄積エネルギーF（x,y）を求めると、F（x-l,y-b）と表すことができ、これは四半平面Ｂ１に含まれる各要素のエネルギー分布に基づくものである。すなわち、上記数１６により求められる任意の点（x,y）における蓄積エネルギーF（x,y）は、ある四半平面全体に含まれる全ての要素のエネルギー分布（Ｅij）に従って求められる。ただし、上記四半平面Ｂ１は矩形Ｚ１ａ（図中斜線で示す範囲）を含むより範囲の広い矩形であるため、任意の点（x,y）における矩形Ｚ１ａからの蓄積エネルギーを求めるためには、数１６をそのまま使用するわけにいかず、当該矩形Ｚ１ａに含まれない範囲内の各要素のエネルギー分布（Ｅij）を加味しない必要がある。
【００４９】
そこで、本実施形態においては、前記矩形Ｚ１ａを形成する複数の任意の大きさからなる複数の四半平面（図７の例ではＢ１,Ｂ２,Ｂ３,Ｂ４の４つの四半平面）の組み合わせを特定し（ステップＳ２３参照）、該特定した複数の四半平面それぞれにおいて任意の点（x,y）における蓄積エネルギーを数１６に従って求め（ステップＳ２４参照）、それぞれの計算結果を加減算して前記矩形Ｚ１ａからの蓄積エネルギーｅ（x,y）を求めるようにしている（ステップＳ２５参照）。ここでは、四半平面Ｂ１のうちｙ座標がｔより大きい範囲の四半平面Ｂ２、四半平面Ｂ１のうちｘ座標がｒより大きい範囲の四半平面Ｂ３を特定し、特定した四半平面毎に任意の点（x,y）における蓄積エネルギーF（x,y）を求める。これら矩形Ｚ１ａを形成する複数の四半平面Ｂ１,Ｂ２,Ｂ３は、図７に示すような関係になる。
【００５０】
上記四半平面Ｂ２,Ｂ３からの蓄積エネルギーを数１６に従って求めると、それぞれ、F（x-l,y-t），F（x-r,y-b）と表すことができる。そこで、図７に示されるような当該矩形Ｚ１ａに含まれない範囲に含まれる各要素のエネルギー分布（Ｅij）を加味しないようにするために、四半平面Ｂ１に基づく蓄積エネルギーF（x-l,y-b）から四半平面Ｂ２,Ｂ３に基づく蓄積エネルギーF（x-l,y-t）、F（x-r,y-b）を減算する。ただし、この減算によると、図７に示されるそれぞれの範囲のうち重なり合った一部範囲（図中において点線で囲んだ範囲）が２度減算されることになり都合が悪い。そのため、調整分として前記一部範囲である四半平面Ｂ４に基づく蓄積エネルギーF（x-r,y-t）を上記減算結果に加算する。このようにして、例えば矩形Ｚ２のエネルギー評価点Ｈ２における矩形Ｚ１ａからの蓄積エネルギーｅ（x,y）が求められる。
【００５１】
矩形Ｚ２のエネルギー評価点Ｈ２における矩形Ｚ１ａからの蓄積エネルギーｅ（x,y）を求める最終的な式の一例を示すと、数１７のように表すことができる。
【数１７】
2014029888.tif 000019
【００５２】
この数１７に示す各項は、数１６に示したような４つの要素の各エネルギー分布（Ｅij）によって求められ、またそれら各要素のエネルギー分布（Ｅij）は数１５に従って初期化によって予め計算されている値である。この数１７に従えば、任意の点（x,y）における蓄積エネルギーF（x,y）を簡単な計算で求めることができる。そうであるならば、補正精度のよい電子ビームの最適照射量を得るために個々の点のグリッド間隔ｓを狭く設定する必要があるような場合であったとしても、精度の高い任意の点（x,y）における蓄積エネルギーF（x,y）を高速につまりは時間をかけることなく短時間に求めることができるようになる。
【００５３】
ところで、矩形からの蓄積エネルギーを求める任意の点がＰＳＦ定義域外の位置（すなわち関数マップ外）にある場合がある。任意の点がＰＳＦ定義域外の位置にあって、ＰＳＦ定義域外のエネルギー分布（Ｅij）を参照する必要がある場合には、ｘやｙをＰＳＦ定義域の境界まで切り上げた（又は切り捨てた）時の蓄積エネルギーの式を用いればよい。図８にその例を示す。
【００５４】
図８は、ＰＳＦ定義域外のエネルギー分布を参照する必要がある場合に用いる蓄積エネルギーの式の一例を示す概念図である。この図８において太線で示す四角枠の外がＰＳＦ定義域外に相当し、前記四角枠を囲む周囲にｘやｙをＰＳＦ定義域の境界まで切り上げた（又は切り捨てた）時の蓄積エネルギーの式を便宜的に示している。この図８に示した各蓄積エネルギーの式F（x,y）によれば、ＰＳＦ定義域外の値を参照する必要がある場合にはＰＳＦ定義域の境界の値を使えばよいことが理解できる。
【００５５】
図３の説明に戻って、ステップＳ２６は前記処理過程において算出された前記「影響範囲」に含まれる各矩形ごとに求められる前記蓄積エネルギーを累算する。ステップＳ２７は、前記「影響範囲」内に含まれる全ての矩形について上記蓄積エネルギー計算を行ったか否かを判定する。「影響範囲」と重なる位置関係にある矩形全てについて上記蓄積エネルギー計算を行っていないと判定した場合には（ステップＳ２７のＮＯ）、ステップＳ２２の処理に戻って前記「影響範囲」内に含まれる他の矩形についても同様に蓄積エネルギーを算出する。他方、「影響範囲」と重なる位置関係にある矩形全てについて上記蓄積エネルギー計算を行ったと判定した場合には（ステップＳ２７のＹＥＳ）、全てのエネルギー評価点について上記処理を行ったか否かを判定する（ステップＳ２８）。
【００５６】
全てのエネルギー評価点について上記処理を行っていないと判定した場合には（ステップＳ２８のＮＯ）、ステップＳ２１の処理へ戻る。全てのエネルギー評価点について上記処理を行ったと判定した場合には（ステップＳ２８のＹＥＳ）、本処理を終了する。このようにして、図形分割後の各矩形ごとに１つずつ決定される評価点ごとに蓄積エネルギーを求める。すなわち、本実施形態では予め求めておいたエネルギー分布（離散的な関数マップＥ）に基づき、組み合わせによって矩形を構成する四半平面毎にエネルギー分布を個々に計算して、それらを累算することによって、エネルギー評価点の蓄積エネルギーを高速な計算スピードで求めることができるようにしている。
【００５７】
図２の説明に戻って、ステップＳ７は、全ての図形Ｚ1〜Ｚ７の（ｐk・ｑk）と（ｒk・ｒk）を合計する。すなわち、ΣｐkｑkとΣｒkｒkを求める。ステップＳ８は、αk＝Σｒkｒk／Σｐkｑk（共役勾配法の処理手順における数８の処理に相当）とする。ステップＳ９は、全ての図形のｘとｒを
ｘk+1＝ｘk＋αkｐk（共役勾配法の処理手順における数９の処理に相当）
ｒk+1＝ｒk−αkｐk（共役勾配法の処理手順における数１０の処理に相当）
とする。
【００５８】
ステップＳ１０は、全ての図形の（ｒk+1・ｒk+1）を合計する（以下では、これをΣｒｒＮｅｘｔと記載）。ステップＳ１１は、ΣｒｒＮｅｘｔが十分に小さいか否かつまりは計算誤差が予め設定された許容値よりも小さくなったか否かを判定する。ΣｒｒＮｅｘｔが十分に小さいと判定した場合には（ステップＳ１１のＹＥＳ）、繰り返し計算を終了して列ベクトルｘを各図形の最適補正量として出力する（ステップＳ１５）。
【００５９】
他方、ΣｒｒＮｅｘｔが十分に小さくないと判定した場合には（ステップＳ１１のＮＯ）、ステップＳ１２及びステップＳ１３及びステップＳ１４の処理を行ったうえで上記ステップＳ６の処理に戻り、上記した各計算処理を繰り返し実行する。ステップＳ１２は、βk＝ΣｒｒＮｅｘｔ／Σｒkｒk（共役勾配法の処理手順における数１１の処理に相当）とする。ステップＳ１３は、全ての図形のｐをｐk+1＝ｒk+1＋βkｐk（共役勾配法の処理手順における数１２の処理に相当）に更新する。ステップＳ１４は、繰り返し回数ｋに「１」を加算する（共役勾配法の処理手順における数１３の処理に相当）。
【００６０】
なお、上記の計算処理を繰り返すか否かを判定する繰り返し条件としては、上記したようなΣｒｒＮｅｘｔの大小を条件とするものに限らない（ステップＳ１１参照）。例えば、シミュレーションなどによって予め収束する繰り返し回数を調査した上で繰り返し回数を設定する方法、繰り返し計算をしても誤差が変動しなくなったことを確認する方法などがある。なお、理論上、繰り返し回数の最大上限回数はｍ回である。
【００６１】
以上のように、本実施形態においては近接効果を補正した電子ビームの最適照射量を求めるのに先立ち算出する必要がある、試料上に描画するパターンを分割した各矩形からの蓄積エネルギーの算出方法に関し、予め複数の各要素毎のエネルギー分布からなる離散的な関数マップを求めておき、この関数マップに従って各矩形からの蓄積エネルギーを求めるようにした。これによれば、予め求めておいた関数マップつまりは各要素のエネルギー分布に従って簡単な計算を行うだけでよいことから、任意の点における各矩形からの蓄積エネルギーを高速につまりは時間をかけることなく求めることができるようになる。また、個々の要素のグリッド間隔を狭く設定した場合であっても、短時間で精度のよい矩形からの蓄積エネルギーを求めることができる、という利点もある。
【００６２】
以上、図面に基づいて実施形態の一例を説明したが、本発明はこれに限定されるものではなく、様々な実施形態が可能であることは言うまでもない。例えば、上述した実施例においては本発明に係る描画方法を可変成形ビーム方式の電子ビーム描画装置に適用した例を示したが、これ以外の方式の描画装置にも適用できる。また、電子ビームの代わりにイオンビームを用いたイオンビーム描画装置に適用することも可能である。さらに、本発明は電子ビーム描画装置の使用目的に限定するものではない。例えば、ウェハ上に直接レジストパターンを形成するという目的以外にも、Ｘ線マスクを作成する際、光ステッパ用マスク、レチクル等を作成する際にも利用可能である。その他、本発明の要旨を逸脱しない範囲で、種々変形して実施することが可能である。
【００６３】
なお、本発明に係る描画方法は、上述したようにして近接効果を補正した電子ビームの最適照射量を上記処理に従って求めるために用いられる各矩形からの蓄積エネルギーを算出する場合に適用することができるだけでなく、例えば、ユーザが描画予測に基づいて電子ビームの最適照射量を求めることのできるようにするために、任意に入力された各矩形に対する照射量に従って電子ビームが照射されたと仮定したときに形成される回路パターンを表示するシミュレーション処理を行う際に用いられる各矩形からの蓄積エネルギーを算出する場合にも適用できることは言うまでもない。つまり、エネルギー分布をシミュレーションする際に用いてよい。
なお、上述した実施例においては共役勾配法を利用して電子ビームの最適照射量を算出する方法を示したがこれに限らず、電子ビームの最適照射量の導出過程において後述するような蓄積エネルギー算出過程を含むものであれば、自己整合法などの従来知られた他のどのような方法を採用してもよい。ただし、上記したように自己整合法などの従来知られた他の方法を利用するよりも共役勾配法を利用した方が、より計算時間を短くすることができるようになる点で有利である。
【符号の説明】
【００６４】
１０…試料室
１１…ターゲット（試料）
１２…試料台
２０…電子光学鏡筒
２１…電子銃
３１…試料台駆動回路部
３２…レーザー測長系
３３…偏向制御回路部
３４…ブランキング制御回路部
３５…可変成形ビーム寸法制御回路部
３６…バッファメモリ及び制御回路部
３７…制御計算部
３８…共益勾配法計算部
３９…蓄積エネルギー計算部
４０…ＣＡＤシステム
Ｈ２…エネルギー評価点
Ｚ１〜Ｚ７…図形（回路パターン）
【図１】
2014029888.tif 000022
【図２】
2014029888.tif 000023
【図３】
2014029888.tif 000024
【図４】
2014029888.tif 000025
【図５】
2014029888.tif 000026
【図６】
2014029888.tif 000027
【図７】
2014029888.tif 000028
【図８】
2014029888.tif 000029
【図９】
2014029888.tif 000030